Accelerating POCS interpolation of 3D irregular seismic data with Graphics Processing Units

Shu-Qin Wang, Xing Gao, Zhen-Xing Yao
Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
Computers & Geosciences, Volume 36, Issue 10, October 2010, Pages 1292-1300 (24 July 2010)


   title={Accelerating POCS interpolation of 3D irregular seismic data with Graphics Processing Units},

   author={Wang, S.Q. and Gao, X. and Yao, Z.X.},

   journal={Computers & Geosciences},





Source Source   



Seismic trace interpolation is necessary for high-resolution imaging when the acquired data are not adequate or when some traces are missing. Projection-onto-convex-sets (POCS) interpolation can gradually recover missing traces with an iterative algorithm, but its computational cost in a 3D CPU-based implementation is too high for practical applications. We present a computing scheme to speed up 3D POCS interpolation with Graphics Processing Units (GPUs). We accelerate the most time-consuming part of the 3D POCS algorithm (i.e. Fourier transforms) by taking advantage of a GPU-based Fourier transform library. Other parts are fine-tuned to maximize the utilization of GPU computing resources. We upload the whole input dataset to the global memory of the GPUs and reuse it until the final result is obtained. This can avoid low-bandwidth data transfer between CPU and GPUs. We minimize the number of intermediate 3D arrays to save GPU global memory by optimizing the algorithm implementation. This allows us to handle a much larger input dataset. When reducing the runtime of our GPU implementation, the coalescing of global memory access and the 3D CUFFT library provide us with the greatest performance improvements. Numerical results show that our scheme is 3X to 29X times faster than the optimized CPU-based implementation, depending on the size of 3D dataset. Our GPU computing scheme allows a significant reduction of computational cost and would facilitate 3D POCS interpolation for practical applications.
No votes yet.
Please wait...

* * *

* * *

HGPU group © 2010-2017 hgpu.org

All rights belong to the respective authors

Contact us: