基于非结构化网格的人工源频率域三维电磁场数值模拟.pdf

第３５卷 第３期 桂 林 理 工 大 学 学 报 Ｖｏｌ３５Ｎｏ３ ２０１５年８月 ＪｏｕｒｎａｌｏｆＧｕｉｌｉｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ Ａｕｇ　２０１５  文章编号：１６７４－９０５７（２０１５）０３－０４２３－０８　　　　　 ｄｏｉ：１０３９６９／ｊ．ｉｓｓｎ１６７４－９０５７２０１５０３００１ 基于非结构化网格的人工源频率域 三维电磁场数值模拟 １ ２ １ １ １ ３ 熊　彬 ，蔡红柱 ，罗天涯 ，李长伟 ，丁彦礼 ，赵建国 （１桂林理工大学 ａ地球科学学院；ｂ广西矿冶与环境科学实验中心，广西 桂林　５４１００４；２ＣｏｌｌｅｇｅｏｆＭｉｎｅｓ＆Ｅａｒｔｈ Ｓｃｉｅｎｃｅｓ，ＵｎｉｖｅｒｓｉｔｙｏｆＵｔａｈ，ＳａｌｔＬａｋｅＣｉｔｙ，ＵＳＡ，８４１１２；３中国石油大学 地球物理与信息工程学院，北京　１０２２４９） 摘　要：借助ＣＯＭＳＯＬ生成完全非结构化三维网格，提出了一种可用于求解非均匀介质中可控源电磁法 三维问题的有限元数值模拟方法。为了回避引入场源带来的奇异性，基于Ｃｏｕｌｏｍｂ规范，导出了关于耦合 二次磁矢量位及电标量位所满足的Ｍａｘｗｅｌｌ方程组；形成了Ｇａｌｅｒｋｉｎ大型线性代数方程；鉴于计算机存贮 空间及巨大的计算量，运用预条件复系数Ｋｒｙｌｏｖ子空间迭代法求解后续形成的大型稀疏对称非正定线性方 程组。主要成果体现在：笛卡尔坐标系网格生成及优良的局部网格加密；引入Ｃｏｕｌｏｍｂ规范和适当的边界 条件保证矢量位唯一；采用滑动最小二乘插值获取高精度二次矢量位旋度及标量位梯度。通过与同类研究 结果的比较，文中所提算法精度得到验证，且该方法具有良好的通用性。 关键词：可控源电磁法；非结构化网格；Ｃｏｕｌｏｍｂ规范；三维正演；有限单元法 中图分类号：Ｐ６３１３２５　　　　　　　　　　 　 　　文献标志码：Ａ 在减少数据解释歧义、降低勘探风险方面， （四面体）来填充二维或三维空间，它消除了结构 可控源电磁法 （ＣＳＥＭ）勘查一直以来都被认为是 化网格中节点的结构性限制，节点和单元的可控性 ［１－２］ 一种重要的地球物理手段 。不可否认的是， 好，能更好地处理边界，适用于模拟不规则计算区 ＣＳＥＭ勘探方法效果的优劣，将直接取决于人们对 域［１４］，在刻画复杂地质构造时有明显优势。非结 ［３］ 野外地质环境电磁响应的深入理解和认识 。随 构化网格在其生成过程中采用一定的准则进行优 着数值技术及计算机性能的快速发展，现在普通 化判断，能生成高质量的网格，很容易控制网格的 的台式工作站已经能容易地处理 ＣＳＥＭ三维正演 大小和节点的密度；此外，它采用随机的数据结构 问题。就当前而言，在电磁问题的正演模拟中， 也有利于进行网格自适应。一旦在边界上指定网格 应用最为广泛的数值方法就是有限单元法和有限 的分布，边界之间就可以自动生成网格，无需分块 ［４］ 差分法 ，尽管其他一些方法也伴随着得到了长 或用户的干预，而且不需要在子域之间传递信息。 足的提高和完善，譬如积分方程法［５－７］、有限体 与此同时，有限单元法能实施完全非结构化网格， ［８］ 积法 、对模型电导率施加某些额外约束的近似 使其单元边界能与地下非均质体的特有几何形态 ［９］ ［１０］ 求解方法 、Ｎｅｕｍａｎｎ级数法 、人工神经网络 达到一致，如偏离的钻孔、流体浸入带以及石油测 ［１１］ 法 ，以及基于低对比度假设的ＢｏｒｎＲｙｔｏｖ型近 井中经常遇到的倾斜地层，正是这种自然支持非 ［１２］ ［１３］ 似求解 ，等等 。 结构化网格的特性，使得在处理具有一定程度几何 完全非结构化三角形网格或四面体网格以 复杂性问题时，有限单元法成为了普遍性的首选 Ｄｅｌａｕｎａｙ法或