一阶声波方程时间四阶精度差分格式的伪谱法求解

在地震波场数值模拟中,伪谱法不产生由空间网格离散引起的数值频散,而常规伪谱法用于求解时间二阶精度差分格式时,则会受到时间差分精度较低的影响而产生数值频散。本文基于一阶声波方程,提出将差分格式的时间差分精度增至四阶,并利用伪谱法求解,从而在避免由空间网格离散引起的数值频散的同时,降低由时间网格离散引起的数值频散。此外,与时间二阶精度差分格式伪谱法相比,时间四阶精度差分格式伪谱法的稳定性条件更为宽松,进而可通过增大时间网格步长提高计算效率。