三维泊松方程数值模拟的多重网格方法

本文简要介绍了多重网格方法的基本思想和原理,然后应用多重网格(MG)方法求解三维泊松方程,网格尺度从17×17×17逐次增加至257×257×257,并与不完全Chelesky共轭梯度法(ICCG)、Gauss直接解法进行比较,结果表明,MG方法计算速度明显优于ICCG、Gauss方法,对于129×129×129网格的三维数值模拟费时43s,比ICCG法快7倍,而对于257×257×257超大型网格的三维数值模拟也仅需412s.     

三维波动方程时空域混合网格有限差分数值模拟方法

常规高阶和时空域高阶有限差分方法广泛应用于三维标量波动方程的数值模拟,这两种差分方法仅利用笛卡尔坐标系中的坐标轴网格点构建三维Laplace差分算子,相应的差分离散波动方程本质上仅具有2阶差分精度,模拟精度低.本文将三维笛卡尔坐标系中非坐标轴网格点分为两类:坐标平面内的非坐标轴网格点和坐标平面外的非坐标轴网格点,系统推...     

三维电磁感应的有限差分法的网格收敛试验

本文用三维有限差分法,在计算电导率异常的电磁响应时作了一次网格收敛试验.该收敛解和其它几种数值解法得到的结果作了比较.这些方法是Hohmann的积分方程法、Lee等的混合法和Gupta等的混合法.我们的结果基本上和Hohmann的以及Gupta等的结果一致.试验表明,有限差分法所用的计算机时间比SANGAM混合法的少.为了缩减迭代过程的计算时间,我们使用了一种新的技巧——变区域迭代法,这种方法的收敛速度和一般的（在全区域迭代的）迭代法的收敛速度相同.     

砂土液化大变形本构模型的三维化及其数值实现

基于砂土液化大变形机理和适用于二维条件的边界面弹塑性本构模型,发展了符合三维应力空间中边界面和剪胀面上的应力映射规则,建立了三维应力空间中砂土液化大变形本构模型.针对模型特点采用半显式的Cutting Plane算法进行应力积分,并采用Pegasus求根算法根据映射规则计算边界面上的应力映射点,在OpenSees开源有限元平台上实现了三维模型的数值化.结合完全耦合的u-p格式有限元单元,对饱和砂土不排水循环扭剪试验进行了模拟,并进行了一个真三维倾斜地基的动力反应分析.计算结果表明模型和所采用的数值算法具有很好的模拟和分析三维条件下砂土液化后大变形的能力.     

三维电磁感应的有限差分法的网格收敛试验

本文用三维有限差分法,在计算电导率异常的电磁响应时作了一次网格收敛试验.该收敛解和其它几种数值解法得到的结果作了比较.这些方法是Hohmann的积分方程法、Lee等的混合法和Gupta等的混合法.我们的结果基本上和Hohmann的以及Gupta等的结果一致.试验表明,有限差分法所用的计算机时间比SANGAM混合法的少.为了缩减迭代过程的计算时间,我们使用了一种新的技巧--变区域迭代法,这种方法的收敛速度和一般的（在全区域迭代的）迭代法的收敛速度相同.     

交错网格高阶差分法三维弹性波数值模拟

本文应用交错网格高阶有限差分方法模拟弹性波在三维各向同性介质中的传播。采用时间上二阶、空间上高阶近似的交错网格高阶差分公式求解三维弹性波位移-应力方程,并在计算边界处应用基于傍轴近似法得到的三维弹性波方程吸收边界条件。在此基础上进行了三维盐丘地质模型的地震波传播数值模拟试算。试算结果表明该方法模拟精度高,在很大程度上减小了数值频散,绕射波更加丰富,而且适用于介质速度具有纵向变化和横向变化的情况。     

基于精细积分法的三维弹性波数值模拟(英文)

波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。     

基于非规则双重网格的三维声波方程模拟

三维声波方程相比二维声波方程能够更好的模拟三维空间的地震波传播,模拟标量近似下的弹性波在三维复杂介质的传播过程.基于非规则网格的正演模拟方法的格子法可以处理很好的刻画起伏地表、速度间断面等复杂构造,但是这类方法需要大量的几何描述来描述网格.本文提出了三维六面体双重网格的格子法来模拟声波方程,一方面该方法继承了格子法能够灵活处理自由表面和速度间断面的特性.另一方面,该方法通过双重网格的实现极大的减少了几何描述文件的大小,可以最大的实现GPU加速,实现粗粒度并行,在节省了几何描述空间的同时达到了很高的加速比.

     

渗透率场敏感系数的数值计算

由渗流微分方程定解问题，利用格林互易定理从理论上导出了渗透率场敏感系数 的计算公式，并借助于数值积分和差分方法给出了渗透率场敏感系数的离散形式. 由Peacem an方程建立了井压渗透率场敏感系数与网格压力渗透率场敏感系数的关系. 对理论模型利用 三维不均匀非稳定渗流场的压强数值解计算了井压对渗透率场的敏感系数，并用敏感系数的 直接计算方法进行了验证，结果表明本文的计算方法与直接计算方法相吻合.     

