摘要

在地震学中,真实的地球表面是一个强阻抗突变界面,可以通过自由表面边界条件对其进行描述。自由表面边界条件在数值模拟中的准确表达是目前地震波有限差分数值模拟的研究热点之一。为解决直接针对自由表面边界条件进行数值近似处理带来的稳定性和精度方面的问题,从平均介质理论出发,结合极限思想,推导了自由表面的另一种新的自适应数学表达方式。本质上,新的自由表面自适应表达方式是对原始弹性波方程中的密度和本构关系的修正,并且在计算过程中不会增加额外的CPU和内存需求。理论分析和数值实验表明,新的自由表面自适应表达方式是泊松比自适应的,对于不同泊松比值的有限差分数值模拟均能给出较为精确的模拟结果,从而可以应用于实际近地表不均匀风化层内的复杂波场的模拟和分析。