滑坡敏感性分析

1、西南、贵州等地基于地理信息系统山地滑坡摘要:本文研究的目的是用地理信息系统的方法(GIS)评价中国西南贵州省周 边地区发生滑坡的可能性。底图是扫描这片区域，单幅滑坡绘制1:50万地形图。 在这片文章中，斜坡、岩性、滑坡单体、构造活动、排水分布和年降水量都将作 为独立的影响因子。因此，六个影响因子的地图准备通过各种权威资源收集信息 并且把它们转换到GIS地图中。在定性地图和梯形模糊数学比重(TFNW)方法 的基础上进行敏感性分析。使用预测图概率，将研究区的滑坡敏感性分为四个类 别:即低、中、高和很高。此外，加权程序显示TFNW是一种有效的方法在滑坡 因子权重分析上。关键字：滑坡灾害空间、地理信息

2、系统、滑坡、模糊数、加权方法 绪论2007年，在中国约有25,000件水文地质的混乱事件造成679人丧生，直接 破坏资产和财产估计约3200万美元根据国土资源部的报告(2008) ( HYPERLINK /)%e3%80%82%e5%9c%a8%e8%b4%b5%e5%b7%9e%ef%bc%8c%e4%b8%ad%e5%9b%bd%e7%9a%84%e4%b8%80%e4%b8%aa%ef%bc%8c%e5%9c%a8%e6%97%a0%e6%95%b0%e7%9a%84%e5%9c%b0%e8%b4%a8%e7%81%be%e5%ae%b3%e4%ba%8b%e4%bb%b6%e4%b8%

3、ad%e5%a4%a7%e7%ba%a6%e6%9c%89 /)。在贵州，中国的一个，在无数的地质灾害事件中大约有 150人丧生和许多基础设施被摧毁。建立滑坡敏感性分析对于土地利用规划和基 础设施的布局非常重要。随之而来的需求是用一种逐渐被重视的地理信息系统(GIS)的方法去预测 已经促使许多随机过程模型发展的事件。地理信息系统是基于计算机的系统，具 有高潜力的存档、处理、分析和显示地理参照数据(Aronoff 1989)的能力并且日 益成为地质灾害分析和风险降低的主要工具(Coppock 1995)。GIS作为基本的分析工具，结合一些合理的模型，采用空间数据管理和操作 对于滑坡危险性区划是非

4、常有效的。因此，许多研究者试图在GIS环境下找出 产生山体滑坡和滑坡的灾害图。滑坡灾害图根据不同的目的和现有的资源，在不 同的尺度编译，从地方到国家，采用多种技术，包括立体航空照片的分析、地貌 域的映射、工程地质勘查和历史档案的索引(Guzzetti et al. 2000)。对于这些技术 的结合经常使用(Galli et al.2008)。至于滑坡灾害分析模型，有许多滑坡危险性评价的研究。统计模型之一的多 元回归模型适用于滑坡危险度区划(Dai et al. 2001;Ohlmacher and Davis 2003)。一 些研究应用概率模型(Parise and Jibson2000;Ba

5、eza and Corominas2001;Clerici et al. 2002;Donati and Turrini2002;Zhou et al. 2002; Lee and Choi2003; Lee et al.2004a)。一些研究采用人工神经网络模型(Pistocchi et al. 2002; Lee et al. 2004b; Lee and Evangelista 2006)。一些研究采用岩土模型和模型的安全系数(Refice and Capolongo2002; Zhou et al. 2003;Carro et al. 2003)。此外，定性地图组合是一种非常流行的方法

6、，用它将各种因素组合获得一个 边坡失稳地图(Stevenson 1977;Saha et al. 2002;Sarkar and Kanungo2004)。这种方 法的区别在于山体滑坡并不需要，因为权重的分配是根据地貌学家的丰富经验 (Galang2004)。尽管如此，这里仍然有一个名为加权线性组合(WLC的)模型，它被称为 是最知名的和最常用的多准则GIS方法(Eastman 1999; Jiang and Eastman 2000;S 5ener et al. 2006;Akgu n andBulut 2007)。本质上，WLC的模型是一种结合各种敏感因子提高滑坡敏感 性分析的方法。这项研

7、究的目的是运用定性地图组合模式在贵州地区采用 GIS 分析滑坡敏感性区划图。每个因素引起的滑坡权重是有模糊简单加权模型实现。 本文的其余部分组织如下。下一节介绍了研究区的基本信息，紧随其后的是各种 偶然因素和定性地图组合模型。这部分“加权方法”描述了一个模糊简单加权模 型。本文以讨论和结论结束。研究区贵州省位于经度 103。36-109。35，纬度24 0 37 -29 0 13 ，覆盖面积176,167km2，海拔高度从140-2,900m(图1)，亚热带型气候。作为山区，该地 区的地质结构非常复杂。它遭受大量的山体活动，如滑坡、泥石流、塌陷和沉降， 主要由从三月到八月的暴雨引起。由于雨水增

