董浩1,3,,
1. 中国地质大学(北京)地球物理与信息技术学院, 北京 100083
2. 四川大学灾后重建与管理学院, 成都 610207
3. 地下信息探测技术与仪器教育部重点实验室, 北京 100083
4. 深地科学与工程教育部重点实验室, 四川大学, 成都 610065
基金项目: 国家高技术研究发展计划（"863"计划）（2014AA06A603），国家自然科学基金青年项目（41504062），深地科学与工程教育部重点实验室（四川大学）开放基金（DESEYU201906）及中央高校基本科研业务费专项资金联合资助
详细信息
作者简介: 郭泽秋, 男, 2017年毕业于中国地质大学(北京), 主要从事电性各向异性介质的大地电磁正反演模拟研究.E-mail:zeqiu_guo@hotmail.com
通讯作者: 董浩, 男, 博士, 讲师, 2013年毕业于中国地质大学(北京), 主要从事大地电磁三维正反演模拟及岩石圈电性结构研究.E-mail:donghao@cugb.edu.cn
中图分类号: P631收稿日期:2018-07-05
修回日期:2019-07-31
上线日期:2019-10-05
A Helmholtz decomposition based pre-condition method for magnetotelluric finite element numerical simulation
GUO ZeQiu1,2,4,,DONG Hao1,3,,
1. School of Geophysics and Information Technology, China University of Geosciences, Beijing 100083, China
2. Institute for Disaster Management and Reconstruction, Sichuan University, Chengdu 610207, China
3. Key Laboratory of Geo-detection of Ministry of Education, Beijing 100083, China
4. Key Laboratory of Deep Earth Science and Engineering(Sichuan University), Ministry of Education, Chengdu 610065, China
More Information
Corresponding author: DONG Hao,E-mail:donghao@cugb.edu.cn
MSC: P631--> Received Date: 05 July 2018
Revised Date: 31 July 2019
Available Online: 05 October 2019
摘要
摘要:本研究针对大地电磁测深法有限元数值模拟中，迭代法求解线性方程组效率较低的问题，利用亥姆霍兹分解原理，将电场矢量双旋度方程的预条件问题转化为基于矢量位的泊松问题和基于标量位的拉普拉斯问题，并在四面体非结构化棱边元离散的情况下，借助节点元辅助网格离散上述预条件问题，进一步利用代数多重网格方法（AMG）实施求解，最终实现预条件算法.利用经典的COMMEMI理论模型进行试算并与前人的积分方程解进行对比，验证了本文数值模拟程序与预条件方法的正确性和可靠性.此外，利用不同自由度规模的实验模型对这一预条件算法的效率进行了测试.结果表明，这一算法可以有效地提升大地电磁测深法棱边有限元数值模拟迭代法的收敛性，计算效率较通用的不完全LU分解预条件算法明显更高；在较大自由度网格（>1000万）数值模拟计算中，其算法效率及内存占用相对直接解法有较大优势，也使小型工作站上利用较大自由度的有限元网格进行大地电磁测深数值模拟计算成为可能.
关键词: 有限单元/
亥姆霍兹分解/
大地电磁/
预条件/
数值模拟
Abstract:To improve the efficiency of iterative methods solving linear equations in finite element simulation of magnetotelluric sounding,we present a pre-conditioning method exploiting the Helmholtz decomposition of vector fields. Helmholtz decomposition is used to transform the pre-conditioned problem of vector curl-curl equation of time-harmonic Maxwell equation into a Poisson problem based on the vector potential and a Laplacian problem with the scalar potential. With tetrahedral unstructured edge element,the preconditioned problem is discretized on auxiliary node element grid and then solved by algebraic multigrid method (AMG). Then,the COMMEMI synthetic model is used to calculate responses,which is compared with those of integral equation method,and verifies the reliability of the simulation program and the pre-condition method. Additionally,the efficiency of the pre-condition algorithm is tested by discretizing the model with a large range of degrees of freedoms (DoFs). The results show that the new algorithm can significantly improve the convergence of iterative methods that applied to magnetotelluric numerical simulation and presents significantly higher efficiency than the common pre-conditioner with incomplete LU decomposition. Under the circumstances of comparably large DoFs (over 10 million),it manifests apparent advantages over direct solver in efficiency and memory consuming,and also makes the magnetotelluric numerical simulation with a finite element mesh of relatively large size possible on a mobile workstation.
Key words:Finite element/
Helmholtz decomposition/
Magnetotellurics/
Pre-conditioning/
Numerical simulation
PDF全文下载地址:
http://www.geophy.cn/data/article/export-pdf?id=dqwlxb_15197
