横向各向同性介质紧致交错网格有限差分波场模拟(英文)

针对有限差分数值模拟的频散问题,本文将交错网格技术和紧致差分格式相结合,推导了横向各向同性介质一阶速度一应力波动方程的紧致交错网格差分格式;对比分析了紧致交错网格差分格式、交错网格差分格式以及紧致差分格式的截断误差主项,并利用Fourier误差分析方法分析了上述三种差分格式的近似精度;在此基础上,分别采用上述三种差分格式进行了波场数值模拟。结果表明,当差分方程阶数相同时,紧致交错网格差分格式截断误差最小,数值频散最弱,差分精度最高,证实了该方法的有效性。     

弹性波在横向各向同性弹性土和饱和土界面上的反射与透射.

根据横向各向同性饱和土体的Biot波动方程,推导出qP1波从饱和土入射到弹性介质的反射系数和透射系数的一般算式,数值分析了各向异性,边界透水条件对反射系数和透射系数的影响．本文的推导思路,适用于qSV波、qP2波从饱和土体入射到弹性土体的情形．      

基于Biot-Squirt方程的波场模拟

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制,对地震波和声波的传播均产生重要影响. Dvorkin和Nur提出了同时包含Biot流动和喷射流动力学机制的统一的BISQ(Biot-Squirt)模型,基于这一模型,尽管有关弹性波在多孔隙介质中的衰减和频散问题已被广泛研究,然而,基于BISQ波传播方程的波场数值模拟至今仍未见报道. 本文从同时包含两种力学机制的孔隙弹性波方程出发,利用FCT有限差分法对含流体孔隙各向同性介质中的地震波和声波进行了数值模拟,并与基于Biot流动的Biot理论之模拟结果进行比较. 数值模拟结果表明:同时包含Biot流动和喷射流动影响的地震波和声波速度比仅包含Biot流动作用的地震波和声波速度慢,慢P波的衰减比根据Biot理论模拟的慢P波衰减更强.     

基于Biot-Squirt方程的波场模拟

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制,对地震波和声波的传播均产生重要影响. Dvorkin和Nur提出了同时包含Biot流动和喷射流动力学机制的统一的BISQ(Biot-Squirt)模型,基于这一模型,尽管有关弹性波在多孔隙介质中的衰减和频散问题已被广泛研究,然而,基于BISQ波传播方程的波场数值模拟至今仍未见报道. 本文从同时包含两种力学机制的孔隙弹性波方程出发,利用FCT有限差分法对含流体孔隙各向同性介质中的地震波和声波进行了数值模拟,并与基于Biot流动的Biot理论之模拟结果进行比较. 数值模拟结果表明:同时包含Biot流动和喷射流动影响的地震波和声波速度比仅包含Biot流动作用的地震波和声波速度慢,慢P波的衰减比根据Biot理论模拟的慢P波衰减更强.     

改进的一阶弹性波方程交错网格高阶差分解法

本文对董良国、马在田等人发表在<地球物理学报>上的一阶弹性波方程交错网格高阶差分解法一文进行了差分格式的改进.发现原文对一阶空间导数采用2N阶差分精度,而对高阶导数却采用简单的中心差分,即对同一方程中的各阶导数的处理存在有的采用高阶格式、有的采用低阶格式的不一致性问题,本文的改进点正是消除了这种不一致性,建立了时间和空间均为四阶精度的差分格式.该方法具有编程简便、易于计算机实现和精度较高的特点.     

基于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.     

一阶弹性波方程交错网格高阶差分解法

提高计算精度和运算效率是所有波场正演方法所追求的目标,本文通过将速度 （应力）对时间的奇数阶高阶寻数转化为应力（速度）对空间的导数,运用时间和空间差分精度 均可达任意阶的高阶差分法,通过交错网格技术,对一阶速度－应力弹性波方程进行了数值求 解．波场快照以及实际模型的正演结果表明,这种求解一阶弹性波方程的高阶差分解法,和 常规的差分法相比网格频散显著减小,精度明显提高,而且可以取较大的空间步长,提高计算 效率。     

二维横各向同性弹性随机介质中的波场特征

本文通过交错网格有限差分正演．模拟了平面地震波在二维横各向同性弹性随机介质模型中的传播及其自激自收时间记录．为研究横各向同性弹性随机介质模型中的波场特征,我们在五个不同的时间区段上,分别计算剖面的三个统计特征(横向中心频率、纵向中心频率、波场能量相对值)．这样,对应每一个横各向同性弹性随机介质模型．均可计算得到15个不同的波场特征量．我们通过在二维横各向同性弹性随机介质中的正演模拟．研究当自相关长度以及介质的各向异性系数变化时,对应的上述波场特征量的变化特点．证实了在随机介质模型中．各向异性系数的变化会引起波场记录上的某些统计特征的变化,归纳得出了若干结论．     

二维横向各向同性介质中人为边界反射的吸收——差分法弹性波场模拟

本文对有限差分法弹性波场三分量正演模拟中遇到的人为边界反射问题进行了研究,分别获得2-D横向各向同性介质中U_x,U_z,与U_y位移分量模拟时的边界吸收条件.从吸收边界上反射系数与入射角之间的关系可以说明,本文提出的边界吸收条件的吸收能力比较强.     

二维横向各向同性介质中人为边界反射的吸收——差分法弹性波场模拟

本文对有限差分法弹性波场三分量正演模拟中遇到的人为边界反射问题进行了研究,分别获得2-D横向各向同性介质中U_x,U_z,与U_y位移分量模拟时的边界吸收条件.从吸收边界上反射系数与入射角之间的关系可以说明,本文提出的边界吸收条件的吸收能力比较强.     

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 …     

