Sequentially ordered single-frequency 2-D acoustic waveform inversion in the Laplace-Fourier domain

被引:38
作者
Shin, Changsoo [1 ]
Koo, Nam-Hyung [2 ]
Cha, Young Ho [1 ]
Park, Keun-Pil [2 ]
机构
[1] Seoul Natl Univ, Dept Energy Resources Engn, Seoul 151742, South Korea
[2] Korea Inst Geosci & Mineral Resources, Petr & Marine Res Div, Taejon 305350, South Korea
关键词
Inverse theory; Seismic tomography; Computational seismology; Acoustic properties; SEISMIC-REFLECTION DATA; ELASTIC INVERSION; GAUSS-NEWTON; PART; TOMOGRAPHY; FIELD; MODEL; STRATEGY; TIME;
D O I
10.1111/j.1365-246X.2010.04540.x
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
P>In the conventional frequency-domain waveform inversion, either multifrequency simultaneous inversion or sequential single-frequency inversion has been implemented. However, most conventional frequency-domain waveform inversion methods fail to recover background velocity when low-frequency information is missing. Recently, new waveform inversion techniques in the Laplace and Laplace-Fourier domain have been proposed to recover background velocity structure from data with insufficient low-frequency information. In such techniques, however, all frequencies are inverted simultaneously, and this requires large computational resources and long computation times. In this paper, we propose a sequentially ordered single-frequency 2-D acoustic waveform inversion using the logarithmic objective function in the Laplace-Fourier domain. Our algorithm sequentially inverts single-frequency data in the Laplace-Fourier domain, thus reducing computational resources. Unlike most conventional waveform inversion methods requiring an initial velocity model close to the true model, we propose a one-step waveform inversion method in seeking to find a final velocity structure from the simple initial model through a hybrid combination of the Laplace domain inversion and the Fourier domain inversion. We adopt and evaluate the multiloop algorithm by modifying the double-loop algorithm commonly used in the conventional frequency-domain waveform inversion. Using the multiloop algorithm repeating loop over frequencies, the quality of the inversion results can be improved and the decision problem of the number of iterations for each frequency can be overcome effectively. Because the sequential order of the Laplace-Fourier frequencies in a 2-D plane should be assigned for inverting Laplace-Fourier frequency data consecutively, we present three different sequential orders of Laplace-Fourier frequencies while considering the multiscale and layer-stripping approach, and we compare the inversion results from the numerical experiments. We applied the sequentially ordered single-frequency 2-D acoustic waveform inversion in the full Laplace-Fourier domain to the synthetic seismic data produced from complex structure model and field data. A realistic model could be recovered in an efficient and robust manner, even using the two-layer homogeneous velocity model as an initial model. The inverted velocity model from the field data was validated by examining the migrated image from the pre-stack depth migration and the flattening of the common-image gathers or by comparing the synthetic shot gather with the real shot gather. The proposed one-step waveform inversion algorithm can be easily extended to the sequential inversion of 3-D acoustic or elastic data in the full Laplace-Fourier domain.
引用
收藏
页码:935 / 950
页数:16
相关论文
共 30 条
[1]  
AMINZADEH F, 1997, SEG EAGE 3D MODELING, V1
[2]   Velocity model-building by 3D frequency-domain, full-waveform inversion of wide-aperture seismic data [J].
Ben-Hadj-Ali, Hafedh ;
Operto, Stephane ;
Virieux, Jean .
GEOPHYSICS, 2008, 73 (05) :VE101-VE117
[3]  
BILLETTE FJ, 2005, P 67 M MADR SPAIN JU
[4]   Applying waveform inversion to wide-angle seismic surveys [J].
Bleibinhaus, Florian ;
Lester, Ryan W. ;
Hole, John. A. .
TECTONOPHYSICS, 2009, 472 (1-4) :238-248
[5]   Full waveform tomography for lithospheric imaging: results from a blind test in a realistic crustal model [J].
Brenders, A. J. ;
Pratt, R. G. .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2007, 168 (01) :133-151
[6]   MULTISCALE SEISMIC WAVE-FORM INVERSION [J].
BUNKS, C ;
SALECK, FM ;
ZALESKI, S ;
CHAVENT, G .
GEOPHYSICS, 1995, 60 (05) :1457-1473
[7]  
HU W, 2007, P 69 M LOND UK JUN 1
[8]   2-D traveltime and waveform inversion for improved seismic imaging: Naga Thrust and Fold Belt, India [J].
Jaiswal, Priyank ;
Zelt, Colin A. ;
Bally, Albert W. ;
Dasgupta, Rahul .
GEOPHYSICAL JOURNAL INTERNATIONAL, 2008, 173 (02) :642-658
[9]   Comparison of scaling methods for waveform inversion [J].
Jang, Ugeun ;
Min, Dong-Joo ;
Shin, Changsoo .
GEOPHYSICAL PROSPECTING, 2009, 57 (01) :49-59
[10]  
Lailly P., 1983, Conference on inverse scattering-Theory and application, P206