有限元算法在声波方程数值模拟中的频散分析.pdf

第 卷 第 期 地 震 学 报 ， ３６ ５ Ｖｏｌ．３６ Ｎｏ．５ 　 　 　 　 　 　 年 月 （ ） ＡＣＴＡＳＥＩＳＭＯＬＯＧＩＣＡＳＩＮＩＣＡ ， ２０１４ ９ ９４４９５５ Ｓｅｔ．２０１４ 　 　 　 ｐ 印兴耀，周建科，吴国忱，梁锴 有限元算法在声波方程数值模拟中的频散分析 地震学报， （）： ．２０１４． ． ３６５ ９４４９５５． ： ／ ｄｏｉ１０．３９６９ ．ｉｓｓｎ．０２５３３７８２．２０１４．０５．０１７． ｊ ， ， ， ＹｉｎＸＹ ＺｈｏｕＪＫ ＷｕＧＣ Ｌｉａｎ Ｋ．２０１４．Ｄｉｓｅｒｓｉｏｎａｎａｌｓｉｓｆｏｒｔｈｅｆｉｎｉｔｅｅｌｅｍｅｎｔａｌｏｒｉｔｈｍｉｎａｃｏｕｓｔｉｃｗａｖｅｅｕａ ｇ ｐ ｙ ｇ ｑ ， （）： ： ／ ｔｉｏｎｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ．犃犮狋犪犛犲犻狊犿狅犾狅犻犮犪犛犻狀犻犮犪 ３６５ ９４４９５５．ｄｏｉ１０．３９６９ ．ｉｓｓｎ．０２５３３７８２．２０１４．０５．０１７． 犵 ｊ 有限元算法在声波方程数值 模拟中的频散分析 印兴耀 周建科 吴国忱 梁 锴 　 　 　 　 （中国山东青岛２６６５８０中国石油大学（华东）地球科学与技术学院） 摘要 针对有限元算法在地震波数值模拟中的数值频散问题，利用集中质量矩阵双线性插值 　 有限元算法，推导了二维声波方程的频散函数 在此基础上采用定量分析方法，对比分析了 ． 网格纵横长度比变化时的入射方向、空间采样间隔、地震波频率以及地层速度对数值频散的 影响 数值算例和模型正演结果表明：当采用集中质量矩阵双线性插值有限元算法时，为了 ． 有效地压制数值频散，在所使用震源子波的峰值频率对应的波长内，采样点数目应不少于 ２０ 个；减小网格长度的纵横比可以有效地抑制入射角（波传播方向与 轴的夹角）较小的地震波 狕 的数值频散；地震波频率越高，传播速度越慢，频散越严重，尤其是当相速度与其所对应的 频率比值小于 ２倍空间采样间隔时，不仅会出现严重的数值频散，还会出现假频现象． 关键词 双线性插值 有限元法 声波方程 集中质量矩阵 数值频散 　　 　 　 　 　 ： ／