复杂介质多尺度变网格地震波数值模拟
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
在应用有限差分法地震波场数值模拟过程中,模拟精度和计算效率是关键问题之一,特别是对于速度变化剧烈的含小尺度地质体的数值模拟尤为重要.为了既能精细刻画介质的局部结构,又能保证模拟的时效性,本文在传统交错网格有限差分算法的基础上,将变网格有限差分方法引入到交错网格高阶差分数值模拟中,对交错网格空间算法进行了改进,并避免了因插值因素影响模拟精度和计算效率.本文采用完全匹配层吸收边界条件实现对可变网格数值模拟的边界处理,根据完全匹配层吸收边界原理和参量的分裂思想,在不同的方向分别加相应的阻尼衰减因子,解决边界的反射问题.此外还分析了数值频散对于模拟结果的影响.文中分别设计了地堑和生物礁两个介质模型对可变交错网格算法的有效性和稳定性进行了验证.数值模拟结果表明,可变交错网格有限差分算法数值模拟能够对地下介质物性的空间变化进行准确刻画,进一步增强数值模拟对复杂介质的适应性,同时为地震数据的波场成像、纵横波联合解释等提供可靠依据.
In the process of seismic wave field numerical simulation using finite difference method,the simulation accuracy and computational efficiency is one of the keys to the problem which is especially important to the numerical simulation of small scale geological body which velocity changes violently.In order to describe the local structure of medium subtly and guarantee the efficiency of the simulation,this article introduces the variable grid finite difference method to the staggered grid high-order finite difference numerical simulation on the basic of the traditional staggered grid finite difference algorithm to improve the staggered grid spatial algorithm and avoid the reduction of the simulation accuracy and computational efficiency caused by the interpolation factor.This article uses the PML absorbing boundary conditions for the boundary treatment of the variable grid numerical simulation that in different directions respectively introduce the corresponding damping factor to solve the problem of boundary reflection which according to the principle of perfectly matched layer absorbing boundary and the parameters of disunion.In addition,this article analyzes the affection of numerical frequency dispersion on the simulation results.In this article,,the graben and organic reef two medium model have been designed to verify the effectiveness and stability of variable staggered grid algorithm.The results show that the variable staggered grid numerical simulation of finite difference algorithm can accurately depict the space variation of underground medium physical properties to further enhance the adaptability of numerical simulation of complex medium,it also can provide reliable basis for wave field imaging and the combined interpretation of p-wave and s-wave.
引文
Alterman Z,Karal F C Jr.1968.Propagation of elastic waves in layered media by finite difference methods[J].Bull.Seism.Soc.Am.,58(1):367-398.
    Cai Q X,He P J,Qin G S,et al.2003.Minimum dispersion algorithm by finite-difference numeric modeling and its applications[J].Oil Geophysical Prospecting(in Chinese),38(3):247-251,262.
    Dong L G.1999.Absorptive boundary condition in elastic-wave numerical modeling[J].Oil Geophysical Prospecting(in Chinese),34(1):45-56.
    Dong L G,Ma Z T,Cao J Z.2000a.A study on stability of the staggered-grid high-order difference method of first-order elastic wave equation[J].Chinese Journal of Geophysics(in Chinese),43(6):856-864.
    Dong L G,Ma Z T,Cao J Z,et al.2000b.A Staggered-grid highorder difference method of one-order elastic wave equation[J].Chinese Journal of Geophysics(in Chinese),43(3):411-419.
    DuQZ,Li B.2009.Seismic wavefield simulation of VTI medium under combined absorption condition[J].Oil Geophysical Prospecting(in Chinese),44(2):173-178.
    Du Q Z,Bai Q Y,Li B.2010.Optimized difference coefficient method seismic wave field numerical simulation in vertical transverse isotropic medium[J].Oil Geophysical Prospecting(in Chinese),45(2):170-176.
    Feng Y J,Yang C C,Wu P.2007.The review of the finitedifference elastic wave motion modeling[J].Progress in Geophysics(in Chinese),22(2):487-491.
    Guo Y H.1994.The pool-formed rule of natural gas and diagenesis of carboniferous reservoir in East Sichuan[J].Journal of Southwestern Petroleum Institute(in Chinese),16(1):1-8.
    Hu Y,Chen H,Su Y,et al.2012.Seismic wave modelling of twophase media based on BISQ model[J].Oil Geophysical Prospecting(in Chinese),47(6):901-907.
    Jastram C,Behle A.1992.Acoustic modelling on a grid of vertically varying spacing[J].Geophysical Prospecting,40(2):157-169.
    Jastram C,Tessmer E.1994.Elastic modelling on a grid with vertically varying spacing[J].Geophysical Prospecting,42(4):357-370.
    Li W J,Wei X C,Ning J R.2009.A new numeric simulating absorbing boundary condition of elastic wave[J].Oil Geophysical Prospecting(in Chinese),44(4):501-507.
    Li Y J.1998.Oil and gas sources and forming time of oil gas pools in Tazhong area,Tarim basin[J].Petroleum Exploration and Development(in Chinese),25(1):11-14.
    Liu Y S,Teng T W,Liu S L,et al.2013.Explicit finite element method with triangle meshes stored by sparse format and its perfectly matched layers absorbing boundary condition[J].Chinese Journal of Geophysics(in Chinese),56(9):3085-3099,doi:10.6038/cjg20130921.
    Madariaga R.1976.Dynamics of an expanding circular fault[J].Bull.Seism.Soc.Am.,66(3):639-666.
    Meng F S,Li Y S,Wang L,et al.2013.Based on variable and staggered grids seismic numerical simulation[J].Progress in Geophysics(in Chinese),28(6):2936-2943,doi:10.6038/pg20130614.
    Moczo P.1989.Finite-difference technique for SH-waves in 2-D media using irregular grids-application to the seismic response problem[J].Geophys.J.Int.,99(2):321-329.
    Pei Z L,Mu Y G.2003.A staggered-grid high-order difference method for modeling seismic wave propagation in inhomogeneous media[J].Journal of the University of Petroleum,China(in Chinese),27(6):17-21.
    Pei Z L.2004.Numerical modeling using staggered-grid high-order finite-difference of elastic wave equation on arbitrary relief surface[J].Oil Geophysical Prospecting(in Chinese),39(6):629-634.
    Sun C Y,Li S J,Ni C K,et al.2008.Wave equation numerical modeling by finite difference method with varying grid spacing[J].Geophysical Prospecting for Petroleum(in Chinese),47(2):123-127.
    Sun C Y,Gong T J,Zhang Y L,et al.2009.Analysis on dispersion and alias in finite-difference solution of wave equation[J].Oil Geophysical Prospecting(in Chinese),44(1):43-48.
    Sun L J,Yin X Y.2011.A finite-difference scheme based on PML boundary condition with high power grid step variation[J].Chinese Journal of Geophysics(in Chinese),54(6):1614-1623,doi:10.3969/j.issn.0001-5733.2011.06.021.
    Virieux J.1984.SH-wave propagation in heterogeneous media;Velocity-stress finite-difference method[J].Geophysics,49(11):1933-1942.
    Virieux J.1986.P-SV wave propagation in heterogeneous media;Velocity-stress finite-difference method[J].Geophysics,51(4):889-901.
    Wang Y G,Xing W J,Xie W X,et al.2007.Study of absorbing boundary condition by perfectly matched layer[J].Journal of China University of Petroleum(in Chinese),31(1):19-24.
    Xue D C,Wang S X.2008.Using combined mass matrix to suppress numeric dispersion[J].Oil Geophysical Prospecting(in Chinese),43(3):318-320.
    Zhang J L,Tian Z P,Wang C X.2008.Elastic wave simulation of P-and S-waves separation by 2-D staggered grid and application[J].Oil Geophysical Prospecting(in Chinese),43(6):717-722.
    Zhu S W,Wei X C.2005.Differential forward modeling of wave equation having irregular grid and any-order precision[J].Oil Geophysical Prospecting(in Chinese),40(2):149-153.
    蔡其新,何佩军,秦广胜,等.2003.有限差分数值模拟的最小频散算法及其应用[J].石油地球物理勘探,38(3):247-251,262.
    董良国.1999.弹性波数值模拟中的吸收边界条件[J].石油地球物理勘探,34(1):45-56.
    董良国,马在田,曹景忠.2000a.一阶弹性波方程交错网格高阶差分解法稳定性研究[J].地球物理学报,43(6):856-864.
    董良国,马在田,曹景忠,等.2000b.一阶弹性波方程交错网格高阶差分解法[J].地球物理学报,43(3):411-419.
    杜启振,李宾.2009.组合吸收边界条件下VTI介质地震波场模拟[J].石油地球物理勘探,44(2):173-17 8.
    杜启振,白清云,李宾.2010.横向各向同性介质优化差分系数法地震波场数值模拟[J].石油地球物理勘探,45(2):170-176.
    冯英杰,杨长春,吴萍.2007.地震波有限差分模拟综述[J].地球物理学进展,22(2):487-491.
    郭一华.1994.川东地区石炭系储层成岩作用和天然气成藏规律[J].西南石油学院学报,16(1):1-8.
    胡英,陈辉,苏云,等.2012.基于双相介质BISQ模型的地震波正演模拟[J].石油地球物理勘探,47(6):901-907.
    李文杰,魏修成,宁俊瑞.2009.一种新的弹性波数值模拟吸收边界条件[J].石油地球物理勘探,44(4):501-507.
    李延钧.1998.塔中地区油气源及成藏时期研究[J].石油勘探与开发,25(1):11-14.
    刘有山,滕吉文,刘少林,等.2013.稀疏存储的显式有限元三角网格地震波数值模拟及其PML.吸收边界条件[J].地球物理学报,56(9):3085-3099,doi:10.6038/cjg20130921.
    孟凡顺,李洋森,王璐,等.2013.基于可变交错网格的地震波数值模拟[J].地球物理学进展,28(6):2936-2943,doi:10.6038/pg20130614.
    裴正林,牟永光.2003.非均匀介质地震波传播交错网格高阶有限差分法模拟[J].石油大学学报(自然科学版),27(6):17-21.
    裴正林.2004.任意起伏地表弹性波方程交错网格高阶有限差分法数值模拟[J].石油地球物理勘探,39(6):629-634.
    孙成禹,李胜军,倪长宽,等.2008.波动方程变网格步长有限差分数值模拟[J].石油物探,47(2):123-127.
    孙成禹,宫同举,张玉亮,等.2009.波动方程有限差分法中的频散与假频分析[J].石油地球物理勘探,44(1):43-48.
    孙林洁,印兴耀.2011.基于PML边界条件的高倍可变网格有限差分数值模拟方法[J].地球物理学报,54(6):1614-1623,doi:10.3969/j.issn.0001-5733.2011.06.021.
    王永刚,邢文军,谢万学,等.2007.完全匹配层吸收边界条件的研究[J].中国石油大学学报(自然科学版),31(1):19-24.
    薛东川,王尚旭.2008.利用组合质量矩阵压制数值频散[J].石油地球物理勘探,43(3):318-320.
    张建磊,田振平,王成祥.2008.二维交错网格纵横波分离的弹性波模拟及应用[J].石油地球物理勘探,43(6):717-722.
    朱生旺,魏修成.2005.波动方程非规则网格任意阶精度差分法正演[J].石油地球物理勘探,40(2):149-153.

版权所有:© 2023 中国地质图书馆 中国地质调查局地学文献中心