伪谱法弹性波场数值模拟中的边界条件

边界条件问题长期以来一直是困扰地震波数值模拟研究者的一个难题,许多人都提出了有效的方法,但是该问题仍然需要研究.本文针对伪谱法地震波场数值模拟的特殊要求,即所有网格点相互耦合,且傅立叶变换的周期性使得方程在网格边界上很难得到一个吸收的替代方程,研究衰减边界条件中衰减系数的变化对边界反射波衰减效果的影响.数值实验结果说明,随着衰减系数的增大,振幅衰减加快,用很少的过渡带网格就可以使边界上的反射波能量很小.但是如果衰减率过大,就会在传播区和过渡带产生干扰波场.因此,在衰减系数的选择上应采取折中办法,即在取某一较小的衰减系数的前提下,设置较少的过渡带网格数使边界的能量反射最少.     

An explicit method for numerical simulation of wave equations

In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with M being a positive integer for regular grids are discussed and illustrated by constructing the second order (M = 1) and the fourth order (M = 2) recursion formulas.     

浅层地震勘探中有限元数值模拟方法的应用研究

针对目前广泛用于解决工程地质问题的浅层地震勘探特点,通过对浅层小尺度介质进行地震波数值模拟,揭示地震波在不同浅层介质模型中的波场特征和传播规律.为浅层地震勘探野外数据采集前合理选择观测系统提供指导性的方案,并能对地震勘探成果解释的合理性给予客观评价.在分析各种数值模拟方法特点的基础上,选择具有精度高、易模拟复杂结构的有...     

Theoretical/numerical model for the transport of non-uniform suspended sediment in open channels

In this paper, we address the transport of multi-disperse suspended sediment mixtures in open channels, via the use of the two-fluid model. To that end, we extend previously developed frameworks for the dilute and non-dilute transport of suspended sediment. Within the scope of the Reynolds-averaged Navier-Stokes (RANS) equations, these modeling frameworks comprise mass and momentum equations for both phases (water and sediment). Here, we calculate the distribution of total volumetric concentration of sediment using two approaches: (1) by considering the mixture as represented by a single size; we call this approach Partial two-fluid model for uniform sediments (PTFMU); and (2) by combining the volumetric concentration of the sediment corresponding to several particle size classes; we call this approach Partial two-fluid model for non-uniform sediments (PTFMNU). In the second approach, we propose a methodology for the computation of the overall velocity of the disperse phase as a function of the velocities of each size class. k-ε type closures to account for the turbulence in the carrier phase (water) are applied. We also consider the coupling between the two phases through the drag force. Velocities of the carrier and disperse phases, and concentrations for each sediment class size are numerically solved by integrating the differential equations over control volumes. In order to validate our models, we compare numerical results to experimental data of Einstein and Chien [H.A. Einstein, N. Chien, Effects of heavy sediment concentration near the bed on velocity and sediment distribution, MRD sediment series report, University of California, Berkley, 1955] and Taggart et al. [W.C. Taggart, C.A. Yermoli, S. Montes, A. Ippen, Effects of sediment size and gradation on concentration profiles for turbulent flow, Massachusetts Institute of Technology, 1972]. Results of mean velocity of the carrier phase are in close agreement with the experimental data. For the prediction of sediment concentrations, we observe that there is a difference in the results using the two approaches mentioned above. We additionally obtain values of the Schmidt number needed to improve the agreement between predictions of the distribution of suspended sediment and the experimental data, and discuss the effect of sediment size and increasing sediment concentration on the values of the Schmidt number.     

三维频率域随钻电磁波测井数值模拟

