基于信息量模型的云南东川泥石流易发性评价

孙 滨，祝传兵，康晓波，叶 雷，刘 益

（1.昆明理工大学国土资源工程学院，云南 昆明 650093；2.自然资源部高原山地地质灾害预报预警与生态保护修复重点实验室，云南 昆明 650093；3.云南省地质环境监测院，云南 昆明 650216；4.云南滇湘矿业咨询有限公司，云南 昆明 650217）

0 引言

泥石流作为我国自然灾害重要类型之一，具有暴发突然、波及范围广、破坏性大等特点。泥石流易发性评价也从最早的定性分析发展到定量分析，一直以来是国内外地学专家学者研究的热点[1-2]，其评价质量的好坏与评价指标的选取、评价模型的确定有着密不可分的关系[3-8]。地形地貌、地质构造及人类活动等是泥石流易发性的主要影响因素，也成为其评价指标选取的重要参考依据[9-10]。常用的泥石流易发性评价模型有专家系统模型（层次分析、专家打分）、数理统计模型（信息量、证据权和逻辑回归等）和机器学习模型（决策树、随机森林、神经网络和支持向量机等）[11-18]，这些模型各有所长，同时也存在一定的缺陷。相对于专家系统模型受人为因素影响明显、机器学习模型存在参数调试较难等问题，信息量模型作为《地质灾害风险调查评价技术要求》（1∶50 000）（试行）推荐的方法具有操作简单、应用广泛和客观性好等优势，能够科学的对各指标区间分级[19]。

东川区地质灾害发育，尤以泥石流分布广泛，危害最为严重，素有“世界泥石流天然博物馆”之称，是全国泥石流危害最严重的地区。据不完全统计，解放以来泥石流已造成人员伤亡300 余人，直接经济损失近4 亿元，并且对生态环境破坏明显，严重制约了东川经济与社会发展。因此，在东川建立科学合理的泥石流易发性评价模型，为泥石流精准防控提供支撑，意义重大[20-21]。

文中在对东川区地质灾害详细调查成果系统梳理基础上，细致分析东川泥石流流域特征、动力学特性与形成机理等，选取坡度、坡向、起伏度、曲率、工程岩组、距水系距离、距断层距离和土地利用类型9 个评价指标，基于信息量模型和GIS 平台技术，以小流域为单元对东川泥石流进行了易发性评价，以期为东川防灾减灾工作提供参考。

1 研究区概况及数据来源

1.1 研究区概况

东川区位于云南省昆明市最北端，面积1 858.79 km2，属侵蚀剥蚀构造地貌，地势呈东西高中间河谷低、南高北低的特征，最大相对高差达3 600 m。区域上属金沙江流域，区域年降水量701～1 163 mm，降水量的85%以上集中在5—10月。区内地层岩性以元古代昆阳群碎屑岩分布最广，约占全区面积的43.68%，古生代碎屑岩夹碳酸盐岩约占19.10%，古生代玄武岩、中生代碎屑岩约占32.36%，新生代黏土岩和松散碎石土等约占4.86%。区内地质构造复杂，主构造线为南北向的小江断裂带。复杂的地形地貌、特殊的地质背景及独特的气候条件导致了东川泥石流频发。

通过收集东川区地质灾害详查、隐患排查和风险普查等数据，结合InSAR 地质灾害隐患早期识别成果和精细化调查野外验证，系统梳理出研究区典型泥石流144 条（图1），以此作为样本数据，开展研究区泥石流易发性定量评价。

图1 研究区泥石流灾害点分布Fig.1 Distribution of debris flow in study area

1.2 数据来源

DEM 数据收集自阿拉斯加卫星设备，制作高程、坡度、坡向、起伏度和曲率5 个评价指标；水系数据收集自OSM 并与DEM 提取水系进行比对制作距水系距离评价指标；地质数据收集自全国地质资料馆东川幅、会理幅、曲靖幅和武定幅20 万地质图[22]制作工程岩组分类和距断层距离2 个评价指标，土地类型数据收集自欧空局官网制作土地利用类型评级指标（表1）。

表1 数据来源及类型Table 1 Data source and types

2 研究方法

2.1 流域单元划分

泥石流易发性常用栅格作为评价单元，虽易于划分和模型计算，但忽略其流域特性，不能有效反映真实的泥石流情况，无法建立合理的评价模型，得出精准的评价结果。文中以收集的12.5 m×12.5 m 的DEM 栅格数据为基础，利用ArcGIS10.8 模型构建器将流域划分过程（填洼→流向→流量→提取河流网络→栅格河网矢量化→盆域分析→栅格转面）模块化，化繁为简，实现自动化生成流域，通过反复调整集水阈值，最终得到最符合实际的流域划分结果。将研究区划分为961 个流域，其中平均流域面积为1.94 km2，最大流域面积为10.62 km2，最小流域面积为0.44 km2。

