Full waveform inversion of time-domain acoustic wave based on first-order approximate instantaneous frequency
XU Yipeng1, REN Zhiming1, LI Zhenchun1, LIU Chang1, HE Zilin1, CHEN Jinmao2
1. SWPI, School of Geosciences, China University of Petroleum(East China), Qingdao, Shandong 266580, China; 2. Shengli Branch of Sinopec Geophysical Corporation, Dongying, Shandong 257000, China
Abstract:A deep velocity model with high accuracy can be obtained by full waveform inversion (FWI) in theory.However,it is difficult to find underground conditions in actual exploration.This is because that,on the one hand,full waveform inversion is strongly dependent on the initial velocity model;and on the other hand,deep signals are weak,and their contribution to the objective function is small.This paper proposes time-domain full waveform inversion based on first-order approximate instantaneous frequency.First,according to the instantaneous phase formula, the objective function and the adjoint source term formula of the first-order approximate instantaneous frequency are derived; Then the initial velocity model for full waveform inversion is established by making full use of the advantage of the instantaneous frequency to highlight low frequency information and deep weak signals.Finally,a gradient attenuation factor is added to further improve the effect of deep inversion and prevent the interference between deep and shallow layers.The feasibility and validity of the method have been proved by model and real data.
任志明.声波和弹性波波动方程有限差分正反演方法研究[D].北京:中国石油大学(北京),2016.REN Zhiming.Research on Finite-Difference Mode-ling and Inversion Methods Based on Acoustic and Elastic Wave Equations[D].China University of Petroleum(Beijing),Beijing,2016.
[2]
罗静蕊,吴如山,高静怀.地震包络反演对局部极小值的抑制特性[J].地球物理学报,2016,59(7):2510-2518.LUO Jingrui,WU Rushan,GAO Jinghuai.Local minima reduction of seismic envelope inversion[J].Chinese Journal of Geophysics,2016,59(7):2510-2518.
[3]
Bozda E,Trampert J,Tromp J.Misfit functions for full waveform inversion based on instantaneous phase and envelope measurements[J].Geophysical Journal International,2011,185(2):845-870.
[4]
Luo J R,Wu R S and Gao F C.Time domain full waveform inversion using instantaneous phase with damping[C].SEG Technical Program Expended Abstracts,2018,37:1472-1476.
[5]
罗静蕊,曹玉玲.应用包络反演为弹性全波形反演构建初始模型[J].油气地球物理,2017,15(2):60-63.LUO Jingrui,CAO Yuling.Applying envelope inversion to build initial model for elastic full waveform inversion[J].Pretroleum Geophysics,2017,15(2):60-63.
[6]
董良国,迟本鑫,陶纪霞,等.声波全波形反演目标函数性态[J].地球物理学报,2013,56(10):3445-3460.DONG Liangguo,CHI Benxin,TAO Jixia,et al.Ob-jective function behavior in acoustic full-waveform inversion[J].Chinese Journal of Geophysics,2013,56(10):3445-3460.
[7]
张海勇.瞬时频率的一种估计方法[J].系统工程与电子技术,2002,20(9):5-6.ZHANG Haiyong.A method for estimating instantaneous frequency[J].Systems Engineering and Electronics,2002,20(9):5-6.
[8]
钟佑明,秦树人,汤宝平.Hilbert-Huang变换中的理论研究[J].振动与冲击,2002,12(4):15-19.ZHONG Youming,QIN Shuren,TANG Baoping.Theoretical research on Hilbert-Huang transformation[J].Journal of Vibration and Shock,2002,12(4):15-19.
[9]
Kosloff D,Baysal E.Forward modeling by a Fourier method[J].Geophysics,1982,47(10):1402-1412.
[10]
Reshef M,Kosloff D,Edwards M,et al.Three-dimension acoustic modeling by the Fourier method[J].Geophysics,1988,53(9):1175-1183.
[11]
程乾生.希尔伯特变换与信号的包络、瞬时相位和瞬时频率[J].石油地球物理勘探,1979,14(3):1-14.CHENG Qiansheng.Hilbert transform and envelope,instantaneous phase and instantaneous frequency of signals[J].Oil Geophysical Prospecting,1979,14(3):1-14.
[12]
贺振华.地震波复数道分析的计算方法[J].石油物探,1980,19(4):64-80.HE Zhenhua.Calculation method of seismic complex trace analysis[J].Geophysical Prospecting for Petroleum,1980,19(4):64-80.
Bartolo L D,Dors C,Mansur W J,et al.A new family of finite-difference schemes to solve the heterogeneous acoustic wave equation[J].Geophysics,2012,77(5):T187-T199.
[15]
Fornberg B.The pseudospectral method:Comparisons with finite differences for the elastic wave equation[J].Geophysics,1987,52(4):483-501.
[16]
成景旺,毛宁波,吕晓春,等.非分裂完全匹配层边界存储时间域全波形反演[J].石油地球物理勘探,2018,53(4):754-764.CHENG Jingwang,MAO Ningbo,LYU Xiaochun,et al.Time domain full waveform inversion with CFS-NPML boundary storage[J].Oil Geophysical Prospecting,2018,53(4):754-764.
[17]
Luo J R,Wu R S.Seismic envelope inversion:reduction of local minima and noise resistance[J].Geophy-sical Prospecting,2015,63(3):597-614.
[18]
Pratt R G.Frequency-domain elastic wave modeling by finite differences:A tool for crosshole seismic imaging[J].Geophysics,1990,55(5):626-632.
[19]
任德玉.地震记录与瞬时频率初探[J].石油地球物理勘探,1980,15(1):7-21.REN Deyu.Preliminary reseacrch on seismic record and instantaneous frequency[J].Oil Geophysical Prospecting,1980,15(1):7-21.
[20]
Wu R S,Luo J Rand Wu B Y.Seismic envelope inversion and modulation signal model[J].Geophysics,2014,79(10):WA13-WA24.
[21]
Alkhalifah T.Scattering-angle based filtering of the waveform inversion gradients[J].Geophysical Journal International,2015,200:363-373.
[22]
史才旺,何兵寿.基于主成分分析和梯度重构的全波形反演[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.