1. 引言
矿井井筒的破裂、井壁的空鼓以及钢筋的移位都给矿井的安全生产带来了严重的威胁,全国已经有多个井筒发生井壁破碎的事故。此类事故不仅影响矿井的正常生产,还会引发更严重的安全隐患,造成巨大的经济损失 [1] [2]。对于井筒损伤的治理工作,通常是采用地面注浆或者井壁直接注浆的方法,但是由于井筒的损伤类型和损伤范围难以确定,导致注浆位置也难以确定,不能有效保证注浆效果,经常出现注浆效果不理想甚至失败的情况 [3] [4] [5]。所以,如何确定井筒的损伤程度以及注浆位置,对于矿井井筒损伤的防治工作至关重要。
鉴于矿井井筒损伤位置较浅、深度较小,且井壁全部由钢筋混凝土浇筑而成,导致传统的矿井物探方法在实际应用中存在一定的局限性 [6],其中矿井瞬变电磁法受关断效应的影响不可避免的在浅部存在盲区,给浅部资料解释带来困难;高密度电法体积效应大,分辨率较低,同时需要良好的接地条件,无法在钢筋混凝土井壁上进行数据采集。而地质雷达对浅部地层介质具有较高分辨率,同时其抗干扰能力强,被广泛应用于结构物无损检测、公路路面塌陷探测等勘察领域 [7] [8] [9]。因此,本文采用地质雷达方法,对北洺河铁矿电梯井井壁的损伤情况进行了详细的检测,对后期的治理工作提供了基础数据,为该矿井的安全生产提供了技术保障。
2. 地质雷达原理
地质雷达(Ground Penetrating Radar,简称GPR)是一种利用高频电磁波(1 MHz~1 GHz)的反射来探测地下不同地层分界面或者目标体的物探方法 [10]。该方法是以宽频带短脉冲形式激发,通过发射天线(T)将高频电磁波送入地下,经地下地层分界面或目标体反射后由接收天线(R)接收,如图1所示。地下地层界面或目标体的反射波的走时为 [11]:
(1)
其中,z是目标体的埋藏深度,v是地下各层介质电磁波的速度,x是发射天线和接收天线之间的距离。

Figure 1. Principle of reflection detection
图1. 反射探测原理
当地下介质中的电磁波速度v为已知时,就可以根据采集到的时间t值求出地层界面或目标体的深度。在剖面测量中,x值一般是固定不变的。图2是对于一个简单的地质模型,雷达探测记录的波形示意图。在波形记录图上各测点均以测线的铅垂方向记录波形,最终构成雷达剖面。

Figure 2. Schematic diagram of the waveform recorded by radar detection
图2. 雷达探测记录的波形示意图
当矿井井筒井壁存在破裂、空鼓以及钢筋的移位时,井壁后的地层完整性也会遭到破坏,此时其介电常数与围岩存在明显差异。在地质雷达测量剖面上,反射波信号幅度会发生变化,导致同相轴不连续或者局部连续,这也是地质雷达探测地下地层分界面或目标体的地球物理基础 [12] [13]。
3. 数据采集
3.1. 地质概况
在矿区内,均为第四系黄土及河床卵石层所覆盖。据钻孔所见,矿区内地层(由新至老)如下 [14]:1)第四系黄土砾石层:上部黄土层厚1~1.5 m。分布于河床两侧。中部砾石层厚度10~143.21 m,在基岩上部有一层厚2~10 m的红色粘土层。2) 石炭、二叠系含煤岩系:主要分布在矿区东端,为呈NNE走向,倾向南东,缓~中倾斜的一套含煤岩系地层,属云驾岭煤田矿区范围。3) 中奥陶系灰岩层:地层由于受后期岩浆岩的侵入影响,岩石大都变质,重结晶作用明显。包括奥陶系中统上马家沟组一段(O22-1)、奥陶系中统上马家沟组二段(O22-2)、奥陶系中统上马家沟组三段(O22-3)、奥陶系中统峰峰组一段(O23-1)以及奥陶系中统峰峰组二段(O23-2)。在垂向上,矿区灰岩上覆有第四系砂砾卵石层、砂质粘土砾石层和底部粘土层,下为燕山期火成岩托底。
矿区内的主要控矿构造为次级的背斜,主要发育在奥陶系中统马家沟石灰岩中部O22-2层位中,在背斜轴部有一宽200~300 m的断裂带,其宽度向上递减,为闪长玢岩、磁铁矿、矽卡岩及角砾岩所充填。到O23-1岩层,由于原岩为角砾状灰岩,岩石疏松破碎,又受后期褶皱影响,岩层产状紊乱,层间滑动剧烈,岩石破碎及糜棱岩化发育。至上部O23-2岩层中,岩层产状渐趋平缓。
矿区含水层主要有奥陶系中统石灰岩含水层、石炭二迭系砂岩薄层灰岩含水层、第四系砂砾石含水层以及闪长岩裂隙含水层。第四系底部粘土层、煤系地层泥页岩和未经风化、构造破坏的完整闪长岩侵入体为矿区主要隔水层。
北洺河铁矿电梯井井筒采用钢筋混凝土井壁结构,井筒长为5.3 m,宽3.15 m,混凝土井壁厚0.3 m。在电梯运行过程中,发现电梯井马头门侧井壁有渗水、局部淋水、钢筋凸起和混凝土脱落现象,如图3所示。为了查明井筒的破裂、井壁的空鼓以及钢筋的移位等情况,为后续加固工作提供依据,需要对电梯井井壁进行地质雷达探测。
3.2. 仪器参数
对于本次井筒损伤探测,使用国产的KON-LD(A)型地质雷达系统。该系统由主机、天线、控制显示单元、数据采集处理软件及配件组成。不同频率的天线,可以满足不同探测深度的要求。由于电梯井井筒内金属架较多,所以为了避免电磁波的多次反射,采用了屏蔽式藕合天线。由于不同天线反映的地质细节深度不同,所以为了能够全面将该矿电梯井井壁的破裂、空鼓以及钢筋的移位探查清楚,综合考虑探测深度和精度,对于同一侧线,分别使用500 MHz和1 GHz天线,采样时窗为100 ns,采样率为1024。
3.3. 探测施工
根据井筒实际情况,在电梯井筒马头门一侧,布置4条测线,测线间距0.5 m,探测起点为−200 m深度处,结束点位−227.9 m深度,测线布置如图4所示。在现场探测时,为了方便施工,工作人员站在电梯轿厢顶,将雷达的发射和接收装置也放在轿厢顶,随后电梯匀速下降,探测点也匀速往下移,进行连续测量。为了方便记录以及解释位置与实际位置对应,在数据采集过程中,每隔5 m进行一次标定,两组天线顺序进行。
4. 数据处理及解释
4.1. 数据处理
由于矿井井筒内金属架、电梯等各种干扰,导致电磁波信号与原始信号发生变化,产生差异,所以在数据处理过程中需要对接收到电磁波反射信号进行处理,数据处理主要包括以下几个方面 [15] [16] [17]:1) 道间平均处理来压制干扰波;2) 滤波处理,用于突出有效波,压制和消除干扰波、多次波;3)信号的增益和衰减处理,用于校正信号损失;4) 速度分析,为时深转换提供速度数据。通过以上处理,可以有效突出异常区域,最大程度的降低电梯井内各种金属干扰的影响,最终实现井筒损伤的精细解释。
4.2. 数据解释
图5~8分别是测线1~测线4不同天线雷达探测的深度剖面图,其中图(a)为500 MHz天线所得结果,反映的是井壁后方4 m范围内的破碎情况;图(b)为1 GHz天线所得的结果,反映的是井壁附近2.5 m范围内空鼓、钢筋的移位情况,图中椭圆圈定的部位是可能存在缺陷的部位。
从雷达图像中波形幅度以及同相轴的连续性来看,在图5(a)中存在多处反射波振幅变弱、同相轴不连续以及反射界面不明显的情况。其中,在−204.3~−204.5 m位置处,深度在0.4~0.9 m处推断有一个小型的裂隙;在−210~−211.5 m的排水口处,深度0.6~1.4 m推断有一裂隙。在−216~−217.5 m位置处,深度在0.4~0.6 m处推断有一小型空洞。在−224~−226.8 m位置处,深度0.4~1.2 m处推断有一小型裂隙群。从图5(b)可以看出,在井壁附近整体缺陷、空鼓和张裂情况较多,特别是在0.3~0.7 m深度,在−203~−205.5 m,−206~−207.5 m,−207.5~−211 m,−214~−214.5 m,−215~−216 m,−217.5~−222 m,−224.5~−227 m范围内反射波同相轴不连续以及反射界面不明显,推断空鼓情况严重。但通过观测,钢筋整体深度分布均匀,没有较大的位移,说明墙面膨胀较为一致。
从雷达图像中波形幅度以及同相轴的连续性来看,在图6(a)中存在多处反射波振幅变弱、同相轴不连续以及反射界面不明显的情况。其中,在−203.5~−204.5 m位置处,深度在0.3~0.7 m处推断有一个小型的裂隙;在−207.5~−214 m位置处,深度0.3~0.7 m推断为一破碎带。在−217.5~−224 m位置处,深度在0.6~0.9 m处推断有一破碎带。在−224~−226.5 m位置处,深度0.2~0.6 m处推断有一小型破碎。从图6(b)可以看出,该测线与测线1情况基本一致:在井壁附近整体缺陷、空鼓和张裂情况较多,特别是在0.3~0.7 m深度空鼓情况严重。其中,在−212~−214.5 m,−215.5~−218 m,−221~−222 m,−224~−227 m范围内反射波同相轴不连续以及反射界面不明显,推断存在空鼓情况。但通过图中可以看出钢筋整体深度分布均匀,没有较大的位移,说明该测线墙面膨胀较一致。

