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

从Maxwell方程出发，开展了三维大地电磁场所满足的边值问题研究，利用加权余量法导出了三维大地电磁有限元方程.介绍了三维矢量有限元六面体网格剖分方式、插值基函数选取，推导了三维大地电磁矢量有限元正演的单元刚度系数矩阵及离散格式.编制了三维矢量有限元大地电磁正演的Matlab程序.三维COMMEMI 3D-1模型的视电阻率曲线与国际通用的标准测试数据能很好地拟合，验证了作者编写的矢量有限元正演程序的正确性.通过对高、低阻异常体的阻抗张量形态分析，说明张量阻抗等值线图能用以大致判断异常体特性，丰富了大地电磁响应特征的表达方式.
Starting from the Maxwell equations, this article studied the boundary conditions of 3D MT. By using the weighted residual method, we derived the three-dimensional MT finite element equation. The three-dimensional vector finite element hexahedral meshing mode was introduced and the basis functions were selected. Then we derived the three-dimensional magnetotelluric vector finite element stiffness coefficient matrix and discrete format. A three-dimensional vector finite element magnetotelluric forward Matlab program was done. The apparent resistivity curve of the dimensional COMMEMI 3D-1 model matches the international standard test data, which proves the correctness of 3D magnetotelluric forward program. With the analysis of high and low resistivity anomalies, it shows that tensor impedance map can roughly determine the anomaly characteristics, which enriches the magnetotelluric response characteristics of expression.