一维波动方程波阻抗反演的同伦方法

文中从地震勘探一维波动方程反问题出发,研究了一种反演地层参数的同伦方法,该方法把非线性方程组的求解转化成常微分方程初值问题的数值求解,从而给出一种稳定的计算速度快、抗噪能力强的全局收敛的反演方法．理论模型和实例试算的结果表明了同伦方法是一种有效的反演算法,特别适用于非线性、多极值的地球物理反演问题,在地球物理非线性反演中具有广泛的应用前景．     

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.     

多孔弹性波方程的多尺度波场模拟

本文研究了二维多孔弹性波方程的多尺度波场数值模拟方法.该多尺度方法可采用较粗的网格计算,同时又能反映细尺度上物性参数的变化信息.文中详细阐述了多尺度模拟方法与算法,并推导了相应的计算格式.基本思想是建立粗细两套网格,在粗网格上,基于有限体积方法计算更新波场;在细网格上,计算多尺度基函数,这基于有限元方法通过求解一个局部化问题得到.对含有随机分布散射体的多孔介质模型进行了数值计算,计算中应用了完全匹配层(PML)吸收边界条件,数值结果验证了本文方法和算法的正确性和有效性.     

Parallel Algorithms for Computational Models of Geophysical Systems

Mathematical models of many geophysical systems are based on the computational processing of large-scale algebraic systems. The most advanced computational tools are based on massively parallel processors. The most effective software for solving partial differential equations in parallel intends to achieve the DDM-paradigm. A set of four algorithms, the DVS-algorithms, which achieve it, and of very general applicability, has recently been developed and here they are explained. Also, their application to problems that frequently occur in Geophysics is illustrated.     

Stochastic partial differential equations in groundwater hydrology

Many problems in hydraulics and hydrology are described by linear, time dependent partial differential equations, linearity being, of course, an assumption based on necessity.Solutions to such equations have been obtained in the past based purely on deterministic consideration. The derivation of such a solution requires that the initial conditions, the boundary conditions, and the parameters contained within the equations be stipulated in exact terms. It is obvious that the solution so derived is a function of these specified, values.There are at least four ways in which randomness enters the problem. i) the random initial value problem; ii) the random boundary value problem; iii) the random forcing problem when the non-homogeneous part becomes random and iv) the random parameter problem.Such randomness is inherent in the environment surrounding the system, the environment being endowed with a large number of degrees of freedom.This paper considers the problem of groundwater flow in a phreatic aquifer fed by rainfall. The goveming equations are linear second order partial differential equations. Explicit form solutions to this randomly forced equation have been derived in well defined regular boundaries. The paper also provides a derivation of low order moment equations. It contains a discussion on the parameter estimation problem for stochastic partial differential equations.     

Coarsening of three-dimensional structured and unstructured grids for subsurface flow

We present a generic, semi-automated algorithm for generating non-uniform coarse grids for modeling subsurface flow. The method is applicable to arbitrary grids and does not impose smoothness constraints on the coarse grid. One therefore avoids conventional smoothing procedures that are commonly used to ensure that the grids obtained with standard coarsening procedures are not too rough. The coarsening algorithm is very simple and essentially involves only two parameters that specify the level of coarsening. Consequently the algorithm allows the user to specify the simulation grid dynamically to fit available computer resources, and, e.g., use the original geomodel as input for flow simulations. This is of great importance since coarse grid-generation is normally the most time-consuming part of an upscaling phase, and therefore the main obstacle that has prevented simulation workflows with user-defined resolution. We apply the coarsening algorithm to a series of two-phase flow problems on both structured (Cartesian) and unstructured grids. The numerical results demonstrate that one consistently obtains significantly more accurate results using the proposed non-uniform coarsening strategy than with corresponding uniform coarse grids with roughly the same number of cells.     

On a monotonicity preserving Eulerian–Lagrangian localized adjoint method for advection–diffusion equations

Eulerian–Lagrangian localized adjoint methods (ELLAMs) provide a general approach to the solution of advection-dominated advection–diffusion equations allowing large time steps while maintaining good accuracy. Moreover, the methods can treat systematically any type of boundary condition and are mass conservative. However, all ELLAMs developed so far suffer from non-physical oscillations and are usually implemented on structured grids. In this paper, we propose a finite volume ELLAM which incorporates a novel correction step rendering the method monotone while maintaining conservation of mass. The method has been implemented on fully unstructured meshes in two space dimensions. Numerical results demonstrate the applicability of the method for problems with highly non-uniform flow fields arising from heterogeneous porous media.     

The method of calculation of the geophysical fields in piecewise homogeneous media

