改进的多道预测算子压制浅水多次波方法

在浅水情况下,由于观测数据中缺少近偏移距信息,水层多次波的压制面临挑战.利用多道预测算子压制水层多次波是浅水环境下压制多次波的重要方法之一,这种方法先从输入数据中估计出多道预测算子,再将预测算子和输入数据做褶积预测出水层相关多次波.然而,估计的多道预测算子很容易受噪声污染,从而影响多次波模型的精度.所以,我们提出了改进的多道预测算子压制浅水多次波方法.该方法先从数据中估计出多道预测算子,并利用估计的算子构建出精确的水层模型;然后,通过计算算子的走时信息、估计振幅信息、合成新算子三个步骤来修正原始的多道预测算子.修正的算子不仅不受噪声影响,还含有精确的走时信息、可靠的振幅信息;最后,该方法用修正的算子来预测多次波,并结合自适应相减,将预测的多次波从输入数据中去除.通过合成数据和实际资料的验证表明,相比于原始的多道预测算子压制多次波方法,改进的方法能够取得更好的压制效果.

     

High-order finite volume WENO schemes for the shallow water equations with dry states

The shallow water equations are used to model flows in rivers and coastal areas, and have wide applications in ocean, hydraulic engineering, and atmospheric modeling. These equations have still water steady state solutions in which the flux gradients are balanced by the source term. It is desirable to develop numerical methods which preserve exactly these steady state solutions. Another main difficulty usually arising from the simulation of dam breaks and flood waves flows is the appearance of dry areas where no water is present. If no special attention is paid, standard numerical methods may fail near dry/wet front and produce non-physical negative water height. A high-order accurate finite volume weighted essentially non-oscillatory (WENO) scheme is proposed in this paper to address these difficulties and to provide an efficient and robust method for solving the shallow water equations. A simple, easy-to-implement positivity-preserving limiter is introduced. One- and two-dimensional numerical examples are provided to verify the positivity-preserving property, well-balanced property, high-order accuracy, and good resolution for smooth and discontinuous solutions.     

Comparison between TVD-MacCormack and ADI-type solvers of the shallow water equations

A total variation diminishing (TVD) modification of the MacCormack scheme is developed for simulating shallow water dynamics on a uniform Cartesian grid. Results obtained using conventional and deviatoric forms of the conservative non-linear shallow water equations (SWEs) are compared for cases where the bed has a varying topography. The comparisons demonstrate that the deviatoric form of the SWEs gives more accurate results than the conventional form, in the absence of numerical balancing of the flux-gradient and source terms. A further comparison is undertaken between the TVD-MacCormack model and an alternating direction implicit (ADI) model for cases involving steep-fronted shallow flows. It is demonstrated that the ADI model is unable to predict trans-critical flows correctly, and artificial viscosity has to be introduced to remove spurious oscillations. The TVD-MacCormack model reproduces all flow regimes accurately. Finally, the TVD-MacCormack model is used to predict a laboratory-scale dyke break undertaken at Delft University of Technology. The predictions agree closely with the experimental data, and are in excellent agreement with results from an alternative Godunov-type model.     

Numerical resolution of well-balanced shallow water equations with complex source terms

This paper presents a well-balanced numerical scheme for simulating frictional shallow flows over complex domains involving wetting and drying. The proposed scheme solves, in a finite volume Godunov-type framework, a set of pre-balanced shallow water equations derived by considering pressure balancing. Non-negative reconstruction of Riemann states and compatible discretization of slope source term produce stable and well-balanced solutions to shallow flow hydrodynamics over complex topography. The friction source term is discretized using a splitting implicit scheme. Limiting value of the friction force is derived to ensure stability. This new numerical scheme is validated against four theoretical benchmark tests and then applied to reproduce a laboratory dam break over a domain with irregular bed profile.     

Positivity-preserving high order well-balanced discontinuous Galerkin methods for the shallow water equations

