Pure qP-wave reverse time migration using low rank finite difference method with spherical difference stencil
YANG Lisheng1, HUANG Jinqiang1,2, GAO Guochao1,2, XIA Peng1,2, HE Yunchuan2, WU Hao1
1. College of Resources and Environmental Engineering, Guizhou University, Guiyang, Guizhou 550025, China; 2. Key Laboratory of Karst Georesoucres and Environment (Guizhou University), Ministry of Education, Guiyang, Guizhou 550025, China
Abstract:Pseudo shear-wave artifacts and numerical instability exist in forward modeling and reverse time migration of quasi-acoustic equations for anisotropic media, producing poor imaging quality. Meanwhile, computational efficiency is an important factor in anisotropic reverse time migration. Therefore, a pure qP-wave reverse time migration method is proposed, where both imaging quality and computational efficiency are taken into account. Firstly, based on the pseudo-analytic solution of the wave equation and the exact anisotropic dispersion relation, the wave field continuation formula of the pseudo-analytic solution of pure qP-wave in VTI and TTI media is constructed respectively, which avoids the derivation of the explicit wave equation and does not need to square the dispersion relation, thus eliminating the pseudo-shear wave artifacts existing in the wave field simulation. Then, by designing a three-dimensional (3D) spherical or two-dimensional (2D) circular difference stencil and solving the difference coefficients appropriate to the model with the help of the low-rank finite difference method, a pure qP-wave field continuation formula based on low-rank finite difference is derived, which reduces the computational complexity by eliminating the need for multiple Fourier transforms during continuation along the time direction. On this basis, the GPU is used to calculate the forward and backward seismic wave fields in parallel, which improves the imaging efficiency of the reverse time migration. The experimental results of typical 2D and 3D models show that the proposed method can eliminate the pseudo shear-wave artifacts, ensure the stability of the calculation process, and have strong adaptability and high computational efficiency in both 2D and 3D anisotropic media.
FLETCHER R P, DU X, FOWLER P J. Reverse time migration in titled transversely isotropic (TTI) media[J]. Geophysics, 2009, 74(6):WCA179-WCA187.
[2]
ALKHALIFAH T.An acoustic wave equation for anisotropic media[J].Geophysics, 2000, 65(4):1239-1250.
[3]
杨富森, 李振春, 王小丹.TTI介质一阶qP波稳定方程波场数值模拟及逆时偏移[J].石油地球物理勘探, 2016, 51(3):497-505.YANG Fusen, LI Zhenchun, WANG Xiaodan.Wave-field numerical simulation and reverse-time migration based on the stable TTI first-order qP-wave equations[J].Oil Geophysical Prospecting, 2016, 51(3):497-505.
[4]
DUVENECK E, BAKKER P M.Stable P-wave modeling for reverse-time migration in tilted TI media[J].Geophysics, 2011, 76(2):S65-S75.
[5]
DU X, BANCROFT J C, LINES L R.Anisotropic reverse-time migration for tilted TI media[J].Geophysical Prospecting, 2007, 55(6):853-869.
[6]
LIU F Q, MORTON S A, JIANG S S, et al.Decoupled wave equation for P- and SV-waves in an acoustic VTI media[C].SEG Technical Program Expanded Abstracts, 2009, 28:2844-2848.
[7]
PESTANA R C, URSIN B, STOFFA P L.Separate P- and SV-wave equations for VTI media[C].SEG Technical Program Expanded Abstracts, 2011, 30:163-167.
[8]
CHU C, MACY B K, ANNO P D.Approximation of pure acoustic seismic wave propagation in TTI media[J].Geophysics, 2011, 76(5):WB97-WB107.
[9]
黄建平, 黄金强, 李振春, 等.TTI介质一阶qP波方程交错网格伪谱法正演模拟[J].石油地球物理勘探, 2016, 51(3):487-496.HUANG Jianping, HUANG Jinqiang, LI Zhenchun, et al.Staggered grid pseudo-spectral numerical simulation of first-order quasi P-wave equation for TTI media[J].Oil Geophysical Prospecting, 2016, 51(3):487-496.
[10]
ZHAN G, PESTANA R C, STOFFA P L.An efficient hybrid pseudo spectral/finite-difference scheme for solving the TTI pure P-wave equation[J].Journal of Geophysics and Engineering, 2013, 10(2):1322-1327.
[11]
张庆朝, 朱国维, 周俊杰, 等.TTI介质qP波伪谱法正演模拟[J].石油地球物理勘探, 2019, 54(2):302-311.ZHANG Qingchao, ZHU Guowei, ZHOU Junjie, et al.qP-wave numerical simulation in TTI media with pseudo-spectral method[J].Oil Geophysical Prospecting, 2019, 54(2):302-311.
[12]
XU S, ZHOU H B.Accurate simulations of pure quasi-P-waves in complex anisotropic media[J].Geophysics, 2014, 79(6):T341-T348.
[13]
LI X Y, ZHU H J.A finite-difference approach for solving pure quasi-P-wave equations in transversely isotropic and orthorhombic media[J].Geophysics, 2018, 83(4):C161-C172.
[14]
徐世刚, 包乾宗, 任志明, 等.结合泊松算法和波场分解成像条件的各向异性优化纯声波方程逆时偏移[J].石油地球物理勘探, 2022, 57(3):613-623.XU Shigang, BAO Qianzong, REN Zhiming, et al.Reverse-time migration of optimized pure acoustic equation in anisotropic media by combining Poisson algorithm and wavefield decomposition imaging condition[J].Oil Geophysical Prospecting, 2022, 57(3):613-623.
[15]
黄金强, 李振春.基于Low-rank分解的复杂TI介质纯qP波正演模拟与逆时偏移[J].地球物理学报, 2017, 60(2):704-721.HUANG Jinqiang, LI Zhenchun.Modeling and reverse time migration of pure quasi-P-wave in complex TI media with a low-rank decomposition[J].Chinese Journal of Geophysics, 2017, 60(2):704-721.
[16]
顾汉明, 张奎涛, 刘春成, 等.基于Low-rank一步法波场延拓的黏声各向异性介质纯qP波正演模拟[J].石油地球物理勘探, 2020, 55(4):733-746.GU Hanming, ZHANG Kuitao, LIU Chuncheng, et al.Low-rank one-step wave extrapolation for pure qP-wave forward modeling in viscoacoustic anisotropic media[J].Oil Geophysical Prospecting, 2020, 55(4):733-746.
[17]
FOMEL S, YING L X, SONG X L.Seismic wave extrapolation using low rank symbol approximation[C].SEG Technical Program Expanded Abstracts, 2010, 29:3092-3096.
[18]
SONG X L, ALKHALIFAH T.Modeling of pesudoacoustic P-waves in orthorhombic media with a low-rank approximation[J].Geophysics, 2013, 78(4):C33-C40.
[19]
SUN J, FOMEL S, YING L X.Low-rank one-step wave extrapolation for reverse time migration[J].Geophysics, 2016, 81(1):S39-S54.
[20]
SONG X L, FOMEL S, YING L X.Low-rank finite-differences and low-rank Fourier finite-differences for seismic wave extrapolation in the acoustic approximation[J].Geophysical Journal International, 2013, 193(2):960-969.
[21]
方刚, FomelS, 杜启振.交错网格Lowrank有限差分及其在逆时偏移中的应用[J].中国石油大学学报(自然科学版), 2014, 38(2):44-51.FANG Gang, FOMEL S, DU Qizhen.Low-rank finite difference on a staggered grid and its application on reverse time migration[J].Journal of China University of Petroleum (Edition of natural Science), 2014, 38(2):44-51.
[22]
袁雨欣, 胡婷, 王之洋, 等.求解二阶解耦弹性波方程的低秩分解法和低秩有限差分法[J].地球物理学报, 2018, 61(8):3324-3333.YUAN Yuxin, HU Ting, WANG Zhiyang, et al.Solving second-order decoupled elastic wave equation using low-rank decomposition and low-rank finite differences[J].Chinese Journal of Geophysics, 2018, 61(8):3324-3333.
[23]
黄金强, 李振春, 江文.TTI介质Low-rank有限差分法高效正演模拟及逆时偏移[J].石油地球物理勘探, 2018, 53(6):1198-1209.HUANG Jinqiang, LI Zhenchun, JIANG Wen.An efficient forward modeling with the low-rank finite-difference algorithm for complex TTI media and its application in reverse time migration[J].Oil Geophysical Prospecting, 2018, 53(6):1198-1209.
[24]
ZHANG Y, ZHANG G.One-step extrapolation method for reverse time migration[J].Geophysics, 2009, 74(4):A29-A33.
[25]
黄金强, 李振春.各向异性介质Low-rank有限差分法纯qP波叠前平面波最小二乘逆时偏移[J].地球物理学报, 2019, 62(8):3106-3129.HUANG Jinqiang, LI Zhenchun.Least-squares reverse time migration of pure qP-wave pre-stack plane-waves based on low-rank finite difference for anisotropic media[J].Chinese Journal of Geophysics, 2019, 62(8):3106-3129.
[26]
包乾宗, 戴雪, 梁雪.一种基于新差分模板和无穷范数的震源波场重建方法[J].石油地球物理勘探, 2022, 57(6):1384-1394.BAO Qianzong, DAI Xue, LIANG Xue.Source wavefield reconstruction based on a new finite-difference stencil and infinity norm[J].Oil Geophysical Prospecting, 2022, 57(6):1384-1394.