gis技术在地质灾害 gis在地质中的应用

地质信息系统技术

一、内容概述

为望城等地区用户提供了全套网页设计制作服务，及望城网站建设行业解决方案。主营业务为做网站、成都网站建设、望城网站设计，以传统方式定制建设网站，并提供域名空间备案等一条龙服务，秉承以专业、用心的态度为用户提供真诚的服务。我们深信只要达到每一位用户的要求，就会得到认可，从而选择与我们长期合作。这样，我们也可以走得更远！

地质信息系统（GIS），产生于20 世纪60 年代。它随着人们对自然资源和环境的规划管理工作的需要以及计算机制图技术的应用而诞生，是一种对大批量空间数据采集、存储、管理、检索、处理和综合分析并以多种形式输出结果的计算机系统。1965 年，W.L.Garrison首先提出了“地质信息系统”这一术语，开创了这一新技术的发展史。此后，美国、加拿大、英国、澳大利亚等国均投入了大量人力、物力和财力，并逐步确立了他们在这一领域里的国际领先地位（黄润秋，2001）。

二、应用范围及应用实例

1.GIS技术在地质灾害信息系统中的应用

随着人口的急剧增长，经济的迅速发展和自然资源的大量消耗，不仅生态环境恶化，而且导致自然灾害（包括地质灾害）频繁发生。美国、印度等国是世界上地质灾害较为严重的国家，地质灾害具有类型多、分布广和成灾强度高的特点。这些地质灾害大部分发生在承灾能力较低的地区，给当地的经济和社会稳定构成了严重的威胁。地质灾害是地质环境质量低劣的表现，它的频发不仅反映了自然地质环境的脆弱性，而且反映了人类工程经济活动与地质环境间矛盾的激化。要使人类工程经济活动与地质环境之间保持较为协调的关系，就必须对地质环境进行评价，以了解不同经济发展过程中区域地质环境的基本态势和变化趋势，为环境管理和城市规划等提供依据，但传统技术手段已不能完全应付迅速反应的地质灾害。地质信息系统作为当前高科技发展的产物，集图形、图像与属性数据管理、处理、分析、输入输出等功能为一体，应是当前地质环境评价与地质灾害预测的强有力工具（赵金平等，2004）。

GIS 技术的产生是计算机技术和信息化发展的共同产物。是管理和研究空间数据的技术系统。可以迅速地获取满足应用需要的信息，能以地图、图形或数据的形式表示处理的结果（曹修定等，2007）。国外尤其是发达国家在GIS应用与地质灾害研究方面已做了很多工作。从20世纪60年代至今，GIS技术的应用也从数据管理、多源数据集数字化输入和绘图输出，到DEM或DTM模型的使用，到GIS结合灾害评价模型的扩展分析，到GIS与决策支持系统（DSS）的集成，到网络GIS，逐步发展深入应用（黄润秋，2001）。

印度Roorkee大学地球科学系的R.P.Gupta和B.C.Joshi（1990）用GIS方法对喜马拉雅山麓的Ramganga Catchment地区进行滑坡灾害危险性分带。该项研究基于多源数据集，如航空像片、MSS磁带数据、MSS图像、假彩色合成图像及各种野外数据，包括地质、构造、地形、土地利用及滑坡分布。以上数据需要进行数字、图像等处理，然后解译绘制出专题平面图，包括地质图（岩性与构造）、滑坡分布图、土地利用图等。这些图件经数字化及有关数据都存储在GIS系统中，找出与滑坡灾害评价相关的因素，如滑坡活动与岩性的关系，滑坡活动与土地利用的关系，不同斜坡类型的滑坡分布情况，滑坡分布与主要断裂带的距离关系。经过统计及经验分析，引入一个滑坡危险系数（LNRF）。LNRF值越大，表示该地滑坡灾害发生的危险性越高。并且对LNRF的3个危险级别分别赋予0、1、2三个权重。考虑到滑坡的发生是多个因素综合作用的结果，故调用GIS的叠加分类模型，将各因素的权重叠加，得到综合图件，图上反映的是每个地区的权重总和。根据给定标准，即可在这张图上勾绘出滑坡灾害危险性分区图。

荷兰ITC的C.J.Van Westen和哥伦比亚IGAC的J.B.Alzate Bonilla（1990）基于GIS对山区地质灾害进行分析。他们在数据采集、整理方面做了大量工作，建立了一套完整的数据库。在此基础上，开发出了分析评价模型，如斜坡稳定性分析模型，其主要功能是计算斜坡稳定的安全系数。另外，两位学者还利用GIS所生成的数字高程模型（DEM），开发出了一部山区落石滚落速率计算模型，并据此绘出了研究区内落石速率分区图（黄润秋，2001）。

美国科罗拉多州立大学Mario Mejia-Navarro和Ellen E.Wohl（1994）在哥伦比亚的麦德林地区，用GIS进行地质灾害和风险评估（姜作勤，2008）。利用GIS对麦德林地区地质灾害进行了分析和研究，重点考虑了基岩和地表地质条件、构造地质条件、气候、地形、地貌单元及其形成作用、土地利用和水文条件等因素。根据各因素的组成成分和灾害之间的对应关系，把每一种因素细分为不同范畴等级，借助于GIS软件（GRASS）的空间信息存储、缓冲区分析、DEM模型及叠加分析等功能，对有关滑坡、洪水和河岸侵蚀等灾害倾向地区进行了灾害分析，并对某一具体事件各构成因素的脆弱性进行评价。

