各向异性介质中波动方程有限元法模拟及其稳定性

被引:21
作者
周辉
徐世浙
刘斌
何樵登
机构
[1] 青岛海洋大学地质地球物理研究所!青岛
[2] 不详
[3] 长春地质学院地球物理系!长春
关键词
各向异性介质; 波动方程; 有限元法; 稳定性条件; 数值模拟;
D O I
暂无
中图分类号
P631.4 [地震勘探];
学科分类号
摘要
利用六面体单元和三线性插值函数,推导出在任意弹性各向异性介质中三维三分量波动方程所满足的有限元方程.给出用时域有限差分求解该有限元方程时,时间、空间采样率和弹性参数应满足的稳定性条件.提出避免每一时间步都求解一大型线性代数方程组的方法.在开始求波场前,先求质量矩阵的逆的因子形式,以后每时间步只需作刚度矩阵的逆的因子形式与向量的乘积运算,可大大提高有限元素法在微机上进行波动方程数值模拟的运算速度.结合两个数值模拟算例,说明该方法的可行性,并与伪谱法作了计算速度的比较,有限元法的计算速度高于伪谱法的计算速度.
引用
收藏
页码:833 / 841
页数:9
相关论文
共 7 条