《三角网格有限元法声波与弹性波模拟频散分析.》.pdfVIP

第 57卷 第 8期 地 球 物 理 学 报 Vo1．57，No．8 2014年 8月 CHINESE JOURNAL OF GE0PHYSICS Aug．，2014 刘少林 ，李小凡，刘有山等．2014．三角网格有限元法声波与弹性波模拟频散分析．地球物理学报 ，57(8)：2620—2630，doi：10 6O38／cjg2O14O821． LiuSL，LiXF，LiuY S，eta1．2014．Dispersionanalysisoftriangle—basedfiniteelementmethodforacousticandelasticwave simulations．ChineseJ．Geophys．(inChinese)，57(8)：2620—2630，doi：10．6O38／cjg2O14O821 三角网格有限元法声波与弹性波模拟频散分析 刘少林 ，李小凡 ，刘有山 ，朱童 ，张美根 1中国科学 院地质与地球物理研究所 ，中国科学院地球深部研究重点实验室 ，北京 100029 2中国石油化工股份有限公司石油物探技术研究院，南京 210014 摘要 本文对声波与弹性波方程进行有限元法离散 ，构造有限元法频散关系的一般特征值 问题 ，分析 了时间离散 格式为 中心差分的三角网格有限元法声波与弹性波模拟的频散特性．比较了三种质量矩阵即分布式质量矩阵、集 中质量矩阵和混合质量矩阵对有限元法频散 的影 响；选取 四种典型三角网格 ，分析 了混合质量矩阵有限元 (MFEM)频散的方 向各向异性 ；数值频散 、方 向各 向异性随插值阶数 的增加逐渐减弱，当空间为三阶插值时，频散 主要表现为随采样率 的变化而几乎无明显方 向各 向异性，其频散幅值也较小．控制其他影响因素不变的情况下， 研究了不同波速 比介质中弹性波的数值频散．最后给出了三角网格MFEM 的数值耗散性． 关键词 有限元法 ；数值频散 ；三角网格 ；声波方程 ；弹性波方程 doi：10．6O38／cjg2O14O821 中图分类号 P631 收稿 日期 2014—03—04，2014-04—08收修定稿 Dispersionanalysisoftriangle-basedfiniteelementmethod foracousticandelasticwavesimulations LIU Shao—Lin ，LIXiao—Fan ，LIU You—Shan ，ZHU Tong ，ZHANG M ei—Gen 1KeyLaboratoryofEarth5DeepInterior，InstituteofGeologyandGeophysics， ChineseAcademyofSciences，Beijing100029，China 2SINOPECGeophysicalResearchInstitute，Nanjing210014，China Abstract Acousticand elasticwave equationsare discretizedby finite elementmethodwith triangles and centraldifference scheme in space and time， respectively． The generalized eigenvalueproblemsare constructed to analyze thenumericaldispersion ofthe finite element method．Threetypesofmassfiniteelementmethod，whichareconsistentfiniteelementmethod (CFEM)，lumpedfiniteel