《地质灾害区划与分区评价与衡量》由会员分享，可在线阅读，更多相关《地质灾害区划与分区评价与衡量（17页珍藏版）》请在金锄头文库上搜索。
1、word第五章 地质灾害区划与分区评价第一节 地质灾害易发区划分与分区评价一、评价思路与方法地质灾害发育现状是其易发性的客观反映，要想准确的进展地质灾害易发性分区，必须依赖遥感解译和野外实际调查工作，本着这一思路，此次调查十分重视对根底地质元素的搜集与分析，野外工作完毕时已形本钱区地质灾害易发性根本轮廓，即初步的定性分区结果；同时考虑到地质环境条件的复杂性，通过对影响地质灾害发育的诸多因素分析，采用半定量方法进展分区计算，作为对定性评价的补充，最后综合两种结果，形本钱区地质灾害的易发性分区图。由于地质灾害易发性的评价结果受到多种因素的影响，而这些因素本身存在着不确定性、模糊性以与各因素之间相互
2、作用的复杂性；如何将复杂的地质因素尽可能的定量化，使分析和评价结果最大限度的符合客观实际情况，是地质工作者广为探讨的问题。层次分析法简称AHP是一种定性和定量相结合的系统化、层次化的分析方法，它不仅适用于存在不确定性和主观信息的情况，还允许以符合逻辑的方式运用经验、洞察力和直觉。由于它在处理复杂的决策问题上的实用性和有效性，因此在各学科模糊综合评判中得以广泛应用。地理信息系统在最近的30多年得到了惊人的开展，广泛应用于资源调查、环境评估、灾害预测等众多领域，借助GIS系统可以完成制图、数字地形分析、空间决策支持、空间分析统计等任务，在GIS平台上进展易发性区划可以在一定程度上防止传统区划工作量
3、大、工作强度大、工作精度不高以与主观影响大的不足。本次工作拟采用基于层次分析法和GIS空间分析统计方法相结合的工作方法对区地质灾害易发性进展评价和区划，主要技术路线和方法如下：1、确定评价单元和评价因子，利用层次分析法确定各因子和各要素的权值。2、对各评价因子指标进展量化，并采用归一化数值变换方法统一量纲。3、在评价指标权值确定和数据归一化的根底上，利用GIS系统的空间分析功能进展数据的空间叠加与统计。4、经统计分析确定易发性区划的分界点，将评价结果分成不同等级。5、在GIS分析成图的根底上综合考虑各种因素，进展修改完善，最终编制工作区地质灾害易发性区划图。二、理论与方法层次分析法是美国运筹学
4、家T.L. Saaty于20世纪70年代中期提出的一种定性与定量相结合的多准如此决策的系统分析方法，其根本原理是把复杂系统分解成目标、准如此、方案等层次，在此根底上进展定性和定量分析的决策。它把人的决策思维过程层次化、数量化、模型化，并用数学手段为分析、决策提供定量的依据，是一种对非定量事件进展定量分析的有效方法，特别是在目标因素结构复杂且缺少必要的数据情况下，需要将决策者的经验判断定量化时该法非常实用。该方法适用于多准如此、多目标或无结构特征的复杂问题的决策分析，它按照各因子相互之间的在支配关系，建立层次结构模型，通过因子的两两比拟，建立判断矩阵，进展层析排序，确定各因子的相对重要性。一层次
5、分析法的主要步骤运用层次分析法建模，大体上可按下面4个步骤进展： 1、建立层次结构模型：综合分统所析系涉与的目标、围、准如此，约束条件，确定综合评价体系中各因索之间的关系，并根据指标的隶属关系进展上下分层排列，形成危险性综合评价体系层次结构。2、将问题中的各个要素划分为不同的层次结构,以框架结构说明各层次之间的从属关系。3、构造判断矩阵：层次分析法的根底主要靠对各层次两两元素相互重要程度差异给出判断，并将这些判断用数值表示出来，形成判断矩阵，即两两比拟矩阵。该矩阵是层次分析法的出发点，也是整个权值确定过程中的关键。在构造判断矩阵的过程中，可以通过向专家发放问卷进展因子间的相对重要程度打分,并用
6、矩阵表示打分结果;4、检验层次分析结果如有误差,需对判断矩阵的元素取值通过专家意见进展调整,从新运算。二评价指标判断矩阵的构建信息是系统分析最根底的数据。任何系统分析都要掌握一定的信息才能进展。层次分析法也需要有相应的信息作为分析的根底，其信息主要来源于人们对不同层次中各个因素之间的相对重要性所做出的判断。通过引入适当的判断标度将这些判断用数值的形式表示出来构成判断知阵。以便比拟本层次各因素与某一因素之间的相对重要性。设B层次中的元素Bi,Bz,Bs,.从与上一层次A中的元素A有关系，如此可以通过判断知阵表示出来，如表5-1-1所示。表5-1-1 判断矩阵示意表AB1B2B3BnB1b11b1
7、2b13b1nB2b21b22b23b2nBnbn1Bn2Bn3bnn表中bij表示对危险性综合评价而言Bi和Bj相对重要程度的数值表示。bij的取值是根据表5-1-2的T.Satty 1-9标度含义来确定。表5-1-2 判断矩阵元需Aij的1-9度标度法标度含 义1两个因子相比拟，两者同样重要3两个因子相比拟，其中一个比另一个稍微重要5两个因子相比拟，其中一个相对另一个来说比拟重要7两个因子相比拟，其中一个相对另一个来说非常重要9两个因袭相比拟，其中一个相对另一个来说极其重要2,4,6,8介于上面两个相邻判断值的中间倒数假如i与j相比拟的判断值为bij,如此j与i比拟的判断值为1/bij三特
8、征值与特征向量计算根据判断矩阵，利用线性代数知识，准确的求出T的最大特征根所对应的特征向量。所求特征向量即为各评价因素的重要性排序，经归一化后即为同一层次相应因素对于上一层次某因素相对重要性的排序权值。本次评价采用和积法进展求解，具体步骤如下：1、将判断矩阵每一列归一化：2、每一列经正规化的判断矩阵按行相加：3、对向量做正规化处理，依次所得到的即为所求特征向量。4、计算判断矩阵的最大特征根max，式中：TAi表示向量TA的第i个元素。四一致性检验为防止其他因素对判断矩阵的干扰，在实际应用中要求判断矩阵满足大体上的一致性，需进展一致性检验。只有通过检验，才能说明判断矩阵在逻辑上是合理的，才能继续
9、对结果进展分析。对判断矩阵进展一致性检验，计算公式：CR=CI/RI 1式中，CR(consistency ratio)为一致性比例。当CR0.10时，认为判断矩阵的一致性是可以承受的，否如此应对判断矩阵作适当修正。CI(consistency index)为一致性指标，按下式计算：CI=(max -n)/(n -1) 2式中：max判断矩阵的最大特征根；n成比照拟因子的个数；RI(random index)随机一致性指标，可查表确定，如表5-1-3所示。表5-1-3 平均随机一致性指数I阶数123456789RI00当R时，就认为判断矩阵具有满意的一致性，否如此就需要重新调整，直到具有满意的
10、一致性为止。三、评价因子选取与敏感性分析影响地质灾害形成的自然因素众多,历史地震地质灾害发生的数量、分布围、活动规模都直接反映了地层岩性、地形地貌、现存新老滑坡以与有关地震动力环境对地震诱发地质灾害的控制作用;此外,土地利用、地下水、植物条件等因素也对震后地质灾害形成起到一定程度的影响。本文采用统计学方法，对研究区山地灾害点与各因子的每个属性进展相对频率组合的定量计算方法，综合市州区震后地质灾害发育情况，本次危险性区划分析中选用了10个影响因子，主要包括：地质构造、地形坡度、海拔高程与水系发育情况等。一地质构造地质构造因素对地质灾害点的发育控制作用十清楚显，在区域地质构造比拟复杂，褶皱比拟强烈
11、，新构造运动比拟活动的地区，地质灾害比拟发育。其影响主要表现在：地质构造决定了地貌形态的分布，对地质灾害发育的临空条件起到间接的控制作用；地质构造带岩石破碎、风化严重，使得边坡的连续性和完整性受到破坏，是地下水最丰富和活动的地区，降低了岩体的抗剪强度；在构造应力作用下，岩体节理、裂隙发育，为崩塌发育提供了条件；活动断层造成地表破裂，岩层结构发生破坏，非活动断层作为地震波的反射界面，可能导致岩体的拉力破坏；断裂构造控制着水系的发育和人类工程活动的分布，对地质灾害的威胁对象起到间接的控制作用研究中，通过GIS软件缓冲区分析和数据统计功能，对研究区灾害点与断裂距离分布关系做了统计：首先，对研究区的断
12、裂做距离缓冲处理，分别得到0-500，500-1000，1000-2000与大于2000米四个缓冲区；然后利用GIS统计功能，对每个缓冲区的灾害数量、缓冲区面积进展统计，计算每个缓冲区灾害点密度。详细数据如表5-1-4所示，灾害点与断裂的分布关系和敏感性关系，如图5-1-1所示。表5-1-4 研究区地质灾害点与断裂距离统计关系断裂缓冲距离(m)2000缓冲区面积(km2)291.83 面积百分比(%)12.40 10.55 22.44 54.61 灾害点数量342842 100灾害点密度(个/km2)0.12 0.11 0.08 0.08 图5-1-1 断裂构造与地质灾害分布关系图二地质灾害频
13、率对市州区按1km*1km网格进展离散，形成2534个空间离散网格，依据灾害点在网格上的分布，进展基于GIS的统计。该计算包括单元面积上灾害发生的频率与地质灾害面积模数比。地质灾害频率比：设第(i,j)单元灾害频率为f(i,j)，单元面积为S(i,j)，单元灾害的频率密度为f(i,j)，整个研究区面积为S，灾害总数为f，总频率密度为f，如此：第(i,j)单元格灾害频率比为： Rf(i,j)=f(i,j)/f其中，f(i,j)=f(i,j)/S(i,j) ；f=f/S。地质灾害面积模数比：设第(i,j)单元灾害体分布面积为Ss(i,j)，单元面积为S(i,j)，单元灾害的面积模数为s(i,j)，
14、整个研究区面积为S，灾害点总面积为s，总面积模数为s，如此：第(i,j)单元格灾害面积模数为：RS(i,j)=s(i,j)/s其中，s(i,j)=s(i,j)/S(i,j) ； s=s/S。经计算可得，单个像元上最大出现地质灾害的频率为8。对空间灾害点的频率分布进展归一，可形成图5-1-2灾害发生频率归一化分布图与5-1-3地质灾害面积模数比归一化图。图5-1-2 地质灾害频率比归一化图图5-1-3地质灾害面积模数比归一化图三坡度与坡度变率利用工作区1：5万DEM数据提取坡度数据。根据前文中的分析，由于工作区滑坡、崩塌灾害主要分布于1060之间的斜坡，10以下斜坡根本不发生滑坡、崩塌等灾害，因
15、此本次评价将60以上斜坡的易发程度定义为1，10以下易发程度定义为0，将坡度数据进展01之间的线性归一化，得到坡度归一化结果图。坡度变化率是对地形根本因子坡度变化情况进展量化的指标，由于斜坡拉应力区的分布与斜坡坡度呈正相关联系，因此随着斜坡坡度变化率增大的斜坡坡脚地带形成的最大剪应力也不断增大，斜坡也就愈容易产生变形破坏。本次通过DEM对全区坡度变化率数据进展提取，然后进展01之间归一化处理之后参与评价。图5-1-4 坡度归一化图图5-1-5 坡度变率归一化四坡向与坡形变率坡形可以利用地表的曲率进展描述和量化，直线形和凸型斜坡在曲率上的表现是曲率0，凹型坡和阶梯型坡的曲率0，因此，可利用ArcGIS平台从DEM数据中提取调查区地表曲率信息，(平面曲率
