摘要
根据贴体网格方法,将粘弹性波有限差分方法推广到含形起伏的模型中。差分计算采用MacCormack格式(Xie and Yao, 1988)。它是一种时间二阶、空间四阶度的格式,隐含对高频成分的耗散,易于消除由非物理因素产生的高频振荡。最后,利用能够精确处理起伏地形的谱元法对新的2D粘弹性有限差分方法进行了验证。将已经立的粘弹性模型,和具有高斯山谷和峡谷的介质模型,同时加载到谱元法程序specfem2D(Komatitsch aTromp, 1999)和有限差分程序中,模拟地震波场。无论是垂直分量还是水平分量波形,二者所的结果都具有较高的一致性,相关系数均达到95%以上。