同样是美国科罗拉多州立大学Mario Mejia-Navarro博士后等人（1996）将GIS技术与决策支持系统（DSS）结合，利用GIS（主要是地质资源分析系统GRASS软件）及工程数学模型建立了自然灾害及风险评估的决策支持系统并应用在科罗拉多州的Glenwood Springs地区（姜作勤等，2001）。应用GIS建立指标数据库，并建立基于GIS的多个控制变量的权重关系式。对泥石流、洪水、地面沉降、由风引起的火灾等灾种进行了灾害敏感性分析、脆弱性分析及风险评估，辅助政府部门做出决策。

美国地质调查局（USGS）已把加强城市地质灾害研究列为21世纪初的重要工作，借助GIS编制美国主要城市地区多种灾害的数字化图件，这种做法与西欧国家的城市地质工作的总趋势一致。其中，美国科罗拉多州格伦伍德斯普林市的城市地质灾害评价项目最具代表性。由于该市位于山区河谷地区，崩滑流地质灾害制约着城市的发展，为此，城市规划部门委托科罗拉多州立大学，开展了GIS地质灾害易损性和风险评价编图研究，最终按14种土地利用适宜性等级，对评价区进行了土地利用区划，圈出了未来城市发展的适宜地段和高风险区，在此基础上建立了城市整体化决策支持系统。

综上所述，可以看出，国外尤其是发达国家将 GIS 应用于地质灾害研究起步较早（表1），研究程度已远远超过我们，此方面的应用也随着GIS技术的自身发展而深入（黄润秋，2001）。

2.GIS在地质矿产勘查中的应用

地质信息系统与现代地球及其相关科学日益增长的需求相适应，以处理地球上任何具有空间方位的海量信息为特征，具定量、定时、定位等优点，近10年来已在地质矿产勘查中得到广泛应用。一个区域各种地质资料（图形、图像、文字、逻辑、数值）的GIS分析实际上代表该区域现阶段较为客观的总认识。目前，野外收集资料、数据建库、GIS分析等尚存在规范化、标准化等问题，GIS本身解决诸多专业性较强地质问题的能力亦不足。但GIS的进一步发展与完善必将使地质矿产勘查进入一个数字化的新时期（周军等，2002）。

GIS因解决地质问题而产生，其雏形可以追溯到20 世纪60 年代。加拿大测量学家R.F.Tomlinson首先于1963年提出地质信息系统这一术语，建成世界上第一个GIS即加拿大GIS（CGIS）一并应用于资源管理与规划。1970～1976年间美国联邦地质调查局建成50多个信息系统并进行综合地质研究，德国在1986 年建成DASCH系统，瑞典、日本等国也陆续建有自己的GIS。GIS的发展与计算机科学的高速发展并行，主要发生在过去的20年中，而近10年来发展更快（周军等，2002）。

表1 国外GIS在地质环境与地质灾害研究中的应用

GIS因解决地质问题而产生，其雏形可以追溯到20 世纪60 年代。加拿大测量学家R.F.Tomlinson首先于1963年提出地质信息系统这一术语，建成世界上第一个GIS即加拿大GIS（CGIS）一并应用于资源管理与规划。1970～1976年间美国联邦地质调查局建成50多个信息系统并进行综合地质研究，德国在1986 年建成DASCH系统，瑞典、日本等国也陆续建有自己的GIS。GIS的发展与计算机科学的高速发展并行，主要发生在过去的20年中，而近10年来发展更快（周军等，2002）。

ArcInfo与ArcView GIS是当前最流行的两个软件包，为美国ESRI（Environmental Systems Research Institute，Inc.）的重要产品，被许多国家官方确定为国土资源、地质、环境等管理、研究的主要地质信息系统。ESRI始建于1969年，由Jack Dansermond和Laura Dangermond用自己平时积蓄的1100美元起步，经过20世纪70年代的艰苦奋斗，1981年推出新型ArcInfo，1986年微机版的PC ArcInfo投入市场，1991 年又一力作ArcView GIS问世。1981年ESRI在其Redlands总部召开首次用户会议，仅18人到场，而1998年的用户大会有来自90个国家的8000多位代表。

ESRI的发展史反映了GIS从无到有、从弱到强、迅速成长壮大的发展历程，也从一个侧面显示出GIS巨大的市场潜力和难以估量的应用价值。

据悉，1995年市场上有报价的GIS 软件已达上千种，但主要占据市场的不过10 余种。除上述提到的ArcInfo与ArcView GIS外，国外的GIS代表作还有MapInfo、ErMapper、Idrisi Endas、Erdas、Genamap、Spans、Tigris等。

GIS已在地质矿产勘查中得到广泛应用，并取得许多瞩目成果。美国、加拿大、澳大利亚早在1985～1989年就将其应用于地质矿产调查和填图。目前，澳大利亚开始利用笔记本电脑以数字形式采集野外地质数据，建立有关数据库，借助ArcInfo与ArcViewGIS编制第二代地质图件。

