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

Shouwei Liu,Huazhong Wang,Shengchang Chen,Xiangning Kong
DOI: https://doi.org/10.6038/cjg20131023
2013-01-01
Abstract:Prestack reverse time migration(RTM) is currently the most precise seismic imaging method. It needs appropriate implementation strategies and efficient computation platforms, because of the needs for large amount of calculation and large storage. Based on high order finite difference RTM, this paper focuses on solving the questions of the lack of video memory and artificial boundary, which must be solved for Graphic Processing Unit (GPU) RTM. Using the domain decomposition technique, we can efficiently realize three dimensional (3D) RTM of any production scale. It eliminates the limitation of GPU memory size. The control equation of the traditional perfectly matched layer (PML) boundary condition is different to that of the inner area. It is not suitable for GPU calculation. This paper realizes the Nearly PML (NPML) boundary condition on GPU. It makes high order finite difference calculation eliminate the needs of conditional branches. And the amount of auxiliary wavefields in boundary area is small. Therefore, wavefields propagation using NPML is efficient on GPU. 3D theoretical data and real field data RTM results prove the correctness of this paper's methods.
What problem does this paper attempt to address?