基于非结构网格剖分的海洋可控源电磁三维正则化反演研究

为更好地处理与解释复杂海底地形条件下测得的海洋可控源电磁数据,本文提出了一种基于非结构网格剖分的频率域海洋可控源电磁数据三维正则化反演方法.该方法首先对海洋地电模型以非结构四面体单元进行离散,然后基于矢量有限元方法获得海洋可控源电磁响应和灵敏度信息,最后采用共轭梯度法求解高斯-牛顿反演方程计算模型修正量.为提高反演的稳定性,通过在反演过程中采用对数转换方法实现反演模型参数的上下限约束.本文分别测试了单测线水平海底地形反演算例和面积性测量的起伏海底地形反演算例.反演结果表明,本文提出的频率域海洋可控源电磁三维反演能够准确地恢复高阻储油层的位置和电阻率信息,且计算效率较高,可用于实测海洋电磁资料的处理与解释.

     

A meshless method to simulate solute transport in heterogeneous porous media

We derive a meshless numerical method based on smoothed particle hydrodynamics (SPH) for the simulation of conservative solute transport in heterogeneous geological formations. We demonstrate that the new proposed scheme is stable, accurate, and conserves global mass. We evaluate the performance of the proposed method versus other popular numerical methods for the simulation of one- and two-dimensional dispersion and two-dimensional advective–dispersive solute transport in heterogeneous porous media under different Pèclet numbers. The results of those benchmarks demonstrate that the proposed scheme has important advantages over other standard methods because of its natural ability to control numerical dispersion and other numerical artifacts. More importantly, while the numerical dispersion affecting traditional numerical methods creates artificial mixing and dilution, the new scheme provides numerical solutions that are “physically correct”, greatly reducing these artifacts.     

Coupling upscaling finite element method for consolidation analysis of heterogeneous saturated porous media

The coupling upscaling finite element method is developed for solving the coupling problems of deformation and consolidation of heterogeneous saturated porous media under external loading conditions. The method couples two kinds of fully developed methodologies together, i.e., the numerical techniques developed for calculating the apparent and effective physical properties of the heterogeneous media and the upscaling techniques developed for simulating the fluid flow and mass transport properties in heterogeneous porous media. Equivalent permeability tensors and equivalent elastic modulus tensors are calculated for every coarse grid block in the coarse-scale model of the heterogeneous saturated porous media. Moreover, an oversampling technique is introduced to improve the calculation accuracy of the equivalent elastic modulus tensors. A numerical integration process is performed over the fine mesh within every coarse grid element to capture the small scale information induced by non-uniform scalar field properties such as density, compressibility, etc. Numerical experiments are carried out to examine the accuracy of the developed method. It shows that the numerical results obtained by the coupling upscaling finite element method on the coarse-scale models fit fairly well with the reference solutions obtained by traditional finite element method on the fine-scale models. Moreover, this method gets more accurate coarse-scale results than the previously developed coupling multiscale finite element method for solving this kind of coupling problems though it cannot recover the fine-scale solutions. At the same time, the method developed reduces dramatically the computing effort in both CPU time and memory for solving the transient problems, and therefore more large and computational-demanding coupling problems can be solved by computers.     

Numerical simulations of non-ergodic solute transport in three-dimensional heterogeneous porous media

Numerical simulations of non-ergodic transport of a non-reactive solute plume by steady-state groundwater flow under a uniform mean velocity, , were conducted in a three-dimensional heterogeneous and statistically isotropic aquifer. The hydraulic conductivity, K(x), is modeled as a random field which is assumed to be log-normally distributed with an exponential covariance. Significant efforts are made to reduce the simulation uncertainties. Ensemble averages of the second spatial moments of the plume and the plume centroid variances were simulated with 1600 Monte Carlo (MC) runs for three variances of log K, _Y²=0.09, 0.23, and 0.46, and a square source normal to of three dimensionless lengths. It is showed that 1600 MC runs are needed to obtain stabilized results in mildly heterogeneous aquifers of _Y²0.5 and that large uncertainty may exist in the simulated results if less MC runs are used, especially for the transverse second spatial moments and the plume centroid variance in transverse directions. The simulated longitudinal second spatial moment and the plume centroid variance in longitudinal direction fit well to the first-order theoretical results while the simulated transverse moments are generally larger than the first-order values. The ergodic condition for the second spatial moments is far from reaching in all cases simulated and transport in transverse directions may reach ergodic condition much slower than that in longitudinal direction.     

Estimates of effective permeability for non-Newtonian fluid flow in randomly heterogeneous porous media

