Abstract:Reverse time migration (RTM) in attenuating media should take absorption and dispersion effects into consideration. The latest proposed viscoacoustic wave equation with decoupled fractional Laplacians facilitates separate amplitude compensation and phase correction in [Formula: see text]-compensated RTM ([Formula: see text]-RTM). However, intensive computation and enormous storage requirements of [Formula: see text]-RTM prevent it from being extended into practical application, especially for large-scale 2D or 3D cases. The emerging graphics processing unit (GPU) computing technology, built around a scalable array of multithreaded streaming multiprocessors, presents an opportunity for greatly accelerating [Formula: see text]-RTM by appropriately exploiting GPUs architectural characteristics. We have developed the cu[Formula: see text]-RTM, a CUDA-based code package that implements [Formula: see text]-RTM based on a set of stable and efficient strategies, such as streamed CUDA fast Fourier transform, checkpointing-assisted time-reversal reconstruction, and adaptive stabilization. The cu[Formula: see text]-RTM code package can run in a multilevel parallelism fashion, either synchronously or asynchronously, to take advantages of all the CPUs and GPUs available, while maintaining impressively good stability and flexibility. We mainly outline the architecture of the cu[Formula: see text]-RTM code package and some program optimization schemes. The speedup ratio on a single GeForce GTX760 GPU card relative to a single core of Intel Core i5-4460 CPU can reach greater than 80 in a large-scale simulation. The strong scaling property of multi-GPU parallelism is demonstrated by performing [Formula: see text]-RTM on a Marmousi model with one to six GPU(s) involved. Finally, we further verified the feasibility and efficiency of the cu[Formula: see text]-RTM on a field data set.

RTM Using Effective Boundary Saving: A Staggered Grid GPU Implementation

Using the Effective Boundary Saving Strategy in GPU-based RTM Programming

Implementation strategy of 3D reverse time migration on GPU/CPU clusters

Improving the memory efficiency of RTM using both Nyquist sampling and DCT based on GPU

CuQ-RTM: A CUDA-based Code Package for Stable and Efficient Q-compensated Reverse Time Migration

GPU Accelerations on the 3D Elastic RTM Method

Simulation and reconstruction for 3D elastic wave using multi-GPU and CUDA-aware MPI

Least-squares Based Rectangular-Grid Cross and Rhombus Stencils for Acoustic Wave Propagation and Reverse Time Migration.

GPU-RRTMG_SW: Accelerating a Shortwave Radiative Transfer Scheme on GPU

Mesh-free Radial-Basis-function-generated Finite Differences and Their Application in Reverse Time Migration

Efficient boundary storage strategies for seismic reverse time migration

A Weighted Runge-Kutta Discontinuous Galerkin Method for Reverse Time Migration

Elastic reverse time migration based on nested triangular mesh in 2D isotropic media

GPU Computing for Longwave Radiation Physics: A RRTM_LW Scheme Case Study

A causal imaging condition for reverse time migration using the Discrete Hilbert transform and its efficient implementation on GPU

Analysis and application of boundary treatment for the computation of reverse-time migration

GPU-based Explicit Time Evolution Method

Wavefield Separation Algorithm of Helmholtz Theory-Based Discontinuous Galerkin Method Using Unstructured Meshes on GPU

Towards Efficient Reverse-time Migration Imaging Computation by Pipeline and Fine-grained Execution Parallelization.

3D prestack reverse time migration based on multiple-GPU collaboration

A domain-decomposition-based parallel approach for 3D geological modeling using radial basis functions interpolation on GPUs