基于模型试验的矿山法地铁隧道支护结构病害检测与识别

基于模型试验的矿山法地铁隧道支护结构病害检测与识别 基于模型试验的矿山法地铁隧道支护结构病害检测与识别

许 越， 于师建*， 张亚栋

(山东科技大学矿业与安全工程学院， 山东 青岛 266590)

摘要：为提高矿山法地铁隧道支护结构病害检测与识别的准确性，采用基于时域有限差分的数值模拟方法，建立常见形状空洞、充水空洞、不同填充介质脱空以及钢格栅拱架、H型钢架下空洞模型，模拟探地雷达波在不同病害中的传播规律，并通过青岛地铁3号线探测实例验证。结果表明： 空洞病害雷达反射波通常呈双曲线形状，反射信号强，但上三角形空洞图像与其他空洞图像有较大差异，表现为倾斜平行的同相轴特征；当空洞或脱空内部含水时，空洞或脱空下部出现明显反射波；钢结构雷达波反射强烈，会对钢结构下方病害探测产生干扰，在实际探测中，应将初期支护与二次衬砌分开探测。

关键词：矿山法隧道； 支护结构病害； 探地雷达； 时域有限差分法； 正演模拟

0 引言

地铁在缓解城市交通压力中发挥着重要的作用，但受隧道设计、施工及运营过程中的一些因素影响，矿山法地铁隧道会出现不同程度的支护结构病害，影响隧道的安全运营和使用寿命[1-3]。探地雷达技术作为一种新的勘探方法，目前广泛应用于地球物理探测方面。利用电磁波在介质传播过程中强度、路径以及波形等的变化可对矿山法地铁隧道支护结构内部进行判断，具有效率高、样本丰富以及无破坏性等特点[4-5]。

在隧道支护结构病害检测和识别方面，许多学者进行了相关研究。吕高等[6]对公路隧道支护结构缺陷进行FDTD正演分析，比较不同填充材料缺陷区的探测成像效果； 张明臣等[7]总结新疆寒区隧道支护结构缺陷特点； 张帆等[8]对隧道支护结构层位特征识别进行模拟研究； 颜培岩等[9]采用探地雷达技术研究暗挖隧道支护结构中钢筋分布特征； 舒志乐等[10]正演模拟隧道支护结构内的方形空洞； 杜良等[11]对隧道支护结构中基本空洞形状进行研究，并建立物理模型进行验证。已有研究多侧重于一种病害，且对三角形、弧形等空洞、空洞或脱空充水以及钢拱架下存在病害的研究较少。本文在前人研究的基础上，依托青岛地铁矿山法隧道支护结构病害检测实践，对矿山法隧道支护结构中常见病害进行分析。

1 正演模拟的时域有限差分法理论基础

基于Maxwell方程的时域有限差分法(FDTD)具有广泛的适用性、通用性、计算时间短和空间占用小等特点，是一种主要的电磁场模拟方法[12]。Maxwell方程组概括了宏观电磁场的基本规律，是支配宏观电磁现象的一组基本方程。FDTD方法是由微分形式的Maxwell旋度方程出发进行差分离散得到的[13]，运用空间离散方式，把带有时间变量的Maxwell旋度方程转化为一组差分方程，并在时间轴上逐步推进求解空间电磁场，然后由电磁问题的初始值及边界条件逐步推进求得以后各时刻的空间电磁场分布[14]。

2 隧道支护结构病害雷达探测模型试验

2.1 模型参数

根据工程实际分别建立模型，使用雷达的子波主频为900 MHz，激励源采用Ricker子波源，空间网格步长为2.5 mm×2.5 mm，模型具体参数见表1。

表1 模型参数

Table 1 Parameters of model

介质电导率/(s/m)相对介电常数电磁波速/(m/ns)混凝土0.0180.120钢筋1063000.001空气010.300砂土1.4×10-430.150岩石0.150.100水0.1×10-2810.033

2.2 支护结构常见形状空洞图像特征

受隧道施工工艺的影响，在初期支护和二次衬砌结构中经常出现近似球形或长方体形的空洞，在二次衬砌模板与模板交界处，也经常出现近似三棱锥形空洞。空洞处支护结构受力不均，易发生开裂、充水、掉块和整体失稳等灾害； 因此，空洞是隧道建设质量和运营维护检测中重点探测的病害之一。

