GPU-Accelerated Seismic Modeling

Efficiently modeling seismic data sets in complex 3D anisotropic media by solving the 3D elastic wave equation is an important challenge in computational geophysics. Using a stress-stiffness formulation on a regular grid, we tested a 3D finite-difference time-domain solver using a second-order temporal and eighth-order spatial accuracy stencil that leverages the massively parallel architecture of graphics processing units (GPUs) to accelerate the computation of key kernels.

This code was developed at the Centre for Petroleum Geosciences at the University of Western Australia. The software is integrated as part of the Madagascar open-source seismic processing suite.

You can check out my GTC 2014 talk about this code here: