Waveform inversion based on adaptive source wavelet extraction and correction for Rayleigh waves
SHAO Guangzhou1, DU Ting1, WU Hua2
1. School of Geological Engineering and Geomatics, Chang'an University, Xi'an, Shaanxi 710054, China; 2. School of Science, Chang'an University, Xi'an, Shaanxi 710064, China
Abstract:Full waveform inversion (FWI) based on Rayleigh waves has made significant advances in recent years. If the receiver coverage is adequate,FWI offers a high resolution in laterally heterogeneous complex media. Source wavelet is one of the key factors affecting the accuracy of waveform inversion. However,the seismic wavelet for field seismic record is unknown,and an improper source wavelet is bound to have a negative effect on inversion results. The conventional solution is extracting and correcting a wavelet with synthesized data. However,wavelet correction is time-consuming because it requires repeated iterations. To solve this problem,this paper proposes a windowed adaptive wavelet extraction method to extract a source wavelet that is in line with the actual situation from the observation records,use it as the initial wavelet for wavelet correction,and thereby achieve the purpose of further optimizing inversion results. Inversion tests on a fault model constructed are conducted with five wavelets,including the actual wavelet,a wrong wavelet,a wrong wavelet combined with wavelet correction,an extracted wavelet,and an extracted wavelet combined with wavelet correction,respectively. The test results show that the wavelet obtained by the proposed method corresponds to the best inversion result,while the wrong wavelet corresponds to the worst one. Although a favorable inversion effect can also be achieved with the extracted wavelet alone,combining wavelet extraction with wavelet correction can save computation time and further improve the accuracy of Rayleigh waveform inversion. More-over,field data application further verifies this conclusion.
SAENGER E H,BOHLEN T. Finite-difference modeling of viscoelastic and anisotropic wave propagation using the rotated staggered grid[J]. Geophysi-cs,2004,69(2):583-591.
[2]
XIA J H,MILLER R D,PARK C B,et al. Comparing shear-wave velocity profiles inverted from multichannel surface wave with borehole measurements[J]. Soil Dynamics and Earthquake Engineering,2002,22(3):181-190.
[3]
LERNER-LAM A L,JORDAN T H. Earth structure from fundamental and higher-mode waveform analysis[J]. Geophysical Journal International,1983,75(3):759-797.
[4]
YOMOGIDA K,AKI K. Amplitude and phase data inversions for phase velocity anomalies in the Pacific Ocean basin[J]. Geophysical Journal International,1987,88(1):161-204.
[5]
SNIEDER R. Large-scale waveform inversions of surface waves for lateral heterogeneity:1.theory and numerical examples[J]. Journal of Geophysical Research—Solid Earth,1988,93(B10):12055-12065.
[6]
DU Z J. Waveform inversion for lateral heterogeneities using multimode surface waves[J]. Geophysical Journal International,2002,149(2):300-312.
[7]
ZENG C,XIA J H,MILLER R D,et al. Feasibility of waveform inversion of Rayleigh waves for shallow shear-wave velocity using a genetic algorithm[J]. Journal of Applied Geophysics,2011,75(4):648-655.
[8]
TRAN K T,MCVAY M. Site characterization using Gauss-Newton inversion of 2-D full seismic waveform in the time domain[J]. Soil Dynamics and Earthquake Engineering,2012,43:16-24.
[9]
GROOS L,SCHÄFER M,FORBRIGER T,et al. The role of attenuation in 2D full-waveform inversion of shallow-seismic body and Rayleigh waves[J]. Geophysics,2014,79(6):R247-R261.
[10]
GROOS L,SCHÄFER M,FORBRIGER T,et al. Application of a complete workflow for 2D elastic full-waveform inversion to recorded shallow-seismic Rayleigh waves[J]. Geophysics,2017,82(2):R109-R117.
[11]
BORISOV D,MODRAK R,GAO F C,et al. 3D elastic full-waveform inversion of surface waves in the presence of irregular topography using an envelope-based misfit function[J]. Geophysics,2018,83(1):R1-R11.
PRATT R G. Seismic waveform inversion in the frequency domain,part 1:theory and verification in a physical scale model[J]. Geophysics,1999,64(3):888-901.
[14]
SUN D,JIAO K,CHENG X,et al. Source wavelet estimation in full waveform inversion[C]. SEG Technical Program Expanded Abstracts,2014,33:1184-1188.
[15]
梁展源,吴国忱,王玉梅. 基于波动方程重建震源子波的三维全波形反演[J]. 石油地球物理勘探,2017,52(6):1200-1207.LIANG Zhanyuan,WU Guochen,WANG Yumei. A 3D full waveform inversion method based on reconstructed source wavelet with wave equation[J]. Oil Geophysical Prospecting,2017,52(6):1200-1207.
[16]
辛天亮,黄建平,解飞,等. 基于数据相似性的不依赖子波的频率域全波形反演[J]. 石油地球物理勘探,2020,55(2):341-350.XIN Tianliang,HUANG Jianping,XIE Fei,et al. Source-independent frequency-domain full waveform inversion based on data similarity[J]. Oil Geophysical Prospecting,2020,55(2):341-350.
[17]
PARK Y,PYUN S. Source wavelet estimation using common mid-point gathers[J]. Geosystem Enginee-ring,2015,18(4):199-204.
[18]
BEHURA J,SNIEDER R. Virtual real source:source signature estimation using seismic interferometry[J]. Geophysics,2013,78(5):Q57-Q68.
[19]
DOU S,AJO-FRANKLIN J B. Full-wavefield inversion of surface waves for mapping embedded low-velocity zones in permafrost[J]. Geophysics,2014,79(6):EN107-EN124.
[20]
PAN Y D,XIA J H,GAO L G. Full waveform inversion of Rayleigh waves for near-surface shear-wave velocity[C]. Near-Surface Asia Pacific Conference,2015,253-256.
[21]
杨涛,张会星,史才旺. 不依赖子波的弹性波混合域全波形反演[J]. 石油地球物理勘探,2019,54(2):348-355.YANG Tao,ZHANG Huixing,SHI Caiwang. Wavelet-independent elastic wave full waveform inversion in hybrid domain[J]. Oil Geophysical Prospecting,2019,54(2):348-355.
[22]
SCHUSTER G T. Seismic Inversion (Investigations in Geophysics)[M]. Society of Exploration Geophy-sicists,Tulsa,Okla,2017,357-363.
[23]
WANG Y,MILLER R D,PETERIE S L,et al. Tunnel detection at Yuma Proving Ground,Arizona,USA—part 1:2D full-waveform inversion experiment[J]. Geophysics,2019,84(1):B95-B105.
[24]
SMITH J A,BORISOV D,CUDNEY H,et al. Tunnel detection at Yuma Proving Ground,Arizona,USA—part 2:3D full-waveform inversion experiments[J]. Geophysics,2019,84(1):B107-B120.