Acceleration of computation speed for elastic wave simulation using a Graphic Processing Unit

Numerical simulation in exploration geophysics provides important insights into subsurface wave propagation phenomena. Although elastic wave simulations take longer to compute than acoustic simulations, an elastic simulator can construct more realistic wavefields including shear components. Therefore, it is suitable for exploration of the responses of elastic bodies. To overcome the long duration of the calculations, we use a Graphic Processing Unit (GPU) to accelerate the elastic wave simulation. Because a GPU has many processors and a wide memory bandwidth, we can use it in a parallelised computing architecture. The GPU board used in this study is an NVIDIA Tesla C1060, which has 240 processors and a 102 GB/s memory bandwidth. Despite the availability of a parallel computing architecture (CUDA), developed by NVIDIA, we must optimise the usage of the different types of memory on the GPU device, and the sequence of calculations, to obtain a significant speedup of the computation. In this study, we simulate two- (2D) and three-dimensional (3D) elastic wave propagation using the Finite-Difference Time-Domain (FDTD) method on GPUs. In the wave propagation simulation, we adopt the staggered-grid method, which is one of the conventional FD schemes, since this method can achieve sufficient accuracy for use in numerical modelling in geophysics. Our simulator optimises the usage of memory on the GPU device to reduce data access times, and uses faster memory as much as possible. This is a key factor in GPU computing. By using one GPU device and optimising its memory usage, we improved the computation time by more than 14 times in the 2D simulation, and over six times in the 3D simulation, compared with one CPU. Furthermore, by using three GPUs, we succeeded in accelerating the 3D simulation 10 times.

[1]  Henri Calandra,et al.  Fast seismic modeling and Reverse Time Migration on a GPU cluster , 2009, 2009 International Conference on High Performance Computing & Simulation.

[2]  Paulius Micikevicius,et al.  3D finite difference computation on GPUs using CUDA , 2009, GPGPU-2.

[3]  Klaus Schulten,et al.  Accelerating Molecular Modeling Applications with GPU Computing , 2009 .

[4]  William J. Dally,et al.  The GPU Computing Era , 2010, IEEE Micro.

[5]  Robert W. Graves,et al.  Simulating seismic wave propagation in 3D elastic media using staggered-grid finite differences , 1996, Bulletin of the Seismological Society of America.

[6]  John E. Stone.  Accelerating Molecular Modeling Applications with GPU Computing , 2009 .

[7]  Miguel Hermanns,et al.  Parallel Programming in Fortran 95 using OpenMP , 2002 .

[8]  Benjamin Edwards,et al.  Development of a Response Spectral Ground‐Motion Prediction Equation (GMPE) for Seismic‐Hazard Analysis from Empirical Fourier Spectral and Duration Models , 2015 .

[9]  Enhua Wu,et al.  Real-time 3D fluid simulation on GPU with complex obstacles , 2004, 12th Pacific Conference on Computer Graphics and Applications, 2004. PG 2004. Proceedings..

[10]  Nader Moussa.  Seismic imaging using GPGPU accelerated reverse time migration , 2009 .

[11]  C. Juhlin.  Imaging of fracture zones in the Finnsjön area, central Sweden, using the seismic reflection method , 1995 .

[12]  J. Virieux.  P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method , 1986 .

[13]  A. Levander.  Fourth-order finite-difference P-SV seismograms , 1988 .

[14]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[15]  John D. Owens,et al.  GPU Computing , 2008, Proceedings of the IEEE.

[16]  V. Červený,et al.  Seismic Ray Theory , 2001, Encyclopedia of Solid Earth Geophysics.

[17]  Liu You,et al.  Real-Time 3D Fluid Simulation on GPU with Complex Obstacles , 2006 .