CMFWI: Coupled Multiscenario Full Waveform Inversion

We present the new method Coupled Multiscenario Full Waveform Inversion (CMFWI) for the solution of the inversion problem arising in seismic applications. As in the case of classical Full Waveform Inversion (FWI), the proposed method is based on seismic reflection signals and it tries to recover the subsoil velocity profile by minimizing a misfit functional between recorded and computed data. CMFWI combines data generated by shooting one source at a time, and shares the effect of this signal with the signals associated with the other sources. Moreover, CMFWI is a new approach with respect to FWI employed with the same sources shot together. We numerically show that CMFWI performs better than FWI. In particular, this comparison focuses on different types of boundary conditions, wave sources, initial guesses of the velocity profile and signal-to-noise ratio.

[1]  P. G. Ciarlet,et al.  Introduction to Numerical Linear Algebra and Optimisation , 1989 .

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

[3]  Yanfei Wang,et al.  Hybrid regularization methods for seismic reflectivity inversion , 2011 .

[4]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[5]  Andrea Walther,et al.  Online Checkpointing for Parallel Adjoint Computation in PDEs: Application to Goal-Oriented Adaptivity and Flow Control , 2006, Euro-Par.

[6]  Andreas Griewank,et al.  Algorithm 799: revolve: an implementation of checkpointing for the reverse or adjoint mode of computational differentiation , 2000, TOMS.

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

[8]  William W. Symes,et al.  The seismic reflection inverse problem , 2009 .

[9]  A. Majda,et al.  Absorbing boundary conditions for the numerical simulation of waves , 1977 .

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

[11]  Kok Lay Teo,et al.  Optimization and control with applications , 2005 .

[12]  Peter M. Bakker,et al.  Stable P-wave modeling for reverse-time migration in tilted TI media , 2011 .

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

[14]  Felix J. Herrmann,et al.  Randomized Full-waveform Inversion: a Dimenstionality-reduction Approach , 2010 .

[15]  Yuzhu Liu,et al.  Improved hybrid iterative optimization method for seismic full waveform inversion , 2013, Applied Geophysics.

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

[17]  Marcos Raydan,et al.  The Barzilai and Borwein Gradient Method for the Large Scale Unconstrained Minimization Problem , 1997, SIAM J. Optim..

[18]  Gilbert Strang,et al.  Computational Science and Engineering , 2007 .

[19]  J. Lions,et al.  Non-homogeneous boundary value problems and applications , 1972 .

[20]  M. Sacchi,et al.  Model parametrization strategies for Newton-based acoustic full waveform inversion , 2018, Journal of Applied Geophysics.

[21]  Gjlles Aubert,et al.  Mathematical problems in image processing , 2001 .

[22]  J. Borwein,et al.  Two-Point Step Size Gradient Methods , 1988 .

[23]  Gerard T. Schuster,et al.  3D Multisource Full-Waveform Inversion Using Dynamic Random Phase Encoding , 2010 .