二维地震波场的五点八阶超紧致有限差分数值模拟

摘要

首先将迎风机制引入五点八阶超紧致有限差分(CCD8)格式,得到五点七阶迎风超紧致(UCCD7)格式,并对两种格式进行数值频散分析和精度分析;其次建立了二阶声波方程的位移场时间四阶离散格式,将五点CCD8格式和五点UCCD7格式分别应用于位移场空间导数的求取,并分析这两种格式的稳定性条件;最后基于PML边界条件,将上述两种格式分别应用于声波方程的均匀介质、水平层状介质及Marmousi模型的数值模拟和波场特征分析及对比.研究结果表明:相较于普通紧致差分,五点CCD8格式具有小截断误差、高模拟精度、低数值频散、高稳定性、所需网格点数少的优点;引入迎风机制后,声波方程的五点UCCD7格式稳定性得到进一步提高;模型试算的结果验证了五点CCD8格式适用于复杂介质的数值模拟,模拟精度和计算效率都高.