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

Numerical simulation of two-dimensional magnetotelluric field based on spectral element method

作　　者：方小姣 李长伟[1,2] 刘剑 FANG Xiaojiao;LI Changwei;LIU Jian(School of Earth Sciences,Guilin University of Technology,Guilin 541006,Guangxi,China;Guangxi Key Laboratory of Hidden Metallic Ore Deposits Exploration,Guilin 541006,Guangxi,China)

机构地区：[1]桂林理工大学地球科学学院,广西桂林541006 [2]广西隐伏金属矿产勘查重点实验室,广西桂林541006

出　　处：《矿产与地质》2019年第5期861-867,共7页Mineral Resources and Geology

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

摘　　要：谱元法(SEM)是基于有限元发展起来的一种数值模拟方法,具备有限元法处理复杂结构的几何灵活性和谱方法的高精度和指数收敛性,将其用于数值模拟计算比传统高阶有限元减少了内存需求和CPU时间。为了提高大地电磁数值模拟的精度和效率,基于Galerkin加权余量法,采用Gauss-Lobatto-Legendre(GLL)正交多项式作为插值基函数,用Intel的Pardiso求解器直接求解线性方程组得到电磁场分布,实现了大地电磁(MT)二维谱元法数值模拟。给出的层状模型和COMMEMI 2D-1模型的数值算例验证了谱元算法和所编写程序的可靠性和稳定性。Spectral Element Method(SEM)is a numerical simulation method on the basis of finite element method.It has an advantage of modeling complex structure with high accuracy and exponential convergence of spectral method.SEM can reduce memory requirement and CPU time compared with traditional high-order finite element method.In order to improve the accuracy and efficiency of magnetotelluric numerical simulation,based on the Galerkin weighted residual method,the Gauss-Lobatto-Legendre(GLL)orthogonal polynomial is used as the interpolation basis function,and the Intel Pardiso solver is used to treat the linear equations directly to obtain the electromagnetic field distribution.It has realized for the two-dimensional spectrum element simulation of earth electromagnetic(MT).The reliability and stability of the spectral element algorithm and program are verified by numerical examples of layered model and COMMEMI 2D-1 model.

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

分 类 号：P631.3[天文地球—地质矿产勘探]