地面物探技术在隧道掘进破碎带区段的应用

2024-03-11朱伟国何雄斌

物探化探计算技术 2024年1期

朱伟国,何雄斌

(广东省地质物探工程勘察院,广州 510800)

0 引言

广佛东环城际珠三角城际轨道交通网是广佛同城化的骨干线路,是推动和实现珠江三角洲区域经济一体化的重要交通基础设施,对提高广州国际综合交通枢纽的地位,促进各个区域经济一体化和提升整体竞争力有着十分重要的意义。全线长约46.75 km,途径广州市番禺区、海珠区、黄埔区、天河区和白云区,共设大石、科学中心等8座车站。沿线需经过地面建筑物密集区、穿越已建成公路、铁路等交通线路,线路经过的地质情况较为复杂,隧道经过洞身地层变化较大,局部呈软硬夹层状,受区内地质构造影响,岩石裂隙普遍发育,局部漏水严重,使得隧洞施工风险较大。为了保证隧道施工安全以及预防地质灾害的发生,有必要开展超前地质预报的工作[1-3]。

超前地质预报主要分为地质调查、地质钻探及物探技术等方式。其中物探技术作为地质预报的重要组成部分需要采用综合物探技术来提高资料解释的可靠性[4-6]。物探技术依据与探测对象的空间位置不同,可分为地面物探技术和隧洞物探技术。地面物探技术依据探测地球物理场的不同,分为浅层地震波法、电磁法及直流电法等几大类[7]。

大太区间隧道在往小里程方向进行TBM掘进施工时,在隧洞开挖至里程K38+823时,掌子面顶部整体出现了塌方,局部区域有大石块掉落现象发生,地表出现了轻微塌陷,对隧洞的施工进展产生了一定的影响。为了探明隧洞上方的不良地质体、断层以及破碎区等地质情况,为隧洞进一步的施工挖掘提供了基础依据,本次应用地面物探技术进行综合探测。针对此次隧道超前预报的地形条件、探测深度要求及周边的干扰情况等因素的考虑,采用了高密度电阻率法和瞬变电磁法两种地面物探技术进行综合探测。实践证明:高密度电法具有成本低、效率高等优点;瞬变电磁法具有灵活、探测深度大等优点;两者相结合能基本探测出地下不良地质体的分布情况,是一种有效的地面物探技术组合。

1 地质概况及地球物理特征

1.1 工程地质概况

拟建大太隧道区间位于广州市白云区太和镇。场区以低山丘陵为主,地形局部陡峭,植被茂密。区内断层和褶皱均较发育,北东向的广从断裂构造组成了本区的主要构造,断裂生成于加里东运动,在海西-印支构造阶段控制着广花凹陷的形成。燕山晚期至喜马拉雅早期对龙归盆地的形成和演变起一定的控制作用,也是区域控岩、控热结构。早期多表现为逆断层,晚期多表现为正断层。由于广从断裂带在该区域与其它方向断裂交接切错[8-9],受区内地质构造影响,使得该区域的岩性变化相对复杂,岩石裂隙普遍较发育。

如图1为隧道地质断面图,从图上可知,隧道所在标高在-10 m～-20 m之间。测区位于石炭系板岩及震旦系片麻岩的转换区域,该区域岩性相对复杂,且岩性风化不均匀。隧洞上方为主要为中风化板岩,前期钻孔揭露该地层由于受区内地质构造影响,岩石裂隙普遍较发育。

图1 隧道地质断面图Fig.1 Geological section of tunnel

1.2 地球物理特征

根据场地收集的钻孔资料表明,测区内第四系冲积层厚度一般小于20 m,土层及岩石的电阻率变化见表1。第四系冲积层视电阻率为30 Ω·m ～100 Ω·m;残积土层视电阻率为50 Ω·m～300 Ω·m;碳质板岩为相对中阻一般为300 Ω·m～1 000 Ω·m;片麻岩视电阻率较高一般为1 000 Ω·m～8 000 Ω·m;不良地质体由于岩石破碎一般为70 Ω·m～240 Ω·m。可见,工区岩土层存在电性差异,为开展高密度电法勘探与瞬变电磁法提供了较好的物性基础。

表1工区电阻率物性参数表Tab.1 Parameter of resistivity for various media

表2 物探异常对应情况表Tab.2 Correspondence of geophysical anomalies

2 方法原理

2.1 高密度电阻率法

