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

摘要

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