三、资料来源

曹修定，阮俊等.2007.GIS技术在地质灾害信息系统中的应用.中国地质灾害与防治学报，18（3）：112～115

黄润秋.2001.面向21世纪地质环境管理及地质灾害评价的信息技术.国土资源科技管理，18：30～34

姜作勤.2008.国内外区域地质调查全过程信息化的现状与特点.地质通报，27（7）：956～964

姜作勤，张明华.2001.野外地质数据采集信息化所涉及的主要技术及其进展.中国地质，28（2）：36～42

赵金平，焦述强.2004.基于GIS的地质环境评价在国外的研究现状.南通工学院学报（自然科学版），3（2）：46～50

周军，梁云.2002.地理信息系统及其在地质矿产勘查中的应用.西安工程学院学报，24（2）：47～50

GIS在地质学中的应用

石油和矿产勘查要求多种数据集进行综合分析。过去对数据存档、检索及迭加分析通常使用图件或表格数据,对比与综合要花费大量时间,遥感与GIS技术则为这些多源勘探数据综合处理提供了现代化手段。

在石油等矿产勘查时,地质学家首先要对各种地质图件、地球物理和地球化学数据、地震剖面以及遥感图像等数据进行综合分析,以便能清楚地了解各种不同数据集之间的关系。

地质数据通常也是由点、线、多边形三种形态构成的。点数据以地球化学分析数据最典型,它与某一特定的取样点有关;线数据可以是一条岩性分界线或一条断裂;多边形数据如某种岩类的出露范围。这些数据,有的采用图件形式,用颜色表示岩石类型(专题图),符号表示地球化学取样点位置,用等值线表示磁场测量值。许多地质数据还以报告、图形或实验室结果表格等形式提供。在GIS中,这些不同的数据集(如地球化学分析数据、航磁调查数据、地震数据、地质图和地形图以及遥感数据)经过数字化、编码、矢量到网格数据转换,产生连续或离散的数据集,存入建立起目标区的地质数据库,图13-1给出了地质地表数据的输入,分析和建库的过程。

在地质数据库中,地质数据按专题内容分层存贮,几何特征以图形图像表达,属性数据则记录在二维关系表中,两者为一对一或一对多的关系。于是,在这个数据模型的基础上,勘探工作区的所有地球物理、地球化学、岩石学及辐射场的数据都可以纳入数据库。一旦工作区的地质数据库被建立,地质学家便可以利用已有的专家(概念)模型来指导数据分析。例如,在石油勘探中,首先利用石油存贮条件与变量之间已知的物理、化学和地质联系来分析数据库提供的数据,对直接或间接与这些联系有关的数据进行分析、处理、生成各种派生数据。表13-1显示某工作区地质数据库中的原始数据和派生数据集。用这些数据所提供的信息来选定油气储藏有利地区。

如将重力和航磁数据叠合,有助于对基底形态的分析。又由于基底形态对沉积盖层构造发育有影响,因而据重力和航磁的一阶、二阶导数可推断出构造的总体特征。又如,基底隆起地区可能影响盖层构造特征,基底凹陷的地区沉积厚度较大,可能成为盆地的沉积中心。

图13-1 地质地表数据处理、分析及建库流程图

背斜构造是重要储油构造。是油气勘探数据库的重要内容。构造的向下延伸范围是一个最有价值的参数,目前的技术水平还难以确定。在数据库中,背斜用多边形表示,并以背斜轴为中心向下延展来定性表达背斜的地下影响范围。

断层对油气的生、储、盖都很重要。断层等密度图与线性体等密度图是用任一网格单元范围内断层/线性体出现的频数来定义的。用邻域分析法计的研究区内围绕每一象元的5×5象元阵列中断层出现的次数。结果图显示出断层/线性体密度。将断层等密度和线性体等密度图进行叠加,合成出一幅描述断裂密度的新图。对盖层断裂密度高值地区进行分析,判明它对区域油气运移和储集的具体作用。

表13-2给出某研究区域模型及其对应的权重,系统据此运行后生成一个新图像。图像的像元值等于各输入的权值求和,将它们进一步分段,便可以表达工作区中油气产出有利性的不同级别,最后圈出高概率产油区。

这种技术方法同样适用于其它矿产勘查、区域成矿预测,工程地质灾害评估与预测等。

GIS技术的引入可能极大改变地质学家的工作模式,使地学工作者面临的对多源地质数据的采集、配准、存储、分析、综合与检索工作,变得形象直观、灵活多样、快速准确,使各种地学模型的生成和发展,在技术上有了主要的支撑系统。

表13-1 原始和派生地质数据

表13-2 模型的输入与数字加权

GIS可以用于地质灾害和土地污染模拟吗？

