Time-domain full waveform inversion with CFS-NPML boundary storage
Cheng Jingwang1, Mao Ningbo1, Lü Xiaochun2, Chang Suoliang3, Yan Hao4, Zhong Hua1
1. Key Laboratory of Exploration Technology for Oil and Gas Resources, Ministry of Education, Yangtze University, Wuhan, Hubei 430100, China; 2. College of Resources and Environment, North China University of Water Resources and Electric Power, Zhengzhou, Henan 450011, China; 3. College of Mining Engineering, Taiyuan University of Technology, Taiyuan, Shanxi 030024, China; 4. Tianjin Branch, CNOOC, Tianjin 300452, China
Abstract:The forward-modeling wavefield storage or reconstruction is a key question in the time-domain full waveform inversion.The strategy of storage replaced by computation based on conventional splitting perfect matched layer (SPML) boundary conditions can effectively reduce the amount of storage.However,conventional SPML boundary conditions require a large number of variables,a large amount of memory,and complex programming.In this paper,an effective boundary storage strategy based on the complex frequency shifted non-splitting perfect matched layer (CFS-NPML) is proposed to reconstruct the forward-modeling wavefield to realize the time-domain full waveform inversion.According to our study,fewer variables in CFS-NPML are needed; the absorption effect is better than conventional SPML boundary conditions when grid points of the boundary absorption layer are less.So less storage capacity is needed than that of conventional SPML boundary conditions in the wavefield reconstruction process.Numerical experiments show that there is almost no difference between the forward-modeling wavefield and the reconstructed wavefield obtained with this storage strategy.The velocity obtained by the full waveform inversion with this storage strategy is consistent with that obtained by the conventional storage wavefield inversion.The wavefield reconstruction does not significantly increase the computer time consuming thanks to parallel computation introduction.
杨勤勇,胡光辉,王立歆.全波形反演研究现状及发展趋势.石油物探,2014,53(1):77-82.Yang Qinyong,Hu Guanghui,Wang Lixin.Research status and development trend of full waveform inversion.GPP,2014,53(1):78-82.
[2]
Virieux J,Operto S.An overview of full waveform inversion in exploration geophysics.Geophysics,2009,74(6):WCC1-WCC26.
[3]
Tarantola A.Inversion of seismic reflection data in the acoustic approximation.Geophysics,1984,49(8):1259-1266.
[4]
Pratt R G.Seismic waveform inversion in the frequency domain,Part 1:Theory and verification in a physical scale model.Geophysics,1999,64(3):888-901.
[5]
Shin C,Cha Y H.Waveform inversion in the Laplace domain.Geophysical Journal International,2008,173(3):922-931.
[6]
Ben-Hadj-Ali H,Operto S,Virieux J.Velocity model building by 3D frequency-domain,full waveform inversion of wide-aperture seismic data.Geophysics,2008,73(5):VE101-VE117
[7]
李媛媛,李振春,张凯等.基于双级并行的弹性波频率域全波形多尺度反演方法.应用地球物理,2015,12(4):545-554.Li Yuanyuan,Li Zhenchun,Zhang Kai et al.Frequency-domain elastic full-waveform multiscale inversion method based on dual-level parallelism.Applied Geophysics,2015,12(4):545-554.
[8]
Han M,Han L G,Liu C C et al.Frequency-domain auto-adapting full wave-form inversion with blended source and frequency-group encoding.Applied Geophysics,2013,10(1):41-52.
[9]
成景旺,吕晓春,顾汉明等.基于柯西分布的频率域全波形反演.石油地球物理勘探,2014,49(5):940-945.Cheng Jingwang,Lü Xiaochun,Gu Hanming et al.Full waveform inversion with Cauchy distribution in the frequency domain.OGP,2014,49(5):940-945.
[10]
张广智,孙昌路,潘新朋等.快速共轭梯度法频率域全波形反演.石油地球物理勘探,2016,51(4):730-737.Zhang Guangzhi,Sun Changlu,Pan Xinpeng et al.Acoustic full waveform inversion with Cauchy in the frequency domain based on fast conjugate gradient method.OGP,2016,51(4):730-737.
[11]
苗永康.基于L-BFGS算法的时间域全波形反演.石油地球物理勘探,2015,50(3):469-474.Miao Yongkang.Full waveform inversion in time domain based on limited-memory BFGS algorithm.OGP,2015,50(3):469-474.
[12]
白璐,韩立国,张盼等.最小平方滤波时间域全波形反演.石油地球物理勘探,2016,51(4):721-729.Bai Lu,Han Liguo,Zhang Pan et al.Time-domain full waveform inversion based on least square filter.OGP,2016,51(4):721-729.
[13]
Sigure L,Etgen J T,Albertin U.3D frequency-domain waveform inversion using time-domian finite-diffe-rence methods.Extended Abstracts of 77th EAGE Conference & Exhibition,2008,F022
[14]
胡光辉,贾春梅,夏洪瑞等.三维声波全波形反演的实现与验证.石油物探,2013,52(4):417-425.Hu Guanghui,Jia Chunmei,Xia Hongrui et al.Implementation and validation of 3D acoustic full waveform inversion.GPP,2013,52(4):417-425.
[15]
冯海新,刘洪,孙军等.基于GPU/CPU和震源随机编码技术的混合域全波形反演.石油物探,2017,56(1):107-115.Feng Haixin,Liu Hong,Sun Jun et al.Hybrid domain full waveform inversion based on GPU/CPU and source random coding technique.GPP,2017,56(1):107-115.
[16]
Symes W W.Reverse time migration with optimal checkpointing.Geophysics,2007,72(5):SM213-SM221.
[17]
Dussaud E,Symes W W,Williamson P et al.Computational strategies for reverse-time migration.SEG Technical Program Expanded Abstracts,2008,27:2267-2271.
[18]
Clapp R G.Reverse time migration:Saving the boun-daries.SEG Technical Program Expanded Abstracts,2008,27:136-144.
[19]
Clapp R G.Reverse time migration with random boundaries.SEG Technical Program Expanded Abstracts,2009,28:2809-2813.
[20]
王保利,高静怀,陈文超等.地震叠前逆时偏移的有效边界储存策略.地球物理学报,2012,55(7):2412-2421.Wang Baoli,Gao Jinghuai,Chen Wenchao et al.Efficient boundary storage strategies for seismic reverse time migration.Chinese Journal of Geophysics,2012,55(7):2412-2421.
[21]
胡昊,刘伊克,常旭等.逆时偏移计算中的边界处理分析及应用.地球物理学报,2013,56(6):2033-2042.Hu Hao,Liu Yike,Chang Xu et al.Analysis and applicaton of boundary treatment for the computation of reverse-time migration.Chinese Journal of Geophy-sics,2013,56(6):2033-2042.
[22]
Qin Z,Lu M,Zheng X et al.The implementation of an improved NPML absorbing boundary condition in elastic wave modeling.Applied Geophysics,2009,6(2):113-121.
Martin R,Komatitsch D,Gedney S D et al.A high-order time and space formulation of the unsplit perfectly matched layer for the seismic wave equation using auxiliary differential equations (ADE-PML).Compu-ter Modeling in Engineering & Sciences,2010,56(1):17-41.
[25]
Komatitsch D,Martin R.An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation.Geophysics,2007,72(5):SM155-SM167.
[26]
Drossaert F H,Giannopoulos A.A nonsplit complex frequency-shifted PML based on recursive integration for FDTD modeling of elastic waves.Geophysics,2007,72(2):T9-T17.
[27]
Zhang Z,Lin Y,Huang L.Full-waveform inversion in the time domain with an energy-weighted gradient.SEG Technical Program Expanded Abstracts,2011,30:2772-2776.
[28]
Boonyasiriwat C,Vaul V,Routh P et al.An efficient multiscale method for time-domain waveform tomography.Geophysics,2009,74(6):WCC59-WCC68.
杨积忠,刘玉柱,董良国.变密度声波方程多参数全波形反演策略.地球物理学报,2014,57(2):628-643.Yang Jizhong,Liu Yuzhu,Dong Liangguo.A multi-parameter full waveform inversion strategy for acoustic media with variable density.Chinese Journal of Geophysics,2014,57(2):628-643.