¹State Key Laboratory of Coal Resources and Mine Safety, China University of Mining and Technology, Beijing, China
²School of Earth and Space Sciences, Peking University, Beijing, China

---

Abstract—The staggered-grid finite difference (FD) method demands significantly computational capability and is inefficient for seismic wave modelling in 2-D viscoelastic media on a single PC. To improve computation speedup, a graphic processing units (GPUs) accelerated method was proposed, for modern GPUs have now become ubiquitous in desktop computers and offer an excellent cost-toperformance- ratio parallelism. The geophysical model is decomposed into subdomains for PML absorbing conditions. The vertex and fragment processing are fully used to solve FD schemes in parallel and the latest updated frames are swapped in Framebuffer Object (FBO) attachments as inputs for the next simulation step. The seismic simulation program running on modern PCs provides significant speedup over a CPU implementation, which makes it possible to simulate realtime complex seismic propagation in high resolution of 2048*2048 gridsizes on low-cost PCs.

Index Terms—seismic, wave propagation, finite difference, viscoelastic media, parallel algorithms, GPU

[PDF]

Cite: Zhangang Wang, Suping Peng, Tao Liu, "GPU accelerated 2-D staggered-grid finite difference seismic modelling," Journal of Software vol. 6, no. 8, pp. 1554-1561, 2011.