当然可以。地理信息系统 (GIS) 是表达处理以及分析与地理分布有关的专业数据的一种技术，它提供了一种快速展示有关地理信息和分析信息的新的手段和平台。 从20世80年代以来，GIS在灾害管理中得到逐步深入的应用：从简单的灾害数据管理、多源数据集数字化输入和绘图输出，到DTM和DEM模型的建立和使用；从GIS结合灾害评价模型的扩展分析；到GIS与决策支持系统的集成；再到WebGIS。GIS的核心是空间数据管理子系统，由空间数据处理和空间数据分析构成。运用GIS所具有的数据采集和提取、转换与编辑、数据集成、数据的重构与转换、查询与检索、空间操作与分析、空间显示和成果输出及数据更新等功能，我们可以根据地质灾害评估的需要，建立以GIS技术为基础的、用于地质灾害评价的空间分析模型，评价结果可以图层的形式显示或者报表、表格形式输出，为专业部门或决策部门提供灾害管理和决策依据。武汉智博创享已经在这方面有不少案例了，可以网上查查看。

地质灾害风险评估方法

滑坡泥石流等地质灾害的不确定性决定了其评估方法采用非确定性分析方法。该类方法是基于地质灾害预测理论的广义系统科学原理，在类比法的基础上发展起来的一类研究方法。随着概率论、数理统计及信息理论、模糊数学理论用于地质灾害预测，目前已形成了多种预测模型，其预测成果可相互对比、检验，从而可使预测成果更具合理性、科学性。目前常用的非确定性分析方法主要有以下几种。

一、参数合成法

参数合成法又称专家经验指数综合评判法。它是最为简单的定量评估方法。该类模型主要是建立在专家丰富的经验基础之上的，通过专家打分法等途径获取专家经验知识，专家选择影响地质灾害的因子并编制成图。根据专家的经验，赋予每个因子一个适当的权重，最后进行加权叠加或合成，形成地质灾害危险性分区图。

它的主要优点是：①可以同时考虑大量的参数；②可以应用于任意比例尺的区域和单体斜坡稳定性评估；③大大降低了隐含规则的使用，定量化程度提高；④整个流程可以在GIS的支持下快速完成，使数据管理标准化，时间短，费用少。主要缺点有：①主观性较强，不同的调查者或专家得出的结果无法进行比较。权值的确定仍含有不同程度的主观性；②隐含的评判规则使结果分析和更新困难；③需要详细的野外调查；④应用于大区域评估时，操作复杂，模型难以推广。

二、数理多元统计模型法

该方法是通过对现有地质灾害及其类似不稳定现象与地质环境条件和作用因素之间的统计规律研究，建立相关的预测模型，从而预测区域地质灾害的危险性。该类模型方法很多，如回归分析、判别分析、聚类分析方法等。

统计分析的前提是已知学习区（训练区）的地质灾害分布情况，根据数理统计理论，建立影响参数和地质灾害发生与否的数学统计模型，在测试区得到验证后，将其应用到地质环境相同或相似的地区，预测研究区的灾害危险性分布规律。因此，统计分析方法评估的结果的可靠度直接取决于测试区原始数据的精度，模型也不能在任何地区推广使用。尽管如此，大量的研究表明，统计分析是目前最为适用的区域地质灾害危险度评估区划方法，它有严格的数理统计理论作基础，数学模型简单易懂，而且与GIS技术能够很好地结合，使庞大的数据得到合理的标准化管理、分析与储存。

多元统计分析中的主成分分析和因子分析方法在环境统计方面有不少成功的应用。将这两种方法结合起来的主成分－因子分析法可以应用于多变量的因子赋权研究（吴聿明，1991）。主成分－因子分析法的主要思想是（应农根，刘幼慈，1987）：在所研究的全部原始变量中将有关信息集中起来，通过探讨相关矩阵的内部依赖结构，将多变量综合成少数彼此互不相关的主成分，以再现原始变量之间的关系，并通过因子荷载矩阵的轴正交或斜交旋转，进一步探索产生这些相关联系的内在原因。

此方法适用于区域地质灾害空间预测研究，对一定地区土地利用、国土开发、城市规划具有宏观指导作用。

三、层次分析法

层次分析法是对一个包括多方面因子而又难以准确量化的复杂系统进行分析评估时，根据各因子之间以及它们与评估目标的相关性，理顺组合方式和层次，据此建立系统评估的结构模型和数学模型；对模型中的各种模糊性因子，根据它们的强度以及对影响对象的控制程度，确定标度指标和作用权重；将这些指标作为基本参数，代入评估模型，逐级进行定量分析并最终取得评估目标。根据地质灾害风险系统组成，大致可通过4个层次的统计分析完成评估工作：以各种要素为主体的基础层统计分析；以危险性、易损性、减灾能力为目的的过渡层分析；以期望损失为目标的准则层分析；以风险度或风险等级为最终目标的目标层分析。

四、模糊与灰色聚类方法

模糊聚类判别法模型以模糊数学理论为基础。由于地质灾害系统的复杂性，用绝对的“非此即彼”不能准确地描述地质灾害系统的客观实际，存在着“亦此亦彼”的模糊现象，不能用1或0二值逻辑来刻画，而需用区间[0,1]的多值（或连续值）逻辑来表达。而模糊数学理论正是适用于地质灾害系统的不确定性，用隶属函数来描述那些边界不清的过渡性问题及受多因素影响的复杂系统的非确定性问题。目前常用的方法有模糊综合评判法、模糊可靠度分析方法及其与层次性原理相结合而派生的模糊层次综合评判法。模糊聚类综合评估的基本步骤是：根据地质灾害风险构成，建立因素集、综合评估集和权重集，确定隶属函数，得到综合评估结果，并进行解释分析。

