ISSN 1009-6248CN 61-1149/P 双月刊

主管单位:中国地质调查局

主办单位:中国地质调查局西安地质调查中心
中国地质学会

    高级检索

    深度随机森林和随机森林算法的滑坡易发性评价对比以汉中市略阳县为例

    贾俊, 毛伊敏, 孟晓捷, 高波, 高满新, 武文英

    贾俊, 毛伊敏, 孟晓捷, 等. 深度随机森林和随机森林算法的滑坡易发性评价对比—以汉中市略阳县为例[J]. 西北地质, 2023, 56(3): 239-249. DOI: 10.12401/j.nwg.2023084
    引用本文: 贾俊, 毛伊敏, 孟晓捷, 等. 深度随机森林和随机森林算法的滑坡易发性评价对比—以汉中市略阳县为例[J]. 西北地质, 2023, 56(3): 239-249. DOI: 10.12401/j.nwg.2023084
    JIA Jun, MAO Yimin, MENG Xiaojie, et al. Comparison of Landslide Susceptibility Evaluation by Deep Random Forest and Random Forest Model: A Case Study of Lueyang County, Hanzhong City[J]. Northwestern Geology, 2023, 56(3): 239-249. DOI: 10.12401/j.nwg.2023084
    Citation: JIA Jun, MAO Yimin, MENG Xiaojie, et al. Comparison of Landslide Susceptibility Evaluation by Deep Random Forest and Random Forest Model: A Case Study of Lueyang County, Hanzhong City[J]. Northwestern Geology, 2023, 56(3): 239-249. DOI: 10.12401/j.nwg.2023084

    深度随机森林和随机森林算法的滑坡易发性评价对比—以汉中市略阳县为例

    基金项目: 中国地质调查局项目“西北典型地区地质灾害调查”(DD20221739),“黄土高原等典型地区地质灾害精细调查与风险管控”(DD20221739)联合资助。
    详细信息
      作者简介:

      贾俊(1985−),男,高级工程师,主要从事地质灾害调查工作。E−mail:geoj@qq.com

      通讯作者:

      孟晓捷(1986−),男,高级工程师,主要从事地质灾害调查工作。E−mail: 270405820@qq.com

    • 中图分类号: P694

    Comparison of Landslide Susceptibility Evaluation by Deep Random Forest and Random Forest Model: A Case Study of Lueyang County, Hanzhong City

    • 摘要:

      针对浅层的机器学习模型泛化能力低而导致其滑坡易发性评价模型预测精度不高的问题,笔者围绕陕西省汉中市略阳县城中心为研究区,采用深度随机森林构建区域地灾易发性评价模型来提升预测精度。依据略阳县滑坡成灾机理研究成果,选取坡度、相对高差、坡向、坡型、工程地质岩组、断裂距离、水系距离、公路铁路距离、植被覆盖等9个因子作为易发性评价指标;将研究区栅格单元按5 m × 5 m进行划分并提取评价因子值,输入深度随机森林评价模型,从而获得研究区易发性评价图。依据评价结果略阳县地质灾害可划分为极高易发区、高易发区、中易发区、低易发区4个等级,面积所占比例分别为5.31%、22.97%、42.11%、29.61%,其划分结果与研究区内地质灾害实际发育情况吻合,合理反映研究区地灾分布的总体特征。深度随机森林的地质灾害易发性预测模型在ROC曲线下面积值(AUC)为91.2%,高于随机森林预测模型的86.3%,表明该模型具有一定的合理性与可行性,可为区域滑坡易发性评价进一步提供新方法。

      Abstract:

      To address the problem of low prediction accuracy of landslide susceptibility evaluation model due to the difficulty of knowledge reuse and generalization of shallow machine learning model, this paper takes Lueyang County, Hanzhong City, Shaanxi Province as the study area and uses deep random forest to build a regional geological disaster susceptibility evaluation model to improve the prediction accuracy. Firstly, based on the research results of landslide mechanism in Lueyang County, nine factors such as slope, relative height difference, slope direction, slope type, engineering geological rock group, fault distance, river system distance, road and railroad distance, and vegetation cover are selected as susceptibility evaluation indexes; secondly, the study area is divided into 5 m × 5 m raster cells and the values of evaluation factors are extracted and input into the depth random forest evaluation model; finally, the susceptibility evaluation map of the study area is obtained. Based on the evaluation results, geological hazards in Lueyang County can be classified into four levels: very high susceptibility, high susceptibility, medium susceptibility, and low susceptibility, with the proportion of area being 5.31%, 22.97%, 42.11%, and 29.61%. The classification results are consistent with the actual development of geological hazards and reasonably reflect the overall characteristics of geological hazard distribution in the study area. In addition, the area under the ROC curve of the geological hazard susceptibility prediction model of deep random forest is 91.2%, which is higher than 86.3% of the random forest prediction model, indicating that the model is reasonable and feasible, and can provide new ideas for the evaluation of regional landslide susceptibility.

    • 昔格达地层(NQx)是出露于中国西南地区的典型晚新生代地层,以四川省盐边县红格乡昔格达村命名,主要分布在金沙江、大渡河、安宁河等干流和支流河谷、断裂谷以及侵蚀洼地中(蒋复初等, 1999; Kong et al., 2009; Deng et al., 2020; 施云云, 2020)(图1)。昔格达地层主要由黏土、粉砂和细砂等细粒组分构成,具有水平层理构造,是一种介于土体与岩体之间的特殊地质体,大多数学者将其归为半成岩(刘惠军等, 2004; 徐则民等, 2011; Du et al., 2020)。该套地层具有不良的工程地质特性,是工程地质问题和地质灾害的有利孕生载体(陈林等, 2015; Ding et al., 2017; Xue et al., 2018; 周平等, 2020)。例如,成昆铁路在隧道施工穿越昔格达地层时,隧道发生拱顶沉降、基地下沉、局部浸湿失稳、边坡滑塌等问题(Zhou et al., 2021);四川汉源中海村昔格达地层滑坡中断了S360省道,造成7人死亡2人失联和8栋民房被毁(徐奕梓等, 2022)。

      图  1  昔格达地层分布示意图
      Figure  1.  Distribution of Xigeda Formation

      目前,已有不少学者对中国四川渡口、攀枝花、西昌、冕宁、汉源、云南寨子村等西南地区的昔格达地层的物质组成、微观结构、物理力学特性等方面开展了大量研究,研究结果表明:①不同地区昔格达地层矿物组成与成分差异性较大,反映出不同地区昔格达地层沉积环境、物源及形成时代等存在差异(王思敬等, 1990; 李小泉, 1996; 孟庆会, 2011)。②昔格达地层黏土岩具有较强胶结结构特征,遇水强度陡然下降的原因之一是土体内黏土矿物结构被破坏(Du et al., 2020; Yang, 2021; Fu et al., 2022; 卢志鹏等, 2022)。③昔格达地层水理特性与其黏粒含量有着密切关系,黏粒含量越高,其液限越大,抗渗性越好(彭盛恩, 1986; 孟庆会, 2011)。④天然条件下昔格达地层具有良好的自稳性,受水影响显著,当含水率超过敏感界线时,强度陡然衰减(彭盛恩, 1986; 文丽娜等, 2005; 张威等, 2011; 安少鹏等, 2013; 周罕等, 2014; Du et al., 2020)。相比川南、攀西等地区的昔格达地层而言,关于川西泸定地区昔格达地层工程地质特性的研究鲜有报道。

      笔者以四川省泸定县海子坪昔格达地层黏土岩为研究对象,基于XRD测试、电镜扫描等测试,分析其物质组成和微观结构特征;通过现场渗透试验、室内力学试验,研究其水理特性与力学特性,并与国内其他地区昔格达地层对比分析,对深化川西地区昔格达地层研究与支撑当地防灾减灾具有一定参考价值。

      四川省泸定县海子坪地处大渡河右岸,共发育七级阶地,昔格达地层发育于T3阶地以上,出露厚度约440 m,为国内现今发现的昔格达地层沉积最厚的剖面(图2图3)。以昔格达地层为基座的T4~T7阶地,拔河高度分别为770 m、820 m、880 m、930 m。在大渡河河谷处发现T1~T3阶地,拔河高度分别为10 m、50 m、100 m,其中T4与T3阶地相对高差达670 m。

      图  2  泸定海子坪昔格达地层分布图(影像据Google Earth)
      Figure  2.  Distribution of Xigeda formation in Haiziping village, Luding county (Image from Google Earth)
      图  3  大渡河海子坪河流阶地剖面图
      Figure  3.  Terrace profile of Dadu River in Haiziping Village

      海子坪昔格达地层由黄色和灰色两种黏土岩互层产出(图4),该套地层底部发育一定厚度的角砾与碎块石层,与下伏基岩呈角度不整合接触。受构造运动影响,昔格达地层整体向西倾斜,倾角约为10°~20°,局部可达30°左右(蒋复初等, 1999)。受植被、风化剥蚀、人类工程活动等因素影响,海子坪昔格达地层露头条件连续性较差,仅在小海子、干海子等多处零星出露。现场调查发现,该套地层中发育较多的浅层滑坡和坡面泥流。在强降雨条件下,部分滑坡失稳后可能顺着沟道运动,并迅速转化为泥石流,对沟口造成一定危害(宋德光等, 2023)。

      图  4  昔格达地层出露及取样剖面
      a. 昔格达组地层取样剖面;b. 黄色土样;c. 灰色土样
      Figure  4.  Xigeda Formation developmental characteristics and sampling profile

      在昔格达地层的成因机制方面,目前已提出的主要有冰湖沉积、断陷湖泊沉积、堰塞湖沉积、泛湖期大湖沉积、山间盆地等(蒋复初等, 1999; 陈智梁等, 2004; 王书兵等, 2006; 王萍等, 2011; 施云云, 2020; 罗璐, 2021),泸定海子坪昔格达地层的成因机制主要有断陷盆地成因和堰塞湖成因两种认识(蒋复初等, 1999; 陈智梁等, 2004)。野外调查表明,大渡河流域泸定-得妥段大型-巨型古滑坡较为发育,且多具有历史堵江特征(吴俊峰等, 2011; Wu et al., 2019),因而认为泸定海子坪昔格达地层极有可能是古滑坡堰塞湖沉积物。在其形成年代学研究方面,不同学者采用古地磁、光释光、ESR、宇生核素等多种测年技术测得的年龄跨度较大(蒋复初等, 1999; 王书兵等, 2006; 王萍等, 2011; 施云云, 2020; 罗璐, 2021)(表1)。

      表  1  泸定海子坪昔格达地层测年方法及结果统计
      Table  1.  The dating methods and statistical results of the Xigeda formation in Haiziping village, Luding county
      测年方法测试年龄(Ma)数据来源
      古地磁4.20~2.60蒋复初等(1999)、
      王书兵等(2006
      光释光1.78~1.13罗璐(2021
      ESR0.70~0.40王萍等(2011
      宇生核素1.04~0.53施云云(2020
      下载: 导出CSV 
      | 显示表格

      在泸定海子坪T4~T6阶地出露的昔格达地层中分别采集黄色和灰色黏土岩样品(图4),并开展X衍射测试(XRD)和电镜扫描(SEM)等测试,分析其矿物成分和微观结构。结果表明,两种土样的矿物成分相同,均以石英和黏土矿物为主,含少量斜长石、方解石、钾长石、菱铁矿和角闪石等矿物(图5图6)。各矿物占比明显不同,黄色黏土岩中的黏土矿物含量比灰色黏土岩高12%,方解石含量少10%,表明黄色黏土岩沉积时期风化作用更强烈,处于暖湿气候环境(罗运利等, 1998)。此外,黄色和灰色黏土岩中黏土矿物成分均以伊利石为主,占比分别高达93%和89%,蒙脱石和绿泥石含量较少(表2)。

      图  5  黄色黏土岩XRD测试分析
      a. XRD分析谱图;b. 矿物成分及含量
      Figure  5.  XRD analysis of yellow claystone sample
      图  6  灰色黏土岩XRD测试分析
      a. XRD分析谱图;b. 矿物成分及含量
      Figure  6.  XRD analysis of gray claystone sample
      表  2  黄色和灰色黏土岩黏土矿物分析结果
      Table  2.  Test results of clay minerals of yellow and gray claystone samples
      序号样品类型黏土矿物检测结果(%)
      蒙脱石
      S
      伊蒙混层
      I/S
      伊利石
      It
      高岭石
      K
      绿泥石
      C
      I/S混层比(%S)
      1黄色黏土岩3934
      2灰色黏土岩5896
      下载: 导出CSV 
      | 显示表格

      根据黄、灰色黏土岩天然状态下的SEM图像可知,黄色黏土岩结构较致密,成分主要为碎屑岩颗粒和填充物。碎屑岩颗粒以石英为主,粒径10~25 μm,沿碎屑岩颗粒边缘分布自生黏土矿物集粒。填充物以伊利石和胶结物为主,含少量绿泥石。伊利石呈片状镶嵌于碎屑粒间孔隙,粒径3~6 μm,绿泥石呈叶片状附着于碎屑岩颗粒边缘。胶结物为黏土矿物黏结,主要表现为填充于粒间孔隙的绒状、针状结构(图7)。灰色黏土岩结构松散,孔隙和微小裂隙较发育,成分主要为粉砂粒碎屑和片状黏土矿物,粉砂粒碎屑以石英为主,粒径5~15 μm。片状黏土矿物主要为伊利石,粒径5~8 μm,黏土矿物颗粒间胶结作用较弱,主要呈片状填充于粒间孔隙(图8)。

      图  7  黄色黏土岩SEM图像
      a. 1 500倍;b. 3 000倍
      Figure  7.  SEM image of yellow claystone sample
      图  8  灰色黏土岩SEM图像
      a. 1 500倍;b. 3 000倍
      Figure  8.  SEM image of gray claystone sample

      基于Image-Pro Plus(IPP)图像分析软件,对黄、灰色黏土岩的微观结构进行定量研究。考虑到SEM图像超过一定放大倍数后存在“失真”现象(梁坤, 2021),选取放大1500倍的SEM图像进行分析(图9图10)。根据测量等效孔径(d)范围,将土样的微观孔隙分成4种孔隙类型:大型孔隙(d>4 μm);中型孔隙(1 μm<d<4 μm);小型孔隙(0.4 μm<d<1 μm);微型孔隙(d<0.4 μm)。结果表明,黄、灰色土样中微观孔隙随机独立分布,连通率较低,灰色黏土岩微观孔隙更为发育(图9b图10b)。黄、灰色黏土岩微观孔隙均以微、小型孔隙为主,中型孔隙次之,大型孔隙并不发育(表3),反映出泸定海子坪昔格达地层静水沉积环境良好,后期受环境动力的扰动影响较小。

      图  9  黄色黏土岩SEM图像和二值图像
      a. SEM图像(1 500倍);b. 二值化图像(白色为孔隙)
      Figure  9.  SEM image and binary image of yellow claystone sample
      图  10  灰色黏土岩SEM图像和二值图像
      a. SEM图像(1 500倍);b. 二值化图像(白色为孔隙)
      Figure  10.  SEM image and binary image of gray claystone sample
      表  3  黄色和灰色黏土岩各孔径级别对应的微观孔隙测量结果
      Table  3.  Measurement results of different pore sizes of yellow and gray claystone samples
      孔径级别(µm)黄色黏土岩灰色黏土岩
      数量 N(个)孔径d(µm)面积A(µm2丰度(C)数量 N(个)孔径d(µm)面积A(µm2丰度(C
      <0.43560.230.040.535120.240.050.53
      0.4~11370.590.240.461780.570.260.51
      1~4221.51.640.45301.241.060.44
      >4000014.115.990.52
       注:CBL的比值,B表示孔隙的短轴长度,L表示孔隙长轴长度,A为平均孔隙面积。
      下载: 导出CSV 
      | 显示表格

      黄、灰色黏土岩的微观结构孔隙丰度(C)分布范围具有较高相似度,孔隙丰度(C)在各个区间都有分布,但主要集中在0.2~0.5与0.8~0.9这两个区间内(图11)。可见泸定海子坪昔格达地层黄、灰色黏土岩的微观孔隙结构均以扁圆状居多,等轴型次之,长条型最少。

      图  11  黄色和灰色黏土岩微观孔隙丰度分布比例
      Figure  11.  Distribution ratio of microscopic porosity of yellow and gray claystone samples

      根据《土工试验方法标准》(GBT 50123-2019),对黄、灰色黏土岩进行了一系列岩土体物理性质测试。结果表明,黄、灰色黏土岩基本物性指标较为接近(表4),其中黄色黏土岩塑性指数为7.8~14.2,灰色黏土岩塑性指数为8.4~12.9,均属粉土–亚黏土。

      表  4  黄色和灰色黏土岩基本物性指标
      Table  4.  Basic physical property indexes of yellow and gray claystone samples
      土样
      编号
      样品
      类型
      含水率
      (%)
      土粒
      比重
      天然密度
      (g/cm3
      干密度
      (g/cm3
      孔隙比液限
      (%)
      塑限
      (%)
      塑性
      指数
      饱和含水
      率(%)
      XGD-01灰色黏土岩3.72.691.861.760.56936.425.910.530.93
      XGD-02黄色黏土岩4.52.701.801.730.55036.324.911.431.91
      下载: 导出CSV 
      | 显示表格

      依照《水利水电工程注水试验规程》(SL 345-2007),在海子坪上多处昔格达地层出露点开展试坑原位双环注水渗透试验(图2),测得的渗透系数为3.62×10−4~7.34×10−4 cm/s,平均渗透系数为5.11×10−4 cm/s,表明该处昔格达黏土岩为中等透水介质(表5)。本次试验结果与其他地区昔格达地层黏土岩渗透特性相比偏大(李小泉, 1996; 宋为广等, 2017)。现场调查发现,泸定海子坪昔格达地层完整性好、无微裂纹或裂缝土块的渗透性较差,但土体内节理裂隙发育,且受扰动后极易形成大量的微小裂隙,导致现场测得的渗透系数偏大。

      表  5  双环注水试验结果
      Table  5.  Results of double-ring water injection tests
      试验点序号渗透系数
      (cm/s)
      平均渗透系数
      (cm/s)
      透水性等级
      S013.62×10−45.11×10−4中等透水
      S027.34×10−4
      S034.36×10−4
      下载: 导出CSV 
      | 显示表格

      现场调查发现,泸定海子坪昔格达地层黄、灰色黏土岩成岩程度较低,取样过程中对土体扰动较大,本研究采用重塑样开展了黄、灰色黏土岩在不同含水率条件下的三轴固结不排水剪切试验。结果表明:在相同围压条件下,土样的峰值强度均随含水率的增加而减小,在试验含水率范围内,土样峰值强度与含水率呈较好的线性关系(图12)。

      图  12  不同围压下峰值强度与含水率关系图
      a. 黄色黏土岩;b. 灰色黏土岩
      Figure  12.  Relationship between peak strength and moisture content under different confining pressures

      两种土样的黏聚力均随含水率的增加而降低,且含水率越高,降幅越大(图13)。当黄色黏土岩含水率由20%增至饱和含水率30.93%时,黏聚力从121.12 kPa降至18.54 kPa,降幅为84.7%;灰色黏土岩含水率由15%增至饱和含水率31.91%时,黏聚力从50.98 kPa降至0.4 kPa,降幅为99.2%。此外,黄色黏土岩黏聚力远高于灰色黏土岩,这也进一步说明了黄色黏土岩黏土矿物间胶结作用更强,结构更为致密;而随着含水率不断增加,黏土矿物间胶结结构被破坏,黏聚力显著降低。

      图  13  黏聚力、内摩擦角与含水率关系
      a. 黄色黏土岩;b. 灰色黏土岩
      Figure  13.  Relationship between cohesion, internal friction angle and moisture content

      当含水率增加时,黄色黏土岩内摩擦角逐渐增大,含水率由20%增至饱和含水率30.93%时,内摩擦角由22.62°增大至29.12°,增幅为28.7%。而灰色黏土岩内摩擦角表现为先降低后增大,但变化幅度不大,不同含水率土样之间内摩擦角的差值在±4°范围内波动。

      根据已有研究(王思敬等, 1990; 孟庆会, 2011)和本次试验结果可知,泸定海子坪昔格达地层黏土岩界限含水率总体低于攀西地区,而塑性指数均远小于攀西地区(图14)。由于土体界限含水率(液限)与土体中的黏土矿物含量呈正相关,推断泸定海子坪昔格达地层黏土岩中的黏土矿物含量比攀西地区的低。

      图  14  不同地区昔格达地层界限含水率对比图
      a. 塑限;b. 液限;c. 塑性指数
      Figure  14.  Correlation of moisture content of Xigeda Formation boundary in different regions

      对比分析不同地区昔格达地层的渗透特性可知,泸定海子坪昔格达地层黏土岩渗透系数的数量级与其他地区粉砂岩类相当(均为10−4cm/s),而明显大于其他地区黏土岩类的渗透系数(主要为10−5~10−7cm/s)(表6)。由此可知,泸定海子坪昔格达地层黏土岩的渗透特性与其他地区昔格达地层的粉砂岩类大致相当,而大于黏土岩类。昔格达地层在不同地区与不同岩性上的渗透特性差异较大,这可能与试验手段、区域地质条件、成岩条件及后期扰动程度等因素密切相关。

      表  6  不同地区昔格达地层渗透特性比较
      Table  6.  Comparison of permeability characteristics of Xigeda Formation in different regions
      地区试验对象试验条件渗透系数(cm/s)透水性等级数据来源
      泸定海子坪粉土-亚黏土现场双环注水3.62×10−4~7.34×10−4中等透水本研究
      攀枝花格里坪亚黏土室内渗透3.7×10−5~7.2×10−4弱透水宋为广等(2017
      粉砂岩室内渗透2.0×10−5~2.78×10−5弱透水左永振等(2016
      攀枝花粟子坪亚黏土室内渗透2.8×10−7~3.3×10−7极弱透水李小泉(1996
      云南龙开口镇黏土岩室内渗透10−8~10−7极弱透水张德强等(2021
      川南地区泥岩现场、室内渗透1.7×10−5~5.25×10−5弱透水钟成等(2012
      砂岩1.6×10−4中等透水
      西昌经久乡黏土岩室内渗透、
      钻孔压水
      1.95×10−7~4.61×10−5弱透水杨碧(2010
      粉砂岩1.47×10−4~6.54×10−4中等透水
      下载: 导出CSV 
      | 显示表格

      对比分析不同地区昔格达地层力学强度指标(张文举, 2003; 黄绍槟等, 2005; 文丽娜等, 2005; 安少鹏等, 2013; 周罕等, 2014; Du et al., 2020)可知,泸定海子坪粉土-亚黏土类黏聚力低于云南寨子村半成岩和攀西地区黏土岩类,而高于攀西地区粉砂岩类;内摩擦角略低于攀枝花新九乡砂岩类,而高于其他地区半成岩、黏土岩及粉砂岩类(图15)。

      图  15  不同地区昔格达地层强度指标与含水率的关系图
      a. 黏聚力;b. 内摩擦角
      Figure  15.  Relationship between strength index and moisture content of Xigeda Formation in different regions

      不同地区的昔格达地层黏聚力均随含水率的增加而降低,但对含水率的敏感性有所不同。随着含水率的增加,泸定海子坪粉土-亚黏土黏聚力的水敏特性最为显著,云南寨子村半成岩次之,攀西地区砂岩类水敏特性最弱(图15a)。此外,其他地区昔格达地层内摩擦角均随含水率的增加而降低,而泸定地区粉土-亚黏土内摩擦角随含水率的增加而增大。随着含水率的增加,泸定海子坪粉土-亚黏土内摩擦角增幅越大,云南寨子村半成岩及攀西地区砂岩类内摩擦角降幅最大,攀西地区黏土岩类内摩擦角降幅最小(图15b)。由上分析可知,不同地区昔格达地层对含水率敏感性存在区域性及岩性上的差异;其中泸定海子坪地区粉土-亚黏土力学特性的水敏特性最为显著,攀西地区黏土岩及云南寨子村半成岩次之,攀西地区砂岩类最小。

      中国西南地区的昔格达地层成岩程度较低,属于一类特殊的硬土软岩,其工程地质特性主要表现为:①昔格达地层的物理力学特性与其黏粒含量有着密切关系,其物理力学性质介于土体和岩体之间。②不同地区昔格达地层水理特性差异较大,这与不同地区昔格达地层成岩程度、赋存环境、物质组成、风化程度及后期受扰动等因素密切相关。③水对昔格达地层强度特性具有显著的弱化效应,这归因于较高的含水率破坏了土体内部黏粒的胶结结构,不同地区昔格达地层的水敏特性有所区别,这与其黏粒含量差异性有关。由此可知,昔格达地层具有良好的孕灾条件,其内部微裂隙较为发育,后期受地震、降雨和工程扰动等因素影响,这些节理裂隙逐渐演化形成地表水入渗的优势通道和控制岩土体稳定性的软弱结构面,易形成崩塌、滑坡和泥石流地质灾害,甚至形成灾害链(铁永波等, 2022)。

      针对昔格达地层遇水软化、易滑致灾的特点,在进行人类工程活动时,须引起足够的重视,针对其控灾因子与致灾模式采取相应的防范措施(冯立等, 2014; 李绵绵等, 2019),如对公路开挖形成的高陡边坡,建议采取挡土墙、排水沟等工程措施:对极端工况诱发下可能形成高位滑坡泥石流的斜坡区段,建议开展详细工程地质调查,对不稳定斜坡开展综合监测,必要时可进行工程加固,设置截排水沟,并在下游沟道设置排导槽与谷坊坝,避免形成高位启动-沟道泥石流的灾害链。

      (1)泸定海子坪昔格达地层黏土岩主要成分为粉细砂与黏土等细粒物质,由黄色和灰色两种黏土岩互层产出,物质组成均以石英与黏土矿物为主,黏土矿物均以伊利石为主。黄色黏土岩的黏土矿物间胶结结构特征明显,而灰色黏土岩的黏土矿物呈薄片散体状,前者比后者结构致密、黏土矿物间胶结作用强。

      (2)黄、灰色黏土岩黏聚力均随含水率的增大而逐渐降低,内摩擦角随含水率的增加而增大。水与黏土矿物作用破坏了黏土矿物间胶结结构,导致黏聚力逐渐丧失,而岩碎屑颗粒相互作用导致内摩擦角增大,且在高含水率条件下增幅更大。

      (3)泸定海子坪昔格达地层黏土岩成岩程度较低,节理裂隙发育,扰动后极易开裂,平均渗透系数为5.11×10−4 cm/s,属中等透水介质,其渗透特性介于攀西地区黏土岩类–砂岩类之间。

      (4)不同地区昔格达地层力学特性均随含水率增加而逐渐减小,其水敏特性存在区域与岩性上的差异,泸定昔格达地层黏土岩遇水强度弱化特性最为显著。泸定海子坪昔格达地层黏土岩黏聚力介于其他地区昔格达地层砂岩类–半成岩及黏土岩之间,而内摩擦角与砂岩类相当。

    • 图  1   研究区位置图

      Figure  1.   Location of researching area

      图  2   研究区灾害示意图

      a. 略阳县谢家坪村阴坡里滑坡;b. 略阳县碾子湾村何硖路旁崩塌

      Figure  2.   Hazard photographs of researching area

      图  3   深度随机森林模型流程结构(据Zhou, 2017

      Figure  3.   Process structure illustration of Deep Random Forest Model

      图  4   滑坡与单评价因子分类图

      Figure  4.   The landslides division and evaluation factors classification

      图  5   略阳县随机森林滑坡易发性评价图

      Figure  5.   Landslide susceptibility assessment map of Lueyang County base on Random Forest Model

      图  6   略阳县深度随机森林滑坡易发性评价图

      Figure  6.   Landslide susceptibility assessment map of Lueyang County base on Deep Random Forest Model

      图  7   深度随机森林与随机森林易发性评价ROC曲线对比图

      Figure  7.   DRF and RF Mode ROC Curve Comparison on Landslide susceptibility assessment

      表  1   各评价因子信息量表

      Table  1   Weighted information values of individual evaluation factors

      因子分级Si(km2Ni(个)Ii
      坡度 <10° 10.535625 8 −0.085244
      10°~20° 13.3961 21 0.639636
      20°~30° 35.77955 46 0.441342
      30°~40° 54.226925 38 −0.165514
      40°~50° 38.6219 12 −0.978836
      >50° 13.119575 12 0.100879
      相对高差 16~175 m 24.930925 53 0.944260
      175~238 m 55.38815 46 0.004353
      238~300 m 49.224425 29 −0.339017
      300~379 m 26.43665 8 −1.005232
      379~605 m 9.6986 1 −2.081904
      坡向 0~45° 23.746850 21 0.067150
      45°~90° 21.137625 23 0.274517
      90°~135° 17.924900 10 −0.393528
      135°~180° 18.869250 14 −0.108399
      180°~225° 20.641625 26 0.420864
      225°~275° 22.018400 18 −0.011429
      275°~315° 20.152475 13 −0.248300
      315°~360° 21.185550 12 −0.378335
      曲率 <−0.5(凹形坡) 11.1099 12 0.267147
      −0.5~0.5
      (直线形坡)
      143.597 119 0.002190
      >0.5(凸形坡) 10.969775 6 −0.413307
      工程地质岩组 坚硬岩组 65.882425 40 −0.308915
      半坚硬岩组 29.980725 13 −0.645528
      软硬相间岩组 58.766425 46 −0.054852
      松散岩组 11.050425 39 1.451170
      距断裂距离 <100 m 31.79 26 −0.010915
      100~200 m 25.87 26 0.195167
      200~500 m 49.79 53 0.252595
      500~1000 m 34.60 23 −0.218224
      >1000 m 23.64 9 −0.775554
      距河流水系
      距离
      <200 m 69.69 86 0.400425
      200~400 m 52.38 31 −0.334384
      400~600 m 27.98 17 −0.308056
      600~800 m 11.37 1 −2.240549
      >800 m 4.28 2 −0.569794
      距公路、铁路
      距离
      <100 m 18.99 58 1.306584
      100~500 m 46.12 41 0.072488
      500~1000 m 38.06 14 −0.809915
      1000~1500 m 26.29 11 −0.681226
      >1500 m 36.23 13 −0.834794
      NDVI −0.41~0.07 1.88 0
      0.07~0.32 7.64 26 1.415100
      0.32~0.52 26.25 56 0.947905
      0.52~0.68 91.74 50 −0.416840
      0.68~0.84 38.18 5 −1.842838
      下载: 导出CSV

      表  2   易发性等级划分与滑坡实际发生比率对比表

      Table  2   Comparison of susceptibility classification and actual landslide occurrence rate

      评价
      方法
      易发性
      等级
      a
      易发分区面积占比(%)
      b
      分区滑坡数量
      c
      滑坡百分比(%)
      c/a
      滑坡发生比率
      随机森林16.3410.730.04
      69.217454.010.78
      11.244029.202.60
      极高3.212216.065.00
      深度
      随机森林
      29.6175.110.17
      42.114129.930.71
      22.976245.261.97
      极高5.312719.713.71
      下载: 导出CSV
    • 黄发明, 殷坤龙, 蒋水华, 等. 基于聚类分析和支持向量机的滑坡易发性评价[J]. 岩石力学与工程地质学报, 2018, 37(1): 156-167

      HUANG Faming, YIN Kunlong, JIANG Shuihua, et al. Landslide susceptibility assessment based on clustering analysis and support vector machine[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(1): 156-167.

      李亭, 田原, 邬伦, 等. 基于随机森林方法的滑坡灾害危险性区划[J]. 地理与地理信息科学, 2014, 30(6): 25-30 doi: 10.3969/j.issn.1672-0504.2014.06.006

      LI Ting, TIAN Yuan, WU Lun, et al. Landslide susceptibility Mapping Using Random Forest[J]. Geography and Geo-Information Science, 2014, 30(6): 25-30. doi: 10.3969/j.issn.1672-0504.2014.06.006

      刘坚, 李树林, 陈涛, 等. 基于优化随机森林模型的滑坡易发性评价[J]. 武汉大学学报·信息科学版, 2018, 43(7): 1085-1091

      LIU Jian, LI Shulin, CHEN Tao, et al. Landslide Susceptibility Assesment Based on Optimized Random Forest Model[J]. Geomatics and Information Science of Wuhan University, 2018, 43(7): 1085-1091.

      孟庆华. 秦岭山区地质灾害风险评估方法研究—以陕西凤县为例[D]. 北京: 中国地质科学院, 2011

      MENG Qinghua. Study on the Methods of Geo-hazards Risk Assessment in Qinling Mountain: A case studyof Feng County, Baoji City, Shaanxi Province[D]. Beijing: Chinese Academy of Geological Sciences, 2011.

      孟晓捷, 张新社, 曾庆铭, 等. 基于加权信息量法的黄土滑坡易发性评价——以1: 5万天水市麦积幅为例[J]. 西北地质, 2022, 55(2): 249-259

      MENG Xiaojie, ZHANG Xinshe, ZENG Qingming, et al. The Susceptibility Evaluation of Loess Landslide Based on Weighted Information Value Method—Taking 1: 50000 Map of Maiji District of Tianshui City As an Example[J]. Northwestern Geology, 2022, 55(2): 249-259.

      邱海军. 区域滑坡崩塌地质灾害特征分析及其易发性和危险性评价研究——以宁强县为例[D]. 西安: 西北大学, 2012: 83

      QIU Haijun. Study on the Regional Landslide Characteristic Analysis and Hazard Assessment: A case study of Ningqiang County[D]. Xi’an: Northwest University, 2012: 83.

      邱维蓉, 吴帮玉, 潘学树, 等. 几种聚类优化的机器学习方法在灵台县滑坡易发性评价中的应用[J]. 西北地质, 2020, 53(1): 222-233

      QIU Weirong, WU Bangyu, PAN Xueshu, et al. Application of Several Cluster-optimization-based Machine Learning Methods in Evaluation of Landslide Susceptibility in Lingtai County[J]. Northwestern Geology, 2020, 53(1): 222-233.

      田春山, 刘希林, 汪佳, 等. 基于CF和Logistic回归模型的广东省地质灾害易发性评价[J]. 水文地质工程地质, 2016, 43(6): 154-170

      TIAN Chunshan, LIU Xilin, WANG Jia. et al. Geohazard susceptibility assessmentbased on CF model and Logistic Regression models in Guangdong[J]. Hydrogeology & Engineering Geology, 2016, 43(6): 154-170.

      王佳运, 毕俊擘, 杨旭东, 等. 山西吕梁山区城镇边坡风险分级与优化[J]. 地质通报, 2020, 39(12): 2004-2012 doi: 10.12097/j.issn.1671-2552.2020.12.014

      WANG Jiayun, BI Junbo, YANG Xudong, et al. Risk classification and optimization of to-wn side slope in Lüliang Mountain[J]. Geological Bulletin of China, 2020, 39(12): 2004-2012. doi: 10.12097/j.issn.1671-2552.2020.12.014

      王杰, 马凤山, 郭捷, 等. 一种改进的区域滑坡危险性评价模型及其应用[J]. 中国地质灾害与防治学报, 2011, 22(2): 14-19 doi: 10.3969/j.issn.1003-8035.2011.02.003

      WANG Jie, MA Fengshan, GUO Jie, et al. An improved regional landslide hazard assessment, model and its application[J]. The Chinese Journal of Geological Hazard and Control, 2011, 22(2): 14-19. doi: 10.3969/j.issn.1003-8035.2011.02.003

      王涛, 吴树仁, 石菊松, 等. 秦岭中部太白县地质灾害发育特征及危险性评估[J]. 地质通报, 2013, 32(12): 1976-1983

      WANG Tao, WU Shuren, SHI Jusong, et al. Case Study Of Landslide Characteristics And Hazard Assessment In Taibai County, Central Qinling Mountains[J]. Geological Bulletinof China, 2013, 32(12): 1976-1983.

      吴常润, 赵冬梅, 刘澄静, 等. 基于GIS和信息量模型的陇川县滑坡易发性评价[J]. 西北地质, 2020, 53(2): 308-320

      WU Changrun, ZHAO Dongmei, LIU Chengjing, et al. Landslide Susceptibility Assessment of Longchuan County Based on GIS and Information Value Model[J]. NorthwesternGeology, 2020, 53(2): 308-320.

      吴润泽, 胡旭东, 梅红波, 等. 基于随机森林的滑坡空间易发性评价: 以三峡库区湖北段为例[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.

      吴孝情, 赖成光, 陈晓宏, 等. 基于随机森林权重的滑坡危险性评价: 以东江流域为例[J]. 自然灾害学报, 2017, 26(5): 119-129

      WU Xiaoqing, LAI Chengguang, CHEN Xiaohong, et al. A landslide hazard assessment based on random forestweight: a case study in the Dongjiang River Basin[J]. Journal of Natural Disarsters, 2017, 26(5): 119-129.

      向喜琼. 区域滑坡地质灾害危险性评价与风险管理[D]. 成都: 成都理工大学, 2005: 22

      XIANG Xiqiong. Regional Landslide Hazard Assessment and Risk Management[D]. Chengdu: Chengdu University of Technology, 2005: 22.

      薛强, 张茂省, 李林等. 基于斜坡单元与信息量法结合的宝塔区黄土滑坡易发性评价[J]. 地质通报, 2015, 34(11): 2108-2115

      XUE Qiang, ZHANG Maosheng, LI Lin. Loess landslide susceptibility evaluation based on slope unit and information value method in Baota District, Yan’an. Geological Bulletin of China, 2015, 34(11): 2108-2115.

      张春山, 韩金良, 孙炜锋, 等. 陕西陇县地质灾害危险性分区评价[J]. 地质通报, 2008, 27(11): 1795-1801 doi: 10.3969/j.issn.1671-2552.2008.11.006

      ZHANG Chunshan, HAN Jinliang, SUN Weifeng, et al. Assessments of geohazard danger zoning in Longxian County, Shaanxi, China[J]. Geological Bulletin of China, 2008, 27(11): 1795-1801. doi: 10.3969/j.issn.1671-2552.2008.11.006

      张茂省, 贾俊, 王毅, 等. 基于人工智能(AI)的地质灾害防控体系建设[J]. 西北地质, 2019a, 52(2): 103-116

      ZHANG Maosheng, JIA Jun, WANG Yi, et. al. Construction of Geological Disarster Prevention and Control System Based on AI[J]. Northwestern Geology, 2019a, 52(2): 103-116.

      张茂省, 薛强, 贾俊, 等. 山区城镇地质灾害调查与风险评价方法及实践[J]. 西北地质, 2019b, 52(2): 125-135 doi: 10.19751/j.cnki.61-1149/p.2019.02.013

      ZHANG Maosheng, XUE Qiang, JIA Jun, et. al. Methods and Practice for the Investigationand Risk Assessment of Geo-hazards in Mountains Towns[J]. Northwestern Geology, 2019b, 52(2): 125-135. doi: 10.19751/j.cnki.61-1149/p.2019.02.013

      张向营, 张春山, 孟华君, 等. 基于Random Forest和AHP的贵德县北部山区滑坡危险性评价[J]. 水文地质工程地质, 2018, 45(4): 142-149

      ZHANG Xiangying, ZHANG Chunshan, MENG Huajun, et al. Landslide hazard evaluation in the northern mountainous area of Guide County based on Random Forest and AHP[J]. Hydrogeology & Engineering Geology, 2018, 45(4): 142-149.

      张永双, 苏生瑞, 吴树仁, 等. 强震区断裂活动与大型滑坡关系研究[J]. 岩石力学与工程学报, 2011, 30(增刊2): 3503-3513

      ZHANG Yongshuang, SU Shengrui, WU Shuren, et al. Research On Relationship Between Fault Movement And Large-Scale Landslide In Intensive Earthquake Region[J]. Chinese Journalof Rock Mechanics and Engineering, 2011, 30(Sup. 2): 3503-3513.

      赵建华, 陈汉林, 杨树峰, 等. 基于决策树算法的滑坡危险性区划评价[J]. 浙江大学学报(理学版), 2004, 31(4): 465-470

      ZHAO Jianhua, CHEN Hanlin, YANG Shufeng, et al. Landslide risk assessment based on decision tree arithmetic[J]. Journal of Zhejiang University (Science Edition), 2004, 31(4): 465-470.

      周静静, 赵法锁, 李辉, 等. 陕西省地质灾害与影响因素相关性研究[J]. 灾害学, 2019, 34(2): 228-234 doi: 10.3969/j.issn.1000-811X.2019.02.041

      ZHOU Jingjing, ZHAO Fasuo, LI Hui, et al. Correlational Research Between Geological Hazards and impact Factors in Shaanxi Province[J]. Journal of Catastrophology, 2019, 34(2): 228-234. doi: 10.3969/j.issn.1000-811X.2019.02.041

      周样样. 陕南地区强降雨条件下突发型地质灾害成因机制研究[D]. 西安: 长安大学, 2013

      ZHOU Yangyang. Study on Formation Mechanism of Abrupt Geological Hazards of Southern Shaanxi Region in Condition of Strong Rainstorm[D]. Xi’an: Chang’an University, 2013.

      Chung C J F, Fabbri A G. Validation of spatial prediction modelsfor landslidehazard mapping[J]. Natural Hazards, 2003, 30(3): 451-472. doi: 10.1023/B:NHAZ.0000007172.62651.2b

      Dou J , Yunus A P , Merghadi A , et al. Different sampling strategies for predicting landslide susceptibilities are deemed less consequential with deep learning[J]. Science of the Total Environment, 2020, 720: 137320. https://doi.org/10.1016/j.scitotenv.2020.137320

      Maher Ibrahim Samee, Biswajeet Pradhana, Saro Lee, et al. Application of convolutional neural networks featuring Bayesian optimization for landslide susceptibility assessment[J]. Catena, 2020, 186: 10424. https://doi.org/10.1016/j.catena.2019.104249

      Saro, Lee, Joo Hyung Ryu, Joong Sun Won, et al. Determination and application of the weights for landslide susceptibility mapping using an artificial neural network[J]. Engineering Geology, 2004, 71(3-4): 289-302. doi: 10.1016/S0013-7952(03)00142-X

      Snoek J, Larochelle H, Adams R P. et al. Practical Bayesian Optimization of Machine Learning Algorithms[J]. Advances in Neural Information Processing Systems, 2012, 4: 1–12.

      Fang Zhice, Wang Yi , Ling Peng, et al. Integration of convolutional neural network and conventional machine learning classifiers for landslide susceptibility mapping[J]. Computer & Geosciencees, 2020, 139: 104470. https://doi.org/10.1016/j.cageo.2020.104470

      Zhou Zhihua, Ji Feng. Deep Forest: Towards An Alternative to Deep Neural Networks[J]. Statistics, 2017, 71(3–4): 289–302.

    图(7)  /  表(2)
    计量
    • 文章访问数:  113
    • HTML全文浏览量:  23
    • PDF下载量:  33
    • 被引次数: 0
    出版历程
    • 收稿日期:  2023-03-14
    • 修回日期:  2023-04-22
    • 网络出版日期:  2023-05-28
    • 刊出日期:  2023-06-19

    目录

    /

    返回文章
    返回