多同性源电极阵列电阻率的地下工程超前探测方法及系统的制作方法
【专利摘要】本发明涉及一种基于多同性源电极阵列电阻率的地下工程超前探测方法，适用于隧道或坑道等地下工程钻爆法或TBM施工过程中的不良地质探测：（1）在工作面上构建阵列式测量电极系；（2）在边墙上组成供电电极系；同时，将另一个供电电极B与测量电极N都置于工作面后远方；（3）利用位于工作面的供电电极A以及远方的供电电极B提供正负电流，利用阵列式测量电极系中的每个测量电极M分别测量电势，并获取工作面及各断面中各测量电极M与测量电极N间的电势差；不断重复本步骤，直至供电电极系移动到预定位置，探测结束；（4）采用基于光滑约束的最小二乘线性反演迭代处理，将测得的电势差数据进行反演，得到工作面前方地质体的三维电阻率图像。
【专利说明】多同性源电极阵列电阻率的地下工程超前探测方法及系统
【技术领域】
[0001]本发明涉及隧道或者坑道等地下工程钻爆法或TBM施工过程中不良地质的探测方法，尤其涉及一种基于多同性源电极阵列电阻率的地下工程超前探测方法及系统。
【背景技术】
[0002]随着水利水电工程和交通工程的建设重点逐步向地形地质条件极端复杂的西部山区和岩溶地区转移，所遇到的工程地质条件异常复杂，施工中面临着突水、突泥等诸多高风险地质灾害。突涌水灾害给施工安全带来了重大灾难和生命财产损失，而预先探明隧道掌子面前方的地质灾害赋存源是规避和防治重大灾害事故的有效途径，因此开展隧道超前预报工作研究显得尤为迫切和重要。
[0003]作为一种重要的工程地球物理探测技术，直流电阻率法以其适应性强、方便快捷、成本低等诸多优点被广泛应用于岩土工程勘察、滑坡地质灾害监测、环境岩土工程评价等领域中。与地面直流电阻率法相比，由于隧道掌子面空间狭小，导致地面的很多探测方法无法应用于隧道突涌水灾害源超前预报。目前直流电阻率法在隧道(巷道)超前预报中大都采用三极法的探测模式，即将供电点A置于掌子面，供电点B置于隧道后方“无穷远”处，测量电极M、N沿着隧道轴向移动测量。
[0004]这一方法主要面临以下问题:①由于测量电极MN在掌子面后方移动测量，易受到掌子面后方异常体的干扰；②由于隧道的特殊环境，隧道三极法的探测模式测线垂直于异常体，导致对掌子面近处的异常体响应敏感，对远处的异常体响应相对不敏感，且供电与接收属于一维的探测模式，分辨率不高。
【发明内容】
[0005]本发明的目的是为了克服上述现有技术的不足，提出了一种动电源点供电、阵列电极接收的基于多同性源电极阵列电阻率的地下工程超前探测方法，适用于隧道或坑道等地下工程钻爆法或TBM施工过程中的不良地质探测。利用多同性源供电的方式来有效压制地下工程工作面后方异常体的干扰，并抑制工作面近处异常体的敏感性，从而相对提高了工作面远处异常体的敏感性，阵列式测量电极的接收模式提高了对异常体的空间分辨率。
[0006]本发明采用的技术方案为:
[0007]一种多同性源电极阵列电阻率的地下工程超前探测方法，它包括以下步骤:
[0008]( I)在地下工程工作面上布置多条平行的测线，在各测线上设置数量相同、间距为d的测量电极M，构成阵列式测量电极系；
[0009](2)沿地下工程工作面轮廓线、在阵列式测量电极系的外围均匀布多个供电电极A，然后以这些供电电极A为出发点在边墙上布置供电测线，并在每条供电测线上再布置若干间隔相同的位于同一工作断面上的供电电极A，从而组成供电电极系；同时，将另一个供电电极B与测量电极N都置于工作面后远方；
[0010](3)利用位于工作面的供电电极A以及远方的供电电极B提供正负电流，利用阵列式测量电极系中的每个测量电极M分别测量电势，并获取各测量电极M与测量电极N间的电势差；然后依次选择工作面后方下一工作断面上的供电电极A提供相同大小的正电流，并继续测量阵列式测量电极系的电势，上位机计算相应电势差；不断重复本步骤，直至供电电极系移动到预定位置，探测结束；
[0011](4)上位机采用基于光滑约束的最小二乘线性反演迭代处理，将步骤(3)中测得的阵列式测量电极系的电势差数据进行反演，得到工作面前方地质体的三维电阻率图像。
[0012]所述步骤(4)中，反演目标函数如下:
[0013]Φ = ( Δ d-Α Δ m)τ ( Δ d-Α Δ m) + λ (CAm)T (CAm) (I)
[0014]式中，Φ为反演目标函数，Ad为观测数据的差值，A为偏导数矩阵，Am为模型参数修正量，C为光滑度矩阵，λ为拉格朗日常数，决定了光滑约束的权重；
[0015]对公式(I)进行变分之后得到具有光滑约束的三维电阻率反演方程:
[0016](ΑΤΑ+λ CtC) Λm = At Λ d (2)
[0017]m(k+1) = m(k)+Am (3)
[0018]式中，m(k)为第k次迭代中的模型参数向量，m(k+1)为下一次反演迭代中的模型参数向量。
[0019]进行反演成像时，首先给模型参数赋初值并进行正演计算，将正演得到的理论观测数据与实际测量数据进行比较，如果方差小于用户的设定值，则输出模型参数；
[0020]否则，计算偏导数矩阵A、光滑度矩阵C、以及理论观测数据与实际测量数据的差值Ad，代入方程(2)求得每次反演迭代中的模型参数修正量Am，再代入方程(3)得到下一次迭代中的模型参数m(k+1)，直至正演的理论观测数据与测量数据的方差小于设定值，则输出模型参数，从而通过反演得到工作面前方地质体的三维电阻率图像。
[0021]一种多同性源电极阵列电阻率的地下工程超前探测系统，它包括布置在地下工程工作面上的多条平行测线，在各测线上设置数量相同、间距为d的测量电极M，构成阵列式测量电极系；
[0022]在地下工程工作面轮廓线以及阵列式测量电极系的外围均匀布多个供电电极A，以这些供电电极A为出发点在边墙上布置供电测线，并在每条供电测线上再布置若干间隔相同的位于同一工作断面上的供电电极A，从而组成供电电极系；
[0023]同时，将另一个供电电极B与测量电极N都置于工作面后远方；
[0024]所述各供电电极均与电源连接，所述各测量电极则与上位机连接。
[0025]所述平行测线有3条，分别为U、L2和L3，它们布置在距离底部j//、I丑和*丑的高度上，其中H为拱顶距离底部的高度。
[0026]所述供电电极B与测量电极N置于工作面后方距离工作面无穷远处，例如100-150m处，此处的100-150m，是无穷远位置，大于5倍的AM即可，IOOm和150m都是大于的，一般认为结果是没有差异的。
[0027]本发明具有以下有益效果:
[0028](I)利用多同性源供电可以有效地压制地下工程工作面后方异常体的干扰，并能够降低掌子面近处异常体的敏感性，从而使得工作面远处异常体的敏感性相对提高；
[0029](2)采用阵列式测量电极的接收模式，有效地提高了对异常体的空间分辨率，增大了探测距离。
[0030](3)采用多个同性源供电电极，与传统三极法相比，降低了掌子面近处异常比，提高了掌子面原处的异常，如表1所示。对于掌子面近处的异常体(4m)，同性源异常比13.1%小于传统定点源异常比16.1%，降低了掌子面近处异常体的敏感性；对于掌子面远处的异常体同性源异常比均大于传统定点源异常，从而提高了对远处异常体的探测能力。
[0031]表1不同探测距离异常对比
【权利要求】
1.一种多同性源电极阵列电阻率的地下工程超前探测方法，其特征是，它包括以下步骤: (1)在地下工程工作面上布置多条平行的测线，在各测线上设置数量相同、间距为d的测量电极M，构成阵列式测量电极系； (2)沿地下工程工作面轮廓线、在阵列式测量电极系的外围均匀布多个供电电极A，然后以这些供电电极A为出发点在边墙上布置供电测线，并在每条供电测线上再布置若干间隔相同的位于同一工作断面上的供电电极A,从而组成供电电极系；同时，将另一个供电电极B与测量电极N都置于工作面后远方； (3)利用位于工作面的供电电极A以及远方的供电电极B提供正负电流，利用阵列式测量电极系中的每个测量电极M分别测量电势，并获取各测量电极M与测量电极N间的电势差；然后依次选择工作面后方下一工作断面上的供电电极A提供相同大小的正电流，并继续测量阵列式测量电极系的电势，上位机计算相应电势差；不断重复本步骤，直至供电电极系移动到预定位置，探测结束； (4)上位机采用基于光滑约束的最小二乘线性反演迭代处理，将步骤(3)中测得的阵列式测量电极系的电势差数据进行反演，得到工作面前方地质体的三维电阻率图像。
2.—种多同性源电极阵列电阻率的地下工程超前探测系统，其特征是，它包括布置在地下工程工作面上的多条平行测线，在各测线上设置数量相同、间距为d的测量电极M，构成阵列式测量电极系； 在地下工程工作面轮廓线以及阵列式测量电极系的外围均匀布多个供电电极A，以这些供电电极A为出发点在边墙上布置供电测线，并在每条供电测线上再布置若干间隔相同的位于同一工作断面上的供电电极A，从而组成供电电极系； 同时，将另一个供电电极B与测量电极N都置于工作面后远方； 所述各供电电极均与电源连接，所述各测量电极则与上位机连接。
3.如权利要求2所述的基于多同性源电极阵列电阻率的地下工程超前探测系统,其特征是，所述平行测线有3条，分别为LpL2和L3,它们布置在距离底部、P和去开的高度上，其中H为拱顶距离底部的高度。
【文档编号】G01V3/06GK103645514SQ201310729095
【公开日】2014年3月19日 申请日期:2013年12月25日 优先权日:2013年12月25日
【发明者】刘斌, 李术才, 聂利超, 杨为民, 王传武, 田明禛, 陈磊, 宋杰, 刘征宇, 李树忱, 王世睿, 李尧 申请人:山东大学