朱文新,何兵寿

• 1:中国海洋大学
• 2:海底科学与探测技术教育部重点实验室

摘要(Abstract)：

数值模拟是研究地震波传播机理的重要工具,也是地震资料采集、处理、和反演领域的重要内容。基于常规矩形网格或交错网格的差分算法存在模拟精度和计算量之间的矛盾,变网格技术可以缓解这一矛盾,但对于TTI介质,常规变网格、变系数差分算法无法实现各分量的交错对应,旋转网格与旋转交错网格技术存在同样的问题。针对上述问题,提出了一种适用于TTI介质正演的新方法:即菱形交错网格剖分法,并给出了基于菱形交错网格的变系数、变网格方法,实现了TTI介质中弹性波场的高精度模拟。

关键词(KeyWords)： TTI介质;弹性波;数值模拟;变网格;菱形网格

Abstract:

Keywords:

基金项目(Foundation): 国家重点研发计划(编号:2018YFC1405900);; 中央高校基本科研业务费专项(编号:201822011);; 国家自然科学基金(编号:41674118,41574105);; 国家重大科技专项(2016ZX05027-002-005)联合资助

作者(Author): 朱文新,何兵寿

参考文献(References)：

• [1]Crampin S. A review of wave motion in anisotropic and cracked elastic-media[J]. Wave motion,1981,3(4):343-391
• [2]Thomsen L. Weak elastic anisotropy[J]. Geophysics,1986,51(10):1954-1966.
• [3]Winterstein D F. Velocity anisotropy terminology for geophysicists[J]. Geophysics,1990,55(8):1070-1088.
• [4]李佳珂，张会星，白冰，等．TTI介质纯准P波一阶压力-速度方程及求解方法[J]．中国煤炭地质，2018,30(7):76-82.
• [5]刘学义，何兵寿．VTI介质中弹性波方程正演的一阶混合吸收边界[J]．中国煤炭地质，2015,27(3):59-63+76.
• [6]张晶，何兵寿，张会星．含直立裂隙介质的弹性波动方程正演模拟[J]．中国煤炭地质，2008,20(1):50-54.
• [7]Moczo P. Finite-difference technique for SH-waves in 2-D media using irregular grids—application to the seismic response problem[J].Geophysical Journal International,2010,99(2):321-329.
• [8]Jastram C,BehleA. Acoustic modeling on a grid of vertically varying spacing[J]. Geophysics Prospecting,1992,40(2):157-169
• [9]Jastram C,Tessmer E. Elastic modeling on a grid with vertically varying spacing[J]. Geophysics Prospecting,1994,42:357-370
• [10]Falk J,Tessmer E,Gajewski D. Tube wave modeling by the finite-difference method with varying grid spacing[J]. Pageoph,1996,148:77-93
• [11]Falk J,Tessmer E,Gajewski D. Efficient finite-difference modelling of seismic waves using locally adjustable time step sizes[J]. Geophys. Prosp,1998,46:603-616
• [12]朱生旺，曲寿利，魏修成，等．变网格有限差分弹性波方程数值模拟方法[J]．石油地球物理勘探，2007,42(6):634-641
• [13]李振春，张慧，张华．一阶弹性波方程的变网格高阶有限差分数值模拟[J]．石油地球物理勘探，2008,43(6):711-717.
• [14]李胜军，孙成禹．声波方程有限差分数值模拟的变网格步长算法[J]．工程地球物理学报，2007,4(3):207-212.
• [15]解闯，宋鹏，谭军，等．声波方程变网格有限差分正演模拟的虚假反射分析[J]．地球物理学进展，2019(2)
• [16]高阳，冯西会，汶小刚，等．三角网格模型剖分方法的研究与应用[J]．中国煤炭地质，2013,25(10):48-52.
• [17]李娜，黄建平，李振春，等．Lebedev网格改进差分系数TTI介质正演模拟方法研究[J]．地球物理学报，2014,57(1):261-269.
• [18]杨绍伟，何兵寿，杨佳佳．弹性波逆时偏移子波拉伸校正[J]．中国煤炭地质，2016,28(2):61-67+77.
• [19]范国苗，何兵寿．P-SV转换波真振幅DMO[J]．中国煤炭地质，2009,21(7):53-56.
• [20]张建敏，何兵寿，唐怀谷．三维TTI介质中的纯准P波方程及求解方法[J]. Applied Geophysics,2017(1):125-132.
• [21]Graves R W. Simulating Seismic Wave Propagation in 3D Elastic Media Using Staggered-Grid Finite Differences[J]. Bulletin of the Seismological Society of America,1996,86(4):1091-1106.
• [22]田雪丰．一阶弹性波交错网格时间高阶差分格式及稳定性分析[J]．中国煤炭地质，2019,31(5):70-78.
• [23]Igel H,Mora P,Riollet B. Anisotropic wave propagation through finite-difference grids[J]. Geophysics,1995,60(4):1203-1216.
• [24]Saenger E H,Gold N,Shapiro S A. Modeling the propagation of elastic waves using a modified finite-difference grid[J]. Wave Motion,2000,31(1):77-92.
• [25]董良国，马在田，曹景忠，等．一阶弹性波方程交错网格高阶差分解法[J]．地球物理学报，2000(03):411-419.
• [26]马光克，李洋森，孙万元，等．可变网格高阶有限差分法逆时偏移研究[J]．石油物探，2016,55(5):728-736.
• [27]董良国，马在田，曹景忠．一阶弹性波方程交错网格高阶差分解法稳定性研究[J]．地球物理学报，2000,43(6):856-864.