灰色聚类综合评估法以灰色系统理论为基础，常用于研究“小样本、贫信息不确定性”问题。在地质灾害预测中，可利用灰色关联分析，评估斜坡稳定性各影响因素的影响程度，可以克服通常数理统计方法作系统分析所导致的缺憾，对样本量和样本的规律性无特殊要求。同样可通过灰色聚类中的灰类白化权函数聚类，在考虑多种影响因素的基础上对各研究单元的危险性状态进行判定，进而完成空间预测中的危险性分区。灰色系统的以灰色模型（GM）为核心的各种预测模型还为分析地质灾害预测中的各种时序数据提供了有效途径，成为目前地质灾害实时跟踪预报的常用方法之一。灰色聚类综合评估的基本步骤是：确定聚类白化数和白化函数，标定聚类权，求聚类系数，构造类向量，求解聚类灰数。

五、信息模型评估法

该类模型的理论基础是信息论。用地质灾害发生过程中熵的减少来表征地质灾害事件产生的可能性，因素组合对某地质灾害事件的确定所带来的不肯定性程度的平均减少量等于该地质灾害系统熵值的变化。认为地质灾害的产生与预测过程中所获取的信息的数量和质量有关,是用信息量来衡量的,信息量越大，表明产生地质灾害的可能性越大。该类模型预测法同统计预测模型一样，适用于中小比例尺区域预测。

信息科学现已成为广泛使用的一门科学，但它的产生却只有短短的半个世纪历史。1948年Shannon 发表的著名论文《通信的数学理论》标志着信息科学的诞生。Shannon把信息定义为“随机事件不确定性的减少”，并把数学统计方法移植到了通信领域，提出了信息量的概念及信息熵的数学公式。信息科学研究的对象是信息，它的重要任务是研究信息的提取、信息传输、信息处理、信息存储等。由于现代自然科学发展的综合整体化趋势，各学科的相互渗透、相互联系，经过几十年的发展，使信息量和信息熵的概念已远远超出了通信领域。信息科学不仅应用于各种自然科学领域，而且已广泛应用在管理、社会等科学领域。

运用信息论方法进行地质学领域的矿床预测研究是由维索奥斯特罗斯卡娅（1968）及恰金（1969）先后提出。赵鹏大等在《矿床统计预测》一书中研究了信息量方法在区域找矿工作中的应用问题。晏同珍、殷坤龙等自1985年起，先后多次在陕南及长江三峡库区探索了信息量方法在区域性滑坡灾害空间预测分区中的应用，并与其他方法（如聚类分析、回归分析、数量化理论方法等）的研究成果进行了比较性研究。艾南山、苗天德（1987）研究了侵蚀流域地貌系统的信息熵问题，他们在斯揣勒的流域面积——高程曲线的基础上构造了侵蚀流域地貌系统的信息熵表达式，并据此作为流域稳定性的一种判定指标。Read J. 和Harr M.（1988）首次将信息熵的概念与斜坡安全系数计算的条分法结合在一起。由于地质灾害预测内容的多样性，所以决定了预测理论和方法的非单一性。晏同珍等（1989）将其概括为三类模型预测法——确定性模型预测法、统计模型预测法、信息模型预测法；前两种模型又可分别称其为“白箱”和“黑箱”模型，而信息模型则是介于两者之间。

地质灾害现象（Y）受多种因素Xi的影响，各种因素所起作用的大小、性质是不相同的。在各种不同的地质环境中，对于地质灾害而言，总会存在一种“最佳因素组合”。因此，对于区域地质灾害预测要综合研究“最佳因素组合”，而不是停留在单个因素上。信息预测的观点认为，地质灾害产生与否是与预测过程中所获取的信息的数量和质量有关，因此可用信息量来衡量：

地质灾害风险评估理论与实践

根据条件概率运算，上式可进一步写成：

地质灾害风险评估理论与实践

式中：I（y,x1x2xn）为因素组合x1x2xn对地质灾害所提供的信息量（bit）；P（y,x1x2xn）为因素x1x2xn组合条件下地质灾害发生的概率；Ix1（y,x2）为因素x1存在时，因素x2对地质灾害提供的信息量（bit）；P（y）为地质灾害发生的概率。

式（2）说明，因素组合x1x2xn对地质灾害所提供的信息量等于因素x1提供的信息量，加上因素x1确定后因素x2对地质灾害提供的信息量，直至因素x1x2xn-1确定后，xn对地质灾害提供的信息量，反映出信息的可加性特征，从而说明区域地质灾害信息预测是充分考虑因素组合的共同影响与作用。

P（y,x1x2xn）和P（y）可用统计概率来表示，各种因素组合对预测地质灾害提供的信息量可正可负，当P（y,x1x2xn）gt；P（y）时，I（y,x1x2xn）gt；0；反之I（y,x1x2xn）lt；0。大于0情况表示因素组合x1x2xn有利于所预测地质灾害的发生，相反情况则表明这些因素组合不利于地质灾害的发生。

