三角谱元法及其在地震正演模拟中的应用-北京玛格泰克科技发展有限.PDFVIP

第 卷 第 期 地 球 物 理 学 报 ， ５７ ４ Ｖｏｌ．５７ Ｎｏ．４ 　 　 　 　 　 年 月 ， ２０１４ ４ ＣＨＩＮＥＳＥ ＪＯＵＲＮＡＬ ＯＦ ＧＥＯＰＨＹＳＩＣＳ Ａ ｒ．２０１４ 　 　 　 ｐ 李琳，刘韬，胡天跃 三角谱元法及其在地震正演模拟中的应用 地球物理学报， （）： ， ： ／ ．２０１４． ． ５７４ １２２４１２３４ｄｏｉ１０．６０３８ ｊｇ ， ， ＬｉＬ ＬｉｕＴ ＨｕＴＹ．２０１４．Ｓｅｃｔｒａｌｅｌｅｍｅｎｔｍｅｔｈｏｄｗｉｔｈｔｒｉａｎｕｌａｒｍｅｓｈａｎｄｉｔｓａ ｌｉｃａｔｉｏｎｉｎｓｅｉｓｍｉｃｍｏｄｅｌｉｎ．犆犺犻狀犲狊犲犑． ｐ ｇ ｐｐ ｇ （ ）， （）： ， ： ／ 犌犲狅犺狊． ｉｎＣｈｉｎｅｓｅ ５７４ １２２４１２３４ｄｏｉ１０．６０３８ 狆狔 ｊｇ 三角谱元法及其在地震正演模拟中的应用 １ ２ １ 李琳 ，刘韬 ，胡天跃 １北京大学地球与空间科学学院，北京 １００８７１ 　 ２中国石化石油勘探开发研究院，北京 １０００８３ 　 摘要 谱元法（ ）是基于有限元（ ）的一种算法，在地震正演模拟中应用广泛，但是大部分研究都是基于四 ＳＥＭ ＦＥＭ 　 边形网格下的谱元法 本文给出了 阶谱元法在三角网格中（ ）的基本原理，包括 形函数的构建，数 ． ２ ＴＳＥＭ Ｌａｒａｎｅ ｇ ｇ 值积分公式的选取 在此基础上，分析了 阶 方法的数值频散特性以及稳定性条件，并引入三角网格下 阶 ． ２ ＴＳＥＭ ３ 有限元方法进行分析对比，数值算例的结果证明 阶 相比于 阶 具有更高的计算精度，以及更宽松的 ２ ＴＳＥＭ ３ ＦＥＭ 稳定条件 最后，本文将 方法应用于中国西部地区的两个含溶洞介质的地质模型中，数值模拟结果表明 ． ＴＳＥＭ ＴＳＥＭ方法能够有效地模拟复杂结构的介质，有助于对地震波场传播特性的认识． 关键词 谱元法；三角网格；频散分析；稳定性分析