高密度电阻率法是以地下被探测目标体与周围介质之间的电性差异为基础,不良地质体因岩性破碎一般表现为低阻,导电能力好,而完整岩石电阻率普遍高,导电能力差,应用人工在地下建立电场,观测和分析人工电场的地下分布规律,在断面上通过电阻率反映不同地质体的导电能力以达到了解岩土分层、基岩起伏形态以及地质构造、软弱夹层等不良地质体分布范围的一种地球物理勘探方法[10-11]。高密度电阻率法实际上是一种快速高效的阵列勘探方法,野外工作时,将数十个电极一次性布设完毕,通过程控式多路电极转换器选择不同的电极组合方式,即可完成野外数据的快速采集[12],见图2。

图2 高密度电法野外工作原理示意图Fig.2 Schematic diagram of high-density resistivity method in the field

2.2 瞬变电磁法

瞬变电磁法是利用不接地回线等方式向地下发送脉冲式一次电磁场,用线圈或接地电极观测由该脉冲电磁场感应的地下涡流产生的二次电磁场的空间和时间分布,来解决有关地质问题的时间域电磁法。瞬变电磁法的探测深度取决于探测信号的采样时刻,早期信号反映了浅表地质异常,晚期信号携带了深层地质信息。不同于频率域电磁法的全程采样方式,瞬变电磁法对数据的采集通常始于发射电流完全关断之后,理论上,瞬变电磁法的这种特性可以避免发射电流激发的一次场响应对目标体涡流产生的二次场响应的干扰,有利于缩小信号的动态范围,提高探测的分辨率[13-14]。

本次采用重庆璀陆探测技术有限公司生产的FCTEM60-1拖曳式高分辨瞬变电磁系统,由瞬变电磁主机、收发一体线圈、数据处理与成像软件组成。采用了“恒压钳位”高速线性关断和无损消互感技术,结合高密度高动态信号采集,具有极强的浅层和高分辨探测能力;拖曳式观测模式极大提高了工作效率;大磁矩发射增强了抗干扰能力,改善了深部探测效果;丰富的线圈配置提供了从浅表到深部,地面、水上、水下和隧道探测的灵活选择[15-17]。

3 数据采集与处理方法

3.1 测线布置

为了探测出掘进位置上方的岩层地质情况,以当前掘进位置为中心,布置2条物探测线, L1和L2,由于线路区域经过陡峭边坡,无法做到与线路完全垂直布设,测线结合实际地形进行了调整,测线布置如图3所示。

图3 工区测线布置图Fig.3 Layout of geophysical prospecting lines in the work area

3.2 参数选择

3.2.1 高密度电阻率法

本次高密度电阻率法采用重庆地质仪器厂生产的DUK-2A型高密度电法仪对场地进行探测,为了兼顾到有效深度和探测精度,选用四极温纳测深装置。本次勘查的参数选择如下:电极距a=3 m,最小隔离系数n=1,最大隔离系数n=20,供电时间为1 s。采用正负方波供电,最高供电电压400 V。

3.2.2 瞬变电磁法

FCTEM60-1拖曳式高分辨瞬变电磁系统是采用收发一体的重叠回线装置。经现场测试采用最优参数如下:测点距为5 m,局部异常地段加密为2.5 m,发送电流60 A,发送频率16 Hz。叠加次数的选择主要与当地的噪声水平有关,理论上叠加次数越大,采集到的衰减曲线信号信噪比越高,经对比试验,本次叠加次数选200次,并重复观测两次取信噪高的数据。

3.3 数据处理

3.3.1 高密度电阻率法

本次高密度电阻率法数据处理软件使用骄佳技术公司的高密度电法数据与极化率层析成像软件Geogiga RTomo6.0。数据处理步骤为:数据格式转换→坏点数据剔除→数据滤波处理→测点高程输入→地形改正→最小二乘法反演→输出反演电阻率断面→物探解释→地质解释剖面输出→完成。

3.3.2 瞬变电磁法

本次瞬变电磁法资料处理采用FCTEM数据处理软件。将原始数据进行预处理后,导入处理软件进行反演。数据处理步骤:在数据处理之前,先建立新工区,输入工区名称,文件保存路径,对工区和测线情况进行简单描述,最后选择数据导入。数据导入成功后,左边显示每一点的衰减曲线,右边显示测线多测道图,如图4。然后选择一维优化反演算法进行反演,先进行初始模型建立,选择高程数据和保存路径,建立初始模型和迭代次数,最后点击确定,进行反演计算。反演结束后,可以查看每个测点的曲线拟合情况,将反演结果导入相应的出图软件格式进行成图。

图4 数据界面示意图Fig.4 Schematic diagram of data interface

4 资料解释与讨论

4.1 L1测线异常总体特征及解释

