An objective function for full-waveform inversion based on frequency-dependent offset-preconditioning

Full-waveform inversion (FWI) is a powerful technique to obtain high-resolution subsurface models, from seismic data. However, FWI is an ill-posed problem, which means that the solution is not unique, and therefore the expert use of the information is required to mitigate the FWI ill-posedness, especially when wide-aperture seismic acquisitions are considered. In this way, we investigate the multiscale frequency-domain FWI by using a weighting operator according to the distances between each source-receiver pair. In this work, we propose a weighting operator that acts on the data misfit as preconditioning of the objective function that depends on the source-receiver distance (offset) and the frequency used during the inversion. The proposed operator emphasizes information from long offsets, especially at low frequencies, and as a consequence improves the update of deep geological structures. To demonstrate the effectiveness of our proposal, we perform numerical simulations on 2D acoustic Marmousi2 case study, which is widely used in seismic imaging tests, considering three different scenarios. In the first two ones, we have used an acquisition geometry with a maximum offset of 4 and 8 km, respectively. In the last one, we have considered all-offsets. The results show that our proposal outperforms similar strategies, for all scenarios, providing more reliable quantitative subsurface models. In fact, our inversion result has the lowest error and the highest similarity to the true model than similar approaches.

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

[2]  S. Ida,et al.  ORBITAL CIRCULARIZATION OF A PLANET ACCRETING DISK GAS: THE FORMATION OF DISTANT JUPITERS IN CIRCULAR ORBITS BASED ON A CORE ACCRETION MODEL , 2014, 1410.8658.

[3]  P. Shive,et al.  An evaluation of Bristow's method for the detection of subsurface cavities , 1990 .

[4]  A. Fichtner Full Seismic Waveform Modelling and Inversion , 2011 .

[5]  L. Sirgue,et al.  The Importance of Low Frequency and Large Offset in Waveform Inversion , 2006 .

[6]  K. Pearson Mathematical Contributions to the Theory of Evolution. III. Regression, Heredity, and Panmixia , 1896 .

[7]  D K Smith,et al.  Numerical Optimization , 2001, J. Oper. Res. Soc..

[8]  Illuminating Santos Basin’s pre-salt with OBN data: Potential and challenges of FWI , 2019, SEG Technical Program Expanded Abstracts 2019.

[9]  Dimitri Komatitsch,et al.  Ultrasonic computed tomography based on full-waveform inversion for bone quantitative imaging , 2017, Physics in medicine and biology.

[10]  Mark Noble,et al.  Robust elastic nonlinear waveform inversion: Application to real data , 1990 .

[11]  G. Corso,et al.  Using fish lateral line sensing to improve seismic acquisition and processing , 2019, PloS one.

[12]  J. Hadamard Sur les problemes aux derive espartielles et leur signification physique , 1902 .

[13]  I. Jones,et al.  High resolution velocity and impedance estimation using refraction and reflection FWI: The Fortuna region, offshore Equatorial Guinea , 2019, SEG Technical Program Expanded Abstracts 2019.

[14]  R. P. Bording,et al.  Applications of seismic travel-time tomography , 1987 .

[15]  A hybrid regularization scheme combining Tikhonov stabilizer with Cauchy function for multiscale frequency domain full-waveform inversion , 2019, SEG Technical Program Expanded Abstracts 2019.

[16]  Jean Virieux,et al.  An overview of full-waveform inversion in exploration geophysics , 2009 .

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

[18]  Tingting Wu,et al.  International Journal of C 2013 Institute for Scientific Numerical Analysis and Modeling Computing and Information an Optimal 9-point Finite Difference Scheme for the Helmholtz Equation with Pml , 2022 .

[19]  Gary Martin,et al.  Marmousi2 An elastic upgrade for Marmousi , 2006 .

[20]  P. Lailly,et al.  Pre-stack inversion of a 1-D medium , 1986, Proceedings of the IEEE.

[21]  F. Mandroux,et al.  Broadband, long-offset, full-azimuth, staggered marine acquisition in the Gulf of Mexico , 2013 .

[22]  E. Haber,et al.  On optimization techniques for solving nonlinear inverse problems , 2000 .

[23]  Jean Virieux,et al.  Crustal seismic imaging from multifold ocean bottom seismometer data by frequency domain full waveform tomography: Application to the eastern Nankai trough , 2006 .

[24]  N. Ricker,et al.  Further developments in the wavelet theory of seismogram structure , 1943 .

[25]  D. J. White,et al.  Two-Dimensional Seismic Refraction Tomography , 1989 .

[26]  G. Corso,et al.  Full-waveform inversion based on Kaniadakis statistics. , 2020, Physical review. E.

[27]  R. Plessix A review of the adjoint-state method for computing the gradient of a functional with geophysical applications , 2006 .

[28]  S. Hanasoge FULL WAVEFORM INVERSION OF SOLAR INTERIOR FLOWS , 2014, 1410.1981.

[29]  Eero P. Simoncelli,et al.  Image quality assessment: from error visibility to structural similarity , 2004, IEEE Transactions on Image Processing.

[30]  Hansruedi Maurer,et al.  Taming the non-linearity problem in GPR full-waveform inversion for high contrast media , 2011 .

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

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

[33]  S. M. Doherty,et al.  Seismic Data Analysis: Processing, Inversion, and Interpretation of Seismic Data , 2000 .

[34]  Pedro Carvalho,et al.  Robust full-waveform inversion using q-statistics , 2020, Physica A: Statistical Mechanics and its Applications.