特定区域内山体滑坡等多灾害预测

２０１７年５月　水利水电快报　ＥＷＲＨＩ　第３８卷第５期　

文章编号：１００６－００８１（２０１７）０５－－００２７－０４　

特定区域内山体滑坡等多灾害预测　

马雅文　Ｇ．加普林达什维利　Ｃ．Ｖ．韦斯腾。　

（１．长江勘测规划设计研究院，湖北武汉４３００１０；２．格鲁吉亚第比利斯国立大学，　第比利斯，０１７９；３．荷兰特温特大学ＩＴＣ学院，恩斯赫德，７５００ＡＥ）　

摘要：滑坡作为山区主要自然灾害之一，给人类生命和财产带来巨大威胁。在历史数据缺乏的情况　下，以格鲁吉亚姆茨赫塔一姆季阿涅季州为研究范围，模拟了该区域的滑坡及其范围。利用ＩＬＷＩＳ软件空　间多标准评价（ＳＭＣＥ）方法评估特定区域范围内土地脆弱性，选择高脆弱性区域作为滑坡起始点。再采用　Ｆｌｏｗ—Ｒ软件模拟了当触发因子较弱、中等和较强３种情况的滑坡范围。结果表明：在充分了解当地环境　的情况下，经验型的滑坡预测可以较好地模拟出滑坡范围；对滑坡范围的模拟比单纯的土地脆弱性分析能　更好地预测这类自然灾害。从格鲁吉亚环保部提供的已有滑坡地点记录看，预测结果与当地有记录的滑　坡灾害基本符合。　关键词：灾害；滑坡；脆弱性评价；空间多标准评价；区域范围；滑坡范围；模型预测分析　中图法分类号：Ｕ４１６．１６３　文献标志码：Ａ　

１研究背景和意义　

自然灾害是影响许多国家经济、社会和环境的　

原因之一。滑坡作为山区主要的自然灾害之一，给　人类生命和财产带来了巨大威胁。广义的滑坡是指　

一切由于重力作用引起的山坡块体运动，还包含泥　石流、岩崩和雪崩等。　

目前，为了降低滑坡带来的危害，在滑坡发生前　及时搬迁是最理想的方法之一。因此对滑坡的预　

测，有助于减轻和降低其带来的损失。然而，由于缺　乏历史数据，通常难以实现对区域规模滑坡的预测。　

本文研究了在缺乏历史数据的情况下，特定区域范　围内山体滑坡范围的预测和模拟方法。利用ＩＬＷＩＳ　

软件空问多标准评价（ＳＭＣＥ）方法评估区域范围内　

土体的脆弱性，选择高脆弱性区域作为滑坡起始点，　利用Ｆｌｏｗ—Ｒ软件模拟了当触发因子（例如地震、　降雨等）较弱、中等和较强３种情况的山体滑坡　

范围。　

收稿日期：２０１７—０３－２０　作者简介：马雅文，女，长江勘测规划设计研究院，工程师　２数据搜集和分析。　

２．１数字高程模型　

通过数字高程模型数据不仅可以获得特定区域　

地形因素（如坡度、角度、坡度曲率），还可以获得水　文参数（如流向、流路、水流网络等）。该研究对比　

了两种数据源生成的数字高程模型，分别是来自　ＡＳＴＥＲ的全球高程模型（ＧＤＥＭ）和格鲁吉亚方面　提供的２０　ｍ精度等高线。　

滑坡范围预测软件Ｆｌｏｗ—Ｒ对区域范围运行　的最小栅格精度是５０　ｍ；ＡＳＴＥＲ　ＧＤＥＭ数据高程精　

度较低（约为３０　Ｉｎ）；ＡＳＴＥＲ数据内插５０　ｒｎ栅格时　

受原格式影响较大，而２０　ｍ精度等高线较好。所以　该研究中所采用的数字高程模型由２０　ｍ精度等高　

线内插获得，栅格大小为５０　ｍ。　

２．２土地覆盖地图　

土地覆盖地图来自１：５００　０００比例尺地貌图。　

・２７・

　２０１７年５月　水利水电快报　ＥＷＲＨＩ　第３８卷第５期　

经过对地貌图数字化，得到矢量化土地覆盖地图。　对土地覆盖地图的分析发现，由于山区地形特殊，山　

地人类活动相对较少，森林和草地等植被所占比重　

较大，人类活动区域如城市、农田、果园、公路铁路等　所占比重较小。土地覆盖也与地形地貌相符，山区　

海拔较高区域，分布有冰山、草地和灌木等；海拔较　低区分布有森林和人类活动区域。　

２．３地质条件　

地质地图来自于１：５００　０００比例尺纸质地图，　

处理方法与土地覆盖地图类似。由于地质种类众　

