GPU-based simulation of seismic wave propagation with hybrid PSM/FDM method

Hybrid PSM/FDM method, approximating the horizontal spatial derivative with pseudospectral method (PSM) while using finite difference method (FDM) to approximate the vertical spatial derivative, is an efficient technique in seismic wave simulation because it makes usage of both the high level of parallelism in FDM and the great accuracy in PSM. With natural parallelism, graphic processing unit (GPU) is very suitable for calculating such kind of problem, and compute unified device architecture (CUDA) drastically simplifies GPU programming. In order to improve the computing efficiency, in this paper, we implemented the 2D hybrid method for seismic wave simulation based on CUDA. The GPU-based hybrid method brings about significant speedup according to numerical tests performed in this paper and proves to be a reliable way for simulating large scale seismic wave propagation.