四川省雅江县上游河段滑坡灾害危险性评价

胡 鹏，杨东强，王荐霖

1.核工业西南勘察设计研究院有限公司，四川 成都610061；2.中国建筑西南勘察设计研究院有限公司，四川 成都610052

0 引言

中国西南地区地势具有西高东低的特点，特别是在一级地貌与二级地貌过渡地带，山高坡陡，滑坡灾害是其最常见、最严重的地灾类型［1-2］.尤其是2020年8月以来，四川省各地区连续性的强降雨使得滑坡灾害频发.因此，为了预测和防治滑坡灾害带来的损失，必须要弄清其诱发因素及成因类型，对该区域滑坡的诱发机率进行系统的分析，优先治理滑坡灾害发生机率大、造成损失严重的地区，才能更有效地减轻、避免各类地质灾害带来生命财产损失［3］.

本研究通过对雅江县上游两岸滑坡灾害分布发育特征的调查，对研究区内滑坡的各个要素及分布规律进行分析统计，提取有关因素，进而建立评价体系，预测评价区内滑坡发生的危险性，并划分危险等级，以期为雅江县流域内滑坡灾害的预警和防治提供依据，为当地的生产活动以及经济发展提供支持.

1 研究区概况

研究区位于四川省甘孜藏族自治州雅江县上游雅砻江流域，整体最大高差约2 500 m.研究区处于西南地区构造活动极为频繁的区域，区域构造属鲜水河断裂带与理塘断裂带之间［4］.两条断裂在第四纪以来活动频繁，加上印支期和喜马拉雅期以来造山运动产生的强烈的河流下切，形成了以中高山为主、呈典型的“V”型深切峡谷的地貌特征.区内灾害隐患点较多，通过野外实地调查共获得滑坡数据114个（图1、2）.

图1 研究区滑坡灾害点分布图Fig.1 Distribution map of landslide disaster points in the study area

2 滑坡灾害危险性评价

2.1 评价方法与原理

目前，区域滑坡危险性评价的方式有很多.相比于其他评价方式，逻辑回归法有着可解决指标间的相互依赖问题、模型清晰、概率值得起推敲等优点；而确定性系数法则可依托于野外实际情况，人为地对逻辑回归得出的结论进行修正［5-6］.故本文选用逻辑回归法与确定性系数法叠加分析，评价研究区滑坡灾害危险性特征.先在研究区内以一定的间距生成3 000个样本点，并随机选取调查滑坡中的100个，共3 100个样本数据作为危险性评价的原始数据；选取研究区内平均分布的14个滑坡数据作为评价成果的验证数据，用以验证评价结果的精确性.主要思路如图3所示.

图2 研究区地层岩性图Fig.2 Stratigraphy of the study area

图3 滑坡灾害危险性评价思路流程图Fig.3 Risk assessment flowchart of landslide disaster

逻辑回归模型的表达式为［7-9］：

式中，xi为关于Pi的自变量，εi为系数.

事件发生的条件概率Pi为：

则事件不发生的概率为：

从而：

式中，α0为常数项；β1…βn为各个评价指标的逻辑回归系数.

确定性系数模型的表达式为［10-13］：

式中，PPa为滑坡数目与评价指标面积的比值；PPs为滑坡总面积与研究区面积的比值.

单个评价因子的确定性权计算公式为：

式中，Wci为评价因子Ci的权重；CFij为评价因子Xi第j个分类的确定性系数.

综合确定性权的计算公式为：

式中，LSP为评价指标的确定性权.

2.2 评价指标体系

依据控制研究区滑坡发育特征的工程环境条件，可将影响滑坡危险性评价的因素分为7个内在指标及1个外在指标（图4）.

图4 滑坡灾害危险性评价指标体系Fig.4 Risk assessment index system of landslide disaster

2.3 滑坡灾害分布特征

2.3.1 坡度

坡度主要影响地面的起伏程度，直接决定斜坡体临空面的发育情况及应力分布特征.通过分析研究区滑坡与地形坡度关系不难看出，坡度为30～45°区间滑坡数量最多，为54个，占滑坡总数的47.37%（图5a、6a）；此坡度区间灾害点密度也最大，为7.06个/100 km2.滑坡的诱发是一个循序渐进的过程，斜坡覆盖层的天然休止角通常为35°左右［14］.当地形坡度过小时，无法诱发滑坡；当地形坡度堆积到天然休止角（30～45°）时，滑坡极易发生了，而过大的地形坡度造成的地质灾害则以崩塌为主.只有当地形坡度在一个合适的区间内时，才容易诱发滑坡.

2.3.2 坡向

斜坡坡向对于斜坡稳定性的影响主要取决于太阳辐射对斜坡树木发育、土体冻融、土壤湿度等要素的影响，进而间接影响滑坡的发生［15-17］.研究区朝向北侧斜坡土体水分不易蒸发，土壤湿度大于南坡，土体的冻胀情况亦强于南坡，从而使得北坡土体裂隙数量多于南坡，故发生滑坡的数量较多（21个，占滑坡总和18.42%），灾害点密度较大，为9.49个/100 km2，（图5b、图7）.

