1887

Abstract

Summary

The reverse time migration (RTM) cross-correlation imaging condition requires to access the forward-propagated and backward-propagated wavefields simultaneously. To reduce the wavefield storage and I/O cost of RTM, we propose a time-space dimension forward wavefield compression and reconstruction method that combines both Nyquist time sampling and spatial wavefield random sampling. Firstly, in the time dimension, the Nyquist sampling law is utilized to reduce the sampling rate, hence only the forward wavefields at Nyquist sampling points are stored. Subsequently, in the space dimension, taking into account the sparse nature of the seismic wavefield, the wavefield snap at each Nyquist sampling point is sampled using the Gaussian random measurement matrix, thereby achieving compressed sampling of the wavefield. Before the cross-correlation operation, the ONSL0 algorithm is utilized to effectively reconstruct the wavefield. Experiments using a 2D layered model and a 3D gradient model demonstrate that our method effectively reduces the memory requirements. In addition, the wavefield reconstruction error using the K-SVD dictionary is smaller than the DCT counterpart, albeit at the cost of a slightly increased computational time. The proposed method can also be used in the gradient calculation in full waveform inversion.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.202510646
2025-06-02
2026-03-13
Loading full text...

Full text loading...

References

  1. CandèsE J, RombergJ, TaoT. [2006]. Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information. IEEE Transactions on information theory, 52(2): 489–509.
    [Google Scholar]
  2. ClappR G. [2009]. Reverse time migration with random boundaries. Seg technical program expanded abstracts 2009. Society of Exploration Geophysicists, 2809–2813.
    [Google Scholar]
  3. GriewankA, WaltherA. [2000]. An implementation of checkpointing for the reverse or adjoint mode of computational differentiation. ACM Transactions on Mathematical Software, 26(1): 19–45.
    [Google Scholar]
  4. LeeD, ChungW. [2022]. Improving the memory efficiency of RTM using both Nyquist sampling and DCT based on GPU, Journal of Geophysics and Engineering, 19(4):706–723.
    [Google Scholar]
  5. Mulder, Plessix. [2004]. A comparison between one-way and two-way wave-equation migration. Geophysics, 69(6): 1491–1504.
    [Google Scholar]
  6. PeiC, ShiL G, LiS H, et al., [2023]. A wavefield reconstruction method using sparse representation and dictionary learning for RTM, Journal of Geophysics and Engineering, 20(5): 946–964.
    [Google Scholar]
  7. SunW, FuL Y. [2013]. Two effective approaches to reduce data storage in reverse time migration. Computers & Geosciences, 56: 69–75.
    [Google Scholar]
  8. SymesW W. [2007]. Reverse time migration with optimal checkpointing. Geophysics, 72(5): 213–221.
    [Google Scholar]
  9. TsaigY, DonohoD L. [2006]. Extensions of compressed sensing. Signal processing, 86(3): 549–571.
    [Google Scholar]
/content/papers/10.3997/2214-4609.202510646
Loading
/content/papers/10.3997/2214-4609.202510646
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error