降雨作用下滑坡灾害监测预警研究*——以遂宁市船山区为例

孙 瑜，韦 玮，王宝亮，李宏俊

(1.四川省核工业地质调查院，四川 成都610061；2.核工业西南勘察设计研究院有限公司，四川 成都610000； 3.云南国土资源职业学院，云南 昆明652501)

四川地跨中国Ⅰ、Ⅱ级地貌台阶，地质环境较为脆弱，属于全国地质灾害发育最严重的省份之一，近年极端天气成频发、多发、高发态势，地质灾害防治工作形势严峻。降雨是导致地质灾害失稳的“元凶”之一[1-2]，尤其对滑坡、崩塌之类的斜坡地质灾害的形成起决定性作用。

国内外众多学者开展研究滑坡和降雨量的关系，已提出不同的滑坡灾害预警预报方法[3-7]。研究发现，滑坡失稳破坏在过程降雨量和降雨强度两个参数中，均存在临界值[8]，当某次降雨发生时，其累计降雨量或降雨强度达到或超过其临界值，极可能导致滑坡灾害成群出现。本文选择四川丘陵区遂宁市船山区为研究区，应用气象-地质耦合模型进行地质灾害监测预警分析。

船山区地处涪江中游，位于四川盆地腹心，境内沟谷迂回，丘包群立，宽缓展布，属于四川盆地典型的因构造剥蚀而成的丘陵景观，涪江两岸大量分布河流堆积阶地。研究区位于新华夏系第三沉降带，地处四川沉降褶皱带之川中褶皱带内，其构造特点表现为以褶皱为主，构造形迹主要为近东西向，褶皱宽阔平缓，区域上表现为排列有序的箕状向斜和鼻状背斜。境内主要出露中生界侏罗系中统上沙溪庙组(J2s)、上统遂宁组(J3s)和蓬莱镇组(J3p)，以及第四系全新统(Q4)地层，基岩主要为近水平状的红色泥岩、砂岩。全区处于四川盆地亚热带湿润季风气候区，四季分明，雨量充沛，气候温和，季风气候显著，具有冬暖春早、夏热秋凉、云雾多、日照少、无霜期长等特点。船山区多年平均气温为17.6℃，每年降雨主要发生在4—10月，占全年降水量的74.6%。

1 降雨与滑坡的关系

根据对2001—2020年船山区降雨资料分析，区内多年平均降雨量约982.8 mm，常年降雨量549.8～1 311.0 mm，最大日降水量50.2～323.7 mm。通过分析该区近年开展的多次地质灾害调查、排查及核查资料，收集当地各行业管理的相关地质灾害统计数据，得到区内滑坡灾害隐患点历年变化情况(表1)。

通过分析历年滑坡失稳破坏时对应的当日降雨量和月累计降雨量数据，可以编绘对应的失稳滑坡数量与日最大降雨量关系图(图1)、失稳滑坡数量与月累计降雨量关系图(图2)。

表1 船山区2001—2020年的降水量和滑坡统计表

图1 日最大降雨量与失稳滑坡数量关系图

图2 月累计降雨量与失稳滑坡发生数量关系图

1.1 滑坡频度与最大单日降雨量的关系

从图1可以看出，船山区日最大降雨量在30～80 mm之间所引起的滑坡数量较多，在日最大降雨量超过190 mm后滑坡失稳数量达到峰值。在一般年份下，滑坡多发生在日降雨量位于20～50 mm之间。位于图1上方的两个点分别对应2013年6月30日和2013年7月1日的两次较大的降雨事件，相应的单日降雨量分别为323.7 mm和197.1 mm，其引发的滑坡数量最多。

1.2 滑坡频度与月累计降雨量的关系

从图2可以看出，船山区滑坡多发生在月累计降雨量超过100 mm时，而累计降雨量达到330 mm时，失稳滑坡数量达到峰值。

参考陈伟[8]的框架：降雨诱发滑坡与前期降雨以及当日降雨关系密切。本文选用滑坡失稳破坏前3 d期降水量(R3)和当日降水量(R0)作为预警预报参数，由于船山区内各雨量站建站时间不一，控制范围有限，并不能完全反映滑坡失稳处局地实际降水量，基于偏安全考虑[8]，将图3中左下角位于滑坡点总数10%的边界作为滑坡发生概率为10%的下包络线，将此下包络线做平行线移至上方包括滑坡点总数90%的位置，此线作为上包络线，根据上、下包络线中间区域滑坡数据分布的离散特征，以包罗尽量多的高频降雨量激发的滑坡点数据为原则，确定滑坡发生概率Pro=50%对应界线[9]，相应数据见表2。

图3 船山区滑坡降雨预警模型

表2 滑坡不同发生概率对应R0、R3值

2 滑坡灾害监测预警设置建议

2.1 不同设计频率最大暴雨量

船山区年最大1/6 h、1 h、6 h、24 h点雨量均值和变差系数主要来源《四川省暴雨统计参数图集(川水发[2010]15号)》，结合Cs=3.5Cv的皮尔逊Ⅲ型曲线的分析，获取不同频率对应的模比系数，从而计算得出不同设计频率对应最大暴雨量(表3、表4)。

表3 不同时间段设计降雨量基础参数

2.2 滑坡预警

滑坡预警应结合滑坡隐患点的基本地质特征、形成机制及发展变化趋势，在全面统筹考虑的基础上，合理确定有效预测参数[10-12]。通常情况下，浅表层位移监测数据是滑坡隐患预测的基本参数，降雨量、地下水动态等监测数据属于滑坡隐患预测的表征参数。同时，滑坡预警亦需全面考虑宏观前兆变形现象。

