时域全波形反演影响因素分析及井数据约束探索
|
摘要
以声波时域全波形反演原理为基础,分析了数据完备性、数据噪声、正问题(包括正问题表征、震源子波、边界及频散等)、初始模型、正则化以及多尺度等因素对时域全波形反演的影响,并提出了相应的实用化建议;探索了加入井数据约束的反演策略,给出了相应的井约束全波形反演误差泛函与梯度计算公式以及CPU并行实现流程,利用复杂地震地质模型分析验证了方法的正确性.
Based on the theory of acoustic full waveform inversion(FWI) in time domain,we have analyzed the effects of factors such as data completness,data noise,modeling problems(modeling problem characterization,the source wavelet,boundary and dispersion),initial model,regularization and multi-scale problems on FWI in time domain,proposed the practical suggestions.And we have explored the inversion strategy with well-data constrained,researched the formula of error function and velocity gradient and implementation process of CPU parallel for well-constrained FWI in time domain.Finally,the correctness of this method is verified by the example of complex seismic geologic model.
Based on the theory of acoustic full waveform inversion(FWI) in time domain,we have analyzed the effects of factors such as data completness,data noise,modeling problems(modeling problem characterization,the source wavelet,boundary and dispersion),initial model,regularization and multi-scale problems on FWI in time domain,proposed the practical suggestions.And we have explored the inversion strategy with well-data constrained,researched the formula of error function and velocity gradient and implementation process of CPU parallel for well-constrained FWI in time domain.Finally,the correctness of this method is verified by the example of complex seismic geologic model.
引文
Bian A F,Yu W H,Zhou H W.2010.Progress in the frequencydomain full waveform inversion method[J].Progress in Geophysics(in Chinese),25(3):982-993,doi:10.3969/j.issn.1004-2903.2010.03.037.
Biondi B,Almomin A.2013.Tomographic full waveform inversion(TFWI)by extending the velocity model along the time-lag axis[C].//SEG Technical Program Expanded Abstracts,1031-1036.
Bunks C,Saleck F M,Zaleski S,et al.1995.Multiscale seismic waveform inversion[J].Geophysics,60(5):1457-1473.
Chen YR,Li Z C,Qin N,et al.2013.Full waveform inversion with wave equation bi-directional illumination optimization[J].Progress in Geophysics(in Chinese),28(6):3015-3021,doi:10.6038/pg2013 0624.
Debayle E,Sambridge M.2004.Inversion of massive surface wave data sets:model construction and resolution assessment[J].J.Geophys.Res.,109(B2),doi:10.1029/2003JB002652.
Dong L G,Chi B X,Tao J X,et al.2013.Objective function behavior in acoustic full-waveform inversion[J].Chinese Journal of Geophysics(in Chinese),56(10):3445-3460,doi:10.6038/cjg20131020.
Gao F X,Liu C,Feng X,et al.2013.Comparisons and analyses of several optimization methods in the application of frequencydomain full waveform inversion[J].Progress in Geophysics(in Chinese),28(4):2060-2068,doi:10.6038/pg20130450.
Lailly P.1983.The seismic inverse problem as a sequence of before stack migrations[A].//Bednar B J,ed.Conference on Inverse Scattering-Theory and Application[M].Tulsa:Society for Industrial and Applied Mathematics(SIAM),206-220.
Liao J P,Liu H X,Zheng G J,et al.2014.Research on multiscale iterative three dimensional acoustic wave full waveform velocity inversion in frequency-space domain-theoretical model[J].Progress in Geophysics(in Chinese),29(1):204-208,doi:10.6038/pg20140l28.
Liu L,Liu H,Zhang H,et al.2013.Full waveform inversion based on modified quasi-Newton equation[J].Chinese Journal of Geophysics(in Chinese),56(7).2447-2451,doi:10.6038/cjg20130730.
Liu Y N,Liu G F,Zang Z F.2014.Multiscale computation problem in seismic wave propagation and imaging[J].Progress in Geophysics(in Chinese),29(1):122-128,doi:10.6038/pg20140116.
Mora P.1987.Nonlinear two-dimensional elastic inversion of multioffset seismic data[J].Geophysics,52(9):1211-1228.
Mora P.1988.Elastic wave-field inversion of reflection and transmission data[J].Geophysics,53(6):750-759.
Nolet G,Montelli R.2005.Optimal parametrization of tomographic models[J].Geophys.J.Int.,161(2):365-372.
Operto S,Virieux J,Amestoy P,et al.2007.3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver:A feasibility study[J].Geophysics,72(5):SM195-SM211.
Plessix R E.2006.A review of the adjoint-state method for computing the gradient of a functional with geophysical applications[J].Geophys.J.Int.,167(2):495-503.
Ren H R.2011.A study of the seismic wave inversion method for imaging in the acoustic media[Ph.D.thesis].Shanghai:Tongji University.
Ren H R,Huang G H,Wang H Z,et al.2013.A research on the Hessian operator in seismic inversion imaging[J].Chinese Journal of Geophysics(in Chinese),56(7):2429-2436,doi:10.6038/cjg20130728.
Shin C,Min D J.2006.Waveform inversion using a logarithmic wavefield[J].Geophysics,71(3):R31-R42.
Sirgue L,Pratt R G.2004.Efficient waveform inversion and imaging:A strategy for selecting temporal frequencies[J].Geophysics,69(1):231-248.
Tang Y X,Lee S W,Baumstein A,et al.2013.Tomographically enhanced full wavefield inversion[C].//SEG Technical Program Expanded Abstracts,1037-1041.
Tape C,Liu Q Y,Maggi A,et al.2010.Seismic tomography of the southern California crust based on spectral-element and adjoint methods[J].Geophys.J.Int.,180(1):433-462.
Tarantola A.1984.Inversion of seismic reflection data in the acoustic approximation[J].Geophysics,49(8):1259-1266.Tarantola A.1987.Inverse problem theory:Methods for data fitting and model parameter estimation[M].2nd ed.University of California:Elsevier.
Wei Z F,Gao H W,Zhang J F.2014.Time-domain full waveform inversion based on an irregular-grid acoustic modeling method[J].Chinese Journal of Geophysics(in Chinese),57(2):586-594,doi:10.6038/cjg20140222.
Yang W Y,Wang X W,Yong X S,et al.2013.The review of seismic Full waveform inversion method[J].Progress in Geophysics(in Chinese),28(2).766-776,doi:10.6038/pg20130225.
卞爱飞,於文辉,周华伟.2010.频率域全波形反演方法研究进展[J].地球物理学进展,25(3):982-993,doi:10.3969/j.issn.1004-2903.2010.03.037.
陈永芮,李振春,秦宁,等.2013.波动方程双向照明优化的全波形反演[J].地球物理学进展,28(6):3015-3021,doi:10.6038/pg20130624.
董良国,迟本鑫,陶纪霞,等.2013.声波全波形反演目标函数性态[J].地球物理学报,56(10):3445-3460,doi:10.6038/cjg20131020.
高凤霞,刘财,冯晅,等.2013.几种优化方法在频率域全波形反演中的应用效果及对比分析研究[J].地球物理学进展,28(4):2060-2068,doi:10.6038/pg20130450.
廖建平,刘和秀,郑桂娟,等.2014.基于多尺度迭代求解三维频率-空间域声波方程的全波形速度反演研究-理论模型[J].地球物理学进展,29(1):204-208,doi:10.6038/pg20140128.
刘璐,刘洪,张衡,等.2013.基于修正拟牛顿公式的全波形反演EJ].地球物理学报,56(7):2447-2451,doi:10.6038/cjg20130730.
刘雅宁,刘国峰,张致付.2014.地震波传播和成像中的多尺度计算问题EJ].地球物理学进展,29(1):122-128,doi:10.6038/pg20140116.
任浩然.2011.声介质地震波反演成像方法研究[博士论文].上海:同济大学.
任浩然,黄光辉,王华忠,等.2013.地震反演成像中的Hessian算子研究[J].地球物理学报,56(7):2429-2436,doi:10.6038/cjg20130728.
魏哲枫,高红伟,张剑锋.2014.基于非规则网格声波正演的时间域全波形反演[J].地球物理学报,57(2);586-594,doi:10.6038/cjg20140222.
杨午阳,王西文,雍学善,等.2013.地震全波形反演方法研究综述[J].地球物理学进展,28(2):766-776,doi:10.6038/pg20130225.
Biondi B,Almomin A.2013.Tomographic full waveform inversion(TFWI)by extending the velocity model along the time-lag axis[C].//SEG Technical Program Expanded Abstracts,1031-1036.
Bunks C,Saleck F M,Zaleski S,et al.1995.Multiscale seismic waveform inversion[J].Geophysics,60(5):1457-1473.
Chen YR,Li Z C,Qin N,et al.2013.Full waveform inversion with wave equation bi-directional illumination optimization[J].Progress in Geophysics(in Chinese),28(6):3015-3021,doi:10.6038/pg2013 0624.
Debayle E,Sambridge M.2004.Inversion of massive surface wave data sets:model construction and resolution assessment[J].J.Geophys.Res.,109(B2),doi:10.1029/2003JB002652.
Dong L G,Chi B X,Tao J X,et al.2013.Objective function behavior in acoustic full-waveform inversion[J].Chinese Journal of Geophysics(in Chinese),56(10):3445-3460,doi:10.6038/cjg20131020.
Gao F X,Liu C,Feng X,et al.2013.Comparisons and analyses of several optimization methods in the application of frequencydomain full waveform inversion[J].Progress in Geophysics(in Chinese),28(4):2060-2068,doi:10.6038/pg20130450.
Lailly P.1983.The seismic inverse problem as a sequence of before stack migrations[A].//Bednar B J,ed.Conference on Inverse Scattering-Theory and Application[M].Tulsa:Society for Industrial and Applied Mathematics(SIAM),206-220.
Liao J P,Liu H X,Zheng G J,et al.2014.Research on multiscale iterative three dimensional acoustic wave full waveform velocity inversion in frequency-space domain-theoretical model[J].Progress in Geophysics(in Chinese),29(1):204-208,doi:10.6038/pg20140l28.
Liu L,Liu H,Zhang H,et al.2013.Full waveform inversion based on modified quasi-Newton equation[J].Chinese Journal of Geophysics(in Chinese),56(7).2447-2451,doi:10.6038/cjg20130730.
Liu Y N,Liu G F,Zang Z F.2014.Multiscale computation problem in seismic wave propagation and imaging[J].Progress in Geophysics(in Chinese),29(1):122-128,doi:10.6038/pg20140116.
Mora P.1987.Nonlinear two-dimensional elastic inversion of multioffset seismic data[J].Geophysics,52(9):1211-1228.
Mora P.1988.Elastic wave-field inversion of reflection and transmission data[J].Geophysics,53(6):750-759.
Nolet G,Montelli R.2005.Optimal parametrization of tomographic models[J].Geophys.J.Int.,161(2):365-372.
Operto S,Virieux J,Amestoy P,et al.2007.3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver:A feasibility study[J].Geophysics,72(5):SM195-SM211.
Plessix R E.2006.A review of the adjoint-state method for computing the gradient of a functional with geophysical applications[J].Geophys.J.Int.,167(2):495-503.
Ren H R.2011.A study of the seismic wave inversion method for imaging in the acoustic media[Ph.D.thesis].Shanghai:Tongji University.
Ren H R,Huang G H,Wang H Z,et al.2013.A research on the Hessian operator in seismic inversion imaging[J].Chinese Journal of Geophysics(in Chinese),56(7):2429-2436,doi:10.6038/cjg20130728.
Shin C,Min D J.2006.Waveform inversion using a logarithmic wavefield[J].Geophysics,71(3):R31-R42.
Sirgue L,Pratt R G.2004.Efficient waveform inversion and imaging:A strategy for selecting temporal frequencies[J].Geophysics,69(1):231-248.
Tang Y X,Lee S W,Baumstein A,et al.2013.Tomographically enhanced full wavefield inversion[C].//SEG Technical Program Expanded Abstracts,1037-1041.
Tape C,Liu Q Y,Maggi A,et al.2010.Seismic tomography of the southern California crust based on spectral-element and adjoint methods[J].Geophys.J.Int.,180(1):433-462.
Tarantola A.1984.Inversion of seismic reflection data in the acoustic approximation[J].Geophysics,49(8):1259-1266.Tarantola A.1987.Inverse problem theory:Methods for data fitting and model parameter estimation[M].2nd ed.University of California:Elsevier.
Wei Z F,Gao H W,Zhang J F.2014.Time-domain full waveform inversion based on an irregular-grid acoustic modeling method[J].Chinese Journal of Geophysics(in Chinese),57(2):586-594,doi:10.6038/cjg20140222.
Yang W Y,Wang X W,Yong X S,et al.2013.The review of seismic Full waveform inversion method[J].Progress in Geophysics(in Chinese),28(2).766-776,doi:10.6038/pg20130225.
卞爱飞,於文辉,周华伟.2010.频率域全波形反演方法研究进展[J].地球物理学进展,25(3):982-993,doi:10.3969/j.issn.1004-2903.2010.03.037.
陈永芮,李振春,秦宁,等.2013.波动方程双向照明优化的全波形反演[J].地球物理学进展,28(6):3015-3021,doi:10.6038/pg20130624.
董良国,迟本鑫,陶纪霞,等.2013.声波全波形反演目标函数性态[J].地球物理学报,56(10):3445-3460,doi:10.6038/cjg20131020.
高凤霞,刘财,冯晅,等.2013.几种优化方法在频率域全波形反演中的应用效果及对比分析研究[J].地球物理学进展,28(4):2060-2068,doi:10.6038/pg20130450.
廖建平,刘和秀,郑桂娟,等.2014.基于多尺度迭代求解三维频率-空间域声波方程的全波形速度反演研究-理论模型[J].地球物理学进展,29(1):204-208,doi:10.6038/pg20140128.
刘璐,刘洪,张衡,等.2013.基于修正拟牛顿公式的全波形反演EJ].地球物理学报,56(7):2447-2451,doi:10.6038/cjg20130730.
刘雅宁,刘国峰,张致付.2014.地震波传播和成像中的多尺度计算问题EJ].地球物理学进展,29(1):122-128,doi:10.6038/pg20140116.
任浩然.2011.声介质地震波反演成像方法研究[博士论文].上海:同济大学.
任浩然,黄光辉,王华忠,等.2013.地震反演成像中的Hessian算子研究[J].地球物理学报,56(7):2429-2436,doi:10.6038/cjg20130728.
魏哲枫,高红伟,张剑锋.2014.基于非规则网格声波正演的时间域全波形反演[J].地球物理学报,57(2);586-594,doi:10.6038/cjg20140222.
杨午阳,王西文,雍学善,等.2013.地震全波形反演方法研究综述[J].地球物理学进展,28(2):766-776,doi:10.6038/pg20130225.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心