非均匀介质弹性波动方程不规则网格有限差分方法.pdf

第21卷第2期 　 计 算 力 学 学 报 　 V o l. 2 1 , N o . 2 2004 年 4 月 Ch inese Journa l of Com puta tiona l M echan ics A p ril 2004 文章编号:(2004) 非均匀介质弹性波动方程的不规则 网格有限差分方法 1, 2 1 2 孙卫涛 , 　杨慧珠 , 　舒继武 ( 1. 清华大学 工程力学系, 北京 100084; 2. 清华大学 计算机科学与技术系, 北京 100084) 摘　要: 从弹性波动方程出发, 提出了一种新的空间不规则网格有限差分方法, 并用于求解非均匀各向异性介质 中的弹性波正演问题。这种方法简单易行, 对于复杂几何结构, 例如低速层、套管井和非平面界面等, 在较细的不 规则网格上进行离散, 计算时间和占用内存更少。与多重网格差分方法相比, 该方法不需要粗、细网格之间的插 值, 所有网格差分计算在同一次空间迭代中完成。具有复杂几何交界面的模型计算, 包括地下透镜体、套管井眼 等, 在确定弹性常数和密度后, 用不规则网格的差分方法更易实现。该方法使用了H igdon 吸收边界条件解决人 工边界反射问题, 引入了新的稳定性条件和网格频散条件, 很好地消除了非物理散射波。理论模型的数值计算表 明, 该方法具有良好的稳定性和计算精度, 在模拟非均匀介质弹性波传播时, 比相同精度的规则网格有限差分 方法计算速度更快。该方法易于推广到非结构网格和三维问题中。 关键词: 有限差分; 不规则网格; 非均匀介质 中图分类号: P 631. 414　　　文献标识码:A 数的介质, 规则网格差分模型通常呈现阶梯形边 1　引　言 界, 非均匀介质弹性常数、密度等物理属性沿折线 理论研究和工程应用中经常遇到求解偏微分 边界进行几何平均, 差分计算也在这样的边界上进 方程的问题, 由于边界条件和介质参数过于复杂, 行, 无法准确模拟真实物理场的分布; 另外, 局部区 大部分方程无法得到解析解, 例如地震勘探中的弹 域的物理参数变化要求精细网格, 导致整个模型网 性波正演问题。有限差分方法是求解双曲型偏微分 格细分, 大大增加计算量。与之相比, 本文提出的不 方程的常用数值算法, 在规则网格上建立差分方程 规则网格差分方法具有格式简单, 运算速度快、占 来近似求解波动方程正演问题的研究数量庞大, 较 用内存较少的优越性, 在参数变化剧烈的局部区域 早见于A lterm an [ 1 ] , Kelly [ 2 ] 和V irieux [ 3, 4 ] 等文献。 和几何特征复杂的边界上, 使用较精细的不规则网 Sho rtley[ 5 ] 首先研究了L ap lace 方程数值解中的不 格, 其他区域使用空间步长较大的规则或不规则网 规则网格差分方法, 给出了矩形不规则网格差分格 [ 6 ] [ 7 ] 格, 很好地解决了区