Multi-scale full waveform inversion based on double diagonal flux correction transport
YANG Ruidong1, HUANG Jianping1, YANG Zhenjie1, MU Xinru1, WANG Zhaozhong2, YANG Shanshan3
1. School of Geosciences, China University of Petroleum (East China), Qingdao, Shandong 266580, China; 2. Institute of Geology and Mineral Resources of Shandong Province, Jinan, Shandong 250100, China; 3. RIPED Tarim Oilfield Company PetroChina, Korla, Xinjiang 841000, China
Abstract:The continuous furthering of oil and gas exploration exposes seismic exploration to increasingly complex issues. Full waveform inversion (FWI) has become a highly promising velocity modeling method due to its high accuracy. Nevertheless, its inversion effect is largely affected by the accuracy of the initial model and the frequency range of the shot record. In actual seismic data processing, the accuracy of the initial velocity model is often not high, despite that an accurate initial velocity model is the key to preventing FWI from falling into local extreme values. Therefore, reducing the depen-dence of FWI on the initial model is vital. For this reason, a multi-scale FWI method based on double-diagonal flux correction is proposed in this paper. Specifically, the two-dimensional time-domain acoustic wave equation is corrected by the flux-corrected transport technology to obtain the corrected wavefields of different frequency components. On this basis, the calculation formulas for FWI gradient and step size are derived, and multi-scale inversion is carried out with the corrected wavefields from low frequency to high frequency. Simulation tests with a modified Marmousi model are conducted to verify the effectiveness of the proposed method. The results show that the proposed method effectively expands the medium- and low-frequency components of the wavefield and reduces the dependence of FWI on the initial model.
杨瑞冬, 黄建平, 杨振杰, 慕鑫茹, 王兆忠, 杨珊珊. 基于双对角通量校正的多尺度全波形反演[J]. 石油地球物理勘探, 2022, 57(5): 1120-1128.
YANG Ruidong, HUANG Jianping, YANG Zhenjie, MU Xinru, WANG Zhaozhong, YANG Shanshan. Multi-scale full waveform inversion based on double diagonal flux correction transport. Oil Geophysical Prospecting, 2022, 57(5): 1120-1128.
TARANTOLA A. Inversion of seismic reflection data in the acoustic approximation[J]. Geophysics,1984,49(8):1259-1266.
[2]
王庆,张建中. 时间域地震全波形反演方法进展[J]. 地球物理学进展,2015,30(6):2797-2806.WANG Qing,ZHANG Jianzhong. Progress in the time-domain full waveform inversion[J]. Progress in Geophysics,2015,30(6):2797-2806.
[3]
卞爱飞,於文辉,周华伟. 频率域全波形反演方法研究进展[J]. 地球物理学进展,2010,25(3):982-993.BIAN Aifei,YU Wenhui,ZHOU Huawei. Progress in the frequency-domain full waveform inversion method[J]. Progress in Geophysics,2010,25(3):982-993.
[4]
MORA P. Nonlinear two-dimensional elastic inversion of multi offset seismic data[J]. Geophysics,1987,52(9):1211-1228.
PRATT R G,SHIN C,HICKS G J. Gauss-Newton and full Newton methods in frequency-space seismic waveform inversion[J]. Geophysical Journal International,1998,133(2):341-362.
[7]
SIRGUE L,PRATT R G. Efficient waveform inversion and imaging:a strategy for selecting temporal frequencies[J]. Geophysics,2004,69(1):231-248.
[8]
CHEN G X,CHEN S C,WU R S. Full waveform inversion in time domain using time-damping filters[C]. SEG Technical Program Expanded Abstracts,2015,34:1451-1455.
[9]
黄建平,崔超,刘梦丽. 基于频率—波数滤波的联合多尺度全波形反演方法及其实现策略[J]. 中国石油大学学报(自然科学版),2018,42(2):50-59.HUANG Jianping,CUI Chao,LIU Mengli. A hybrid multi-scale full waveform inversion method based on frequency-wavenumber filter and its implementation strategies[J]. Journal of China University of Petroleum (Edition of Natural Science),2018,42(2):50-59.
[10]
SHIN C,CHA Y H. Waveform inversion in the Laplace domain[J]. Geophysical Journal International,2008,173(3):922-931.
[11]
SHIN C,HA W. A comparison between the behavior of objective function for waveform inversion in the frequency and Laplace domains[J]. Geophysics,2008,73(5):VE119-VE133.
[12]
WU R S,LUO J R,WU B Y. Ultra-low-frequency information in seismic data and envelope inversion[C]. SEG Technical Program Expanded Abstracts,2013,32:3078-3082.
[13]
胡勇,韩立国,许卓,等. 基于精确震源函数的解调包络多尺度全波形反演[J]. 地球物理学报,2017,60(3):1088-1105.HU Yong,HAN Liguo,XU Zhuo,et al. Demodulation envelope multi-scale full waveform inversion based on precise seismic source function[J]. Chinese Journal of Geophysics,2017,60(3):1088-1105.
[14]
包乾宗,陈俊霓,吴浩. 基于地震数据包络的多尺度全波形反演方法[J]. 石油物探,2018,57(4):584-591.BAO Qianzong,CHEN Junni,WU Hao. Multi-scale full waveform inversion based on logarithmic envelope of seismic data[J]. Geophysical Prospecting for Petroleum,2018,57(4):584-591.
[15]
ZHANG Pan,HAN Liguo,ZHOU Yan,et al. Passive-source multitaper-spectral method based low-frequency data reconstruction for active seismic sources[J]. Applied Geophysics,2015,12(4):585-597.
[16]
GUO Y D,HUANG J P,CUI C,et al. Direct building background velocity field by plane-wave multi-source multi-scale full-waveform inversion[J]. Exploration Geophysics,2020,51(3):327-340.
[17]
BORIS J P,BOOK D L. Flux-corrected transport. Ⅰ:SHASTA,a fluid transport algorithm that works[J]. Journal of Computational Physics,1973,11(1):38-69.
[18]
BOOK D L,BORIS J P,HAIN K. Flux-corrected transport Ⅱ:Generalizations of the method [J]. Journal of Computational Physics,1975,18(3):248-283.
[19]
FEI T,LARNER K. Elimination of numerical dispersion in finite-difference modeling and migration by flux-corrected transport[J]. Geophysics,1995,60(6):1830-1842.
[20]
吴国忱,王华忠. 波场模拟中的数值频散分析与校正策略[J]. 地球物理学进展,2005,20(1):58-65.WU Guochen,WANG Huazhong. Analysis of nume-rical dispersion in wave-field simulation[J]. Progress in Geophysics,2005,20(1):58-65.
[21]
陈可洋. 地震波数值模拟中优化的通量校正传输方法[J]. 内陆地震,2012,26(2):169-179.CHEN Keyang. An optimized flux-corrected transport method in the seismic wave numerical simulation[J]. Inland Earthquake,2012,26(2):169-179.
[22]
赵威,李伟波,桂志先,等. 基于波动方程有限差分数值模拟及FCT消频散分析[J]. 中国锰业,2018,36(2):187-191.ZHAO Wei,LI Weibo,GUI Zhixian,et al. Finite difference numerical simulation of wave equation and its dispersion analysis of FCT elimination[J]. China’s Manganese Industry,2018,36(2):187-191.
[23]
KALITA M,ALKHALIFAH T. Flux-corrected transport for full-waveform inversion[J]. Geophysical Journal International,2019,217(3):2147-2164.
[24]
张彬彬,张军华,吴永亭. 地震数据低频信号保护与拓频方法研究[J]. 地球物理学进展,2019,34(3):1139-1144.ZHANG Binbin,ZHANG Junhua,WU Yongting. Research on protection and extension for seismic low frequencies[J]. Progress in Geophysics,2019,34(3):1139-1144.
[25]
廖建平,刘和秀,郑桂娟,等. 基于多尺度迭代求解三维频率—空间域声波方程的全波形速度反演研究——理论模型[J]. 地球物理学进展,2014,29(1):204-208.LIAO Jianping,LIU Hexiu,ZHENG Guijuan,et al. Research on multiscale iterative three dimensional acoustic wave full waveform velocity inversion in frequency-space domain:theoretical model[J]. Progress in Geophysics,2014,29(1):204-208.
[26]
张军华,张在金,张彬彬,等. 地震低频信号对关键处理环节的影响分析[J]. 石油地球物理勘探,2016,51(1):54-62.ZHANG Junhua,ZHANG Zaijin,ZHANG Binbin,et al. Low frequency signal influences on key seismic data processing procedures[J]. Oil Geophysical Prospecting,2016,51(1):54-62.
[27]
PLESSIX R E. A review of the adjoint-state method for computing the gradient of a functional with geophysical applications[J]. Geophysical Journal International,2006,167(2):495-503.
[28]
邹志辉,张翊孟,卞爱飞,等. 常规检波器低频数据的评价与恢复及其在地震成像中的应用[J]. 石油地球物理勘探,2016,51(5):841-849.ZOU Zhihui,ZHANG Yimeng,BIAN Aifei,et al. Low-frequency evaluation and recovery of conventional geophone data and applications in seismic imaging[J]. Oil Geophysical Prospecting,2016,51(5):841-849.
[29]
芮拥军,邹志辉,王胜阁,等. 基于高低频联测的常规检波器数据低频振幅和相位同时恢复方法[J]. 石油地球物理勘探,2017,52(4):631-643.RUI Yongjun,ZOU Zhihui,WANG Shengge,et al. Low-frequency amplitude and phase simultaneous recovery of conventional geophones data based on joint deployment of high- and low-frequency geophones[J]. Oil Geophysical Prospecting,2017,52(4):631-643.
[30]
廖建平,王华忠,李伟波,等. 三维微分方程滤波法及线性噪声压制[J]. 石油物探,2007,46(6):550-556.LIAO Jianping,WANG Huazhong,LI Weibo,et al. 3D partial differential equation filtering and linear noise attenuation[J]. Geophysical Prospecting for Petroleum,2007,46(6):550-556.