Using the rapid expansion method for accurate time-stepping in modeling and reverse-time migration

Reverse-time migration is based on seismic forward modeling algorithms, where spatial derivatives usually are calculated by finite differences or by the Fourier method. Time integration in general is done by finite-difference time stepping of low orders. If the spatial derivatives are calculated by high-order methods and time stepping is based on low-order methods, there is an imbalance that might require that the time-step size needs to be very small to avoid numerical dispersion. As a result, computing times increase. Using the rapid expansion method (REM) avoids numerical dispersion if the number of expansion terms is chosen properly. Comparisons with analytical solutions show that the REM is preferable, especially at larger propagation times. For reverse-time migration, the REM needs to be applied in a time-stepping manner. This is necessary because the original implementation based on very large time spans requires that the source term is separable in space and time. This is not appropriate for reverse-time migration where the sources have different time histories. In reverse-time migration, it might be desirable to use the Poynting vector information to estimate opening angles to improve the quality of the image. In the solution of the wave equation, this requires that one calculates not only the pressure wavefield but also its time derivative. The rapid expansion method can be extended easily to provide this time derivative with negligible extra cost.

[1]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[2]  Paul G. Richards,et al.  Quantitative Seismology: Theory and Methods , 1980 .

[3]  Edip Baysal,et al.  Forward modeling by a Fourier method , 1982 .

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

[5]  H. Tal-Ezer,et al.  Spectral methods in time for hyperbolic equations , 1986 .

[6]  M. A. Dablain,et al.  The application of high-order differencing to the scalar wave equation , 1986 .

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

[8]  Z. Koren,et al.  AN ACCURATE SCHEME FOR SEISMIC FORWARD MODELLING , 1987 .

[9]  Alfred Behle,et al.  NUMERICAL SOLUTION OF THE ACOUSTIC AND ELASTIC WAVE EQUATIONS BY A NEW RAPID EXPANSION METHOD1 , 1989 .

[10]  C. Edward.  High-order (space And Time) Finite-difference Modeling of the Elastic Wave Equation , 1990 .

[11]  Robert J. Geller,et al.  Optimally accurate second order time-domain finite difference scheme for computing synthetic seismograms in 2-D and 3-D media , 2000 .

[12]  K. Yoon,et al.  Challenges In Reverse-time Migration , 2004 .

[13]  Kurt J. Marfurt,et al.  Reverse-Time Migration using the Poynting Vector , 2006 .

[14]  Yu Zhang,et al.  Two-step Explicit Marching Method For Reverse Time Migration , 2008 .

[15]  P. Stoffa,et al.  Rapid Expansion Method (REM) For Time-stepping In Reverse Time Migration (RTM) , 2009 .

[16]  Sverre Brandsberg-Dahl,et al.  The Pseudo-analytical Method: Application of Pseudo-Laplacians to Acoustic And Acoustic Anisotropic Wave Propagation , 2009 .

[17]  Paul L. Stoffa,et al.  Time evolution of the wave equation using rapid expansion method , 2010 .