High Order Finite Difference Forward Simulation of 2D and 3D Seismic Wave Field Based on GPU

This paper mainly introduces the seismic acoustic wave filed of high order finite difference forward and the problem of GPU computing. Through the numerical simulation technology, the forward problem of seismic wave propagation is realized. On the theoretical basis of the wave equation and with the utilization of Taylor series expansion, the wave equation's finite difference and its discrete form can be got. The 2D and 3D GPU forward procedures are written by C++ and CUDA, using the shared memory to prompt communication speed between different GPU threads, and improve the problem of 3D model sharing memory capacity limits on finite difference order. Different size models for 2D and 3D GPU and CPU program to calculate are established, running the program and getting consumption of time in order to compare the computational efficiency of the two procedures. The Test shows that, both in 2D and 3D model, the time consumption of GPU forward program is far less than that of CPU forward program, and with data amount increasing, the speedup is more and more significant. The test result proves that the efficiency of GPU program is better than that of CPU program.