Shallow water equations with a non-flat bottom topography have been widely used to model flows in rivers and coastal areas. An important difficulty arising in these simulations is the appearance of dry areas where no water is present, as standard numerical methods may fail in the presence of these areas. These equations also have still water steady state solutions in which the flux gradients are nonzero but exactly balanced by the source term. In this paper we propose a high order discontinuous Galerkin method which can maintain the still water steady state exactly, and at the same time can preserve the non-negativity of the water height without loss of mass conservation. A simple positivity-preserving limiter, valid under suitable CFL condition, will be introduced in one dimension and then extended to two dimensions with rectangular meshes. Numerical tests are performed to verify the positivity-preserving property, well-balanced property, high order accuracy, and good resolution for smooth and discontinuous solutions.     

A unified formulation for the three-dimensional shallow water equations using orthogonal co-ordinates: theory and application

In this paper, the formulations of the primitive equations for shallow water flow in various horizontal co-ordinate systems and the associated finite difference grid options used in shallow water flow modelling are reviewed. It is observed that horizontal co-ordinate transformations do not affect the chosen co-ordinate system and representation in the vertical, and are the same for the three- and two-dimensional cases. A systematic derivation of the equations in tensor notation is presented, resulting in a unified formulation for the shallow water equations that covers all orthogonal horizontal grid types of practical interest. This includes spherical curvilinear orthogonal co-ordinate systems on the globe. Computational efficiency can be achieved in a single computer code. Furthermore, a single numerical algorithmic code implementation satisfies. All co-ordinate system specific metrics are determined as part of a computer-aided model grid design, which supports all four orthogonal grid types. Existing intuitive grid design and visual interpretation is conserved by appropriate conformal mappings, which conserve spherical orthogonality in planar representation. A spherical curvilinear co-ordinate solution of wind driven steady channel flow applying a strongly distorted grid is shown to give good agreement with a regular spherical co-ordinate model approach and the solution based on a β-plane approximation. Especially designed spherical curvilinear boundary fitted model grids are shown for typhoon surge propagation in the South China Sea and for ocean-driven flows through Malacca Straits. By using spherical curvilinear grids the number of grid points in these single model grid applications is reduced by a factor of 50–100 in comparison with regular spherical grids that have the same horizontal resolution in the area of interest. The spherical curvilinear approach combines the advantages of the various grid approaches, while the overall computational effort remains acceptable for very large model domains.     

南海北部深水盆地浅水流的地球物理特性及识别

浅水流(Shallow Water Flow,SWF)是深水环境海底浅部地层中超压的砂体流动,是对深水钻井最具破坏力的一种地质灾害,严重制约深水油气开发.为有效预测和防治浅水流,需要对浅水流地球物理特性进行研究,并在研究区内加以识别.本文借鉴国外主要深水盆地对浅水流问题的研究经验,对南海北部深水盆地潜在浅水流区域采取以属性判定、超压分析为主,振幅识别为辅的方法进行预测.精细层序地层学解释发现,南海北部深水盆地存在上新世以来的古珠江深水水道沉积体系和第四纪水道,这些水道砂体疏松未固结、孔隙度大、有效应力低、几乎表现出流体特性.基于遗传算法的混合反演方法发现,研究区存在典型的AVO响应,横波速度极低,低频特征明显,振幅强度弱,连续性较好,存在极性反转,高泊松比和高纵横波速度比.研究结果表明,南海北部陆坡具备浅水流发生的潜在条件,深水水道发育区为潜在的浅水流危险区,浅水流具有独特的地震响应特征,泊松比高达0.49,纵横波速度比约为3.5~9或更高,SWF层位对地震属性的敏感度VP/VS>AVO响应>泊松比.     

Balancing the source terms in a SPH model for solving the shallow water equations

A shallow flow generally features complex hydrodynamics induced by complicated domain topography and geometry. A numerical scheme with well-balanced flux and source term gradients is therefore essential before a shallow flow model can be applied to simulate real-world problems. The issue of source term balancing has been exhaustively investigated in grid-based numerical approaches, e.g. discontinuous Galerkin finite element methods and finite volume Godunov-type methods. In recent years, a relatively new computational method, smooth particle hydrodynamics (SPH), has started to gain popularity in solving the shallow water equations (SWEs). However, the well-balanced problem has not been fully investigated and resolved in the context of SPH. This work aims to discuss the well-balanced problem caused by a standard SPH discretization to the SWEs with slope source terms and derive a corrected SPH algorithm that is able to preserve the solution of lake at rest. In order to enhance the shock capturing capability of the resulting SPH model, the Monotone Upwind-centered Scheme for Conservation Laws (MUSCL) is also explored and applied to enable Riemann solver based artificial viscosity. The new SPH model is validated against several idealized benchmark tests and a real-world dam-break case and promising results are obtained.     