8、加了空隙水压力，土体强度降低， 最后导致滑坡。因此，每年在这个地区降雨诱发的地质灾害频繁发生，严重威胁 人民生命财产安全并且妨碍经济发展。XlUiAN|ALegendAltifuda1 M 25W-li0M J 3*H-I5O 3 3544-2WC图1研究区的分布与海拔高度偶然因素和定性地图组合模型滑坡灾害是发生在一段指定时间内的一个特定区域的潜在破坏性现象 (Varnes1984)。这个定义既包括了空间又包括了时间的尺度。这种危险性的区划 是指将陆地表面划分区域，并根据实际或潜在的土石流或是斜坡上的一些大块物 质将这些区域分类(Varnes1984)。滑坡危险性区划的最终地图显示了相对滑坡灾

9、 害的空间分布也常被称为敏感性地图(Brabb1984; Hansen 1984)。定性地图组合原理的描述为：可能发生山体滑坡作为一个确定的因果关系的 若干因素的结果(Varnes1978)。每个地图的因子进行加权然后分成子类，然后将 权重分配给每个子类。该地区上覆盖不同的加权因子和子类比重允许为每一个相 对不稳定的地图得分求和，这些得分可以将地图分为不同的不稳定类。偶然因素为了评估潜在的滑坡，需要识别和分析偶然因素。在这篇文章中使用了一下 参数:斜坡、岩性、堆积密度、构造活动、年降水量和分布。一个常用的因子经 常用在边坡失稳的分析中即植被，但不会用在这个模型中主要是因为最常见的植 被数据较少

10、。数字高程模型（DEM ）和岩层，地质构造和排水管网在研究区域获得了全 国基金会地质信息系统（NFGIS-M50DB）所提供的基本地图和空间数据库。滑坡 资料是由土地资源部提供的，而降水资料是由贵州省气象局提供。这是一幅1:50 万的地图。由于这些层由连续数据组成，它们需要为离散子类从新分类。这些子类是有 些主观，但都在该地区以往研究的基础之上，如表1所示。子类分类的主观性始 终存在并将继续成为边坡失稳区划项目中最主观的方面。然而，在一些不同的研 究领域使用一些一致性可能有助于减少主观性的影响。基本因素和触发因素都会 使边坡产生灾害预警图，如表1.这篇文章中只使用了基本因素生成斜坡敏感性的 区

11、划图。表1滑坡危险性评价因子因素子类1234基本因 素斜坡（0）缓25中 25-35陡 35-45很陡45岩性G1G2G3G4构造活动非常弱弱中强渗透强度弱中强很强滑坡范围（#/100km2）10年降水量（mm/y）1,0001,000-1,4001,400-2,0002,000触发因 素五天的总降水量（mm）200第二天降水预测（mm）150斜坡最初，数字高程模型（DEM）对这个区域的研究是建立在使用ArcGIS9.0软 件的地形图。在1:50万的地形图中，等高线被用来生成一个数字高程模型。坡 角的计算采用DEM，斜坡输出的度数分为四个子类（表1）。岩性岩性的资料来自1:50万的NFGIS-

12、M50DB图，它是通过实地调查得到的。四 大类型的岩石和沉积物存在：（1）碳酸盐岩占区域的70%，包括层状和大块的石 灰石，石英质的灰岩和泥灰岩；（2）岩浆岩占区域面积的7%；（3）变质岩占区 域面积的11%，包括砂岩，粉砂岩，凝灰岩；（4）海陆沉积岩占区域面积的11%， 由粘土，粉质粘土，细砂，粗砂，砾石和鹅卵石组成。根据造成滑坡的可能性， 将岩石分为四个子类（表1）并且易滑岩性图如图2所示。图2易滑岩性图（1）G1:组成松散堆积物的粘土，砾石土，卵石土或砂卵石土；粉砂岩，页 岩，泥灰岩，凝灰岩，千枚岩或粘土岩；变质砂岩，粉砂岩或变质凝灰岩。（2）G2:中厚层的砂岩，页岩，泥灰岩，凝灰岩，千

