基于非结构化网格的二维大地电磁自适应有限元正演模拟研究

摘要

大地电磁测深法资料解释中正演和反演起着重要作用.正演是反演的基础,速度快高精度的正演为反演的实用化奠定基础.因此对正演算法的研究一直是MT数值模拟中研究的重点.有限元法是MT正演模拟中通常使用的方法.但有限法实现大地电磁正演模拟的过程中,利用规则化的网格模拟复杂模型,如带地形和断层构造,会产生较大的几何离散误差.模型的初始网格剖分程度也影响着正演模拟的精度,且同一套网格并不能满足所有频点的计算精度.因此,本文采用非结构化的四边形网格对研究区域进行网格剖分,从而能够模拟复杂的地电模型.并与自适应有限元相结合,从一个粗网格出发,通过网格单元的后验误差估计值来指导网格的局部加密以此获取较高精度的正演结果.通过对带地形的模型进行模拟计算,验证了该方法的有效性.