ISSN 1003-8035 CN 11-2852/P
  • 中国科技核心期刊
  • CSCD收录期刊
  • Caj-cd规范获奖期刊
  • Scopus 收录期刊
  • DOAJ 收录期刊
  • GeoRef收录期刊
欢迎扫码关注“i环境微平台”

基于滑坡破坏模式分析的易发性评价以三峡库区首段泄滩河左岸为例

朱宇航, 黄海峰, 殷坤龙, 郭子正, 郭飞, 赖鹏

朱宇航,黄海峰,殷坤龙,等. 基于滑坡破坏模式分析的易发性评价−以三峡库区首段泄滩河左岸为例[J]. 中国地质灾害与防治学报,2023,34(2): 156-166. DOI: 10.16031/j.cnki.issn.1003-8035.202112035
引用本文: 朱宇航,黄海峰,殷坤龙,等. 基于滑坡破坏模式分析的易发性评价−以三峡库区首段泄滩河左岸为例[J]. 中国地质灾害与防治学报,2023,34(2): 156-166. DOI: 10.16031/j.cnki.issn.1003-8035.202112035
ZHU Yuhang,HUANG Haifeng,YIN Kunlong,et al. Evaluation of landslide susceptibility based on landslide failure mode analysis: A case study of the left bank of Xietan River in the first section of Three Gorges Reservoir[J]. The Chinese Journal of Geological Hazard and Control,2023,34(2): 156-166. DOI: 10.16031/j.cnki.issn.1003-8035.202112035
Citation: ZHU Yuhang,HUANG Haifeng,YIN Kunlong,et al. Evaluation of landslide susceptibility based on landslide failure mode analysis: A case study of the left bank of Xietan River in the first section of Three Gorges Reservoir[J]. The Chinese Journal of Geological Hazard and Control,2023,34(2): 156-166. DOI: 10.16031/j.cnki.issn.1003-8035.202112035

基于滑坡破坏模式分析的易发性评价——以三峡库区首段泄滩河左岸为例

基金项目: 国家自然科学基金项目 (42107489);水电工程智能视觉监测湖北省重点实验室开放基金项目(2020SDSJ02);三峡库区地质灾害教育部重点实验室开放基金项目(2020KDZ09)
详细信息
    作者简介:

    朱宇航(1999-),男,硕士,主要从事地质灾害风险评价方面的研究工作。E-mail:cugzyh@cug.edu.cn

    通讯作者:

    黄海峰(1978-),男,博士,教授,主要从事地质灾害监测防治等方面的教学与科研工作。E-mail:hhf@ctgu.edu.cn

  • 中图分类号: P642.22

