高密度点法在隧道围岩分类勘探中的应用.docxVIP

高密度点法在隧道围岩分类勘探中的应用 近年来，高密度电法的应用越来越广泛。该方法的显著利益在于它拥有巨大的数据来源、高精度的勘探和广泛的应用。例如堤坝病害检测、岩溶勘探、地质分层以及在铁路勘探中都有很好运用。高密度电法作为一种电性勘探方法, 不可避免的受到地形条件的制约。在大多数勘探运用中, 由于都在较平地形上进行, 故可以不考虑地形影响。然而在地形起伏较大时, 如铁路、公路隧道勘探中, 它导致视电阻率曲线畸变严重。有时淹没地质体的信息, 有时歪曲有用信息, 有时呈现假信息, 使人们不能作出正确的推断解释, 降低了电阻率法的解释精度, 甚至导致错误的解释。因而, 研究地形影响及其改正方法, 是提高山区高密度电法勘探地质效果的一个关键。高密度地改方法有多种, 精度较高的是有限元法, 其次是如有限差分法、角域法等。有限元法计算量大, 一般在运用中不可能划分太多单元格, 因而也导致精度降低, 特别是如果要进行三维地改时, 其速度和效果更成问题。角域法相对精度低, 但突出优点是计算速度快, 特别适合野外工程勘探的要求。但在角域法地改计算中, 也因计算方法不同而导致精度差异, 本文将着重讨论这一点。 另外, 在铁路隧道物探围岩勘探中, 经常由于地质构造 (如断层) 有不同产状, 导致划分的围岩类别错误, 这时辅以高密度电法查明构造情况 (产状) , 可以大大提高围岩类别划分准确度。这也是我们要研究的内容。 1 复杂地形的修正曲线 我们借助数学公式解析单元地形的点源畸变电场, 再将单元地形的视电阻率畸变曲线组合叠加, 就能近似地取得复杂地形的改正曲线。此方法的基本思想为: (1) 复杂地形可以由一些简单的规则地形单元组成。只要能算出单元地形 (“角域”) 的畸变电场, 就能通过简便的组合叠加, 近似地算出复杂地形畸变。 (2) 点源场角域地形的电位是可以求解的。计算出它的电位畸变值, 再通过组合叠加获得复杂二度体地形的点源场畸变值。 1.1 rc点及规划法上的电位 研究地形对点源场的畸变影响问题的难度比线源场要大, 仅对最简单的二度体角域地形才能求得完全解。角域只有一个拐点, 两侧坡长“无限”, 如图1所示。角域顶角为?=π+β+α, 当?π时为无限长山脊, ?π时为无限长山谷。同时规定左坡的水平坐标或径向坐标取负值;右坡的取正值。角域下面均匀介质的电阻率为ρ0。 以拐点为原点建立圆柱坐标系 (r, θ, z) , 点电源在 (rc, 0, 0) , 测量点在 (rp, θ, z) 。 电位方程和边界条件如下: {?2V?r2+1r?V?r+1r2?2V?θ2+?2V?z2=0?V?θ|θ=0=0,?V?θ|θ=?=0V在rc点有奇性limr→∞V=0(1)????????????????2V?r2+1r?V?r+1r2?2V?θ2+?2V?z2=0?V?θ∣∣θ=0=0,?V?θ∣∣θ=?=0V在rc点有奇性limr→∞V=0(1) 用分离变量法解方程, 得电位公式: V=Ιρ0π?√|rcrp|{12Q-12(ξ)+∞∑m=1cosmπθ?Qmπ?-12(ξ)}(2)V=Iρ0π?|rcrp|√{12Q?12(ξ)+∑m=1∞cosmπθ?Qmπ??12(ξ)}(2) 其中ξ=|r2c+r2p+z22rcrp|ξ=∣∣r2c+r2p+z22rcrp∣∣。当rc与rp的符号相同时, θ=0;否则θ=?。Qv(ξ) 为第二类勒让德函数, V是它的阶。 在上述电位公式中, 由于Qv(1) 不存在, 不允许ξ=1。因此, 凡在rc, rp和z三者满足r2c+r2p+z22rcrp=±1r2c+r2p+z22rcrp=±1的电位上, 电位V有奇性。然而, 实际上V仅在电源点 (rc) 上, 即rp=rc, z=0情况时才有奇性, 其他情况均应有解, 如下: 当电源点与测量点分别在两侧山坡上, 而且两者到顶点的距离相等, 即rc=-rp, z情况时的电位公式为: V(rp,??0)=Ιρ0π?rc{π4+∞∑m=1(-1)m[φ(m+12)-φ(mπ?+12)]}(3)V(rp,??0)=Iρ0π?rc{π4+∑m=1∞(?1)m[φ(m+12)?φ(mπ?+12)]}(3) 式中?(x)=limn→∞[lnn-1x-1x+1-1x+2?-1x+n] 另外, 当电源点或测量点有一个在角域顶点上时, 电位公式为: V=Ιρ02?r′(4) 式中:r′——是顶点到测量点或电源点的距离。当z=0时, r′等于|rp|或|rc|。 利用 (2) 、 (3) 、 (4) 中的公式即可计算角域电位。 1.2 种装置类型的设计 通过上面的论述, 不难得到高密度的地形改正公式。设电极的间距为dx, 则三种装置类型, 即α, β, γ装置的地改计算公式为: