Regional landslide hazard assessment using the IV-RF coupling model and critical monthly average rainfall threshold:A case study from Fuling District, Chongqing
-
摘要:
提高降雨型滑坡易发性预测精度和构建合适的降雨阈值模型对区域滑坡危险性评价具有重要意义。以重庆市涪陵区为例,采用信息量模型、BP神经网络模型、随机森林模型、信息量-BP神经网络耦合模型和信息量-随机森林耦合模型进行区域滑坡易发性评价,对比不同模型下的接受者操作特征曲线、曲线下方面积和易发性分布规律。提出滑坡临界月平均降雨阈值模型,反演出不同时间概率下的临界月平均降雨阈值。将易发性结果与时间概率等级进行耦合得到区域滑坡危险性评价结果并随机选取30次滑坡事件与4次典型滑坡事件进一步验证了评价精度。研究结果表明:信息量和机器学习模型进行耦合,弥补了机器学习在前期数据输入和非样本选择的缺点,提升了单一机器学习模型的预测精度,其中信息量-随机森林耦合模型预测精度最高;随机选取的30例滑坡样本中,有20例滑坡(占67%)位于发生时间概率50%以上区域,验证了临界月平均降雨阈值模型的精度;随机选取的4例典型滑坡样本中,时间概率等级基本为P4或P5,且位置均位于高危险区与极高危险区中,与现场调查结果基本一致,说明基于信息量-随机森林耦合模型和临界月平均降雨阈值的区域滑坡危险性评价结果准确且可靠。
Abstract:Improving the accuracy of susceptibility prediction for rainfall-induced landslides and establishing suitable rainfall threshold models are of great significance for regional landslide hazard assessment. Taking Fuling District of Chongqing as a case study, the information value model, BP neural network model, random forest model, information value-BP neural network coupled model, and information value-random forest coupled model were used to evaluate regional landslide susceptibility. By comparing the receiver operating characteristic (ROC) curves, area under the curve (AUC), and susceptibility distribution patterns of different models, a critical monthly average rainfall threshold model for landslides is proposed, and critical monthly average rainfall thresholds for different temporal probabilities were inferred. The susceptibility results were coupled with temporal probability levels to produce regional landslide hazard assessment results. The evaluation accuracy is further validated with 30 randomly selected landslide events and 4 typical landslide cases. The results show that coupling the Information Value and machine learning models compensates for the shortcomings of machine learning in early data input and non-sample selection, enhancing the predictive accuracy of single machine learning models. Among these, the information value-random forest coupled model exhibits the highest predictive accuracy; of the 30 randomly selected landslide samples, 20 cases (67%) occurred in areas with a temporal probability of over 50%, validating the accuracy of the critical monthly average rainfall threshold model. The 4 typical landslide samples selected randomly were primarily in the P4 or P5 temporal probability levels and were located in high to very high-risk areas, aligning well with field survey results. This indicates that the regional landslide hazard assessment based on the information value-random forest coupled model and the critical monthly average rainfall threshold is accurate and reliable.
-
0. 引言
随着全球气候变暖和极端天气事件的增多,暴雨引发的地质灾害在全球范围内发生频率显著提升,并造成了严重后果[1 − 4]。尤其是在山地和丘陵地区,降雨诱发的浅层滑坡不仅频率增加,其破坏性也日益显著[5 − 8]。这类滑坡通常具有较高流动性,在短时间内沿坡面快速滑移扩展,造成大面积地表破坏,对区域环境和人类活动构成了重大威胁。浅层滑坡不仅常导致道路阻断和基础设施受损,甚至直接危及人类生命安全[9 − 12]。目前,研究者主要围绕滑坡的编目建立[13 − 14]、空间分布[15]、滑坡评价[16 − 19]、诱发机制[20 − 24]和监测预警[25 − 30]等方向开展了广泛研究。其中,详细且精确的滑坡数据编目是所有研究的重要基础。滑坡数据编目的主要方法包括基于遥感影像的目视解译法、自动化提取法以及基于地面调查的人工标注法等。目视解译法具有精度高的优势,但耗时较长,且对解译者经验依赖较大[31 − 33];自动化提取法则能够处理大范围区域的滑坡识别,但在复杂地形条件下容易受到误差干扰[34 − 36];地面调查法虽然精确,但难以覆盖大范围区域,因此在时效性上存在较大局限性。
2024年4月发生在中国广东省韶关市的极端暴雨事件在短时间内引发了大量浅层滑坡,并造成区域内广泛的基础设施破坏[37]。尽管已有研究关注到了这次事件并初步对滑坡分布进行了解译。但是受限于时效性,滑坡编目的完整性和精确度有待进一步提高。因此,本研究利用降雨前后的高分辨率卫星影像,结合GIS平台进行目视解译,同时辅以实地调查数据进行验证,最终构建了此次事件诱发的详细滑坡编目库。考虑到降雨型滑坡的失稳过程不仅由单一因素决定,而是多种因素综合作用的结果,本研究在分析滑坡分布时,纳入了地形、地质构造等九类环境控制因素。此外,还通过数理统计分析探索了各个控制因素与滑坡分布的关系。研究结果不仅为强降雨诱发滑坡的易发性研究提供了重要参考,也为区域内滑坡防治措施的制定提供了科学依据。
1. 研究区及事件概况
2024年4月中旬以来,广东省大部分地区遭受持续性暴雨侵袭,特别是在4月16日—22日期间,暴雨和特大暴雨覆盖了包括韶关在内的多个区域,极端降水过程给当地造成了严重影响。据广东省气象服务中心发布,4月1日—30日,广东全省平均雨量497.4 mm,打破4月雨量历史纪录,较常年同期(176.8 mm)显著偏多181%。其中韶关市4月平均雨量是常年3.96倍,为增幅最大(https://weibo.com/2015316631/OcokI8rmB)。图1展示了4月以来韶关市区气象站点数据,其中多日降雨量超过100 mm,最大单日降雨量接近200 mm,累计降雨量更是接近
1000 mm。广东省韶关市武江区江湾镇位于粤北地区,地处华南褶皱系的南岭构造带内,区域地形以山地和丘陵为主(图2),整体呈现四周高、中间低的地貌特征。区域内河网发育,水系分布广泛。地质构造复杂,基岩主要由花岗岩和片麻岩组成,上覆地层则以第四纪松散堆积物和风化残积土为主[38 − 39]。受长期构造活动和强降水影响,该地区成为地质灾害高发区,特别是在暴雨等极端天气条件下,滑坡和泥石流等灾害频率显著增加[40]。在2024年4月的极端暴雨事件中,江湾镇作为粤北暴雨中心之一,降雨量远超常年均值,极端降水引发了大量山体滑坡灾害,导致多条道路中断,房屋倒塌,基础设施严重受损,给区域经济和居民生活带来了巨大影响。
2. 数据与方法
已有研究表明,降雨诱发的浅层滑坡在遥感影像上通常表现为显著的色调变化和纹理特征差异。本研究利用高分辨率光学遥感卫星影像,采用人工目视解译的方法提取滑坡数据。具体而言,基于Planet卫星3月和4月的月度合成影像(空间分辨率为3 m),通过对比降雨前后多时相影像,能够精准识别此次暴雨触发的滑坡位置与边界,同时有效剔除非此次降雨事件引发的滑坡,从而确保数据的高精度和时效性。
如图3所示,图3(a)(b)展示了同一区域在降雨前后的遥感影像,滑坡区域在降雨后显现出明显的地表裸露特征,其边界以黄色线条标注。类似地,图3(c)(d)为另一处区域的影像对比,也清晰反映了降雨导致的滑坡发生范围。为进一步验证遥感解译结果,本研究结合了灾害发生后的现场调查,如图3(e)(f)所示。图3(e)拍摄于现场滑坡点,图3(f)记录了另一处滑坡现场,照片显示了滑坡坡面破坏和周边环境的具体情况。通过实地观测的滑坡特征,对遥感影像上降雨滑坡的解译标准进行调整优化,尽可能保证数据的可靠性和准确性。
影响因子的选择方面,我们考虑了相关的研究成果和经验,包括已有文献中对降雨触发滑坡分布规律和易发性的分析[41 − 45],同时也参考了研究区域的具体情况[41, 46 − 48]。对所选因子进行了相关性分析,参考Pearson相关系数和统计意义对高度相关且冗余的因子进行剔除。例如研究区范围较小,所涉及的岩性较为单一(花岗岩残积土为主),就本研究范围而言,其统计价值相对有限,因此后续未进一步分析。值得一提的是,尽管岩性在统计分析层面暂未凸显关键作用,但鉴于具有较强代表性,因此该区域是开展物理模型试验的理想研究区。最终我们选取高程、坡度、坡向、剖面曲率、相对坡位指数、地形湿度指数、土地利用类型、距离道路距离和距离河流距离纳入分析范围。
对于影响因子的计算,本研究采用了30 m分辨率的数字高程模型(digital elevation model,DEM)数据,其具体来源为Copernicus DEM(https://doi.org/10.5270/ESA-c5d3d65),直接反映了高程信息。基于GIS软件的表面分析功能计算DEM中的坡度、坡向和剖面曲率。基于地形分析功能,计算获得相对坡位指数和地形湿度指数。土地利用类型数据来源于中国科学院发布的30 m分辨率GLC_FCS30D数据(https://zenodo.org/records/8239305),本文对其进行了分类和重分类,以方便后续分析。水系和道路数据则基于自然资源部发布的全国1∶100万基础地理信息(https://www.webmap.cn/)。本文对水系和道路进行距离分析,确定各研究区各点位置到水体及道路的实际距离,得到距离数据。上述因子最终均以30 m分辨率的栅格TIFF 格式导出并进行后续分析。具体研究区的因子分布如图4所示,其中图4(g)的土地利用类型在表1中具体列出,提供了各影响因子的地理位置和分布情况。
表 1 土地利用类型编号Table 1. Land use type codes编号 具体类型 10 雨养耕地 11 草本植被覆盖 20 灌溉耕地 51 开阔常绿阔叶林 52 封闭常绿阔叶林 61 开阔落叶阔叶林 62 封闭落叶阔叶林 71 开阔常绿针叶林 72 封闭常绿针叶林 120 灌木地 121 常青灌木地 180 湿地 190 不透水面 210 水体 尽管部分因子分辨率与遥感影像分辨率(3 m)存在差异,但二者的适用目的并不相同:高分辨率影像用于精准解译滑坡边界与位置,确定滑坡区域。而因子数据旨在提取区域尺度的控制因素(如坡度、地形湿度指数等),揭示滑坡分布与地形因子的普遍规律。因子分辨率的精度可能会在一定程度上导致结果的偏差。但受限于数据获取限制与难度,文中所采用的数据已是目前研究区公开可获取的高分辨率的因子数据。此外,本研究关注的是区域尺度滑坡分布与地形因子的统计规律,而非单点精度,因此分析结果仍具有代表性和可靠性。
3. 结果与分析
3.1 滑坡数据库及空间分布
解译结果表明,此次降雨事件在研究区范围内共计触发
6310 处浅层滑坡。图5展示了2024年4月韶关暴雨事件后,在研究区内解译的滑坡分布与密度情况。滑坡密度图以栅格单元为统计单元,1 km为搜索半径,计算得到滑坡的数量密度。直观呈现了滑坡的空间分布特征和密集区域。从滑坡数量地理分布来看,呈现显著的空间聚集性。在图5标注的1、2、3、4四个区域,滑坡密度明显高于周围区域,达到了每平方千米分布有近150个滑坡。这些高密度区域整体呈现北东—南西的趋势。图5(b)(c)子图为图5(a)中标注的1和2区域的放大图,图3(b)(d)为图5(a)中标注的3和4区域的放大图,展示了滑坡在遥感影像上的具体形态特征。结合实地调查发现,上述高密度区域均位于省道S520(江湾段)两侧。
此次降雨事件触发的滑坡的总面积达到5.85 km2,其中最大面积为
22368.6 m2。图6展示了滑坡面积-滑坡密度图。并用对数正态分布函数进行拟合。随着滑坡面积的增加,滑坡密度先增加然后逐渐减少。这意味着小型规模的滑坡更为常见,集中在102~103 m2,而大型滑坡则相对较少。具体的统计结果表明,有3853 处滑坡面积在102~103 m2,占总滑坡数量的66.7%,这些滑坡的总面积达到1.72 km2,占滑坡总面积的29.9%。3.2 几何形态特征分析
在滑坡流动性分析中,滑坡的几何形态参数是反映其运动特征的关键指标[49 − 51]。根据国际工程地质与环境协会(IAEG)滑坡术语专委会(C37-Landslide Nomenclature)工作组的研究成果,选择开源测量工具ALPA对滑坡几何形态参数进行提取[52]。图7(a)展示了降雨滑坡面积与滑行距离之间的关系。通过频率直方图显示,最大滑行距离为311.62 m,大部分滑坡的滑行距离在75 m以下,占比80.84%,共有
4672 处。其中,滑行距离在25~50 m之间的滑坡有2137 处,占比37%。散点图的横坐标为滑坡面积(A)纵坐标(L),数据点通过颜色变化表示密度。采用幂律拟合模型得到的曲线(式1),其R2值为0.78,清晰地展示了滑坡面积与滑行距离之间的正相关趋势,并提供95%置信区间和预测区间。直方图的纵坐标与散点图保持一致,进一步支持所观察到的趋势。$$ L=\mathrm{lg}A^{3.72} $$ (1) 式中:L——滑行距离/m;
A——滑坡面积/m2。
图7b 显示了滑坡面积(A)与滑坡长宽比(ε)之间的关系。直方图表明,ε值最高可达16.49,主要集中在1~4,占比67.23%,共有3 885处滑坡。尽管两者之间存在一定相关性,但由于数据复杂性,无法用单一数学模型准确拟合。因此,研究采用包络线(式2)描述总体趋势,绝大多数滑坡散点位于包络线以下,高密度区域集中在面积102~103 m2,ε值在1~4,表明这类滑坡更为常见。随着面积增大,长宽比通常增高,这反映出较大滑坡往往形状更为延展;同时,数据的高度分散性提示,不同滑坡之间的几何变异性源自多种因素,如地形特征、地质构造和环境条件。
$$ \varepsilon=4.05\mathrm{lg}A-4.97 $$ (2) 式中:ε——滑坡长宽比;
A——滑坡面积/m2。
图7(c)展示了滑坡(L)与高差(H)之间的关系。大部分滑坡的高差为20 m以下,包括
3183 处滑坡,占比55.08%,最大滑坡高程达到132.49 m。采取线性模型得到的拟合线表明滑行距离与高差呈显著正相关(式3),拟合模型的R2为0.63。这表明随着高差的增加,滑坡的滑行距离也趋于增加。这意味着,较高的滑坡起始点往往导致更远的滑行距离。这种关系与滑坡的动力学过程相关,其中高差提供了更多的势能,从而增加了滑坡体的动能和流动性。$$ H=0.41L+1.19 $$ (3) 式中:H——滑坡高差/m;
A——滑坡面积/m2。
图7(d)展示了滑坡的H/L与滑坡面积(A)之间的关系。在滑坡动力学研究中,H/L是衡量滑坡能量转化与运动特性的重要指标之一。从理论上讲,当H/L值较小时,意味着在相对较小的高差下能够实现较长的滑行距离,这表明滑坡体在运动过程中能够更有效地将重力势能转化为动能,并且在运动过程中受到的阻力相对较小。因此在一些研究中,该指标也被用于反映滑坡的流动性[53 − 56]。此次事件触发的滑坡的H/L主要分布在0.2~0.6,包括
3563 处滑坡,占比61.65%,呈现出区域浅层滑坡的强流动性特征。在散点图中,大部分滑坡均位于我们得到的包络线之下(式4)。在这之中,滑坡的H/L与面积之间存在一定的负相关关系。这意味着,随着滑坡面积的增加,H/L有减小的趋势。
$$ H/L=-0.63\mathrm{lg}A+2.67 $$ (4) 式中:H——滑坡高差/m;
L——滑行距离/m;
A——滑坡面积/m2。
3.3 控制因素分析
3.3.1 滑坡数量和滑坡面积
滑坡数量和滑坡总面积在不同高程区间内的分布存在显著差异,见图8(a)。滑坡主要集中在海拔150~350 m,其中200~250 m区间内的滑坡数量(
1369 个)和滑坡总面积(约1.446 3×106 m2)达到峰值。从总体趋势来看,滑坡数量和总面积在150~350 m呈现快速增加,随后在350 m以上逐渐减少。在500 m以上的高海拔区域,滑坡数量和总面积显著下降。坡度是影响滑坡发生的重要地形因子,直接决定了地表重力分量、降雨汇流速度以及土壤稳定性等关键特性。一般来说,坡度较大的区域更容易因降雨诱发滑坡,因为陡峭的坡面使得抗剪力不足以抵抗重力和孔隙水压力的增大[57 − 58]。从统计分析结果来看,见图8(b),滑坡主要集中在坡度16°~28°,其中20°~24°区间的滑坡数量(
1228 个)和滑坡总面积(约1.183 1×106 m2)均达到峰值。本文分析该现象主要由以下原因导致:从图4(b)可以看出,在研究区范围内,高坡度区间面积占比非常小,尽管在理论上应存在更高的滑坡概率,但实际的滑坡数量和面积受限于其实际面积的局限。此外,在其他广东省区域的降雨滑坡研究中,包括韶关市[37]、揭西县[14]、陆河县[59]、龙川县[10, 60]的研究中,都存在滑坡主要分布于坡度45°以下范围的结论。综合现有研究,研究区基岩主要由花岗岩和片麻岩组成,上覆地层以第四纪松散堆积物和风化残积土为主,其他区域具有相似的岩性组成。关于此类岩性的滑坡失稳机制研究中提到,降雨过程中,由于降雨强度远高于花岗岩残土渗透系数,强降雨主要影响边坡浅表层,未入渗的雨水产生地表径流,通过浅表面滑动破坏。在斜坡浅表饱和区产生饱和渗流场,并演变成分布不均匀的特征,即在斜坡中下部饱和区厚度大于上部饱和区的厚度,因此坡体中下部变形破坏较多,对应坡度较缓的区域[60 − 61]。此外广东省对坡度大于25°的区域规定了林木采伐限制,也是该现象的支持原因之一[7]。已有研究表明,山区迎风坡地带往往具有更高的降雨量[62 − 64],因此能够接收更多降雨的水动力影响以及其坡面汇水能力增强,因此相对其他坡向,迎风坡滑坡分布更为集中[33, 65]。本研究中,滑坡数量和面积在东坡(E)和东南坡(SE)上达到峰值,其中东坡滑坡数量为
1082 个,总面积为96.78×104 m2,东南坡滑坡数量为1026 个,总面积为95.97×104 m2。此外,南坡(S)滑坡数量为976个,总面积略高于东坡和东南坡,达到97.20×104 m2。相较之下,北坡(N)、西北坡(NW)和东北坡(NE)的滑坡数量和总面积显著较少,见图8(c)。这同样与周边其他地区的研究具有相似的趋势[10, 14]。广东省属于亚热带季风气候,副热带高压、暖湿夏季风、北方冷空气和台风是影响广东省整体降雨情况的关键因素[64]。4月为华南前汛期阶段[66],南海夏季风还未暴发,因此偏南的暖湿气流是该地区的主导风向[67 − 68]。因此该现象符合暴雨期间山体迎风坡向滑坡分布远多于背风坡向的规律。3.3.2 剖面曲率
剖面曲率用于描述地表在沿坡面垂直于坡度方向的曲率特征。它反映了地表在沿坡面的切线方向上的弯曲程度,通常用此曲率来表征流体在重力作用下在表面上的向下加速和减速,因此可以反映降雨地表径流的汇集或分散程度[69 − 70]。剖面曲率为0~0.005内滑坡总面积最大,约为200.65×104 m2,而滑坡数量在−0.005至0区间最高,为2 057个,见图8(d)。在更典型的凹形地形中(剖面曲率为−0.01至−0.005)中,滑坡数量较多,共计904个,但其滑坡总面积较小,约66.92×104 m2。这表明凹形地形更易形成小规模滑坡,可能是由于汇水作用增强了局部水动力条件,导致坡体局部失稳而触发较多小型滑坡。而在典型凸形地形(剖面曲率为0.005至0.01)中,滑坡数量减少,为734个,但滑坡总面积更大,约为94.54×104 m2,说明凸形地形虽然滑坡数量较少,但倾向于形成规模更大的滑坡。这可能是因为凸形地形的水分分散条件较好,土壤水分更易沿坡体流失,从而需要更大的外界触发条件(如强降雨或高坡度)才能导致滑坡,而一旦触发,滑坡的规模则更大。
3.3.3 地形湿度指数
地形湿度指数主要用于评估土壤湿度的空间分布影响[71 − 72]。从统计结果来看,见图8(e),地形湿度指数值在4至5的区间内,滑坡数量和总面积均达到峰值,分别为2 695个和约283.64×104m2。这表明该区间是滑坡的高发区。
3.3.4 相对坡位指数
相对坡位指数通常是滑坡发生的起点位置。图8(f)表明,滑坡主要集中在相对坡位指数值接近0的区间(−0.2~0.2),对应斜坡中部区域。其中相对坡位指数为0~0.2的区间滑坡总面积最大,约为1.338 2×106 m2,滑坡数量为
1301 个;而在相对坡位指数为−0.2~0的区间滑坡数量最多,为1340 个,总面积为1.181 5×106 m2。3.3.5 距河流距离
河流在降水诱发滑坡中扮演着重要角色,其主要影响机制包括水动力侵蚀作用、地下水位变化以及河流汇流对边坡稳定性的影响。距离河流的远近直接决定了河流对坡体侵蚀和潜在滑坡影响的强弱[73 − 74]。从统计结果来看,见图8(g),滑坡数量和滑坡总面积随距离河流的增加而呈现显著递减趋势。距离河流100 m以内的区域内滑坡最为集中,滑坡数量为
1396 个,总面积约为1.277 8×106 m2。这表明河流附近的边坡因受侵蚀和冲刷作用影响较大,坡体稳定性显著降低,成为滑坡的高发区域。3.3.6 距道路距离
道路是影响滑坡发生的重要人为因素,其主要作用包括对坡体的直接切割、排水系统的改变以及地质结构的扰动。道路施工和交通荷载会显著降低坡体的稳定性。图8(h)显示,在距离道路100 m内的区域,滑坡数量为904个,总面积约为78.16×104 m2,距离道路100~200 m的区域滑坡数量有所减少,为740个,但滑坡总面积略有增加,达到约81.71×104 m2。随着距离的增加,滑坡数量和面积逐步减少。尽管研究区土地利用类型丰富,但是实际上滑坡主要集中分布的土地利用类型较少。因此我们对原有细分的类型进行了合并。图8(i)显示,阔叶林区域滑坡最多,数量达到
4102 个,总面积约为3.802 9×106 m2,占总滑坡分布的绝大多数。此外耕地区域滑坡数量为378个,总面积约为32.36×104 m2。这也反映了此次降雨时间对农业经济的严重影响。4. 讨论
近年来,随着全球气候变暖和极端天气事件的日益频繁,强降雨诱发的地质灾害正呈现出更高的发生频率和更大的破坏性[3, 75 − 76]。相关研究中,滑坡编目数据库的完善与标准化是未来的重要研究方向。在我国西南山区,研究者对2017年6月8日,重庆西部平行岭谷地区的暴雨事件进行分析,得到滑坡487处[15]。得到2023年7月4日重庆万州暴雨触发滑坡946处[77]。在黄土高原地区,研究表明2013年天水地区强降雨过程中触发黄土滑坡
54000 处[78 − 79]。在东南沿海地区,相关研究区域涉及了福建[80 − 81]、广东[10, 14, 37]、浙江[65, 82 − 83]和安徽[84 − 85]多省。图9展示了部分具有明确滑坡数量和面积的研究。降雨滑坡的规模效应在本研究中有所体现,但尚需更广泛的数据支撑。图9中的拟合趋势线表明,滑坡数量与滑坡面积之间存在明显的线性正相关关系。但是区域之间的差异性依然存在。这种差异可能与区域降雨强度、地形地貌特征以及土地利用类型等控制因素密切相关。未来研究应结合遥感解译与地质调查,进一步揭示不同区域滑坡触发的主控因素及其对滑坡规模的影响。例如量化降雨强度-历时与滑坡规模的关系,将是下一步研究的重要方向。气候变化对极端降雨事件的驱动作用需要更多关注。近年来,东南沿海地区受台风和梅雨季节的影响,极端降雨事件愈发频繁[87 − 88]。台风-暴雨-滑坡灾害链效应也逐步加剧,图中多个东南沿海省份的事件的时间分布也体现了这一趋势。这表明未来可能需要更多地将气象预报与滑坡灾害防控相结合,构建动态监测和实时预警体系,从而减少极端天气带来的损失。此外,尽管当前降雨触发滑坡数据库的研究逐步丰富,但在滑坡定义、编目精度和分类方法上可能存在一定的差异。为提高研究的对比性和适用性,未来应推动滑坡数据库的标准化和开放共享,从而为更高效的滑坡灾害预测模型提供基础。
5. 结论
(1)滑坡空间分布特征。本次暴雨诱发的浅层滑坡总计
6310 处,滑坡分布密度在北东—南西向的某些区域显著升高,表现出强烈的空间聚集性。滑坡总面积达5.85 km2,单个滑坡面积以102~103 m2的小规模为主,占比66.7%。大面积滑坡相对较少,主要分布在局部高坡度区域。(2)滑坡形态特征。滑坡的几何形态特征表明,滑行距离与滑坡面积呈幂律相关,滑坡长宽比主要集中在1~4,占比67.23%。滑坡高差大多数小于20 m,但高差较大的滑坡通常具有更长的滑行距离。滑坡的H/L主要分布在0.2~0.6,反映了浅层滑坡的高流动性特征,且滑坡规模增大时H/L呈下降趋势。
(3)控制因素分析。滑坡的发生受到多种地形因子和环境条件的综合影响。高程(150~350 m)和坡度(16°~28°)是滑坡的高发区间,其中高程200~250 m和坡度20°~24°的滑坡数量和面积均达到峰值。东坡和东南坡的滑坡数量和面积显著高于其他坡向,这与坡面的降雨汇水能力和水动力条件密切相关。地形湿度指数在4~5区间时滑坡最为集中,表明中等湿度条件下坡体更易失稳。滑坡数量随距离河流和道路的增加显著递减。距河流100 m以内区域滑坡数量最多,为
1396 处,显示出河流侵蚀作用对边坡稳定性的影响显著。距道路100 m以内的滑坡数量为904处,表明道路施工引发的边坡失稳问题需特别关注。本研究通过目视精确解译构建了此次降雨事件触发的详细滑坡编目,统计并揭示了自然因素和人为活动对滑坡触发的影响,为滑坡灾害风险评估提供了重要支持。未来研究需进一步结合动态降雨强度—历时模型与地质条件,深入揭示极端天气背景下滑坡的触发机制与空间规律。同时,建议加强区域性滑坡数据库的标准化建设,并发展基于气象数据的实时监测与预警系统,以应对气候变化带来的地质灾害挑战。
-
表 1 基于易发性与时间概率等级的区域滑坡危险性评价表
Table 1 Regional landslide hazard assessment table based on susceptibility and temporal probability levels
时间等级 易发性 极低
易发性低
易发性中
易发性高
易发性极高
易发性P1(0<P(x)≤P1) 极低
危险性极低
危险性极低
危险性极低
危险性低
危险性P2(P1<P(x)≤P2) 极低
危险性极低
危险性低
危险性低
危险性中
危险性P3(P2<P(x)≤P3) 极低
危险性低
危险性中
危险性中
危险性高
危险性P4(P3<P(x)≤P4) 极低
危险性低
危险性中
危险性高
危险性极高
危险性P5(P4<P(x)≤1) 低
危险性中
危险性高
危险性极高
危险性极高
危险性表 2 评价因子分级结果
Table 2 The grading results of assessment factors
评价
因子分级 分级面积/km2 Si/S
(×100)滑坡面积/km2 Ni/N
(×100)I 坡度/(°) 0~10 997.6 33.9 1.8 22.8 −0.4 10~20 996.8 33.9 3.7 46.8 0.3 20~30 629.7 21.4 2.0 25.3 0.2 30~40 246.3 8.4 0.4 5.1 −0.5 >40 71.8 2.4 0.0 0.0 −5.3 坡向/
(°)北(337.5~22.5) 168.2 5.7 0.3 3.8 −0.4 东北(22.5~67.5) 67.7 2.3 1.5 19.0 2.1 东(67.5~112.5) 320.3 10.9 1.0 12.7 0.2 东南(112.5~157.5) 337.5 11.5 0.8 10.1 −0.1 南(157.5~202.5) 294.5 10.0 0.8 10.1 0.0 西南(202.5~247.5) 313.1 10.6 0.8 10.1 0.0 西(247.5~292.5) 352.2 12.0 1.3 16.5 0.3 西北(292.5~337.5) 400.8 13.6 1.1 13.9 0.0 平面(−1) 687.7 23.4 0.3 3.8 −1.8 曲率 <−9 100.9 3.4 0.0 0.0 −5.6 −9~−6 172.7 5.9 0.1 1.3 −1.5 −6~−3 1201.1 40.8 3.4 43.0 0.1 −3~0 550.7 18.7 0.3 3.8 −1.6 0~3 550.7 18.7 4.1 51.9 1.0 3~6 176.0 6.0 0.0 0.0 −6.2 6~12 92.9 3.2 0.0 0.0 −5.5 >12 97.7 3.3 0.0 0.0 −5.6 高程/
m<200 129.4 4.4 1.9 24.1 1.7 200~300 377.1 12.8 2.6 32.9 0.9 300~400 449.6 15.3 1.6 20.3 0.3 400~500 417.2 14.2 0.7 8.9 −0.5 500~600 377.3 12.8 0.4 5.1 −0.9 600~700 435.3 14.8 0.4 5.1 −1.1 700~800 375.7 12.8 0.2 2.5 −1.6 800~900 124.1 4.2 0.1 1.3 −1.2 900~ 1000 72.7 2.5 0.0 0.0 −5.3 > 1000 184.1 6.3 0.0 0.0 −6.2 地形湿度指数 <10 558.1 19.0 0.9 11.4 −0.5 10~20 825.4 28.1 2.2 27.8 0.0 20~30 449.1 15.3 1.6 20.3 0.3 30~40 215.4 7.3 1.2 15.2 0.7 40~50 108.0 3.7 0.5 6.3 0.5 >50 786.5 26.7 1.5 19.0 −0.3 地表粗
糙度<1.1 2346.2 79.7 6.7 84.8 0.1 1.1~1.2 400.9 13.6 1.0 12.7 −0.1 1.2~1.3 120.8 4.1 0.2 2.5 −0.5 1.3~1.4 55.9 1.9 0.0 0.0 −5.0 1.4~1.5 15.8 0.5 0.0 0.0 −3.7 >1.5 2.7 0.1 0.0 0.0 −2.0 地层 梁山组+栖霞组+
茅口组并层26.1 0.9 0.0 0.1 −1.9 大冶组与嘉陵江
组并层533.8 18.1 0.3 3.8 −1.6 巴东组 212.6 7.2 1.0 12.7 0.6 地层 吴家坪组与长兴
组并层48.8 1.7 0.0 0.0 −4.9 韩家店组 9.0 0.3 0.0 0.0 −3.2 须家河组 137.7 4.7 0.4 5.1 0.1 珍珠冲组 122.2 4.2 0.7 8.9 0.8 自流井组 85.3 2.9 0.7 8.9 1.1 新田沟组 85.4 2.9 0.2 2.5 −0.1 沙溪庙组 841.7 28.6 2.4 30.4 0.1 龙马溪组与小河坝组并层 2.2 0.1 0.0 0.0 −1.8 蓬莱镇组 367.7 12.5 0.0 0.5 −3.2 遂宁组 469.8 16.0 2.2 27.8 0.6 岩层倾
角/(°)<10 632.0 21.5 1.8 22.8 0.1 10~20 1096.4 37.3 3.9 49.4 0.3 20~30 443.6 15.1 1.1 13.9 −0.1 30~40 364.9 12.4 0.6 7.6 −0.5 40~50 342.2 11.6 0.5 6.3 −0.6 50~60 58.2 2.0 0.0 0.0 −5.1 >60 5.2 0.2 0.0 0.0 −2.6 岩层倾
向/(°)北(337.5~22.5) 376.3 12.8 1.0 12.9 0.0 东北(22.5~67.5) 319.5 10.8 0.7 8.5 −0.2 东(67.5~112.5) 367.2 12.5 1.0 12.5 0.0 东南(112.5~157.5) 364.8 12.4 0.8 9.6 −0.3 南(157.5~202.5) 305.1 10.3 0.5 6.6 −0.4 西南(202.5~247.5) 313.5 10.6 1.5 18.4 0.6 西(247.5~292.5) 427.1 14.5 1.5 18.7 0.3 西北(292.5~337.5) 476.2 16.1 1.0 12.8 −0.2 距断层距离/km <1 101.2 3.4 0.5 6.3 0.6 1~2 111.9 3.8 0.6 7.6 0.7 2~3 128.4 4.4 0.2 2.5 −0.5 3~4 141.3 4.8 0.4 5.1 0.1 4~5 154.4 5.2 0.3 3.8 −0.3 >5 2305.2 78.3 5.9 74.7 0.0 距水系距离/m 0~200 173.6 5.9 2.6 32.4 1.7 200~400 102.1 3.5 1.2 14.9 1.5 400~600 104.7 3.6 0.2 2.5 −0.4 600~800 94.5 3.2 0.2 2.9 −0.1 800~ 1000 102.9 3.5 0.2 2.9 −0.2 1000 ~1200 84.8 2.9 0.1 0.9 −1.2 1200 ~1400 100.2 3.4 0.1 0.8 −1.4 1400 ~1600 90.9 3.1 0.7 8.9 1.1 > 1600 2087.6 71.0 2.7 33.9 −0.7 距道路距离/m 0~200 310.1 10.5 2.1 28.5 1.0 200~400 197.5 6.7 0.5 6.6 0.0 400~600 177.1 6.0 0.8 11.2 0.6 600~800 144.4 4.9 0.4 5.1 0.0 800~ 1000 140.5 4.8 0.2 2.2 −0.8 1000 ~1200 110.4 3.8 0.2 2.9 −0.3 1200 ~1400 122.2 4.2 0.3 4.0 0.0 > 1400 1739.1 59.1 3.4 44.9 −0.3 土层厚
度/m0~2.5 2503.5 85.1 4.0 50.6 −0.5 2.5~5 272.4 9.3 1.3 16.5 0.6 5~7.5 30.5 1.0 0.9 11.4 2.4 7.5~10 69.8 2.4 0.9 11.4 1.6 10~12.5 18.6 0.6 0.2 2.5 1.4 12.5~15 19.6 0.7 0.2 2.5 1.3 15~17.5 13.8 0.5 0.4 5.1 2.4 17.5~20 9.3 0.3 0.0 0.0 −3.2 >20 5.2 0.2 0.0 0.0 −2.6 归一化植被指数 <−0.1 53.1 1.8 0.6 7.6 1.4 −0.1~0 15.8 0.5 0.2 2.5 1.5 0~0.1 24.1 0.8 0.1 1.3 0.4 0.1~0.2 62.9 2.1 0.1 1.3 −0.5 0.2~0.3 135.6 4.6 0.5 6.3 0.3 0.3~0.4 720.1 24.5 1.7 21.5 −0.1 0.4~0.5 1624.4 55.2 3.7 46.8 −0.2 >0.5 306.1 10.4 1.0 12.7 0.2 土地利用类别 耕地 1452.9 49.4 5.3 67.5 0.3 林地 726.6 24.7 0.4 4.7 −1.7 草地 512.1 17.4 1.3 16.6 −0.1 灌木地 67.1 2.3 0.1 1.2 −0.7 湿地 0.3 0.0 0.0 0.1 2.1 水体 88.7 3.0 0.0 0.1 −3.1 人造地表 94.4 3.2 0.8 9.8 1.1 表 3 雨量站月平均降雨量
Table 3 Average monthly rainfall of rainfall stations
雨量站 月平均降雨量/mm 长寿 95.6 涪陵 91.4 丰都 91.6 凤来 106.4 雨台山 81.9 睦和 87.5 双河口 104.2 大木 117.6 河图 75.6 武陵山 89.3 龙潭 102.5 表 4 时间概率等级与月平均降雨量关系表
Table 4 Relationship between temporal probability levels and monthly average rainfall
时间概率等级 时间概率 月平均降雨量/mm P1 0<P(x)≤0.1 <78 P2 0.1<P(x)≤0.25 78~87 P3 0.25<P(x)≤0.5 87~96 P4 0.5<P(x)≤0.75 96~106 P5 0.75<P(x)≤1 >106 表 5 30次滑坡样本发生时间概率等级统计表
Table 5 Temporal probability levels of occurrence for 30 landslide samples
时间概率等级 P1 P2 P3 P4 P5 样本数 2 2 6 12 8 表 6 4次典型滑坡信息
Table 6 Information on four typical landslides
滑坡事件 日期 月平均降雨量/mm 时间概率等级 滑坡1 2020-06-29 98.7 P4 滑坡2 2018-08-01 103.4 P4 滑坡3 2020-06-27 118.2 P5 滑坡4 2016-06-02 121.3 P5 -
[1] 谢家龙,李远耀,王宁涛,等. 考虑库水位及降雨联合作用的云阳县区域滑坡危险性评价[J]. 长江科学院院报,2021,38(12):72 − 81. [XIE Jialong,LI Yuanyao,WANG Ningtao,et al. Assessment of regional landslide hazard in Yunyang County considering the combined effect of reservoir water level and rainfall[J]. Journal of Yangtze River Scientific Research Institute,2021,38(12):72 − 81. (in Chinese with English abstract)] XIE Jialong, LI Yuanyao, WANG Ningtao, et al. Assessment of regional landslide hazard in Yunyang County considering the combined effect of reservoir water level and rainfall[J]. Journal of Yangtze River Scientific Research Institute, 2021, 38(12): 72 − 81. (in Chinese with English abstract)
[2] LIN Qigen,LIMA P,STEGER S,et al. National-scale data-driven rainfall induced landslide susceptibility mapping for China by accounting for incomplete landslide data[J]. Geoscience Frontiers,2021,12(6):101248. DOI: 10.1016/j.gsf.2021.101248
[3] LI Langping,LAN Hengxing,GUO Changbao,et al. A modified frequency ratio method for landslide susceptibility assessment[J]. Landslides,2017,14(2):727 − 741. DOI: 10.1007/s10346-016-0771-x
[4] HUANG Faming,YIN Kunlong,HUANG Jinsong,et al. Landslide susceptibility mapping based on self-organizing-map network and extreme learning machine[J]. Engineering Geology,2017,223:11 − 22. DOI: 10.1016/j.enggeo.2017.04.013
[5] HUANG Faming,ZHANG Jing,ZHOU Chuangbing,et al. A deep learning algorithm using a fully connected sparse autoencoder neural network for landslide susceptibility prediction[J]. Landslides,2020,17(1):217 − 229. DOI: 10.1007/s10346-019-01274-9
[6] GUO Zizheng,YIN Kunlong,GUI Lei,et al. Regional rainfall warning system for landslides with creep deformation in Three Gorges using a statistical black box model[J]. Scientific Reports,2019,9:8962. DOI: 10.1038/s41598-019-45403-9
[7] NGUYEN B Q V,KIM Y T. Regional-scale landslide risk assessment on Mt. Umyeon using risk index estimation[J]. Landslides,2021,18(7):2547 − 2564. DOI: 10.1007/s10346-021-01622-8
[8] HUANG Faming,CAO Zhongshan,GUO Jianfei,et al. Comparisons of heuristic,general statistical and machine learning models for landslide susceptibility prediction and mapping[J]. Catena,2020,191:104580. DOI: 10.1016/j.catena.2020.104580
[9] HE Qian,JIANG Ziyu,WANG Ming,et al. Landslide and wildfire susceptibility assessment in southeast Asia using ensemble machine learning methods[J]. Remote Sensing,2021,13(8):1572. DOI: 10.3390/rs13081572
[10] WANG Zitao,LIU Qimeng,LIU Yu. Mapping landslide susceptibility using machine learning algorithms and GIS:A case study in Shexian County,Anhui Province,China[J]. Symmetry,2020,12(12):1954. DOI: 10.3390/sym12121954
[11] 周超,甘露露,王悦,等. 综合非滑坡样本选取指数与异质集成机器学习的区域滑坡易发性建模[J]. 地球信息科学学报,2023,25(8):1570 − 1585. [ZHOU Chao,GAN Lulu,WANG Yue,et al. Landslide susceptibility prediction based on non-landslide samples selection and heterogeneous ensemble machine learning[J]. Journal of Geo-Information Science,2023,25(8):1570 − 1585. (in Chinese with English abstract)] ZHOU Chao, GAN Lulu, WANG Yue, et al. Landslide susceptibility prediction based on non-landslide samples selection and heterogeneous ensemble machine learning[J]. Journal of Geo-Information Science, 2023, 25(8): 1570 − 1585. (in Chinese with English abstract)
[12] 贾雨霏,魏文豪,陈稳,等. 基于SOM-I-SVM耦合模型的滑坡易发性评价[J]. 水文地质工程地质,2023,50(3):125 − 137. [JIA Yufei,WEI Wenhao,CHEN Wen,et al. Landslide susceptibility assessment based on the SOM-I-SVM model[J]. Hydrogeology & Engineering Geology,2023,50(3):125 − 137. (in Chinese with English abstract)] JIA Yufei, WEI Wenhao, CHEN Wen, et al. Landslide susceptibility assessment based on the SOM-I-SVM model[J]. Hydrogeology & Engineering Geology, 2023, 50(3): 125 − 137. (in Chinese with English abstract)
[13] 李国营,刘平,张凯,等. 量纲统一在滑坡易发性评价中的影响分析[J]. 水文地质工程地质,2024,51(3):118 − 129. [LI Guoying,LIU Ping,ZHANG Kai,et al. Analysis of the influence of dimensional unity in landslide susceptibility assessment[J]. Hydrogeology & Engineering Geology,2024,51(3):118 − 129. (in Chinese with English abstract)] LI Guoying, LIU Ping, ZHANG Kai, et al. Analysis of the influence of dimensional unity in landslide susceptibility assessment[J]. Hydrogeology & Engineering Geology, 2024, 51(3): 118 − 129. (in Chinese with English abstract)
[14] 易靖松,王峰,程英建,等. 高山峡谷区地质灾害危险性评价——以四川省阿坝县为例[J]. 中国地质灾害与防治学报,2022,33(3):134 − 142. [YI Jingsong,WANG Feng,CHENG Yingjian,et al. Study on the risk assessment of geological disasters in alpine valley area:A case study in Aba County,Sichuan Province[J]. The Chinese Journal of Geological Hazard and Control,2022,33(3):134 − 142. (in Chinese with English abstract)] YI Jingsong, WANG Feng, CHENG Yingjian, et al. Study on the risk assessment of geological disasters in alpine valley area: A case study in Aba County, Sichuan Province[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33(3): 134 − 142. (in Chinese with English abstract)
[15] 陈洪凯,魏来,谭玲. 降雨型滑坡经验性降雨阈值研究综述[J]. 重庆交通大学学报(自然科学版),2012,31(5):990 − 996. [CHEN Hongkai,WEI Lai,TAN Ling. Review of research on empirical rainfall threshold of rainfall-induced landslide[J]. Journal of Chongqing Jiaotong University (Natural Science),2012,31(5):990 − 996. (in Chinese with English abstract)] CHEN Hongkai, WEI Lai, TAN Ling. Review of research on empirical rainfall threshold of rainfall-induced landslide[J]. Journal of Chongqing Jiaotong University (Natural Science), 2012, 31(5): 990 − 996. (in Chinese with English abstract)
[16] 许强. 对滑坡监测预警相关问题的认识与思考[J]. 工程地质学报,2020,28(2):360 − 374. [XU Qiang. Understanding the landslide monitoring and early warning:Consideration to practical issues[J]. Journal of Engineering Geology,2020,28(2):360 − 374. (in Chinese with English abstract)] XU Qiang. Understanding the landslide monitoring and early warning: Consideration to practical issues[J]. Journal of Engineering Geology, 2020, 28(2): 360 − 374. (in Chinese with English abstract)
[17] ROSI A,PETERNEL T,JEMEC-AUFLIČ M,et al. Rainfall thresholds for rainfall-induced landslides in Slovenia[J]. Landslides,2016,13(6):1571 − 1577. DOI: 10.1007/s10346-016-0733-3
[18] 赵伯驹,李宁,幸夫诚,等. 基于I-CF模型的四川德格县滑坡危险性评价与区划[J]. 中国地质灾害与防治学报,2023,34(5):32 − 42. [ZHAO Boju,LI Ning,XING Fucheng,et al. Landslide geological hazard assessment based on the I-CF model of Dege County in Sichuan Province[J]. The Chinese Journal of Geological Hazard and Control,2023,34(5):32 − 42. (in Chinese with English abstract)] ZHAO Boju, LI Ning, XING Fucheng, et al. Landslide geological hazard assessment based on the I-CF model of Dege County in Sichuan Province[J]. The Chinese Journal of Geological Hazard and Control, 2023, 34(5): 32 − 42. (in Chinese with English abstract)
[19] PRADHAN B,LEE S,BUCHROITHNER M F. A GIS-based back-propagation neural network model and its cross-application and validation for landslide susceptibility analyses[J]. Computers,Environment and Urban Systems,2010,34(3):216 − 235. DOI: 10.1016/j.compenvurbsys.2009.12.004
[20] L-B STATISTICS,BREIMAN L. Random forests[J]. Machine Learning,2001:5 − 32.
[21] 黄发明,曹中山,姚池,等. 基于决策树和有效降雨强度的滑坡危险性预警[J]. 浙江大学学报(工学版),2021,55(3):472 − 482. [HUANG Faming,CAO Zhongshan,YAO Chi,et al. Landslides hazard warning based on decision tree and effective rainfall intensity[J]. Journal of Zhejiang University (Engineering Science),2021,55(3):472 − 482. (in Chinese with English abstract)] HUANG Faming, CAO Zhongshan, YAO Chi, et al. Landslides hazard warning based on decision tree and effective rainfall intensity[J]. Journal of Zhejiang University (Engineering Science), 2021, 55(3): 472 − 482. (in Chinese with English abstract)
[22] 谭玉敏,郭栋,白冰心,等. 基于信息量模型的涪陵区地质灾害易发性评价[J]. 地球信息科学学报,2015,17(12):1554 − 1562. [TAN Yumin,GUO Dong,BAI Bingxin,et al. Geological hazard risk assessment based on information quantity model in Fuling District,Chongqing City,China[J]. Journal of Geo-Information Science,2015,17(12):1554 − 1562. (in Chinese with English abstract)] TAN Yumin, GUO Dong, BAI Bingxin, et al. Geological hazard risk assessment based on information quantity model in Fuling District, Chongqing City, China[J]. Journal of Geo-Information Science, 2015, 17(12): 1554 − 1562. (in Chinese with English abstract)
[23] 柳依莎,杨华. 基于信息量模型的地质灾害危险性评价研究——以重庆市涪陵区为例[J]. 重庆师范大学学报(自然科学版),2012,29(4):34 − 40. [LIU Yisha,YANG Hua. On the information content model of geological hazard assessment in the area of Fuling of Chongqing[J]. Journal of Chongqing Normal University (Natural Science),2012,29(4):34 − 40. (in Chinese with English abstract)] LIU Yisha, YANG Hua. On the information content model of geological hazard assessment in the area of Fuling of Chongqing[J]. Journal of Chongqing Normal University (Natural Science), 2012, 29(4): 34 − 40. (in Chinese with English abstract)