通过引入高阶精度的近似解析离散算子,给出了一种求解声波方程的八阶NAD方法.数值误差分析和计算效率结果显示,与四阶LWC方法和八阶LWC方法相比,八阶NAD方法具有高数值精度、高计算效率和低存储量.应用NAD方法模拟地震波在复杂非均匀3层介质和Marmousi模型中的传播,数值结果表明该方法能有效压制数值频散,具有较强的地震波模拟适用能力.
We first present the eighth-order NAD method for solving the acoustic equation by introducing the nearly analytic discrete operator with high-order accuracy.Then the analysis of numerical errors and results of computational efficiency show the eighth-order NAD method has high numerical accuracy, high computational ef- ficiency and low computing storage, compared against the fourth-order Lax-Wendroff correction (LWC) scheme and the eighth-order LWC scheme.At last,we apply the eighth-order NAD method to model the seismic wave- field propagating in complex heterogeneous three-layer medium and Marmousi model.The numerical results illus- trate that the new method can effectively suppress numerical dispersion.And the eighth-order NAD method has strong application of seismic wave simulation.