区域地质灾害预测是在对研究区域网格单元划分的基础上进行的，根据不同地区具体的地质、地形条件，采用相应的网格形状和网格大小，进一步结合区域地质灾害分布图开展信息统计分析。假定某区域内共划分成N个单元，已经发生地质灾害的单元为N0个。具相同因素x1x2xn组合的单元共M个，而在这些单元中有地质灾害的单元数为M0个。按照统计概率代表先验概率的原理，式（1），因素x1x2xn在该地区内对地质灾害提供的信息量为：

地质灾害风险评估理论与实践

如果采用面积比来计算信息量值，则式（3）可表示成：

地质灾害风险评估理论与实践

式中：A为区域内单元总面积；A0为已经发生地质灾害的单元面积之和；S为具相同因素x1x2xn组合的单元总面积；S0为具相同因素x1x2xn组合单元中发生地质灾害的单元面积之和。

一般情况下，由于作用于地质灾害的因素很多，相应的因素组合状态也特别多，样本统计数量往往受到限制，故采用简化的单因素信息量模型的分步计算，再综合叠加分析相应的信息量模型改写为：

地质灾害风险评估理论与实践

式中：I为预测区某单元信息量预测值；Si为因素xi所占单元总面积；S0i为因素xi单元中发生地质灾害的单元面积之和。

六、实证权重法

实证权重法（Weights of evidence，）是加拿大数学地质学家Agterberg等（1989）提出的一种基于二值（存在或不存在）图像的地学统计方法，是在假设条件独立的前提下，基于贝叶斯定理（Bayesian’rule）的一种定量预测方法。Bonham-Carter等（1990）和Harris等（2001）都先后应用WOE方法来预测矿产的远景分布。通过对已知成矿情况网格单元的预测因子和响应因子之间的统计分析，计算出权重，然后对各待预测网格单元的各预测因子进行加权综合，最后，通过确定每一单元响应因子出现的概率大小便可得到不同级别的成矿远景区。

Van Westen进一步将模型应用到灾害危险性评估领域。数据驱动权重模拟方法的主要原理是利用滑坡历史分布数据，建立滑坡分布与各影响因子之间的统计关系，即根据在各影响因子不同类别中滑坡分布的统计情况来确定各影响因子对滑坡灾害的贡献率（权重）大小。这种采用数据进行权重确定的方法被称为数据驱动模型。与专家知识模型相比，权重的确定更加科学和可靠，避免了专家的主观性所带来的不确定性。最后，利用另一时期的滑坡分布历史数据对评估结果进行检验和成功率预测，调整不合理的边界，使评估结果更加具有可信度。基于统计学的Bayesian方法的数据驱动权重模型所采用的统计方法更加严谨，充分考虑了滑坡影响因素之间的关系，以及各影响因素与滑坡灾害的关系；并进行影响因素的独立性分析，找出最关键的影响因子。在此基础上计算各影响因素的权重。

七、非线性模型预测法

非线性模型预测法又称BP神经网络法，是把一组样本的输入输出问题变为一个非线性优化问题而建立的预测模型。

鉴于地质灾害系统具有复杂性特点，很难用简单的线性方程表达，因此使一批非线性预测模型迅速发展起来。如分形理论就是通过研究地质灾害系统的自相似性来对地质灾害的运动规律进行研究。易顺民应用分形理论研究了区域性滑坡灾害活动的自相似结构特征，发现在地质灾害活动的高潮期到来前有明显的降维。吴中如、黄国明等依据分形理论提出了滑坡变形失稳判据及滑坡蠕滑的相空间模型，是地质灾害时间预报的一种全新思路。自组织理论探索地质灾害复杂系统如何从无序进化到有序的自组织过程；突变理论主要从定量的角度描述非线性系统在临界失稳时的突变行为，为地质灾害时间预报提供了一种新途径；分形理论则从几何的角度探讨系统内各个层次间的自相似性，应用在地质灾害过程描述及过程预报中，化复杂为简单，化定性为定量；混沌动力学探讨非线性地质灾害系统在其演化过程中的不可逆性和演化行为对初值的敏感性。

人工神经网络（Artificial Neural Network，简称ANN）是由大量与自然神经细胞类似的人工神经元广泛互连而成的网络。网络的信息处理由神经元之间的相互作用来实现，知识与信息的存贮表现为网络元件互连间分布式的物理联系，网络的学习和识别决定于各神经元连接权系的动态演化过程。人工神经网络是一个超大规模非线性连续时间自适应信息处理系统。目前人工神经网络的应用已渗透到许多领域，为学习识别和计算提供了新的现代途径。

人工神经网络使用比较方便，它的信息处理过程同人脑一样，是一个黑箱，如图1-6所示。在实际应用中，和人们打交道的只是它表层的输入和输出，而内部信息处理过程是看不到的。对于不懂神经网络内部原理的人，也可将自己的问题交给这种网络进行解决，只要把你的例子让它学习一段时间，它就可以解决与之有关的问题。这正符合地质灾害预测理论的基本原理和思路。

图1-6 神经网络信息处理示意图

根据人工神经网络对生物神经系统的不同组织层次和抽象层次的模拟，人工神经网络可以分为多种类型。目前已有40余种人工神经网络模型。引用于地质灾害预测评估的多层前馈神经网络模型（Back Propagation,简称BP模型）是目前应用最广泛、发展最成熟的一种神经网络模型，如图1-7所示，它是按层次结构构造的，包括一个输入层、一个输出层和一个或多个隐含层。

