“…Originated from early 1970s, the modern seismic migration method initially tried to back propagate the reflected seismic wave observed on the surface to its starting location through mathematical method, and gain the image which can indicate the location of the subsurface reflectivity (Claerbout, 1971(Claerbout, , 1976Ma et al, 1898;Robein, 2010;Stolt, 1986;Bleistein et al, 1999;Li et al, 2011). After 40 years development, seismic migration has developed from post-stack to pre-stack migration, from two-dimension to three-dimension, from time domain migration to depth domain migration, from isotropic medium to anisotropic medium, from scalar acoustic wave equation to vector elastic wave equation (Du, 2012;Chang, 1994), from conventional migration without iteration to true amplitude migration and iterative least squares migration (Zhang, 2013a;Yao, 2012aYao, , 2012bDong, 2012;Gray, 1997;Sun et al, 2002;Cui et al, 2004;Liu et al, 2006;Zhang, 2006;Li, 2008;Zhou et al, 2014), and it also tried to realize velocity and impedance inversion (Zhang, 2013b(Zhang, , 2013c. The calculation of seismic migration can be divided into two broad categories: 1) is based on Kirchhoff integral, 2) is based on solving differential equation.…”