建立模型模拟隧道支护结构中常见的空洞形状，模型大小为2.7 m×0.65 m，由混凝土和空气2种介质组成，不同形状空洞病害模型及其正演扫描结果见图1。可以看出： 1)圆形空洞雷达图像呈明显的双曲线形状，空洞的上下边界处出现2个明显的反射波峰，且能量衰减较快； 2)下三角形空洞的雷达图像也表现为双曲线形状，但并不完整，在右侧出现缺失，空洞内部反射波较多； 3)弧形空洞雷达图像与下三角形空洞雷达图像相似，但线型比较完整； 4)上三角形与下三角形的雷达图像有较大的差异，主要因为上三角形斜边为第1反射界面，当电磁波传播至此界面时发生第1次反射，反射波能量较强，且由于倾斜的原因，使得每一道扫描的电磁波传播至界面的时间近似相等，在倾斜上形成近似平行的同相轴，即在水平距离2.1～2.3 m 形成了倾斜平行的同相轴，但偏离模型中上三角形空洞的位置(水平距离为2.3～2.5 m)。三角形上顶点位置出现双曲线图像，其位置与模型相对应。左侧顶点两侧也出现了一定程度的双曲线图像。从正演扫描图中还可以看到，同相轴加强的位置一般为空洞缺陷的开始位置。

(a) 模型图

(b) 正演扫描图

图1 不同形状空洞病害模型及其正演扫描图

Fig. 1 Models of different shaped hollows and their forward modeling

2.3 支护结构空洞含水及在钢筋影响下的图像特征

隧道支护结构空洞很多是充水的，在工程检测中，准确区分含水空洞和不含水空洞有着很重要的作用。建立包含6个矩形空洞的模型，左侧2个模拟不同大小的空洞，中间2个模拟半充水和完全充水空洞，右侧2个模拟不同钢筋网格间距下的空洞。左侧第1个空洞大小为0.15 m×0.05 m，其余空洞大小为0.15 m×0.1 m。支护结构空洞充水模型及其正演扫描结果见图2。

(a) 模型图

(b) 正演扫描图

图2 支护结构空洞充水模型及其正演扫描图

Fig. 2 Models of hollows contained with water and their forward modeling

从图2可以看出： 1)左侧2个空洞处存在明显的反射波，其雷达图像与图1(b)弧形空洞雷达图像相似，都表现为较明显的双曲线形状，且随着空洞厚度的增加，高频部分出现了分层情况，表明上下界面的反射波逐渐分开； 2)半充水空洞雷达图像双曲线顶点位置下部反射波强度强于上部的反射波，说明空洞下部出现了一个介电常数差异更大的反射界面，即为模型中的水-空气反射界面，由于垂直距离较近，反射波的子波相互叠加，不能确定充水量大小； 3)完全充水空洞雷达图像上不仅在洞顶出现了反射波，而且在洞底下方也出现了明显的杂乱反射波，这是因为水的介电常数为81，其电磁波传播速度为0.033 m/ns，小于空气中的电磁波传播速度，天线接收反射波和折射波的走时较大； 4)分析右侧2个不同钢筋网格间距下空洞的雷达图像可以得知，钢筋网的雷达图像呈明显的弧形，钢筋的间距越小，图像排列越紧凑，也越明显。此外，钢筋网会对下方空洞的分辨造成一定影响，钢筋间距越减小，这种影响越强烈。在实际探测时，往往需要采取一些方法(如预测反褶积法、一致性消除法等)来消除钢筋的多次波影响。

2.4 不同填充介质脱空图像特征

连续的空洞会形成脱空，脱空病害主要出现在支护结构层与层之间，而且脱空层里面的填充体并不完全为空气，会出现泥土甚至水等介质，直接影响支护结构的稳定性[2]。建立模型模拟二次衬砌和初期支护之间的脱空，设置3个连续的含不同介质的脱空层，厚度均为10 cm，其模型及正演扫描结果见图3。从雷达图像上可明显看到脱空层的位置，砂土脱空层位置处出现明显的水平反射波图像。由于水对电磁波的强烈吸收作用，含水脱空层上部位置处的反射波能量比其他填充体的反射波能量小，但在含水脱空层的下部出现多次反射波。由于空气的介电常数与砂土的基本相同，因此其雷达图像与砂土脱空层雷达图像基本类似。

(a) 模型图

(b) 正演扫描图

