摘要
交错网格有限差分已经被广泛地应用到地震波数值模拟与偏移成像中.传统的交错网格有限差分空间导数的高阶差分系数一般是通过Taylor级数展开法求取的,而用该方法确定的差分系数来求解空间导数时,一般只是在一个较小的频率范围内才能取得比较高的精度.本文针对交错网格有限差分的空间一阶导数的求解,从其频散关系出发,采用最小二乘方法推导一种在更大频率范围内都能满足该导数数值解精度要求的任意偶数阶精度差分系数,然后分析该差分系数求解空间一阶导数的数值频散,最后采用该最小二乘交错网格有限差分法对弹性波场进行数值模拟.同时与传统的Taylor级数展开交错网格有限差分的数值频散及波场模拟的精度和计算效率进行比较.数值频散和弹性波场模拟分析表明:在相同的离散参数和几乎不改变计算效率的情况下,最小二乘交错网格有限差分能有效地压制数值频散,并比传统的Taylor级数展开交错网格有限差分具有更高的数值模拟精度.