THE APPLICATION OF PML BOUNDARY CODITIONS IN THE SIMULATION OF SEISMIC WAVE BY THE HIGH-ORDER STAGGERED-GRID FINITE DIFFERENCES AT THE TRANSVERSELY ISOTROPIC MEDIA

In the realm of the numerical simulation, finite difference method and finite element method are more intuitive and effective than other simulation methods. In the process of simulating seismic wave propagation, the finite differences method is widely used because of its high computational efficiency and the advantage of the algorithm is more efficient. With the demand of precision, more and more researchers have proposed more effective methods of finite differences, such as the high-order staggered-grid finite differences method, which can restore the actual process of wave propagation on the premise of ensuring accuracy and improving the efficiency of operation. In the past numerical simulation of seismic wave field, different models of isotropic medium are mostly used, but it is difficult to reflect the true layer situation. With the research demand of natural seismology and seismic exploration, the research on anisotropic media is more and more extensive. Transversely isotropic(TI)media can well simulate the seismic wave propagation in the formation medium, such as gas-bearing sandstone, mudstone, shale et al., the character of TI media is reflected by introducing the Thomsen parameters to reflect its weak anisotropy of vertical direction by using Thomson parameter. Therefore, studying the process of seismic wave propagation in TI media can restore the true information of the formation to the greatest extent, and provide a more reliable simulation basis for the numerical simulation of seismic wave propagation. In the geodynamic simulation and the numerical simulation of the seismic wave field, under the limited influence of the calculation area, if no boundary conditions are added, a strong artificial boundary reflection will be generated, which greatly reduces the validity of the simulation. In order to minimize the influence of model boundaries on the reflection of seismic waves, it is often necessary to introduce absorbing boundary conditions. At present, there are three types of absorption boundary conditions: one-way wave absorption boundary, attenuation absorption boundary, and perfectly matched layer(PML)absorption boundary. In terms of numerical simulation of seismic waves, the boundary absorption effect of PML is stronger than the first two, which is currently the most commonly used method, and it also represents the cutting-edge development direction of absorption boundary technology. The perfectly matched layer absorbing boundary is effectively applied to eliminating the reflective waves from model boundaries, but for transversely isotropic medium, the effect of the absorbing is not very well. For this reason, the elastic dynamic wave equations in transversely isotropic media are derived, and we describe a second-order accurate time, tenth-order accurate space, formulation of the Madariaga-Virieux staggered-grid finite difference methods with the perfectly matched layer(PML)are given. In addition, we have established vertical transversely isotropic(VTI)media and arbitrary inclined tilted transversely isotropic(TTI)media models, using a uniform half-space velocity model and a two-layer velocity model, respectively. By combining the actual geoscience background, we set the corresponding parameters and simulation conditions in order to make our model more research-oriented. When setting model parameters, different PML thickness, incident angle, source frequency and velocity layer models were transformed to verify the inhibition of boundary reflection effect by PML absorption boundary layer. The implementations of this simulation show that the formula is correct and for the transversely isotropic(TI)media of any angular symmetry axis, when the thickness of the PML layer reaches a certain value, the seismic wave reflection effect generated by the artificial boundary can be well suppressed, and the absorption effect of PML is not subject to changes in incident angle and wave frequency. Therefore, the results of our study indicate that our research method can be used to simulate the propagation process of seismic waves in the transversely isotropic(TI)media without being affected by the reflected waves at the model boundary to restore the actual formation information and more valuable geological research.     

双相各向同性介质弹性波高精度波场分离数值模拟方法

为了便于研究双相介质固流相混合弹性波场中纵横波波场的传播规律,提出了基于交错网格的Biot双相各向同性介质弹性波动方程高精度波场分离正演数值模拟方法.采用高阶交错网格有限差分法来构建一阶双曲型双相各向同性介质弹性波动方程正演算子实现波场正演,并在每一步递推过程中,分别计算出同相和流相分量相应的散度场(纯纵波场)和旋度场...     

准各向异性粘弹介质地震波的数字仿真

本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。     

准各向异性粘弹介质地震波的数字仿真

本文论述准各向异性粘弹波动方程有限元数值解地震波响应的原理与计算方法,并给出地震衰减因子与方程中粘滞系数的计算关系式。列出了所研制的地面记录地震剖面与井中记录垂直地震剖面,以及波场时间切片图、散度与旋度图、应变能与动能的时间切片图。本方法适于任意分层形态且层间物性参数为强间断的实际地震剖面的仿真,旨在研究正问题以实现反问题。     

弹性波数值模拟的非规则网格差分法

基于应力、速度混合变量弹性波方程及任意四边形网格差分算子,给出了交错计算应力及速度的非规则网格弹性波应力一速度差分法该方法融合了有限元法能适应复杂形状边界及差分法无需计算刚度阵的特点,具有较高的计算精度,所需计算机存储空间较少,计算效率也很高．基于积分平衡方程引入了任意形状自由表面的边界条件,且通过局部滤波改善了自由表面边界条件的稳定性,使得该方法可应用于考虑地表形状影响的地震波数值模拟     

弹性波传播数值模拟的区域分裂法

本文针对弹性波波场传播的数学模型,提出了一种基于有限元和有限差分耦合的区域分裂方法．它有灵活的网格剖分,克服了单纯用差分方法对区域的依赖性,可以适用于曲边地表;达到同样精度所需的计算量比有限元方法小;并易于实现并行计算．数值实验表明了算法的有效性．