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.

Optical Image Real-time MTF Compensation Method Based on CUDA

Extreme-Quality Computational Imaging via Degradation Framework

Non-blind optical degradation correction via frequency self-adaptive and finetune tactics

Image Restoration of Hybrid Time Delay and Integration Camera System with Residual Motion

Unified Optical Distortion Correction Method for Imaging Systems Using A Concise Geometrical Transformation Model

A flexible and accurate total variation and cascaded denoisers-based image reconstruction algorithm for hyperspectrally compressed ultrafast photography

Optimization on FFT of small size in real-time image restoration based on GPU

High-fidelity image reconstruction for compressed ultrafast photography via an augmented-Lagrangian and deep-learning hybrid algorithm

Monte Carlo-based fluorescence molecular tomography reconstruction method accelerated by a cluster of graphic processing units

GPU-accelerated Iterative Method for FD-OCT Image Reconstruction with an Image-Level Cross-Domain Regularizer.

Fast iterative method empowered by GPU acceleration for Fourier-domain optical coherence tomography image reconstruction

Realization of CUDA-based Real-Time Registration and Target Localization for High-Resolution Video Images

Image Reconstruction by Mumford–Shah Regularization for Low-Dose CT with Multi-Gpu Acceleration

Sophisticated Deep Learning with On-Chip Optical Diffractive Tensor Processing

Distortion correction algorithm for UAV remote sensing image based on CUDA

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

Acceleration of Early-Photon Fluorescence Molecular Tomography with Graphics Processing Units

Image quality enhancement using original lens via optical computing.

Optical Computing System for Fast Non-uniform Image Deblurring

High-performance real-world optical computing trained by in situ model-free optimization

High-performance real-world optical computing trained by in situ gradient-based model-free optimization