A high accuracy seismic modeling and RTM method on GPU

Abstract

The article discusses a high accuracy seismic modeling and reverse time migration (RTM) method. High-order finitedifference (FD) equation is used for modeling and RTM with fourth-order in time domain, eighth-order in x, y direction of space domain and sixteenth-order in z direction. This high-order FD solution can effectively reduce numerical dispersion. PML absorbing boundary conditions is used to effectively solve artificial boundary reflection problem. Thereby, the accuracy of modeling and RTM are significantly improved. By the advantage of GPU, the computation efficiency is largely increased. Research of theoretical data and real data shows this method is very applicable to the modeling and imaging for the media with complex surface and complex underground structure.