在现有的普通计算机的内存和运算速度的条件下,三维频率域电磁波测井数值模拟非常困难,为了研究随钻仪器在复杂测井环境中的响应规律,从Maxwell电磁响应方程出发,针对大斜度井井眼和侵入剖面的几何特点,采用新的网格划分方法,并应用基于交错网格的有限差分得到了三维频率域电磁响应差分计算格式,采用改进的ICCG（不完全乔尔斯基共轭梯度）方法,对一维变带宽存储的大型复稀疏矩阵进行了求解,得到了随钻电磁波测井响应.计算结果表明：一维变带宽存储方法很好地解决了大型稀疏矩阵的存储问题;改进的ICCG方法得到的结果真实可靠并大大提高了求解效率;随钻电阻率曲线随着井斜角度增大会出现“极化角”,低阻侵入会使“极化角”弱化;相位电阻率曲线受围岩影响较小,受侵入影响较大.     

Large-scale numerical simulation of groundwater flow and solute transport in discretely-fractured crystalline bedrock

A large-scale fluid flow and solute transport model was developed for the crystalline bedrock at Olkiluoto Island, Finland, which is considered as potential deep geological repository for spent nuclear fuel. Site characterization showed that the main flow pathways in the low-permeability crystalline bedrock on the island are 13 subhorizontal fracture zones. Compared to other sites investigated in the context of deep disposal of spent nuclear fuel, most deep boreholes drilled at Olkiluoto are not packed-off but are instead left open. These open boreholes intersect the main fracture zones and create hydraulic connections between them, thus modifying groundwater flow. The combined impact of fracture zones and open boreholes on groundwater flow is simulated at the scale of the island. The modeling approach couples a geomodel that represents the fracture zones and boreholes with a numerical model that simulates fluid flow and solute transport. The geometry of the fracture zones that are intersected by boreholes is complex, and the 3D geomodel was therefore constructed with a tetrahedral mesh. The geomodel was imported into the numerical model to simulate a pumping test conducted on Olkiluoto Island. The pumping test simulation demonstrates that fracture-borehole intersections must be accurately discretized, because they strongly control groundwater flow. The tetrahedral mesh provides an accurate representation of these intersections. The calibrated flow model was then used for illustrative scenarios of radionuclide migration to show the impact of fracture zones on solute transport once the boreholes were backfilled. These mass transport simulations constitute base cases for future predictive analyses and sensitivity studies, since they represent key processes to take into consideration for repository performance assessment.     

Numerical simulation of elastic wave propagation based on the transversely isotropic BISQ equation

Introduction More real models are being developed by the modern seismology. As we all know, the earth is not a simple elastic body. Oil and gas reservoir, ground surface, seashore zone, sea bottom layer, etc, are porous solid media with fluids. It has been confirmed that there are two main fluid flow mechanisms in these media (Dvorkin, Nur, 1993), i.e., the Biot flow mechanism (Biot, 1956, 1962) based on the macroscopic property and the Squirt-flow mechanism (Mavko, Nur, 1979) based on the …     

Convergence of realization-based statistics to model-based statistics for the LU unconditional simulation algorithm: some numerical tests

 In geostatistics, stochastic simulation is often used either as an improved interpolation algorithm or as a measure of the spatial uncertainty. Hence, it is crucial to assess how fast realization-based statistics converge towards model-based statistics (i.e. histogram, variogram) since in theory such a match is guaranteed only on average over a number of realizations. This can be strongly affected by the random number generator being used. Moreover, the usual assumption of independence among simulated realizations of a random process may be affected by the random number generator used. Simulation results, obtained by using three different random number generators implemented in Geostatistical Software Library (GSLib), are compared. Some practical aspects are pointed out and some suggestions are given to users of the unconditional LU simulation method.     

代县地电阻率降雨影响及数值模拟

多年观测数据表明,代县地电阻率年变受降雨影响比较显著,通过褶积滤波法获取降雨量对视电阻率的影响值,定量分析降雨与视电阻率的相关性。在此基础上,利用代县地电阻率电测深数据及相关地质资料,建立三维有限元模型,模拟降雨对地电阻率的影响,进一步确定降雨对视电阻率影响的物理机制。结果表明:降雨是造成代县地电阻率年变的主要因素,且与视电阻率呈正相关性;降雨造成表层介质含水饱和度发生变化,使得相应电阻率下降近10倍,从而引起地电阻率年变幅度发生近1/10的改变。     

