0.   引言

滑坡是一种常见的地质灾害，边坡不稳定造成的滑坡问题己成为全球性的三大自然灾害之一^[1]，它是指在一定地形、地质条件下，由于各种自然或人为因素破坏了原有的力学平衡条件，使得边坡处于不稳定状态，在自重或其他作用下，沿一定的软弱带作整体的、缓慢的、间歇的甚至是突发的向前向下移动的不良地质现象^[2]，如2010年7月27日四川汉源二蛮山在地震、强降雨等影响下发生滑坡，损失十分惨重。

滑坡稳定性评价是滑坡灾害研究的热点问题，也是难点^[3]。滑坡发生的原因有环境条件和诱发条件，受滑坡体特征、地层岩性、地质构造、地形地貌、地震和降雨等多个因素影响，深入系统地研究主导因素和诱发因素与滑坡的发生和空间分布规律之间的关系，对防灾减灾工作有重要的应用价值。

GIS在边坡灾害方面的应用始于20世纪70年代后期，20世纪80年代，特别是进入20世纪90年代以来其应用日趋广泛。近年来，随着人们对边坡防灾的重视和边坡灾害研究理论与方法的逐步完善，GIS在边坡稳定性分析中的应用快速发展，目前GIS在边坡灾害研究中的应用主要有斜坡信息管理、滑坡与动静态环境因素相关分析、滑坡灾害空间预测系统和斜坡灾害预警系统等4个方面^[4~9]。本文以汉源县二蛮山滑坡体为研究对象，通过野外调查和室内实验，分析了该滑坡的基本特征和形成机理。在此基础上，应用GIS手段进行数据采集与建库、显示与编辑、存储与管理以及空间分析，利用内置的VBA（Visual Basic for Application）开发平台对其进行二次开发、结构定制和功能扩展，开发滑坡稳定性评价模块，从而实现灾害预测预报功能，希望能有效预防灾害发生，提高预防及减少地质灾害的能力，特别对余震、降雨比较频繁的区域有着重要意义。

1.   研究区概况

2010年7月27日，受极端天气和持续强降雨的影响，四川汉源县二蛮山山体突发滑坡，滑坡长1.6 km，高620 m，滑坡体约120×10⁴ m³，造成居住于左侧斜坡上的双合村一组5户20人在滑坡-碎屑流的高速运动过程中被掩埋。二蛮山沟谷原为一高频泥石流沟，堆积于沟谷中段的松散物质于主滑坡结束约30 min后，在重力影响下形成第二次滑坡，滑坡体到达万工集镇，造成部分房屋被掩埋，91户（涉及391人）房屋倒塌，414户1338人房屋受到影响，1500人被迫紧急转移^[10]。本次滑坡具有突发性、隐蔽性强、危害性大等特点，值得关注重视。

二蛮山位于万工集镇背后，地理坐标为北纬29°18′，东经102°42′，总体呈东西向展布；三面环水，流经研究区的有白岩河和大渡河，对二蛮山岩土体侧蚀严重；西南侧为一高频泥石流大沟，滑坡源位于大沟上段右侧斜坡上，万工集镇处于沟口部位的堆积扇上，而双合村一组村民居住地处于大沟中下段的左侧斜坡上（见图 1）。研究区属亚热气候带，干热少雨，暴雨季节（6—9月）降水量占全年的70.6%，且多以夜间阵雨和暴雨为主，是滑坡、泥石流活动的高峰期，而地下水类型主要是松散堆积层孔隙水，其次是基岩裂隙水和碳酸盐岩溶隙水^[11]。本次滑坡在成灾机理、动力特征以及成灾模式等方面比较复杂，据相关文献，二蛮山滑坡分为主滑坡与二次滑坡两部分（见图 2），主滑坡又可分为滑源区、携卷滑塌区、主堆积区，二次滑坡又分为滑源区和堆积区^[12]。

图  1  滑坡体概况

Figure  1.  The situation of landslide body

下载: 全尺寸图片 幻灯片

图  2  滑坡体特征

Figure  2.  The features of landslide body

下载: 全尺寸图片 幻灯片

2.   滑坡机理分析

“7.27”汉源二蛮山滑坡泥石流地质灾害是内外力共同耦合作用的结果，受该地区特殊的地质构造、地形地貌、岩性、地震和降雨等多因素影响，其中地震和降雨起主导控制作用；但在对具体区域进行评价时，往往因客观条件限制，并不需要所有的因素都参与评价，而应根据研究区的实际情况具体选择。

2.1   地形地貌条件