Weakly nonlinear shallow water magnetohydrodynamic waves

We consider an electrically conducting rotating fluid governed by the shallow water magnetohydrodynamic equations with no diffusion. We use an a priori asymptotic technique (the method of geometric optics or ray method) to study weakly nonlinear hydromagnetic waves. These waves are intermediate in length in the following sense: they are much longer than the fluid depth but much shorter than the radius of the earth. The time scale for the waves is much longer than that of the free surface oscillations and the approximation varies on an even longer timescale. The waves we are considering are studied in the beta plane approximation for an ambient magnetic field parallel to the equator which varies in the direction perpendicular to the equator. The leading order approximation gives a dispersion relation for the waves, which are generally found to be confined to bands about the equator as well as in bands at higher and lower latitudes. At the next order of approximation, a conservation law is found for the wave amplitude. We also obtain an equation governing the behavior of the leading order mean azimuthal velocity which is forced to grow linearly with time.     

A weighted surface-depth gradient method for the numerical integration of the 2D shallow water equations with topography

A finite volume MUSCL scheme for the numerical integration of 2D shallow water equations is presented. In the framework of the SLIC scheme, the proposed weighted surface-depth gradient method (WSDGM) computes intercell water depths through a weighted average of DGM and SGM reconstructions, in which the weight function depends on the local Froude number. This combination makes the scheme capable of performing a robust tracking of wet/dry fronts and, together with an unsplit centered discretization of the bed slope source term, of maintaining the static condition on non-flat topographies (C-property). A correction of the numerical fluxes in the computational cells with water depth smaller than a fixed tolerance enables a drastic reduction of the mass error in the presence of wetting and drying fronts. The effectiveness and robustness of the proposed scheme are assessed by comparing numerical results with analytical and reference solutions of a set of test cases. Moreover, to show the capability of the numerical model on field-scale applications, the results of a dam-break scenario are presented.     

Ekman motion in shallow open sea in the presence of time-harmonic variation of water depth

The time variation of the wind-induced flow in a homogeneous unbounded sea region has been analytically investigated. The time-dependent Ekman solution in a homogeneous, shallow open sea has been further extended by taking into account the time variation in water depth which might be caused by either tidal motion or change in the mixed layer thickness. The solution approach taken in this study is based on the Galerkin-eigenfunction method originally developed by Heaps [1972. On the numerical solution of three-dimensional hydrodynamical equations for tides and storm surges. Memoires de la Societe Royale des Sciences de Liege Serie 6(2), 143–180]. A series of calculations have been made with emphasis on the influence of the time variation in water depth upon the build-up of Ekman spirals in the presence of oscillatory variations in water depth. It has been found that two oscillations contribute to the wind drift current, the inertial oscillation and the depth-variation-induced oscillation; the inertial oscillation decays with time, but the depth-variation-induced oscillation remains undamped despite the presence of bottom friction. The presence of time-harmonic variation produces peculiar forms of hodograph with a curled or circular pattern according to the angular frequency of the water depth variation.     

Identification of wind stress on shallow water surfaces by optimal smoothing

Using the state space approach, an on-line filter procedure for combined wind stress identification and tidal flow forecasting is developed. The stochastic dynamic approach is based on the linear twodimensional shallow water equations. Using a finite difference scheme, a system representation of the model is obtained. To account for uncertainties, the system is embedded into a stochastic environment. By employing a Kalman filter, the on-line measurements of the water-level available can be used to identify and predict the shallow water flow. Because it takes a certain time before a fluctuation in the wind stress can be noticed in the water-level measurements, an optimal fixed-lag smoother is used to identify the stress.     

基于改善水质的浅水湖泊引调水模式的评价指标

