Elastic wavefield separation methods and their applications in reverse time migration
Wu Xiao1,3, Liu Yang1,2,3, Cai Xiaohui4
1. State Key Laboratory of Petroleum Resources and Prospecting, China University of Petroleum(Beijing), Beijing 102249, China; 2. Faculty of Petroleum, Karamay Campus of China University of Petroleum(Beijing), Karamay, Xinjiang 834000, China; 3. CNPC Key Laboratory of Geophysical Prospecting, China University of Petroleum(Beijing), Beijing 102249, China; 4. Institute of Geotechnical Engineering, Nanjing University of Technology, Nanjing, Jiangsu 210009, China
Abstract:The elastic wavefield separation is a key technique for elastic reverse time migration (RTM),which is mainly based on Helmholtz decomposition.In this paper,we analyze and summarize elastic wavefield separation methods,and classify them into two kinds:one is direct separation and the other is indirect separation.The former obtains P-and S-wavefields by directly calculating the divergence and curl of the original wavefield.This kind includes the space-domain methods and the wavenumber-domain methods.The latter obtains P-and S-wavefields at the next time using recursively the separated P-and S-wavefields of the previous time.We further introduce the basic principles of three space-domain methods,a wavenumber-domain method and an indirect separation method.We adopt a homogenous model,a layered model and the Marmousi model to demonstrate the effectiveness and efficiency of these five methods.The two conventional space-domain methods have the least computational complexity,but the separation results do not preserve the amplitude.The indirect separation method is effective and the calculation is moderate,but there are energy anomalies at velocity interfaces.The improved space-domain method has the best separation effect (but usually needs larger computational cost),followed by the wavenumber-domain method.Finally,the influences of these five methods on the elastic RTM results are tested with the Marmousi model,which indicate both the improved space-domain method and the wavenumber-domain method have the satisfactory imaging results.
杨午阳,张厚柱,撒利明等.逆时偏移有关问题探讨.石油地球物理勘探,2016,51(6):1251-1262.Yang Wuyang,Zhang Houzhu,Sa Liming et al.Some issues about reverse time migration.OGP,2016,51(6):1251-1262.
[2]
Baysal E,Kosloff D and Sherwood J.Reverse time migration.Geophysics,1983,48(11):1514-1524.
[3]
Chang W F and McMechan G A.Elastic reverse-time migration.Geophysics,1987,52(10):1365-1375.
[4]
丁亮,刘洋.逆时偏移成像技术研究进展.地球物理学进展,2011,26(3):1085-1100.Ding Liang and Liu Yang.Progress in reverse time migration imaging.Progress in Geophysics,2011,26(3):1085-1100.
[5]
Morse P M and Feshbach H.Methods of Theoretical Physics.McGraw-Hill Book Company,1953.
[6]
Sun R,Chow J and Chen K.Phase correction in separating P-and S-waves in elastic data.Geophysics,2001,66(5):1515-1518.
[7]
Sun R,McMechan G A,Hsiao H H et al.Separating P and S-waves in prestack 3D elastic seismograms using divergence and curl.Geophysics,2004,69(1):286-297.
[8]
Du Q,Zhu T and Ba J.Polarity reversal correction for elastic reverse time migration.Geophysics,2012,77(2):S31-S41.
[9]
Aki K and Richards P G.Quantitative Seismology:Theory and Methods.W H Freeman & Co,San Francisco,1980.
[10]
Zhu H.Elastic wavefield separation based on the Helmholtz decomposition.Geophysics,2017,82(2):S173-S183.
[11]
Zhang Q and McMechan G A.2D and 3D elastic wavefield vector decomposition in the wavenumber domain for VTI media.Geophysics,2010,75(3):D13-D26.
[12]
陈志文,贺振华.二维弹性波动方程F-K分离.石油物探,1992,31(2):47-57.Cheng Zhiwen and He Zhenhua.F-K decomposition of 2-D elastic wave equation.GPP,1992,31(2):47-57.
[13]
井涌泉,高红伟,王维红.二维各向同性介质P波和S波分离方法研究.地球物理学进展,2008,23(5):1412-1416.Jing Yongquan,Gao Hongwei and Wang Weihong.The study of P-waves and S-waves separation in two-dimensional isotropic media.Progress in Geophysics,2008,23(5):1412-1416.
[14]
马德堂,朱光明.弹性波波场P波和S波分解的数值模拟.石油地球物理勘探,2003,38(5):482-486.Ma Detang and Zhu Guangming.Numerical modeling of P-wave and S-wave separation in elastic wavefield.OGP,2003,38(5):482-486.
[15]
李振春,张华,刘庆敏等.弹性波交错网格高阶有限差分法波场分离数值模拟.石油地球物理勘探,2007,42(5):510-515.Li Zhenchun,Zhang Hua,Liu Qingmin et al.Numeric simulation of elastic wavefield separation by staggering grid high-order finite-difference algorithm.OGP,2007,42(5):510-515.
[16]
张建磊,田振平,王成祥.二维交错网格纵横波分离的弹性波模拟及应用.石油地球物理勘探,2008,43(6):717-722.Zhang Jianlei,Tian Zhenping and Wang Chengxiang.Elastic wave simulation of P-and S-waves separation by 2-D staggered grid and application.OGP,2008,43(6):717-722.
[17]
刘洪林,陈可洋,杨薇等.高阶交错网格有限差分法纵横波波场分离数值模拟.地球物理学进展,2010,25(3):877-884.Liu Honglin,Chen Keyang,Yang Wei et al. Numerical modeling of P-and S-wave field separation with high-order finite-difference scheme.Progress in Geophysics,2010,25(3):877-884.
[18]
唐小平,白超英,刘宽厚.伪谱法分离波动方程弹性波模拟.石油地球物理勘探,2012,47(1):19-26.Tang Xiaoping,Bai Chaoying and Liu Kuanhou.Ela-stic wavefield simulation using separated equations through pseudo-spectral method.OGP,2012,47(1):19-26.
[19]
Xiao X and Leaney W S.Local vertical seismic profiling (VSP) elastic reverse-time migration and migration resolution:Salt-flank imaging with transmitted P-to-S waves.Geophysics,2010,75(2):S35-S49.
[20]
Wang W,McMechan G A and Zhang Q.Comparison of two algorithms for isotropic elastic P and S vector decomposition.Geophysics,2015,80(4):T147-T160.
[21]
Liu Y and Sen M K.Scalar wave equation modeling with time-space domain dispersion-relation-based staggered-grid finite-difference schemes.Bulletin of the Seismological Society of America,2011,101(1):141-159.
[22]
Liu Y.Optimal staggered-grid finite-difference schemes based on least-squares for wave equation modelling.Geophysical Journal International,2014, 197(2):1033-1047.
[23]
刘洋.波动方程时空域有限差分数值解及吸收边界条件研究进展.石油地球物理勘探,2014,49(1):35-46.Liu Yang.Research progress on time-space domain finite-difference numerical solution and absorption boundary conditions of wave equation.OGP,2014,49(1):35-46.
[24]
董良国.弹性波数值模拟中的吸收边界条件.石油地球物理勘探,1999,34(1):46-56.Dong Liangguo.Absorptive boundary condition in elastic-wave numerical modeling.OGP,1999,34(1):46-56.
[25]
Du Q,Guo C F and Zhao Q.Vector-based elastic reverse time migration based on scalar imaging condition.Geophysics,2017,82(2):S111-S127.
[26]
张伟,石颖.矢量分离纵横波场的弹性波逆时偏移.地球物理学进展,2017,32(4):1728-1734.Zhang Wei and Shi Ying.Elastic reverse time migration based on vector decomposition of P-and S-wavefields.Progress in Geophysics,2017,32(4):1728-1734.
[27]
王鹏飞,何兵寿.基于行波分离的三维弹性波矢量场点积互相关成像条件.石油地球物理勘探,2017,52(3):477-483.Wang Pengfei and He Bingshou. Vector field dot product cross-correlation imaging based on 3D elastic wave separation.OGP,2017,52(3):477-483.
[28]
Balch A H and Erdemir C.Sign-change correction for prestack migration of P-S converted wave reflections.Geophysical Prospecting,1994,42(6):637-663.