0 引 言
寒区高铁隧道口边坡受冻融循环作用易失稳变形
[1,2],直接威胁行车安全并增加运维成本
[3]。为此,边坡监测技术正从单点人工观测向多源自动化感知演进
[4,5]:早期主要依赖全球导航卫星系统(Global navigation satellite system,GNSS)和全站仪进行表面位移监测
[6],虽操作简便,但受限于人工频次和空间覆盖范围,对边坡快速变形的捕捉率不足50%,预警效能有限
[7]。随后,引入倾角仪、应变计及TDR等接触式传感器,实现了连续自动观测,但在寒区低温和冻融环境下稳定性不足,误报率高达10%~15%
[8]。近年来,为突破覆盖范围限制,遥感技术(如InSAR、激光扫描、无人机摄影)被应用于大范围形变识别
[9-11],毫米波雷达结合视频识别技术则实现了寒区坡面的实时监测
[12,13],国际研究趋势也聚焦于通过多源数据融合降低单一传感器的局限性
[14,15]。
然而,现有技术在实际应用中仍存在显著局限。GNSS与全站仪获取的数据点位分散、信息孤立
[16];接触式传感器在低温环境下易中断和误报;遥感及无人机监测虽覆盖范围广,但受气象和地形干扰,实时性与可靠性不足
[17]。部分工程仍依赖人工记录和数据处理,不仅效率低下,还易导致信息滞后和判断失误
[18,19]。2022年贵广高铁D1862次列车在暴雨中因泥石流侵入隧道口导致脱轨事故
[20],便暴露了现有依赖单一传感器的监测体系在潜在风险识别与实时预警能力上的不足。
多源数据融合理论上可提升监测精度,通过整合地面与空基传感器信息增强对复杂地形边坡的监测能力
[21]。例如,毫米波雷达结合视频识别可实时捕捉坡面异动,而InSAR与无人机影像则能提供大尺度形变信息,实现跨时间尺度的数据互补。但在实际应用中仍面临诸多挑战:数据来源异质,格式与精度差异大,传统融合方法处理困难;寒区低温和复杂气象易导致传感器性能下降,造成监测结果不准确或延迟;实时性与可靠性不足,难以支撑连续预警。因此,亟须构建高精度、智能化、实时化的边坡监测体系。
本文构建了集“感知-传输-分析-预警”于一体的多源融合监测系统,依托北斗高精度定位、毫米波雷达和InSAR等技术,开发数字化预警平台,实现边坡风险的智能识别与动态预警,并通过利民隧道口工程验证其可行性。
1 工程地质与隐患特征
如
图1所示,利民隧道(K68+361.61~K69+210.57)位于黑龙江省哈尔滨市阿城区,洞口边坡陡峭(高差>15 m)。该边坡基岩风化、节理发育,覆盖层松散(粉质黏土-碎石混杂),含局部软弱带。工程区属寒温带季风气候(1991-2020年年均气温为3.6 ℃,历史极低温达-35 ℃)
[22],年均降雨量523 mm,显著低于降雨型滑坡高发区域降雨量水平。现场勘测显示,该区域冻融作用显著(最大季节冻深>1.5 m),主要活动于3~5月解冻期及10~11月冻结期,是诱因滑坡灾害的主要因素。其灾害机制体现为以下3类耦合效应:
(1)冻融-水热耦合效应:当土体温度处于-5° C~0° C的相变阈区,且质量含水率>25%时,易诱发浅层滑移或坡脚剪切变形。
(2)春融期水文激增效应:融雪叠加降雨导致土体饱和,抗剪强度骤降,进而引发既有滑裂带扩展。
(3)列车动荷载扰动效应:加剧冻融期土体结构重组与孔隙水压力波动。
综合来看,滑坡灾害高敏感时段为春融期、强降雨过程后及高频运营时段。
2 多源监测预警系统构建
2.1 监测指标体系
为精准识别寒区高铁隧道口边坡在冻融、水热及动荷载耦合作用下的复杂失稳演化过程,需构建具有针对性、可量化、可联动特征的监测指标体系。
基于利民隧道工程地质条件与隐患特征,将监测对象划分为3类,即结构变形类、环境扰动类、异常响应类。据此构建多源信息指标体系(见
表1),为边坡安全状态的多维度、全过程动态感知与智能研判提供支撑。
2.2 基于有限元风险识别的监测部署
2.2.1 边坡有限元模型构建与失稳区域识别
基于利民隧道口边坡原始CAD图,提取特征点坐标,利用Surfer软件插值拟合边坡面坐标。结合工程实际情况选取关键点,并依据断面图确定土层厚度,建立边坡三维有限元模型(见
图2,其中
X为隧道轴向,
Y为垂直轴向,
Z为竖直向)。结合地质资料确定土体类型后,通过ANSYS计算分析
XY、
YZ、
XZ方向的剪切弹性应变,识别出区域1和2为潜在滑坡失稳区(见
图3)。
2.2.2 “点-线-面-体”融合感知部署策略
基于有限元分析识别的潜在滑坡失稳区(区域1、区域2,见
图3),本文提出以下分层监测部署策略(见
图4):
(1) 点监测:在区域1、区域2的关键位置(
图3中菱形/矩形标记点)部署北斗位移监测终端、测斜仪、温湿度传感器及雨量计(分别对应
图4中设备1~4),获取关键点位的实时状态信息。
(2) 线监测:通过分析区域1、区域2内北斗位移监测终端的点位数据,生成位移趋势线,评估潜在滑移面的线性变形特征。
(3) 面监测:结合区域1、区域2内点监测设备连线所勾勒的范围,并利用区域1、区域4部署的毫米波雷达/视频监控设备(分别对应
图4中设备6、7),综合分析坡体表面的形变状态。
(4) 体监测:综合区域1、区域4内毫米波雷达/视频监控数据及高分遥感影像(分别对应
图4中数据6、7),评估包含区域1、区域2在内的边坡整体稳定性状态。
2.2.3 多维度监测设备部署方案
结构变形监测:设备布设于区域1、区域2的坡顶、坡面、坡体。基准站点位需对空条件良好(视场内高度角≥15 °,且无遮挡)。
环境扰动监测:温湿度传感器布设于区域1、区域2(各布设2台),雨量计布设于隧道口区域(布设1台)。
异常响应监测:依托毫米波雷达、InSAR雷达遥感卫星、高清摄像头与高分数据影像,识别边坡特定点位/区域的树木位移、危岩落石、异物侵限、防护栅栏状态及坡面整体稳定性。该监测模块可与雷达数据联动校验,用于春融期或暴雨期间的应急响应。监测设备部署实景如
图5所示。
虽然多源传感设备的部署需投入一定成本,但通过结合边坡潜在失稳特征的精准布点、设备功能的综合利用及针对寒区环境的适应性设计,可在保障监测效果的同时合理控制成本。
2.3 预警模型与预警等级划分
2.3.1 预警模型构建
为克服冻融期突发形变实时监测难题,本文构建了基于北斗位移触发、图像复核与多源融合的三级自动化预警模型(见
图6),实现滑坡/崩塌的早期识别与快速响应。具体流程如下:
(1)一级触发(位移判定):实时分析北斗位移数据。基于莱因达准则(3σ)分析历史同期同位置数据,确定位移统计阈值。结合位移速率阈值(见
表2)综合判定:若位移量及位移速率均未超限,则终止流程;若位移量或位移速率任意一项超限,则触发二级判定。
(2)二级识别(图像复核):调用毫米波雷达图像,结合栅栏状态、侵限记录,检测活动边界、结构变形或异物侵限。若无异常,则判定误触发并终止流程;若确认异常,则进入三级判定
[23]。
(3)三级预警(多源融合):此阶段融合毫米波雷达形变信息、InSAR图斑、光学影像、视频监控与北斗位移曲线等多源数据进行综合判定。其中,InSAR数据重点关注变形速率>5 mm/年且持续增加的区域;裂缝监测依据扩展速率预警,当扩展速率>1 mm/d时,对应一级预警条件;环境因子评估则结合
表2中的寒区特征参数(如地温、含水率、裂缝扩展速率)展开。若多源数据均未超限,则终止流程;若确认存在滑移面扩展、持续变形或动态异常(符合
表2高等级预警位移速率,并伴随对应寒区特征),则触发报警。
(4)风险响应:当触发一级(红色)预警后,系统自动启动应急联动机制,向平台及责任单位发送预警(支持人工确认),同时提升监测频率强化跟踪。具体响应措施按预警等级执行(见
表2)。
本模型通过自动化决策克服人工流程的固有延迟,显著提升冻融期瞬时形变监测预警能力。
2.3.2 预警等级划分与阈值确定
本文采用我国现行边坡变形四级预警体系(绿/黄/橙/红对应初始变形至加速临滑阶段)。预警等级划分以位移速率为核心指标(见
表2,位移速率阈值范围为<0.5、0.5~1.5、1.5~2.5、>2.5 mm/d),其设定依据如下:
(1)历史数据分析:基于利民隧道口边坡冻融失稳案例库,统计各变形阶段位移速率特征。
(2)规范文献参考:综合边坡工程规范及文献建议值
[24,25]。
(3)强度折减法验证:按理论失稳位移值的70%校核阈值。
需要注意的是,位移速率阈值的应用需紧密结合寒区特殊的冻融环境进行风险评估。针对利民隧道工程,冻融循环(地温波动、水分迁移及冻胀融沉)显著影响变形演化。在相同位移速率下,若出现特定寒区环境(如冻胀期地温<-5 ℃、融沉期地温>0 ℃、高含水率、裂缝加速扩展),则边坡失稳风险将显著升高。
表2所列寒区特征参数(地温、含水率、裂缝扩展速率)源于历史数据统计,用于辅助判断位移速率对应的实际风险。例如,当位移速率达到二级预警(1.5~2.5 mm/d)时,若同步出现融沉环境(地温>0 ℃)或裂缝扩展迹象(裂缝扩展速率<1 mm/d),则采取高于该等级的响应措施,因其实际风险已逼近更高等级。
综上,通过位移速率指标论证及与寒区冻融参数的关联分析,本文构建了适用于本工程的寒区边坡四级预警模型,详见
表2所列预警等级、位移速率阈值、寒区特征参数及响应机制。
3 数字化预警平台
3.1 平台架构设计与核心实现
针对寒区高铁隧道口边坡具有的强动态性与突发性特征,本文开发的数字化预警平台构建了“天-空-地”多源集成监测预警系统,整合InSAR(天基)、毫米波雷达(空基)及地面传感器(地基)等设备(部署见2.2节),实现空间、参数与功能互补。
该平台采用“感知层-传输层-分析层-应用层”4层架构(见
图7)。
感知层集成北斗、传感器、雷达、摄像头、InSAR/高分影像,覆盖结构变形、环境扰动、异常响应监测。
传输层采用4G/5G/Wi-Fi/微波等多制式冗余通信,结合北斗CORS统一坐标,通过数据加密、断点续传、误码校验保障寒区数据传输的可靠性。
分析层/应用层通过自主开发的上位机软件实现,其核心功能包括实时多源数据可视化(图表/三维模型)、历史趋势回溯、视频异常识别、设备管理及预警规则配置。该系统集成5大模块,即系统管理、感知数据处理、数据管理、影像监控、雷达/北斗/高分数据分析(含稳定性评估算法);支持台式机与移动APP远程访问及权限管理,满足高铁工务运维需求。
3.2 核心功能模块
3.2.1 监测数据分析
寒区高铁隧道口边坡北斗位移数据的预处理采用多级滤波:首先,剔除原始报文中经纬度高程为0或位置异常的无效解,仅保留差分定位状态为4的有效数据;其次,数据经4G传输后,采用两级均值算法构建基准值(见式(1)),即先对2 h数据剔除20%极值后取均值,再对12组均值进行平均运算,最终将得到的稳定基准值存储于Flash中。
为提升数据的可靠性,采用正弦拟合模型去除温度、含水率等环境数据的周期性噪声。2020-2021年监测期内,1~4号设备的位移数据符合高斯分布(
R²>0.93)。基于莱因达3
σ准则,计算得出各设备高程位移与水平位移临界值(
μ+3
σ),见
表3。其中,2号设备的高程位移临界值最高(18.33 mm),3号设备的水平位移临界值最低(3.88 mm)。该阈值量化了寒区高铁隧道口边坡预警标准,有效降低了误报与漏报风险。
在融合层面,智能监测预警系统实现了结构数据与环境扰动的时序对齐与空间投影映射。基于正弦拟合模型,建立温度与位移响应的关联函数,即:
式中:d(t)为边坡位移(如高程位移或水平位移);为冻融循环频率;为位移与温度变化的时间滞后量;A为位移振幅,反映冻融最大响应;B为基线位移,反映长期累积趋势。其中,A和B通过实测数据的最小二乘回归联合确定。工程上,A值越大表明边坡对冻融变形更敏感,B值的变化可揭示边坡整体稳定性的演化趋势。
2020年11月至2021年5月监测数据显示,大气温度在-27~17 ℃之间波动,地表温度滞后大气7~10天且呈正弦变化;冻胀前期(11~12月)与消融期(3~5月)温度变化较快,稳定冻结期(1~2月)温度波动平缓。边坡高程位移与温度显著相关:负温区间,高程位移在-5~0 mm范围浮动;温度回升后,高程位移由-5 mm增至5 mm,其趋势与温度拟合曲线(
R²=0.941)一致。该规律验证了热胀冷缩与冻胀力对边坡变形的驱动作用,并揭示了冻融界面迁移及饱和带变化的影响机制(见
图8)。
3.2.2 三维可视化展示
智能监测预警系统集成多源数据空间分析功能(见
图9),具体如下:
(1)高分/高清影像地貌识别(见
图9(a)、(b))。
(3)InSAR大范围形变场解析(见
图9(d))。
这些功能为预警决策提供空间基准,技术人员可手动调用影像数据,辅助风险研判。
数字化平台通过多源数据自动分析与智能预警,替代人工巡检,显著降低运维工时与人力成本;其精准预警推送使现场处置由“全面排查”转向“靶向操作”,为寒区高铁隧道口边坡监测提供可推广方案,赋能同类工程。
3.2.3 智能预警发布
北斗定位监测作为智能预警系统的核心模块,通过专属界面(见
图10)实时高精度监测边坡地表位移。以1号北斗毫米终端为例,当高程位移从-40.0 mm变为50.0 mm时,系统基于预设阈值(如位移速率>0.5 mm/d触发初级预警)动态分析位移趋势。
当位移曲线界面(见
图10)检测到数据异常时,系统立即激活预警机制,弹出窗口显示异常设备、预警内容及状态判定,并依据
表2规则确定预警等级。在预警信息界面(见
图10)中,不同预警等级以不同颜色标识(如红色代表一级高风险),同步触发对应级别的报警弹窗(见
图11)。
当预警触发后(见
图11),技术人员通过以下流程辅助决策:
(1)在预警界面(见
图10)记录异常设备位置信息(如“1号北斗毫米终端”)。
(2)返回“边坡稳定性分析”界面点选设备的位置,切换至三维可视化系统(见
图9)。
(3)调用该点位的毫米波雷达实时数据(见
图9(c)),复核边坡微变形状态。
(4)结合北斗位移时序曲线(见
图10)与InSAR历史形变(见
图9(d))评估风险趋势。
此流程通过多源数据空间关联提升研判效率(案例验证见4.2节)。
4 工程验证与案例分析
4.1 工程验证
为验证寒区高铁隧道口边坡多源融合智能监测预警系统在冻融循环阶段的适用性,选取利民隧道2021年4~6月3组典型监测数据进行分析。该时段覆盖春末冻融启动、夏季冻土消融及雨季高温耦合等特征工况,通过InSAR广域形变监测与北斗局部监测的协同验证,完整呈现系统“感知-分析-预警”流程的有效性。
图12(研究区4~6月InSAR形变分布)中红色五星集中区(高风险区)与边缘黄色区(潜在缓慢变形带)为局部监测点部署及预警等级判定提供宏观依据。
位置/时段:InSAR高风险区(127.172°E附近,3号监测区),2021年4月27日-5月4日(冻融循环启动阶段)。
监测触发:InSAR监测形变量达0.03~0.05 m(超二级预警阈值0.02 m);北斗水平位移速率达17~41 mm/d(超二级预警阈值1.5~2.5 mm/d),触发二级识别。
验证分析:毫米波雷达与高清影像未检出结构滑移或异物侵限,但捕获表层冻胀性微小隆起(幅值<5 mm),排除误触发。融合分析显示,地温为-0.5~0 ℃(冻融相变临界区),含水率为0.45%~0.54%,结合春末冻土未完全融化特征,判定冻胀力为主导变形机制(符合冻融-水热耦合效应理论)。
预警决策:综合判定为加速变形阶段,触发二级(橙色)预警。
位置/时段:InSAR潜在缓慢变形带(形变量为0.01~0.03 m,3号监测区外围),2021年5月14~21日(冻土消融后期)。
监测触发:InSAR形变量为0.01~0.03 m(超三级预警阈值0.01 m);北斗水平位移速率达15~28 mm/d(处三级预警阈值0.5~1.5 mm/d上限),触发三级识别。
验证分析:毫米波雷达与高清影像显示坡面无裂缝扩展,但检测到冻土消融性沉降(幅值3~8 mm),与InSAR形变趋势一致。融合分析显示,地温为3.58~4.95 ℃(>0 ℃,冻胀力消退),含水率为0.54%(受初雨及冻土残余水影响),虽未达饱和,但受“春融期水文激增效应”滞后作用导致持续变形。
预警决策:判定为等速变形阶段,触发三级(黄色)预警。
位置/时段:InSAR覆盖区外围(2号监测区),2021年6月1~8日(高温与雨季叠加阶段)。
监测触发:北斗监测显示2021年6月5日高程位移速率短时突增至-19.01 mm/d(超一级预警阈值2.5 mm/d),而6月8日累计水平位移达55.09 mm(突发性特征);两者分别代表不同日期的局部异常峰值。综合InSAR广域形变、地温含水率及雷达影像等多源信息,触发一级识别。
验证分析:环境参数显示地温为9.42 ℃(>5 ℃),含水率为15.4%。经多终端协同验证确认,位移突变呈局部孤立性(非区域性失稳)。融合研判表明,高温环境加剧土体弱化,叠加雨季潜在的饱和风险及“春融期水文激增效应”滞后影响,符合一级预警核心指标。
预警决策:触发一级(红色)预警。
4.2 冻融型滑坡前兆案例验证
2020~2021年监测期内,数字化预警平台成功捕捉一起冻融型滑坡前兆事件,验证了智能监测预警系统对高风险场景的精准识别能力。2号设备监测数据(见
图14)显示,2021年2月18~21日,水平位移由14.7 mm增至16.8 mm,高程位移由-25.4 mm降至-35.1 mm;同期土体温度从-5.47 ℃升至-5.42 ℃(冻土相变活跃区),含水率由0.32%升至0.35%。
关键预警特征如下:
(1)实时水平位移速率连续3日>2 mm/d。
(2)高程位移剧变(3日累计-9.7 mm)。
(3)土体温度>-5 ℃,且含水率达预警阈值(0.35%)。
系统据此于2月21日触发一级滑坡预警,响应流程如下:
(1)定位异常设备(2号北斗毫米终端,见
图10)。
(2)加载毫米波雷达数据(见
图9(c)),检测到该区域微变形≥1.2 mm。
(3)关联北斗位移加速趋势(见
图10)及InSAR沉降形变速率。
(4)综合判定为一级风险,启动应急响应。
现场验证:雷达监测点对应区域出现径向裂缝及土体松动,系统自动加密监测至1次/h。
系统响应时效(理论最坏值)如下:
(1)感知层:北斗阈值识别延迟≤1 s(采样率为1 Hz,见
表1)。
(2)传输层:4G网络延迟均值为3 s(见3.1节)。
(4)全链路响应时间≤8 s,远低于24 h人工巡检周期。
5 结 论
(1)融合北斗定位、InSAR遥感、毫米波雷达、温湿传感及图像识别技术,构建覆盖“结构变形-环境扰动-异常响应”的三维指标体系;提出“点-线-面-体”协同部署策略,实现位移场(精度±3 mm)、温度场(±0.5 ℃)、含水率(±1.5%)及裂缝演化的动态感知,突破寒区复杂环境时空监测瓶颈。
(2)建立基于位移速率核心指标、环境因子阈值与图像扰动判据的智能预警模型,形成“位移触发→图像复核→多源融合”的分级响应闭环;通过耦合位移统计阈值(莱因达准则)、温湿临界参数(见
表2)及裂缝扩展判据,为灾害孕育至临滑全过程分级预警提供有力支撑。
(3)研发数字化预警平台,集成多源数据处理、三维可视化与应急联动功能,实现监测-分析-响应闭环管理;利民隧道工程验证结果表明,本文提出的智能监测预警系统具备寒区冻融期高适应性(-30 ℃冰雪环境下设备可用率达100%)与可靠性(精准捕获2021年2月滑坡前兆,预警逻辑与形变特征吻合度>92%);多源模块协同响应时间≤8 s,显著优于24 h人工巡检时效。
(4)当前局限:①多源数据权重的动态优化需结合地质力学模型进一步深化;②单工程试点约束统计泛化能力;③极端天气下光学影像识别误差率>15%。未来方向:①开发地质特征驱动的自适应融合算法;②拓展多气候区(极寒/湿热)长周期验证;③构建环境鲁棒性增强的图像识别模型。
辽宁省教育厅科学研究项目(LJKZ0507)
辽宁省教育厅科学研究项目(LJKQZ20222462)