Reverse Time Migration Using the Pseudospectral Time-Domain Algorithm

We propose a pre-stack reverse time migration (RTM) seismic imaging method using the pseudospectral time-domain (PSTD) algorithm. Traditional pseudospectral method uses the fast Fourier transform (FFT) algorithm to calculate the spatial derivatives, but is limited by the wraparound effect due to the periodicity assumed in the FFT. The PSTD algorithm combines the pseudospectral method with a perfectly matched layer (PML) for acoustic waves. PML is a highly effective absorbing boundary condition that can eliminate the wraparound effect. It enables a wide application of the pseudospectral method to complex models. RTM based on the PSTD algorithm has advantages in the computational efficiency compared to traditional methods such as the second-order and high order finite difference time-domain (FDTD) methods. In this work, we implement the PSTD algorithm for acoustic wave equation based RTM. By applying the PSTD-RTM method to various seismic models and comparing it with RTM based on the eighth-order FDTD method, we find that PSTD-RTM method has better performance and saves more than 50% memory. The method is suitable for parallel computation, and has been accelerated by general purpose graphics processing unit.

[1]  Woojae Seong,et al.  Pseudospectral time-domain method for pressure-release rough surface scattering using a surface transformation and image method , 2009, Numerical Algorithms.

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

[3]  Paul L. Stoffa,et al.  Rapid expansion and pseudo spectral implementation for reverse time migration in VTI media , 2012 .

[4]  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 .

[5]  Tsili Wang,et al.  Finite-difference modeling of elastic wave propagation: A nonsplitting perfectly matched layer approach , 2003 .

[6]  Gang Zhao,et al.  MULTIDOMAIN PSEUDOSPECTRAL TIME-DOMAIN (PSTD) METHOD FOR ACOUSTIC WAVES IN LOSSY MEDIA , 2004 .

[7]  Olivier Bou Matar,et al.  Pseudo-Spectral simulation of 1D nonlinear propagation in heterogeneous elastic media , 2015 .

[8]  Jianguo Zhao,et al.  An unsplit complex-frequency-shifted PML based on matched Z-transform for FDTD modelling of seismic wave equations , 2012 .

[9]  Wei Zhang,et al.  Unsplit complex frequency-shifted PML implementation using auxiliary differential equations for seismic wave modeling , 2010 .

[10]  W. D. Roeck,et al.  Flow and geometrical effects on radiated noise from exhaust pipes computed by the Fourier pseudospectral time-domain method , 2015 .

[11]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[12]  George A. McMechan,et al.  Imaging conditions for prestack reverse-time migration , 2008 .

[13]  Bradley E. Treeby,et al.  Quantifying numerical errors in the simulation of transcranial ultrasound using pseudospectral methods , 2014, 2014 IEEE International Ultrasonics Symposium.

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

[15]  Qing Huo Liu,et al.  Large-scale simulations of electromagnetic and acoustic measurements using the pseudospectral time-domain (PSTD) algorithm , 1999, IEEE Trans. Geosci. Remote. Sens..

[16]  E. Baysal,et al.  Reverse time migration , 1983 .

[17]  Carlos Spa,et al.  Numerical Absorbing Boundary Conditions Based on a Damped Wave Equation for Pseudospectral Time-Domain Acoustic Simulations , 2014, TheScientificWorldJournal.

[18]  G. H. F. Gardner,et al.  FORMATION VELOCITY AND DENSITY—THE DIAGNOSTIC BASICS FOR STRATIGRAPHIC TRAPS , 1974 .

[19]  Qing Huo Liu,et al.  The perfectly matched layer for acoustic waves in absorptive media , 1997 .

[20]  A. Giannopoulos,et al.  A nonsplit complex frequency-shifted PML based on recursive integration for FDTD modeling of elastic waves , 2007 .

[21]  S. Gedney,et al.  An Auxiliary Differential Equation Formulation for the Complex-Frequency Shifted PML , 2010, IEEE Transactions on Antennas and Propagation.

[22]  Zhiming Li,et al.  Multistep reverse time migration , 2009 .

[23]  José Escolano,et al.  Impedance boundary conditions for pseudo-spectral time-domain methods in room acoustics , 2010 .

[24]  Larry Lines,et al.  Analysis of higher-order, finite-difference schemes in 3-D reverse-time migration , 1996 .

[25]  Paul L. Stoffa,et al.  An efficient hybrid pseudospectral/finite-difference scheme for solving the TTI pure P-wave equation , 2013 .

[26]  José Escolano,et al.  Semi-empirical boundary conditions for the linearized acoustic Euler equations using Pseudo-Spectral Time-Domain methods , 2011 .

[27]  Aria Abubakar,et al.  Application of the nearly perfectly matched layer in acoustic wave modeling , 2007 .

[28]  Qing Huo Liu,et al.  The PSTD Method with the 4th-Order Time Integration for 3D TAT Reconstruction of a Breast Model , 2014 .

[29]  Thomas Hagstrom,et al.  A New Construction of Perfectly Matched Layers for Hyperbolic Systems with Applications to the Linearized Euler Equations , 2003 .

[30]  B. Fornberg.  The pseudospectral method: Comparisons with finite differences for the elastic wave equation , 1987 .

[31]  Xu Wang,et al.  Computing prestack Kirchhoff time migration on general purpose GPU , 2011, Comput. Geosci..

[32]  Changsoo Shin,et al.  Acceleration of stable TTI P-wave reverse-time migration with GPUs , 2013, Comput. Geosci..

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

[34]  Qing Huo Liu,et al.  Direct $Z$-Transform Implementation of the CFS-PML Based on Memory-Minimized Method , 2015, IEEE Transactions on Microwave Theory and Techniques.

[35]  Qing Huo Liu,et al.  PERFECTLY MATCHED LAYERS FOR ELASTODYNAMICS: A NEW ABSORBING BOUNDARY CONDITION , 1996 .

[36]  J. Forssén,et al.  The extended Fourier pseudospectral time-domain method for atmospheric sound propagation. , 2010, The Journal of the Acoustical Society of America.

[37]  John C. Bancroft,et al.  Anisotropic reverse‐time migration for tilted TI media , 2007 .

[38]  Qing Huo Liu,et al.  The pseudospectral time-domain (PSTD) algorithm for acoustic waves in absorptive media. , 1998, IEEE transactions on ultrasonics, ferroelectrics, and frequency control.

[39]  Roland Martin,et al.  A High-Order Time and Space Formulation of the Unsplit Perfectly Matched Layer for the Seismic Wave Equation Using Auxiliary Differential Equations (ADE-PML) , 2010 .

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

[41]  Qing Huo Liu,et al.  A Multidomain PSTD Method for 3D Elastic Wave Equations , 2004 .

[42]  Manuel Kindelan,et al.  On the construction and efficiency of staggered numerical differentiators for the wave equation , 1990 .

[43]  Qing Huo Liu,et al.  The PSTD algorithm: A time-domain method requiring only two cells per wavelength , 1997 .

[44]  D. Komatitsch,et al.  An unsplit convolutional perfectly matched layer improved at grazing incidence for the seismic wave equation , 2007 .

[45]  Hong Liu,et al.  Accelerating finite difference wavefield-continuation depth migration by GPU , 2012, Applied Geophysics.

[46]  Haohuan Fu,et al.  Selecting the right hardware for reverse time migration , 2010 .

[47]  John Etgen,et al.  Computational methods for large-scale 3D acoustic finite-difference modeling: A tutorial , 2007 .