基于正交贴体网格的黏弹介质地震波模拟

常规的有限差分地震波数值模拟采用笛卡尔坐标系的规则网格对计算区域进行网格剖分,在模拟起伏地表下的地震波场时,不仅不利于实现自由边界条件,而且由于阶梯状网格近似,在网格角点处容易产生虚假散射波,从而影响模拟精度。为此,将计算流体力学中的正交贴体网格生成技术引入起伏地表下黏弹介质的网格剖分中,采用优化的同位网格有限差分法计算曲线坐标系的一阶速度—应力方程,并通过选择性滤波法消除同位网格差分产生的格点振荡。正交贴体网格能够精确地描述起伏地表,网格的正交性也使得实施自由边界条件时无需做复杂的坐标转换和插值运算。数值算例表明:该方法得到的数值解与解析解完全吻合;与规则网格有限差分法模拟结果相比,在相同网格间距条件下,该方法能够有效消除阶梯状网格引起的虚假散射波,从而提高数值模拟精度。此外,起伏地表两层和三层黏弹介质模型模拟结果表明,该方法对复杂模型同样有较强的适用性。