“…As one of the most popular methods of numerical solution for partial differential equations, the finite-difference method (FDM) has been widely utilized in seismic modeling [15][16][17][18] and migration [19][20][21]. To improve the accuracy and stability of a FDM in numerical modeling, many methods have been developed, including difference schemes of variable grid [22,23], irregular grid [24,25], standard staggered grid [26,27], rotated staggered grid [28], variable time step [29], and implicit methods [30].…”