大地电磁三维矢量有限元正演模拟-2019年文档.docxVIP

大地电磁三维矢量有限元正演模拟 ：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. 大地电磁（MT是以电离层激发的天然交变电磁场为场源， 在地表观测相互正交的电场、 磁场分量来获取地电构造信息的一 种重要地球物理勘探方法 [1].MT 不需要庞大的发射源设备，只 需采用比较轻便的接收设备，野外工作方便、成本低，被广泛应 用于地壳和上地幔电性结构的研究， 在石油天然气勘探、 矿产资 源勘探、 工程与环境普查等领域， 发挥着举足轻重的作用 [2-9]. 可以预见，三维 MT勘探技术是地球物理中深层领域的研究热点 及今后MT的发展趋势，而三维 MT正演是理解MT勘探物理现象 并认识地质体电磁响应规律的有效手段，显然尤其重要 尽管矢量FEM拥有诸多优点，但在地球物理的电磁法正演领 域中，其应用并不多见，尚需要进一步完善 . 目前的研究主要包 括：Yoshimura等[10]开展了矢量FEM的MT响应数值模拟，并 将矢量FEM的计算结果与交错网格 FDM的计算结果进行了对比; Mitsuhata等[11]利用矢量FEM和节点FEM耦合的方法对三维 MT数值模拟；Nam [12]采用不规则六面体矢量 FEM直接计算电 场，研究了起伏地形下 MT的电阻率和相位的变化规律；刘长生 等[13] 将完全非结构化四面体单元引入到矢量有限元中， 实现了 三维大地电磁h-型自适应矢量有限元正演；王烨[14]开展了高 频率大地电磁法矢量有限元正演， 并采用改进的威尔金森方法求 解大型病态方程组， 提高了迭代速度； 顾观文等 [15] 开展了矢量 有限元法MT 三维地形数值模拟，研究了地形起伏下三维阻抗张 量的变化规律；杨军等 [16] 采用非结构四面体单元的三维矢量 FEM实现了海洋可控源电磁数值模拟；苏晓波等 [17]采用规则六 面体单元的三维矢量 FEM实现了大地电磁数值模拟，并对网格剖 分的重要性进行了研究 . 在前人基础上，作者推导了三维大地电磁矢量 FEM正演的离 散形式，应用矢量FEM算法计算了三维COMMEMI 3D-国际模型 [18]的MT视电阻率及模型张量阻抗，研究了高低阻异常体的电 磁响应特性，有效地指导了 MT的资料解释. 图5( a)，( b)分别为f=0.1 Hz时XY模式与YX模式下 矢量FEM正演视电阻率曲线.分析图5( a)，( b)可知，两幅 图中的矢量FEM曲线与COMMEMf提供的数据都能够很好地吻 合，说明无论是在低频还是高频部分， 应用矢量FEM开展三维大 地电磁正演，都具有较高的精度，同时也验证了矢量FEM算法及 程序的正确