13、枚岩或粘土岩；变质 砂岩，粉砂岩或变质凝灰岩。（3）G3：厚层状砂岩，含砾泥岩，页岩，泥灰岩，凝灰岩，泥页岩。（4）G4:特厚层砂岩，砾岩，石灰岩，白云岩，硅质岩，磷矿体。构造活动岩石结构复杂，通常形成断裂构造，断层破碎或者挤压褶皱，它们容易导致 滑坡事件。山体爆破或地震触发滑坡的可能性由距离的拉长而减小。根据构造活 动和断层的距离贵州这片区域可以分为四个子类（表1）。水系分布滑坡位置和分水岭之间有一个密切的关系，集水区的面积与土壤可以吸收和渗透的水量有 关。更大的区域与更强的渗透性和更高的滑坡几率有关，流域的长度暗示了流域面积的大小， 地下水的集中程度和沉淀物的多少。水的功能，如湖泊，池塘和

14、大江大河已在这一研究中获得。贵州省的水系分布图从1:50万的NFGIS-M50DB图中获得（图3）,根据水系的大小和 距离，贵州省这区域的河流可分为四个子类（表1）。*5,uWwhwn图3水系分布图滑坡点从2002年3月到2004年9月的山体滑坡清单图的研究领域由贵州省国土资源部编制，在这 片研究区中的最新的滑坡清单是2006年9月通过实地调查建立的。因此，精致的滑坡清单 地图如图4。新的滑坡点显示为2,392处，包括325处泥石流和130处崩塌。根据滑坡密度 把它分为四个子类（表1）。图4滑坡点年降水量降水量的空间分布在确定交叉区的降雨量和判断危害的可能性方面是一个重要的因素。历史 数据表明

15、，最多最严重的地质灾害几乎都发生在雨季，每年的3月到八月在这个地方。贵州 省的年降雨量分布图如图5。根据年降雨量的降雨强度可分为四个子类（表1）。图5年降雨量分布图QMC模型这些分析在光栅环境中进行，因为网格在功能分析上既简单又易于使用，如叠加。因此，所 有的矢量数据将被转换为栅格地图。一旦因素被确定，只要对各种因素的影响有很好的了解，地图定性结合的方法非常简单并且 可以轻松的应用到任何地方。此方法的测绘单位将是网格单元。这种模式的基本步骤是：（1）各因素图分成离散的子类（已做）。（2）各子因素的权重和总体因素值。加权处理过程已经在“加权方法”中详细描述过了。一旦权重被确定，基于子类的权重并根

16、 据子类划分网格，这样数据层被重新归类到新的层次。是一个关于斜坡子类权重的重要概念。（3）一个不稳定的数字地图只能使用各种加权方法覆盖（Van Westen1994）,最后的成果图会 对每一个栅格，每一个得分在分配的权重基础上有一个范围。使用栅格计算器，重新分类数据层划分为算术交叉重叠和原始数据的不稳定性，S，每一个 栅格计算值。对于每一个格子，k,是每一个子类因子的特殊组合，即不稳定评分，S j在该 单元格中的计算公式为（1）。S j = 26 wk .Wj（1）i =1这里i是一个因子，ij是因子i的子类。式的结果S j如果接近1表示很有可能发生滑坡然而如果数字接近于。则表示这种可 能性很

17、小。这样，概率图有一个数值的连续规模获得。将地图分为敏感性的子类还是有必要 的。Ayalew and Yamagishi(2005已在考虑根据分位数，天然缺口，等间隔和标准差将敏感性分为 四个系统。他们已经发现以分位数为基础的分类系统有一个缺点，就是它将相差甚远的值分 到同一个类中。天然缺口在数据值存在大跳跃时分类效果更好，使用等间隔的方法有未必是 有效的因为它强调一类相对于其他类。作为比较的结果，他们选用了标准差的方法，因为它 使用的生成类的方法。在这篇文章中，用标准差方法分成四个敏感性类，它们是很高，高，中，低。根据这一敏感 性地图(图6)，总面积的9.24%发现是低敏感性区域，中和高敏感

18、性区域分别占34.78%和 41.53%，非常高易发区占研究区的14.45%。图6贵州省滑坡敏感性地图加权方法在QMC和WLC模型中，一个首要的问题是分给每个因素的权重分别是多少。有许多技术， 例如统计指数法(Wi)(Van Westen1997)权重因子(WF) (C , evik and Topal2003;O ztekin and Topal2005)还有层次分析法(AHP) (Saaty 1983; Saaty and Vargas 2001; Ayalew et al. 2005;Sener et al. 2006; Akgun and Bulut 2007)去找权重。这项研究提出了

19、梯形模糊数的加权方法分配给 每个因素和子类权重。模糊集理论在宇宙坐标系中讨论X，一个X的模糊子集A的定义是一个隶属函数旦A(x)中的每个兀素x映射在X中的一个区间0,1的实数。功能函数目(X)标志着隶属函数A中的x。当旦(x)AA值较大时，x的每个值在A中较大(Keufmann and Gupta 1991)。一个模糊集A = （a,b,c,d）在实数R上，avbvcvd,被称为梯形模糊数（Keufmann and Gupta1991）如果它的隶属函数如（2）所示。x - a7,a x b, b 一 a（2）1, b x c,X - c7,c x d, d - c0, x aord 45)A4