二蛮山山体地处高山峡谷区，属高中山地形，坡向南西，沟谷纵坡坡度平均27°，滑坡区为一险峻高陡斜坡，稳定性差，地形复杂，上陡下缓，区内地形高差悬殊大（相对高差近千米），处于高势能区，使滑坡具有高速启动的位能条件，为滑坡的发生提供了基本环境条件，从而引起滑坡。

2.2   地层岩性

本区主要由二叠系峨眉山玄武岩和梁山组、阳新组构成。玄武岩岩体内柱状节理发育，表层岩体风化严重，塬面破碎；梁山组以钙质、炭质黏土岩为主；阳新组以生物碎屑灰岩、白云岩为主，区内主要是生物碎屑灰岩，灰岩表层风化溶蚀强烈，陡峻且凹凸不平，多发育溶蚀沟槽。滑坡源主要位于沟谷左岸的玄武岩斜坡上，陡峭崎岖，加之“5.12”汶川地震使斜坡产生了一定的变形、破坏，造成土质疏松，地形破碎，为滑坡提供了物质基础，也为地下水渗入提供了良好的通道^[12]。

2.3   降雨的诱发作用

2010年7月23—27日，汉源县出现连续强降雨天气过程，降雨主要集中在24—25日，累计降雨量达到163 mm，尤其7月25日小时降雨量达到50 mm（见图 3），超过短时特大暴雨标准（≥50 mm/h）。在7月23日降雨过程中，滑源区斜坡便出现局部开裂，滑坡量较小；随后在7月25日二蛮山大沟内暴发了特大滑坡泥石流灾害，一部分滑坡体沿着山体快速向下移动，另一部分抛向了山腰左侧村民居住的区域；至7月27日滑源区斜坡出现剧烈变形和崩滑，致使沟槽及两旁碎屑物质缓缓滑向沟口堆积扇上的万工集镇。强降雨造成土壤含水饱和，岩土软化，从而降低了岩层间的抗滑力，增大了下滑力，并对滑源区斜坡坡脚严重掏蚀，导致斜坡失稳，是诱发此次特大山洪泥石流灾害的直接原因。

图  3  滑坡前48 h逐小时和累计降雨量

Figure  3.  Hourly and accumulated rainfall before 48 hours of landslide occuring

下载: 全尺寸图片 幻灯片

3.   数据与方法

3.1   数据源分析

本研究以SRTM（Shuttle Radar Topography Mission，航天飞机雷达地形测绘）为基础数据，此数据带有空间特征和属性特征，包括高程、坡度、坡向等信息，可免费使用且广泛应用于地质研究中，具有现势性好、分辨率高和精度高等特点；同时结合研究区地形图、降雨和地质等数据进行分析。本研究所用的软件平台是ArcGIS 10.1，为美国Esri公司研发产品，具有空间数据采集、存储与管理、查询检索、专题制图、数据输出、空间分析以及二次开发等功能。

3.2   研究方法

本文主要研究思路是对SRTM-DEM数据进行预处理，利用ArcGIS的空间分析功能提取高程、坡向和坡度等派生数据，建立针对研究区的坡度、坡向和高程等评价指标图形库，并建立基础属性数据库。GIS与模糊数学评价方法相结合，并基于ArcGIS平台开发出评价插件，输入滑坡相关的具体预测参数，如高程、最大月（日）降雨量、平均月（日）降水量等，运行模型，计算每一个评价单元的所有指标，得出潜在危险系数；在此基础上，使用ArcGIS的重分类功能对危险系数进行分级，合理划分灾害危险区划范围，进而专题制图绘制出研究区域危险区划评价图，对灾害发生的可能性做出迅速评价；最后与现场实地考察和其他研究结果对比，进行成果检验分析。本文评价体系包括评价单元选取、评价指标的选取及量化、选定评价模型、危险性分析与分级和评价结果制图等内容。

3.2.1   评价单元选取

本文根据滑坡机理和地貌特征，将斜坡共划分为约200个单元（见图 4）。

图  4  评价单元划分图

Figure  4.  A map showing division of the evaluation unit

下载: 全尺寸图片 幻灯片

3.2.2   构建评价指标体系

结合现场调查研究和室内实验，对灾害发育条件、特征及机理进行分析，选取了7种主要成灾因子，即滑坡体特征、地质构造、地震、坡度、坡向（见图 5）、高程（见图 6）及降雨量。但是由于条件限制，地震资料收集难度大，因此本文主要研究降雨诱发型滑坡，深入系统地研究降雨与滑坡发生和空间分布规律之间的关系。