多、分布零散，再加上本文针对大范围内的灾害进行　模拟，因此按地质形成年代进行了区分。　

２．４水网和路网　

水网数据来自数字高程模型，利用软件ＩＬＷＩＳ　

处理获得。这里水网并不仅指河流区域，而是泛指　

地形均为沟壑状的区域。泥石流等灾害发生在沟壑　中。本文设定集水范围为０．５　ｋｍ　时，水网形成。　

统计数据显示，陡峭地区的道路极大程度地造　成了滑坡灾害，因此在本研究中，考虑坡度不小于　

２０。的道路。　

２．５　ＭＯＤＩＳ雪域覆盖数据　

由于本文也分析雪崩灾害，因此下载了ＭＯＤＩＳ　

雪域覆盖数据（ＭＯＤ１０Ａ２）。ＭＯＤ１０Ａ２包含每８　ｄ　内最大雪域面积。该研究时下载了目标区域有雪月　

份近５　ａ来的雪域情况。选取每月中云层覆盖最少　的地图，转化为布尔型数据，然后对所有年份和月份　

的雪域地图进行叠加，选取数值较大的作为雪域　范围　

３研究方法　

３．１土地脆弱性评价　

土地脆弱性地图是指定性或定量得出自然灾害　

空问分布可能性的地图，包含对已有或者潜在的滑　坡灾害程度的描述。　

土地脆弱性地图的生成，通常需要目标区域已　

有滑坡的频率和影响因素多年的记录。通常情况　

下，要获取这种记录费时费力，现在世界上有较完整　记录的多为发达国家，如法国和意大利；对于发展中　

・２８．　国家来说，记录通常不全甚至没有。因此，在缺乏历　

史数据的情况下，无法利用统计模型生成土地脆弱　性地图。　另一种方法是基于经验的地图生成方法，具体　

做法是由专家进行打分。空间多标准评价（ＳＭＣＥ）　就是当客观信息不完全时使用的基于专家知识的模　

型。在这种模型里，基础数据（例如坡度、地质、路　网等）方面的专家经验弥补了测量数据的缺乏。为　

了获得专家经验，需要数值化基础数据和合理的结　合方法。　

ＳＭＣＥ方法之一是在ＧＩＳ的环境下，将基础数　据地图（以下称为要素地图）标准化为从０～１，然后　

对每个要素地图赋予一定的权重进行加和，得到数　值化的土地脆弱性地图。　

基于已有文献和已知数据，该研究中分别选出　

４种自然灾害的要素地图，然后赋予标准化的数值。　３．１．１　山体滑坡灾害要素选择和标准化　

到现在为止，还没有统一山体滑坡要素地图的　标准，不同文献里面选择不同的要素，且大部分研究　都是基于研究区域的特性和已有的数据。Ｆ．Ｎａｄｉｍ　

等选择地形、地质、植被覆盖和土地湿度来进行山体　滑坡预测　，Ａ．Ｙａｌｃｉｎ选择地质、天气、土地覆盖、　坡度进行研究　Ｊ。经过对目标区域的了解，选择以　

下５种数据作为山体滑坡的影响要素。　

（１）地形是影响山体滑坡的最重要因素。坡度　

长度、坡度曲率、坡度方向和陡峭程度都影响着滑坡　的程度。统计数据显示，当坡度为２Ｏ。～３Ｏ。时，滑　

坡发生频率最高，标准化值为１，坡度更大和更小时　都依次降低。　

（２）地貌。在植被覆盖丰富的区域，由于植物　根系固定，发生滑坡概率较小，而缺少植被覆盖的区　域容易风化，因此，发生滑坡概率较大。该研究赋予　

铁路、居民点等人为活动较多的区域标准化值为１，　

将湖泊、冰川等区域标准化值为零，其他依据土地覆　盖取０～１之问的值。　

（３）岩层和地质结构也影响滑坡，已有文献表　明风化程度和岩石种类对山体滑坡有影响，然而由　

于不同区域土地和岩性不同，无法将其他区域的经　验用于目标区域。目标区域岩石的时间分布为第四　

纪、新第三纪、早第三纪、白垩纪和侏罗纪。每个纪　

约有１０～３０种岩层。因此，研究时邀请了一位目标　区域地质专家和两位非目标区域地质专家分别打　

分。结果显示，非目标区域地质专家打分相对集中，　

与地层关系较一致；目标区域地质专家打分差异性　马雅文等特定区域内山体滑坡等多灾害预测　

更大，且与地方已有的山体滑坡更一致。该研究集　

合地层、目标区域专家打分和非目标区域专家打分　（其中目标区域专家打分权重更大）来标准化岩层　影响地图。　