东海海域潜在地震海啸的数值模拟初步研究

针对东海冲绳海槽地区的地震地质背景,对东海海域潜在的地震海啸进行了预研究. 假设了冲绳海槽在发生8.5级大地震,断层错距高达15 m的极端地震情况引发的海啸对中国东部沿海地区的影响. 初步数值模拟结果表明,该地震引发的海啸的最大初始波高为4.3 m,4小时左右传至浙江沿岸,近岸各处波高为1——2 m,其中局部地区波高为2.4 m;约7——8小时靠近上海海岸线（若震源在中冲绳海槽地区,海啸传到上海最快大约7小时）,近岸波高约为1 m. 近岸区域地形变化复杂,海岛密布,局部地形条件可能会很大地影响实际各地点海啸波高,加上海啸在岸边爬高及港湾效应,估计波高还会升高. 给出了冲绳海槽南、中、北部发生潜在地震海啸的传播等时图. 笔者在东海设置了3个地震及海啸监测站,基于海啸模拟结果绘制了监测站处的海啸随时间演化曲线,分析了预研究成果对海啸预警可能发挥的作用.      

物理地震危险性分析中的地震动模拟方法研究进展

随着计算机计算能力的大幅提升和地震科学方法认知的不断丰富和完善,众多的学者开始用基于物理机理真实断层模型和速度结构模型的地震动数值模拟来评估未来可能发生的破坏性地震的灾害特征。本文初步汇总了地震动数值模拟在地震危险性分析中的研究现状、计算方法和研究实验项目,阐述了地震动数值模拟在物理地震危险性分析及应用中的发展趋势。从抗震设防的角度上来看：基于地震动数值模拟的地震危险性分析将有助于提前做好防震救灾预警。本文综述了地震动数值模拟在物理地震危险性分析中的研究和应用,该方法的不断完善会使其在未来具有巨大的潜力。     

完全匹配层在数值模拟孔隙介质声波传播中的应用

The nonsplitting perfectly matched layer （NPML） absorbing boundary condition （ABC） was first provided by Wang and Tang （2003） for the finite-difference simulation of elastic wave propagation in solids. In this paper, the method is developed to extend the NPML to simulating elastic wave propagation in poroelastic media. Biot＇s equations are discretized and approximated to a staggered-grid by applying a fourth-order accurate central difference in space and a second-order accurate central difference in time. A cylindrical twolayer seismic model and a borehole model are chosen to validate the effectiveness of the NPML. The results show that the numerical solutions agree well with the solutions of the discrete wavenumber （DW） method.     

Finite element equations and numerical simulation of elastic wave propagation in two—phase anisotropic media

Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotopic media are derived in this paper.Numerical solution of finite element equations is given.Finally,properties of elastic wave propagation are observed and analyzed through FEM modeling.     

脉动风场模拟技术的研究与进展

通过对脉动风场模拟技术近20年来国内外研究成果的总结,阐述了脉动风场模拟研究的发展历程。首先简要介绍了近地风场的特性和脉动风场模拟技术的研究现状,然后从脉动风场模拟技术的研究方法出发,分别介绍了谐波合成法、线性滤波法和小波方法在风场模拟技术中的研究发展现状和所面临的问题,最后对风场模拟技术的未来趋势进行了展望。     

Using the extended finite element method for simulation of transient well leakage in multilayer aquifers

