Three-dimensional frequency-domain full-waveform inversion with an iterative solver

With the acquisition of wide-aperture seismic data sets, full-waveform inversion is an attractive method for deriving velocity models. Three-dimensional implementations require an efficient solver for the wave equation. Computing 3D time-harmonic responses with a frequency-domain solver is complicated because a large linear system with negative and positive eigenvalues must be solved. Time-domain schemes are an alternative. Nevertheless, existing frequency-domain iterative solvers with an efficient preconditioner are a viable option when full-waveform inversion is formulated in the frequency domain. An iterative solver with a multigrid preconditioner is competitive because of a high-order spatial discretization. Numerical examples illustrated the efficiency of the iterative solvers. Three dimensional full-waveform inversion was then studied in the context of deep-water ocean-bottom seismometer acquisition. Three dimensional synthetic data inversion results showed the behavior of full-waveform inversion with respect to the initial model and the minimum frequency available in the data set. Results on a 3D real ocean-bottom seismometer data set demonstrated the relevance of full-waveform inversion, especially to image the shallow part of the model.

[1]  Walter I. Futterman,et al.  Dispersive body waves , 1962 .

[2]  G. Chavent,et al.  Identification de la Non-Linearité D'Une équation Parabolique Quasilineaire , 1974 .

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

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

[5]  A. Tarantola,et al.  Two‐dimensional nonlinear inversion of seismic waveforms: Numerical results , 1986 .

[6]  Albert Tarantola,et al.  Linearized inversion of multioffset seismic reflection data in the ω-k domain , 1986 .

[7]  Peter Mora,et al.  Inversion = Migration + Tomography , 1988, Shell Conference.

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

[9]  William W. Symes,et al.  Velocity inversion by differential semblance optimization , 1991 .

[10]  Jorge Nocedal,et al.  A Limited Memory Algorithm for Bound Constrained Optimization , 1995, SIAM J. Sci. Comput..

[11]  Guy Chavent,et al.  Determination of background velocities by multiple migration fitting , 1995 .

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

[13]  R. G. Pratt,et al.  Two‐dimensional velocity models from wide‐angle seismic data by wavefield inversion , 1996 .

[14]  Guy Chavent,et al.  Waveform Inversion of Reflection Seismic Data for Kinematic Parameters by Local Optimization , 1998, SIAM J. Sci. Comput..

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

[16]  G. Chavent,et al.  An optimal true-amplitude least-squares prestack depth-migration operator , 1999 .

[17]  G. Schuster,et al.  Least-squares migration of incomplete reflection data , 1999 .

[18]  Guy Chavent,et al.  Migration‐based traveltime waveform inversion of 2-D simple structures: A synthetic example , 2001 .

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

[20]  Graham. J. Hicks,et al.  Arbitrary source and receiver positioning in finite‐difference schemes using Kaiser windowed sinc functions , 2002 .

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

[22]  Jean Virieux,et al.  Multiscale imaging of complex structures from multifold wide-aperture seismic data by frequency-domain full-waveform tomography: application to a thrust belt , 2004 .

[23]  Jean Virieux,et al.  Quantitative imaging of complex structures from dense wide‐aperture seismic data by multiscale traveltime and waveform inversions: a case study , 2004 .

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

[25]  Cornelis Vuik,et al.  A Novel Multigrid Based Preconditioner For Heterogeneous Helmholtz Problems , 2005, SIAM J. Sci. Comput..

[26]  Cornelis Vuik,et al.  A new iterative solver for the time-harmonic wave equation , 2006 .

[27]  S. Operto,et al.  3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study , 2007 .

[28]  R. G. Pratt,et al.  Full waveform tomography for lithospheric imaging: results from a blind test in a realistic crustal model , 2007 .

[29]  Cornelis Vuik,et al.  A parallel multigrid-based preconditioner for the 3D heterogeneous high-frequency Helmholtz equation , 2007, J. Comput. Phys..

[30]  René-Édouard Plessix,et al.  A Helmholtz iterative solver for 3D seismic-imaging problems , 2007 .

[31]  William W. Symes,et al.  Reverse time migration with optimal checkpointing , 2007 .

[32]  Xiaoye S. Li,et al.  Frequency response modelling of seismic waves using finite difference time domain with phase sensitive detection (TD–PSD) , 2007 .

[33]  William W. Symes,et al.  Migration velocity analysis and waveform inversion , 2008 .

[34]  René-Édouard Plessix,et al.  Resistivity imaging with controlled-source electromagnetic data: depth and data weighting , 2008 .

[35]  Wim A. Mulder,et al.  Exploring some issues in acoustic full waveform inversion , 2008 .

[36]  Reinhard Nabben,et al.  Multilevel Projection-Based Nested Krylov Iteration for Boundary Value Problems , 2008, SIAM J. Sci. Comput..

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

[38]  Use of refraction, reflection, and wave-equation-based tomography for imaging beneath shallow gas: A Trinidad field data example , 2008 .

[39]  Jean Virieux,et al.  Velocity model building by 3D frequency-domain, full-waveform inversion of wide-aperture seismic data , 2008 .