沧榆高速顶管区地质雷达探测分析.pdf

中 国 煤 炭 地 质COAL GEOLOGY OF CHINA第 35 卷 8 期2023 年 8 月Vol.35 No.8Aug.2023沧榆高速顶管区地质雷达探测分析卢斌荧（陕西省交通规划设计研究院有限公司，陕西西安，710065）摘要：在高速公路运营阶段，利用地质雷达探测并分析次生灾害的案例并不多见，尤其是在探测公路顶管穿越区的次生灾害方面。在道路顶管施工中，不可避免出现次生灾害，影响公路的正常运营。为便捷的识别路基以下的次生地质灾害对路网的影响，本研究应用地质雷达技术，以特定的顶管为中心，两侧各15米，分布多条测线，探测道路病害，经过治理后复测，对探测出的病害路段再次治理后探测。关键词：地质雷达；次生灾害探测；公路顶管穿越区中图分类号：P642文献标识码：ADetection and Analysis of GPR in Cangyu Highway Pipe Crossing AreasLU Binying（Shaanxi Provincial Transport Planning Design and Research Institute Co.，Ltd.，Xi an，Shaanxi 710065）Abstract:During the operation phase of highways，there are not many cases of using ground penetrating radar to detect and analyzesecondary disasters，especially in detecting secondary disasters area where highway top pipes pass through.During road jackingconstruction,secondary hazards inevitably occur,affecting the routine operation of the highway.In order to conveniently identify theimpact of secondary geological hazards below the roadbed on the road network,this study applies GPR technology to detect roaddiseases by distributing multiple survey lines centered on a specific pipe jacking,15 meters on each side,and after treatment and re-testing,and detecting the detected diseased sections after treatment again.Keywords:ground penetrating radar(GPR);secondary disaster detection;highway pipe-crossing area地质雷达是工程技术人员探测地下目标的一种常用手段。20世纪 50年代，地质雷达开始用于冰层、冰川厚度的探测，因其携带方便、快速准确等的优势，广泛应用工程验收等无损检测，在公路隧道探测、岩土工程勘察等领域，能够直观、实时的显示工程地质情况。随着计算机技术的突破性发展，地质雷达性能越来越强大，地质雷达探测在理论研究、硬件设备制造、探测原理、数据的解译和分析处理方面取得了重大突破，天线覆盖范围越来越广，探测深度和精度得到了极大拓展。随着技术的不断更新发展，其应用领域更加广泛，社会经济效益显著。在高速公路勘察设计阶段，为隧道的围岩分级，判明裂隙发育区、富水区的位置提供了参考，也为探明隧道施工区的地质情况，钻孔的布设、断层的分布提供了基础数据。然而，在高速公路运营阶段，地质雷达探测次生灾害治理应用的案例并不多见。地质雷达（GPR），是一种利用高频无线电磁波来确定介质内部物质分布规律的地球物理探测方法，在地下传播过程中，通过对电磁波回波的研究和分析，就可以获得地下目标体的几何形态与空间位置1-3。姜化冰对深度不同、尺寸各异的方形空洞进行正演模拟，分析各种因素空洞目标雷达图谱的变化4；李世念等利用基于三维时域有限差分方法原理的GprMax 3D软件，充气型和充水型空洞三维探地雷达正演模拟，对数据体不同方向的图像特征进行分析，并结合工程实例，验证了正演模拟结果的可靠性5；李政用Matlab语言编写模拟程序，对隧道衬砌检测和超前预报中的典型地电模型进行正演模拟，取得了正演图像，提高图像的解释水平6；韩俊涛建立了常见铁路路基病害模型和检测的图谱判识基本规律，圆形空洞呈双曲线特征，矩形空洞在边界上呈现一段水平圆弧7；刘胜峰通过不同高度、圆形和方形两种规则空洞雷达探测实验，结论是雷达对水平尺度探测效果良好，但是垂直大小的精度逊于水平方向8；陈婕通过GprMax计算原理，建立了不同大小，不同深度、不同形状的不同填充介质的混凝土空洞模型，总结出其成像的规作者简介：卢斌荧（1972），男，博士，高级工程师，研究方向：遥感地质、地理信息系统。E-mail：收稿日期：2023-05-06责任编辑：张燕生doi:10.3969/j.issn.1674-1803.2023.08.12文章编号：1674-1803（2023）08-0070-058 期卢斌荧：沧榆高速顶管区地质雷达探测分析律9。韩佳明总结出地下圆形空洞雷达扫描图像对称轴左侧为单调递减的连续凹函数，厚度由不断增大至稳定10；对称轴右侧为单调递增的连续凹函数，厚度由稳定至不断减小；随着埋深与半径增大，曲线曲率减小，曲线图像趋于缓和，曲线的特征体现出张开弧度增大的趋势。1研究背景1.1研究区概况本项目位于榆林市，场地地形高低起伏，地面高程在1 149.801 342.80m。场地地貌单元主要分为风砂滩地地貌和黄土梁地地貌，风砂滩地地貌，表层风积砂较厚，低洼处常年积水形成海子，滩地大小不一，多曲折相连，黄土梁地多被风沙半覆盖，梁地大小不一。沿线地表植被稀疏，主要生长为黄蒿、杨树等沙生植物。在区域范围内地基土自上而下分层描述如下：层细砂（Q4eol）：褐黄色，含云母片，稍湿饱和，松散中密状态，以稍密状态为主，摇震反应中等，无光泽反应，干强度低，韧性低。层厚 5.508.50m，层底深度为 5.508.50m，土质均匀，稍湿，松散-中密，以稍密状态为主，中压缩性，工程性能一般，可选做管道和工作坑地基持力层。层黄土（Q3eol）：黄褐色红褐色，饱和，可塑软塑，以可塑状态为主，针状孔隙发育，含钙质结核，局部含粉土夹层，偶见蜗牛碎片，无湿陷性，属中压缩 性 土。层 厚 4.505.30m，层 底 深 度 为 10.6013.00m，土质均匀，稍湿，坚硬硬塑，以硬塑状态为主，中压缩性土，无湿陷性，工程性能一般，可选做管道和工作坑地基持力层。层粉土（Q3eol）：褐黄色，含云母片，饱和，中密，含钙质结核，局部钙质结核富集，摇震反应中等，无光泽反应，干强度低，韧性低。未揭穿，最大揭露厚度11.40m，土质均匀，稍湿，中密密实转台，以中密状态为主，中压缩性，工程性能一般，可选作管道和工作坑地基持力层。地层情况的复杂程度，影响到了地质雷达数据分析的精准程度，本研究充分收集研究区的地层情况，便于后续研判。1.2区域灾害治理实验概况本次实验根据探测结果提供的数据作为灾害治理设计的依据，对探测出的缺陷（脱空和不密实）区域进行治理，治理后运行一段时间后再次探测，对仍存在灾害缺陷，进一步治理，直至完全探测不出缺陷为止，满足公路运营的需要。2数据来源和技术方法2.1数据来源本研究收集到高速公路和顶管的施工阶段的资料，包括地质勘察资料、路基和路面结构图、顶管施工图，详细整理了路面及其以下的包括路面结构层、路基、底基层的结构，底基层以下的地质钻探资料。这些资料都作为后续分析的参考资料。2.2地质雷达检测原理地质雷达技术（Ground Penetrating Radar，简称GPR）利用主频为106109Hz波段的电磁波，以宽频带短脉冲的形式，由地面通过地质雷达天线发射器发送到地下，电磁波经地下目的体或地层反射界面后返回，被地质雷达天线接收器所接收，通过对收到的雷达信号进行处理和解译，达到探测地下目标体的目的11-12。地表下的介质比较复杂，电磁波吸收和反射的存在差异，通过对反射波形的分析，可以建立地下介 质 的 结 构 模 型（地 质 雷 达 检 测 原 理 示 意 图见图1）。图1 地质雷达检测原理示意Figure 1 Schematic of GPR detection principle71中国煤炭地质第 35 卷电磁波在介质中传播，电场强度分量关系式如下：Ex=E0e-kx（1）式中：Ex为传播距离Z=x处的电场强度；E0为传播距离Z=0处的电场强度；k为传播系数，其表达式为k=+j，其中为衰减常数，表示每单位距离衰减程度的常数；j为电流密度；为相移常数，表示每单位距离落后的相位11。=221+()2-1（2）=221+()2+1（3）式中：为电磁波的角频率（=2f，f为电磁波中心频率）；为介质的磁导率；为介质的电导率；为介质的介电常数。实际计算中，通常用相对磁导率（r=/0）和相对介电常数（r=/0）计算介质的衰减常数和相移常数。=20022+()02-（4）=20022+()02+（5）式中：0=410-7H/m，为空气磁导率：0=1/3610-9F/m，为空气介电常数。在介质中传播的平面电磁波，在遇到不同的波阻抗（）界面时将发生反射和透射，产生反射波和透射波。垂直极化平面波的反射系数：R=2cosi-1cosn2cosi+1cosn（6）平行极化平面波的反射系数：R=1cosi-2cosn1cosi+2cosn（7）式中：i为电磁波入射角；n为电磁波反射角；1和2为两个不同介质的波阻抗。（1）高阻介质波阻抗公式如下式表示：=（8）在垂直入射情况下，1=n=0，则：R=1-21+2（9）R=2-11+2（10）上式表明，介质介电常数的不同与反射信号相关。对于高阻介质，反射信号的大小主要取决于介质介电常数的差异。（2）低阻介质波阻抗公式如下式表示：=s2wm（11）在垂直入射情况下，i=n=0，则：R=1-21+2（12）R=2-11+2（13）上式表明，对于低阻介质，反射信号的大小主要取决于介质电导率的差异11。介质的复杂程度与反射信号表现出高关联性，这种特性可以进行地表下的地质分析。3数据采集与分析3.1天线选取与设备标定本次检测采用地质雷 SIR-4000 主机及配套100MHz和400MHz屏蔽天线，连续采集，主机时窗选用200ns、AD采样16位、1 024扫描样点数、128扫描线/s、9节点自动增益、3道水平平滑；垂向高通滤波50MHz、低通滤波300MHz、3道水平平滑。首先对检测目标进行介电常数标定，检测前对路基的介电常数做现场标定，每处顶管位置标定1次，每次采集3个数据，根据已知顶管位置计算出路基介电常数值，取平均值为该段路基的介电常数。经现场标定，本次检测段路基介电常数取值为8。3.2测线布设根据探测需要，在顶管下穿区域沿车道方向布设测线，测线布置覆盖顶管穿越区域（地质雷达检测测线布设示意图见图2）。3.3地质雷达数据处理本 次 的 地 质 雷 达 数 据 采 用 数 据 处 理 软 件RADAN7进行处理。处理流程：数据输入数据编辑能量均衡数值滤波时深转换图形编辑输出剖面图（地质雷达主要处理方法及功能见表1）。地质雷达在检测过程中，由于受周围环境电磁噪声、仪器自身噪声等因素的影响，数据记录中除目标体信息外，还包含其他许多干扰信息。这些干扰信息降低了数据的信噪比，往往会掩盖真实异常，728 期卢斌荧：沧榆高速顶管区地质雷达探测分析图2 地质雷达检测测线布设示意图Figure 2 Schematic diagram of GPR detection line layout甚至造成假异常，使检测结果不准确，因此在利用地质雷达数据进行检测结果解释之前，需要进行数据处理来压制干扰，提高数据的信噪比。对检测结果进行处理可以提高信噪比，但是要得到检测结果，还需要结合现场调查情况及相关构筑物结构信息和现场记录，尽量剔除假异常，使地质雷达检测信息和目标体真实情况相对应，进而实现地质雷达数据的解释。地质雷达检测数据的处理，需要先找出雷达波图谱中的反射波同相轴不连续、产生弯曲、反射波能量强、回波振幅反应较强的区域，再根据雷达波的相位、频率和幅值变化等因素综合分析判定，最终确定病害分布情况。表1 地质雷达主要处理方法及功能Table 1 Main processing methods and functions of GPR处理对象波形较平的噪声高频噪声（像雪花一样的波形）展宽了的波形绕射波（压缩了的双曲线）、校正反射层位置增强低振幅波的能见度分辨细微异常特征生成清晰的显示文件方法水平高通滤波、垂直高通滤波、空间滤波、背景修复水平低通滤波、垂直低通滤波、空间滤波反褶积运算偏移算法函数、改变显示增益Hilbert变换、空间滤波算法函数、层位拾取、静校正3.4实验过程本次实验段落为3条煤矿疏干水综合利用项目下穿高速工程管线区域（下穿深度一般大于10m，不含路基填高），对下穿段高速公路道路土体密实度进行探测，防止施工引起道路路基土体松散、空洞，对行车安全造成安全隐患，及时掌握探测段内高速公路路基密实度情况，以便及时采取相应措施，快速处理以消除安全隐患。经过第一次探测，发现仅1号管线周围有5处脱空和1处不密实（表2），其余两条管线区域正常；对1号管线区域发现的不密实和脱空区域经过专业设计和施工治理，进行第二次探测和治理（表3）；对于仍然存在不密实的区域进行于第三次探测，探测结果未发现不密实现象。实验数据比较发现：第一次探测出5处灾害缺陷，经过治理后第二次探测依然有三处，分别是测线237 不密实区域经过治理后缺陷大小由 2.70m2.58m缩小为 1.13m0.65m，缩小至 10%，测线 242的不密实区域经过治理后的缺陷大小由 4.80m1.72m缩小为1.91m0.5m，缩小至11%，测线245的脱空的区域经过治理后缺陷大小由1.33m2.19m缩小为 1.13m1.58m，缩小至 60%（由脱空变为不密实）。对这些不密实的区域治理后进行第三次探测，没有发现不密实等灾害缺陷，公路营运正常，实验数据真实有效。73中国煤炭地质第 35 卷表2 第一次地质雷达探测典型病害缺陷统计Table 2 Statistics of typical disease detected by GPR at the first time序号12345线路沧榆高速测线编号230237242245246缺陷类型不密实不密实不密实脱空不密实严重程度一般一般一般一般一般可信程度疑似疑似疑似疑似疑似起始位置9.956.605.257.9310.8终止位置10.769.3010.059.2615.71起始深度/m4.052.681.961.221.67终止深度/m7.765.263.683.414.72缺陷大小/mm0.813.712.702.584.801.721.332.195.533.04表3 治理后第一次地质雷达探测病害缺陷统计Table 3 Statistics of disease detected by GPR at the first time after treatment序号123线路沧榆高速测线编号311333323缺陷类型不密实不密实不密实严重程度一般一般一般可信程度疑似疑似疑似起始位置5.0113.8212.45终止位置6.9214.9513.58起始深度/m1.251.131.30终止深度/m1.751.782.88缺陷大小/mm1.910.501.130.651.131.583.5实验结果判定对地质雷达检测数据进行处理可以提高信噪比，但是要得到检测结果，还需要结合现场调查情况及相关构筑物结构信息和现场记录，尽量剔除假异常，使地质雷达检测信息和目标体真实情况相对应，进而实现地质雷达数据的准确解释。主要判定特征应符合下列要求：1）密实：信号幅度较弱，甚至没有界面反射信号。2）不密实：界面的强反射信号同相轴呈绕射弧形，且不连续，较分散。3）脱空：界面反射信号强，三振相明显，在其下部仍有强反射界面信号，两组信号时程差较大。本次探测以特定的顶管为中心，两侧各15m，分布多条测线，经两次探测后治理，第三次探测无异常，时间间隔分别为1个月、2个月和4个月，通过近一年的运营观察，历经夏天高温和冬天严寒的运行实践，未发现明显的路基沉降的变化，证明此方法可在高速公路顶管穿越区指导公路运营养护部门对道路次生灾害的检测，快速确定治理方案，保证公路的正常运行。4结论通过地质雷达法对顶管施工区域高速公路路基检测，可成功排除道路存在的安全隐患，保障高速公路运营安全。经过施工前后地质雷达检测数据分析对比，发现地质雷达法在道路病害检测中发挥出了明显的优势，总结如下：1）地质雷达法应用于运营阶段的高速公路道路探测，无损、快捷、对交通运营影响较小。2）地质雷达法对高速公路道路病害探测效果显著，通过病害整治前后地质雷达数据对比，发现道路病害范围明显缩小及无异常，证明地质雷达在道路病害探测中数据的准确性及方法的可靠性。3）地质雷达在探测过程中干扰因素较多，对数据处理、分析、解译的要求较高，需要详细收集道路结构、地层地质资料，以及现场实际情况，有助于提高数据的准确性。4）地质雷达应用于高速公路病害探测，能有效排除道路安全隐患，保障高速公路运营安全，在今后的勘察探测中，将会发挥出更大的价值。参考文献：1刘澜波，钱荣毅探地雷达：浅表地球物理科学技术中的重要工具 J 地球物理学报，2015，58（8）：2606-2617.2陈义群，肖柏勋论探地雷达现状与发展 J 工程地球物理学报，2005，2（2）：149-155.3胡群芳，郑泽昊，刘海，等三维探地雷达在城市市政管线渗漏探测中的应用 J 同济大学学报：自然科学版，2020，48（7）：972-981.4姜化冰探地雷达地下空洞目标探测研究 D 黑龙江哈尔滨：哈尔滨工业大学，20165李世念，王秀荣，林恬，等基于GprMax的道路空洞三维探地雷达正演数值模拟 J 中国地质灾害与防治学报，2020，31（3）：132-138.6李政探地雷达在公路隧道中的应用研究 D 广西南宁：广西大学，2014.7韩俊涛铁路路基病害的地质雷达正演模拟与探测方法研究D 北京：北京交通大学，20188刘胜峰地质雷达应用于公路隧道衬砌无损检测的实验研究D 湖南长沙：长沙理工大学，2007.9陈婕基于GprMax3.0的探地雷达混凝土空洞目标正演模拟与试验分析 D 湖南长沙：长沙理工大学，2018.10韩佳明，刘明明，牛宇凯，等.地下圆形空洞地质雷达成像机理J.科学技术与工程，2022，22（23）：10130-10137.11刘英利.地质雷达在工程物探中的应用研究 D 四川成都：成都理工大学，2008.12彭珍发.西商高速公路建设中的地质雷达技术应用与研究D 四川成都：成都理工大学，2011.74