高密度电阻率法在边坡地形校正中的应用.docxVIP

高密度电阻率法在边坡地形校正中的应用 在起伏地形下进行电法测量时，获得的数据不可避免地受到地形起伏的影响，视场的信息量发生变化，存在有用地质信息的异常剖面。人们无法正确地作出正确的判断，也无法降低电法解释的准确性，也无法做出错误的解释。因此，地形变化是改善起伏地形电法调查地质影响的重要关键。地形改正一般采用“比值法”, 即用实测视电阻率除以“地形校正因子 (地形模型视电阻率/地形模型电阻率) ”, 从而得到改正后或平地形下的视电阻率。精度较高的是有限元法, 其次是有限差分法、边界单元法、角域法等。有限元法计算量大, 一般在运用中不可能划分太多单元格, 因而也导致精度降低, 特别是如果要进行三维地改时, 其速度和效果更成问题。角域法相对精度低, 但优点是计算速度快, 适合工程勘探的要求。 小砭沟位于延安市枣园镇, 交通较为便利, 面积约为65万m2, 共建房屋200栋。区内滑坡地质灾害较多, 基岩面起伏, 地质结构复杂, 不良地质体的成因、分布、性状具有相当的随机性, 对工程的影响也具有相当的不确定性, 并且隐伏于地下不易查明。因此, 在前期勘察基础上, 需采用地球物理勘探进一步查明滑坡体的饱和带和基岩层的分布情况, 为处治设计和施工提供工程地质依据。 区内地貌单元属黄土冲沟地貌, 四周属黄土梁峁与冲沟相间地貌, 地形影响严重。本次物探工作采用重庆地质仪器厂的DUK-2A高密度电阻率法测量系统, 数据采集采用温纳装置, 电极数60, 点距5m, 拟采用角域地改方法进行地形校正, 并结合钻孔资料分析地改效果。 1 边坡地形异常视电阻率计算 借助数学公式解析点源条件下单元角域地形的电源畸变电场, 采用Fortran语言编写有关程序, 并计算边坡地形上的纯地形异常视电阻率值。进行角域地形改正基本思想为:①通过数学方法求取简单角域地形的电位公式, 进而求取视电阻率畸变值;②通过组合叠加得复杂角域地形上的视电阻率畸变值, 得到纯地形影响;③利用比值法进行地形改正。 1.1 角域地形上的电位 所谓角域, 是指由两个半无限大平面所限的二度区域, 在区域中充满电阻率为ρ的均匀且各向同性导电介质。两半平面之张角ue788即为角域之顶角。两侧面与水平面之夹角为侧面之倾角, 分别用β和α表示左右两侧面之倾角, 规定斜坡在水平轴下面, 坡角取负值, 在上面取正值。则 显然, 当Φπ时, 角域即为坡长无限之二度山脊;当Φπ时, 角域即为坡长无限之二度山谷。角域是二度地形的基本单元, 在电阻率法的地形改正中很有用处。 现有一角域二度体, 角域的张角为Φ, 在角域以下有一电阻率为ρ的均匀介质, 地面以上为空气, 电阻率为无穷大。在地面的A点有一电流强度为I的点电源供电, A到坐标原点的距离为r0。为了求出电位的分布, 首先建立坐标系, 将x轴取为过供电点A而垂直于地形两斜面的交线, x轴与地形两侧面交线的交点为O, y轴取为过O点而垂直于地形的水平面, z轴为地形两侧面的交线。坐标轴正向如图1所示。 这样, 角域地形上M点电位解析解为 式中, ξ= (r2+r02+z2) / (2rr0) , Qγ(ξ) 为第二类勒让德函数。式 (2) 本身不完整, 存在以下问题: 第一, Qγ(1) 不存在, 故在此时解有奇异性, 不能用式 (2) 计算, 使用 当测量电极M在顶点时, 由电极互换原理得 角域地形点电源电场电位由式 (2) ~式 (5) 这几个电位公式及其附属公式计算。 为了突出地形影响, 用U/U0表示电位畸变值U (τ) (U0为背景场) , 在讨论z=0时, 电位畸变值是非常方便的。将式 (2) 除以U0=Iρ0/ (2π|r-r0|) 后, 得到 根据电极互换原理, 有U (τ) =U (1/τ) , 在计算U (τ) 时只考虑-1≤τ1的情况。 1.2 角域上点源场的视电阻率畸变值 实际地形是连绵起伏的。实践表明, 可像图2那样用以折代曲的办法拟合连续起伏的二度地形。将各个角域Φ1、Φ2、Φ3、Φ4的异常相乘叠加起来, 得到实际地形异常的近似值。无论地形形态如何, 地形异常主要发生在坡度改变点及其附近, 因此可以把角域看作二度地形的最基本单元。 当几个角域上供电电极、测量电极的实际水平位置相互对齐后, 复合地形上点源场的视电阻率畸变值可以由乘积叠加公式 决定。其中, ρs/ρ是连续有限地形上某一点的地形异常值, (ρs/ρ)1… (ρs/ρ)n是各个角域在相应点的地形异常值;β1…βn为各角域在相应点的坡度, 根据这种组合叠加的方法可以得到各种复杂地形的视电阻率异常, 得到纯地形影响。 1.3 地表电阻率剖面 当在地形起伏的斜坡上进行测量时, 由于观测点不在水平位置, 使地下电场的分布发生了变化, 与水平地面情况相比, 此时测得的电阻率包含了地