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

被引:22
作者
刘长生 [1 ,2 ]
汤井田 [1 ]
任政勇 [1 ,3 ]
冯德山 [1 ]
机构
[1] 中南大学信息物理工程学院
[2] 长沙航空职业技术学院计算机系
[3] 瑞士联邦理工学院(EIH)地球物理系
基金
湖南省自然科学基金;
关键词
MT三维正演; 非结构化网格; 矢量有限元; 残差型后验误差; h-型自适应有限元;
D O I
暂无
中图分类号
P631.325 [];
学科分类号
0818 ; 081801 ; 081802 ;
摘要
基于能够模拟复杂模型的非结构化网格,提出基于矢量单元的三维自适应有限元大地电磁模拟算法。其过程是:利用残差型的后验误差算子初步估算粗网格上的单元误差,通过加密误差超过限定的单元,生成新的网格;对新的网格重复上一步过程,从而得到更加精确的数值结果;重复迭代过程直到计算结果的精度达到预定要求为止,从而生成最优化的有限元网格;基于COMMEMI 3D-1 MT模型的数值模拟,验证本文算法的正确性。研究结果表明:通过自适应的网格加密和迭代求解过程,本文算法可以产生迭代收敛的数值结果,计算结果具有较高的精度。
引用
收藏
页码:1855 / 1859
页数:5
相关论文
共 7 条
[1]  
基于矢量有限元的高频大地电磁法三维数值模拟.[D].王烨.中南大学.2008, 03
[2]   三维地电断面电阻率测深有限元数值模拟 [J].
阮百尧 ;
熊彬 ;
徐世浙 .
地球科学, 2001, (01) :73-77
[3]   Three-dimensional magnetotellurics modeling using edge-based finite-element unstructured meshes [J].
Liu Changsheng ;
Ren Zhengyong ;
Tang Jingtian ;
Yan Yan .
APPLIED GEOPHYSICS, 2008, 5 (03) :170-180
[4]  
Methods for modelling electromagnetic fields Results from COMMEMI—the international project on the comparison of modelling methods for electromagnetic induction.[J].M.S. Zhdanov;I.M. Varentsov;J.T. Weaver;N.G. Golubev;V.A. Krylov.Journal of Applied Geophysics.1997, 3
[5]  
电子计算机在电法勘探中的应用.[M].罗延钟;张桂青编著;.武汉地质学院出版社.1987,
[6]  
有限元素在电法测井中的应用.[M].李大潜编;.石油工业出版社.1980,
[7]  
地球物理中的有限单元法..徐世浙;.科学出版社.1994,