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

来源期刊：矿产与地质2019年第5期

论文作者：方小姣 李长伟 刘剑

文章页码：861 - 867

关键词：谱元法;大地电磁;GLL正交多项式;二维正演;

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