图1-7 BP网络模型

实际上，BP模型是把一组样本的输入输出问题变为一个非线性优化问题。我们可以把这种模型看成一个从输入到输出的映射，这个映射是高度非线性的。如果输入节点数为n，输出节点数为m，则神经网络表示的是从n维欧氏空间到m维欧氏空间的映射。

在预测识别过程中，标准样本的选择是否得当，是预测是否成功的关键。一般来说，学习样本最好能涵盖预测对象的所有状态，具有广泛的代表性。在确定网络结构时，一般来讲，一个隐层的三层BP模型已可进行任意精度模拟任何连续函数。隐含层结点数目过少，不能有效地映射输入层和输出层之间的关系；过多，收敛速度过慢。因此，中间层结点数目的选取，需经过反复演算训练，才能得出较为理想的节点数。在计算过程中，为了提高效率，可以适当降低输入结点的数目，减少训练样本的维数，以增加网络的稳定性，同时还可以通过增加冲量项法或者自适应调节学习率、共轭梯度法等方法提高迭代收敛速度。

BP模型运用到地质灾害危险性区划中，可以通过样本区的标准样本的学习建立相应预测网络，从而推广到预测区进行预测。网络的输入层的变量对应于影响地质灾害产生的主要影响因素，变量可以是二态变量，也可以是具体的观测数据。当然由于各变量存在单位或数量级的差异，必须把变量数据经过正规化或标准化处理。输出层对应的是地质灾害预测等级（极高、高、中等、低、极低）的划分，或是危险程度的具体数值表达，如稳定性系数、破坏概率等，这就要求样本区的研究精度较高，指标细化程度较高。

八、地质灾害风险分析与GIS技术

地理信息系统（GIS）是集计算机科学、信息科学、现代地理学、遥感测绘学、环境科学、城市科学、空间科学、管理科学和现代通讯技术于一体的一门新兴学科。具体而言，GIS是指对各种地理信息及其载体（文字、数据、图表、专题图等）进行输入、存储、检索、修改、量测、运算、分析、输出等的技术系统。GIS的主要功能有采集、存储、管理、分析、输出各种数据、数据维护和更新、区域空间分析以及多因素综合分析和动态监测等。GIS不仅可以像传统的数据库管理系统（DBMS）那样管理数字和文字（属性）信息，而且还可以管理空间（图形）信息；它可以使用各种空间分析的方法，对多种不同的信息进行综合分析，寻找空间实体间的相互关系，分析和处理一定区域内分布的现象和过程。当代地理信息系统正向能够提供丰富、全面的空间分析功能的智能化GIS的方向发展。智能化的GIS具有强大的空间建模功能，能够构建各种具有专业性、综合性、集成性的地学分析模型来完成具体的实际工作，解决以前只有靠地学专家才能解决的问题。

GIS把各种与空间信息相关的技术与学科有机地融合在一起，并与不同数据源的空间与非空间数据相结合，通过空间操作与模型分析，提供对规划、管理、决策有用的信息产品。GIS为我们提供了一种认识和理解地学信息的新方式，GIS强大的空间分析功能和空间数据库管理能力为我们研究区域地质灾害提供了一个科学、便捷的崭新途径。

作为数字地球的核心技术之一，GIS经过将近40年的发展，已经成为一种日益成熟的空间数据处理技术和方法。它提供了一种认识和理解地学信息的新方式，已广泛应用于国土资源调查、环境质量评估、区域规划设计、公共设施管理等方面。在地质灾害研究领域，GIS技术的应用已从最初的数据管理、多源数据采集数字化输入和绘图输出，到数字高程模型、数字地面模型的使用、GIS 结合灾害评估模型的扩展分析、GIS与决策支持系统的集成、GIS虚拟现实技术的应用等，并逐步发展与深入应用。

各种地质灾害都是在地球表层一定空间范围和一定时间限度内发生的，尽管不同种类的地质灾害之间、同一种类的地质灾害的不同个体之间大都形态各异，形成机理也是千差万别，但它们都是灾害孕育环境与触发因子共同作用的结果，而这些都与空间信息密切相关，利用GIS技术不仅可以对各种地质灾害及其相关信息进行管理，而且可以从不同空间和时间的尺度上分析地质灾害的发生与环境因素之间的统计关系，评估各种地质灾害的发生概率和可能的灾害后果。地质灾害危险性区划图属于一种综合图件，而且具有一定时段内的静态特点，因此需要不断更新；尤其是有新的地质灾害发生的时候，更应及时修订。由于GIS技术的空间分析、制图功能和可视化的特点，所以GIS技术在地质灾害区划研究方面正得到快速发展，以GIS软件为技术平台的地质灾害的危险性、易损性和风险评估的系统研究逐步成为本领域研究的发展方向，并有可能在不远的未来与网络技术相结合。

国外尤其是发达国家，对GIS技术应用于地质灾害领域的研究已做了很多工作。从20世纪80年代至今，GIS技术的应用已从数据管理、多源数据采集、数据化输入和绘图输出，到数字高程模型、数字地面模型的使用、GIS结合灾害评估模型的扩展分析、GIS与决策支持系统（DSS）的集成、GIS虚拟现实技术的使用，都得到不断的发展和广泛的应用。在滑坡灾害研究领域，GIS技术的应用已经比较成熟，主要体现在以下几个方面：

