三维电磁扩散场数值模拟及磁化效应的影响-物理学报.PDFVIP

物 理 学 报   Acta  Phys.  Sin.   Vol. 68, No. 3 (2019)    030201     三维电磁扩散场数值模拟及磁化效应的影响* 1) 1)† 1)2) 李志旋     岳明鑫     周官群 1) (中国科学技术大学地球和空间科学学院, 合肥　230026) 2) (安徽惠洲地质安全研究院股份有限公司, 合肥　231202) (2018 年8 月21 日收到; 2018 年12 月3 日收到修改稿) 采用矢量有限元法实现了三维电磁扩散场数值模拟, 并成功将其应用在大地电磁的正演研究中. 为灵活 精确地拟合起伏地形和地下不规则构造, 采用由不规则四面体单元组成的非结构化网格, 可根据模型设计的 需要调整网格的大小. 引入了基于二次场理论, 将解析的一次场从总场中扣除, 直接计算二次场, 使得误差仅 局限于相对较小的二次场, 以提高总场计算精度. 常规的节点有限元法不满足电性分界面上法向电场不连续 和无源区单元内电流密度无散, 违反麦克斯韦方程组. 为克服节点有限元法的弊端, 使用矢量有限元法求解基 于二次电场的偏微分方程. 另外, 在算法设计中, 考虑了磁导率参数的变化, 可以模拟磁导率不均匀的模型. 通过与COMMEMI模型已发表的结果对比, 证明了本文算法的正确性和精确性. 为突显非结构网格优势, 计算 了椭球异常体模型和任意地形模型的MT响应, 并详细讨论了地形和磁化效应对三维数值模拟结果的影响. 关键词：电磁扩散, 矢量有限元法, 数值模拟, 磁化效应, 大地电磁 PACS ：02.60.Cb, 02.70.Dh, 91.25.Qi, 91.25.wb 　DOI: 10.7498/aps.68 1   引　言 剖分, 对于复杂的模型仍难以剖分. 非结构化网格 的采用, 使得这一问题迎刃而解, 拓展了有限元模 电磁扩散场数值模拟方法中, 积分方程法、有 拟的应用范围[11−13]. 另外, 直接计算二次场以提高 限差分法和有限元法是三种经典并广泛使用的方 计算精度, 基于二次场的算法其一次场可解析地求 [1−3] 法. 积分方程法 受并矢格林函数的复杂性限制, 得, 不存在误差, 虽然二次场由数值法解得, 存在 [4−6] 仅适用于简单规则模型的模拟. 有限差分法 的 误差, 但二次场相对于一次场较小, 误差比重不大, 优点是易于实现, 其方程的离散形式决定了该方法 从而整体上提高了总场的精度. 只能模拟可四边形或六面体剖分的计算域, 对于复 可从两个角度来求解三维电磁场问题: 其一, 杂模型, 如起伏地形和不规则地下结构, 采取台阶 基于势 (场的矢势和标势)求解, 然后对势微分求 状近似误差大, 且剖分难实现. 相对而言, 采用非