一阶弹性波方程的变网格高阶有限差分数值模拟

摘要

使用可变网格的有限差分法进行地震模拟有许多独特的优点,主要表现为对地质模型的离散化更为合理,在低速带和复杂构造区域,可将局部网格划分得相对精细些,不仅提高了模拟精度,消除了因采样不足导致的频散现象,而且可以减少计算机内存需求,保持模型计算的灵活性.本文提出一种新的基于高阶交错网格技术的弹性波数值模拟方法,通过改麦网格的空间步长实现了局部网格加密技术,弥补了常规网格的缺陷和不足.试算结果表明,本文提出的算法稳定性较好,且能够提高模拟精度,减少计算时间,提高计算效率.