Two effective approaches to reduce data storage in reverse time migration

Prestack reverse time migration (RTM) requires extensive data storage since it computes wavefields in forward time and accesses wavefields in reverse order. We first review several successful schemes that have been proposed to reduce data storage, but require more computational redundancies. We propose two effective strategies to reduce data storage during RTM. The first strategy is based on the Nyquist sampling theorem, which involves no extra computational cost. The fact is that the time sampling intervals required by numerical algorithms or given by field records is generally several times smaller than that satisfied by the Nyquist sampling theorem. Therefore, we can correlate the source wavefields with the receiver wavefields at the Nyquist time step, which helps decrease storage of time history. The second strategy is based on a lossless compression algorithm, which is widely used in computer science and information theory. The compression approach reduces storage significantly at a little computational cost. Numerical examples show that the two proposed strategies are effective and efficient. Two approaches are proposed to reduce data storages in reverse time migration.The Nyquist scheme uses an imaging condition to reduce data storages.The Nyquist scheme introduces with no extra computational costs.The compression scheme employs compression algorithm to reduce data storages.The two approaches are proved to be effective via numerical examples.

[1]  D. Huffman A Method for the Construction of Minimum-Redundancy Codes , 1952 .

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

[3]  G. McMechan,et al.  Viscoelastic true-amplitude prestack reverse-time depth migration , 2008 .

[4]  Rui Yan,et al.  An angle-domain imaging condition for elastic reverse time migration and its application to angle gather extraction , 2012 .

[5]  Phil D. Anno,et al.  Seismic Data Compression And Regularization Via Wave Packets , 2010 .

[6]  Yu Zhang,et al.  A Stable TTI Reverse Time Migration And Its Implementation , 2009 .

[7]  Peter Moczo,et al.  Efficient Methods to Simulate Planar Free Surface in the 3D 4th-Order Staggered-Grid Finite-Difference Schemes , 2002 .

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

[9]  Alexey Romanenko,et al.  3D wave-packet decomposition implemented on GPUs [Электронный ресурс] , 2011 .

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

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

[12]  Paul Sava,et al.  Isotropic angle-domain elastic reverse-time migration , 2008 .

[13]  Giovanni Motta,et al.  Handbook of Data Compression , 2009 .

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

[15]  Necati Gulunay,et al.  Seismic trace interpolation in the Fourier transform domain , 2003 .

[16]  Hervé Chauris,et al.  Tips and tricks for Finite difference and i/o-less FWI , 2011 .

[17]  Yang Liu,et al.  Seislet transform and seislet frame , 2010 .

[18]  Tariq Alkhalifah,et al.  Acoustic approximations for processing in transversely isotropic media , 1998 .

[19]  Jon F. Claerbout,et al.  Imaging the Earth's Interior , 1985 .

[20]  Tariq Alkhalifah,et al.  An acoustic wave equation for anisotropic media , 2000 .

[21]  N. Whitmore Iterative Depth Migration By Backward Time Propagation , 1983 .

[22]  Robin P. Fletcher,et al.  Time-varying boundary conditions in simulation of seismic wave propagation , 2011 .

[23]  Bing Zhou,et al.  NUMERICAL SEISMOGRAM COMPUTATIONS FOR INHOMOGENEOUS MEDIA USING A SHORT, VARIABLE LENGTH CONVOLUTIONAL DIFFERENTIATOR1 , 1993 .

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

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

[26]  Abraham Lempel,et al.  A universal algorithm for sequential data compression , 1977, IEEE Trans. Inf. Theory.