（４）与水网的距离也影响了滑坡灾害。该研究　

取０．５　ｋｍ　为最小汇水面积，利用Ｓｔｒａｈｌｅｒ等级分配　水网等级。根据已有滑坡统计发现，一级水网和二　级水网对滑坡影响最大，因此，一、二级水网赋值为　１，其他为零。图１为水网的Ｓｔｒａｈｌｅｒ等级划分。　

图１　Ｓｔｒａｈｌｅｒ等级划分　（５）道路建设对山区地质灾害也有影响，尤其　是在坡度较大的区域。研究中考虑修建在坡度大于　

２０。的道路。　

３．１．２泥石流灾害要素选择和标准化　

由于发生速度快、范围大，泥石流是山区最危险　的自然灾害中之一。沉淀物、降水和坡度３个因素　

对泥石流起关键作用。相比山体滑坡，泥石流对水　的程度更为依赖。在该研究中，取和山体滑坡相同　

的要素预测泥石流土地脆弱性。区别在于对不同的　要素赋予不同的权重。此内容在下一节描述。　３．１．３岩崩灾害要素选择和标准化　

岩崩通常发生在坡度陡峭、植被少、岩石风化程　

度高的区域。因此选择坡度、土地覆盖和地质作为　

岩崩脆弱性的要素。　（１）通常考虑坡度和岩崩关系时，取一个坡度　

阈值作为是否易发生岩崩的分界线。加拿大铁路统　计数据显示，坡度越高，岩崩可能性越大。研究中根　

据此原理标准化坡度对岩崩的影响程度。　（２）土地覆盖也影响岩崩发生的可能性。植被　

覆盖密度越大，发生岩崩的可能性越小，森林区域岩　

崩可能性比草地要小。岩石裸露越少，发生岩崩的　可能性越小，例如冰川、水库发生岩崩可能性为零。　

（３）地质要素对岩石性质和风化程度有影响，　

在研究中考虑了目标区域的断层和岩石的性质来标　

准化地质要素。　

３．１．４雪崩灾害要素选择和标准化　雪崩的要素包括了地形、土地覆盖、道路影响、　ＭＯＤＩＳ雪域覆盖和高程。其中前３个要素用于预　

测土地脆弱性，后两个要素为布尔数值，用于剔除无　降雪和高程不足以发生雪崩的区域。地形因素用于　

得出坡度、坡向、横向面曲率和汇水面积，与土地覆　盖和道路影响要素共同确定雪崩土地脆弱性。　

３．２土地脆弱性地图和滑坡范围预测　

３．２．１滑坡起始点预测　

研究采用ＳＭＣＥ法生成土地脆弱性地图，该法　

是一个经常用于自然资源和灾害管理的ＧＩＳ应用。　其主要原理是先将各要素进行标准化（这在前一节　

已经完成），然后分别赋予权重，根据各要素和相应　

权重组合综合得出结果。　采用层次分析法（Ａｎａｌｙｔｉｃ　ｈｉｅｒａｒｃｈｙ　ｐｒｏｃｅｓｓ）获　

得各要素的权重。由于缺乏历史数据，本文进行了　敏感性分析，以保证权重的合理性。　

采用ＩＬＷＩＳ软件中ＳＭＣＥ模型建模，将各要素　和相应权重组合，得到各类滑坡脆弱性地图。之后　

选取脆弱性高的土地区域作为滑坡起始点。利用　Ｆｌｏｗ—Ｒ软件模拟滑坡范围　３］。Ｆｌｏｗ—Ｒ中参数的　

设置则依据在历史数据充足的区域试验和验证后获　得。为了模拟真实场景，将触发因子分为较弱、中等　

和较强３类，不同触发因子触发的脆弱性阈值不同，　

触发因子较强时，阈值较小，滑坡起始点数量较多。　

３．２．２滑坡范围预测　

预测滑坡范围Ｆｌｏｗ—Ｒ软件是一个以栅格数据　

为基础的基于经验模型。滑坡范围由双重算法决　定：基于地形确定滑坡方向的算法和基于能量转换　

确定滑坡长度的算法　。　

对滑坡方向的算法确定了滑坡从起始点向四周　８个栅格的方向。Ｆｌｏｗ—Ｒ中包含对不同地质灾害　

适宜的不同算法，以多方向算法为例。　当ｔａ　＞０时，　

：　（１）　∑（ｔａ　）　＝１　’。　式中，　，　为滑坡方向（１，２，…，８），　为向第ｉ个栅　格运行的滑动物的量（１，２，…，８），　为中间栅格与　

第ｉ个栅格之间坡度正切值，　为变量。　

滑坡长度的算法取决于滑动物的动能。由于滑　坡总量为正，动能的计算取决于单位栅格内的滑动