三角网格有限元法波动模拟的数值频散及稳定性研究

三角网格有限元法能够准确模拟复杂构造和复杂介质条件下的地震波场,数值频散和稳定性条件是地震波数值模拟中参数选择的主要依据.基于均匀的线性三角网格单元,根据结构刚度矩阵的组装原理以及平面波理论,推导了集中质量矩阵下两种网格结构的声波频散函数以及稳定性条件,并对数值频散特性以及稳定性进行了详细研究:三角网格单元中波动的数值频散除了受到空间采样间隔、单元网格纵横比和波传播方向等常规因素的影响外,还受到网格布局的影响,过锐或过钝的三角单元会对波动数值频散产生不良的影响,不同类型的单元网格、单元纵横比对应着不同的稳定性条件,正三角单元中的波动具有较好的数值频散特性,其数值各向异性(频散随波传播方向的变化)效应最弱,稳定性条件也较为宽松.最后通过数值模拟直观地验证了以上分析结果,为有限元正演三角网格的剖分和参数的设置提供一定的理论依据.