20、0.2430.044岩性B0.161G1B0.4190.067G2B20.3370.054G3B30.2070.033G4B40.0370.007构造活动C0.182非常弱q0.0690.013弱C20.1710.031中C30.3390.062强C40.4210.076渗透强度D0.040弱D0.0420.002中D20.1980.008强D30.2840.011很强D40.4760.019滑坡范围 (#/100km2)E0.23510E40.2950.116年降水量(mm/y)F0.2002,000F40.3700.074基于上述概念模型的TFNW的程序如下：（1）确定决策者的重要性（或可

21、靠度），假设有一个由k个专家组成的委员会，他们负责决 定滑坡危害因素或子类的权重，如表1所示。专家的重要性（或可靠度）是R（t=1,2,.,k）, 当R七e （0,1）并且Z气=1。每个专家的重要性（或可靠度）为w, （t=1,2, .,k），被最终决策者所分配。因此，R,的重要性由式（3）定义。如果R1=R2=Rk，专家小组被称 为同质小组，否则，专家小组被称为异构（非均质）小组。R =-（t=1,2,.,k）（3）*寸乙wti=1（2）介绍语言加权变量（表2）给专家评估滑坡危害因素或子类（表1）的重要性并且计算单个因素或子类的模糊加权综合。例如，让w广（a并，七七,d/ ,（j=1,2,3

22、,4）,（t=1,2,.,k）,是专家t赋予斜坡（因素A）子类A j （表3）的语言权重。A j的模糊加权总量wj由k个 专家组成的委员会评估，可由式（4）定义。w =(a ,b ,c ,d )=( R aj j j j jt jtt=1&b ,2rc ,&d )t jtt jtt jtt =1t = 1t = 1(4)（3）为了减少个别滑坡危险因素或子类的模糊权重的模糊性，wj的模糊性记为d（wj）,因此， 用式子（5）表示。d(w )=(a +b +c +d )/4(j=1,2,3,4)(5)j j j j j（4）计算规范化的权重和建立权重向量，子类A j的规范化权重记为w广用式子（6）

23、表示。d （w）w .=（j=1,2,3,4）（6）J 尸 d （w ）jj=1如果尸w. =1，则A j子类的权重向量可以表示为式子（7）。j=1w=（w1 ,w2 ,w3 ,w4）（7）例子在这一部分，A j子类的权重向量，作为一个例子，计算过程按上述程序。（1）由三个专家组成的委员会，D1,D2和D3是因素A （斜坡）子类给出的权重，四个选 择属性可以考虑：（a）缓（A）；（b）中（A2）；（c）陡（A3）和（d）非常陡（A4）。所 有的专家都认为是同等重要的，即R1=R2=R3=1/3，因此，这个委员会是一个同质群体。（2）使用语言变量和它们各自的权重模糊数（表2）来评估每个子类A j

24、 （表3）的重要权 重。根据评估值（表4），计算每个子类的模糊权重（表4）用式子（4）。表4子类A j的模糊权重子类专家模糊加权总汇D1D2D3缓（A）VL(0,0,0,3)VL(0,0,0,3)L(0,3,3,5)(0,1,1,11/3)中（A2）H(5,7,7,10)H(5,7,7,10)VH(7,10,10,10)(17/3,8,8,10)陡(A3)VH(7,10,10,10)H(5,7,7,10)VH(7,10,10,10)(19/3,9,9,10)非常陡（A4）M(3,5,5,7)M(3,5,5,7)H(5,7,7,10)(11/3,17/3,17/3,8)（3）计算聚合模糊权重值的模糊价值（表5）使用式子（5）和各子类的普遍价值（表5）用式子（6）。这些是权重向量所产生的权重值w=（w,w2，w3, w4 ）=（0.060,0.334,0.363,0.243）。第四类比第三类分配的权重要低，这是由于这个事实，无论是崩积层或风化粘土可以在450 的斜坡上站立。子类的普遍权重B .，C j，D j，E j和Fj可以通过重复计算第（1）到（3）步得到。同理， 因素的普遍权重A,B,C,D,E和F可以通过同样的步骤得到，全部的普遍权重值如表3所示。如式子（1）所示，在网格k中不稳定评分S j相当于所有子类所