Source-independent frequency-domain full waveform inversion based on data similarity
XIN Tianliang1, HUANG Jianping1, XIE Fei1,2, ZHOU Bing3, LU Zizhuo1
1. School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China; 2. East China Branch, Sinopec, Nanjing, Jiangsu 210011, China; 3. CNOOC(China) Tianjin Company, Tianjin 300452, China
Abstract:In terms of single-frequency data features in frequency domain,we formulate an objective function (SOD) in complex number field,which is similar to the function describing data similarity in real number domain,to mitigate the non-linearity of inversion.This method is practical because amplitude matching between synthetic data and observed data becomes less important.The source wavelet is naturally avoided in the objective function; besides,the reference trace is not required.But for standard trace normalized (STN) full waveform inversion,the convergence may slow down or even suspend if the reference trace is not properly estimated.Model tests show that ① SOD inversion is independent of source wavelet,and it is superior to additional two source-independent algorithms (average trace normalized (ATN) and STN) in the accuracy of inversion; ② SOD inversion is based on data similarity,and it is tolerant to energy discrepancies between different data sets; thus,the inversion is less dependent on the initial velocity model; ③ random noises may account for similar proportions in the numerator and denominator through multiply and add operations.This means that the objective function of SOD is robust; ④ SOD could be integrated with source encoding to effectively improve computational efficiency.
Tarantola A.Inversion of seismic reflection data in the acoustic approximation[J].Geophysics,1984,49(8):1259-1266.
[2]
Pratt R G.Frequency-domain elastic wave modeling by finite differences:A tool for cross-hole seismic imaging[J].Geophysics,1990,55(5):626-632.
[3]
Pratt R G.Inverse theory applied to multi-source cross-hole tomography,Part 2:Elastic wave-equation method[J].Geophysical Prospecting,1990,38(3):311-329.
[4]
Pratt R G,Worthington M H.Inverse theory applied to multi-source cross-hole tomography,Part 1:Acoustic wave-equation method[J].Geophysical Prospecting,1990,38(3):287-310.
[5]
李福元,韦成龙,邓桂林,等.从海洋地震资料直达波提取震源子波[J].石油地球物理勘探,2019,54(3):512-521.LI Fuyuan,WEI Chenglong,DENG Guilin,et al.Source signature extraction from marine seismic direct waves[J].Oil Geophysical Prospecting,2019,54(3):512-521.
[6]
Song Z M,Williamson P R,Pratt R G.Frequency domain acoustic-wave modeling and inversion of crosshole data:Part II-Inversion method,synthetic experiments and real-data results[J].Geophysics,1995,60(3):796-809.
[7]
Pratt R G.Seismic waveform inversion in the frequency domain,Part 1:Theory and verification in a physical scale model[J].Geophysics,1999,64(3):888-901.
[8]
Zhou B,Greenhalgh S A.Crosshole seismic inversion with normalized full-waveform amplitude data[J].Geophysics,2003,68(4):1320-1330.
[9]
Lee K H,Kim H J.Source-independent full-waveform inversion of seismic data[J].Geophysics,2003,68(6):2010-2015.
[10]
Choi Y,Shin C,Min D J,et al.Efficient calculation of the steepest descent direction for source-independent seismic waveform inversion:An amplitude approach[J].Journal of Computational Physics,2005,208(2):455-468.
[11]
Greenhalgh S A,Zhou B.Surface seismic imaging by multi-frequency amplitude inversion[J].Exploration Geophysics,2003,34(4):217-224.
[12]
Xu K,Greenhalgh S A,Wang M Y.Comparison of source-independent methods of elastic waveform inversion[J].Geophysics,2006,71(6):R91-R100.
[13]
Choi Y,Min D J.Source-independent elastic waveform inversion using a logarithmic wavefield[J].Journal of Applied Geophysics,2012,76(10):13-22.
[14]
Choi Y,Alkhalifah T.Source-independent time-domain waveform inversion using convolved wavefields:Application to the encoded multisource waveform inversion[J].Geophysics,2011,76(5):R125-R134.
[15]
敖瑞德,董良国,迟本鑫.不依赖子波、基于包络的FWI初始模型建立方法研究[J].地球物理学报,2015,58(6):1998-2010.AO Ruide,DONG Liangguo,CHi Benxin.Source-independent envelope-based FWI to build an initial model[J].Chinese Journal of Geophysics,2015,58(6):1998-2010.
[16]
Zhang Q,Zhou H,Li Q,et al.Robust source-inde-pendent elastic full-waveform inversion in the time domain[J].Geophysics,2016,81(2):R29-R44.
[17]
杨涛,张会星,史才旺.不依赖子波的弹性波混合域全波形反演[J].石油地球物理勘探,2019,54(2):348-355.YANG Tao,ZHANG Huixing and SHI Caiwang.Wavelet-independent elastic wave full waveform inversion in hybrid domain[J].Oil Geophysical Prospecting,2019,54(2):348-355.
[18]
Gao F,Williamson P,Pratt R G.A new objective function for full waveform inversion:Differential sem-blance in the data domain[C].SEG Technical Program Expanded Abstracts,2014,33:1178-1183.
[19]
Liao J P,Guo Z W,Liu H X,et al.Application of frequency-dependent traveltime tomography to 2D crosswell seismic field data[J].Journal of Environmental & Engineering Geophysics,2017,22(4):421-426.
[20]
Liao J P,Wang H Z,Ma Z T.2-D elastic wave modeling with frequency-space 25-point finite-difference operators[J].Applied Geophysics,2009,6(3):259-266.
[21]
Pratt R G,Shin C,Hick G J.Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion[J].Geophysical Journal International,1998,133(2):341-362.
[22]
Pratt R G,Symes W.Semblance and differential semblance optimisation for waveform tomography:A frequency domain implementation[C].Sub-basalt Imaging Conference Expanded Abstracts,2002,183-184.
[23]
成景旺,吕晓春,顾汉明,等.基于柯西分布的频率域全波形反演[J].石油地球物理勘探,2014,49(5):940-945.CHENG Jingwang,LYU Xiaochun,GU Hanming,et al.Full waveform inversion with Cauchy diatribution in the frequency domain[J].Oil Geophysical Prospecting,2014,49(5):940-945.
[24]
周聪,刘江平,罗银河,等.二维频率域全波场有限差分数值模拟方法[J].石油地球物理勘探,2014,49(2):278-287.ZHOU Cong,LIU Jiangping,LUO Yinhe,et al.2D full-wavefield modeling in frequency domain using finite-difference[J].Oil Geophysical Prospecting,2014,49(2):278-287.
[25]
Jo C H,Shin C,Suh J H.An optimal 9-point,finite difference,frequency-space,2-D scalar wave extrapolator[J].Geophysics,1996,61(2):529-537.
[26]
Godwin J,Sava P.Simultaneous Source Imaging by Amplitude Encoding[R].CWP Report,2010,21-42.
[27]
解飞,黄建平,李振春.一种基于高效迭代解法的频率域全波形反演[J].地球物理学报,2018,61(8):1998-2010.XIE Fei,HUANG Jianping,LI Zhenchun.An efficient full waveform inversion in frequency domain based on wavefield iteration method[J].Chinese Journal of Geophysics,2018,61(8):1998-2010.