2.2 信息量模型基本原理

信息量模型的理论基础是信息论[23]，运用概率论和数理统计的方法以信息熵的概念来分析各种评价指标作用下泥石流易发性的模型。通过各种评价指标与泥石流灾害点空间叠加分析，依托GIS 平台计算其单个指标信息量，然后进行多个指标的加权叠加得到综合信息量，从而建立泥石流易发性评价模型，其信息量值越大，表明易发性越高。

式中：I——各种评价指标加权的总信息量，可作为泥石流易发性指数；

Nj——单个评价指标特定分级区间内含有泥石流的个数；

N——泥石流总数量；

Sj——单个评价指标特定分级区间内栅格数；

S——总栅格数。

3 研究区泥石流易发性评价

3.1 评价指标选取原则

泥石流的形成影响因素众多，其演化是一个复杂的多因素耦合作用的地质过程。东川自1961年中科院建立野外观测站至今，对东川泥石流动力地貌过程与区域演化规律、运动学与动力学特征、流体物理力学与流变特性等方面取得了一系列国际先进水平成果。本文在系统分析东川泥石流触发因素基础上，结合大量专家学者泥石流易发性评价研究成果进行了评价指标的优选，选取坡度、坡向、起伏度、曲率、工程岩组、距水系距离、距断层距离和土地利用类型9 个评价指标。

3.2 评价流程

首先通过GIS 平台得到各评价指标的状态分级[24-25]，然后对其分级区间进行重分类，与泥石流灾害点图层进行空间叠加分析统计，计算各评价指标图层的信息量值；采用栅格转面-空间连接-面转栅格实现将各评价指标的状态分级信息量值赋值到栅格图层中，运用空间分析工具叠加各评价指标的信息量栅格图层获取总信息量值，并按流域单元划分总信息量栅格图层，以子流域总信息量的平均值作为该子流域的信息量值，并按自然间断法对其进行重分类，实现研究区的泥石流易发性分区。

3.3 评价指标状态分级

因各评价指标的状态分级对信息量模型精度影响较大，科学合理划分其状态分级至关重要。文中在处理各评价指标状态分级时首先采用自然间断法一般将其细分为10～15 级，与泥石流灾害点图层进行叠加分析，得出其信息量值，然后优先将信息量值为零的分级与相邻分级合并，接着将信息量值相近的分级与相邻分级合并，最后计算归并后的分级信息值，从而实现最优各评价指标状态分级（图2、表2）。

表2 各因素状态信息量表Table 2 Weighted information values of each factor

（1）高程

以研究区12.5 m 精度的DEM 栅格数据为基础，最高点为拱王山雪岭，海拔4 344 m，最低点位于金沙江与小江交汇处，海拔660 m，高差大于3 600 m。以660～4 344 m 为区间，将高程分为5 级，见图2（a）。泥石流主要分布在660～3 000 m，占泥石流总数量的84.34%，见图3（a）。

（2）坡度

从0°起，以10°为间隔，将坡度分为5 级，见图2（b）。泥石流点在坡度各分级均有分布，主要集中在20°～40°，见图3（b）。

（3）坡向

将坡向划为10 个分级，见图2（c）。泥石流在南、西南、东南方向比例为34.57%，在北、北西、北东方向为37.93%，可见泥石流在北坡比例大于南坡，见图3（c）。

图2 各评价因子状态分级图Fig.2 State classification of the factors

（4）地形起伏度

地形起伏度可以直观的反映山体的相对高差，是划分地貌类型的一个重要指标，同时还能有效地体现人类活动与地质灾害发育程度的相关性。借助ArcGIS 平台提取出研究区地形起伏度在0～441 m，分为0～20 m平坦起伏、20～40 m 微起伏、40～60 m 小起伏，＞60 m较大起伏4 个分级，见图2（d）。泥石流在微起伏占比最高，为40.84%，见图3（d）。

（5）曲率

曲率主要是用来反映地形弯曲程度的指标，将曲率划为4 个分级，见图2（e）。泥石流主要分布在曲率-1～2，占泥石流总数量的69.52%，见图3（e）。

（6）工程岩组分类

以研究区20 万地质图为基础，根据地层岩性的工程地质特性，将研究区工程岩组划分为4 大类，见图2（f）。泥石流主要分布在较软岩组内，占泥石流总数量的56.47%，见图3（f）。

（7）距水系距离

研究区水系发育，分布有小江流域和普渡河流域，均属金沙江水系。以研究区内河流为中心作200 m、400 m、＞400 m 缓冲区分析，见图2（g）。结果表明，在距河流400 m 范围内泥石流点分布集中，占泥石流总数量的78.08%，见图3（g）。

（8）距断层距离

研究区位于川滇菱形块体内部小江断裂与普渡河断裂夹持地带，断层发育。以研究区1∶20 万构造纲要图为基础，以断层为中心作1 km、2 km、3 km、＞3 km 缓冲区分析，见图2（h）。分析表明泥石流主要位于断层0～2 km 范围内，占比84.03%，见图3（h）。

