[发明专利]多同性源电极阵列电阻率的地下工程超前探测方法及系统

说明书

技术领域

本发明涉及隧道或者坑道等地下工程钻爆法或TBM施工过程中不良地质的探测方法，尤其涉及一种基于多同性源电极阵列电阻率的地下工程超前探测方法及系统。

背景技术

随着水利水电工程和交通工程的建设重点逐步向地形地质条件极端复杂的西部山区和岩溶地区转移，所遇到的工程地质条件异常复杂，施工中面临着突水、突泥等诸多高风险地质灾害。突涌水灾害给施工安全带来了重大灾难和生命财产损失，而预先探明隧道掌子面前方的地质灾害赋存源是规避和防治重大灾害事故的有效途径，因此开展隧道超前预报工作研究显得尤为迫切和重要。

作为一种重要的工程地球物理探测技术，直流电阻率法以其适应性强、方便快捷、成本低等诸多优点被广泛应用于岩土工程勘察、滑坡地质灾害监测、环境岩土工程评价等领域中。与地面直流电阻率法相比，由于隧道掌子面空间狭小，导致地面的很多探测方法无法应用于隧道突涌水灾害源超前预报。目前直流电阻率法在隧道（巷道）超前预报中大都采用三极法的探测模式，即将供电点A置于掌子面，供电点B置于隧道后方“无穷远”处，测量电极M、N沿着隧道轴向移动测量。

这一方法主要面临以下问题：①由于测量电极MN在掌子面后方移动测量，易受到掌子面后方异常体的干扰；②由于隧道的特殊环境，隧道三极法的探测模式测线垂直于异常体，导致对掌子面近处的异常体响应敏感，对远处的异常体响应相对不敏感，且供电与接收属于一维的探测模式，分辨率不高。

发明内容

本发明的目的是为了克服上述现有技术的不足，提出了一种动电源点供电、阵列电极接收的基于多同性源电极阵列电阻率的地下工程超前探测方法，适用于隧道或坑道等地下工程钻爆法或TBM施工过程中的不良地质探测。利用多同性源供电的方式来有效压制地下工程工作面后方异常体的干扰，并抑制工作面近处异常体的敏感性，从而相对提高了工作面远处异常体的敏感性，阵列式测量电极的接收模式提高了对异常体的空间分辨率。

本发明采用的技术方案为：

一种多同性源电极阵列电阻率的地下工程超前探测方法，它包括以下步骤：

（1）在地下工程工作面上布置多条平行的测线，在各测线上设置数量相同、间距为d的测量电极M，构成阵列式测量电极系；

（2）沿地下工程工作面轮廓线、在阵列式测量电极系的外围均匀布多个供电电极A，然后以这些供电电极A为出发点在边墙上布置供电测线，并在每条供电测线上再布置若干间隔相同的位于同一工作断面上的供电电极A，从而组成供电电极系；同时，将另一个供电电极B与测量电极N都置于工作面后远方；

（3）利用位于工作面的供电电极A以及远方的供电电极B提供正负电流，利用阵列式测量电极系中的每个测量电极M分别测量电势，并获取各测量电极M与测量电极N间的电势差；然后依次选择工作面后方下一工作断面上的供电电极A提供相同大小的正电流，并继续测量阵列式测量电极系的电势，上位机计算相应电势差；不断重复本步骤，直至供电电极系移动到预定位置，探测结束；

（4）上位机采用基于光滑约束的最小二乘线性反演迭代处理，将步骤（3）中测得的阵列式测量电极系的电势差数据进行反演，得到工作面前方地质体的三维电阻率图像。

所述步骤（4）中，反演目标函数如下：

Φ＝(Δd-AΔm)^T(Δd-AΔm)+λ(CΔm)^T(CΔm)    （1）

式中，Φ为反演目标函数，Δd为观测数据的差值，A为偏导数矩阵，Δm为模型参数修正量，C为光滑度矩阵，λ为拉格朗日常数，决定了光滑约束的权重；

对公式（1）进行变分之后得到具有光滑约束的三维电阻率反演方程：

(A^TA+λC^TC)Δm＝A^TΔd    （2）

m^(k+1)＝m^(k)+Δm    （3）

式中，m^(k)为第k次迭代中的模型参数向量，m^(k+1)为下一次反演迭代中的模型参数向量。

进行反演成像时，首先给模型参数赋初值并进行正演计算，将正演得到的理论观测数据与实际测量数据进行比较，如果方差小于用户的设定值，则输出模型参数；