复杂地形条件下直流电阻率异常三维数值模拟研究-中国科技论文在线.PDF

第 39 卷 　第 4 期 地质与勘探 Vol. 39 　No. 4 2003 年 7 月 　　　　　　　　　　　　GEOLOGY AND PROSPECTING　　　　　　　　　　　　 July ,2003 技术·方法 复杂地形条件下直流电阻率异常三维数值模拟研究 1 2 1 熊 　彬 , 阮百尧 , 罗延钟 ( 1. 中国地质大学地球物理系 ,武汉 　430074 ;2. 桂林工学院 ,桂林 　541004) [摘 　要]复杂地形条件下的电阻率异常解释是十分困难的 : ①地形的起伏改变了地下电流场 ,使其 分布规律不同于水平地电条件下的电流场 ; ②在实际的观测过程中 ,测点至场源的距离与地表水平情况 的差异直接影响到了正常的解释工作。针对这一实际工作中遇到的问题 ,运用四面体单元的三维有限 元数值方法对复杂条件下的电阻率异常进行了理论研究 ,文中选用了几例典型地质模型。结果表明 ,采 用四面体单元的有限元法研究复杂条件的影响是可行的 ,有效的。 [ 关键词]复杂地形 　直流电阻率 　三维 　四面体单元 ( ) [ 中图分类号]P631. 3 　[文献标识码]A 　[文章编号]0495 - 5331 2003 04 - 0060 - 05 含 1681 个未知数的方程组。后来 , 徐世浙等[13 ,14 ] 0 　引 言 ( ) 1985 ,1991 开始用边界单元法研究这个问题 , 陈伯 在山区进行电阻率法勘探时[1 ] ,复杂的地形条 舫等[15 ] ( 1998) 又采用有限差分的方法计算了三维 件是个重要的干扰因素。它往往会导致视电阻率曲 地形影响的电磁感应。 线发生严重的畸变 ,产生形似高、低阻体的假异常 ; 众所周知 ,有限单元法最适合于物性复杂分布 有时甚至会掩没地质体的有用异常 ;有时又歪曲有 的地球物理问题 ,而有限差分和边界元法在这一领 用异常 ,造成地质体异常形态的变更和位置的移动。 域又各有其局限性。本文将在前期数值模拟工 而且 ,我国的金属矿又多位于山区[2 ] ,因此对复杂条 作[16 ,17 ] 的基础上 ,给出对具有复杂地形条件的地下 件下的电阻率异常规律进行研究具有重要的理论及 异常体进行有限元数值模拟的具体办法。 实际意义 ,是提高电阻率法地质效果的重要问题之 一。 1 　边值问题 以往的文献中有不少用数值模拟的办法来研究 三维构造中双点电源电场中电位的边值问题归 复杂条件下的电场变化规律 ,但这其中大都是对二 纳如下 : 维问题进行计算 ,或者是仅仅就纯地形的情况进行 (