（9）土地利用类型

土地利用数据采用欧空局（ESA）发布的2020年10 m 分辨率的全球土地利用（World Cover）数据。研究区土地利用主要包括林地、灌木、草地、耕地、建筑用地、裸地/稀疏植被区和开阔水域7 类，见图2（i）。泥石流主要集中分布在林地、草地和耕地内，占比88.82%，见图3（i）。

3.4 基于GIS 的信息量计算和易发性评价

评价指标体系建立完成后，基于GIS 平台将各评价指标状态分级（图2）与泥石流灾害点进行空间叠加分析，得出各评价指标状态分级区间内泥石流分布数量（图3），根据信息量模型公式进行各评价指标信息量值计算，然后通过栅格叠加分析计算总信息量值，以自然间断法将其划分为4 个等级，得到研究区泥石流易发性分区图，见图4（a）。低易发区[-14.15，-5.60]、中易发区[-5.60，-2.72]、高易发区[-2.72，0.42]及极高易发区[0.42，8.10]，占东川区全域栅格比例为18.49%、32.43%、33.30%和15.78%，其灾点占栅格频率比分别为0.038、0.150、0.542、4.843，极高易发区和高易发区内灾点频率比最高，符合客观事实。

图3 各指标分级分区面积和泥石流数量相关性统计图Fig.3 Statistical charts of correlation between the area of factor and the number of debris flow points

3.5 基于小流域的信息量法易发性评价

目前泥石流易发性定量评价多采用直接以栅格单元为评价单元，往往会出现单条泥石流流域内存在不同的易发性分区，忽略了泥石流单元的整体性，与实际的环境不相符合，对指导泥石流的精准防控不利。为真实反映泥石流的环境演变，解决单条泥石流流域内分区差异性，文中将基于栅格单元计算得来的泥石流易发性总信息量值以流域为单元进行划分，然后统计每个子流域范围内信息量值的平均值，作为研究区子流域单元信息量值，通过自然间断法将其划分为4 个等级：低易发区[-7.85，-3.97]、中易发区[-3.97，-2.25]、高易发区[-2.25，-0.31]及极高易发区[-0.31，3.39]，得到研究区泥石流易发性分区图，见图4（b）。

图4 泥石流易发性分区图Fig.4 Zoning map of debris-flows susceptibility

4 结果分析与精度评价

4.1 易发性结果分析

研究区泥石流极高易发区和高易发区主要集中在小江流域沿岸的阿旺镇、铜都镇、拖布卡镇一带分布，占全域面积的37.27%，其中著名的蒋家沟、大小白泥沟、大桥河沟以及城市后山4 条沟均位于极高易发区，与历史泥石流灾害点吻合度较高。

极高易发区和高易发区内地质环境特征典型，主要表现在以下几点：（1）位于小江活动断裂附近，构造极为发育，地震频发；（2）工程岩组以较软岩组为主，受小江断裂等影响岩石极为破碎，为泥石流发生提供了物质来源；（3）地形陡峻，相对高差大，地貌以侵蚀剥蚀为主；（4）距离水系越近泥石流易发性越高；（5）主要集中在小江河谷一线，人类工程活动强烈，对地形地貌破坏作用明显。

4.2 易发性精度评价

文中采用深度学习二分类中常用评价模ROC 型曲线进行检验，构建以易发性分区累计百分比为横轴，泥石流灾点数累计百分比为纵轴，曲线下面积（AUC）表示易发性预测成功率，其值越接近1，代表准确率越高[26-28]。本次采用的信息量模型AUC为0.876，准确度较高（图5）。

图5 ROC 曲线Fig.5 The curve of receiver operating characteristic

5 结论

（1）构建了云南东川泥石流易发性定量评价指标体系，对评价指标状态分级进行优化处理后发现：泥石流主要分布在高程660～3 000 m、坡度20°～40°，曲率-1～2，地形集中在微起伏地形，工程岩组主要以较软岩组为主，距河流400 m 内，断层在0～2 km，土地利用类型以林地、草地和耕地为主，而坡向对其影响不大。

（2）东川泥石流极高易发区和高易发区主要集中在小江流域沿岸的阿旺镇、铜都镇、拖布卡镇一带分布，占全域面积的37.27%。

（3）通过信息量模型可以很好建立东川泥石流易发性定量预测模型，ROC 曲线检验模型AUC=0.876，准确度较高，建模结果与历史泥石流灾害点吻合度较高，较好地揭示了研究区泥石流易发性特征，为东川防灾减灾工作提供参考。

文中在评价指标选取和状态分级过程中仍存在一些不足之处。评价指标选取和状态分级合并主观性较强，后续可引入类似ROC 曲线等验证模型进行定量化评价，提高其合理性和科学性。