基于非结构化网格的三维大地电磁自适应矢量有限元模拟.pdfVIP

第 41 卷第 5 期 中南大学学报( 自然科学版) Vol.41 No.5 2010 年 10 月 Journal of Central South University (Science and Technology) Oct. 2010 基于非结构化网格的三维大地电磁自适应矢量有限元模拟 1, 2 1 1, 3 1 刘长生 ，汤井田 ，任政勇 ，冯德山 (1. 中南大学 信息物理工程学院，湖南 长沙，410083； 2. 长沙航空职业技术学院 计算机系，湖南 长沙，410014； 3. 瑞士联邦理工学院(EIH) 地球物理系，瑞士 苏黎士，CH8092) 摘 要：基于能够模拟复杂模型的非结构化网格，提出基于矢量单元的三维自适应有限元大地电磁模拟算法。其 过程是：利用残差型的后验误差算子初步估算粗网格上的单元误差，通过加密误差超过限定的单元，生成新的网 格；对新的网格重复上一步过程，从而得到更加精确的数值结果；重复迭代过程直到计算结果的精度达到预定要 求为止，从而生成最优化的有限元网格；基于 COMMEMI 3D-1 MT 模型的数值模拟，验证本文算法的正确性。 研究结果表明：通过自适应的网格加密和迭代求解过程，本文算法可以产生迭代收敛的数值结果，计算结果具有 较高的精度。 关键词：MT 三维正演；非结构化网格；矢量有限元；残差型后验误差；h-型自适应有限元 中图分类号：P631.322 文献标志码：A 文章编号：1672−7207(2010)05−1855−05 Three-dimension magnetotellurics modeling by adaptive edge finite-element using unstructured meshes 1, 2 1 1, 3 1 LIU Chang-sheng , TANG Jing-tian , REN Zheng-yong , FENG De-shan (1. School of Info-physics and Geometrics Engineering, Central South University, Changsha 410083, China; 2. Department of Computer Science, Changsha Aeronautical Vocational and Technical College, Changsha 410014, China; 3. Institute of Geophysics, ETH-Swiss Federal Institute of Technology, Zurich CH8092, Switzerland) Abstract: Based on unstructured grid which can simulate complex geoelectric models, the vector edge-based finite-elements were proposed. An adaptive finite element three-dimensional magnetotelluric mode