A 3D full waveform inversion method based on reconstructed source wavelet with wave equation
Liang Zhanyuan1, Wu Guochen1, Wang Yumei2
1. School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China; 2. Research Institute of Geophysics, Shengli Oilfield Branch Co., SINOPEC, Dongying, Shandong 257022, China
Abstract:We propose in this paper a 3D full waveform inversion method based on reconstructed source wavelet with wave equation.We resolve the wave equation with direct waves obtained from near offsets of prestack shot gathers as boundary conditions and take time series at the source point as source wavelet.This method uses the principle of seismic source wave field reconstruction to simulate the inverse process of the wave propagation along the surface,and extracts direct waves at near offsets,which avoids the interference such as refraction wave,reflection wave and diving wave.So the method has less dependence on the initial model.For tests,we select the 3D SEG/EAGE overthrust model and field data to analyze the phase difference of 3D seismic wavelet based on the source wavefield reconstruction theory of wave equation.Test results show that there is only phase difference between the reconstructed wavelet and the input wavelet.After the phase adjustment,the reconstructed wavelet is used for the 3D full waveform inversion; the inversion velocity is close to the real one.Results are more reliable for channels.The waveform and phase between observed records and simulated records have also good correspondence,which validates the feasibility and applicability of the proposed method.
梁展源, 吴国忱, 王玉梅. 基于波动方程重建震源子波的三维全波形反演[J]. 石油地球物理勘探, 2017, 52(6): 1200-1207.
Liang Zhanyuan, Wu Guochen, Wang Yumei. A 3D full waveform inversion method based on reconstructed source wavelet with wave equation. Oil Geophysical Prospecting, 2017, 52(6): 1200-1207.
Lailly P.The seismic inverse problem as a sequence of before stack migrations.Conference on Inverse Scattering:Theory and Application.Society for Industrial and Applied Mathematics,Philadelphia,PA,1983,206-220.
[2]
Tarantola A.Inversion of seismic reflection data in the acoustic approximation.Geophysics,1984,49(8):1259-1266.
[3]
王毓玮,董良国,黄超等.降低弹性波全波形反演强烈非线性的分步反演策略.石油地球物理勘探,2016,51(2):288-294.Wang Yuwei,Dong Liangguo,Huang Chao et al.A multi-step strategy for mitigating severe nonlinearity in elastic full-waveform inversion.OGP,2016,51(2):288-294.
[4]
李志晔,李振春,刘玉金等.基于波场分离的层析波形反演方法.石油地球物理勘探,2016,51(2):295-300.Li Zhiye,Li Zhenchun,Liu Yujin et al.Tomographic waveform inversion with wave-field decomposition.OGP,2016,51(2):295-300.
[5]
Virieux J,Operto S.An overview of full-waveform inversion in exploration geophysics.Geophysics,2009,74(6):WCC1-WCC26.
[6]
付继有,李振春,杨国权等.声介质波动方程反射旅行时反演方法.石油地球物理勘探,2015,50(6):1134-1140.Fu Jiyou,Li Zhenchun,Yang Guoquan et al.Wave equation reflection travel-time inversion in acoustic media.OGP,2015,50(6):1134-1140.
[7]
Maurer H,Greenhalgh S,Latzel S.Frequency and spatial sampling strategies for crosshole seismic waveform spectral inversion experiments.Geophysics,2009,74(6):C79-C89.
[8]
Krebs J R,Anderson J E,Hinkley D et al.Fast full-wavefield seismic inversion using encoded sources.Geophysics,2009,74(6):C177-C188.
[9]
Habashy T,Abubakar A,Pan G et al.Source-receiver compression scheme for full-waveform seismic inversion.Geophysics,2011,76(4):R95-R108.
[10]
Ben-Hadj-Ali H,Operto S,Virieux J.An efficient frequency-domain full waveform inversion method using simultaneous encoded sources.Geophysics,2011,76(4):R109-R124.
[11]
Li X,Aravkin A,van Leeuwen T et al.Fast randomized full-waveform inversion with compressive sensing.Geophysics,2012,77(3):A13-A17.
[12]
Gauthier O,Virieux J,Tarantola A.Two-dimensional nonlinear inversion of seismic waveforms:Numerical results.Geophysics,1986,51(7):1387-1403.
[13]
Feng Bo,Wang Huazhong,Tian Lixin et al.A strategy for source wavefield reconstruction in reverse time migration.SEG Technical Program Expanded Abstracts,2011,20:3164-3168.
[14]
陈乐寿,王光锷.大地电磁测深法.北京:地质出版社,1984.
[15]
杨培杰,印兴耀.地震子波提取方法综述.石油地球物理勘探,2008,43(1):123-128.Yang Peijie,Yin Xingyao.Summary of seismic wavelet pick-up.OGP,2008,43(1):123-128.
[16]
Lee K H,Kim H J.Source-independent full-waveform inversion of seismic data.Geophysics,2003,68(6):2010-2015.
[17]
Xu K,Greenhalgh S A,Wang M Y.Comparison of source-independent methods of elastic waveform inversion.Geophysics,2006,71(6):R91-R100.
[18]
Choi Y,Min D J.Source-independent elastic waveform inversion using a logarithmic wavefield.Journal of Applied Geophysics,2012,76:13-22.
[19]
Zhou B,Greenhalgh S A.Crosshole seismic inversion with normalized full-waveform amplitude data.Geophysics,2003,68(4):1320-1330.
[20]
Choi Y,Alkhalifah T.Source-independent time-domain waveform inversion using convolved wavefields:Application to the encoded multisource waveform inversion.Geophysics,2011,76(5):R125-R134.
[21]
Song Z M,Williamson P R,Pratt R G.Frequency-domain acoustic-wave modeling and inversion of crosshole data:Part Ⅱ-Inversion method,synthetic experiments and real-data results.Geophysics,1995,60(3):796-809.
[22]
敖瑞德,董良国,迟本鑫.不依赖子波、基于包络的FWI初始模型建立方法研究.地球物理学报,2015,58(6):1998-2010.Ao Ruide,Dong Liangguo,Chi Benxin.Source-independent envelope-based FWI to build an initial model.Chinese Journal of Geophysics,2015,58(6):1998-2010.
[23]
Zhou H,Amundsen L,Zhang G.Fundamental issues in full waveform inversion.SEG Technical Program Expanded Abstracts,2012,31.
[24]
Warner M,Ratcliffe A,Nangoo T et al.Anisotropic 3D full-waveform inversion.Geophysics,2013,78(2):R59-R80.
[25]
Ma Y,Hale D.Quasi-Newton full-waveform inversion with a projected Hessian matrix.Geophysics,2012,77(5):R207-R216.