基于有限体积法直流电阻率25维数值模拟.pdfVIP

第３４卷 第３期 桂 林 理 工 大 学 学 报 Ｖｏｌ３４Ｎｏ３ ２０１４年８月 ＪｏｕｒｎａｌｏｆＧｕｉｌｉｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ Ａｕｇ　２０１４  文章编号：１６７４－９０５７（２０１４）０３－０４１１－０５　　　　　 ｄｏｉ：１０３９６９／ｊ．ｉｓｓｎ１６７４－９０５７２０１４０３００１ 基于有限体积法直流电阻率２５维数值模拟 １，２ １ １ ２ １ 熊　彬 ，刘云龙 ，李长伟 ，ＣＡＩＨｏｎｇｚｈｕ，罗天涯 （１桂林理工大学 地球科学学院，广西 桂林　５４１００４； ２ＤｅｐａｒｔｍｅｎｔｏｆＧｅｏｌｏｇｙａｎｄＧｅｏｐｈｙｓｉｃｓ，ＵｎｉｖｅｒｓｉｔｙｏｆＵｔａｈ，Ｕｔａｈ８４１１２，ＵＳＡ） 摘　要：为进一步解决传统直流电法正演计算中所存在的问题，给出了基于有限体积法的直流电阻率法 ２５维数值模拟算法。针对数值模拟中至为关键的场源奇异性及边界效应，介绍了一种右端项校正技术。 通过构建若干典型的地电断面模型，计算地表电位及视电阻率，并对其异常响应进行分析。算例结果表 明，通过校正方程组右端项后产生的源矢量不再是包含近源处突变电位的矢量，它能够补偿正演算符因边 界取得不够远而丢失的信息，进而减少边界及源奇异性对正演模拟的误差，验证了本文算法的可行性。 关键词：直流电阻率；有限体积法；２５维；数值模拟 中图分类号：Ｐ６３１３　　　 　　　　　 　　　　　　文献标志码：Ａ ［１６］ 直流电阻率法已被广泛的应用于寻找金属、 的正反演问题；高增等 （２０１０）研究了基于有 非金属矿床、工程地质勘查、环境监察等领域。 限体积的感应测井正演模拟响应。 开展快速精确有效的电阻率响应数值模拟对资料 本文在前人工作的基础上，对直流电阻率中 的采集与地质解释至关重要［１－３］。现阶段，直流 的２５维问题进行了研究。为了提高近源处模拟 电阻率数值模拟方法主要集中在有限单元法［４－７］、 精度及消除边界效应，采用了一种右端项校正的 ［８］ ［９－１０］ 有限差分法 、边界单元法 。有限体积法是 方法，对典型地电构造的表层电位及视电阻率异 在有限单元与有限差分的基础上发展起来的一种 常特征进行了讨论。 数值方法，它既吸收了有限差分的便捷性，同时 １　边值问题 又融合了有限元网格剖分的灵活性，并且在整个 计算区域遵循积分守恒原理。边界处理也较为简 假设地下导电性分布沿地质构造走向方向 （ｙ 单，计算效率高［１１］，因此在地球物理正演计算方 方向）无变化，与之相应的波数域中直流电场２５ ［１２］ ［１７］ 法研究中逐渐被关注。Ｄｅｙ等 在１９７９年首先将 维边值问题 为 有限体积法应用到电阻率数值模拟并取得良好的 ２ 　 ·［（ｘ，ｚ） Ｕ（ｘ，ｋ，ｚ）］－ｋ （ｘ，ｚ）Ｕ（ｘ，ｋ，ｚ）  σ  σ 效果；罗延钟等［１３］ （１９８７）采用三角网格剖分， 　　 ＝－Ｉ（Ａ）， ； （１ａ） δ ∈Ω 利用有限体积法计算了地形对视电阻率的影响