The extended finite element (XFEM) is applied to the problem of transient leakage from abandoned or free-flowing artesian wells in perforated aquifer-aquitard systems. To more accurately capture the singularities in potentiometric head at the wells, the standard linear finite element basis is locally augmented with asymptotic analytical solutions which enable more accurate calculations of leakage rates between aquifers. Highly accurate flux estimates are obtained without the need for higher mesh resolution near wells. Simulations are carried out to test both the accuracy and convergence properties of the XFEM implementation, and the XFEM results are compared to those of a high-resolution standard finite element model. It is seen that for the type of singularity-driven problem posed here, the standard FEM is unable to resolve leakage rates without very fine discretization, but that the XFEM performs robustly with fewer degrees of freedom. The impact of aquifer geometric heterogeneity on leakage rates is assessed and seen to be an important factor in determining total leakage. It is demonstrated that the XFEM may be a valuable tool in many water resources applications where small-scale effects can impact global system behavior.     

活断层上均匀盖层中氡浓度分布的数值模拟及反演拟合

Based on the convection and diffusion mechanisms of radon migration, in this paper we deduce the two-dimensional differential equation for radon transportation in the overburden above active fault zones with an unlimited extension along the strike. Making use of the finite difference method, the radon concentration distribution in the overburden above active faults is calculated and modeled. The active fault zone parameters, such as the depth and the width of the fault zone, and the value of radon concentration, can be inverted from the measured radon concentration curve. These realize quantitative interpretation for radon concentration anomalies. The inversion results are in good agreement with the actual fault zone parameters.     

基于柱坐标系的隧道空间溶洞探测全波场数值模拟

溶洞作为碳酸盐岩地区常见的一类局部不良地质体,其波场特征常常与形态、填充物、空间方位等因素相关,而受此影响地震类方法对溶洞的探测效果往往不理想.因此,要想实现地震类方法对溶洞的准确探测,首先需要明确不同状态下溶洞的波场特征.针对上述问题,考虑到实际隧道空间形态,本文首次采用柱坐标系下的一阶速度-应力弹性波动方程对隧道空间进行数值模拟,通过通量矫正技术压制频散,同时对已开挖段隧道与围岩的边界采用自由边界条件处理.对不同半径、不同填充物、不同方位的溶洞以及复杂溶洞模型的数值模拟结果及波场特征分析表明:本文提出的方法可实现隧道空间全波场高精度数值模拟,波场特征符合波的运动学与动力学特征;不同情况下溶洞波场特征的差异可有效的指导溶洞的探测与识别.     

Parallel numerical simulation of two-phase flow model in porous media using distributed and shared memory architectures

A two-phase (water and oil) flow model in a homogeneous porous media is studied, considering immiscible and incompressible displacement. This model is numerically solved using the Finite Volume Method (FVM) and we compare four numerical schemes for the approximation of fluxes on the faces of the discrete volumes. We describe briefly how to obtain the mathematical and computational models applying axiomatic formulations and generic programming. Two strategies of parallelization are implemented in order to reduce the execution time. We study distributed (cluster of CPUs) and shared (Graphics Processing Units) memory architectures. A performance comparison of these two architectures is done along with an analysis of the four numerical schemes, for a water-flooding five-spot pattern model.     

上海地区地壳应变能力数值模拟

上海市综合深井项目第2阶段将建设9—10个综合深井观测站。利用蠕变方程,通过ANSYS有限元软件进行数值分析,模拟上海市9个观测点在当前地壳运动水平下的形变信息,为综合深井第2阶段建设提供理论基础。     

Mechanical model of precursory source proces for strong earthquakes in western Yun－nanarea——conceptandnumericalsimula┐tion

Ｍｅｃｈａｎｉｃａｌｍｏｄｅｌｏｆｐｒｅｃｕｒｓｏｒｙｓｏｕｒｃｅｐｒｏ┐ｃｅｓｆｏｒｓｔｒｏｎｇｅａｒｔｈｑｕａｋｅｓｉｎｗｅｓｔｅｒｎＹｕｎ┐ｎａｎａｒｅａ——ｃｏｎｃｅｐｔａｎｄｎｕｍｅｒｉｃａｌｓｉｍｕｌａ┐ｔｉｏｎＹＯＵ－ＪＩＮＳＵ（苏有...