各向异性介质弹性波传播的三维不规则网格有限差分方法

提出一种新的三维空间不规则网格有限差分方法，模拟具有地形构造的非均匀各向异性介质中弹性波传播过程.该方法通过具有二阶时间精度和四阶空间精度的不规则交错网格差分算子来近似一阶弹性波动方程，与多重网格不同，无需在精细网格和粗糙网格间进行插值，所有网格点上的计算在同一次空间迭代中完成.针对具有复杂物性参数和复杂几何特征的地层结构，使用精细不规则网格处理粗糙界面、断层和空间界面等复杂几何构造,理论分析和数值算例表明，该方法不但节省了大量计算机内存和计算时间，而且具有令人满意的稳定性和精度.