基于谱元法的二维大地电磁数值模拟

论文摘要

谱元法（SEM）是基于有限元发展起来的一种数值模拟方法,具备有限元法处理复杂结构的几何灵活性和谱方法的高精度和指数收敛性,将其用于数值模拟计算比传统高阶有限元减少了内存需求和CPU时间。为了提高大地电磁数值模拟的精度和效率,基于Galerkin加权余量法,采用Gauss-Lobatto-Legendre（GLL）正交多项式作为插值基函数,用Intel的Pardiso求解器直接求解线性方程组得到电磁场分布,实现了大地电磁（MT）二维谱元法数值模拟。给出的层状模型和COMMEMI 2D-1模型的数值算例验证了谱元算法和所编写程序的可靠性和稳定性。

论文目录

• 0 引言
• 1 二维大地电磁正演理论
• 2 MT二维谱元法的实现
•   2.1 插值多项式及GLL数值积分
•   2.2 单元分析
•   2.3 视电阻率及相位的求解
• 3 数值计算及分析
• 4 结论

文章来源

类型: 期刊论文

作者: 方小姣,李长伟,刘剑

关键词: 谱元法,大地电磁,正交多项式,二维正演

来源: 矿产与地质 2019年05期

年度: 2019

分类: 工程科技Ⅰ辑,基础科学

专业: 地质学,地球物理学,矿业工程

单位: 桂林理工大学地球科学学院,广西隐伏金属矿产勘查重点实验室

基金: 国家自然科学基金（编号：41464002）,广西自然科学基金（编号：2013GXNSFAA019277）,桂林理工大学博士科研启动基金（编号：GUTQDJJ2011038）,广西隐伏金属矿产勘查重点实验室建设项目（编号：14-A-01,-05,12-071-20）共同资助

分类号: P631.325

页码: 861-867

总页数: 7

文件大小: 241K

下载量: 40