滑坡是我国高发、频发的自然灾害之一,属地质灾害链的重要组成部分^[1],作为典型的非线性动力学过程,滑坡的形成机理受到区域地质环境、多相介质融合、多场耦合条件等复杂因素的控制。根据《全国地质灾害防治“十三五”规划》^[2],受到地质地貌条件、极端天气气候、降水时空分布不均等因素的影响,我国滑坡灾害防治工作所面临的形势依然严峻。作为滑坡灾害的研究热点之一,易发性评价可综合反映出滑坡的威胁区域和潜在破坏能力,是灾害危险性评价、风险评价的基础,通过滑坡灾害的易发性分析,可探讨滑坡在不同条件下的形成机理,厘清灾害的易发性区域。

近年来,国内外研究人员基于不同的地质环境条件,构建并使用多样化的模型来分析滑坡灾害的影响因子与易发性的关系:Ionut Sandric利用证据权法,对罗马尼亚国布雷扎镇进行滑坡敏感性分析,认为证据权法具有较高的精度^[3]; Binh Thai Pham采用误差降低剪枝的装袋算法、多重提升算法、旋转森林算法、随机子空间算法和误差降低剪枝算法分别对印度国北阿坎德邦比托拉格尔地区进行滑坡易发性评价,认为基于误差降低剪枝的装袋模型精确度最高^[4]; Mukhiddin Juliev采用SI,FR,CF指标对乌兹别克斯坦国博斯坦利克地区进行滑坡敏感性评价,划分出5个敏感性区域^[5]; 林齐根利用Newmark位移模型,建立了基于物理机理的汶川县地震滑坡易发性评估模型^[6]; 夏辉分别以SVM模型和ANN模型对巫山县进行滑坡易发性评价,认为SVM模型的预测能力略优于ANN模型^[7]; 赵东梅应用maxent模型对哈尼梯田核心区进行了滑坡易发性评价^[8]; 蒋德明结合确定性系数概率模型和敏感性指数,对泸水县的滑坡孕育因素进行分析^[9]。这些单类模型的应用在一定程度上克服了概率统计理论的局限性,对于不同耦合模型在滑坡易发性评价中的探讨较为缺乏。

应用于地质灾害评价的信息量模型,是基于熵的概念来分析因素,其概率^[10]通过信息熵的增减来表示。BP神经网络适用于求解内部运行机制较为复杂的科学问题,通过计算影响地质灾害发生的各因素的信息量,经单因素或加权叠加^[11]后,由神经网络训练带有已知经验的数据集来自动提取“合理的”求解规则,当实际输出与期望值的误差达到目标后,再通过求得的贡献率来量化影响因素,同时避免了专家打分或赋值等人为因素所造成的干扰。信息量法与BP神经网络已在泥石流评价^[12]、生态风险评价^[13]、库区减灾^[14]、煤矿生产^[15]、核物理技术^[16]、地壳活动分析等^[17]领域得到应用,评价模型稳健。湘西州是武陵山区土质滑坡的密集发育区,由于孕灾环境、诱发机制的不确定性与模糊性,迫切需要进行土质滑坡的易发性评价。基于此,本文利用GIS技术与信息量、BP神经网络结合的模型探讨土质滑坡的影响因素和空间分布特征,以期为武陵山区的防灾减灾政策提供科学依据。

研究区位于湖南省西北部,下辖吉首市、泸溪县、凤凰县、花垣县、保靖县、古丈县、永顺县和龙山县。区内气候属亚热带季风性湿润气候,多年平均气温16.2 ℃,多年平均相对湿度78%～80%,水系主要属沅水与澧水,主体河流为沅江及其支流酉水、武水,干流长度大于5 km、流域面积在10 km²以上的河流440余条。研究区属我国第二级阶梯云贵高原的东北部,构造上位于扬子地台东南缘,毗邻华南加里东褶皱区^[18],进入新构造运动期以来,研究区表现为断块间歇性差异隆升,平均抬升速率为0.77 mm/a^[19]。

先通过信息量法将9个因子特定属性段的值,结合ArcGIS的数据处理模块,将属性值转化为影响因子的信息量,再将其转换为数值矩阵后加载至Matlab R2013。参与运算的影响因子作为9个输入层,影响因子的贡献率作为9个输出层,结合公式(3),设置隐含层神经元数为8,输出节点为1,迭代次数为10⁴次,目标误差为10^-7,构建BP神经网络进行训练。最后基于训练结束时BP神经网络的输出结果,求得各影响因子对土质滑坡的贡献比例(图2—3)。

图2 模型构建框架

(1)高程。高程是土质滑坡发生的基础条件,不同的高程对滑坡运动的影响也存在差异。根据研究尺度与土质滑坡点的实际情况,从350 m起,以300 m为步长,根据ArcGIS中的自然断点法进行分级,得到土质滑坡的栅格比例与信息量值。发生土质滑坡的高程多为350～650 m。