图  5  坡向图

Figure  5.  A map of aspect

下载: 全尺寸图片 幻灯片

图  6  高程图（黑色为阴影效果）

Figure  6.  A map of elevation

下载: 全尺寸图片 幻灯片

3.2.3   评价模型

边坡稳定性分析评价受制于诸多因素，具有耦合性、复杂性、不确定性等特点。本文运用模糊数学综合评判方法加以分析和判别，把这些具有模糊概念的自然属性量化，通过数学函数使其转变为定量计算^[13]。

模型计算公式为：

式中：H——安全系数；F_坡度，F_坡向，F_高程，F_岩，F_雨——单元地形坡度、坡向、高程、单元岩性及单元降雨量危险性等级；K——效应增大修正系数，R₁，R₂，R₃，R₄，R₅——单元地形坡度、坡向、高程、单元岩性及单元降雨量的影响系数（即权重系数），本文中的取值分别为0.225，0.114，0.286，0.175和0.200。

降雨量是本文主要的分析变量。利用ANSYS对二蛮山滑坡坡体进行建模，FLAC3D模拟降雨滑坡过程，得出引起滑坡的降雨临界值，进而确定降雨对滑坡的影响系数（即参数建议值）^[14]。针对本区的特点，本文按高程把研究区划分为0.8～1.0 km，1.0～1.2 km，1.2～1.4 km，1.4～1.6 km和大于1.6 km等5个区域，并结合层次分析法（AHP）将坡度划分为0°~10°、10°~20°、20°~30°、30°~60°和＞60°等5个等级，将坡向划分为平面、SW和S、NW和W、SE和E以及NE和E等5种类型，将岩性划分为土层、灰岩和玄武岩等3种类型，输入不同降雨量进行预测。

3.2.4   评价结果及分析

本文基于ArcGIS平台，利用内置的VBA（Visual Basic for Application）进行二次开发、结构定制和功能扩展，开发滑坡稳定性评价模块（见图 7），从而实现灾害预测预报功能。

图  7  评价插件运行界面

Figure  7.  The running interface of evaluation plugin

下载: 全尺寸图片 幻灯片

输入相应降雨预测参数，可对不同降雨量进行分析。使用ArcGIS重分类功能对危险系数进行分级（见表 1），合理划分灾害危险区划范围，通过专题制图绘制出研究区域危险性评价图，将研究区域划分为低、中和高3个不同危险性等级的区域，并可给出各区域占研究区总面积的比例。

表  1  划分标准与面积比例

Table  1.  The standard and proportion scale

下载: 导出CSV 

| 显示表格

按照以上危险性分区原则和方法，结合已有研究基础，对二蛮山山体进行地质灾害危险区划分，结果如图 8所示。由图 8可以看出，滑体集中在北坡中部，其中高度危险区（红色区域）面积约占7.5%，主要集中在大沟上段右侧斜坡，高程在1600 m左右，与二蛮山滑坡实际灾情基本吻合，需要密切关注留意；中度危险区（黄色区域）面积约占13.7%，大部分处于高程1300 m左右，分布于高度危险区（红色区域）周围，容易受影响形成次级滑坡灾害；低风险区（绿色区域）面积占79.8%。

图  8  降雨量为6 mm时危险性评价图

Figure  8.  A risk assessment map for the 6 mm rainfall

下载: 全尺寸图片 幻灯片

4.   结论与认识

本文以汉源县二蛮山滑坡为例，在不考虑地震的作用下，主要研究降雨对滑坡产生的影响，并在此基础上进行滑坡危险性预测。

安全系数（H）和危险区所占面积比例的分析结果均表明，该区滑坡灾害主要发生在高程大于1.4 km的地段。研究区约21.2%的面积处于30°—50°的坡度地段，随着高度和坡度的增大，发生滑坡的可能性也增加。

采用GIS与模糊数学综合评判方法相结合进行滑坡稳定性评价，其评价结果与实地考察分析和极限平衡法计算结果基本相符，较好地反映了滑坡体实际状态。但是由于滑坡泥石流灾害问题的复杂性，仍然存在着许多有待于进一步完善的地方，例如前期评价因子选取、区间划分、区划结果的准确性、各影响因子参数建议值、不同研究区域评价预测模型选取等问题还有待于更加深入的研究。

本研究可为震后滑坡体稳定性评价提供参考，对做好强余震预报和次生灾害防治工作具有重要意义，尤其对受大地震影响严重、余震频繁、降雨充沛的地区具有参考价值。