Comparison of Frequency-Selection Strategies for 2D Frequency-Domain Acoustic Waveform Inversion

Several frequency-selection strategies have been used to obtain global minimum solutions in waveform inversion. One strategy, called the discretization method, is to discretize frequencies with a large sampling interval to minimize redundancy in wavenumber information. Another method, the grouping method, groups frequencies with redundancy in wavenumber information. The grouping method can be carried out in two ways. With the first method, the minimum frequency is fixed and the maximum frequency is gradually extended upward (i.e., the overlap-grouping method). Under the second method, frequencies are not overlapped across the groups and waveform inversion proceeds from lower to higher frequency groups (i.e., the individual-grouping method). In this study, we compare these three frequency-selection strategies using both synthetic and real data examples based on logarithmic waveform inversion. Numerical examples for synthetic and real field data demonstrate that the three frequency-selection methods provide solutions closer to the global minimum compared to solutions resulting from simultaneously performed waveform inversion, and that the individual-grouping method yields slightly better resolution for the velocity models than the other methods, particularly for the deeper part. These results may imply that using either too small or too large data sets at every stage slightly deteriorates inversion results, and that grouping data in appropriately sized aggregations improves inversion results.

[1]  M. Toksöz,et al.  Diffraction tomography and multisource holography applied to seismic imaging , 1987 .

[2]  Z. M. Song,et al.  Frequency-domain acoustic-wave modeling and inversion of crosshole data; Part II, Inversion method, synthetic experiments and real-data results , 1995 .

[3]  Cengiz Esmersoy The backpropagated field approach to multidimensional velocity inversion , 1986 .

[4]  Changsoo Shin,et al.  Comparison of waveform inversion, part 1: conventional wavefield vs logarithmic wavefield , 2007 .

[5]  Gerard T. Schuster,et al.  An efficient multiscale method for time-domain waveform tomography , 2009 .

[6]  K. Marfurt Accuracy of finite-difference and finite-element modeling of the scalar and elastic wave equations , 1984 .

[7]  R. Gerhard Pratt,et al.  Efficient waveform inversion and imaging: A strategy for selecting temporal frequencies , 2004 .

[8]  R. Shipp,et al.  Two-dimensional full wavefield inversion of wide-aperture marine seismic streamer data , 2002 .

[9]  A. Devaney Geophysical Diffraction Tomography , 1984, IEEE Transactions on Geoscience and Remote Sensing.

[10]  R. Pratt Seismic waveform inversion in the frequency domain; Part 1, Theory and verification in a physical scale model , 1999 .

[11]  Roelof Versteeg,et al.  The Marmousi experience; velocity model determination on a synthetic complex data set , 1994 .

[12]  C. Bunks,et al.  Multiscale seismic waveform inversion , 1995 .

[13]  Changsoo Shin,et al.  Waveform inversion using a logarithmic wavefield , 2006 .

[14]  Bernard C. Levy,et al.  Variable Background Born Inversion by Wavefield Backpropagation , 1988 .

[15]  C. Shin,et al.  Waveform inversion in the Laplace domain , 2008 .

[16]  Changsoo Shin,et al.  Waveform inversion using a back-propagation algorithm and a Huber function norm , 2009 .

[17]  B. Ursin,et al.  Seismic velocity analysis in the scattering‐angle/azimuth domain , 2003 .

[18]  A. Tarantola Inversion of seismic reflection data in the acoustic approximation , 1984 .

[19]  C. Shin,et al.  Improved amplitude preservation for prestack depth migration by inverse scattering theory , 2001 .

[20]  S. Kim,et al.  Multifrontal Solver Combined with Graph Partitioners , 1999 .

[21]  A. J. Berkhout Multidimensional linearized inversion and seismic migration , 1984 .