An understanding of the interplay between non-Newtonian effects in porous media flow and field-scale domain heterogeneity is of great importance in several engineering and geological applications. Here we present a simplified approach to the derivation of an effective permeability for flow of a purely viscous power–law fluid with flow behavior index n in a randomly heterogeneous porous domain subject to a uniform pressure gradient. A standard form of the flow law generalizing the Darcy’s law to non-Newtonian fluids is adopted, with the permeability coefficient being the only source of randomness. The natural logarithm of the permeability is considered a spatially homogeneous and correlated Gaussian random field. Under the ergodic hypothesis, an effective permeability is first derived for two limit 1-D flow geometries: flow parallel to permeability variation (serial-type layers), and flow transverse to permeability variation (parallel-type layers). The effective permeability of a 2-D or 3-D isotropic domain is conjectured to be a power average of 1-D results, generalizing results valid for Newtonian fluids under the validity of Darcy’s law; the conjecture is validated comparing our results with previous literature findings. The conjecture is then extended, allowing the exponents of the power averaging to be functions of the flow behavior index. For Newtonian flow, novel expressions for the effective permeability reduce to those derived in the past. The effective permeability is shown to be a function of flow dimensionality, domain heterogeneity, and flow behavior index. The impact of heterogeneity is significant, especially for shear-thinning fluids with a low flow behavior index, which tend to exhibit channeling behavior.     

Modelling of groundwater flow in heterogeneous porous media by finite element method

The HySuf‐FEM code (Hydrodynamic of Subsurface Flow by Finite Element Method) is proposed in this article in order to estimate the spatial variability of the transmissivity values of the Berrechid aquifer (Morocco). The calibration of the model is based on the hydraulic head, hydraulic conductivity and recharge. Three numerical tests are used to validate the model and verify its convergence. The first test case consists in using the steady analytical solution of the Poisson equation. In the second, the model has been compared with the hydrogeological system which is characterized by an unconfined monolayer (isotropic layer) and computed by using PMWIN‐MODFLOW software. The third test case is based on the comparison between the results of HySuf‐FEM and the multiple cell balance method in the aquifer system with natural boundaries case. Good agreement between the Hydrodynamic of Subsurface Flow, the numerical tests and the spatial distribution of the thickening of the hydrogeological system is deduced from the analysis and the interpretations of hydrogeological wells. Copyright © 2011 John Wiley & Sons, Ltd.     

基于BISQ模型的三维双相各向异性介质数值模拟

Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.     

Probabilistic collocation and lagrangian sampling for advective tracer transport in randomly heterogeneous porous media

The Karhunen-Loeve (KL) decomposition and the polynomial chaos (PC) expansion are elegant and efficient tools for uncertainty propagation in porous media. Over recent years, KL/PC-based frameworks have successfully been applied in several contributions for the flow problem in the subsurface context. It was also shown, however, that the accurate solution of the transport problem with KL/PC techniques is more challenging. We propose a framework that utilizes KL/PC in combination with sparse Smolyak quadrature for the flow problem only. In a subsequent step, a Lagrangian sampling technique is used for transport. The flow field samples are calculated based on a PC expansion derived from the solutions at relatively few quadrature points. To increase the computational efficiency of the PC-based flow field sampling, a new reduction method is applied. For advection dominated transport scenarios, where a Lagrangian approach is applicable, the proposed PC/Monte Carlo method (PCMCM) is very efficient and avoids accuracy problems that arise when applying KL/PC techniques to both flow and transport. The applicability of PCMCM is demonstrated for transport simulations in multivariate Gaussian log-conductivity fields that are unconditional and conditional on conductivity measurements.     

Effect of permeability variations on solute transport in highly heterogeneous porous media

The effect of aquifer heterogeneity on flow and solute transport in two-dimensional isotropic porous media was analyzed using the Monte Carlo method. The two-dimensional logarithmic permeability (ln K) was assumed to be a non-stationary random field with its increments being a truncated fractional Lévy motion (fLm). The permeability fields were generated using the modified successive random additions (SRA) algorithm code SRA3DC [http://www.iamg.org/CGEditor/index.htm]. The velocity and concentration fields were computed respectively for two-dimensional flow and transport with a pulse input using the finite difference codes of MODFLOW 2000 and MT3DMS. Two fLm control parameters, namely the width parameter (C) and the Lévy index (α), were varied systematically to examine their effect on the resulting permeability, flow velocity and concentration fields. We also computed the first- and second-spatial moments, the dilution index, as well as the breakthrough curves at different control planes with the corresponding concentration fields. In addition, the derived breakthrough curves were fitted using the continuous time random walk (CTRW) and the traditional advection-dispersion equation (ADE). Results indicated that larger C and smaller α both led to more heterogeneous permeability and velocity fields. The Lévy-stable distribution of increments in ln K resulted in a Lévy-stable distribution of increments in logarithm of the velocity (ln v). Both larger C and smaller α created sharper leading edges and wider tailing edges of solute plumes. Furthermore, a relatively larger amount of solute still remained in the domain after a relatively longer time transport for smaller α values. The dilution indices were smaller than unity and increased as C increased and α decreased. The solute plume and its second-spatial moments increased as C increased and α decreased, while the first-spatial moments of the solute plume were independent of C and α values. The longitudinal macrodispersivity was scale-dependent and increased as a power law function of time. Increasing C and decreasing α both resulted in an increase in longitudinal macrodispersivity. The transport in such highly heterogeneous media was slightly non-Gaussian with its derived breakthrough curves being slightly better fitted by the CTRW than the ADE, especially in the early arrivals and late-time tails.