波动方程有限差分法频散衰减方法研究

波动方程式已被广泛用于正演和偏移中波场的数值模拟,但是,空间和时间上的微分近似值的精确才能确保有限差分法结果的精确。传统上,取得精确的数值需要通过使用相对更密集的有限差分网格或很长的有限差分算子;否则,数值误差(即数值频散)将会在很大程度上污染有效信号;然而,以上两种方法将导致大量增加的计算成本。这里提出一种简单和低计算成本的改进的波动方程式用来压制数值频散,当改进的波动方程式数值求解采用有限差分法时,数值频散相对于标准的波动方程式有明显的改善,且几乎没有增加额外的计算。二维数值试验表明,修改后的波动方程式极大地提高了图像质量,而且修改后的计算增量可忽略不计。