2.5维电磁响应有限元模拟与波数取值的研究.pdfVIP

维普资讯 第30卷 第2期 物探化探计算技术 2008年3月 文章编号：1o01一l749(2008)O2 l35—1O 2．5维电磁响应的有限元模拟与波数取值研究 沈金松 ，孙文博 (1．中国石油大学 资源与信息学院地球物理系，北京 102249； 2．中国石油集团公司 物探重点实验室，北京 102249) 摘 要：利用有限元法实现了任意方向偶极子源在二维介质中频率域电磁响应的数值模拟，研 究 了波数取值对模拟结果的影响。通过对构造走 向的Fourier变换 ，将全三维电磁 问题 ，转化为 一 系列二维问题 ，并在波数域求解，极大地减小了计算工作量，导出了波数域耦合适用于二维电 性介质中任何方向电或磁偶极子响应计算的电磁场方程。针对每个给定的波数，上述耦合电磁 场方程用等参有限元方法在 一z平面内求解。采用Fourier逆变换，将波数域解积分，得到空间 域电磁场。针对电磁模拟计算中，源点的奇异性，采用具有一定面积的伪 函数表达源电流分 布，使数值解精度得以提高。另外，采用等参有限元，使地下复杂地质体得到准确表达。利用不 同波数值对均匀介质与层状介质的模拟结果与解析解的对比，验证了算法的正确性与精度。利 用层状介质模型的解析解与数值计算结果的对比，分析了波数的优化取值范围及取值点数对数 值模拟结果的影响，考察了算法对非均匀介质的适应性。 关键词：2．5维；波数范围；电磁响应模拟；等参有限元；伪 函数 中图分类号：P631．325 文献标识码：A 必须考虑地下，目标体的不规则形态及电性分布的 0 前言 非均匀性。同时，还应考虑这些问题对传统的模拟 方法 (如有限差分方法)存在巨大困难。另外，在 电磁测量能用于地层 电阻率研究，对油气勘 实际工作中，有时还不得不将源置于非均匀体附 探，固体矿勘查，地下水等地下资源探测，以及垃圾 近，此时，必须考虑源对非均匀体的影响 引。 等地下污染的评价，都具有重要的价值 J。地面、 作者在文中讨论的2．5维电磁问题，是指偶极 井中，及井地等不同探测方式的广泛应用所获取的 子源在二维地下介质中的响应。对于2．5维电磁 大量电磁测量数据，迫切需要有效的处理解释方 问题，国内、外均有学者开展研究 ，但由于源处理及 法与相应的软件。对于一维及二维电磁处理解释 地下复杂介质有效表达困难 ，公开发表的结果还十 方法已有较成熟的软件 J，2．5维与三维处理解 分有限。StoyerandGreenfield_1利用有限差分方 释方法逐渐为业界关注，也取得了一些有意义的成 法计算了频率域磁偶极子源在二维介质中的响应， 果_4J。对电磁数据的三维解释，尽管有最 自然和 对规则异常体得到较好的结果；Leppin 利用隐 合理的处理方法，但由于计算资源的限制，目前只 格式有限差分，计算了方形电流线圈的时间域电磁 限于简单模型的构造与粗剖分的实现 71。对于2．5 响应，对有限尺寸源得到了较好的结果；Leeand 维处理，主要障碍在于对有限尺寸源 的精确处 Morrison ]利用有限元方法计算了层状介质 中磁 理 】，以及地下介质复杂形态和电性非均匀性的 偶极子源的电磁响应 ；Sugeng等 采用等参单元 有效表达 J。为了模拟复杂介质，在数值计算中， 的有限元法计算了方形线圈的二次电磁场，避开了 收稿 日期：2007—05—09 改回 日期：2007—10—24