适于声波方程数值模拟的时间-空间域有限差分算子改进线性算法

摘要

有限差分法广泛应用于地震波场的数值延拓,确定合适的有限差分算子以减小数值频散是有限差分法的一个重要研究内容.近年来为了进一步抑制数值频散和增加时间步长,新的有限差分模板得到了应用,对于此,前人使用泰勒展开方法和最小二乘方法确定有限差分算子系数.本文在以前工作的基础上,使用改进的线性方法确定新模板的有限差分系数,并与传统模板线性方法进行对比;通过频散分析和正演模拟验证出新模板线性方法能够更好地保持频散关系,在相同的精度下效率提高了一倍,从而说明了改进的线性方法的有效性.