Efficient 3D elastic full-waveform inversion using wavefield reconstruction methods

ABSTRACTIn reverse time migration (RTM) or full-waveform inversion (FWI), forward and reverse time propagating wavefields are crosscorrelated in time to form either the image condition in RTM or the misfit gradient in FWI. The crosscorrelation condition requires both fields to be available at the same time instants. For large-scale 3D problems, it is not possible, in practice, to store snapshots of the wavefields during forward modeling due to extreme storage requirements. We have developed an approximate wavefield reconstruction method that uses particle velocity field recordings on the boundaries to reconstruct the forward wavefields during the computation of the reverse time wavefields. The method is computationally effective and requires less storage than similar methods. We have compared the reconstruction method to a boundary reconstruction method that uses particle velocity and stress fields at the boundaries and to the optimal checkpointing method. We have tested the methods on a 2D vertical trans...

[1]  Andreas Fichtner,et al.  The adjoint method in seismology – I. Theory , 2006 .

[2]  Changsoo Shin,et al.  Efficient calculation of a partial-derivative wavefield using reciprocity for seismic imaging and inversion , 2001 .

[3]  J. Virieux.  P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[4]  Robin M. Weiss,et al.  Solving 3D anisotropic elastic wave equations on parallel GPU devices , 2013 .

[5]  Tieyuan Zhu,et al.  Time-reverse modelling of acoustic wave propagation in attenuating media , 2013 .

[6]  Tieyuan Zhu,et al.  Viscoelastic time-reversal imaging , 2015 .

[7]  J. Etgen,et al.  An overview of depth imaging in exploration geophysics , 2009 .

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

[9]  R. Clapp.  Reverse time migration with random boundaries , 2009 .

[10]  William W. Symes,et al.  Computational Strategies For Reverse-time Migration , 2008 .

[11]  Lianjie Huang,et al.  Reducing the computer memory requirement for 3D reverse-time migration with a boundary-wavefield extrapolation method , 2014 .

[12]  L. Amundsen,et al.  Introduction to Petroleum Seismology , 2005 .

[13]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[14]  Andreas Griewank,et al.  Achieving logarithmic growth of temporal and spatial complexity in reverse automatic differentiation , 1992 .

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

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

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

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

[19]  G. McMechan.  MIGRATION BY EXTRAPOLATION OF TIME‐DEPENDENT BOUNDARY VALUES* , 1983 .

[21]  R. Mittet.  Implementation of the Kirchhoff integral for elastic waves in staggered‐grid modeling schemes , 1994 .

[22]  P. Mora.  Nonlinear two-dimensional elastic inversion of multioffset seismic data , 1987 .

[23]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

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

[25]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[26]  Lijian Tan,et al.  Time-reversal checkpointing methods for RTM and FWI , 2012 .

[27]  Olav Holberg,et al.  COMPUTATIONAL ASPECTS OF THE CHOICE OF OPERATOR AND SAMPLING INTERVAL FOR NUMERICAL DIFFERENTIATION IN LARGE-SCALE SIMULATION OF WAVE PHENOMENA* , 1987 .

[28]  Lianjie Huang,et al.  A wave-energy-based precondition approach to full-waveform inversion in the time domain , 2012 .