Acoustic prestack reverse time migration based on the symplectic method
Feng,Lanlan,Yang,Dinghui,Zhou,Yanjie
DOI: https://doi.org/10.6038/cjg20141227
2014-01-01
Abstract:Prestack reverse time migration (RTM) is currently one of the most accurate seismic migration imaging methods. An important step of its implementation is to solve the full-wave equation, so the numerical algorithm of solving the full wave equation quickly and effectively is crucial for reverse-time migration. The nearly analytic symplectic partitioned Runge-Kutta (NSPRK) method is a high-efficiency and high-precision numerical symplectic difference method developed in recent years, which can effectively suppress the numerical dispersion when a coarse grid is used. So it can increase computation efficiency and save computer memory demand. In this article, we develop a sixth-order NSPRK method based on the fourth-order NSPRK method, analyze the stability and numerical dispersion of the new sixth-order NSPRK method, compare its calculation efficiency and conduct wave field simulation. Then we apply it to the acoustic prestack reverse time migration, obtaining a sixth-order symplectic prestack reverse time migration method with low numerical dispersion, which can be applied to wavefield extrapolation with large steps and compute in long time. We get migration images of the Sigsbee2B model, and compare them to the fourth-order NSPRK method, the sixth-order conventional finite difference method, and the fourth-order Lax-Wendroff correction method (LWC). Numerical results show that the method based on the NSPRK prestack reverse time migration can get better imaging results than the fourth-order NSPRK method, the sixth-order conventional finite difference method, and the fourth-order LWC prestack reverse time migration method, especially in the case of coarse meshes.