CuQ-RTM: A CUDA-based code package for stable and efficient Q-compensated reverse time migration

Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. F1-F15 ◽  
Author(s):  
Yufeng Wang ◽  
Hui Zhou ◽  
Xuebin Zhao ◽  
Qingchen Zhang ◽  
Poru Zhao ◽  
...  
Keyword(s):  
Large Scale ◽  
Gpu Computing ◽  
Graphics Processing Unit ◽  
Processing Unit ◽  
Program Optimization ◽  
Adaptive Stabilization ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Time Migration

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.

Fast seismic modeling and reverse time migration on a graphics processing unit cluster

2011 ◽  
Vol 24 (7) ◽  
pp. 739-750 ◽  
Author(s):  
Rached Abdelkhalek ◽  
Henri Calandra ◽  
Olivier Coulaud ◽  
Guillaume Latu ◽  
Jean Roman
Keyword(s):  
Graphics Processing Unit ◽  
Processing Unit ◽  
Seismic Modeling ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Migration ◽  
Graphics Processing

Enhanced prestack depth imaging of wide-azimuth data from the Gulf of Mexico: A case history

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB79-WB86 ◽  
Author(s):  
Xuening Ma ◽  
Bin Wang ◽  
Cristina Reta-Tang ◽  
Wilfred Whiteside ◽  
Zhiming Li
Keyword(s):  
Image Quality ◽  
Gulf Of Mexico ◽  
Large Scale ◽  
Velocity Model ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Delayed Imaging ◽  
Data Set ◽  
Depth Imaging ◽  
Time Migration

We present a case study of enhanced imaging of wide-azimuth data from the Gulf of Mexico utilizing recent technologies; and we discuss the resulting improvements in image quality, especially in subsalt areas, relative to previous results. The input seismic data sets are taken from many large-scale wide-azimuth surveys and conventional narrow-azimuth surveys located in the Mississippi Canyon and Atwater Valley areas. In the course of developing the enhanced wide azimuth processing flow, the following three key steps are found to have the most impact on improving subsalt imaging: (1) 3D true azimuth surface-related multiple elimination (SRME) to remove multiple energy, in particular, complex multiples beneath salt; (2) reverse-time migration (RTM) based delayed imaging time (DIT) scans to update the complex subsalt velocity model; and (3) tilted transverse isotropic (TTI) RTM to improve image quality. Our research focuses on the depth imaging aspects of the project, with particular emphasis on the application of the DIT scanning technique. The DIT-scan technique further improves the accuracy of the subsalt velocity model after conventional ray-based subsalt tomography has been performed. We also demonstrate the uplift obtained by acquiring a wide-azimuth data set relative to a standard narrow-azimuth data set, and how orthogonal wide-azimuth is able to enhance the subsalt illumination.

scholarly journals Reverse time migration of 3D vertical seismic profile data

Geophysics ◽  
2016 ◽  
Vol 81 (1) ◽  
pp. S31-S38 ◽  
Author(s):  
Ying Shi ◽  
Yanghua Wang
Keyword(s):  
Parallel Implementation ◽  
Processing Unit ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Time Step ◽  
Absorbing Boundary ◽  
Data Set ◽  
Time Migration ◽  
Wavefield Simulation ◽  
Vsp Data

Reverse time migration (RTM) has shown increasing advantages in handling seismic images of complex subsurface media, but it has not been used widely in 3D seismic data due to the large storage and computation requirements. Our prime objective was to develop an RTM strategy that was applicable to 3D vertical seismic profiling (VSP) data. The strategy consists of two aspects: storage saving and calculation acceleration. First, we determined the use of the random boundary condition (RBC) to save the storage in wavefield simulation. An absorbing boundary such as the perfect matching layer boundary is often used in RTM, but it has a high memory demand for storing the source wavefield. RBC is a nonabsorbing boundary and only stores the source wavefield at the two maximum time steps, then repropagates the source wavefield backwards at every time step, and hence, it significantly reduces the memory requirement. Second, we examined the use of the graphic processing unit (GPU) parallelization technique to accelerate the computation. RBC needs to simulate the source wavefield twice and doubles the computation. Thus, it is very necessary to realize the RTM algorithm by GPU, especially for a 3D VSP data set. GPU and central processing unit (CPU) collaborated parallel implementation can greatly reduce the computation time, where the CPU performs serial code, and the GPU performs parallel code. Because RBC does not need the same huge amount of storage as an absorbing boundary, RTM becomes practically applicable for 3D VSP imaging.