（1）建立基于GIS的滑坡灾害信息管理系统。如Keane James M.（1992）, BaharIrwan（1998）, Bliss Norman B.（1998）等将GIS运用到滑坡灾害历史数据的管理及预测成果成图表征中。

（2）GIS技术与各种评估模型结合运用到滑坡危险性预测中。如Matula（1987），Lekkas E.（1995）, Randall（1998）, Dhakal Amod Sagar（1999）等利用GIS的空间分析功能与预测模型的结合，完成了滑坡预测因素的空间叠加，进行滑坡危险性预测，得出相应的预测分区图和滑坡敏感性图。

（3）进行基于GIS的滑坡灾害风险分析预测与管理。如 Ellene（1994）,Leroi（1996）,Bunza（1996）, Castaneda Oscar E.（1998）, Atkinson（1998）, Michael（2000）, Aleotti（2000）等从影响滑坡灾害风险的因素出发，利用GIS的空间分析功能进行因素叠加，实现风险评估并结合GIS的信息管理功能，对灾害信息进行管理，最终进行管理决策，大到防灾减灾的目的。目前，国外在滑坡灾害预测领域已基本实现了RS与GIS的紧密结合，个别项目已达到了3S技术的结合。

国内基于GIS技术开展地质灾害评估工作起步较晚，目前还没有成熟实用的地质灾害预测评估的GIS系统。姜云、王兰生（1994）在山区城市地面岩体稳定性管理与控制中应用了GIS技术，以重庆市为典型研究对象，对地面岩体变形破坏进行了时空预测预报；同时，通过分析城市地质环境对土地工程利用的制约关系，应用GIS的信息存储、查询、空间叠加运算及DEM模型等功能，做出地力等级划分，并编制了斜坡稳定性综合评估分区图。雷明堂、蒋小珍等（1994）将GIS技术运用在岩溶塌陷评估中，完成了研究区岩溶塌陷危险度评估及分区。成都理工学院（1998）和中国地质环境监测院及国土资源部长江三峡地质灾害防治指挥部合作进行了“地质灾害信息系统及防治决策支持系统”开发试验工作，初步建立了一个全国地质灾害调查与综合评估系统。中国国土资源经济研究院、中国地质大学、中国地质科学院岩溶地质研究所、国土资源部实物地质资料中心（2002）联合开展了“全国地质灾害风险区划”项目攻关，利用国产软件MAPGIS，对全国小比例尺滑坡、泥石流、岩溶塌陷地质灾害进行了基于GIS的风险评估（包括地质灾害危险性评估、易损性评估和风险性区划）。朱良峰等在国产版权的MAPGIS软件平台上，开发了一套地质灾害风险评估系统RISKANLY。这套基于GIS技术的地质灾害风险分析不仅方法上可行，而且技术上先进，代表着地质灾害风险分析的发展方向。当然，无论是地质灾害的危险性分析模型，还是区域社会经济易损性分析模型，都有待于实践中的进一步研究与发展，这显然是应该随着人类对地质灾害本质属性认识的逐渐深化而不断发展的。

随着我国社会经济的迅速发展和城市化进程的加快，崩塌、滑坡、泥石流、地面塌陷等地质灾害破坏的广度与深度也在迅速增大，需要更加关注地质灾害的区域时空预测研究。与地质灾害有关的相关因素很多且成因复杂，都与空间信息密切相关，因此，利用GIS技术不仅可以对地质灾害相关的各种空间信息进行管理，而且可以从不同的空间和时间尺度上分析地质灾害的发生与环境因素之间的统计关系，评估地质灾害的发生风险和可能的灾害范围。因此，基于GIS的地质灾害风险评估与区划将会在未来我国的社会经济发展中起着重要的作用。

九、小结

地质灾害风险评估涉及两个重要的方面：一是地质灾害发生的可能性问题，二是人类自身、社会及环境等对象对地质灾害的抵御能力问题。因此，地质灾害的定义采用国际上的geological hazard一词。本书遵循科学性、通用性的原则，结合国内近年来在地质灾害风险评估领域已初步形成的有代表性的术语表达方式，在联合国教科文组织提出的统一定义的基础上，对地质灾害风险评估所涉及的基本术语定义如下：

（1）危险度H（Hazard）。特定地区范围内某种潜在的地质灾害现象在一定时期内发生的概率。

（2）易损性V（Vulnerability）。某种地质灾害现象以一定的强度发生而对承灾体可能造成的损失程度，易损性可以用0-1来表示，0表示无损失，1表示完全损失。

（3）承灾体E（Element at risk）。特定区域内受地质灾害威胁的各种对象，包括人口、财产、经济活动、公共设施、土地、资源、环境等。

（4）风险度R（Risk）。承灾体可能受到各种地质灾害现象袭击而造成的直接和间接经济损失、人员伤亡、环境破坏等。风险等于危险性、易损性、承灾体价值三者的乘积。

风险度（R）=危险度（H）×易损度（V）×承灾体价值（E）