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

基于自组织特征映射网络-随机森林模型的滑坡易发性评价以江西大余县为例

何书, 鲜木斯艳·阿布迪克依木, 胡萌, 陈康

何书,鲜木斯艳·阿布迪克依木,胡萌,等. 基于自组织特征映射网络-随机森林模型的滑坡易发性评价−以江西大余县为例[J]. 中国地质灾害与防治学报,2022,33(1): 132-140. DOI: 10.16031/j.cnki.issn.1003-8035.2022.01-16
引用本文: 何书,鲜木斯艳·阿布迪克依木,胡萌,等. 基于自组织特征映射网络-随机森林模型的滑坡易发性评价−以江西大余县为例[J]. 中国地质灾害与防治学报,2022,33(1): 132-140. DOI: 10.16031/j.cnki.issn.1003-8035.2022.01-16
HE Shu, ABUDIKEYIMU XMSY, HU Meng, et al. Evaluation on landslide susceptibility based on self-organizing feature map network and random forest model:A case study of Dayu County of Jiangxi Province[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33(1): 132-140. DOI: 10.16031/j.cnki.issn.1003-8035.2022.01-16
Citation: HE Shu, ABUDIKEYIMU XMSY, HU Meng, et al. Evaluation on landslide susceptibility based on self-organizing feature map network and random forest model:A case study of Dayu County of Jiangxi Province[J]. The Chinese Journal of Geological Hazard and Control, 2022, 33(1): 132-140. DOI: 10.16031/j.cnki.issn.1003-8035.2022.01-16

基于自组织特征映射网络-随机森林模型的滑坡易发性评价——以江西大余县为例

基金项目: 江西省教育厅科学技术研究项目(GJJ180436);江西省自然科学基金面上项目(20171BAB203029)
详细信息
    作者简介:

    何 书(1978-),男,贵州遵义人,博士,副教授,主要从事地质灾害防治与环境地质方面研究与教学工作。 E-mail:769844918@qq.com

  • 中图分类号: P642.2

Evaluation on landslide susceptibility based on self-organizing feature map network and random forest model:A case study of Dayu County of Jiangxi Province

  • 摘要: 为深入探讨评价单元和非滑坡样本选取对滑坡易发性预测的影响,构建了一种基于自组织特征映射网络-随机森林模型的滑坡易发性评价模型。该模型针对栅格单元和斜坡单元在滑坡易发性评价中的不足,结合栅格单元和斜坡单元的相互关系,提出了滑坡易发性指数的优化计算方法。在此基础上,基于随机森林Tree Bagger分类器构建滑坡易发性评价模型,通过对比分析自组织特征映射网络和随机方法选取非滑坡样本对评价结果的影响,探讨自组织特征映射网络、随机森林和自组织特征映射网络-随机森林三种评价模型的有效性;将评价模型应用于大余县滑坡易发性评价。结果显示,随机森林模型和自组织特征映射网络-随机森林模型的预测精度较高,分别达到91.19%和94.94%,成功率曲线的AUC值分别为0.822和0.849,表明自组织特征映射网络-随机森林模型具有更高的预测率和成功率, 自组织特征映射网络聚类的预测精度虽然有限,但作为非滑坡样本的选择方法,能够有效提高随机森林模型的评价精度。
    Abstract: In order to further explore the influence of evaluation units and non-landslide sample selection methods on landslide susceptibility prediction, a landslide susceptibility evaluation model is established based on self-organizing feature map network and random forest model in this paper. According to the relationship between grid units and slope units, an optimized calculation method of landslide susceptibility index is proposed. Aiming at the deficiencies of grid units and slope units in the evaluation of landslide susceptibility, this model proposes an optimized calculation method for landslide susceptibility index based on the relationship between grid cells and slope cells. On this basis, a landslide susceptibility evaluation model was established based on the random forest Tree Bagger classifier. By comparing and analyzing the influence of self-organizing feature map network and random non-landslide sample selection methods on the evaluation results, the effectiveness of the three evaluation models of self-organizing feature map network, random forest and self-organizing feature map network -random forest were discussed. The evaluation model has been applied to the landslide susceptibility evaluation in Dayu County. The results show that the prediction accuracy of random forest and self-organizing feature map network-random forest is higher, reaching 91.19% and 94.94% respectively, and the AUC of success rate curve was 0.822 and 0.849 respectively. It shows that self-organizing feature map network-random forest has higher prediction rate and success rate, although the prediction accuracy of self-organizing feature map network clustering is limited, it can effectively improve the evaluation accuracy of random forest model as the basis for selecting non landslide samples.
  • 中国是地质灾害高发的国家,各类地质灾害给人们的生命财产造成了巨大的损失[1]。滑坡作为一种主要的地质灾害,具有隐蔽性强、危害性大、突发性高等特点,广泛分布于中国山区和峡谷地带[2-4]。近年来,高山峡谷滑坡频频发生。因此,对高山峡谷地区的滑坡灾害开展早期识别,能够为防灾减灾事业及政府部门决策提供一种有效的手段。

    地表形变是反应当前坡体稳定性及运动状态最直接的物理量,因此,监测地表形变可以为探测滑坡等地质灾害的隐患点提供重要信息[1]。传统的地表形变监测方法主要采用精密水准测量、全球定位系统 (GPS) 等。但这些监测方法存在变形监测工作量大、费时、花费大、测点难以保存等缺陷,同时,传统的监测方法无法对大区域滑坡形变进行探测[5-6]。与传统地表形变相比,合成孔径雷达干涉测量 (Interferometric Synthetic Aperture Radar, InSAR) 技术是近几年发展起来的一种新型大地测量手段。具有覆盖范围广、穿透云层、全天候运作、精度高的特点,理论上可以获得非常精确的数字高程模型和毫米量级的地表形变信息,已被成功用于滑坡灾害监测[7-9]。目前,利用InSAR技术对滑坡早期识别的研究取得了一些成功范例:张路等[2]采用自主研发的相干散射体时序InSAR方法,成功识别出了17处持续变形中的不稳定坡体。Dai等[10]利用InSAR对潜在的滑坡进行早期识别,10个潜在滑坡区域被早期发现,结果表明,利用InSAR可以作为获取滑坡早期识别的有效手段。韩守富等[11]证明了InSAR技术在黄土高原地区地质灾害早期识别方面的适用性和准确性, 可以应用于黄土高原地区地质灾害隐患识别预警。冯文凯等[12]利用SBAS- InSAR技术对金沙江流域沃达村滑坡进行地表形变监测,表明SBAS- InSAR技术在复杂山区地质灾害监测预警领域有较为广阔的应用前景,为类似老滑坡监测预警提供了新的思路与借鉴。戴可人等[13]利用时间序列InSAR技术对雅砻江流域雅江县—木里县段的高山峡谷区域进行了滑坡灾害隐患广域早期识别,成功探测到8处隐患区域。以上方法能够有效识别滑坡灾害,但对于深切割高山峡谷区的滑坡早期识别,仅利用SAR单轨道数据监测会致使SAR成像几何畸变造成部分滑坡不能识别,只有通过升降轨数据结合的方式才能全面准确的对滑坡灾害进行早期识别。

    文章利用SBAS-InSAR技术,采用升降轨数据结合的方式对深切割高山峡谷区的滑坡进行早期解译识别,并根据识别结果,对不同潜在滑坡类型进行了分析与讨论,为高山峡谷地区防灾减灾及滑坡灾害早期识别提供一种更为全面的方法。

    小基线集(SBAS-InSAR)是在差分InSAR基础上发展起来的一种新的时间序列分析方法,能够降低相位噪声和误差[14]

    假定在时间$ {t}_{1} $$ {{t}}_{\text{s}} $内获取同一地区的S幅SAR影像,然后根据干涉组合条件,在短基线距的条件下形成N幅干涉条纹图,且满足:

    $$ \frac{S}{2}\ll N\ll \frac{S(S-1)}{2} $$ (1)

    $ t_{\text{A}} $$ t_{\text{B}} $$ t_{\text{A}} < t_{\text{B}} $)时刻两幅影像生成的第${{i}}({{i}}= $$ \mathrm{1,2},\cdots ,N)$幅干涉图,在去除平地及地形相位影响后,第$ x $个像素的干涉相位可表示为[14]

    $$ \begin{split} {\varphi ^i}_{{\rm{A,B}}}\left( {x,r} \right) =& {\varphi ^i}_{{\rm{def}}}\left( {x,r} \right) + {\varphi ^i}_{{\rm{topo}}}\left( {x,r} \right) + \\ &{\varphi ^i}_{{\rm{aps}}}\left( {x,r} \right) + {\varphi ^i}_{{\rm{orb}}}\left( {x,r} \right) + {\varphi ^i}_{{\rm{noise}}}\left( {x,r} \right) \end{split}$$ (2)

    式中:r——斜距;

    ${\varphi ^i}_{{\rm{def}}}\left( {x,r} \right)$——$t_{\text{A}} — t_{\text{B}}$对应的斜距向变形;

    ${\varphi ^i}_{{\rm{topo}}}\left( {x,r} \right)$——地形相位误差;

    ${\varphi ^i}_{{\rm{aps}}}\left( {x,r} \right)$——大气相位误差;

    ${\varphi ^i}_{{\rm{orb}}}\left( {x,r} \right)$——基线轨道引起的相位误差;

    ${\varphi ^i}_{{\rm{noise}}}\left( {x,r} \right)$——噪声误差。

    假定tk时刻和tk+1时刻不同干涉图间的形变速率为${{{v}}_{{{k,k + 1}}}}$,则$t_{\text{A}}— t_{\text{B}}$间的累积形变可表示为式:

    $$ \delta {\varphi _{{\rm{def}}}}\left( {x,r} \right) = \frac{{4\pi}}{\lambda }\sum\limits_{j = {t_\rm{{{A}}}}}^{{t_\rm{B}} - 1} {\left( {{t_{{k + 1}}} - {t_{k}}} \right){v_{k,k + 1}}} $$ (3)

    $ N $幅干涉条纹图进行三维时空相位解缠即可求出不同SAR获取时间的形变速率。

    文中利用SBAS-InSAR技术对深切割高山峡谷区滑坡灾害早期识别,方法的重点有两个关键:①针对深切割高山峡谷区的滑坡如何识别;②如何保证识别结果的准确性。

    对于深切割高山峡谷区的滑坡早期识别,由于地形条件等因素的影响,仅利用SAR单轨道数据监测会致使SAR成像几何畸变造成部分滑坡不能识别。为此,在技术上采用升降轨数据结合的方式,对深切割高山峡谷区滑坡灾害进行早期识别。弥补了SAR单轨道数据识别存在的不全面、不准确等弊端。

    为保证识别结果的准确性,仅凭形变监测结果无法区分是否为潜在的滑坡导致形变,为此,引入高分辨率光学影像对形变区域进行辅助识别,结合形变范围、高程、坡度、植被覆盖和坡体是否具有滑坡特征等进行识别,避免过度依赖形变结果导致的误判等问题。同时高分辨率光学影像能起到验证识别结果的作用,早期的滑坡一般会有一定的滑动痕迹,但不代表所有早期滑坡都有滑动痕迹。综上,引入高分辨率光学影像等作为辅助识别能有效保证识别结果的准确性。

    东川小江流域为世界典型暴雨泥石流区,被称为“泥石流的天然博物馆”,该地区泥石流的发生往往与滑坡密切相关,是属于暴雨型滑坡泥石流,即先滑坡,再经暴雨冲刷后形成泥石流[15]。文中选取小江沿线两侧高山峡谷作为研究区(图1),以小江为界,河谷凹陷,形成“V”字型,属于典型的深切割高山峡谷,东侧为牯牛寨山,最高峰海拔4017.3 m;西部为拱王山,最高峰海拔4344.1 m。该区域地势陡峭,其独特的地形和地质构造,致使局部区域暴雨多、土质松软、水土流失严重,导致该区域内地质灾害频发。

    图  1  研究区位置
    Figure  1.  Location of study area

    研究数据是从欧空局(European Space Agency, ESA)下载的40景C波段Sentinel-1A升降轨影像,时间跨度为2018年4月25日—2019年4月20日,极化方式为VV,成像模式为IW。数据参数如表1所示。为了提高影像轨道精度,引入了POD精密定轨星历数据,使用JAXA(日本宇宙航空研究开发机构)提供的ALOS WORLD 3D 30 m空间分辨率的数字高程模型(Digital Elevation Model, DEM),用于去除地形相位影响。拼接后的DEM如图2所示。

    表  1  Sentinenl-1A数据参数
    Table  1.  Sentinenl-1A data parameters
    轨道
    方向
    成像
    模式
    波段波长入射角
    /(°)
    极化
    方式
    距离
    分辨率/m
    方位
    分辨率/m
    重访
    周期/d
    升轨IWC5.6334.17VV52012
    降轨IWC5.6339.35VV52012
    下载: 导出CSV 
    | 显示表格
    图  2  研究区DEM
    Figure  2.  Digital elevation model of study area

    选取时间跨度2018年4月25日—2019年4月20日的Sentinel-1A斜距单视复数(Single Look Complex,SLC)影像,其中升降轨数据各20景。利用SARscape软件进行处理,升轨和降轨数据选取日期为2018-08-11和2018-08-13的影像作为超级主影像,通过设置临界基线和时间基线,生成103和100对像对,设置多视数为1∶4可以较好地抑制斑点噪声,采用Minimum Cost Flow解缠方法和Goldstein滤波方法做干涉工作流,最终生成干涉图,调整删除不理想的数据,在研究区生成较为理想的部分升降轨干涉图如图3所示。从图3可以看出,针对山区地区,运用升降轨数据处理后得到的干涉图相干性较为理想。

    图  3  部分升降轨干涉图
    Figure  3.  Part of the lifting rail interference diagram

    经过轨道精炼和重去平,估算和去除残余的恒定相位和解缠后还存在的相位坡道,进行第一次反演、第二次反演,最后对序列信息进行地理编码获得研究区域2018年4月25日—2019年4月20日雷达视线方向(LOS)的形变速率图。如图4所示,正值表示形变朝着卫星的方向,负值表示形变沿着远离卫星的方向。图4(a)为采用升轨数据获取的地表形变速率图,图4(b)表示采用降轨数据获取的地表形变速率图。从图4可以看出,采用升轨获取的形变主要分布在小江东侧,最大形变速率为−162.074 mm/a,降轨获取的形变主要分布在小江西侧西北方向,最大形变速率为−120.425 mm/a。致使不同轨道得到的形变结果不同,是由于升轨数据飞行方向大致从南到北,雷达视线(LOS)方向位于右侧,能够很好地将峡谷两侧由西向东的地表形变监测出来,相反,降轨数据的飞行方向与之相反,能够将峡谷两侧由东向西的地表形变监测出来。为此,利用不同轨道数据可以互补,使得监测结果更为准确全面,能够避免单一轨道带来的几何畸变等问题。

    图  4  升降轨获得的地表形变速率图
    Figure  4.  Surface deformation rate map obtained by lifting rail

    利用InSAR技术在植被覆盖地区进行地形测量非常困难,这主要是因为电磁场和(或)散射体的物理特性随时间的变化而造成的,植被覆盖度高地区,失相干严重,导致监测地面形变信息的能力较差,对形变监测精度的影响较大[16-17]。为考虑监测形变结果的有效性,引入归一化植被指数(Normalized Difference Vegetation Index,NDVI)对研究区植被覆盖数据进行分析。如图5所示,NDVI结果被限定在[−1,1],负值表示地面覆盖为云、水、雪等,对可见光高反射,0表示有岩石或裸土等,正值表示有植被覆盖,且随覆盖度增大而增大。对比图4图5可以发现,升降轨获取的部分形变区处于植被覆盖较高的地区,这正是由于植被造成失相干严重,致使这些区域监测结果不准确,在分析和识别滑坡时,需把其剔除,避免导致滑坡识别的误判。

    图  5  研究区NDVI
    Figure  5.  NDVI in the study area

    深切割高山峡谷地区同样存在非滑坡原因的沉降变形,为保证获取的形变信息为滑坡灾害,仅凭形变监测结果无法区分是否为潜在的滑坡导致形变,为此,引入高分辨率光学影像对形变区域进行辅助识别,结合形变范围、高程、坡度、植被覆盖和坡体是否具有滑坡特征等进行识别,避免过度依赖形变结果导致的误判等问题。由于不同轨道获取的形变数据能够对研究区不同方向的滑坡进行早期识别,分别对不同轨道获取的数据结合光学影像解译,升轨数据解译结果如图6所示,黑色区域为滑坡形变区域,红色虚线区域为对应光学影像滑坡灾害隐患,黄色箭头代表滑坡方向,共解译出11处潜在滑坡,分别用${\rm{H}}01\text{、}{\rm{H}}02\text{、}\cdots\cdots\text{、}{\rm{H}}11$进行编号,其中H01、H02、H08和H09位于现存泥石流沟处,H01位于尖山沟老尖山、土岩子一带,H02位于蒋家沟一带,H08、H09位于大白泥沟和小白泥沟一带,其他潜在滑坡有一定的滑坡痕迹。识别出来的滑坡区域平均最大形变速率超过−40 mm/a,最大形变速率−120 mm/a,位于H06岩子脚村上方。根据滑坡是否直接威胁附近周围村落,对11处潜在滑坡进行危险划分,共识别出4处高风险滑坡,分别是H06、H09、H10和H11,H06坡体对岩子脚村构成威胁,H09坡体对大村和鲁纳窝村构成威胁,H10对小多红村和妥托村构成威胁,H11对小凹子村、厂上村、大麦地、小村子等构成威胁,这些区域一旦发生滑坡,极有可能对当地群众造成生命及财产威胁。

    图  6  升轨数据滑坡灾害识别结果
    Figure  6.  Landslide disaster identification results of lifting data

    降轨数据解译结果如图7所示,共解译出7处潜在滑坡,分别用${\rm{H}}12、{\rm{H}}13\text{、}\cdots\cdots\text{、}{\rm{H}}18$进行编号,滑坡区域平均形变速率均超过−40 mm/a,最大形变区域位于H18大白泥沟和小白泥沟连接处,其最大形变速率为−94.11 mm/a。根据滑坡是否直接威胁附近周围村落,对7处潜在滑坡区进行危险划分,共识别出1处滑坡高风险区,编号为H16,该坡体对姑庄村和新寨田村构成威胁。对比图6图7发现,不同轨道识别出来的潜在滑坡分布不同,利用升轨获取的滑坡主要分布在小江右侧峡谷中部,降轨获取的滑坡主要分布在小江左侧西北方向,其中不同轨道识别到的潜在滑坡有相同位置,但形变区域不一致,例如,H01、H13共处于尖山沟一带,H02、H15共处于蒋家沟一带,H08、H18共处于大小白泥沟一带。由此证明,利用升降轨结合的方式能够有效识别深切割高山峡谷地区不同坡度方向存在的潜在滑坡,避免了单一轨道存在的识别结果不准确,不全面等问题。

    图  7  降轨数据滑坡灾害识别结果
    Figure  7.  Landslide disaster identification results of rail descent data

    经过对比光学影像发现,研究区识别出来的潜在滑坡可定义为三种类型,分别是处于泥石流区域的潜在滑坡(H09)、有滑坡痕迹的潜在滑坡(H06)和无滑坡痕迹的潜在滑坡(H16)。为有效了解不同类型潜在滑坡的形变趋势,分别对三种典型的潜在滑坡进行分析。

    H09属于典型的泥石流区域潜在滑坡,位于小白泥沟,该区域为暴雨型泥石流多发区。利用获取的形变速率值与三维光学影像叠置分析,其形变速率如图8所示,存在A、B和C三个潜在滑坡区,呈V字型分布,滑动方向由坡面向下滑动,在每个潜在滑坡上方能明显的观察到断裂滑动痕迹,最大形变速率为−103.013 mm/a,位于潜在滑坡A面上方。从图8中可以看到,B、C面分别位于小白泥沟两侧,其潜在滑坡属于长期滑坡,潜在滑坡上端形变速率大,下端由于滑坡堆积物导致一定的抬升。但A面潜在滑坡还未形成真正意义的滑坡,一旦A面滑坡,极有可能威胁大村和鲁纳窝村民的生命及财产安全。为更有效全面地分析该种类型潜在滑坡的特点,选取A面进行详细分析,可以看到A面呈现三处不均匀形变,有明显的拉张裂缝且有滑坡痕迹,其中P1位于整个潜在滑坡顶端,P2位于潜在滑坡中部,P3位于潜在滑坡右侧。三个点均有可能在后期的降雨等因素下发展为新的潜在滑坡区。图9为3处潜在滑坡特征点形变量与当月降雨量的关系图,可以看到特征点的形变量与降雨量有一定的相关关系。研究区总体位于小江断裂带沿线,该区域断裂带宽5~20 km,呈现弱剪切强挤压活动特征[18],以挤压穹起隆升变形为主[19],从形变时序图9可以看出,三个特征点在2018年5月先抬升后沉降,抬升可能是由于断裂带挤压活动致使,后经过雨水冲刷导致后期慢慢滑动沉降。P1、P2和P3点的平均形变量分别为−10.571 mm、−22.564 mm、−19.516 mm,年形变量分别为−49.063 mm、−58.740 mm、−62.635 mm。该区域月均降雨量为102 mm。随着时间的推移,形变量逐渐增加,到一定临界值便形成真正意义上的滑坡。

    图  8  H09潜在滑坡形变速率图
    Figure  8.  Deformation rate diagram of H09 potential landslide
    图  9  P1—P3时间序列曲线与降雨量
    Figure  9.  Time series curve of P1—P3 and rainfall

    H06属于典型的有滑坡痕迹的潜在滑坡,位于牯牛山岩子脚上方。图10为该潜在滑坡形变速率与三维影像的叠置图,从图10中可以白色虚线内两处潜在滑坡有明显的滑动痕迹,滑动方向为坡顶沿着坡脚移动。最大形变速率位于P5附近,达到−127.093 mm/a。选取特征点P4、P5形变量与月降雨量构建关系图(图11)。随着时间推移,形变量不断增大。特征点的形变量与月降雨量有一定的相关关系,从形变时序图11可以看出,两个特征点在2018年5月先抬升后沉降,该潜在滑坡位于小江断裂带沿线,可能是由于断裂带挤压活动致使其抬升。当月降雨量增多加剧潜在滑坡的形变,这是由于降雨致使表层松散土体随着雨水滚落流失。P4、P5点的平均累积形变量分别为−37.562 mm、−42.054 mm,年形变速率分别为−105.903 mm/a、−112.469 mm/a。该潜在滑坡区最高形变点高程为3606 m,地势陡峭,从年形变量可以看出该区域存在很大安全隐患,一旦发生滑坡,岩子脚村村民将面临极大的生命及财产安全。

    图  10  H06潜在滑坡形变速率图
    Figure  10.  Deformation rate diagram of H06 potential landslide
    图  11  P4—P5时间序列曲线与降雨量
    Figure  11.  P4—P5 time series curve and rainfall

    H16属于典型的无滑坡痕迹的潜在滑坡,位于姑庄村和新寨田上方。图12为该潜在滑坡区形变速率与三维影像的叠置图,从图12中可以看到,位于潜在滑坡上方有明显的形变,最大形变速率为−80.141 mm/a。通过光学影像看到山体表面并无明显的滑动痕迹。该潜在滑坡区最高高程为2304 m,坡度大,滑动方向大致沿着山沟滑动。选取特征点P6、P7形变量与月降雨量构建关系图,如图13所示,形变与降雨有一定关系。从形变时序图13可以看出,P6、P7两个特征点总是先抬升后沉降再抬升这一反复过程,在2018年5月可能是由于断裂带挤压活动致使其抬升,其他月份抬升微小可能是由于雨水冲刷导致土地堆积抬升,再经雨水冲刷致使微小沉降。P6、P7平均累积形变量分别为−29.219 mm、−21.281 mm,年形变速率分别为−63.788 mm/a、−53.110 mm/a。两个特征点形变曲线大致一直,从潜在滑坡上方到中部形变速率逐渐减小。该潜在滑坡一旦发生滑坡,将直接威胁姑庄村和新寨田两个村。

    图  12  H16潜在滑坡形变速率图
    Figure  12.  Deformation rate diagram of H16 potential landslide
    图  13  P6—P7时间序列曲线与降雨量
    Figure  13.  P6—P7 time series curve and rainfall

    文中基于SBAS-InSAR技术,采用Sentinel-1升降轨数据结合互补的方式对东川小江沿线两侧深切割高山峡谷区滑坡灾害进行早期识别实验,得出以下结论:

    (1)升轨数据识别出的滑坡主要分布在小江右侧,降轨反之,这表明相比单一轨道获取的识别结果,利用升降轨结合的方式能够更全面的监测和识别高山峡谷滑坡,避免单轨道对高山峡谷区滑坡进行早期识别存在SAR成像几何畸变造等问题。

    (2)在高山峡谷地区植被覆盖度过高容易致使失相干,同时,有部分地区并不属于滑坡导致的形变,为此,引入光学影像,结合形变范围、高程、坡度、植被覆盖和坡体是否具有滑坡特征等进行识别,避免过度依赖形变结果导致的误判等问题。

    (3)通过H09、H06、H16 等3个典型潜在滑坡灾害分析,可以看到滑坡的形成和降雨量有一定的关系。

    (4)文中共识别出18处潜在滑坡区,其中H06、H09、H10、H11、H16等5处为滑坡高风险区,证明利用该方法可作为高山峡谷区滑坡灾害识别的有效手段,但目前还有一定的不足,比如在高山峡谷地区SAR数据量过少致使部分滑坡并未获取到形变量等问题。

  • 图  1   SOM-随机森林模型流程

    Figure  1.   Flow chart of SOM-vandom forest model

    图  2   研究区地理位置及滑坡空间分布图

    Figure  2.   Location of the study area and landslide distribution

    图  3   各指标特征空间分布

    Figure  3.   Spatial distribution of the factors

    图  4   滑坡易发性分区结果

    Figure  4.   Zoning map of landslide susceptibility

    图  5   模型误差与决策树数量关系

    Figure  5.   Relationship between model error and decision tree quantity

    图  6   滑坡易发性预测成功率曲线

    Figure  6.   Success rate curves of landslide susceptibility indexes

    表  1   评价指标体系

    Table  1   Evaluation index system

    滑坡影响因素分类及分级
    高程(1)[130,200);(2)[200,300);(3)[300,400);(4)[400,500);(5)[500,600);(6)[600,700);(7)[700,800);(8)[800,900);(9)[900,1346);
    植被归一化指数(1)[−1,0.2);(2)[0.2,0.4);(3)[0.4,0.6);(4)[0.6,0.8);(5)[0.8,1.0)
    土地利用类型(1)矿山工程用地;(2)乔木;(3)耕地和荒地;(4)城镇用地
    坡度/(◦)(1)[0,10);(2)[10,20);(3)[20,30);(4)[30,40);(5)[40,54.6)
    总曲率(1)[−6.3,−0.82);(2)[−0.82,−0.27);(3)[−0.27,−0.18);(4)[−0.18,0.88);(5)[0.88,6.44)
    岩土类型(1)C:白云质灰岩类;(2)D:砂岩类;(3)Q:松散沉积物;(4)H:变余砂岩类;
    (5)K:红色砂砾岩类;(6)R:岩浆岩类;(7)Z:板岩千枚岩
    道路密度/(km∙km-2)(1)[0,0.393);(2)[0.393,0.786);(3)[0.786,1.178);(4)[1.178,1.571);(5)[1.571,1.964)
    距道路的距离/m(1)[0,50);(2)[50,100);(3)[100,150);(4)[150,200);(5)≥200
    距水系的距离/m(1)[0,50);(2)[50,100);(3)[100,150);(4)[150,200);(5)≥200
    距断层的距离/m(1)[0,200);(2)[200,400);(3)[400,600);(4)[600,800);(5)≥800
    下载: 导出CSV

    表  2   不同滑坡易发性分区的滑坡频率

    Table  2   Landslide frequency in different landslide susceptibility zones

    序号滑坡
    易发性
    等级
    SOM 随机森林 SOM-随机森林
    分区面积/km2滑坡频率
    /(个·km−2
    比例/% 分区面积/km2滑坡频率
    /(个·km-2
    比例/% 分区面积/km2滑坡频率
    /(个·km−2
    比例/%
    1372.450.00811 561.5200 545.460.00000
    2较低495.010.06267.76 113.700.01763.94 122.060.00821.84
    3中等288.020.180522.37 137.650.02184.88 139.190.01443.22
    4较高61.030.278534.5 152.830.091620.49 152.120.118326.53
    5122.530.277534.37 373.340.316170.70 380.210.305168.41
    下载: 导出CSV
  • [1] 黄发明, 叶舟, 姚池, 等. 滑坡易发性预测不确定性: 环境因子不同属性区间划分和不同数据驱动模型的影响[J]. 地球科学,2020,45(12):4535 − 4549. [HUANG Faming, YE Zhou, YAO Chi, et al. Uncertainties of landslide susceptibility prediction: different attribute interval divisions of environmental factors and different data-based models[J]. Earth Science,2020,45(12):4535 − 4549. (in Chinese with English abstract)
    [2]

    PEETHAMBARAN B, ANBALAGAN R, KANUNGO D P, et al. A comparative evaluation of supervised machine learning algorithms for township level landslide susceptibility zonation in parts of Indian Himalayas[J]. Catena,2020,195:104751. DOI: 10.1016/j.catena.2020.104751

    [3] 周超, 殷坤龙, 曹颖, 等. 基于集成学习与径向基神经网络耦合模型的三峡库区滑坡易发性评价[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)
    [4]

    BAHARVAND S, RAHNAMARAD J, SOORI S, et al. Landslide susceptibility zoning in a catchment of Zagros mountains using fuzzy logic and GIS[J]. Environmental Earth Sciences,2020,79(10):1 − 10.

    [5] 杨永刚, 殷坤龙, 赵海燕, 等. 基于C5.0决策树-快速聚类模型的万州区库岸段乡镇滑坡易发性区划[J]. 地质科技情报,2019,38(6):189 − 197. [YANG Yonggang, YIN Kunlong, ZHAO Haiyan, et al. Landslide susceptibility evaluation for township units of bank section in Wanzhou district based on C5.0 decision tree and K-means cluster model[J]. Geological Science and Technology Information,2019,38(6):189 − 197. (in Chinese with English abstract)
    [6] 武雪玲, 沈少青, 牛瑞卿. GIS支持下应用PSO-SVM模型预测滑坡易发性[J]. 武汉大学学报(信息科学版),2016,41(5):665 − 671. [WU Xueling, SHEN Shaoqing, NIU Ruiqing. Landslide susceptibility prediction using GIS and PSO-SVM[J]. Geomatics and Information Science of Wuhan University,2016,41(5):665 − 671. (in Chinese with English abstract)
    [7] 夏辉, 殷坤龙, 梁鑫, 等. 基于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)
    [8] 刘坚, 李树林, 陈涛. 基于优化随机森林模型的滑坡易发性评价[J]. 武汉大学学报(信息科学版),2018,43(7):1085 − 1091. [LIU Jian, LI Shulin, CHEN Tao. Landslide susceptibility assesment based on optimized random forest model[J]. Geomatics and Information Science of Wuhan University,2018,43(7):1085 − 1091. (in Chinese with English abstract)
    [9] 杜国梁,杨志华,袁颖,等. 基于逻辑回归-信息量的川藏交通廊道滑坡易发性评价[J].水文地质工程地质,2021,48(5):102-111.

    DU Guoliang, YANG Zhihua, YUAN Ying, et al. Landslide susceptibility mapping in the Sichuan-Tibet traffic corridor using logistic regression-information value method[J]. Hydrogeology & Engineering,2021,48(5):102-111.(in Chinese with English)

    [10] 连志鹏, 徐勇, 付圣, 等. 采用多模型融合方法评价滑坡灾害易发性: 以湖北省五峰县为例[J]. 地质科技通报,2020,39(3):178 − 186. [LIAN Zhipeng, XU Yong, FU Sheng, et al. Landslide susceptibility assessment based on multi-model fusion method: A case study in Wufeng County, Hubei Province[J]. Bulletin of Geological Science and Technology,2020,39(3):178 − 186. (in Chinese with English abstract)
    [11] 郭子正, 殷坤龙, 付圣, 等. 基于GIS与WOE-BP模型的滑坡易发性评价[J]. 地球科学,2019,44(11):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(11):4299 − 4312. (in Chinese with English abstract)
    [12] 鲜木斯艳·阿布迪克依木, 何书. 基于MIV-BP神经网络的滑坡易发性空间预测[J]. 人民长江,2019,50(12):140 − 144. [ABUDIKEYIMU XMSY, HE Shu. Spatial prediction on landslide vulnerability based on MIV-BP neural network[J]. Yangtze River,2019,50(12):140 − 144. (in Chinese with English abstract)
    [13] 李文娟, 邵海. 基于遥感影像多尺度分割与地质因子评价的滑坡易发性区划[J]. 中国地质灾害与防治学报,2021,32(2):94 − 99. [LI Wenjuan, SHAO Hai. Landslide susceptibility assessment based on multi-scale segmentation of remote sensing and geological factor evaluation[J]. The Chinese Journal of Geological Hazard and Control,2021,32(2):94 − 99. (in Chinese with English abstract)
    [14]

    PARYANI S, NESHAT A, JAVADI S, et al. Comparative performance of new hybrid ANFIS models in landslide susceptibility mapping[J]. Natural Hazards,2020,103(2):1961 − 1988. DOI: 10.1007/s11069-020-04067-9

    [15] 谈树成, 赵娟娟, 杨林, 等. 基于GIS和信息量-快速聚类模型的滑坡易发性研究: 以云南省福贡县为例[J]. 云南大学学报(自然科学版),2018,40(6):1148 − 1158. [TAN Shucheng, ZHAO Juanjuan, YANG Lin, et al. Evaluation of landslide susceptibility based on GIS and the information value & the K-means cluster model—Take Fugong County as an Example[J]. Journal of Yunnan University (Natural Sciences Edition),2018,40(6):1148 − 1158. (in Chinese with English abstract)
    [16] 孙长明, 马润勇, 尚合欣, 等. 基于滑坡分类的西宁市滑坡易发性评价[J]. 水文地质工程地质,2020,47(3):173 − 181. [SUN Changming, MA Runyong, SHANG Hexin, et al. Landslide susceptibility assessment in Xining based on landslide classification[J]. Hydrogeology & Engineering Geology,2020,47(3):173 − 181. (in Chinese with English abstract)
    [17] 郭子正, 殷坤龙, 黄发明, 等. 基于滑坡分类和加权频率比模型的滑坡易发性评价[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)
    [18] 黄发明, 殷坤龙, 蒋水华, 等. 基于聚类分析和支持向量机的滑坡易发性评价[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. (in Chinese with English abstract)
    [19] 李颖, 应保胜, 容芷君, 等. 基于SOM的产品设计结构模块划分及其评价[J]. 武汉科技大学学报,2018,41(4):301 − 306. [LI Ying, YING Baosheng, RONG Zhijun, et al. SOM-based product design structure clustering and its evaluation[J]. Journal of Wuhan University of Science and Technology,2018,41(4):301 − 306. (in Chinese with English abstract)
    [20] 王小川, 史峰, 郁磊. MATLAB神经网络43个案例分析[M]. 北京: 北京航空航天大学出版社, 2013.

    WANG Xiaochuan, SHI Feng, YU Lei. Analysis of 43 cases of neural network in matlab [M]. Beijing: Beijing University of Aeronautics & Astronautics Press, 2013. (in Chinese)

    [21] 林荣福, 刘纪平, 徐胜华, 等. 随机森林赋权信息量的滑坡易发性评价方法[J]. 测绘科学,2020,45(21):131 − 138. [LIN Rongfu, LIU Jiping, XU Shenghua, et al. Evaluation method of landslide susceptibility based on random forest weighted information[J]. Science of Surveying and Mapping,2020,45(21):131 − 138. (in Chinese with English abstract)
    [22] 李亭, 田原, 邬伦, 等. 基于随机森林方法的滑坡灾害危险性区划[J]. 地理与地理信息科学,2014,30(6):25 − 30. [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. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-0504.2014.06.006
    [23] 颜阁, 梁收运, 赵红亮. 基于GIS的斜坡单元划分方法改进与实现[J]. 地理科学,2017,37(11):1764 − 1770. [YAN Ge, LIANG Shouyun, ZHAO Hongliang. An approach to improving slope unit division using GIS technique[J]. Scientia Geographica Sinica,2017,37(11):1764 − 1770. (in Chinese with English abstract)
    [24] 罗路广, 裴向军, 黄润秋. 强震山区地震滑坡发生概率研究: 以九寨沟国家地质公园为例[J]. 岩石力学与工程学报,2020,39(10):2079 − 2093. [LUO Luguang, PEI Xiangjun, HUANG Runqiu. Earthquake-triggered landslide occurrence probability in strong seismically mountainous areas: A case study of Jiuzhaigou National Geopark[J]. Chinese Journal of Rock Mechanics and Engineering,2020,39(10):2079 − 2093. (in Chinese with English abstract)
    [25] 黄启乐, 陈伟, 唐绪波, 等. 区域地质灾害评价中斜坡单元划分方法研究[J]. 自然灾害学报,2017,26(5):157 − 164. [HUANG Qile, CHEN Wei, TANG Xubo, et al. Study on the method of slope unit zoning in regional geo-hazards risk assessment[J]. Journal of Natural Disasters,2017,26(5):157 − 164. (in Chinese with English abstract)
    [26] 程温鸣, 彭令, 牛瑞卿. 基于粗糙集理论的滑坡易发性评价: 以三峡库区秭归县境内为例[J]. 中南大学学报(自然科学版),2013,44(3):1083 − 1090. [CHENG Wenming, PENG Ling, NIU Ruiqing. Landslide susceptibility assessment based on rough set theory: Taking Zigui County territory in Three Gorges Reservoir for example[J]. Journal of Central South University (Science and Technology),2013,44(3):1083 − 1090. (in Chinese with English abstract)
    [27] 宫清华, 黄光庆, 张冬良, 等. 基于斜坡单元的浅层滑坡风险区划: 以华南松岗河小流域为例[J]. 安全与环境学报,2017,17(2):615 − 620. [GONG Qinghua, HUANG Guangqing, ZHANG Dongliang, et al. On the application of the slope element unit method to the shallow landslide risk-prevention mapping with a case study sample of Songgang watershed[J]. Journal of Safety and Environment,2017,17(2):615 − 620. (in Chinese with English abstract)
    [28] 薛强, 张茂省, 高波. 斜坡单元支持下基于土体含水率的陕西省清涧县城区黄土滑坡危险性评价[J]. 中国地质,2020,47(6):1904 − 1914. [XUE Qiang, ZHANG Maosheng, GAO Bo. Hazard assessment of loess landslide based on soil moisture content and supported by slope unit in Qingjian City, Shaanxi Province[J]. Geology in China,2020,47(6):1904 − 1914. (in Chinese with English abstract) DOI: 10.12029/gc20200624
    [29]

    PHAM B T, TIEN BUI D, PRAKASH I, et al. Rotation forest fuzzy rule-based classifier ensemble for spatial prediction of landslides using GIS[J]. Natural Hazards,2016,83(1):97 − 127. DOI: 10.1007/s11069-016-2304-2

    [30] 武雪玲, 任福, 牛瑞卿. 多源数据支持下的三峡库区滑坡灾害空间智能预测[J]. 武汉大学学报(信息科学版),2013,38(8):963 − 968. [WU Xueling, REN Fu, NIU Ruiqing. Spatial intelligent prediction of landslide hazard based on multi-source data in Three Gorges Reservoir area[J]. Geomatics and Information Science of Wuhan University,2013,38(8):963 − 968. (in Chinese with English abstract)
  • 期刊类型引用(6)

    1. 翟越,谢怡帆,汪铁楠,王铭,张韵生. 基于改进突变理论的岩体崩塌危险性评价. 甘肃科学学报. 2023(02): 63-70 . 百度学术
    2. 鲁晓,祁生文,郑博文,郭忻怡,李永超,郭松峰,邹宇,唐凤娇,姚翔龙,宋帅华,马丽娜,张琳鑫,刘方翠,罗光明,梁宁,台大平. 川藏交通廊道崩滑灾害分布及其危险性评价. 工程地质学报. 2023(03): 718-735 . 百度学术
    3. 张广泽,陈国庆,王哲威,王栋,裴本灿. 川藏铁路雅安至昌都段高陡危岩发育特征与快速评价. 工程科学与技术. 2022(02): 30-40 . 百度学术
    4. 贺攀,郭荣昌,张蕊,余岭燕. 基于不同评价因子组合的铁路沿线滑坡危险性评价. 兰州交通大学学报. 2022(05): 34-41 . 百度学术
    5. 王文沛,殷跃平,胡卸文,李滨,刘明学,祁小博. 碎屑流冲击下桩梁组合结构拦挡效果及受力特征研究. 地质力学学报. 2022(06): 1081-1089 . 百度学术
    6. 戚利荣,王家鼎,张登飞,张永双,李贞孝,孙嘉兴,马剑飞. 冻融循环作用下花岗岩损伤的宏微观尺度研究. 水文地质工程地质. 2021(05): 65-73 . 百度学术

    其他类型引用(5)

图(6)  /  表(2)
计量
  • 文章访问数:  489
  • HTML全文浏览量:  216
  • PDF下载量:  256
  • 被引次数: 11
出版历程
  • 收稿日期:  2021-04-14
  • 修回日期:  2021-06-21
  • 网络出版日期:  2022-02-13
  • 刊出日期:  2022-02-24

目录

/

返回文章
返回