黏弹TTI介质中旋转交错网格高阶有限差分数值模拟
详细信息 本馆镜像全文    |  推荐本文 | | 获取馆网全文
摘要
以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.
On the basis of Carcione′s theories of viscoelasticity and anisotropy,two-dimensional,three-component,first-order velocity-stress wave equations of viscoelastic tilted transversely isotropic(viscoelastic TTI) media are presented and a rotated staggered grid any-order finite-difference scheme is used to numerically solve the equations.Equations of the perfectly matched layer(PML) are derived for the wave equations in viscoelastic TTI media and the rotated staggered grid any-order finite-difference scheme is also used to solve these equations.Results of numerical modeling indicate that the modeling precision is high and the absorbing boundary condition is good in the viscoelastic TTI media,and high-precision snapshots of wave field and synthetic seismograms can be obtained,and they can reflect the characteristic of viscoelasticity and anisotropy.
引文
[1]Carcione J M.Wave propagation in anisotropic linearviscoelastic media:Theory and simulated wavefields.GeophysicalJournal International,1990,101(3):739-750.
    [2]Carcione J M.Constitutive model and wave equations forlinear,viscoelastic,anisotropic media.Geophysics,1995,60(2):537-548.
    [3]Carcione J M.Staggered mesh for the anisotropic andviscoelastic wave equation.Geophysics,1999,64(6):1863-1866.
    [4]Zhang Z J,Wang G J,Harris J M.Multi-componentwavefield simulation in viscous extensively dilatancy anisotropicmedia.Physics of the Earth and Planetary Interiors,1999,114(1-2):25-38.
    [5]张忠杰,滕吉文,贺振华.EDA介质中地震波速度、衰减与品质因子方位异性研究.中国科学(E辑),1999,29(6):569-574.Zhang Z J,Teng J W,He Z H.Azimuthal anisotropy ofseismic velocity,attenuation and Q value in viscous EDAmedia.Science in China(Series E)(in Chinese),1999,29(6):569-574
    [6]杜启振,杨慧珠.方位各向异性黏弹性介质波场有限元模拟.物理学报,2003,52(8):2010-2014.Du Q Z,Yang H Z.Finite-element methods for viscoelasticand azimuthally anisotropic media.Acta Physica Sinica(inChinese),2003,52(8):2010-2014.
    [7]杜启振.各向异性黏弹性介质伪谱法波场模拟.物理学报,2004,53(12):4428-4434.Du Q Z.Wavefield forward modeling with the pseudo-spectral method in viscoelastic and azimuthally anisotropicmedia.Acta Physica Sinica(in Chinese),2004,53(12):4428-4434.
    [8]杜启振,王延光,付水华.方位各向异性黏弹性介质波场数值模拟.地球物理学进展,2006,21(2):502-504.Du Q Z,Wang Y G,Fu S H.Wavefield forward modelingwith the pseudo-spectral method in viscoelastic and azimuthallyanisotropic media.Progress in Geophysics(in Chinese),2006,21(2):502-504.
    [9]郭智奇,刘财,杨宝俊等.黏弹各向异性介质中地震波场模拟与特征.地球物理学进展,2007,22(3):804-810.Guo Z Q,Liu C,Yang B J.Seismic wave fields modeling andfeature in viscoelastic anisotropic media.Progress in Geophysics(in Chinese),2007,22(3):804-810.
    [10]谭未一,赵兵,张中杰.黏弹性VTI介质地震波模拟特征分析.地球物理学进展,2007,22(4):1305-1311.Tan W Y,Zhao B,Zhang Z J.The analysis of the wave fieldcharacteristics in viscoelastic VTI medium.Progress inGeophysics(in Chinese),2007,22(4):1305-1311.
    [11]李军,苏云,李录明.各向异性黏弹性介质波场数值模拟.大庆石油学院学报,2009,33(6):38-42.Li J,Su Y,Li L M.Wavefield numeric simulation inanisotropic viscoelastic media.Journal of Daqing PetroleumInstitute(in Chinese),2009,33(6):38-42.
    [12]刘洋,李承楚,牟永光.具有倾斜对称轴的横向各向同性介质中的弹性波.石油地球物理勘探,1998,33(2):161-169.Liu Y,Li C C,Mou Y G.Elastic wave in transverseisotropic media with dipping symmetric axis.Oil GeophysicalProspecting(in Chinese),1998,33(2):161-169.
    [13]Qin Y L,Zhang Z J,Li S L.CDP mapping in tiltedtransversely isotropic(TTI)media.Part I:Method andeffectiveness.Geophysical Prospecting,2003,51(4):315-324.
    [14]Qin Y L,Zhang Z J,Xu S H.CDP mapping in tiltedtransversely isotropic(TTI)media.Part II:Velocityanalysis by combining CDP mapping with a geneticalgorithm.Geophysical Prospecting,2003,51(4):325-332.
    [15]Virieux J.P-SV wave propagation in heterogeneous media:velocity-stress finite-difference method.Geophysics,1986,51(4):889-901.
    [16]冯英杰,杨长春,吴萍.地震波有限差分模拟综述.地球物理学进展,2007,22(2):487-491.Feng Y J,Yang C C,Wu P.The review of the finite-difference elastic wave motion modeling.Progress inGeophysics(in Chinese),2007,22(2):487-491.
    [17]Saenger E H,Gold N,Shapiro S A.Modeling thepropagation of elastic waves using a modified finite-differencegrid.Wave Motion,2000,31(1):77-92.
    [18]Saenger E H,Bohlen T.Finite-difference modeling ofviscoelastic and anisotropic wave propagation using therotated staggered grid.Geophysics,2004,69(2):583-591.
    [19]陈浩,王秀明,赵海波.旋转交错网格有限差分及其完全匹配层吸收边界条件.科学通报,2006,51(17):1985-1994.Chen H,Wang X M,Zhao H B.Rotated staggered gridfinite-difference and the perfectly matched layer.ChineseScience Bulletin(in Chinese),2006,51(17):1985-1994.
    [20]张鲁新,符力耘,裴正林.不分裂卷积完全匹配层与旋转交错网格有限差分在孔隙弹性介质模拟中的应用.地球物理学报,2010,53(10):2470-2483.Zhang L X,Fu L Y,Pei Z L.Finite difference modeling ofBiot's poroelastic equations with unsplit convolutional PMLand rotated staggered grid.Chinese J.Geophys.(in Chinese),2010,53(10):2470-2483.
    [21]杜启振,孙瑞艳,张强.组合边界条件下二维三分量TTI介质波场数值模拟.石油地球物理勘探,2011,46(2):187-195.Du Q Z,Sun R Y,Zhang Q.Numerical simulation of three-component elastic wavefield in 2DTTI media in the conditionof the combined boundary.Oil Geophysical Prospecting(inChinese),2011,46(2):187-195.
    [22]Liu Y,Sen M K.An implicit staggered-grid finite-differencemethod for seismic modelling.Geophysical JournalInternational,2009,179(1):459-474.
    [23]Igel H,Mora P,Riollet B.Anisotropic wave propagationthrough finite-difference grids.Geophysics,1995,60(4):1203-1216.
    [24]Berenger J P.A perfectly matched layer for the absorption ofelectromagnetic waves.Journal of Computational Physics,1994,114(2):185-200.
    [25]Collino F,Tsogka C.Application of the perfectly matchedabsorbing layer model to the linear elastodynamic problem inanisotropic heterogeneous media.Geophysics,2001,66(1):294-307
    [26]Zhang Z J,Teng J W,Wan Z C.Establishment ofparameteric equations of seismic wave fronts in 2-D anisotropicmedia.Chinese Science Bulletin,1994,39(22):1890-1894.
    [27]裴正林,王尚旭.任意倾斜各向异性介质中弹性波波场交错网格高阶有限差分法模拟.地震学报,2005,27(4):441~451.Pei Z L,Wang S X.A staggered-grid high-order finite-differencemodeling for elastic wave field in arbitrary title aniotropic media.Acta Seismologica Sinica(in Chinese),2005,27(4):441~451.

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