Reverse Time Migration of Decoupled qP- and qSV-wave in TTI Media with the Optimal Quadratic Approximation

MU Xinru1,2 HUANG Jianping1,2 LI Zhenchun1,2 GUO Xu1,2 YONG Peng1,2 YUAN Shuangqi1,2

(1.School of Geosciences, China University of Petroleum (East China), Qingdao, Shandong, China 266580)
(2.Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao, Shandong, China 266000)

【Abstract】In the seismic data processing, if the influence of anisotropy and the symmetry inclination of tilted transversely isotropic (TTI) media changes are ignored, the imaging position of targets will be deviated, and even the amplitude imbalance may be caused by the caustic problem. The pseudo-shear-wave interference exists in the forward modeling and reverse time migration (RTM) when using the conventional TI coupled pseudo-acoustic wave equations, and the numerical simulation is unstable when the anisotropy parameter ε < δ. We proposed a new approach to solve the problem. By referring to the experience of predecessors, we used the idea of the optimal quadratic approximation (OQA) to approximate the accurate phase velocity formula in VTI medium, and approximate phase velocity formulas of qP-wave and qSV-wave with relatively high accuracy are obtained respectively. Then, based on the approximate phase velocity formula, the decoupled qP-wave and qSV-wave equations in TTI media were derived. We used a hybrid finite-difference/pseudo-spectral method to solve the wave equations to realize efficient and stable numerical simulation. Finally, high-precision, high-efficiency, and stable TTI media forward modeling and RTM were established. Numerical simulation results show that the proposed method can achieve migration imaging in complex TTI media.

【Keywords】 reverse time migration (RTM); TTI media; pure qP-wave; pure qSV-wave; optimal quadratic approximation (OQA); hybrid finite-difference/pseudo-spectral method;


【Funds】 Special Fund of Taishan Scholars for Youth Talent Project (Tsqn20161016) National Oil and Gas Major Project—“Scattered-wave Seismic Imaging Technology Based on Multiple Scattering Theory” (2016ZX05014-008-005) National Oil and Gas Major Project—“Joint Imaging based on Full-Waveform Inversion and Least-Square Migration for Thin Interbeds” (2016ZX05002-005-002)