图3 不同填充介质脱空模型及其正演扫描图

Fig. 3 Models of support structure gaps filled with different dielectrics and their forward modeling

2.5 钢格栅拱架下空洞图像特征

钢格栅拱架通常由4根花纹钢焊接而成，是初期支护中主要的受力部分，但在混凝土喷射过程中，由于受到钢筋的阻力及混凝土的流动性影响，易在拱架后方出现空洞。因此，建立模型模拟初期支护中的钢格栅拱架以及其后方的矩形、三角形空洞。模型大小为2.5 m×0.5 m，设置3个钢拱架，拱架之间间距为75 cm。组成钢拱架的钢筋直径为22 mm，间距为15 cm。前2个拱架下面分别设置矩形空洞和三角形空洞； 第3个不设空洞，作为对照。钢格栅拱架下方空洞模型及其正演扫描结果见图4。可以看出： 钢筋检测图像为典型的双曲线形状，矩形空洞处出现了同相轴水平平行的反射波，且两侧有双曲线图像； 三角形空洞处出现了同相轴倾斜平行的反射波图像。虽然2个空洞的反射图像受到钢筋的干扰，但仍然在扫描图中很具特征地显现出来。

(a) 模型图

(b) 正演扫描图

图4 钢格栅拱架下方空洞模型及其正演扫描图

Fig. 4 Models of hollows behind steel grating arch and their forward modeling

2.6 H型钢架后方空洞图像特征

当隧道开挖遇到V级围岩，且围岩位于挤压强烈的断裂带内时，采用H型钢架进行初期支护效果较好，同时悬挂钢筋网，H型钢架间隔50～100 cm。建立模型的大小为2.2 m×0.6 m，2个H型钢架间距为80 cm，钢架后方为间距20 cm的钢筋网，由于挂网钢筋的阻力，在钢架之间形成一定大小的喷射混凝土空洞，形状为弧形，其模型及其正演扫描结果见图5。可以看出： H型钢架的图像表现为双曲线形状，顶端出现平直的反射波，下方反射波振荡叠加； 钢筋网表现为典型的双曲线图像，间隔为20cm，与模型预设值相同； 弧形空洞位置处图像与图1(b)中弧形空洞图像相似，但由于钢筋反射波的屏蔽作用，空洞底部无法准确识别。

(a) 模型图

(b) 正演扫描图

图5 H型钢架下方空洞模型及其正演扫描图

Fig. 5 Models of hollows behind H-shaped steel arch and their forward modeling

3 工程实例

在对矿山法地铁隧道支护结构病害检测FDTD正演分析的基础上，结合青岛地铁3号线某区间的检测实例，分析实际工程中的支护结构病害雷达检测图像。青岛地铁3号线是青岛第1条地铁线路，西起青岛火车站，终到青岛北站，是连接青岛主城区的主要交通干线，总长约24.8 km，有车站22座。为了保证隧道的施工安全以及质量，在施工过程中采用探地雷达对支护结构质量进行检测。采用美国GSSI公司研制的SIR-30E型探地雷达主机，仪器为新一代双通道、高分辨率且体积轻便的专业型地质透视仪，天线为900 MHz的单体屏蔽天线，探测深度可达0.7 m，可以满足检测的需求。

3.1 仪器参数设定

结合现场实际情况对仪器参数进行设置，工作中选用的技术参数如下： 1)天线频率,中心频率为900 MHz的收发一体式屏蔽天线； 2)工作方式,连续剖面法； 3)测量方式,测量轮长度测量； 4)采样点数,512； 5)数据方式,波动成像、灰度显示。

3.2 数据处理

现场数据采集后，使用GSSI公司提供的雷达图像专业处理软件Radan7进行处理和分析。Radan7软件处理数据的主要流程为： 原始数据—调整测量方向—核定测量标记—切除废段记录—水平距离均衡—0点校正—数字滤波—偏移处理—能量均衡—文件标注—时深转换—输出雷达深度剖面图。在实际处理数据的过程中，并不是严格按照上述步骤进行，而宜根据实际测量的数据进行选择。

3.3 检测结果分析