图7 滑坡数量与坡向关系图Fig.7 The diagram of landslide numbers by slope aspects

2.3.3 拔河高度

斜坡的拔河高度反映了河流对河谷底部的切割程度，反映了斜坡临空面的发育情况及坡体的应力分布特征［18］.当拔河高度过低时，临空高度不足，但其受到人类工程活动的影响较大；当拔河高度过高时，人类工程活动较少，坡顶地形坡度较大，崩塌为主要的地质灾害方式.研究区拔河高度600～900 m范围的滑坡数目最多，为36个，占总和的36.98%；拔河高度300～600 m范围灾害点密度最大，为8.07个/100 km2，占总数30.70%（图5c）.滑坡的数量与拔河高度呈现近似正态分布的关系（图6b）.

2.3.4 工程地质岩组

依据研究区地层岩性、岩体结构面发育情况、岩体风化情况、岩体强度的不同，可将区内地层分为较软岩类、较软岩夹较硬岩类、较硬岩夹较软岩类、较硬岩类4类.通过图5d可知，研究区滑坡主要分布于软硬互层岩类的边坡中.当斜坡岩层性质一致时，岩层的内摩擦角、内聚力与斜坡稳定性系数呈近似线性关系；而斜坡岩层性质不一致时，这种线性关系便发生改变，雨水极易沿岩体结构面入渗坡体，软岩受雨水软化后便会形成软弱带，斜坡极易沿着这层软弱带产生滑动.当软岩厚度较小时，岩层内聚力及内摩擦角的细微变化都会对斜坡稳定系数产生较大影响.因而，较硬岩夹较软岩地区相对于其他区域更容易引发滑坡［19］（图6c）.研究区较硬岩夹较软岩类地区滑坡个数72个，占滑坡总和的63.16%，灾害点密度为0.10个/100 km2.

2.3.5 断裂带

研究区范围内断裂带数量虽较少，但通过查阅区域地质资料可知，距雅江县城下游约10 km处，断裂带较发育，受区域断裂影响较大.研究区断裂带附近滑坡分布相对密集（图5e），滑坡与断裂带距离在1000～2000 m范围灾害点密度最大，为18.09个/100 km2.

2.3.6 岸坡结构

虽然研究区岸坡结构以横向坡、斜向坡为主，其中斜向坡分布滑坡数目最多，但顺向坡结构松散，抗风化能力低，受扰动后极易引发滑坡，灾害点密度最大，为9.48个/100 km2，占滑坡总和的22.81%（图5f、6d）.

2.3.7 河流水系

研究区属深切峡谷地貌，河流的侧蚀作用及下蚀作用不断冲刷岸坡，致使岸坡岩体软化，临空面发育，斜坡稳定性降低，容易造成新滑坡的发生及老滑坡的复活.研究区内滑坡沿雅砻江干流、支流两岸呈明显的沿河流线性分布特征，其中距河流水系0～200 m范围滑坡数目最多，为35个，占35.07%；距河流水系200～400 m范围灾害点密度最大，为10.04个/100 km2，占滑坡总和的29.82%（图5g、6e）.

2.3.8 降雨

雨水入渗，增大了坡体的下滑力，同时研究区地层以较硬岩夹较软岩类为主，所夹较软岩极易受到雨水的侵蚀，导致承载能力降低，形成滑动面.研究区年平均降雨量主要为780 mm（图5h），该范围内滑坡的数目最多，为90个，占滑坡总和的78.95%，灾害点密度最大，为8.47个/100 km2（图6f）.

图5 研究区滑坡危险性评价图Fig.5 Risk assessment maps of landslide disaster in the study area

图6 评价指标与滑坡数量关系Fig.6 Relationship between evaluation indexes and amount of landslide

2.4 滑坡灾害结构组成特征

经野外实地调查，研究区滑坡可分为土质滑坡及岩质滑坡两大类.

2.4.1 土质滑坡这类滑坡多发生于沿道路两侧及已有滑坡前缘产生的二次滑动，滑坡规模通常相对较小，滑源区高程较低，受人类活动的影响较大，危害形式以堵塞道路为主.滑坡堆积体主要为残坡积、滑坡堆积块碎石土，碎石粒径通常在3～5 cm.

2.4.2 岩质滑坡

这类滑坡是研究区滑坡的主要类型，滑坡的稳定性受“锁固段”控制，滑坡形态以长条形、箕形为主，规模通常较大，危害也较大.依据滑坡堆积体组成结构特征，可将其分成完全解体滑坡和不完全解体滑坡两类.

完全解体滑坡多发生在拔河高度较高地带，滑坡滑动区域长度通常较长，整体落差较大.滑体处在滑源区时有着较大的势能，当“锁固段”破裂滑动产生时，这些势能全部转化为动能并顺着结构面释放，使滑体在滑动时完全解体，呈高速流动状态.堆积体形态为块石堆积，滑坡形态呈长条形、扇形（图8）.

图8 完全解体滑坡堆积体形态特征Fig.8 Morphological characteristics of completely disintegrated landslide

