Abstract:We propose in this paper a three-dimensional full waveform inversion based on the first-order velocity-stress acoustic wave equation.First based on the time-domain three-dimensional first-order velocity-stress acoustic wave equation,an adjoint equation of the time-domain three-dimensional first-order velocity-stress acoustic wave equation and a corresponding gradient equation of P-wave velocity are derived with the adjoint-state method based on the perturbation theory.Then the staggered grid finite difference method is adopted to obtain the synthetic forward wavefield and the backward adjoint wavefield in the time domain.After that the gradient of the P-wave velocity is calculated.On this basis,the conjugate gradient method is used to update P-wave velocity model.Finally the 3D full waveform inversion method based on the first-order velocity-stress acoustic wave equation is implemented.Because this method is based on the first-order velocity-stress acoustic wave equation,the staggered grid finite difference method can be easily applied,and the synthetic wavefield and adjoint wavefield are directly and quickly extrapolated based on GPU in the time domain.The feasibility and effectiveness of the proposed method are demonstrated by model tests.
段文胜,李飞,黄录忠等.OVT域宽方位层析速度建模与深度域成像.石油地球物理勘探,2016,51(3):521-528.Duan Wensheng,Li Fei,Huang Luzhong et al.Wide azimuth tomographic velocity model building and depth imaging in the OVT domain.OGP,2016,51(3):521-528.
[2]
杨午阳,张厚柱,撒利明等.逆时偏移关键问题探讨.石油地球物理勘探,2016,51(6):1251-1262.Yang Wuyang,Zhang Houzhu,Sa Liming et al.Some issues about reverse time migration.OGP,2016,51(6):1251-1262.
[3]
魏哲枫,高红伟,张剑锋.基于非规则网格声波正演的时间域全波形反演.地球物理学报,2014,57(2):586-594.Wei Zhefeng,Gao Hongwei and Zhang Jianfeng.Time domain full waveform inversion based on an irregular grid acoustic modeling method.Chinese Journal of Geophysics,2014,57(2):586-594.
[4]
杨午阳,王西文,雍学善等.地震全波形反演方法研究综述.地球物理学进展,2013,28(2):766-776.Yang Wuyang,Wang Xiwen,Yong Xueshan et al.The review of seismic full waveform inversion method.Progress in Geophysics,2013,28(2):766-776.
[5]
石玉梅,张研,姚逢昌等.基于声学全波形反演的油气藏地震成像方法.地球物理学报,2014,57(2):607-617.Shi Yumei,Zhang Yan,Yao Fengchang et al.Metho-dology of seismic imaging for hydrocarbon reservoirs based on acoustic full waveform inversion.Chinese Journal of Geophysics,2014,57(2):607-617.
[6]
Tarantola A.Inversion of seismic reflection data in the acoustic approximation.Geophysics,1984,49(8):1259-1266.
[7]
Mora P.Nonlinear two-dimensional elastic inversion of multioffset seismic data.Geophysics,1987,52(9):1211-1228.
[8]
Warner M,Ratcliffe A,Nangoo T.Anisotropic 3D full-waveform inversion.Geophysics,2013,78(2):R59-R80.
[9]
Plessix R E,Baeten G,Maag J W et al.Full waveform inversion and distance separated simultaneous swee-ping:A study with a land seismic data set.Geophysical Prospecting,2012,60(4):733-747.
[10]
Mei J W,Ahmed S,Searle A et al.Application of full waveform inversion on an Alaska land 3D survey.SEG Technical Program Expanded Abstracts,2014,33:981-986.
[11]
Plessix R E.A review of the adjoint-state method for computing the gradient of a functional with geophysical applications.Geophysical Journal International,2006,167(2):495-503.
[12]
Castellanos C,Etienne V,Hu G H et al.Algorithmic and methodological developments towards full waveform inversion in 3D elastic media.SEG Technical Program Expanded Abstracts,2011,30:2793-2798.
[13]
胡光辉,贾春梅,夏洪瑞等.三维声波全波形反演的实现与验证.石油物探,2013,52(4):417-425.Hu Guanghui,Jia Chunmei,Xia Hongrui et al.Implementation and validation of 3D acoustic full waveform inversion.GPP,2013,52(4):417-425.
[14]
Ben-Hadj-Ali H,Operto S and Virieux J.An efficient frequency-domain full waveform inversion method using simultaneous encoded sources.Geophysics,2011,76(4):R109-R124.
[15]
张广智,孙昌路,潘新朋等.快速共轭梯度法频率域声波全波形反演.石油地球物理勘探,2016,51(4):730-737.Zhang Guangzhi,Sun Changlu,Pan Xinpeng et al.Acoustic full waveform inversion in the frequency domain based on fast conjugate gradient method.OGP,2016,51(4):730-737.
[16]
Shin C,Cha Y H.Waveform inversion in the Laplace domain.Geophysical Journal International,2008,173(3):922-931.
[17]
Virieux J and Operto S.An overview of full-waveform inversion in exploration geophysics.Geophysics,2013,78(4):R151-R166.
[18]
Shin C,Cha Y H.Waveform inversion in the Laplace-Fourier domain.Geophysical Journal International,2009,177(4):1067-1079.
[19]
Vavryuk V and Kühn D.Moment tensor inversion of waveforms:a two-step time-frequency approach.Geophysical Journal International,2012,190(3):1761-1776.
[20]
Kim Y,Shin C,Calandra H et al.An algorithm for 3D acoustic time-Laplace-Fourier-domain hybrid full waveform inversion.Geophysics,2013,78(4):R151-R166.
[21]
Zou K,Langlo L T,Rønholt G et al.Full-waveform inversion in a shallow water environment:A North Sea 3D towed-streamer data example.SEG Technical Program Expanded Abstracts,2013,32:919-923.
[22]
Mothi S,Schwarz K and Zhu H F.Impact of full-azimuth and long-offset acquisition on full waveform inversion in deep water Gulf of Mexico.SEG Technical Program Expanded Abstracts,2013,32:924-928.
Fichtner A.Full seismic waveform modelling and inversion.Springer,Heidelberg/Dordrecht/London/New York,2011.
[26]
李海山,杨午阳,雍学善.二维时间域黏声波全波形反演.石油地球物理勘探,2018,53(1):87-94.Li Haishan,Yang Wuyang,Yong Xueshan.Visco-acoustic wave full waveform inversion in the 2D time domain.OGP,2018,53(1):87-94.
[27]
Shi Z J,Shen J.Step-size estimation for unconstrained optimization methods.Computational & Applied Mathematics,2005,24(3):399-416.
[28]
冯海新,刘洪,孙军等.基于GPU/CPU和震源随机编码技术的混合域全波形反演.石油物探,2017,56(1):107-115.Feng Haixin,Liu Hong,Sun Jun et al.Hybrid domain full waveform inversion based on GPU/CPU and source random coding technique.GPP,2017,56(1):107-115.
[29]
Köhn D.Time Domain 2D Elastic Full Waveform Tomography[D].Kiel University,Germany,2011.
[30]
Choi Y and Alkhalifah T.Source-independent time-domain waveform inversion using convolved wavefields:Application to the encoded multisource waveform inversion.Geophysics,2011,76(5):R125-R134.
[31]
Asnaashari A,Brossier R,Garambois S et al.Regula-rized seismic full waveform inversion with prior model information.Geophysics,2013,78(2):R25-R36.