Elastic full waveform inversion based on frequency-division encoding
SHAO Xiangqi1,2, HE Bingshou1,2, SHI Caiwang3
1. Key Laboratory of Submarine Geoscience and Propecting Technology, Ministry of Education, Ocean University of China, Ocean University of China, Qingdao, Shandong 266100, China; 2. Laboratory of Marine Mineral Resources Evaluation and Detection, Qingdao, Shandong 266100, China; 3. Department of Earth and Space Science, South University of Science and Technology, Shenzhen, Guangdong 518000, China
Abstract:The calculation amount and efficiency limit the application of full waveform inversion technology.Source coding technology can effectively reduce the calculation amount of full waveform inversion and improve the inversion efficiency.The traditional source coding technology requires that each gun has the same receiver array,and it has some shortcomings in the adaptability of the recording system,so it can not be directly applied to the recording system with rolling array.In this paper,an elastic full wave inversion method based on frequency-division coding is proposed.Its basic idea is:In forward process,multiple shot records can be extended at the same time,but each shot record is assigned with a harmonic source of different frequency,so the wave fields forward extended at the same time do not coincide with each other in the frequency spectrum.In the process of forward and residual backward propagation of the shot point wave field,phase sensitivity detection is used to separate the single-frequency wave field of each gun to realize the complete separation of the wave fields of each gun,and then the waveform inversion gradient is constructed to update the model.Compared with the traditional source encoding method,the simultaneous interpreting method can not only be speeded up in inversion,but also be applied to the recording system with rolling array in a wider field.
Tarantola A.Inversion of seismic reflection data in the acoustic approximation[J].Geophysics,1984,49(8):1259-1266.
[2]
Bunks C,Saleck F M,Zaleski S,et al.Multiscale seismic waveform inversion[J].Geophysics,1995,60(5):1457-1473.
[3]
Pratt R G.Seismic waveform inversion in the frequency domain,Part1:Theory and verification in a physical scale model[J].Geophysics,1999,64(3):888-901.
[4]
Sirgue L,Pratt R G.Efficient waveform inversion and imaging:A strategy for selecting temporal frequencies[J].Geophysics,2004,69(1):231-248.
[5]
胡光辉,李熙盛,郭丽,等.构造约束全波形反演及其海上资料应用[J].石油物探,2018,57(4):592-596.HU Guanghui,LI Xisheng,GUO Li,et al.Structure-constrained full waveform inversion and its application in marine seismic data[J].Geophysical Prospecting for Petroleum,2018,57(4):592-596.
[6]
Pratt R G,Shin C,Hicks.Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion[J].Geophysical Journal International,1998,133(2):341-362.
[7]
Sheen D H,Tuncay K,Baag C E,et al.Time domain Gaussâ-Newton seismic waveform inversion in elastic media[J].Geophysical Jornal of the Royal Astronomical Society,2006,167(3):1373-1384.
Ma Y,Hale D.Quasi-Newton full-waveform inversion with a projected Hessian matrix[J].Geophysics,2012,77(5):R207-R216.
[10]
刘璐,刘洪,张衡,等.基于修正拟牛顿公式的全波形反演[J].地球物理学报,2013,56(7):2447-2451.LIU Lu,LIU Hong,ZHANG Heng,et al.Full waveform inversion based on modified quasi-Newton equation[J].Chinese Journal of Geophysics,2013,56(7):2447-2451.
[11]
韩璇颖,印兴耀,曹丹平,等.基于分段快速模拟退火的零偏VSP全波形反演[J].石油物探,2019,58(1):103-111.HAN Xuanying,YIN Xingyao,CAO Danping,et al.Zero-offset VSP velocity inversion with FWI using segmented fast simulated annealing[J].Geophysical Prospecting for Petroleum,2019,58(1):103-111.
[12]
Mao J,Wu R S,Wang B L.Multiscale full waveform inversion using GPU[C].SEG Technical Program Expanded Abstracts,2012,31:1-7.
[13]
Shin J,Ha W,Jun H,et al.3D Laplace-domain full waveform inversion using a single GPU card[J].Computers & Geosciences,2014,67(6):1-13.
[14]
王拓.基于GPU加速的混合域全波形反演技术研究[D].北京:中国石油大学(北京),2016.WANG Tuo.Method Studty of Hybrid Time-frequency Domain FWI Based GPU Acceleration[D].China University of Petroleum(Beijing),Beijing,2016.
[15]
张光超.基于GPU并行计算的时间域弹性波全波形反演方法研究[D].北京:中国石油大学(北京),2016.ZHANG Guangchao.Time-domain Elastic Full Waveform Inversion with GPU Parallel Computing[D].China University of Petroleum(Beijing),Beijing,2016.
[16]
Romero L A,Ghiglia D C,Ober C C,et al.Phase encoding of shot records in prestack migration[J].Geophysics,2000,65(2):426-436.
[17]
Krebs J R,Anderson J E,Hinkley D,et al.Fast full-wave field seismic inversion using encoded sources[J].Geophysics,2009,74(6):WCC177-WCC188.
[18]
Tao Y,Sen M K.Frequency-domain full waveform inversion with plane-wave data[J].Geophysics,2013,78(1):R13-R23.
[19]
Ben-Hadj-Ali H,Operto S,Virieux J.An efficient frequency-domain full waveform inversion method using simultaneous encoded sources[J].Geophysics,2011,76(4):R109-R124.
[20]
Han M,Han L G,Liu C C,et al.Frequency-domain auto-adapting full waveform inversion with blended source and frequency-group encoding[J].Applied Geophysics,2013,10(1):41-52.
[21]
Díaz E,Guitton A.Fast full waveform inversion with random shot decimation[C].SEG Technical Program Expanded Abstracts,2011,30:2804-2808.
[22]
Ha W,Shin C.Efficient Laplace-domain full waveform inversion using a cyclic shot subsampling method[J].Geophysics,2013,78(2):R37-R46.
[23]
Shi C W,He B S.Multiscale full-waveform inversion based on shot subsampling[J].Applied Geophysics,2018,15(2):261-270.
[24]
Liu C C,Han M,Han L G,et al.Application of principal component analysis for frequency-domain full waveform inversion[C].SEG Technical Program Expanded Abstracts,2012,31:1-5.
[25]
段超然,韩立国.主成分分析波场重构反演与全波形反演联合速度重构[J].石油地球物理勘探,2016,51(6):1134-1140.DUAN Chaoran,HAN Liguo.A joint velocity reconstruction method:principal component analysis based wavefield-reconstruction inversion conbined with full waveform inversion[J].Oil Geophysical Prospecting,2016,51(6):1134-1140.
[26]
Nihei K T,Li X.Frequency response modelling of seismic waves using finite difference time domain with phase sensitive detection (TD-PSD)[J].Geophysical Journal of the Royal Astronomical Society,2007,169(3):1069-1078.
[27]
Huang Y S,Schuster G T.Full waveform inversion with multisource frequency selection of marine streamer data[J].Geophysical Prospecting,2018,66(7):1243-1257.
[28]
史才旺,何兵寿.基于主成分分析和梯度重构的全波形反演[J].石油地球物理勘探,2018,53(1):95-104.SHI Caiwang,HE Bingshou.Full waveform inversion based on principal component analysis and gradient reconstruction[J].Oil Geophysical Prospecting,2018,53(1):95-104.
[29]
Choi Y,Alkhalifah T.Application of multi-source waveform inversion to marine streamer data using the global correlation norm[J].Geophysical Prospecting,2012,60(4):748-758.
[30]
夏冬明,宋鹏,谭军,等.基于归一化局部互相关算法的相位编码全波形反演[J].中国海洋大学学报(自然科学版),2017,47(7):105-111.XIA Dongming,SONG Peng,TAN Jun,et al.The phase-encoded full waveform inversion based on normalized local cross-correlation algorithm[J].Periodical of Ocean University of China,2017,47(7):105-111.
[31]
Huang Y S,Schuster G T.Multisource least-squares migration of marine streamer and land data with frequency-division encoding[J].Geophysical Prospecting,2012,60(4):663-680.
[32]
Zhang Q C,Mao W J,Zhou H.Hybrid-domain simultaneous-source full waveform inversion without crosstalk noise[J].Geophysical Journal International,2018,215(3):1659-1681.
[33]
Choi Y,Alkhalifah T.Source-independent time-domain waveform inversion using convolved wavefields:App-lication to the encoded multisource waveform inversion[J].Geophysics,2011,76(5):R125-R134.
[34]
杨涛,张会星,史才旺.不依赖子波的弹性波混合域全波形反演[J].石油地球物理勘探,2019,54(2):348-355.YANG Tao,ZHANG Huixing,SHI Caiwang.Wavelet-independent elastic wave full waveform inversion in hybrid domain[J].Oil Geophysical Prospecting,2019,54(2):348-355.
[35]
辛天亮,黄建平,解飞,等.基于数据相似性的不依赖子波的频率域全波形反演[J].石油地球物理勘探,2020,55(2):341-350.XIN Tianliang,HUANG Jianping,XIE Fei,et al.Source-independent frequency-domain full wavefrom inversion based on data similarity[J].Oil Geophysical Prospecting,2020,55(2):341-350.