选取一些具有代表性的检测图像进行分析。图6(a)为初期支护质量较好的雷达探测图像，雷达剖面显示地层反射波能量均匀，无超挖、欠挖和空洞形成的异常反射波，混凝土与围岩密实性好，无脱空现象； 图6(b)为典型的初期支护不密实雷达探测图像，反射波同相轴发生畸变，呈现多个细密的弧形，表示混凝土局部不均匀变化、密实性较差； 图6(c)为典型的空洞缺陷图像，根据上文空洞模型试验结果可知，此处为圆形空洞或者较小的矩形空洞，且空洞中没有充水； 图6(d)为钢格栅结构雷达图像，表现为典型的双曲线反射图像，两两一组，个别钢格栅间距较大或较小，但整体均匀，符合设计要求。

(a) 完整初期支护图像

(b) 初期支护不密实图像

(c) 空洞缺陷图像

(d) 钢格栅结构图像

图6 典型雷达检测图像

Fig. 6 Typical radar reflection images

4 结论与建议

在空洞病害的模型试验中，同相轴加强的位置点一般为空洞缺陷的开始位置。空洞探地雷达图像基本呈双曲线形，其中： 矩形空洞图像中存在水平形状波形； 三角形空洞表现为连续倾斜同相轴波形特征； 充水空洞以及脱空层中介质含水时，由于水对电磁波的吸收作用会产生能量较弱的反射波，但是在扫描图下方会出现反射波的多次振荡；当空洞位于钢格栅拱架以及H型钢架下方时，受钢结构的屏蔽影响，会影响空洞形态和大小的判定，因此应尽量在二次衬砌施工前对初期支护质量进行检测。在探地雷达数据分析方面，建议今后结合数值模拟的图谱特征，编写数据处理程序，实现探测数据的反演，得到支护结构内部的图像，而不仅仅是通过图像特征去判断支护结构的完整情况，从而更好地指导隧道支护结构质量评估。

参考文献(References)：

[1] 张顶立,张素磊,房倩,等.铁路运营隧道支护结构背后接触状态及其分析[J].岩石力学与工程学报,,32(2): 217-224.(ZHANG Dingli, ZHANG Sulei, FANG Qian, et al. Study of contact state behind tunnel lining in process of railway operation and its analysis[J].Chinese Journal of Rock Mechanics and Engineering, , 32(2): 217-224.(in Chinese))

[2] 周文轩,李菁,郭辉.探地雷达在隧道衬砌无损检测中的应用[J].施工技术,,43(11): 85-89. (ZHOU Wenxuan, LI Jing, GUO Hui. The application of ground penetrating radar in nondestructive detecting of tunnel lining[J]. Construction Technology, , 43(11): 85-89. (in Chinese))

[3] 王红良,张运良,安永林,等.地质雷达无损检测在隧道二次衬砌中的应用[J].铁道工程学报,(11): 73-76.(WANG Hongliang, ZHANG Yunliang, AN Yonglin, et al. The application of nondestructive testing using radar in tunnel lining[J]. Journal of Railway Engineering Society, (11): 73-76.(in Chinese))

[4] 谢勇勇,廖红建,昝月稳.探地雷达检测铁路路基病害的二维正演模拟[J].浙江大学学报(工学版)，,44(10): 1907-1911.(XIE Yongyong, LIAO Hongjian, ZAN Yuewen. Two-dimensional forward simulation of railway roadbed disease based on ground penetrating radar technique[J].Journal of Zhengjiang University(Engineering Science), , 44(10): 1907-1911.(in Chinese))

[5] 陈承申.探地雷达二维有限元正演模拟[D].长沙: 中南大学,.(CHEN Chengshen. Two-dimension finite element forward simulation for ground penetrating radar model[D].Changsha: Central South University,.(in Chinese))

[6] 吕高,李宁,刘新星,等.公路隧道衬砌缺陷几何形态及填充物FDTD正演分析[J]. 岩石力学与工程学报,,33(7): 1415-1423.(LYU Gao,LI Ning,LIU Xinxing, et al. FDTD forward modeling of geometric shape and fillings of lining defects of highway tunnel[J].Chinese Journal of Rock Mechanics and Engineering, , 33(7): 1415-1423.(in Chinese))

[7] 张明臣,魏燧,蒋斌松,等.地质雷达在寒区隧道衬砌质量检测中的应用研究[J].现代隧道技术,,53(1): 187-191，201.(ZHANG Mingchen, WEI Sui, JIANG Binsong,et al. On the application of ground penetrating radar for quality detection of tunnel linings in cold regions[J].Modern Tunnelling Technology, , 53(1): 187-191，201.(in Chinese))