Leveraging the accelerated processing units for seismic imaging: A performance and power efficiency comparison against CPUs and GPUs

2017 ◽  
Vol 32 (6) ◽  
pp. 819-837 ◽  
Author(s):  
Issam Said ◽  
Pierre Fortin ◽  
Jean–Luc Lamotte ◽  
Henri Calandra
Keyword(s):  
Graphics Processing Units ◽  
Power Efficiency ◽  
Seismic Imaging ◽  
Graphics Processing Unit ◽  
Processing Unit ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Single Node ◽  
Time Migration ◽  
Graphics Processing

Oil and gas companies rely on high performance computing to process seismic imaging algorithms such as reverse time migration. Graphics processing units are used to accelerate reverse time migration, but these deployments suffer from limitations such as the lack of high graphics processing unit memory capacity, frequent CPU-GPU communications that may be bottlenecked by the PCI bus transfer rate, and high power consumptions. Recently, AMD has launched the Accelerated Processing Unit (APU): a processor that merges a CPU and a graphics processing unit on the same die featuring a unified CPU-GPU memory. In this paper, we explore how efficiently may the APU be applicable to reverse time migration. Using OpenCL (along with MPI and OpenMP), a CPU/APU/GPU comparative study is conducted on a single node for the 3D acoustic reverse time migration, and then extended on up to 16 nodes. We show the relevance of overlapping the I/O and MPI communications with the computations for the APU and graphics processing unit clusters, that performance results of APUs range between those of CPUs and those of graphics processing units, and that the APU power efficiency is greater than or equal to the graphics processing unit one.

Application of RTM 3D angle gathers to wide-azimuth data subsalt imaging

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB175-WB182 ◽  
Author(s):  
Yan Huang ◽  
Bing Bai ◽  
Haiyong Quan ◽  
Tony Huang ◽  
Sheng Xu ◽  
...  
Keyword(s):  
Model Updating ◽  
Velocity Model ◽  
Azimuth Angle ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Additional Information ◽  
Time Migration ◽  
Reflection Angle ◽  
Subsalt Imaging

The availability of wide-azimuth data and the use of reverse time migration (RTM) have dramatically increased the capabilities of imaging complex subsalt geology. With these improvements, the current obstacle for creating accurate subsalt images now lies in the velocity model. One of the challenges is to generate common image gathers that take full advantage of the additional information provided by wide-azimuth data and the additional accuracy provided by RTM for velocity model updating. A solution is to generate 3D angle domain common image gathers from RTM, which are indexed by subsurface reflection angle and subsurface azimuth angle. We apply these 3D angle gathers to subsalt tomography with the result that there were improvements in velocity updating with a wide-azimuth data set in the Gulf of Mexico.

Use of RTM full 3D subsurface angle gathers for subsalt velocity update and image optimization: Case study at Shenzi field

Geophysics ◽  
2011 ◽  
Vol 76 (5) ◽  
pp. WB27-WB39 ◽  
Author(s):  
Zheng-Zheng Zhou ◽  
Michael Howard ◽  
Cheryl Mifflin
Keyword(s):  
Model Building ◽  
Transverse Isotropy ◽  
Velocity Model ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Edge Diffraction ◽  
Time Migration ◽  
Image Optimization ◽  
Free Generation

Various reverse time migration (RTM) angle gather generation techniques have been developed to address poor subsalt data quality and multiarrival induced problems in gathers from Kirchhoff migration. But these techniques introduce new problems, such as inaccuracies in 2D subsurface angle gathers and edge diffraction artifacts in 3D subsurface angle gathers. The unique rich-azimuth data set acquired over the Shenzi field in the Gulf of Mexico enabled the generally artifact-free generation of 3D subsurface angle gathers. Using this data set, we carried out suprasalt tomography and salt model building steps and then produced 3D angle gathers to update the subsalt velocity. We used tilted transverse isotropy RTM with extended image condition to generate full 3D subsurface offset domain common image gathers, which were subsequently converted to 3D angle gathers. The angle gathers were substacked along the subsurface azimuth axis into azimuth sectors. Residual moveout analysis was carried out, and ray-based tomography was used to update velocities. The updated velocity model resulted in improved imaging of the subsalt section. We also applied residual moveout and selective stacking to 3D angle gathers from the final migration to produce an optimized stack image.

An efficient local operator-based Q-compensated reverse time migration algorithm with multistage optimization