Evaluation of landslide susceptibility based on landslide failure mode analysis: A case study of the left bank of Xietan River in the first section of Three Gorges Reservoir

  • 摘要: 三峡库区首段发育有大量岩质滑坡,其中很多灾害点极具隐蔽性且目前并未被查明。文中以三峡库区首段泄滩河左岸为研究区,以区内唯一破坏的卡门子湾顺层岩质滑坡为例,在分析其成因机制的基础上归纳总结了该地区顺层岩质滑坡的破坏模式,并以此确定了高程、坡度、坡向、起伏度、平面曲率、剖面曲率、地层岩性、距河流距离及距道路距离共9个评价指标因子及疑似滑坡隐患点,将这些灾害隐患点作为滑坡样本,运用ALSA模型开展研究区的滑坡易发性分区,最后采用ROC曲线及现场复查等方法验证评价结果的可靠性。预测结果表明:研究区内顺层岩质滑坡的极高易发区和较高易发区大致呈面状分布,主要集中在岩性为侏罗系中统上沙溪庙组紫红色泥岩夹砂岩和西北坡向的近库岸地区。现场验证发现易发分区结果与滑坡破坏模式分布规律较吻合,表明基于滑坡破坏模式选择滑坡样本得到的滑坡易发性结果在整体上也能反映研究区滑坡概率空间分布规律,在缺乏准确滑坡样本时可作为一种替补方案。上述研究结果为基于滑坡破坏模式选取滑坡样本开展易发性评价工作提供了理论支持和科学依据。
    Abstract: There are a large number of rock landslide disasters developed in the first section of the Three Gorges Reservoir area, many of which are very hidden and have not been identified. In this paper, taking the left bank of Xietan River in the first section of the Three Gorges Reservoir as the study area, taking the only bedding rock landslide in Kamenziwan as an example, the failure mode of bedding rock landslide in this area is summarized on the basis of analyzing its genesis mechanism. Nine evaluation index factors, including elevation, slope aspect, slope, relief, plane curvature, profile curvature, formation lithology, distance from river and distance from road, as well as suspected hidden danger points of landslide disaster are determined. These hidden danger points are taken as landslide samples. Automatic Landside Susceptibility Assessment Model (ALSA) was used to carry out landslide Susceptibility zoning in the study area. Finally, ROC curve and field review were used to verify the reliability of the evaluation results. The prediction results show that the extremely high and highly prone areas of bedding rock landslides in the study area are distributed in a plane shape, mainly concentrated in the middle Jurassic Upper Shaximiao Formation purplish red mudstone intercalated sandstone, and the northwest slope direction near the reservoir bank area. Field verification shows that the results of prone zoning are consistent with the distribution law of landslide failure mode, indicating that the landslide susceptibility results obtained by selecting landslide samples based on landslide failure mode can also reflect the spatial distribution law of landslide probability in the study area on the whole, and can be used as a substitute scheme in the absence of accurate landslide samples. The above research results provide theoretical support and scientific basis for selecting landslide samples to carry out vulnerability assessment based on landslide failure mode.
  • 滑坡是我国发育最广泛的地质灾害之一,尤其是21世纪以来,随着人类工程活动的开展,滑坡发生频率上升,造成了巨大的人员伤亡和经济损失。对滑坡的风险评估和预测是目前防灾减灾的重要途径之一。随着计算机技术的发展,目前对滑坡风险评估逐渐由定性向定量转变,在定量评价中,堆积层厚度对滑坡识别具有重要意义,并且是一个重要的控制因素,能否对其进行准确预测是定量评价的关键。

    目前,国内外对于土壤厚度的预测研究较为成熟,主要有基于物理机制、基于点采样插值、基于遥感的方法和基于环境相关推理的方法[1]。李岩等[2]采用多种插值方法对草原栗钙土层厚度空间分布进行预测,并选出适合研究区的最优插值方法。王桂林等[3]基于粒子群优化对协同克里金法进行了优化,结果表明该法对提高土层厚度的预测精度有较大作用。Scarphone等[4]采用广义线性模型(GLM)、随机森林(RF)、GLM残差克里格法和RF残差克里格法预测土壤厚度,结果表明GLM残差克里格法是预测浅层土壤厚度的最佳方法。Pradhan等[5]基于相对起伏度开发了一种空间土壤厚度模型,结果表明该模型的预测能力要优于多元线性回归模型和S模型。Li等[6]采用多元线性回归(MLR)、支持向量回归(SVR)、随机森林(RF)和极端梯度推进(XGBoost)4种常用机器学习方法分别对河南省土壤厚度进行预测。

    堆积层厚度预测与土壤厚度预测有相似之处,区别在于堆积层厚度预测需要更多地考虑地质成因,国内对堆积层厚度的研究较少,多以空间插值为主,但对外界因素考虑不足,难以满足精度要求。刘磊等[7]根据已知区域堆积层厚度与其标志因子之间的关系,建立了基于随机森林算法的堆积层厚度模型,估算了整个万州主城区的堆积层厚度的空间分布。段贵娟等[8]以水平相关距离为基础数据,采用普通克里金法对南海区域的厚度分布进行了空间预测,预测效果良好。邵晨灿等[9]通过分析各成因类型土层厚度的影响因素,采用决策树模型建立影响因素与土层厚度之间的关系,推测了贵州省瓮安县的土层厚度分布。虽然目前国内外对于堆积层厚度的研究较少,但是对于土壤厚度的预测研究较为成熟,对堆积层厚度的研究可以借鉴土壤学中基于环境相关推理的方法,该方法假设相同的土壤景观单元发育相同的土壤类型,如Wang等[10]通过采用模糊C-均值法和模糊隶属函数构建土壤景观模型,从而预测中国三峡秭归地区土壤深度的空间变化。Akumu等[11]开发了一种基于GIS模糊逻辑建模方法的土壤深度建模和制图技术,制图精度达到94%,优于随机森林模型。Cheng等[12]提出了一种从单个土壤多边形中挖掘土壤-环境关系的方法,对传统土壤图进行了更新。基于环境相关推理法应用到堆积层厚度中为:相同地质环境下的堆积层厚度相似,可通过已知点的堆积层厚度与影响因子之间的关系预测整个区域的堆积层厚度。

    三峡库区万州区铁峰乡地质灾害以滑坡为主,探索区内堆积层厚度对于滑坡稳定性评价显得尤为重要。鉴于此,本文以铁峰乡为研究对象,在实地调查的基础上,首先总结了区域堆积层的成因模式,并据此添加区域堆积层厚度控制点;其次定量分析了堆积层厚度与影响因子之间的相关性,并基于Apriori算法挖掘并分析影响因子与堆积层厚度的关联准则;最后,选用三种机器学习模型对研究区堆积层厚度进行空间预测,并对预测结果进行对比分析。研究结果揭示了区域第四系堆积层的成因机制,为机器学习在堆积层厚度预测领域奠定了基础。

    数据挖掘流程主要由数据采集,数据预处理,数据离散化和关联准则挖掘四步组成。本文采用两步聚类算法及Apriori算法,通过数据挖掘流程处理堆积层厚度数据,得到堆积层厚度与影响因子间的关联准则。

    (1)数据采集:将堆积层厚度进行收集归纳,获得基础数据。

    (2)数据预处理:对已有数据进行审核、筛选、排序等处理措施。

    (3)数据离散化:由于Apriori算法只能处理离散型变量,故采用两步聚类法将堆积层厚度离散化。

    (4)关联准则挖掘:应用Apriori算法对离散后的变量进行关联规则集合选择,生成有效的关联规则。

    通过两步聚类算法对堆积层厚度类型进行转化,将分类所得的变量采用Apriori算法计算,挖掘堆积层厚度的关联准则。

    Apriori算法由Agrawal和Srikant提出[13],该算法首先生成高于最小支持度的频繁项目集,在第一步产生的频繁项目集中生成高于最小可信度的关联准则。

    频繁项目集是指对包含项目a的项集T,其支持度大于或等于用户指定的支持度阈值,即:

    |T(a)||T|minsupp (1)

    包含k个项目的频繁项目集称为频繁k项集,记为Lk。如图1上层的项目集ababcabcd,当满足最小支持度时均为频繁k项集。

    图  1  Apriori算法示意图
    Figure  1.  Schematic diagram of Apriori algorithm

    Apriori算法实现过程详见图1。首先,搜索出产生长度为1的频繁项集L1L1又用于产生长度为2的频繁项目集L2,如此循环,搜索所有的频繁项目集。

    从频繁项目集中产生简单的关联准则,按置信度大于置信度阈值的条件,选择出有效规则集合。对每个频繁项目集L,计算L所有非空子集L的置信度,如果CL'LL'大于用户指定的置信度阈值,即

    CL(LL)=|T(L)||T(L)|minconf (2)

    则生成关联准则L(LL)

    支持向量回归机模型(Support vector regression,SVR)由美国学者Vapnik[14]提出,主要用于解决非线性问题的回归,其基本思想是:将样本数据分为两部分,分别为训练样本和测试样本,通过事先选择的非线性映射将输入向量(训练样本)映射到高维度特征空间,利用其有限的样本数据进行训练计算,求得达到最佳拟合效果时在这个空间中构造出的最优决策函数模型,对比测试项以分析模型的测试效果。该模型以统计训练理论为基础,实现小样本条件下的统计规律和学习,特点为结构简单,且性能比传统的BP神经网络有明显的优势。因此,本模型在堆积层厚度预测方面具有一定的优越性。

    在SVM模型的参数优化方法中,灰狼算法(Grey wolf optimization,GWO)较粒子群算法(Particles swarm optimization,PSO)、遗传算法(Genetic algorithm,GA)等传统启发式算法而言,在研究组合式优化问题上已被证明具有显著优势,能够有效地解决多种参数的联合寻优问题。此外,GWO-SVM模型更适合小样本数据,且具有更好的学习效果。

    GWO优化SVM的基本思想是将SVM的惩罚因子C、核函数参数g和不敏感参数ε作为猎物,通过不断更新狼群位置,对目标猎物的最优解进行定位。利用得到的最优解构建SVM预测模型。本文设置的参数范围为:C=[0.01,100]g=[0.01,100]ε = [0.001,1]

    采用GWO-SVM耦合模型,对堆积层厚度进行预测,具体预测流程如下:

    (1)输入训练集和测试集。

    (2)初始化狼群数量、迭代次数,设置惩罚因子C、核函数参数g和不敏感参数ε的范围。

    (3)GWO以Cgε为猎物进行优化,达到最大迭代次数时输出最优值。

    (4)利用寻优得到的参数构建预测模型,利用该模型进行堆积层厚度预测。

    铁峰乡位于重庆市万州区北部山区(图2),长江上游,东邻九池乡,南靠石龙乡、石和乡,西连两汇乡,北接开州区,总面积为51.46 km2。属亚热带季风气候,四季分明,雨量充沛。地势总体上南北两侧高,中间海拔低,以构造侵蚀中低山-河谷地貌为主。区内出露地层有侏罗系中统沙溪庙组和新沟组、侏罗系下统自流井组和珍珠冲组、三叠系须家河组和巴东组。第四系覆盖物广,主要为残坡积物、崩滑堆积物和河流冲洪积物。研究区位于铁峰山背斜北西翼核部,背斜全长约114 km,走向50°~70°,轴线呈凸向NNW的弧形,轴部地层为三叠系巴东组,两翼为侏罗系中、下统。

    图  2  万州区位置及铁峰乡滑坡分布图
    Figure  2.  Location of Wanzhou District and distribution of landslides in Tiefeng Township

    铁峰乡发育的地质灾害类型包括堆积层滑坡、顺层岩质滑坡和岩质崩塌。根据地质灾害调查资料,研究区共发育历史滑坡灾害14处(堆积层滑坡12处,岩质滑坡2处),新生滑坡灾害12处(堆积层滑坡9处,岩质滑坡3处),堆积层滑坡主要分布于研究区南侧的顺向坡区段,目前多处于蠕滑变形状态;顺层岩质滑坡主要发育于自流井组和珍珠冲组地层内,以软弱夹层为主要的控滑结构,其稳定状态受公路开挖的影响极大。研究区崩塌灾害点主要分布于研究区斜交坡、北侧逆向坡结构的自然高陡边坡,以及公路开挖形成的高切坡之上,总计87处,以楔形体滑移式崩塌为主,占崩塌灾害点总量的56.3%,其次为坠落式崩塌,占比32.2%。

    本次野外调查共确定了192个堆积层厚度数据控制点。对于全区而言,这些控制点不足以完成铁峰乡的堆积层厚度预测。因此,需在野外调查的基础上分析堆积层的成因模式,并据此添加一定量的堆积层控制点。根据铁峰乡坡体结构和岩层组合关系,对堆积层的成因进行分析总结,将研究区第四系堆积层分为残坡积和崩滑堆积两大类,崩滑堆积又可分为顺向坡远程滑动堆积、顺向坡近距离滑动堆积、顺向坡侧向崩塌堆积和软硬相间岩层崩滑堆积四种成因模式。

    指地表岩石风化后残留在原地的堆积物(图3),可区分为机械风化作用和化学风化作用。

    图  3  残坡积示意图
    Figure  3.  Schematic diagram of residual slope

    (1)坡顶及原地残积物

    由于重力和流水的搬运,宏观上沟谷及山麓地区多为冲、洪积物分布。残积物范围较为局限,山顶及未受扰动的风化原地堆积较为明显。常见的差异为较软的岩性风化堆积厚度较大,较硬的岩性厚度较薄。

    (2)斜坡上部残积物

    主要为风化堆积的残积物,未受重力作用及流水侵蚀而发生崩滑运动。该处第四系堆积形态与原岩较为相似,层状结构与地层一致。多位于发生过崩滑区域的上部,且相对于中下部厚度较薄。

    (3)软硬相间岩性残积物

    由于岩性差异形成的风化差异,不同的岩性处风化厚度不同,堆积形态也有所差异。整体上表现为原生堆积物在软岩处风化强烈,地形较缓且堆积较厚;在硬质岩处风化较弱,地形较陡而堆积较薄。

    指因崩塌滑坡运动造成的堆积物,原始堆积物形态发生改变。可分为顺向坡远程滑动堆积、顺向坡近距离滑动堆积、顺向坡侧向崩塌堆积和软硬相间岩层崩滑堆积四种成因模式。

    (1)顺向坡远程滑动堆积

    原岩风化形成的堆积物发生远程滑动从而改变坡面原始形态。从纵剖面上来看,斜坡后部地形较缓,且堆积体较薄,主要为残留的滑坡堆积物和较薄的残坡积物;中前部由于滑坡下滑堆积,厚度较大,且易形成局部平台,主要物质为滑坡堆积物。若滑坡运动距离较大,滑坡堆积物可能覆盖坡底冲沟内的冲洪积物,甚至冲入河道,改变河道形态(图4)。

    图  4  顺向坡远程滑动堆积示意图
    Figure  4.  Schematic diagram of long-distance sliding accumulation along the slope

    (2)顺向坡近距离滑动堆积

    原岩风化形成的堆积物在内外地质作用的影响下沿软弱层面发生失稳破坏,改变原始地貌形态。原始地貌平面上为正常侵蚀堆积的山坡,沟谷之间过渡平衡,而发生滑动后横剖面中部可见明显的非正常侵蚀,堆积物下移而导致中间地形下错明显,平面上可见明显的洼地。剖面上原地形为自然过渡的缓坡,滑动部分导致前缘鼓胀、堆积较厚;后缘发生滑动出现地形下陷。坡顶为基岩,上部为风化后未完全滑动的残积物,中下部为风化后滑动堆积的滑坡堆积物,厚度较大,前部可能存在被滑坡堆积所部分覆盖的冲洪积物(图5)。

    图  5  顺向坡近距离滑动堆积示意图
    Figure  5.  Schematic diagram of close range sliding accumulation along the slope

    (3)顺向坡侧向崩塌堆积

    顺向坡两侧冲沟陡崖在侧向侵蚀的作用下发生崩塌,崩积物堆积于陡崖下部,从而改变原始地貌。近崖处地形较陡,堆积层厚度较薄;随陡崖向外地形逐渐变缓,堆积厚度增厚;由于崩积物运动能量有限,远处堆积厚度逐渐变薄直到停止运动。由于崩坡积物粒径差异大,其形成的堆积结构特点明显,粒径较大的崩积物运动距离较大,多堆积于中前部,形成整体上前粗后细的堆积特点(图6)。

    图  6  顺向坡侧向崩塌堆积示意图
    Figure  6.  Schematic diagram of lateral collapse and accumulation along the slope

    (4)软硬相间岩层崩滑堆积

    软硬相间地层且临空条件较好,存在明显的差异风化特征,导致地形上出现了多级陡缓相间的差异变化,且不同位置堆积厚度及范围不同。平台堆积物厚度与陡崖的剖面形态相关,可分为三种类型:①坡下坡上陡崖高度一致,侧向的侵蚀范围相当,堆积厚度一定,以崩积物为主;②坡下陡崖高于坡下时,坡下堆积范围较坡上大,且厚度较大,以崩积物为主;③当侧向侵蚀同时存在滑坡时,坡上以崩积物为主,坡下滑动区则下伏滑坡堆积物,上覆新生崩积物(图7)。

    图  7  软弱相间岩层分布模式示意图
    Figure  7.  Schematic diagram of the distribution mode of soft and weak alternating rock formation

    在研究区野外调查获得的192个堆积层厚度数据点的基础上,结合上述堆积层成因模式,添加包含0点和非0点在内的堆积层厚度控制点,累计共得到1211个数据样本点。根据堆积层成因模式和堆积特点,本文考虑了坡度、坡向、流向、地形起伏、剖面曲率、斜坡结构等13个影响因子。采用ArcGIS软件对研究区的1∶10000地形图进行栅格划分,并对上述13个因子进行提取并分级,统计各调查点的堆积层厚度,分析堆积层厚度的分布规律。本文选取坡度、坡向、流向、高程、地形起伏、剖面曲率、斜坡结构六个因子进行分析(图8)。

    图  8  各因子统计图
    Figure  8.  Statistical chart of each factor

    (1)坡度

    从坡度因子图可以看出,研究区坡度主要集中在10°~40°范围内,该坡度范围内堆积层平均厚度为2.5~3.0 m。从堆积层平均厚度可以看出,堆积层厚度与坡度成反比,在坡度大于50°时,堆积层厚度接近0,而在0~40°范围内,堆积层厚度大于1.0 m,这表明在坡度低的地方,土层更容易堆积保留,从而形成一定的堆积层厚度,见图8(a)。

    (2)坡向

    从坡向因子图中可以看出,研究区内整体呈现出越接近于南北向的地方土层厚度越大的规律,在坡向近东西的区域内,堆积层的厚度普遍较小,通常小于0.5 m,在坡向为南北向的区域内堆积层的平均厚度达到最大值,通常在1.5~2 m。此规律与研究区内的地形地貌相吻合,研究区内坡向基本都是沿南北向分布,土体的崩滑运动主要也是沿这个方向发展,从而使土体在此方向上的堆积厚度大于其他方向,见图8(b)。

    (3)流向

    从流向因子图可以看出,堆积层平均厚度在大体上呈现出越容易向北汇水的区域堆积层厚度越厚,这与坡向因子相似,当堆积层点向北汇水时,该点更可能位于顺向坡内,而顺向坡更易发生崩滑运动,从而造成向北汇水点的堆积层厚度更厚,见图8(c)。

    (4)地形起伏

    从地形起伏因子可以看出,堆积层厚度分布规律与坡度因子完全一致,与地形起伏程度呈反比,地形起伏越小的区域内,土体越容易堆积,堆积层平均厚度在2.5~3 m之内,而当地形起伏达到4 m时,堆积层平均厚度仅为0.5 m,地形起伏超过5 m时,堆积层难以形成,平均厚度接近0。该规律与研究区地形相符合,由于研究区位于山区,地形起伏较大的区域往往也是坡度较大的区域,受地形限制,土体难以堆积,见图8(d)。

    (5)剖面曲率

    从剖面曲率因子可以看出,堆积层调查点主要集中在0~10、60~70、>70三个分级内。堆积层厚度分布呈现出与剖面曲率呈反比的规律,当剖面曲率位于0~10之内时,堆积层平均厚度达2.5~3 m,而当剖面曲率达到70时,厚度仅有0.5 m。剖面曲率大代表地表高程变化率大,即地势急剧变化,因此土层难以堆积保留,见图8(e)。

    (6)斜坡结构

    从斜坡结构因子可以看出,研究区主要为斜交坡和顺向坡。从堆积层厚度上看,土体在顺向坡中更容易堆积,从而形成一定的堆积层,顺向坡中的堆积层厚度在2.5~3 m之间,而斜交坡中堆积层厚度仅为1 m左右。这与地形作用息息相关,在顺向坡中,山体上部的土体由于崩滑等一系列作用从而向下移动,在顺向坡中往下堆积,从而在下部容易形成一定厚度的堆积层,见图8(f)。

    采用GRA(Grey relation analysis)和MIC(Maximal information coefficient)两种算法来分析上述13个影响因子与堆积层厚度的相关性,计算结果如表1所示。根据计算结果,各因子与堆积层厚度的GRA值差异不大,但是对于MIC算法,因子1、2、3、6、10、12与堆积层厚度的MIC值均大于0.2,表示具有较强的相关性。因此,本文选取以上6个因子(坡度、坡向、流向、地形起伏、剖面曲率、斜坡结构)作为关联前项,堆积层厚度作为关联后项,根据已知样本点数据挖掘两者的关联准则。

    表  1  各因子与堆积层厚度相关性指标
    Table  1.  Correlation index between each factor and accumulation layer thickness
    因子GRAMIC
    因子1坡度0.69230.6977
    因子2坡向0.68780.5285
    因子3流向0.86730.4696
    因子4高程0.70500.0752
    因子5汇水条件0.63490.0417
    因子6地形起伏0.68560.6933
    因子7地层0.75460.0486
    因子8距水系距离0.63660.0525
    因子9距陡崖距离0.60630.0637
    因子10剖面曲率0.69530.2141
    因子11平面曲率0.67830.0784
    因子12斜坡结构0.63640.3271
    因子13植被覆盖0.64400.0943
    下载: 导出CSV 
    | 显示表格

    由于数据挖掘中Apriori算法仅能处理离散型变量,因此,需要对影响因子进行聚类。采用两步聚类算法对坡度、坡向、流向、地形起伏、剖面曲率、斜坡结构6个因子进行聚类,如表2所示。堆积层厚度按照(0~2 m)、(2~4 m)、(4~6 m)、(6~8 m)、(>8 m)划分为HD1、HD2、HD3、HD4、HD5五类。

    表  2  各因子属性指标
    Table  2.  Attribute indicators of each factor
    因子属性类别
    因子1坡度0~20°F11
    20°~40°F12
    >40°F13
    因子2坡向337.5°~67.5°F21
    67.5°~157.5°F22
    157.5°~247.5°F23
    247.5°~337.5°F24
    因子3流向1&2F31
    4&8F32
    16&32F33
    64&128F34
    因子6地形起伏0~2F61
    2~4F62
    >4F63
    因子10剖面曲率0~20F101
    20~40F102
    40~60F103
    >60F104
    因子12斜坡结构顺向坡F121
    逆向坡F122
    斜交坡F123
    下载: 导出CSV 
    | 显示表格

    根据已知样本点数据,以上述六个影响因子作为关联前项,堆积层厚度作为关联后项,设定支持度阈值为0,置信度阈值为100%,以前项最大为3生成堆积层厚度与影响因子的关联准则。共生成关联准则193条(部分准则如表3所示)。从关联准则可以看出,当堆积层厚度大于2 m时,坡度小于40°、地形起伏小于4 m,而当堆积层厚度小于2 m时,坡度大于40°、地形起伏大于4 m,说明坡度和地形起伏是造成堆积层厚度空间分布差异的主要因素。

    表  3  堆积层厚度关联准则
    Table  3.  Association criterion of accumulation layer thickness
    规则 ID规则支持度/% 置信度/%
    1YZ3=F31 & YZ6=F62 & YZ2=F21HD50.08100.00
    2YZ3=F31 & YZ2=F21 & YZ1=F12HD50.08100.00
    3YZ2=F22 & YZ3=F34 & YZ1=F11HD30.08100.00
    4YZ2=F22 & YZ3=F34 & YZ6=F61HD30.08100.00
    5YZ3=F34 & YZ1=F12 & YZ2=F24HD20.17100.00
    6YZ2=F23 & YZ10=F102 & YZ1=F11HD20.08100.00
    7YZ10=F102 & YZ2=F22 & YZ1=F11HD20.08100.00
    8YZ3=F31 & YZ6=F63HD118.58100.00
    9YZ3=F31 & YZ1=F13 & YZ6=F63HD117.18100.00
    10YZ2=F22 & YZ1=F13 & YZ6=F63HD117.01100.00
    下载: 导出CSV 
    | 显示表格

    研究区共划分了486097个栅格点,其中1211个堆积层调查点,484886个预测栅格点。从1211个调查点中,随机选取1011个点作为训练样本,剩余200个作为测试样本,484886个未知堆积层厚度的栅格点作为预测样本,分别采用PSO-SVM、GA-SVM、GWO-SVM 3种机器学习模型对研究区堆积层厚度分布开展训练及预测。采用PSO、GA、GWO三种算法对SVM模型进行三参数寻优,各种优化方法的搜索区间如下:惩罚因子C=[0.01,100];核函数参数g=[0.01,100];不敏感参数ε = [0.001,1],其中GWO算法的寻优结果为:C=25.61,g=0.054,ε=0.131。各模型测试结果及精度如图9表4所示。

    图  9  各模型测试结果对比(测试样本)
    Figure  9.  Comparison of test results of each model (test samples)
    表  4  各模型预测精度对比
    Table  4.  Comparison of prediction accuracy of each model
    预测模型MSERMSEQC平均厚度/m
    PSO-SVM401.221.420.470.791.49
    GA-SVM198.571.000.230.761.76
    GWO-SVM152.240.870.190.831.83
    下载: 导出CSV 
    | 显示表格

    采用上述3种模型对区域堆积层厚度进行预测,预测结果如图10所示。根据图9表4可知,GWO-SVM模型的测试结果与实际最为吻合,预测结果MSE为152.24,RMSE为0.87,精度最佳。另外,采用GWO-SVM的全区堆积层厚度预测结果平均值为1.83 m,与现场实际调查的结果较为相符(根据调查点数据,堆积层平均厚度为1.43 m)。因此,选取GWO-SVM模型得到的区域堆积层厚度分布图为最终结果,如图10(c)所示。

    图  10  各模型的区域堆积层厚度预测结果
    Figure  10.  Prediction results of regional accumulation layer thickness of each model

    图10(c)可知,颜色较深区域为堆积层厚度较薄地区,对应于研究区坡度较陡以及陡崖带附近,与上述坡度越大堆积层厚度越小的规律相符;从近陡崖处到远离陡崖的位置,堆积层厚度呈现出薄-厚-薄的规律;总体沿南北方向呈条带状分布,与研究区地形地貌相契合,在坡脚处堆积较厚。

    (1)万州区铁峰乡第四系堆积层可分为残坡积和崩滑堆积两类,崩滑堆积有四种成因模式,分别为:顺向坡远程滑动堆积、顺向坡近距离滑动堆积、顺向坡侧向崩塌堆积和软硬相间岩层崩滑堆积。

    (2)堆积层厚度与六种因子表现出较强的相关性:坡向、坡度、流向、地形起伏、剖面曲率和斜坡结构。其中坡度和地形起伏是造成堆积层厚度空间分布差异的主要因素。

    (3)采用三种机器学习模型对铁峰乡堆积层厚度进行预测,结果显示GWO-SVM模型的预测结果与实际最为吻合,表明该方法可较好地应用于区域第四系堆积层厚度预测,为滑坡灾害识别提供更为可靠的依据。

  • 图  1   研究区位置图

    Figure  1.   Location map of the study area

    图  2   泄滩河左岸现场调查图(摄于2020年8月)

    Figure  2.   Field geotechnical photo of the left bank of Xietan River (August 2020)

    图  3   卡门子湾滑坡全貌及优势结构面

    Figure  3.   Arial view and dominant structural plane of Kamenziwan landslide

    图  4   卡门子湾滑坡结构面赤平投影

    Figure  4.   Equated projection of Kamenziwan landslide section

    图  5   卡门子湾滑坡破坏前后库水位和降雨量

    Figure  5.   Summery of reservoir water level and rainfall before and after Kamenziwan landslide failure

    图  6   卡门子湾滑坡破坏模式概化示意图

    Figure  6.   Generalized schematic diagram of failure mode of Kamenziwan landslide

    图  7   泄滩河北段粉砂岩山脊夹泥岩沟谷微地貌示意图

    Figure  7.   The terrain and watershed partition of siltstone ridge and mudstone valley in Xietan River north section

    图  8   易发性评价指标图层

    Figure  8.   Layer of susceptibility assessment index

    图  9   疑似滑坡区为训练样本的滑坡易发性分区图

    Figure  9.   Distribution map of landslide susceptibility area

    图  10   易发性结果ROC曲线

    Figure  10.   ROC curve of susceptibility results

    图  11   现场验证位置及照片(摄于2021年1月)

    Figure  11.   Field review verification photos (January 2021)

    表  1   卡门子湾滑坡破坏模式总结表

    Table  1   Summary table of failure mode of Kamenziwan landslide

    岸坡类型缓倾切向坡
    破坏模式视倾向顺层牵引式滑坡
    孕灾(六面体)
    结构面
    斜坡表面“上陡下缓前临空”的台阶状折线地形
    底部滑带左侧顺层,右侧切层
    前缘剪出口河流冲刷侵蚀形成临空面,
    剪出口高程在145~175 m
    后缘边界一组陡倾结构面切割,出露岩层切面
    左侧边界多组结构面相互切割形成阶梯状边界
    右侧边界岩层面
    边界特征两个约束边界(右、后)+两个自由边界(左、前)
    物质组成
    条件
    滑体块裂岩体
    滑带中后部由三组结构面及岩层面形成阶梯状
    滑带,前部岩层溃曲形成缓倾结构面
    滑床侏罗系中统沙溪庙组(J2s)上部灰绿色砂岩
    夹泥岩,下部紫红色泥岩夹砂岩
    下载: 导出CSV

    表  2   卡门子湾滑坡区易发性结果分析表

    Table  2   Summey table of landslide susceptibility results at Kamenziwan landslide area

    以疑似滑坡区为样本的滑坡易发性评价结果
    卡门子湾滑坡区栅格数各分区占比/%极高及较高易发区占比/%
    极高易发598626.682.8
    高易发1264056.2
    中易发375016.7
    低易发1130.5
    极低易发00.0
    总计22489100
    下载: 导出CSV
  • [1] 殷坤龙,朱良峰. 滑坡灾害空间区划及GIS应用研究[J]. 地学前缘,2001,8(2):279 − 284. [YIN Kunlong,ZHU Liangfeng. Landslide hazard zonation and application of GIS[J]. Earth Science Frontiers,2001,8(2):279 − 284. (in Chinese with English abstract) DOI: 10.3321/j.issn:1005-2321.2001.02.010
    [2] 黄波林,殷跃平,李滨,等. 三峡工程库区岩溶岸坡岩体劣化及其灾变效应[J]. 水文地质工程地质,2020,47(4):51 − 61. [HUANG Bolin,YIN Yueping,LI Bin,et al. Rock mass deterioration and its catastrophic effect of karst bank slope in the Three Gorges Project Reservoir area[J]. Hydrogeology & Engineering Geology,2020,47(4):51 − 61. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.202003055
    [3] 卫童瑶,殷跃平,高杨,等. 三峡库区巫山县塔坪H1滑坡变形机制[J]. 水文地质工程地质,2020,47(4):73 − 81. [WEI Tongyao,YIN Yueping,GAO Yang,et al. Deformation mechanism of the taping H1 landslide in Wushan County in the Three Gorges Reservoir area[J]. Hydrogeology & Engineering Geology,2020,47(4):73 − 81. (in Chinese with English abstract) DOI: 10.16030/j.cnki.issn.1000-3665.202003043
    [4] 郭子正,殷坤龙,付圣,等. 基于GIS与WOE-BP模型的滑坡易发性评价[J]. 地球科学,2019,44(12):4299 − 4312. [GUO Zizheng,YIN Kunlong,FU Sheng,et al. Evaluation of landslide susceptibility based on GIS and WOE-BP model[J]. Earth Science,2019,44(12):4299 − 4312. (in Chinese with English abstract)
    [5] 张钟远,邓明国,徐世光,等. 镇康县滑坡易发性评价模型对比研究[J]. 岩石力学与工程学报,2022,41(1):157 − 171. [ZHANG Zhongyuan,DENG Mingguo,XU Shiguang,et al. Comparison of landslide susceptibility assessment models in Zhenkang County,Yunnan Province,China[J]. Chinese Journal of Rock Mechanics and Engineering,2022,41(1):157 − 171. (in Chinese with English abstract)
    [6]

    MOHAMMADY M,POURGHASEMI H R,PRADHAN B. Landslide susceptibility mapping at Golestan Province,Iran:a comparison between frequency ratio,Dempster-Shafer,and weights-of-evidence models[J]. Journal of Asian Earth Sciences,2012,61:221 − 236. DOI: 10.1016/j.jseaes.2012.10.005

    [7] 王世宝,庄建琦,樊宏宇,等. 基于频率比与集成学习的滑坡易发性评价:以金沙江上游巴塘—德格河段为例[J]. 工程地质学报,2022,30(3):817 − 828. [WANG Shibao,ZHUANG Jianqi,FAN Hongyu,et al. Evaluation of landslide susceptibility based on frequency ratio and ensemble learning:Taking the Batang-Dege section in the upstream of Jinsha River as an example[J]. Journal of Engineering Geology,2022,30(3):817 − 828. (in Chinese with English abstract)
    [8]

    DU Juan,GLADE T,WOLDAI T,et al. Landslide susceptibility assessment based on an incomplete landslide inventory in the Jilong Valley,Tibet,Chinese Himalayas[J]. Engineering Geology,2020,270:105572. DOI: 10.1016/j.enggeo.2020.105572

    [9]

    STEGER S,BRENNING A,BELL R,et al. Exploring discrepancies between quantitative validation results and the geomorphic plausibility of statistical landslide susceptibility maps[J]. Geomorphology,2016,262:8 − 23. DOI: 10.1016/j.geomorph.2016.03.015

    [10]

    MEINHARDT M,FINK M,TÜNSCHEL H. Landslide susceptibility analysis in central Vietnam based on an incomplete landslide inventory:Comparison of a new method to calculate weighting factors by means of bivariate statistics[J]. Geomorphology,2015,234:80 − 97. DOI: 10.1016/j.geomorph.2014.12.042

    [11]

    ABEDINI M,TULABI S. Assessing LNRF,FR,and AHP models in landslide susceptibility mapping index:A comparative study of Nojian watershed in Lorestan Province,Iran[J]. Environmental Earth Sciences,2018,77(11):405. DOI: 10.1007/s12665-018-7524-1

    [12] 夏辉,殷坤龙,梁鑫,等. 基于SVM-ANN模型的滑坡易发性评价—以三峡库区巫山县为例[J]. 中国地质灾害与防治学报,2018,29(5):13 − 19. [XIA Hui,YIN Kunlong,LIANG Xin,et al. Landslide susceptibility assessment based on SVM-ANN Models:A case stualy for Wushan County in the Three Gorges Reservoir[J]. The Chinese Journal of Geological Hazard and Control,2018,29(5):13 − 19. (in Chinese with English abstract)
    [13] 殷坤龙,柳源. 滑坡灾害区划系统研究[J]. 中国地质灾害与防治学报,2000,11(4):28 − 32. [YIN Kunlong,LIU Yuan. Systematic studies on landslide hazard zonation[J]. The Chinese Journal of Geological Hazard and Control,2000,11(4):28 − 32. (in Chinese with English abstract) DOI: 10.3969/j.issn.1003-8035.2000.04.007
    [14] 周天伦,曾超,范晨,等. 基于快速聚类-信息量模型的汶川及周边两县滑坡易发性评价[J]. 中国地质灾害与防治学报,2021,32(5):137 − 150. [ZHOU Tianlun,ZENG Chao,FAN Chen,et al. Landslide susceptibility assessment based on K-means cluster information model in Wenchuan and two neighboring counties,China[J]. The Chinese Journal of Geological Hazard and Control,2021,32(5):137 − 150. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2021.05-17
    [15] 石菊松,张永双,董诚,等. 基于GIS技术的巴东新城区滑坡灾害危险性区划[J]. 地球学报,2005,26(3):275 − 282. [SHI Jusong,ZHANG Yongshuang,DONG Cheng,et al. GIS-based landslide hazard zonation of the new Badong County site[J]. Acta Geosicientia Sinica,2005,26(3):275 − 282. (in Chinese with English abstract) DOI: 10.3321/j.issn:1006-3021.2005.03.014
    [16] 郭子正,殷坤龙,黄发明,等. 基于滑坡分类和加权频率比模型的滑坡易发性评价[J]. 岩石力学与工程学报,2019,38(2):287 − 300. [GUO Zizheng,YIN Kunlong,HUANG Faming,et al. Evaluation of landslide susceptibility based on landslide classification and weighted frequency ratio model[J]. Chinese Journal of Rock Mechanics and Engineering,2019,38(2):287 − 300. (in Chinese with English abstract) DOI: 10.13722/j.cnki.jrme.2018.0838
    [17]

    MA Shuyue,QIU Haijun,HU Sheng,et al. Quantitative assessment of landslide susceptibility on the Loess Plateau in China[J]. Physical Geography,2020,41(6):489 − 516. DOI: 10.1080/02723646.2019.1674559

    [18] 罗路广,裴向军,黄润秋,等. GIS支持下CF与Logistic回归模型耦合的九寨沟景区滑坡易发性评价[J]. 工程地质学报,2021,29(2):526 − 535. [LUO Luguang,PEI Xiangjun,HUANG Runqiu,et al. Landslide susceptibility assessment in Jiuzhaigou scenic area with GIS based on certainty factor and Logistic regression model[J]. Journal of Engineering Geology,2021,29(2):526 − 535. (in Chinese with English abstract) DOI: 10.13544/j.cnki.jeg.2019-202
    [19] 罗路广,裴向军,崔圣华,等. 九寨沟地震滑坡易发性评价因子组合选取研究[J]. 岩石力学与工程学报,2021,40(11):2306 − 2319. [LUO Luguang,PEI Xiangjun,CUI Shenghua,et al. Combined selection of susceptibility assessment factors for Jiuzhaigou earthquake-induced landslides[J]. Chinese Journal of Rock Mechanics and Engineering,2021,40(11):2306 − 2319. (in Chinese with English abstract) DOI: 10.13722/j.cnki.jrme.2021.0198
    [20]

    HONG Haoyuan,CHEN Wei,XU Chong,et al. Rainfall-induced landslide susceptibility assessment at the Chongren area (China) using frequency ratio,certainty factor,and index of entropy[J]. Geocarto International,2017,32(2):139 − 154.

    [21] 刘璐瑶,高惠瑛. 基于证据权与Logistic回归模型耦合的滑坡易发性评价[J]. 工程地质学报,2021:1 − 11 [2021-05-13]. [LIU Luyao,GAO Huiying. landslide susceptibility assessment based on coupling of woe model and Logistic regression model[J]. Journal of Engineering Geology,2021:1 − 11 [2021-05-13]. (in Chinese with English abstract) DOI: 10.13544/j.cnki.jeg.2020-482
    [22] 李郎平,兰恒星,郭长宝,等. 基于改进频率比法的川藏铁路沿线及邻区地质灾害易发性分区评价[J]. 现代地质,2017,31(5):911 − 929. [LI Langping,LAN Hengxing,GUO Changbao,et al. Geohazard susceptibility assessment along the Sichuan-Tibet railway and its adjacent area using an improved frequency ratio method[J]. Geoscience,2017,31(5):911 − 929. (in Chinese with English abstract) DOI: 10.3969/j.issn.1000-8527.2017.05.004
    [23]

    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

    [24]

    ZHANG Yixing,LAN Hengxing,LI Langping,et al. Optimizing the frequency ratio method for landslide susceptibility assessment:A case study of the Caiyuan Basin in the southeast mountainous area of China[J]. Journal of Mountain Science,2020,17(2):340 − 357. DOI: 10.1007/s11629-019-5702-6

    [25]

    HE Keqiang,WANG Shangqing,DU Wen,et al. Dynamic features and effects of rainfall on landslides in the Three Gorges Reservoir region,China:Using the Xintan landslide and the large Huangya landslide as the examples[J]. Environmental Earth Sciences,2010,59(6):1267. DOI: 10.1007/s12665-009-0114-5

    [26]

    TANG Huiming,WASOWSKI J,JUANG C H. Geohazards in the three Gorges Reservoir Area,China - Lessons learned from decades of research[J]. Engineering Geology,2019,261:105267. DOI: 10.1016/j.enggeo.2019.105267

    [27]

    YIN Yueping,HUANG Bolin,ZHANG Quan,et al. Research on recently occurred reservoir-induced Kamenziwan rockslide in Three Gorges Reservoir,China[J]. Landslides,2020,17(8):1935 − 1949. DOI: 10.1007/s10346-020-01394-7

    [28]

    Regional level landslide inventory maps of the Shyok River watershed, Northern Pakistan[J]. Bulletin of Engineering Geology and the Environment, 2016, 75(2): 563 − 574.

    [29] 吴润泽,胡旭东,梅红波,等. 基于随机森林的滑坡空间易发性评价—以三峡库区湖北段为例[J]. 地球科学,2021,46(1):321 − 330. [WU Runze,HU Xudong,MEI Hongbo,et al. Spatial susceptibility assessment of landslides based on random forest:A case study from Hubei section in the Three Gorges Reservoir area[J]. Earth Science,2021,46(1):321 − 330. (in Chinese with English abstract)
    [30] 周超,殷坤龙,曹颖,等. 基于集成学习与径向基神经网络耦合模型的三峡库区滑坡易发性评价[J]. 地球科学,2020,45(6):1865 − 1876. [ZHOU Chao,YIN Kunlong,CAO Ying,et al. Landslide susceptibility assessment by applying the coupling method of radial basis neural network and adaboost:A case study from the Three Gorges Reservoir area[J]. Earth Science,2020,45(6):1865 − 1876. (in Chinese with English abstract)
    [31] 闫举生,谭建民. 基于不同因子分级法的滑坡易发性评价—以湖北远安县为例[J]. 中国地质灾害与防治学报,2019,30(1):52 − 60. [YAN Jusheng,TAN Jianmin. Landslide susceptibility assessment based on different factor classification methods:A case study in Yuanan County of Hubei Province[J]. The Chinese Journal of Geological Hazard and Control,2019,30(1):52 − 60. (in Chinese with English abstract)
    [32]

    LIU Zhongqiang,GILBERT G,CEPEDA J M,et al. Modelling of shallow landslides with machine learning algorithms[J]. Geoscience Frontiers,2021,12(1):385 − 393. DOI: 10.1016/j.gsf.2020.04.014

    [33]

    LAN H X,ZHOU C H,WANG L J,et al. Landslide hazard spatial analysis and prediction using GIS in the Xiaojiang watershed,Yunnan,China[J]. Engineering Geology,2004,76(1/2):109 − 128.

    [34] 周越,曾昭发,唐海燕,等. 公路勘察中滑坡体的地球物理特征与分析—以张榆线公路勘察为例[J]. 吉林大学学报(地球科学版),2021,51(2):638 − 644. [ZHOU Yue,ZENG Zhaofa,TANG Haiyan,et al. Geophysical characteristics of landslide body in highway reconnaissance: A case study in highway prospecting of Zhangyu line[J]. Journal of Jilin University (Earth Science Edition),2021,51(2):638 − 644. (in Chinese with English abstract)
  • 期刊类型引用(2)

    1. 徐开慧. 差异性边界约束边坡稳定性分析及治理措施研究. 江苏建筑. 2024(02): 112-115 . 百度学术
    2. 陆烨,黄越川. 基于DIC分析的含石量对碎石土边坡稳定性影响. 中国地质灾害与防治学报. 2023(04): 49-57 . 本站查看

    其他类型引用(0)

图(11)  /  表(2)
计量
  • 文章访问数:  5258
  • HTML全文浏览量:  2204
  • PDF下载量:  472
  • 被引次数: 2
出版历程
  • 收稿日期:  2021-12-27
  • 修回日期:  2022-04-05
  • 录用日期:  2022-04-07
  • 网络出版日期:  2023-02-19
  • 刊出日期:  2023-04-24

目录

/

返回文章
返回