大地电磁三维矢量有限元正演模拟

摘要

从Maxwell方程出发,开展了三维大地电磁场所满足的边值问题研究,利用加权余量法导出了三维大地电磁有限元方程.介绍了三维矢量有限元六面体网格剖分方式、插值基函数选取,推导了三维大地电磁矢量有限元正演的单元刚度系数矩阵及离散格式.编制了三维矢量有限元大地电磁正演的Matlab程序.三维COMMEMI 3D-1模型的视电阻率曲线与国际通用的标准测试数据能很好地拟合,验证了作者编写的矢量有限元正演程序的正确性.通过对高、低阻异常体的阻抗张量形态分析,说明张量阻抗等值线图能用以大致判断异常体特性,丰富了大地电磁响应特征的表达方式.