Download this article


    [1] Thomsen L. Weak elastic anisotropy [J]. Geophysics, 1986, 51 (10): 1954–1966.

    [2] Wu Guochen. Seismic Wave Propagation and Imaging in Anisotropic Media [M]. Petroleum University Press, Dongyong, Shandong, 2006 (in Chinese).

    [3] Du X, Fletcher R P, and Fowler P J. A new pseudo-acoustic wave equation for VTI media [C]. Extended Abstracts of 70th EAGE Conference & Exhibition, 2008, H033.

    [4] DUAN Xinyi, LI Zhenchun, HUANG Jianping, et al. A prestack Gaussian beam depth migration in common-shot domain for anisotropic media [J]. Geophysical Prospecting for Petroleum, 2014, 53 (5): 579–586 (in Chinese).

    [5] HUANG Jinqiang, LI Zhenchun. Modeling and reverse time migration of pure quasi-P-waves in complex TI media with a low-rank decomposition [J]. Chinese Journal of Geophysics, 2017, 60 (2): 704–721 (in Chinese).

    [6] Crampin S. Anisotropy in exploration seismics [J]. First Break, 1984, doi: 10.3997/1365-2397.1984006.

    [7] Winterstein D F. Velocity anisotropy terminology for geophysicists [J]. Geophysics, 1990, 55 (8): 1070–1088.

    [8] DU Qizhen, QIN Tong. Multicomponent prestack reverse-time migration of elastic waves in transverse isotropic medium [J]. Chinese Journal of Geophysics, 2009, 52 (3): 801–807 (in Chinese).

    [9] Yong P, Huang J P, Li Z C, et al. Elastic-wave reverse-time migration based on decoupled elastic-wave equations and inner-product imaging condition [J]. Journal of Geophysics & Engineering, 2016, 13 (6): 953–963.

    [10] Zhou H, Zhang G, and Bloor R. An anisotropic acoustic wave equation for VTI media [C]. Extended Abstracts of 68th EAGE Conference & Exhibition, 2006, H033.

    [11] Du X, Bancroft J C, and Lines L R. Anisotropic reverse-time migration for tilted TI media [J]. Geophysical Prospecting, 2007, 55 (6): 853–869.

    [12] Duveneck E, Bakker P M. Stable P-wave modeling for reverse-time migration in tilted TI media [J]. Geophysics, 2011, 76 (2): S65–S75.

    [13] WANG Weiguo, XIONG Shuijin, XU Huaning, et al. Reverse-time migration using pseudo-spectral method for tilted TI media [J]. Oil Geophysical Prospecting, 2012, 47 (4): 566–572 (in Chinese).

    [14] CHENG Jiubing, KANG Wei, WANG Tengfei. Description of qP-wave propagation in anisotropic media, Part Ⅰ: Pseudo-pure-mode wave equations [J]. Chinese Journal of Geophysics, 2013, 56 (10): 3474–3486 (in Chinese).

    [15] LIU Wenqing, WANG Xiwen, WANG Yuchao, et al. A regularized qP-wave equation for TTI media and its application to reverse time migration [J]. Chinese Journal of Geophysics, 2016, 59 (3): 1059–1069 (in Chinese).

    [16] LI Zhenchun, YANG Fusen, WANG Xiaodan. High-precision numerical simulation of first-order qP-waves in the transversely isotropic (TI) medium optimized by the LS-RSGFD method [J]. Chinese Journal of Geo-physics, 2016, 59 (4): 1477–1490 (in Chinese).

    [17] HU Shuhua, WANG Yuchao, LIU Wenqing, et al. Pure quasi-P wave stable simulation in complex TTI media [J]. Oil Geophysical Prospecting, 2018, 53 (2): 280–287 (in Chinese).

    [18] Tsvankin I. P-wave signatures and notation for transversely isotropic media: An overview [J]. Geophysics, 1996, 61 (2): 467–483.

    [19] Alkhalifah T. Acoustic approximations for processing in transversely isotropic media [J]. Geophysics, 1998, 63 (2): 623–631.

    [20] Alkhalifah T. An acoustic wave equation for aniso-tropic media [J]. Geophysics, 2000, 65 (4): 1239–1250.

    [21] Klíe H and W Toro. A new acoustic wave equation for modeling in anisotropic media [C]. SEG Technical Program Expanded Abstracts, 2001, 20: 1171–1174.

    [22] Hestholm S. Acoustic VTI modeling using high-order finite differences [C]. SEG Technical Program Expanded Abstracts, 2007, 26: 139–143.

    [23] Duveneck E, Milcik P, Bakker P M, et al. Acoustic VTI wave equations and their application for anisotropic reverse-time migration [C]. SEG Technical Program Expanded Abstracts, 2008, 27: 2186–2190.

    [24] Zhang Y and Zhang H. A stable TTI reverse time migration and its implementation [C]. SEG Technical Program Expanded Abstracts, 2009, 28: 2794–2798.

    [25] 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 (in Chinese).

    [26] ZHANG Yan, WU Guochen. Methods of removing pseudo SV-wave artifacts in TTI media qP-wave reverse-time migration [J]. Chinese Journal of Geophy-sics, 2013, 56 (6): 2065–2076 (in Chinese).

    [27] Zhang H, Zhang G, and Zhang Y. Removing S-wave noise in TTI reverse time migration [C]. SEG Technical Program Expanded Abstracts, 2009, 28: 2849–2853.

    [28] LI Bo, LI Min, LIU Hongwei, et al. Stability of reverse time migration in TTI media [J]. Chinese Journal of Geophysics, 2012, 55 (4): 1366–1375 (in Chinese).

    [29] Chu C, Macy B K, and Anno P D. Approximation of pure acoustic seismic wave propagation in TTI media [J]. Geophysics, 2011, 76 (5): WB97-WB107.

    [30] Zhan G, Pestana R C, and Stoffa P L. Decoupled equations for reverse time migration in tilted transversely isotropic media [J]. Geophysics, 2012, 77 (2): 37–45.

    [31] Song G, Huang R, Tian J, et al. A new qP-wave equation for 2D VTI media [C]. SEG Technical Program Expanded Abstracts, 2016, 35: 3977–3981.

    [32] Mu X R, Huang J P, Huang J Q, et al. A precise si-mulation of pure qP-wave in VTI media [C]. Extended Abstracts of 80th EAGE Conference & Exhibition, 2018, doi: 10.3997/2214-4609.201800970.

    [33] Li X and Zhu H. A finite-difference approach for solving pure quasi-P wave equations in transversely isotropic and orthorhombic media [J]. Geophysics, 2018, 83 (4): C161–C172.

    [34] Zhan G, Pestana R C, and Stoffa P L. An efficient hybrid pseudo spectral/finite difference scheme for solving TTI pure P wave equation [C]. 13rd International Congress of the Brazilian Geophysics Society & EXPOGEF, Rio de Janeior, Brazil, 2013, 1322–1327.

    [35] Kosloff D D, Baysal E. Forward modeling by a Fourier method [J]. Geophysics, 1982, 47 (10): 1402–1412.

    [36] Yang J, Zhang S, Zhu H. Isotropic elastic wavefields decomposition using fast Poisson solvers [C]. SEG Technical Program Expanded Abstracts, 2017, 36: 4716–4720.

    [37] Baysal E, Dan D K, Sherwood J W C. Reverse time migration [J]. Geophysics, 1983, 48 (11): 1514–1524.

    [38] LIU Hongwei, LI Bo, LIU Hong, et al. The algorithmof high order finite difference pre-stack reverse time migration and GPU implementation [J]. Chinese Journal of Geophysics, 2010, 53 (7): 1725–1733 (in Chinese).

    [39] LI Zhenchun, LI Qingyang, HUANG Jianping, et al. A stable and high-precision dual-variable grid forward modeling and reverse time migration method [J].Geophysical Prospecting for Petroleum, 2014, 53 (2): 127–136 (in Chinese).

    [40] Feng B, Wang H Z. Reverse time migration with source wavefield reconstruction strategy [J]. Journal of Geophysics & Engineering, 2011, 9 (1): 69–74.

    [41] WANG Baoli, GAO Jinghuai, CHEN Wenchao, et al. Efficient boundary storage strategies for seismic reverse time migration [J]. Chinese Journal of Geophy-sics, 2012, 55 (7): 2412–2421 (in Chinese).

This Article


CN: 13-1095/TE

Vol 54, No. 06, Pages 1280-1292+1174

December 2019


Article Outline


  • 0 Introduction
  • 1 Methodology
  • 2 Model tests
  • 3 Conclusion
  • References