不完全解体滑坡堆积体形态多为箕形，滑体主要沿岩体结构面以整体下挫的形式或沿基覆界线产生滑动，后缘可见圈椅状地形、双沟同源等滑坡的典型特征.这类滑坡堆积体前缘多有“假基岩”出露，表面往往有一定厚度的覆盖层，由残坡积、崩坡积和滑坡堆积组成（图9）.

图9 不完全解体滑坡堆积体形态特征Fig.9 Morphological characteristics of incompletely disintegrated landslide

2.5 危险性评价结果

通过2.3节的分析可知，滑坡的诱发是由多个因素引起的，且每个因素之间都有着内在的联系.逻辑回归分析的核心在于通过对比各个因素之间的内在联系从而得出回归系数，可使评价结果的准确性得到保障.其具体公式如下［20-21］：

式中，Rij为滑坡灾害面积Szij与其对应的评价指标面积Sij的比值；Xij为评价指标的归一化值；ij为第i个评价指标的第j个子集；m为子集个数；

依据公式（9）（10）可求得研究区各个评价指标分级的归一化值（表1）.如表1中坡度0～15°范围内滑坡面积与坡度为0～15°区域面积的比值即为Rij，Rij与研究区面积的比值即为Xij.

表1 评价指标归一化值计算表Table 1 Calculation results for normalized value of evaluation indexes

随后将各评价指标的归一化值Xij赋值在3 100个样本点上，再将其导入SPSS软件计算各个评价指标的回归系数.回归系数值越大，说明其对滑坡诱发的影响程度越高.计算结果见表2.

依据公式（6）（7）计算各个评价指标分级的确定性系数权值（表3）.

将所得回归系数（表2）与各评价指标分级的确定性系数权值（表3）分别代入ArcGIS软件中，加权求和分别得出研究区滑坡发生的概率P1、P2，随后将P1与P2相叠加，进而得到研究区滑坡灾害危险性成果图（图10）.

图10 滑坡灾害危险性成果图Fig.10 Risk zoning map of landslide disaster in the study area

表2 回归系数计算结果Table 2 Calculation results of regression coefficient

表3 各个评价指标的确定性系数权值Table 3 Certainty factor weights of evaluation indexes

3 评价结果验证与分析

为了确保评价结果的准确性，利用以下两种方法对其进行检验.

1）查验点验证

研究区共调查滑坡灾害114个，随机选取其中的100个滑坡点进行评价分析，剩余14个点作为查验点，用以验证成果的精确性.14个查验点在危险性分区图中的分布如表4所示.

表4 查验点分布统计表Table 4 Risk zoning and number of landslide disaster inspection points

2）ROC曲线验证

本研究选用ROC检验曲线用以验证评价结果的精确性.ROC曲线下面积越大，表明评价结果的准确程度越高.ROC曲线下面积为85.48%（图11），评价结果精确性较高.

图11 ROC检验曲线Fig.11 ROCtest curve of risk zoning

评价结果表明：研究区高危险区、极高危险区主要位于呷拉镇一带山脊的两侧，分布滑坡74个，占总数的65.91%，灾害点密度为8.00个/100 km2.该区域范围内年平均降雨量为780 mm，地层以较硬岩夹较软岩为主，岸坡结构多为顺向坡，且分布有一条断裂带，皆为易引发滑坡的致灾因素.评价结果（表5）与实际野外调查结果相符.

表5 滑坡灾害危险性评价结果统计表Table 5 Risk evaluation results of landslide disaster

4 结论

依据雅江县上游河段滑坡的野外调查结果，结合DEM高程数据，分析研究区内滑坡灾害的发育分布特点，并基于逻辑回归与确定性系数的叠加分析，对研究区滑坡危险性进行分区评价，取得如下成果.

1）雅江县上游河段滑坡主要分布在海拔2 500～3 000 m、拔河高度600～900 m、坡度30～45°、距离河流水系0～200 m范围内较硬岩夹较软岩类地带.

2）雅江县上游河段滑坡以岩质滑坡为主，岩质滑坡又可进而分为完全解体滑坡和不完全解体滑坡两类.完全解体滑坡多发生在拔河高度相对较高、岩体结构破碎地段，滑体多以碎屑流的方式滑动，滑坡形态多呈长条形，堆积体为块石堆积；不完全解体滑坡形态多呈箕形，滑体沿岩体结构面整体下挫或沿基覆界线产生滑动.

3）雅江县上游河段高危险区、极高危险区面积占总面积的46.75%，该区内滑坡数目为滑坡总和的65.91%.说明该区域内滑坡分布相对聚集，危害程度较高，与实际野外调查结果相符.

4）雅江县上游河段极高危险区、高危险区主要位于呷拉镇一带山脊的两侧；中等危险区主要位于呷拉镇、瓦多乡斜坡的中部，普巴绒乡上游河段；低危险区、极低危险区主要位于普巴绒乡、瓦多乡坡脚河流的两岸，木绒乡上游河段.

5）经查验点及ROC曲线验证，危险性评价成果有较高的精确度（85.48%），能够作为研究区防灾减灾工作及河谷两岸开发利用的合理方案依据.