Least-squares reverse time migration with a hybrid stochastic conjugate gradient
Li Chuang1,2, Huang Jianping1,2, Li Zhenchun1,2, Wang Rongrong3, Sun Miaomiao14
1. School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China; 2. Laboratory for Marine Mineral Resources, Pilot National Laboratory for Marine Science and Technology, Qingdao, Shandong 266580, China; 3. Institute of Electronics, Chinese Academy of Sciences, Beijing 100190, China; 4. Shengli Branch, Petroleum Engineering Geophysical Limited Corporation, SINOPEC, Dongying, Shandong 257088, China
Abstract:To improve the efficiency of the least-squares reverse time migration (LSRTM),a hybrid stochastic conjugate gradient method is derived with the stochastic optimization theory.The correlation-based shot sampling method is proposed to reduce the amount of migrated data set,which also reduces stochastic characteristics of the gradient.After sampling shot data,the stochastic gradient descent iteration and the conjugate gradient iteration are used alternatively to update the migration imaging.The conjugate gradient iteration improves the convergence efficiency and image quality of LSRTM.Numerical tests on the synthetic data of SEG/EAGE rugged topography model and 2D field data from the exploration area M show that the proposed method can improve the computational efficiency of standard LSRTM,and has higher image quality and faster convergence compared with the stochastic gradient descent method.
陈生昌,张博.基于波场垂向和水平方向外推的高陡构造偏移.地球物理学报,2012,55(4):1300-1306.Chen Shengchang,Zhang Bo.Steeply-dipping structures migration based on the wave field vertical and horizontal extrapolation.Chinese Journal of Geophy-sics,2012,55(4):1300-1306.
[3]
Tang Y.Target-oriented wave-equation least-squares migration/inversion with phase-encoded Hessian.Geophysics,2009,74(6):WCA95-WCA107.
[4]
Tarantola A.Inversion of seismic reflection data in the acoustic approximation.Geophysics,1984,49(8):1259-1266.
[5]
Nemeth T,Wu C J,Schuster G T.Least-squares migration of incomplete reflection data.Geophysics,1999,64(1):208-221.
[6]
Kuehl H,Sacchi M D.Least-squares wave-equation migration for AVP/AVA inversion.Geophysics,2003,68(1):262-273.
[7]
周华敏,陈生昌,任浩然等.基于照明补偿的单程波最小二乘偏移.地球物理学报,2014,57(8):2644-2655.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.
[8]
Dai W,Fowler P,Schuster G T.Multi-source least-squares reverse time migration.Geophysical Prospecting,2012,60(4):681-695.
[9]
Dutta G and Schuster G T.Sparse least squares reverse time migration using seislets.SEG Technical Program Expanded Abstracts,2015,34:4232-4237.
[10]
Hou J and Symes W W.Accelerating extended least-squares migration with weighted conjugate gradient iteration.Geophysics,2016,81(4):S165-S179.
[11]
Li C,Huang J P,Li Z C et al.Regularized least-squares migration of simultaneous-source seismic data with adaptive singular spectrum analysis.Petroleum Science,2017,14(1):61-74.
[12]
Li C,Huang J P,Li Z C et al.Plane-wave least-squares reverse time migration with a preconditioned stocha-stic conjugate gradient method.Geophysics,2018,83(1):S33-S46.
[13]
Romero L A,Ghiglia D C,Ober C C et al.Phase encoding of shot records in prestack migration.Geophy-sics,2000,65(2):426-436.
[14]
Hu J,Wang H,Fang Z et al.Efficient amplitude encoding least-squares reverse time migration using cosine basis.Geophysical Prospecting,2016,64(6):1483-1497.
[15]
Schuster G T,X Wang,Y Huang et al.Theory of multisource crosstalk reduction by phase-encoded statics.Geophysical Journal International,2011,184(3):1289-1303.
[16]
Li C,Huang J P and Li Z C.A hybrid encoding method for seismic data migration.Electronic Journal of Geotechnical Engineering,2017,22(11):4343-4357.
[17]
李娜.基于Huber范数的多震源最小二乘逆时偏移.石油地球物理勘探,2017,52(5):941-947.Li Na.Multi-source least-squares reverse time migration based on Huber norm.OGP,2017,52(5):941-947
[18]
李振春,郭振波,田坤.黏声介质最小平方逆时偏移.地球物理学报,2014,57(1):214-228.Li Zhenchun,Guo Zhenbo,Tian Kun.Least-squares reverse time migration in visco-acoustic media.Chinese Journal of Geophysics,2014,57(1):214-228.
[19]
Dai W and Schuster G T.Plane-wave least-squares reverse-time migration.Geophysics,2013,78(4):5165-5177.
[20]
Huang J P,Li C,Wang R R et al.Plane-wave least-squares reverse time migration for rugged topography.Journal of Earth Science,2015,26(4):471-480.
[21]
Huang Y,Schuster G T.Multisource least-squares migration of marine streamer and land data with fre-quency-division encoding.Geophysical Prospecting,2012,60(4):663-680.
[22]
黄建平,孙郧松,李振春等.一种基于分频编码的最小二乘裂步偏移方法.石油地球物理勘探,2014,49(4):702-707.Huang Jianping,Sun Yunsong,Li Zhenchun et al.Least-squares split-step migration based on frequency-division encoding.OGP,2014,49(4):702-707.
[23]
Leeuwen T V,Aravkin A Y and Herrmann F J.Seismic waveform inversion by stochastic optimization.International Journal of Geophysics,2011,1-18.
[24]
Leeuwen T V and Herrmann F J.Fast waveform inversion without source-encoding.Geophysical Prospecting,2013,61(S1):10-19.
[25]
Ha W and Shin C.Efficient Laplace-domain full waveform inversion using a cyclic shot subsampling me-thod.Geophysics,2013,78(2):R37-R46.
[26]
Moghaddam P P,Keers H,Herrmann F J et al.A new optimization approach for source-encoding full-wave-form inversion.Geophysics,2013,78(3):R125-R132.
[27]
李闯,黄建平,李振春等.预条件最小二乘逆时偏移方法.石油地球物理勘探,2016,51(3):513-520.Li Chuang,Huang Jianping,Li Zhenchun et al.Preconditioned least-squares reverse time migration.OGP,2016,51(3):513-520.
[28]
李振春,李闯,黄建平等.基于先验模型约束的最小二乘逆时偏移方法.石油地球物理勘探,2016,51(4):738-744.Li Zhenchun,Li Chuang,Huang Jianping et al.Regularized least-squares reverse time migration with prior mode.OGP,2016,51(4):738-744.