"Acceleration of seismic data processing with wave-packet decomposition using NVidia CUDA"
Nikitin V.V., Duchkov A.A., and Andersson F.

Seismic data are characterized by multidimensionality, large data sizes, and irregular structures. In this paper we consider an optimal decomposition of seismic data using the basis of Gaussian wave packets. We implemented and optimized a number of fast algorithms for forward and inverse transforms for three-dimensional seismic data decomposition. The algorithms implemented on GPU demonstrate 2-6x speedup compared to 20-core CPUs. The programs were tested on synthetic seismic data sets: data reconstruction by Gaussian wave-packet coefficients, data compression, denoising, and interpolation in the case of missing traces.

Keywords: Gaussian wave packets, Fourier transform, CUDA.

  • Nikitin V.V. – MAX IV Laboratory, Lund University, Fotongatan 2, 22592, Lund, Sweden; Ph.D., Postdoctoral Researcher, e-mail: viktor.nikitin@maxiv.lu.se
  • Duchkov A.A. – Trofimuk Institute of Petroleum Geology and Geophysics, Siberian Branch of Russian Academy of Sciences; prospekt Koptyuga 3, Novosibirsk, 630090, Russia; Ph.D., Head of Laboratory, e-mail: duchkovaa@ipgg.sbras.ru
  • Andersson F. – Center for Mathematical Sciences, Lund University, Solvegatan 18, 22100, Lund, Sweden; Ph.D., Associate Professor, e-mail: fredrik.andersson@math.lth.se