有限单元法在电阻率法地形改正中的应用 - 矿产勘查.pdf

第２卷　第３期 矿 产 勘 查 Ｖｏｌ．２　Ｎｏ．３ ２０１１年 ５月 　　　　　　　　　　　　　　　ＭＩＮＥＲＡＬＥＸＰＬＯＲＡＴＩＯＮ　　　　　　　　　　　　　　　 Ｍａｙ，２０１１ 有限单元法在电阻率法地形改正中的应用 １ ２ １ 梁芳敏 ，欧东新 ，魏继祖 （１．河南省有色金属地质矿产局第六地质大队物化探院，洛阳　４７１００２； ２．桂林理工大学地球科学学院，桂林　５４１００４） 摘　要　文章应用有限单元法进行点源二维电场数值模拟计算。在正演过程中采用三角单元剖分网格，有利于对 复杂地形的模拟。电导率（ ）和电位在相邻节点间线性变化，使得电导率和电位与实际地质构造更接近、 σ 与实际地电场变化规律更易吻合；在点源二维电场正演基础上，将一定地形条件下的二维地电模型的视 电阻率响应，分解为畸变场和正常场两部分，采用比值法消除由地形起伏引起的视电阻率畸变，从而达到 地形校正的目的。并通过设计几个模型及实际算例，验证该方法对地形校正的效果。 关键词　有限单元法　数值模拟　视电阻率　地形校正 中图分类号：ＴＰ３０１；Ｐ６３１．３　文献标识码：Ａ　　　文章编号：１６７４－７８０１（２０１１）０３－０２９３－０５ 将三维问题变成二维问题，再应用有限单元法对积 ０　引言 分区域进行三角单元剖分，经过单元分析、全体单元 在电法勘探的实际工作中，常常受到地形因 相加合成： 素的干扰影响，如在山区及丘陵地带的起伏地形 １ Ｔ Ｔ Ｆ（Ｕ）＝ Ｆ（Ｕ）＝ ＵＫＵ＋ＵＰ （１） 等，使视电阻率曲线发生畸变。这些干扰使视电  ｅ ２ 其中Ｋ ＝Ｋ ＋Ｋ ＋Ｋ ，Ｋ＝ Ｋ，是总体系数矩阵。 阻率曲线变化复杂，或者使异常曲线特征点变得 ｅ １ｅ ２ｅ ３ｅ  ｅ 不明显。由于地形影响因素的存在，直接影响和 当双电流源供电时 降低了电法勘探的地质效果，忽略这些影响将使 １ １ Ｔ Ｐ＝ ０ ０…０＋ Ｉ　 － Ｉ　０…０　０ ， 推断解释工作难以进行，甚至导致错误的地质结 ( ２ ２ ) 论。因此，弄清楚地形对电剖面电阻率异常曲线 对（１．１）式求变分并令其为零： Ｔ Ｔ Ｔ Ｆ（Ｕ）＝ ＵＫＵ－ ＵＰ＝ Ｕ（ＫＵ－Ｐ）＝０ 的影响，并采取有效的措施即地形改正，对电法实 δ δ δ δ 　　由于 Ｕ ０，所以， 测数据进行预处理是电法资料解释前非常重要的 δ ≠ 工作。计算点源二维电场的地形影响的数值方法 ＫＵ－Ｐ＝０　或ＫＵ＝Ｐ