(2)坡度。坡度反映了地表单元的陡缓程度,是分析土质滑坡形成条件、易发性的重要基础数据。利用ArcGIS中Surface Analysis的坡度工具Slope提取研究区的坡度数据,最容易发生土质滑坡的坡度为10°～20°。

(3)地形起伏度。地形起伏度是定量描述地貌形态的重要参数,起源于前苏联科学院地理所提出的“地形切割深度”概念,目前已成为划分地貌类型的重要指标。地形起伏度可在宏观上描述地表起伏,由研究区内海拔最高处与最低处之差确定^[22],计算公式为:

H=H_max-H_min(3)

式中:H为地形起伏度; H_max为最高海拔; H_min为最低海拔。

(4)岩石坚硬程度。研究区出露地层自中元古界冷家溪群至新生界第四系,缺失石炭系、新近系、古近系。在一个区域内,可供滑坡形成的固体物质量首先取决于地质因素,软弱或软硬均有分布的岩层易遭到破坏与侵蚀,有利于滑坡的形成^[23]。根据湖南省1:250万地质图和《岩土工程勘察规范》^[24],对研究区的岩石坚硬程度进行划分(表2),在较硬岩间发生土质滑坡的比例最高,次为软岩。

图3 各影响因子分级与信息量统计

(5)沟谷密度。作为描述地表被水道切割破碎程度的重要因素,沟谷密度与地面破碎程度成正比,在地表径流易形成的情况下,土壤侵蚀逐渐加剧并不断诱发滑坡^[25]。沟谷密度的计算公式为:

D_s=(∑L)/A(4)

式中:D_s为沟谷密度; ∑L为研究区内的沟谷长度; A为研究区面积。

经ArcGIS提取沟谷、矢量化、剔除伪沟谷后,得到研究区的沟谷密度。

(6)年降雨量。在降雨期间或者降雨过后,岩土体内的空隙水压力升高,降雨入渗在裂隙处形成局部饱水带^[26],使得潜在滑动面上的有效应力与抗剪强度降低^[27],诱发岩土体失稳与变形^[28]。年降雨量越集中,对土壤黏性、含水量、抗剪强度等性质的影响越显著,年降雨量在1 370～1 420 mm所引发的土质滑坡最多。

(7)地貌类型。研究区地貌类型分为侵蚀溶蚀型中山峰丛洼地地貌(I)、侵蚀溶蚀型低山溶丘谷地地貌(Ⅱ)、河谷侵蚀型堆积地貌(Ⅲ)、侵蚀构造型中低山地貌(Ⅳ)和溶蚀构造型高山台地峡谷地貌(Ⅴ),侵蚀构造型中低山地貌中最易发生土质滑坡。

(8)土壤类型。研究区地带性土壤、耕作土壤均有分布(表3),成土母质复杂。粗骨土分布于沟深坡陡的丘陵地带,土层浅薄,岩石碎屑或砾石含量高; 红壤呈带状分布于研究区东部和北部,范围较广; 酸性紫色土作为初育土壤,母岩疏松,易于崩解。

(9)土壤侵蚀强度。根据目前我国土壤侵蚀强度的划分标准^[29](表4)和区内强烈的岩溶发育程度,水力侵蚀为研究区土壤侵蚀最主要的类型。强度水力侵蚀最易诱发土质滑坡。

表2 研究区岩石坚硬程度

表3 研究区土壤类型

(1)本文基于信息量法与BP神经网络建立了湘西州土质滑坡易发性评价模型,根据ROC检验结果,AUC值达0.76,有73%的土质滑坡点位于高易发区,密度与发生比率随易发性等级的升高而增加,与灾害等级的划分原则相符、与实际情况吻合。通过量化9个影响因子对土质滑坡的贡献率,认为岩石坚硬程度、地形起伏度、年降雨量是土质滑坡的主控因子,地貌类型、坡度、土壤侵蚀强度对土质滑坡的发育产生较大影响。

(2)湘西州土质滑坡表现为残坡积松散物较厚、边坡稳定性差的特征,高易发区占湘西州面积的24.87%,主要分布于地形切割强烈、软硬岩层相间、降水丰沛的西北部和东部地区; 中易发区占湘西州面积的41.79%,分布于湘西州南部以侵蚀溶蚀型峰丛洼地、侵蚀溶蚀型溶丘谷地为主的地貌区; 低易发区占湘西州面积的33.34%,主要分布于地势平坦、土壤侵蚀强度低的中部、西南部地区。

(3)湘西州是湖南省中低山区滑坡的密集分布区,未来的防治重点应侧重于交通干线两侧以及峡谷人口密集区。针对相同类型的滑坡灾害区划研究,未来可利用数据挖掘、机器学习算法的多源观测数据模型,来提高预测、判断和区划的精度,最大程度的管理好滑坡区域。