(1)预警模型。滑坡预警模型建立和预警判据确定，建议参照以下方法和原则：

①合理建立适宜的滑坡地质模型，辅以科学的监测方法和监测内容，确定有针对性的滑坡预警模型。

②积极建立工程地质类比分析、决策树、层次分析、逻辑回归、人工神经网络、支持向量机和随机森林等数学模型，进行多因子和多模型的综合评判，保证滑坡预警的准确度。

③滑坡预警模型建立后，积极利用已获取历史上的滑坡失稳监测资料，进行反演分析，并检验模型的适宜性和有效性，最后确定相应的滑坡灾害预警判据。

④通过对宏观地表变形迹象、自动化监测设备获取的数据分析，以及区域地质灾害气象预报的综合研判，参考地表位移、位移变形速率和位移加速度等监测结果，从而确定滑坡预警模型，最后确定四级预警的相关指标。

⑤通过不同类型预警指标获取相应参数来确定预警等级存在差异性，预警结果应结合上述参数综合判定。

滑坡地质灾害预警模型可见图4。

(2)预警方案。监测阶段主要结合滑坡诱发因素和变形失稳的运动特点，采用地表裂缝位移监测、降雨量监测[13-16]。按坡体变形破坏的变形特点、发展阶段及可能发生的时间，可将滑坡预警监测概括为四类[17]：Ⅰ级(警报)、Ⅱ级(警戒)、Ⅲ级(警示)、Ⅳ级(注意)。基于斜坡变形演化三阶段理论进行预警分级(图5)。

表4 不同设计频率最大暴雨量计算表

图4 滑坡地质灾害预警模型

图5 滑坡变形的三阶段演化图

①Ⅰ级(警报)：滑坡处于临界变形状态，如图5所示DE段，短期临滑前变形征兆特征极为显著，在数小时或数日内发生失稳的概率极大。由警报降雨量、变形位移、变形速度、变形加速度等指标确定。降雨量达到警报雨量阈值，同时参考位移参数变化情况确定。采取措施：组织应急抢险队伍，第一时间撤离危险区和影响区的所有人员，并对相关区域进行交通管制。

②Ⅱ级(警戒)：滑坡处于初始加速变形阶段，如图5所示CD段，局部地段出现滑坡宏观变形前兆特征，在数日内或数周内发生破坏的概率较大。由警报降雨量、变形位移、变形速度、变形加速度等指标确定。降雨量达到相应的警报雨量阈值，并结合地表位移监测指标变化情况综合确定。采取措施：增大所有监测项目的专业监测频率，启动相应的应急预案，划定滑坡灾害的危险区和影响区，组织危险区内的所有人员迅速撤离避让，并对滑坡危险区进行交通管制。

③Ⅲ级(警示)：滑坡处于匀速变形阶段，如图5所示BC段，坡体出现明显变形，在数月或1年内发生破坏的概率很大。由警报降雨量、变形位移、变形速度、变形加速度等指标确定。降雨量达到相应的警报雨量阈值，并结合地表位移监测指标变化情况综合确定。采取措施：组织专业监测人员进行专业加密监测。

④Ⅳ级(注意)：滑坡处于初始变形阶段，如图5所示AB段，1年内发生破坏的可能性不大。由警报降雨量、变形位移、变形速度、变形加速度等指标确定。降雨量达到相应的警报雨量阈值，并结合地表位移监测指标变化情况综合确定。采取措施：加强群测群防工作，并制定应急预案。

根据滑坡变形各类相关监测数据，进行综合统计分析，当监测参数达到某一预警等级值或情况时，按“就高不就低”的原则，赓即启动相应的预警预案。根据船山区实际情况，降雨量作为首要监测指标，当小时降雨量、日降雨量、累计降雨量达到预警值时，即可确定相应预警等级；此外，地表裂缝位移监测作为次要监测指标，同时考虑雨水等的影响；在上述基础上，结合滑坡隐患点群测群防监测员获取的实时情况相应启动适宜的预警等级。综合上述情况，当地政府应及时组织专家、政府管理人员进行会商，启动相应的滑坡应急预案，重点是优先级应急预案。

(3)预警阈值。根据表3至表5获取的降雨相关数据，结合船山区地质灾害风险动态管控机制，提出相应的滑坡灾害隐患点监测预警阈值，按照“一点一策”的原则调整单点的降雨阈值。滑坡预警等级同时结合“船山区近年开展的地质灾害自动化实时专业监测体系建设项目”监测预警阈值参数确定(表5)。

3 结论和讨论

(1)遂宁市船山区内发育的滑坡失稳与降雨关系密切，尤其是滑坡失稳破坏前3 d期降水量(R3)和当日降水量(R0)，可以作为滑坡失稳的主要判识参数。

(2)基于20年历史降雨数据分析，初步建立滑坡降雨相关预警模型，获取滑坡在不同发生概率时对应的R0、R3值，结合不同设计频率最大暴雨量计算结果，综合考虑从地质成因机制分析得出的滑坡不同变形阶段的破坏特征，提出船山区滑坡监测预警模型、预警方案和预警阈值，有助于指导当地滑坡灾害防御和主动避灾防灾。

(3)鉴于滑坡失稳受多因素控制，基于气象降雨模型的监测预警未全面纳入全部滑坡影响因子，上述监测预警方案在接下来的地质灾害防治实践中，仍需结合每年实际新增的滑坡灾害隐患情况，进行相应修正，提高预警精度，才能达到有效服务于船山区滑坡灾害预警的目的。

表5 滑坡监测预警等级表划分