三角网格有限元法声波与弹性波模拟频散分析-地球物理学报.PDFVIP

第 卷 第 期 地 球 物 理 学 报 ， ５７ ８ Ｖｏｌ．５７ Ｎｏ．８ 　 　 　 　 　 年 月 ， ２０１４ ８ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ａｕ．２０１４ 　 　 　 ｇ 刘少林，李小凡，刘有山等 三角网格有限元法声波与弹性波模拟频散分析 地球物理学报， （）： ， ： ．２０１４． ． ５７８ ２６２０２６３０ｄｏｉ１０． ／ ６０３８ ｊｇ ， ， ， ＬｉｕＳＬ ＬｉＸＦ ＬｉｕＹＳ ｅｔａｌ．２０１４．Ｄｉｓｅｒｓｉｏｎａｎａｌｓｉｓｏｆｔｒｉａｎｌｅｂａｓｅｄｆｉｎｉｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄｆｏｒａｃｏｕｓｔｉｃａｎｄｅｌａｓｔｉｃｗａｖｅ ｐ ｙ ｇ （ ）， （）： ， ： ／ ｓｉｍｕｌａｔｉｏｎｓ．犆犺犻狀犲狊犲犑．犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ５７８ ２６２０２６３０ｄｏｉ１０．６０３８ 狆狔 ｊｇ 三角网格有限元法声波与弹性波模拟频散分析 １ １ １ ２ １ 刘少林 ，李小凡 ，刘有山 ，朱童 ，张美根 １中国科学院地质与地球物理研究所，中国科学院地球深部研究重点实验室，北京　１０００２９ ２中国石油化工股份有限公司石油物探技术研究院，南京 ２１００１４ 　 摘要 本文对声波与弹性波方程进行有限元法离散，构造有限元法频散关系的一般特征值问题，分析了时间离散 　 格式为中心差分的三角网格有限元法声波与弹性波模拟的频散特性 比较了三种质量矩阵即分布式质量矩阵、集 ． 中质量矩阵和混合质量矩阵对有限元法频散的影响；选取四种典型三角网格，分析了混合质量矩阵有限元 （ ）频散的方向各向异性；数值频散、方向各向异性随插值阶数的增加逐渐减弱，当空间为三阶插值时，频散 ＭＦＥＭ 主要表现为随采样率的变化而几乎无明显方向各向异性，其频散幅值也较小 控制其他影响因素不变的情况下， ． 研究了不同波速比介质中弹性波的数值频散 最后给出了三角网格 的数值耗散性 ． ＭＦＥＭ ． 关键词 有限元法；数值频散；三角网格；声波方程；弹性波方程 　 ： ／ 中图分类号 收稿日期 ， 收修定稿 ｄｏｉ１０．６０３８ Ｐ６３１ ２０１４０３０４２０１