不依赖子波、基于包络的FWI初始模型建立方法研究
|
摘要
地震全波形反演(FWI)从理论走向实际面临着诸多难题,其中之一就是需要一个较高精度的初始模型,另一个难题就是需要一个较为精确的震源子波,初始模型和震源子波的准确程度严重影响着全波形反演的最终结果.为此,本文提出了不依赖子波、基于包络的FWI初始模型建立的方法,建立了相应的目标函数,推导出了反演的梯度,给出了伴随震源的表达式,理论上分析了不依赖子波FWI的可行性.在数值试验中,讨论了参考道的选取方式,通过分析归一化目标函数收敛速率,认为近偏移距参考道优于远偏移距参考道,在地震数据含干扰噪音时,平均道作为参考道要优于最小偏移距参考道.通过包络、包络对数、包络平方三种目标函数反演结果的比较,发现包络对数目标函数对深层的反演效果最好.通过不同子波的试验进一步验证了本方法的正确性.
Full waveform inversion(FWI)is facing many difficulties in practice.One of the issues is how to construct a relatively precise initial model.Another one is the treatment of source wavelet.Both seriously affect the final inversion result.To tackle this problem,we propose a source-independent envelope-based FWI approach to build an initial model for the traditional FWI.On one hand,the cross-convolved wavefields,including the convolution of observed wavefields with a reference trace from the modeled wavefields and the convolution of modeled wavefields with a reference trace from the observed wavefields,can eliminate source wavelet effect in the misfit function of waveform inversion.On the other hand,the envelope of wavefields contains abundant low-frequency components,which are crucial for initial model inversion.Therefore,we use the envelope of cross-convolved wavefields to establish misfit functions,the feasibility of which is proved by theoretical analysis.The expressions of the corresponding gradient and adjoint source are derived.Furthermore,we use a Marmousi synthetic data example to test our method.The proposed method is tested by using different reference traces,misfit functions andsource wavelets.(1)The normalized misfit of long offset reference trace converges more slowly than that of short-offset reference trace for noise-free data and the average trace as reference trace is better than the minimum-offset reference trace for noisy data.(2)The logarithmic envelope is the best and the squared envelope is the worst misfit function to invert the deep part of the model.The initial model built by our method choosing the logarithmic misfit function and minimum-offset reference trace is much better than that of the traditional envelope-based FWI when the source wavelet is incorrect.(3)Using other kinds of source wavelets,such as 20 Hz Ricker wavelet without low frequency below 5 Hz and 12 Hz Gauss wavelet,both of which are quite different from the true 8Hz Ricker wavelet,the inversion result of our method can provide agood initial model for the traditional FWI.In addition,the initial model can be used in seismic depth migration,and the imaging result is close to the RTM result with true velocity model and much better than that with gradient initial model.In conclusion,first,we propose a source-independent envelope-based FWI approach to build an initial model for the traditional FWI without any prior information when the source wavelet is unknown and even the data lack low frequency.The initial model built by our method can also be used as the velocity model for pre-stack depth migration.Second,the short-offset reference trace is better than long-offset reference trace for noise-free data and average trace as reference trace is better than the minimum-offset reference trace for noisy data.Third,the logarithmic envelope misfit function is the optimal misfit function.Finally,by changing the source wavelet,it is further proved that our method is source-independent.
Full waveform inversion(FWI)is facing many difficulties in practice.One of the issues is how to construct a relatively precise initial model.Another one is the treatment of source wavelet.Both seriously affect the final inversion result.To tackle this problem,we propose a source-independent envelope-based FWI approach to build an initial model for the traditional FWI.On one hand,the cross-convolved wavefields,including the convolution of observed wavefields with a reference trace from the modeled wavefields and the convolution of modeled wavefields with a reference trace from the observed wavefields,can eliminate source wavelet effect in the misfit function of waveform inversion.On the other hand,the envelope of wavefields contains abundant low-frequency components,which are crucial for initial model inversion.Therefore,we use the envelope of cross-convolved wavefields to establish misfit functions,the feasibility of which is proved by theoretical analysis.The expressions of the corresponding gradient and adjoint source are derived.Furthermore,we use a Marmousi synthetic data example to test our method.The proposed method is tested by using different reference traces,misfit functions andsource wavelets.(1)The normalized misfit of long offset reference trace converges more slowly than that of short-offset reference trace for noise-free data and the average trace as reference trace is better than the minimum-offset reference trace for noisy data.(2)The logarithmic envelope is the best and the squared envelope is the worst misfit function to invert the deep part of the model.The initial model built by our method choosing the logarithmic misfit function and minimum-offset reference trace is much better than that of the traditional envelope-based FWI when the source wavelet is incorrect.(3)Using other kinds of source wavelets,such as 20 Hz Ricker wavelet without low frequency below 5 Hz and 12 Hz Gauss wavelet,both of which are quite different from the true 8Hz Ricker wavelet,the inversion result of our method can provide agood initial model for the traditional FWI.In addition,the initial model can be used in seismic depth migration,and the imaging result is close to the RTM result with true velocity model and much better than that with gradient initial model.In conclusion,first,we propose a source-independent envelope-based FWI approach to build an initial model for the traditional FWI without any prior information when the source wavelet is unknown and even the data lack low frequency.The initial model built by our method can also be used as the velocity model for pre-stack depth migration.Second,the short-offset reference trace is better than long-offset reference trace for noise-free data and average trace as reference trace is better than the minimum-offset reference trace for noisy data.Third,the logarithmic envelope misfit function is the optimal misfit function.Finally,by changing the source wavelet,it is further proved that our method is source-independent.
引文
BozdagˇE,Trampert J,Tromp J.2011.Misfit functions for fullwaveform inversion based on instantaneous phase and envelopemeasurements.Geophysical Journal International,185(2):845-870.
Brossier R,Operto S,Virieux J.2009.Seismic imaging of complexonshore structures by 2D elastic frequency-domain full-waveforminversion.Geophysics,74(6):WCC105-WCC118.
Bunks C,Saleck F M,Zaleski S,et al.1995.Multiscale seismicwaveform inversion.Geophysics,60(5):1457-1473.
Chi B X,Dong L G,Liu Y Z.2013.Full waveform inversion basedon envelope objective function.∥75th Annual InternationalConference and Exhibition,EAGE,Extended Abstracts.
Chi B X,Dong L G,Liu Y Z.2014.Full waveform inversion methodusing envelope objective function without low frequency data.Journal of Applied Geophysics,109:36-46.
Choi Y,Shin C,Min D J,et al.2005.Efficient calculation of thesteepest descent direction for source-independent seismic waveforminversion:An amplitude approach.Journal of ComputationalPhysics,208(2):455-468.
Choi Y,Alkhalifah T.2011.Source-independent time-domainwaveform inversion using convolved wavefields:Application tothe encoded multisource waveform inversion.Geophysics,76(5):R125-R134.
Dong L G,Ma Z T,Cao J Z,et al.2000.A staggered-grid highorder difference method of one-order elasticwave equation.Chinese Journal of Geophysics(in Chinese),43(3):411-419,doi:10.3321/j.issn:0001-5733.2000.03.015.
Lee K H,Kim H J.2003.Source-independent full-waveform
inversion of seismic data.Geophysics,68(6):2010-2015.
Liu B L,Dong L G,Chi B X.2014.Method of Source-independentFull Waveform Inversion[Ph.D.thesis](in Chinese).Shanghai:Tongji University.
Pratt R G,Shin C,Hick G J.1998.Gauss-Newton and full Newtonmethods in frequency-space seismic waveform inversion.GeophysicalJournal International,133(2):341-362.
Pratt R G.1999.Seismic waveform inversion in the frequencydomain,Part 1:Theory and verification in a physical scalemodel.Geophysics,64(3):888-901.
Sirgue L,Pratt R G.2004.Efficient waveform inversion andimaging:A strategy for selecting temporal frequencies.Geophysics,69(1):231-248.
Song Z M,Williamson P R,Pratt R G.1995.Frequency-domainacoustic-wave modeling and inversion of crosshole data:Part IIInversion method,synthetic experiments and real-data results.Geophysics,60(3):796-809.
Tarantola A.1984.Inversion of seismic reflection data in theacoustic approximation.Geophysics,49(8):1259-1266.
Virieux J,Operto S.2009.An overview of full-waveform inversionin exploration geophysics.Geophysics,74(6):WCC1-WCC26.
Wu R S,Luo J R,Wu B Y.2014.Seismic envelope inversion andmodulation signal model.Geophysics,79(3):WA13-WA24.
Zhou B,Greenhalgh S A.2003.Crosshole seismic inversion withnormalized full-waveform amplitude data.Geophysics,68(4):1320-1330.
董良国,马在田,曹景忠等.2000.一阶弹性波方程交错网格高阶差分解法.地球物理学报,43(3):411-419,doi:10.3321/j.issn:0001-5733.2000.03.015.
刘宝龙,董良国,迟本鑫.2014.不依赖于子波的全波形反演方法研究[硕士论文].上海:同济大学.
Brossier R,Operto S,Virieux J.2009.Seismic imaging of complexonshore structures by 2D elastic frequency-domain full-waveforminversion.Geophysics,74(6):WCC105-WCC118.
Bunks C,Saleck F M,Zaleski S,et al.1995.Multiscale seismicwaveform inversion.Geophysics,60(5):1457-1473.
Chi B X,Dong L G,Liu Y Z.2013.Full waveform inversion basedon envelope objective function.∥75th Annual InternationalConference and Exhibition,EAGE,Extended Abstracts.
Chi B X,Dong L G,Liu Y Z.2014.Full waveform inversion methodusing envelope objective function without low frequency data.Journal of Applied Geophysics,109:36-46.
Choi Y,Shin C,Min D J,et al.2005.Efficient calculation of thesteepest descent direction for source-independent seismic waveforminversion:An amplitude approach.Journal of ComputationalPhysics,208(2):455-468.
Choi Y,Alkhalifah T.2011.Source-independent time-domainwaveform inversion using convolved wavefields:Application tothe encoded multisource waveform inversion.Geophysics,76(5):R125-R134.
Dong L G,Ma Z T,Cao J Z,et al.2000.A staggered-grid highorder difference method of one-order elasticwave equation.Chinese Journal of Geophysics(in Chinese),43(3):411-419,doi:10.3321/j.issn:0001-5733.2000.03.015.
Lee K H,Kim H J.2003.Source-independent full-waveform
inversion of seismic data.Geophysics,68(6):2010-2015.
Liu B L,Dong L G,Chi B X.2014.Method of Source-independentFull Waveform Inversion[Ph.D.thesis](in Chinese).Shanghai:Tongji University.
Pratt R G,Shin C,Hick G J.1998.Gauss-Newton and full Newtonmethods in frequency-space seismic waveform inversion.GeophysicalJournal International,133(2):341-362.
Pratt R G.1999.Seismic waveform inversion in the frequencydomain,Part 1:Theory and verification in a physical scalemodel.Geophysics,64(3):888-901.
Sirgue L,Pratt R G.2004.Efficient waveform inversion andimaging:A strategy for selecting temporal frequencies.Geophysics,69(1):231-248.
Song Z M,Williamson P R,Pratt R G.1995.Frequency-domainacoustic-wave modeling and inversion of crosshole data:Part IIInversion method,synthetic experiments and real-data results.Geophysics,60(3):796-809.
Tarantola A.1984.Inversion of seismic reflection data in theacoustic approximation.Geophysics,49(8):1259-1266.
Virieux J,Operto S.2009.An overview of full-waveform inversionin exploration geophysics.Geophysics,74(6):WCC1-WCC26.
Wu R S,Luo J R,Wu B Y.2014.Seismic envelope inversion andmodulation signal model.Geophysics,79(3):WA13-WA24.
Zhou B,Greenhalgh S A.2003.Crosshole seismic inversion withnormalized full-waveform amplitude data.Geophysics,68(4):1320-1330.
董良国,马在田,曹景忠等.2000.一阶弹性波方程交错网格高阶差分解法.地球物理学报,43(3):411-419,doi:10.3321/j.issn:0001-5733.2000.03.015.
刘宝龙,董良国,迟本鑫.2014.不依赖于子波的全波形反演方法研究[硕士论文].上海:同济大学.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心