Conjugate gradient and cross-correlation based least-square reverse time migration and its application
Sun Xiao-Dong1,2, Ge Zhong-Hui1, and Li Zhen-Chun1
1. China University of Petroleum (Hua Dong), Qingdao 266580, China.
2. Laboratory for Marine Mineral Resource, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266071, China.
Abstract Although conventional reverse time migration can be perfectly applied to structural imaging it lacks the capability of enabling detailed delineation of a lithological reservoir due to irregular illumination. To obtain reliable reflectivity of the subsurface it is necessary to solve the imaging problem using inversion. The least-square reverse time migration (LSRTM) (also known as linearized reflectivity inversion) aims to obtain relatively high-resolution amplitude preserving imaging by including the inverse of the Hessian matrix. In practice, the conjugate gradient algorithm is proven to be an efficient iterative method for enabling use of LSRTM. The velocity gradient can be derived from a cross-correlation between observed data and simulated data, making LSRTM independent of wavelet signature and thus more robust in practice. Tests on synthetic and marine data show that LSRTM has good potential for use in reservoir description and four-dimensional (4D) seismic images compared to traditional RTM and Fourier finite difference (FFD) migration. This paper investigates the first order approximation of LSRTM, which is also known as the linear Born approximation. However, for more complex geological structures a higher order approximation should be considered to improve imaging quality.
This research is sponsored by The National Natural Science Fund (No. 41574098) and Sinopec Geophysical Key Laboratory Open Fund (No. wtyjy-wx2016-04-2).
Cite this article:
. Conjugate gradient and cross-correlation based least-square reverse time migration and its application[J]. APPLIED GEOPHYSICS, 2017, 14(3): 381-386.
[1]
Dai,W., Boonyasiriwat, C., and Schuster, Ge. T.,2010, 3D Multi-source Least-squares Reverse Time Migration: 80th Ann. Internat. Mtg, Soc. Expl. Geophys., SEG, Expanded Abstracts, 3120−3124.
[2]
Dai, W., Huang, Y. S., and Schuster, G. T., 2013, Least-squares reverse time migration of marine data with frequency-selection encoding: 83th Ann. Internat. Mtg, Soc. Expl. Geophys., SEG, Expanded Abstracts, 3231−3236.
[3]
Dai,W., and Schuster,G. T., 2013, Plane-wave least-squares reverse-time migration:Geophysics,78(4),165−177.
[4]
Dai,W., and Schuster, J., 2009, Least-squares migration of simultaneous sources data with a deblurring filter: 79th Ann. Internat. Mtg, Soc. Expl. Geophys., SEG, Expanded Abstracts, 2990−2994.
[5]
Dai,Wei., Wang, X., and Schuster, G. T., 2011,Least-squares migration of multisource data with a deblurring filter:Geophysics,76(5), 135−146.
[6]
Du, Q. Z., Guo, C. F., Zhao, Q., Gong, X. F., Wang, C. X., and Li, X. Y., 2017, Vector-based elastic reverse time migration based on scalar imaging condition: Geohysics, 82(2), 111−127.
Dutta,G., and Schuster, G. T., 2014, Attenuation compensation for least-squares reverse time migration using the viscoacoustic-wave equation: Geophysics,79(6), 251−262.
[9]
Erlangga,Y. A., Vuik,K., Oosterlee, K., Plessix, R. E., and Mulder, W. A.,2004, A robust iterative solver for the two-way wave equation based on a complex shifted-Laplace preconditioner: 74th Ann. Internat. Mtg, Soc. Expl. Geophys., SEG, Expanded Abstracts, 1897−1900.
[10]
He, R., You, J. C., Liu, B., Wang, Y. C., Deng, S. C., and Zhang, F. Q., 2017, High-order generalized screen propagator migration based on particle swarm optimization: Applied Geophysics, 14(1), 64−72.
[11]
He, Y. Y., Hu, T. Y., He, C., and Tan, Y. Y., 2016, P-wave attenuation anisotropy in TI media and its application in fracture parameters inversion: Applied Geophysics, 13(4), 649− 656.
[12]
Li, C., Huang, J. P., Li, Z. C., and Wang, R. R., 2017, Preconditioned prestack plane-wave least squares reverse time migration with singular spectrum constraint: Applied Geophysics, 14(1), 73−86.
[13]
Ren, H R, Huang, G. H, Wang, H. Z, et al., 2013, A research on the Hessian operator in seismic inversion imaging: Chinese J. Geophys, 56(7), 2429−2436.
[14]
Ren, H. R.,Wang, H. Z., and Huang, G. H., 2012, Theoretical analysis and comparison of seismic wave inversion and imaging methods:Lithologic Reservoirs, 24(5), 12−18.
[15]
Ren, H. R.,Wang, H. Z., andHuang, G, H., 2012, Analysis of the basic problem of seismic wave inversion:Lithologic Reservoirs, 24(6), 1−9.
Schuster,G. T., Dai,W., Zhan, G., and Boonyasiriwat, C., 2010, Theory of Multisource Crosstalk Reduction by Phase-Encoded Statics: 80th Ann. Internat. Mtg, Soc. Expl. Geophys., SEG, Expanded Abstracts, 3110−3114.
[18]
Wang, Y. B., Zheng, Y. K., Xue,Q. F., Chang,X., Tong, W., and Luo,Y., 2017,Reverse time migration of multiples: Reducing migration artifacts using the wavefield decomposition imaging condition:Geophysics,82(4), 307−314.
Zhang, D. L.,Gerard, T., and Schuster, G. T., 2014, Least-squares reverse time migration of multiples: Geophysics, 79(1), 11−21.
[21]
Zhang,Y., Duan, L., and Xie, Y., 2013, A stable and practical implementation of least-squares reverse time migration: 83th Ann. Internat. Mtg, Soc. Expl. Geophys., SEG, Expanded Abstracts, 3716−3720
[22]
Zong, Z. Y., Yin, X. Y., and Li, K., 2016, Joint AVO inversion in the time and frequency domain with Bayesian interference: Applied Geophysics, 13(4), 631−640.
2017, Vol. 14 