Geophysics ◽  
2018 ◽  
Vol 83 (3) ◽  
pp. S249-S259 ◽  
Author(s):  
Tong Zhou ◽  
Wenyi Hu ◽  
Jieyuan Ning
Keyword(s):  
Seismic Attenuation ◽  
Depth Information ◽  
Frequency Noise ◽  
Processing Unit ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Pass Filter ◽  
Low Pass Filter ◽  
Time Migration ◽  
Imaging Results

Most existing [Formula: see text]-compensated reverse time migration ([Formula: see text]-RTM) algorithms are based on pseudospectral methods. Because of the global nature of pseudospectral operators, these methods are not ideal for efficient parallelization, implying that they may suffer from high computational cost and inefficient memory usage for large-scale industrial problems. In this work, we reported a novel [Formula: see text]-RTM algorithm — the multistage optimized [Formula: see text]-RTM method. This [Formula: see text]-RTM algorithm uses a finite-difference method to compensate the amplitude and the phase simultaneously by uniquely combining two techniques: (1) a negative [Formula: see text] method for amplitude compensation and (2) a multistage dispersion optimization technique for phase correction. To prevent high-frequency noise from growing exponentially and ruining the imaging results, we apply a finite impulse response low-pass filter using the Kaiser window. The theoretical analyses and numerical experiments demonstrate that this [Formula: see text]-RTM algorithm precisely recovers the decayed amplitude and corrects the distorted phase caused by seismic attenuation effects, and hence produces higher resolution subsurface images with the correct structural depth information. This new method performs best in the frequency range of 10–70 Hz. Compared with pseudospectral [Formula: see text]-RTM methods, this [Formula: see text]-RTM approach offers nearly identical imaging quality. Based on local numerical differential operators, this [Formula: see text]-RTM method is very suitable for parallel computing and graphic processing unit implementation, an important feature for large 3D seismic surveys.

Acoustic and elastic numerical wave simulations by recursive spatial derivative operators

Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. T167-T174 ◽  
Author(s):  
Dan Kosloff ◽  
Reynam C. Pestana ◽  
Hillel Tal-Ezer
Keyword(s):  
Synthetic Data ◽  
Analytic Solutions ◽  
Numerical Dispersion ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Data Set ◽  
Lamb’S Problem ◽  
Dynamic Elasticity ◽  
Time Migration ◽  
Recursive Operators

A new scheme for the calculation of spatial derivatives has been developed. The technique is based on recursive derivative operators that are generated by an [Formula: see text] fit in the spectral domain. The use of recursive operators enables us to extend acoustic and elastic wave simulations to shorter wavelengths. The method is applied to the numerical solution of the 2D acoustic wave equation and to the solution of the equations of 2D dynamic elasticity in an isotropic medium. An example of reverse-time migration of a synthetic data set shows that the numerical dispersion can be significantly reduced with respect to schemes that are based on finite differences. The method is tested for the solutions of the equations of dynamic elasticity by comparing numerical and analytic solutions to Lamb’s problem.

Reverse Time Migration Using the Pseudospectral Time-Domain Algorithm

2016 ◽  
Vol 24 (02) ◽  
pp. 1650005 ◽  
Author(s):  
Jiangang Xie ◽  
Zichao Guo ◽  
Hai Liu ◽  
Qing Huo Liu
Keyword(s):  
Time Domain ◽  
Acoustic Waves ◽  
Fdtd Method ◽  
Pseudospectral Method ◽  
Imaging Method ◽  
Processing Unit ◽  
Reverse Time ◽  
Reverse Time Migration ◽  
Eighth Order ◽  
Time Migration

We propose a pre-stack reverse time migration (RTM) seismic imaging method using the pseudospectral time-domain (PSTD) algorithm. Traditional pseudospectral method uses the fast Fourier transform (FFT) algorithm to calculate the spatial derivatives, but is limited by the wraparound effect due to the periodicity assumed in the FFT. The PSTD algorithm combines the pseudospectral method with a perfectly matched layer (PML) for acoustic waves. PML is a highly effective absorbing boundary condition that can eliminate the wraparound effect. It enables a wide application of the pseudospectral method to complex models. RTM based on the PSTD algorithm has advantages in the computational efficiency compared to traditional methods such as the second-order and high order finite difference time-domain (FDTD) methods. In this work, we implement the PSTD algorithm for acoustic wave equation based RTM. By applying the PSTD-RTM method to various seismic models and comparing it with RTM based on the eighth-order FDTD method, we find that PSTD-RTM method has better performance and saves more than 50% memory. The method is suitable for parallel computation, and has been accelerated by general purpose graphics processing unit.

Sign in / Sign up

Export Citation Format

Share Document