Regularized least-squares reverse time migration with prior model
Li Zhenchun1, Li Chuang1, Huang Jianping1, Wang Rongrong2
1. School of Geosiences, China University of Petroleum (East China), Qingdao, Shandong 266580, China;
2. College of Information and Control Engineering, China University of Petroleum (East China), Qingdao, Shandong 266580, China
Abstract:Least-squares reverse time migration (LSRTM) converges slowly, and sometimes drops into local extremum because of the ill-posedness of the inversion problem. On the other hand, the influence of irregular geometry and the absorption of underground layers result in some blank areas of illumination, in which the structures cannot be imaged by LSRTM. To solve these problems, the paper presents the theory of regularized least-squares reverse time migration (RLSRTM) with prior model. The prior model is constructed from logging data and incorporated into LSRTM as the regularization constraint. And dynamic regularization parameters and preconditioned regularization-term gradients are proposed to ensure better constraint. Based on numerical tests on a sparse Marmousi model, the following observation are obtained: 1LSRTM can suppress the migration artifacts and compensate energy in deep part compared with reverse time migration (RTM), but the compensation to uneven illumination is limited and the structures in blank illumination areas cannot be imaged; 2RLSRTM with preconditioned regularization-term gradient can further compensate energy of structures with poor illumination, produce more clear images of the boundary of anticlines and other layers in deep part, and even recover some information in blank illumination areas; 3RLSRTM without preconditioned regularization-term gradient produces images with some blurry boundaries and false structures. Therefore, RLSRTM with preconditioned regularization-term gradient can accelerate the convergence, improve the resolution and amplitude preservation of the images, ensure the stability of the inversion, and has reduced sensitivity to low signal-to-noise ratio shot data compared with LSRTM and RLSRTM without preconditioned regularization-term gradient.
Huang Jianping,Li Zhenchun,Kong Xue et al.A study of least-squares migration imaging method for fractured-type carbonate resivior.Chinese Journal of Geophysics,2013, 56(5):1716-1725.
[3]
Kuehl H,Sachhi M D.Split-step WKBJ least-square migration/inversion of incomplete data.5th SEGJ International Symposium Imaging Technology, 2001.
[4]
Kuehl H,Sacchi M D.Least-squares wave-equation migration for AVP/AVA inversion.Geophysics,2003,68(1): 262-273.
Zhou Huamin,Chen Shengchang,Ren Haoran et al. One-way wave equation least-squares migration based on illumination compensation. Chinese Journal of Geophysics, 2014, 57(8): 2644-2655.
[6]
Li C,Huang J P,Li Z C.Plane-wave least square reverse time migration for rugged topography.SEG Technical Program Expanded Abstracts,2014,33: 3742-3746.
[7]
Huang J P,Li C,Wang R R et al.Plane-wave least square reverse time migration for rugged topography.Journal of Earth Science,2015, 26(4):471-480.
[8]
Li C,Huang J P,Li Z C.Application of plane-wave least square migration in fault block reservoirs:a case study.76th EAGE Conference & Exhibition, 2014.
Huang Jianping,Cao Xiaoli,Li Zhenchun et al.Application of least-square reverse time migration in the high resolution imaging of near-surface.OGP,2014, 49(1):107-112.
[10]
Dai W,Fowler P,Schuster G T.Multi-source least-squares reverse time migration.Geophysical Prospecting,2012,60(4): 681-695.
Huang Jianping, Li Chuang, Li Qingyang et al. Least-square reverse time migration with static plane-wave encoding. Chinese Journal of Geophysics,2015,58 (6):2046-2056.
[13]
Dai W, Schuster G T. Plane-wave least-squares reverse-time migration.Geophysics, 2013, 78(4): 5165-5177.