以浅水湖泊引清调水的目的为基础,结合生态水力学要求,兼考虑经济性因子,尝试性地建立了浅水湖泊调水方案评价指标体系.评价指标包括湖泊水质改善效果指标、水力条件指标和经济效益指标.以玄武湖调水模式为例,采用二维水质水量模型对玄武湖不同的引水规模、引水方式、引水口和出水口流量分配的引调水方案进行模拟计算,统计出各个评价指标的结果,综合分析评价出相对最优的引调水方案.该评价方法可推广到其它类似的浅水湖泊引清调水模式的评价中.     

Experimental investigation of flow pattern and sediment deposition in rectangular shallow reservoirs

This paper reports the experimental investigation of flow pattern, preferential regions of deposition and trap efficiency as a function of the length of rectangular shallow reservoirs. Four flow patterns were identified (from longer to shorter reservoirs): an asymmetric flow with two reattachment points, an asymmetric flow with one reattachment point, an unstable flow, and a symmetric flow without any reattachment point. Using dye visualizations, the median value and the temporal variability of the reattachment lengths were precisely measured for the asymmetric flows. For each stable flow, sediment tests with plastic particles were carried out. The regions of deposition on the bed of the reservoir were clearly a function of the flow pattern. The transition from an asymmetric flow pattern to a symmetric flow pattern was responsible for an abrupt decrease of the trap efficiency; a number of regression laws were discussed to take it into account.     

Well-balanced high-order centered schemes on unstructured meshes for shallow water equations with fixed and mobile bed

In this paper, we study the numerical approximation of the two-dimensional morphodynamic model governed by the shallow water equations and bed-load transport following a coupled solution strategy. The resulting system of governing equations contains non-conservative products and it is solved simultaneously within each time step. The numerical solution is obtained using a new high-order accurate centered scheme of the finite volume type on unstructured meshes, which is an extension of the one-dimensional PRICE-C scheme recently proposed in Canestrelli et al. (2009) [5]. The resulting first-order accurate centered method is then extended to high order of accuracy in space via a high order WENO reconstruction technique and in time via a local continuous space–time Galerkin predictor method. The scheme is applied to the shallow water equations and the well-balanced properties of the method are investigated. Finally, we apply the new scheme to different test cases with both fixed and movable bed. An attractive future of the proposed method is that it is particularly suitable for engineering applications since it allows practitioners to adopt the most suitable sediment transport formula which better fits the field data.     

Application of third generation shallow water wave models in a tidal environment

Wave modeling was performed in the German Bight of the North Sea during November 2002, using the spectral wave models, namely the K-model and Simulating WAves Nearshore (SWAN), both developed for applications in environments of shallow water depths. These models mainly differ with respect to their dissipation source term expressions and in exclusion or inclusion of nonlinear wave–wave interactions. The K-model uses nonlinear dissipation and bottom dissipation, and neglects quadruplet wave–wave interaction whereas, SWAN includes, besides bottom dissipation, dissipation by white-capping and depth induced wave breaking and triad wave–wave interaction. The boundary spectra were extracted from the WAM model results of a North Sea hindcast of the HIPOCAS project, wind fields, tidal current and water level variations from the results of models used in the Belawatt project. The purpose of this study was to test the performance of both wave models to see whether they were able to predict near-shore wave conditions accurately. The runs were performed with and without tidal current and level variations to determine their effect on the waves. Comparisons of model results with buoy measurements show that taking into account tides and currents improve the spectral shape especially in areas of high current speeds. Whereas SWAN performed better in terms of spectral shape, especially in case of two peaked spectra, the K-model showed better results in terms of integrated parameters.Responsible Editor: Hans Burchard     

Lattice Boltzmann simulations of the transient shallow water flows

A two-dimensional lattice Boltzmann model (LBM) is presented for transient shallow water flows. The model is based on the shallow water equations coupled with the large eddy simulation model. In order to obtain accurate results efficiently, a multi-block lattice scheme is applied at the area where a local finer grid is needed for strong change in physical variables. The model is verified by applying to five cases with transient processes: (a) a tidal wave over steps; (b) a perturbation over a submerged hump; (c) partial dam break flow; (d) circular dam break flow; (e) interaction between a dam break surge and four square cylinders. The objectives of this study are to validate the two-dimensional LBM in transient flow simulation and provide the detailed transient processes in shallow water flows.