Abstract: We gain the eight-order NAD-RK method based on the acoustic wave equation of two-dimension.This method uses the nearly analytic discretization method to conduct eight-order discretization on high order partial derivatives of the space,and employs the four-order Runge-Kutta method to conduct four-order discretization on temporal derivatives.The method is applied to wave-field numerical simulations from the two-layer acoustic and three-layer acoustic models in the 2-D case.The paper compares the method against the eighth-order LAX-Wendroff correction (LWC) and the eighth-order staggered-grid (SG) finite-difference methods.These results illustrate that the eight-order NAD-RK method has advantages such as weak numerical dispersion and high computational simulation effect.