Abstract:The conventional staggered grid finite difference method is a numerical simulation method which is often used to study seismic wave propagation.However, different components of each variable in the staggered grid is staggered defined.For undefined points, their variables need to be interpolated, which reduces the accuracy of simulation.Therefore, the 2D 3C velocity-stress elastic wave equation of the TTI media is deduced on the basis of previous studies and the Lebedev grid is used to process the high accuracy differential separation.This scheme can avoid errors caused by the wave field interpolation in the conventional staggered grid when dealing with anisotropic media and it improves the modulated precision.Meanwhile, the multi-axis perfectly matched layer (MPML) method is introduced into this method.The forward modeling of the single layer TTI medium and the complex TTI medium with the lens is carried out.The results indicate as follows:A.The Lebedev grid does not need to interpolate the wavefield when the elastic wave equation of the anisotropic medium is discrete.Compared with the conventional staggered grid method, the simulation accuracy is higher.The test shows that P-wave, fast S-wave and slow S-wave can be clearly observed when the TTI medium is simulated and the polarization direction of fast S-wave and slow S-wave is opposite.The velocity characteristics of the three kinds of waves observed in single gun record also conform to the propagation regulation of seismic wavefield; B.After introducing MPML into this method, the phenomenon of boundary reflection is effectively suppressed without affecting the simulation effect.
Tsvankin I,Thomsen L.Inversion of reflection traveltimes for transverse isotropy.Geophysics,1995,60(4):1095-1107.
[2]
Alterman Z,Karal F C.Propagation of elastic waves in layered media by finite difference methods.Bulletin of the Seismological Society of America,1968,58(1):367-398.
[3]
Virieux J.SH wave propagation in heterogeneous media:velocity-stress finite-difference method.Geophysics,1984,49(11):1933-1942.
[4]
Virieux J.P-SV wave propagation in heterogeneous media:velocity-stress finite-difference method.Geophysics,1986,51(4):889-901.
[5]
Dablain M A.The application of high-order differencing to scalar wave equation.Geophysics,1986,51(1):54-66.
[6]
Graves R W.Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences.Bulletin of the Seismological Society of America,1996,86(4):1091-1106.
[7]
董良国,马在田,曹景忠等.一阶弹性波方程交错网格高阶差分解法.地球物理学报,2000,43(3):411-419.Dong Liangguo,Ma Zaitian,Cao Jingzhong et al.Study on stability of staggered grid high order difference method for first order elastic wave equation.Chinese Journal of Geophysics,2000,43(3):411-419.
[8]
裴正林.三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟.中国石油大学学报(自然科学版),2004,28(5):23-29.Pei Zhenglin.Numerical simulation of elastic wave equation with staggered grid high-order finite difference method in three dimensional anisotropic medium.Journal of China University of Petroleum(Edition of Natural Science),2004,28(5):23-29.
[9]
裴正林,董玉珊,彭苏萍.裂隙煤层弹性波场方位各向异性特征数值模拟研究.石油地球物理勘探,2007,42(6):665-672.Pei Zhenglin,Dong Yushan,Peng Suping.Numerical simulation of azimuthal anisotropy of elastic wave field in fractured seams.OGP,2007,42(6):665-672.
[10]
裴正林.层状各向异性介质中横波分裂和再分裂数值模拟.石油地球物理勘探,2006,41(1):17-25.Pei Zhenglin.Numerical simulation of splitting and splitting of transverse wave in layered anisotropic media.OGP,2006,41(1):17-25.
[11]
李东,董守华,赵小翠等.煤储层双相EDA介质的地震波场模拟.地球物理学进展,2011,26(2):654-663.Li Dong,Dong Shouhua,Zhao Xiaocui et al.Simulation of seismic wave field in coal reservoir with dual phase EDA medium.Progress in Geophysics,2011,26(2):654-663.
镇晶晶,刘洋,李敏.几种裂缝模型的波场传播特征比较.石油地球物理勘探,2012,47(6):908-917.Zhen Jingjing,Liu Yang,Li Min.Comparison of wave field characteristics among the fracture rock models.OGP,2012,47(6):908-917.
[14]
Lebedev V I.Difference analogues of orthogonal decompositions,basic differential operators and some boundary problems of mathematical physics I.USSR Computational Mathematics and Mathematical Phy-sics,1964,4(3):69-92.
[15]
Bernth H,Chapman C.A comparison of finite-difference grids for anisotropic elastic modeling.SEG Technical Program Expanded Abstracts,2010,29:2916-2920.
[16]
Lisitsa V,Vishnevskiy D.Lebedev scheme for the numerical simulation of wave propagation in 3D anisotropic elasticity.Geophysical Prospecting,2010,58(4):619-635.
[17]
林朋,彭苏萍,卢勇旭等.基于旋转交错网格的双相各向异性介质二维三分量波场模拟.煤炭学报,2016,41(5):1203-1211.Lin Peng,Peng Suping,Lu Yongxu et al.Study on 2D/3C wave propagation in two-phase anisotropic media using the rotated staggered-grid method.Journal of China Coal Society,2010,58(4):619-635.
[18]
李娜,黄建平,李振春等.Lebedev网格改进差分系数TTI介质正演模拟方法研究.地球物理学报,2014,57(1):261-269.Li Na,Huang Jianping,Li Zhenchun et al.Research on Lebedev mesh modified differential coefficient TTI forward modeling method.Chinese Journal of Geophysics,2014,57(1):261-269.
[19]
李娜,李振春,黄建平等.Lebedev网格与标准交错网格耦合机制下的复杂各向异性正演模拟.石油地球物理勘探,2014,49(1):121-131.Li Na,Li Zhenchun,Huang Jianping et al.Complex anisotropic forward modeling based on Lebedev mesh and standard staggered mesh coupling mechanism.OGP,2014,49(1):121-131.
[20]
黄建平,杨宇,李振春等.基于M-PML边界的Lebedev网格起伏地表正演模拟方法及稳定性分析.中国石油大学学报(自然科学版),2016,40(4):47-56.Huang Jianping,Yang Yu,Li Zhenchun et al.Forward modeling method and stability analysis of Lebedev grid based on M-PML boundary.Journal of China University of Petroleum(Edition of Natural Science),2016,40(4):47-56.
[21]
杨宇,黄建平,雷建设等.Lebedev网格黏弹性介质起伏地表正演模拟.石油地球物理勘探,2016,51(4):698-706.Yang Yu,Huang Jianping,Lei Jianshe et al.Forward modeling of viscoelastic medium in Lebedev lattice.OGP,2016,51(4):698-706.
[22]
杨庆节,刘财,耿美霞等.交错网格任意阶导数有限差分格式及差分系数推导.吉林大学学报,2014,44(1):375-385.Yang Qingjie,Liu Cai,Geng Meixia et al.Staggered grid finite difference scheme and coefficients deduction of any number of derivatives.Journal of Jilin University(Earth Science Edition),2014,44(1):375-385.
[23]
王洋,刘洪,张衡等.一种全局优化的隐式交错网格有限差分算法及其在弹性波数值模拟中的应用.地球物理学报.2015,58(7):2508-2524.Wang Yang,Liu Hong,Zhang Heng et al.A global optimized implicit staggered grid finite difference algorithm and its application in numerical simulation of elastic wave.Chinese Journal of Geophysics,2015,58(7):2508-2524.
[24]
尹志恒,狄帮让,魏建新等.裂隙参数对纵波能量衰减影响的物理模型研究.石油地球物理勘探,2012,47(5):728-734.Yin Zhiheng,Di Bangrang,Wei Jianxin et al.P-wave attenuation by fracture parameter on physical models.OGP,2012,47(5):728-734.
[25]
沈金松,詹林森,马超.裂缝等效介质模型对裂缝结构和充填介质参数的适应性.吉林大学学报:地球科学版,2013,43(3):993-1003.Shen Jinsong,Zhan Linsen,Ma Chao.The applicability of fracture effective model to diffeerent parameters of fracture geometric strutures and media filled in fracture pores.Journal of Jinlin University(Earth Science Edition),2013,43(3):993-1003.
[26]
姚振岸,孙成禹,谢俊法等.黏弹TTI介质旋转交错网格微地震波场模拟.石油地球物理勘探,2017,52(2):253-263.Yao Zhen'an,Sun Chengyu,Xie Junfa et al.Simulation of micro-seismic wave field in viscous TTI medium rotating staggered.OGP,2017,52(2):253-263.