[8] 张帆,张立,廖红建,等.基于空气耦合式探地雷达探测的衬砌层位特征识别的正演模拟[J].隧道建设,,35(2): 121-127. (ZHANG Fan, ZHANG Li, LIAO Hongjian, et al. Forward simulation for identification of tunnel lining layer characteristics based on detecting by air-coupled ground penetrating radar[J].Tunnel Construction, , 35(2): 121-127.(in Chinese))

[9] 颜培岩,曾新霞,刘柳,等.探地雷达技术在暗挖隧道钢筋间距检测中的正演模拟与应用[J]. 特种结构，,32(6): 90-95. (YAN Peiyan, ZENG Xinxia, LIU Liu, et al. Forward simulation and application of ground penetrating radar technology in detection of the rebar space

of the bored tunnel[J]. Special Structures, , 32(6): 90-95.(in Chinese))

[10] 舒志乐,刘新荣,朱成红,等. 隧道衬砌空洞探地雷达三维探测模型试验研究[J]. 岩土力学,,32(增刊1): 551-558.(SHU Zhile, LIU Xinrong, ZHU Chenghong,et al. Study of model test about 3D GPR detection of tunnel lining cavity[J]. Rock and Soil Mechanics, , 32(S1): 551-558. (in Chinese))

[11] 杜良,刘树才,梁棋念,等. 隧道检测中基本形状空洞探地雷达图像特征[J].地球物理学进展,,30(6): 2892-2896.(DU Liang, LIU Shucai, LIANG Qinian, et al. GPR image feature of basic shape cavity in tunnel inspecting [J]. Progress in Geophysics, , 30(6): 2892-2896. (in Chinese))

[12] 冯德山,戴前伟,何继善，等.探地雷达GPR正演模拟的时域有限差分实现[J].地球物理学进展,,21(2): 630-636.(FENG Deshan, DAI Qianwei, HE Jishan,et al. Finite difference time domain method of GPR forward simulation [J].Progress in Geophysics, , 21(2): 630-636. (in Chinese))

[13] Goldman D. Ground-penetrating radar simulation in engineering and archaeology[J].Geophysics,1994,59: 224-232．

[14] Yee K S. Numerical solution of initial boundary value problem involving Maxwell equation in isotropic media[J]. IEEE Trans Antennas Propagate, 1996, 14(3): 302-307．

Detection and Recognition of Support Structure Diseases of Metro Tunnel Constructed by Mining Method Based on Forward Modeling

XU Yue, YU Shijian*, ZHANG Yadong

(College of Mining and Safety Engineering, Shandong University of Science and Technology, Qingdao 266590, Shandong, China)

Abstract： The service life of Metro tunnel constructed by mining method would be affected by cracks and corrosion of and hollows behind support structure. In order to improve the detection and recognition accuracy of support structure disease of Metro tunnel constructed by mining method, models of common shaped hollows, hollows contained with water, support structure gaps filled with different dielectrics and hollows behind steel grating arch and H-shaped steel arch are established by finite-difference time-domain method(FDTD). The propagation rules of ground penetrating radar wave in tunnel support structure are analyzed; the reflected results of different support structure diseases are analyzed; and a case study is made on engineering detection of Line No. 3 of Qingdao Metro. The results show that: 1) The radar reflection waves of hollows usually show as hyperbolic curves and their signals are strong. The radar images of triangle hollows which show as parallel inclined lineup are greatly different from those of other shaped hollows. 2) When the hollows and support structure gaps contained with water, the radar reflection waves are clearly visible at bottom of hollows and support structure gaps. 3) The steel structures of support structure have great influence on the detection and recognition of tunnel structure diseases; as a result, the diseases detection of primary support and secondary lining should be carried out respectively.

Keywords： tunnel constructed by mining method; tunnel support structure disease; ground penetrating radar; finite-difference time-domain method; forward modeling

收稿日期：-05-19;

修回日期：-12-13

第一作者简介：许越(1990—)，男，内蒙古兴安盟人，山东科技大学安全技术及工程专业在读硕士，研究方向为灾害预测与防治。E-mail： xuyuechn@。*通讯作者： 于师建，E-mail： ysj7179@。

DOI：10.3973/j.issn.1672-741X..02.009

中图分类号：U 45 文献标志码： A

文献标志码：A

文章编号： 1672-741X()02-0179-06