Figure 5. The depth profile of Line 1 (a) the result of 500 MHz; (b) the result of 1 GHz
图5. 测线1深度剖面图(a) 500 MHz天线探测结果;(b) 1 GHz天线探测结果

Figure 6. The depth profile of Line 2 (a) the result of 500 MHz; (b) the result of 1 GHz
图6. 测线2深度剖面图(a) 500 MHz天线探测结果;(b) 1 GHz天线探测结果
从雷达图像中波形幅度以及同相轴的连续性来看,在图7(a)中存在多处反射波振幅变弱、同相轴不连续以及反射界面不明显的情况。其中,在−200.5~−201.5 m位置处,深度在0.3~0.7 m处推断有一个小型的裂隙;在−203~−204 m位置处,深度0.3~0.7 m处推断为一破碎带;在−209.5~−212.5 m位置处,深度在0.2~0.7 m处推断有一破碎带;在−213.7~−214.9 m、−215.4~−220 m和−221.5~−224.5 m位置处,深度0.3~0.9 m处推断有一小型破碎。从图7(b)中可以看出,该测线井壁附近整体空鼓、张裂情况较小,主要集中在0.3~0.7 m深度。在−201~−201.5 m,−214~−216 m以及−224~−227 m处反射波同相轴不连续以及反射界面不明显,推断以上位置存在空鼓情况。与测线1和测线2一致,该测线范围内钢筋整体深度分布均匀,没有较大的位移,说明墙面膨胀较一致。

Figure 7. The depth profile of Line 3 (a) the result of 500 MHz; (b) the result of 1 GHz
图7. 测线3深度剖面图(a) 500 MHz天线探测结果;(b) 1 GHz天线探测结果
从雷达图像中波形幅度以及同相轴的连续性来看,在图8(a)中存在多处反射波振幅变弱、同相轴不连续以及反射界面不明显的情况。其中,在−201.5~−202.3 m位置处,深度在0.4~0.8 m处推断有一个小型的空洞;在−203~−204 m位置处,深度0.3~0.7 m处推断为一破碎带;在−206.7~−208.5 m位置处,深度在0.4~0.7 m处推断有一破碎带;在−213.5~−216 m位置处,深度0.2~0.6 m处推断有一小型破碎;在−217.3~−218 m、−222.5~−223 m位置处,深度0.4~0.7 m处推断有一小型破碎;在−224.6~−225.5 m位置处,深度0.2~0.7 m处推断有一小型破碎;在−218~−222 m位置处,深度1.4~2.4 m处推断有一小型破碎区。从图8(b)中可以看出,该测线井壁附近整体空鼓、张裂情况也比较小,同样主要集中在0.3~0.7 m深度。其中,−203.9~−204.9 m,−210.2~−211 m,−226~−227.9处反射界面不明显,推断存在空鼓情况;−226~−227.9,30~40 cm处反射波同相轴不连续,推断该位置存在破碎情况。该测线范围内钢筋整体深度分布均匀,没有较大的位移。

Figure 8. The depth profile of Line 4 (a) the result of 500 MHz; (b) the result of 1 GHz
图8. 测线4深度剖面图(a) 500 MHz天线探测结果;(b) 1 GHz天线探测结果
综合4条测线探测结果分析发现电梯井马头门一侧井壁深部整体性较好,井壁后方破碎和裂隙发育较浅,最深为2.4 m,大部分在1 m以内,推断是由于地下水浸泡引起膨胀变形导致。井壁浅部附近(0.3~0.7 m)的空鼓和张裂情况在井壁右侧比较严重,越靠近左侧,浅部损伤越轻。由于破碎和裂隙发育较浅,且发育面积较小,可以直接采用锚网支护加固的方法进行治理。
5. 结论与建议
5.1. 结论
采用地质雷达探测技术对电梯井井筒损伤探测表明,损伤位置主要集中于深部(0.4~1.2 m)范围内,岩石破碎和裂隙发育情况较为严重,而浅部虽然存在一定的空鼓和张裂,但钢筋分布较为均匀,没有明显移位。通过分析认为井壁后方围岩遇水膨胀变形是造成井筒损伤的主要原因。
从整体结果来看,对于铁矿电梯井井筒损伤的探测,采用地质雷达能够满足探测精度的要求,探测结果达到了预期目的。同时,本次地质雷达探测结果,为后续井壁治理工作提供了指导和依据。
5.2. 建议
1) 由于井筒内干扰因素较多,对于数据处理和解释也带来一定的困难,这也是后续需要进一步研究解决的问题。
2) 随着开采水平的不断下移,回采错动界限逐渐向电梯井靠近,受爆破冲击及地压作用,电梯井井壁后期产生变形、开裂的几率相应增加。因此,后期需加强对井壁的监测工作。
基金项目
国家十三五重点研发计划项目(编号:2016YFC0801600)。