Recurrent algorithms for solving 3D direct problems on the distribution of geophysical fields in the piecewise-homogeneous media are built, based on methods of integral transforms, integral representations, and integral equations. The algorithms allow the recalculation of Green’s functions of one medium into Green’s function of another medium. The computational procedure is applicable to problems with equations of elliptic, parabolic, hyperbolic, and mixed types. The algorithms can be split and implemented on computing clusters and multiprocessor complexes. Inverse problems and algorithms of searching for boundaries of piecewise-homogeneous media as extremals of the Tikhonov regulating functional, approximated by the spline functions, are formulated.     

A mixed-dimensional finite volume method for two-phase flow in fractured porous media

We present a vertex-centered finite volume method for the fully coupled, fully implicit discretization of two-phase flow in fractured porous media. Fractures are discretely modeled as lower dimensional elements. The method works on unstructured, locally refined grids and on parallel computers with distributed memory. An implicit time discretization is employed and the nonlinear systems of equations are solved with a parallel Newton-multigrid method. Results from two-dimensional and three-dimensional simulations are presented.     

一种反演地幔横向密度差异的方法

重力位球谐系数及板块运动速度是两种独立的地球物理观测资料,但二者同时又是地幔内横向密度差异及与之相联系的地幔流动的结果,从而,中长波重力异常与板块运动的观测提供了探索地幔横向非均匀的一条途径。本文根据适合于低雷诺流动的最小粘滞损耗原理,利用张量场及矢量场的标量分解与传播子解的数学方法,将观测的重力位球谐系数和板块运动速度作为约束条件,发展了一种反演地幔横向密度差异的方法。     

Unstructured Gridding for MODFLOW from Prior Groundwater Flow Models: A New Paradigm

This work introduces a new unstructured gridding approach relying on feedback from a previous groundwater flow model. All cells in a relatively coarse model using a rectilinear grid are recursively subdivided following a cell wise specific discharge-based indicator to generate quadtree, octree or Voronoï grids. This technique leverages the full potential of the latest MODFLOW engines. The suitability of this approach is demonstrated on challenging single and multilayered heterogeneous formations. The proposed method is straightforward to implement in existing software packages. It supports iterative updating of groundwater flow models from the legacy rectilinear to unstructured grids.     

Multilayer shallow water flow using lattice Boltzmann method with high performance computing

A multilayer lattice Boltzmann (LB) model is introduced to solve three-dimensional wind-driven shallow water flow problems. The multilayer LB model avoids the expensive Navier–Stokes equations and obtains stratified horizontal flow velocities as vertical velocities are relatively small and the flow is still within the shallow water regime. A single relaxation time BGK method is used to solve each layer coupled by the vertical viscosity forcing term. To increase solution stability, an implicit step is suggested to obtain flow velocities. The main advantage of using the LBM is that after selecting appropriate equilibrium distribution functions, the LB algorithm is only slightly modified for each layer and retains all the simplicities of the LBM within the high performance computing (HPC) environment. The performance of the parallel LB model for the multilayer shallow water equations is investigated on CPU-based HPC environments using OpenMP. We found that the explicit loop control with cache optimization in LBM gives better performance on execution time, speedup and efficiency than the implicit loop control as the number of processors increases. Numerical examples are presented to verify the multilayer LB model against analytical solutions. We demonstrate the model’s capability of calculating lateral and vertical distributions of velocities for wind-driven circulation over non-uniform bathymetry.     

A mass conservative scheme for simulating shallow flows over variable topographies using unstructured grid

Most available numerical methods face problems, in the presence of variable topographies, due to the imbalance between the source and flux terms. Treatments for this problem generally work well for structured grids, but most of them are not directly applicable for unstructured grids. On the other hand, despite of their good performance for discontinuous flows, most available numerical schemes (such as HLL flux and ENO schemes) induce a high level of numerical diffusion in simulating recirculating flows. A numerical method for simulating shallow recirculating flows over a variable topography on unstructured grids is presented. This mass conservative approach can simulate different flow conditions including recirculating, transcritical and discontinuous flows over variable topographies without upwinding of source terms and with a low level of numerical diffusion. Different numerical tests cases are presented to show the performance of the scheme for some challenging problems.     

A two-grid method for coupled free flow with porous media flow

This paper presents a two-grid method for solving systems of partial differential equations modelling incompressible free flow coupled with porous media flow. This work considers both the coupled Stokes and Darcy as well as the coupled Navier-Stokes and Darcy problems. The numerical schemes proposed are based on combinations of the continuous finite element method and the discontinuous Galerkin method. Numerical errors and convergence rates for solutions obtained from the two-grid method are presented. CPU times for the two-grid algorithm are shown to be significantly less than those obtained by solving the fully coupled problem.     

Simulation of dam-break flow with grid adaptation