图5(a)为L1线高密度电法二维反演解释剖面图,在点号1 120～1 140下部存在条带状低阻异常,视电阻率在200 Ω·m～400 Ω·m左右,异常两侧为相对高阻,电阻率在1 000 Ω·m以上,推测该低阻异常为断裂破碎带的反映,编号Fw1;在点号1 170～1 195下部存在条带状低阻异常,视电阻率为200 Ω·m ～400 Ω·m之间,两侧为相对高阻,电阻率均在1 000 Ω·m以上,推测为断裂破碎带的反映,编号Fw2;点号1 110～1 170下部分布低阻异常,视电阻率小于50 Ω·m,推测该区域岩性相对破碎,为基岩风化不均所致。隧道经过位置在点号1 137～1 172之间,浅部为相对高阻,推测为该区域基岩埋深稍浅,上部岩性较完整,下部存在相对低阻,表明隧道上方的基岩岩性存在风化不均匀的现象。

(a)L1线高密度电法二维反演剖面图 (b)L1线瞬变电磁反演剖面图图5 L1线综合物探反演剖面图Fig.5 comprehensive geophysical inversion profile of line L1

图5(b)为L1线瞬变电磁法反演剖面图,测线区域的岩性层状分布比较明显,可分为覆盖土层及强、中风化基岩层,表层电阻率在100 Ω·m～300 Ω·m,平均厚度约为15 m,与前期钻孔揭露的情况基本一致。岩层电阻率在300 Ω.m～1 100 Ω·m,局部存在相对低阻,推测为断裂破碎及岩性不均匀的反映。在点号1 120～1 350段对应标高75 m的范围,电阻率相对较低,两侧电阻率表现为相对高阻,推断该区域受断裂影响,基岩较为破碎,编号为Fw1;在点号1170～1185段对应标高40 m～65 m的范围,电阻率相对较低,两侧电阻率表现为相对高阻,推断该区域受断裂影响,基岩较为破碎,编号为Fw2;隧道范围对应测点号1 135～1 175之间,上方的电阻率总体变化平缓,隧道上方电阻率显示较为完整,仅在点号1 143～1 153附近对应标高22 m～34 m之间、点号1 167～1 178附近对应标高6 m～48 m之间,存在相对低阻,推测为基岩破碎区域。

4.2 L2测线异常总体特征及解释

如图6(a)为L2线高密度电法二维反演剖面图,测线区域视电阻率高、低阻异常相间,等值线平缓变化,浅部视电阻率在50 Ω·m ～200 Ω·m左右,主要反映覆盖土层,中深部电阻率大于1 000 Ω·m,主要反映基岩。高低阻相间则反映岩石风化不均。设计隧道经过水平位置在1 132～1 166号点之间,浅部的电阻率相对低,反映风化土层,下部为相对高阻,反映岩石。从剖面图可知测线区域存在三个明显的等值线下凹或局部低阻异常区,推测为岩体塌陷或节理裂隙集中发育区。

如图6(b)从瞬变电磁二维反演电阻率剖面图上看,表层电阻率在100 Ω.m～300 Ω·m,岩层电阻率在300 Ω·m～1 100 Ω·m。在测线范围标高在10 m～40 m之间电阻率变化比较大,推测该区域的岩性分布不均;隧道范围对应在L2线的点号1 135～1 165之间,在点号1 120～1 130附近对应标高13 m～23 m之间、点号1 164～1 175附近对应标高6 m～22 m之间,存在相对低阻,推测为基岩破碎区域。由于基岩破碎区域距离隧道顶板较近,不排除由于隧道掘进过程中,震动引起该区域的破碎基岩产生掉块现象,因此对该区域的异常应予以重视。

4.3 物探异常综合解释与分析

通过高密度电阻率法及瞬变电磁法均可得到测线L1、L2地层结构、断裂破碎带、节理裂隙发育区的分布情况。从两种方法的反演剖面可知,高密度电阻率法和瞬变电磁法反映覆盖土层平均厚度基本一致,断裂破碎带的空间展布也能一一对应。

结合物探异常对应情况表可知,高密度电法异常G1与瞬变电磁法异常T1、T2对应、高密度电法异常G2与瞬变电磁法异常T3对应、高密度电法异常G3与瞬变电磁法异常T4、T5对应、高密度电法异常G4与瞬变电磁法异常T7对应。其中高密度电法异常G3、瞬变电磁法异常T5与已知塌陷区域水平位置相吻合,证明了两种方法均能较好反映出地面沉降区域。

从异常特征上看,高密度电法反映出来的低阻异常范围较大,且受场地限制的情况下无法外延,使得深部的地质情况存在缺失,而反映同一区域的瞬变电磁法异常则更为精细,且有效探测深度更大。所以瞬变电磁法能够弥补高密度电法探测深度不足的问题。