The unsteady free surface flow caused by sudden collapse of a dam produces discontinuities in the flow variables. As the flow surges downstream, it forms a moving bore front with steep gradients of water height and velocity. In the numerical simulation of this flow, proper grid distribution can play a crucial part in the prediction and resolution of the solutions. The use of presently available numerical schemes to solve this problem on a uniform course grid system fails to resolve the characteristic flow features and hence do a poor job in simulating this flow. In this paper, an adaptive grid which adjusts itself as the solution evolves is used for a better resolution of the flow properties. Rai and Anderson's¹² method is used to determine the grid speed; however, a different partial differential equation based on the conservative principle of grid arc lengths for clustering grids in one-dimensional flow is used along with the St. Venant equations to numerically simulate the flow. Both the subcritical and the supercritical flows under extreme boundary conditions are solved using this technique. With a specified number of grid points, this provides better quality solutions as compared to those obtained with uniformly distributed grids.     

A NEW MATHEMATICAL MODEL FOR2-D FLOW

A new mathematical model for 2-D flow is formulated with accurate satisfac tion of boundary conditions in conjunction with square or any grids,so that it may also yield accurate results when the domain of computation is in the shape of a strip as occurring in rivers. The basic equations are split into three sets of compo- nent equations of which two sets may be transformed to the same form. This model is a part or a model devised for 2-D flows with sediment.     

Practical aspects of stochastic dynamic tidal modelling

Kalman filtering for stochastic dynamic tidal models, is a hyperbolic filtering problem. The questions of observability and stability of the filter as well as the effects of the finite difference approximation on the filter performance are studied. The degradation of the performance of the filter, in case an erroneous filter model is used, is investigated. In this paper we discuss these various practical aspects of the application of Kalman filtering for tidal flow identification problems. Filters are derived on the basis of the linear shallow water equations. Analytical methods are used to study the performance of the filters under a variety of circumstances.     

非线性地球物理反演方法:回顾与展望

自20世纪90年代以来，非线性地球物理反演方法已走向成熟，这些方法包括线性化迭代法，仿真淬火法，遗传算法及联合反演方法等，线性化迭代反演基于开放物理系统状态发生相变的原理，要进一步改善模型参数化的技术和迭代过程的自组织；仿真淬火法与遗传算法基于自然过程的指数率或生物演化的优生率，可以相互结合以提高解估计的分辨率与置信度；联合反演要结合岩石物理性质的统计规律，才能取得兼容地质与综合方法的应用效果，地球动力学中的反问题不仅涉及偏微分方程系数项的求解，而且涉及初始条件或初始边界的求解，对地球动力学作用过程研究特别重要。     

Practical implementation of the fractional flow approach to multi-phase flow simulation

Fractional flow formulations of the multi-phase flow equations exhibit several attractive attributes for numerical simulations. The governing equations are a saturation equation having an advection diffusion form, for which characteristic methods are suited, and a global pressure equation whose form is elliptic. The fractional flow approach to the governing equations is compared with other approaches and the implication of equation form for numerical methods discussed. The fractional flow equations are solved with a modified method of characteristics for the saturation equation and a finite element method for the pressure equation. An iterative algorithm for determination of the general boundary conditions is implemented. Comparisons are made with a numerical method based on the two-pressure formulation of the governing equations. While the fractional flow approach is attractive for model problems, the performance of numerical methods based on these equations is relatively poor when the method is applied to general boundary conditions. We expect similar difficulties with the fractional flow approach for more general problems involving heterogenous material properties and multiple spatial dimensions.     

A Stable and Efficient Numerical Algorithm for Unconfined Aquifer Analysis

The nonlinearity of equations governing flow in unconfined aquifers poses challenges for numerical models, particularly in field-scale applications. Existing methods are often unstable, do not converge, or require extremely fine grids and small time steps. Standard modeling procedures such as automated model calibration and Monte Carlo uncertainty analysis typically require thousands of model runs. Stable and efficient model performance is essential to these analyses. We propose a new method that offers improvements in stability and efficiency and is relatively tolerant of coarse grids. It applies a strategy similar to that in the MODFLOW code to the solution of Richard's equation with a grid-dependent pressure/saturation relationship. The method imposes a contrast between horizontal and vertical permeability in gridblocks containing the water table, does not require "dry" cells to convert to inactive cells, and allows recharge to flow through relatively dry cells to the water table. We establish the accuracy of the method by comparison to an analytical solution for radial flow to a well in an unconfined aquifer with delayed yield. Using a suite of test problems, we demonstrate the efficiencies gained in speed and accuracy over two-phase simulations, and improved stability when compared to MODFLOW. The advantages for applications to transient unconfined aquifer analysis are clearly demonstrated by our examples. We also demonstrate applicability to mixed vadose zone/saturated zone applications, including transport, and find that the method shows great promise for these types of problem as well.