十年网站开发经验 + 多家企业客户 + 靠谱的建站团队
量身定制 + 运营维护+专业推广+无忧售后,网站问题一站解决
RS+GIS即遥感结合地理信息技术求取滑坡规模。在获得滑坡后的OrbView立体像对数据后,制作了能较精确反映滑坡后地形的5m栅格的DEM,以下简称为滑后DEM,然后由以下步骤求取滑坡规模。
创新互联自2013年起,先为西充等服务建站,西充等地企业,进行企业商务咨询服务。为西充企业网站制作PC+手机+微官网三网同步一站式服务解决您的所有建站问题。
(一)滑后图像配准
以离滑坡时间最近的2004-09-11QUICKBIRD图像作为滑坡后形态的信息源,使其与滑后DEM精确配准,该DEM反映滑坡后地面各点的高程,该项工作在MAPGIS的图像处理镶嵌配准功能中完成。如图3⁃9。
图3⁃9滑坡与DEM精确配准
(二)求滑坡前DEM
如有同样精度的滑坡前后DEM,则求其高程差后即可方便地再求得滑坡体积。但本项目实际情况是:在滑坡范围只有1:5万数字地形制作的20m栅格DEM反映滑前地形,由于在高山峡谷地区,且地形图制作年代已较长,在滑坡体上只有很少的点并是经圆滑后的数据,不能准确反映滑前地形。本研究提出用滑后DEM制作滑前DEM,其步骤如下:①将滑前图像与滑后DEM配准,主要是边界部分配准,如图3⁃9,形成高精度正射影像;②将滑后DEM沿垂直滑动方向分为若干剖面,本次划分25条剖面,剖面上每一点的DEM值为各剖面上该点滑后的高程值,其中红点为该剖面的高程变化最大点,其与滑坡边界套合较好;紫色点为剖面上的最低点,与滑动面套合较好,如图3⁃10;③求各剖面的滑前地形线,将滑后DEM套合在滑前正射图像上,参考原滑前20m栅格DEM,根据图像上的滑前地形,以一定间隔逐点修正各剖面滑前DEM,获得各剖面的滑前DEM,并与滑后DEM套合。这样,各剖面的上界线(图3⁃10中各条剖面的虚线)为滑前DEM值,下底线(各剖面线的实线)为滑后DEM值;④在MAPGIS平台求得各剖面的面积ΔS,如表3⁃4;⑤以公式
ΔSdy计算滑坡体积,1和n为剖面编号,S为每一块剖面的面积,dy为每一个剖面的间隔。间隔越小,求取的滑坡体积越准确。根据以上步骤求得:2004-06-22滑坡的体积为63.7×104m3,滑体平均宽(沿河)150m,平均长(顺坡)180m,平均厚24m。
图3⁃10 求取滑坡前后地形示意
表3-4 2004-06-22滑坡各剖面面积
本计算值与上述滑坡规模估算值比较减小了9.3×104m3。平均厚度值减小了6m。RS+GIS方法是将滑坡分为25块剖面来求规模的,以各剖面上两侧高程变化最大点作为滑坡边界,以剖面上最低位置作为滑面位置,滑前地形也是逐条剖面逐点修正的,故有较可靠的依据,其结果也更准确。但前一种估算方法较简单,两者相比,精度可达86%左右,需要时还是可以采用的。
B.1 填图
在地图上表示滑坡灾害是一种既实用又方便的手段,它可以表示关于滑坡的多种信息,并可以以不同的详细程度进行组合。当灾害图件与土地利用图件并行使用时,它们又成为很有价值的规划工具。通常,滑坡灾害填图可以分成三个阶段来实现。第一阶段,进行区域性填图,或初勘填图,主要是统合各种已有的数据,并确定可能构成问题的地区。这种区域性(有时称为“小比例尺”)填图常由省、州或联邦政府地质调查局来完成。第二阶段是社区范围填图(有时称为“中比例尺”),主要针对有问题的复杂地区进行较详细的地表和地下填图。第三阶段,进行具体场地的大比例尺填图。如果资源和经费有限,可以越过区域性填图阶段,直接把注意力放在那些已知的、有问题的地区。
◆三种类型的综合性填图
这里介绍三种类型的综合性填图:①区域性填图;②社区范围填图;③具体场地的详细填图。
区域性填图
区域性填图或称初勘填图,是为了给区域规划提供基础数据,主要提供今后的两个阶段的填图所需要的基本信息,设定将来填图工作的优先顺序。
这种类型的图件主要是简单的历史滑坡发生记录,或滑坡易发生地可能性分布图。通过这类图件,可以在区域上确定构成滑坡灾害问题的区域,了解滑坡形成的条件。该类图件主要关注那些可能发生滑坡运动的地质单元或环境。这类填图主要依赖于图像解析(对航空照片的地质解译),初勘阶段的野外填图,以及对所有相关地质资料的收集和整理。填图比例尺为1:1万~1:400万,甚至更小。
社区范围填图
这种类型的填图既要在三维空间上确定滑坡发生的可能性,又要考虑滑坡发生的条件和原因。在此阶段,必须制定有关土地利用、分区及建筑的指南,对下一阶段的详细填图也必须提出详细的建议。地下部分的勘测也必须在此阶段进行,从而提供附有剖面图的地图。填图比例尺为1:1000~1:1万。
具体场地的详细填图具体场地的详细填图涉及该场地具体问题的确定、分析和解决方案,常常要求精确到针对某个居民点。这一阶段的填图往往由计划进行场地开发的土地所有者,委托给私人咨询公司实施,内容一般包括进行详细钻孔记录的钻探、取样、实验室实验分析,以获取设计和建设时必须的资料。这个阶段的填图比例尺可以有变化,但多为1:500或1:600。
◆滑坡分布图的三个重要准则
对于设计师和一般公众来说,滑坡分布图有三种类型:①滑坡编录图(landslide inventories);②滑坡易发性评价图(landslide susceptibility maps);③滑坡灾害危险性分区图(landslide hazard maps)。
滑坡编录图
滑坡编录图标明那些曾经毁于滑坡过程的地区(图B1)。这类图件的详细程度可能有所不同。简单图件中可能只显示滑坡发生过的较广泛的地区;而详细图件中可能对每一处滑坡都进行详细描述并进行分类:显示滑落崖的位置,区分沉陷区和堆积区,区分正在活动的滑坡和已经停止活动的滑坡,滑坡开始活动的地质年代,滑坡的运动速度,以及其他有关滑坡体厚度和物质成分的资料及数据。
简单的滑坡编录图给出整个区域的滑坡分布概观,并指出在哪些地区需要进行更详细的调查研究。详细的滑坡编录图能够提供某地区不同滑坡过程的解释,能够用来对滑坡易发生地区的开发行为进行限制或制止,并且为滑坡防治工程措施的设计提供建议。它们还能够为其他派生图件的准备提供良好的基础,这些派生图件包括以确定土地利用为目的的边坡稳定性图和滑坡灾害危险性分级图等。一种途径是使用航空照片,对其中一些被选择的地区(滑坡发生可能性大的地区)进行现场调查和确认,然后以编号的形式在地图上标识所有滑坡信息。在地图上显示的滑坡信息,必须包括以下的部分或全部内容:活动状态,滑坡确定的准确程度,边坡运动的主要类型,滑坡体的估计厚度,滑坡体物质的类型,以及活动的
日期或时期。
在美国,区域性地图最常用的比例尺是1:2.4万(在加拿大为1:5万),这是因为美国地质调查局广泛提供的高质量地形图采用了这一比例尺,并且与之相匹配的航空照片也很容易找到。其他的在美国常用的地图比例尺为1:5万、1:10万和1:25万。
滑坡易发性评价图
滑坡易发性评价图是基于历史滑坡目录图的滑坡图件,它描述那些具有滑坡潜在危险性的地区(图B2)。通过对那些可能引起滑坡的主要因素(如陡峭的边坡,在被饱和或受扰动后将大幅度失去强度的软弱的地质单元,排水不良的岩体或土体)与过去的滑坡分布进行关联分析,而得到将来的滑坡易发性评价结果。这些图件仅表示边坡的相对稳定性,它们并不作出绝对的滑坡发生危险性预测。
滑坡易发性评价图可以被认为是从历史滑坡编录图中派生出来的,因为历史滑坡编录图是滑坡易发性评价图的最基本的图件。例如,如果将一张地质图叠加在一张显示有既存滑坡分布的历史滑坡编录图上,就可以确定特定的易于产生滑坡的地质单元。通过外推这种信息,可以预测具有潜在滑坡可能性的其他地区。一些更复杂的图件可以包括更多的信息,例如边坡坡角和排水条件等。
滑坡灾害危险性分区图
滑坡灾害危险性分区图显示灾害过程所造成的危险性的空间分布(图B3):在哪些地区曾经发生过滑坡,现在正在活动的滑坡,而最重要的是将来在哪些地方可能发生滑坡。对于某个特定的地区,灾害危险性分区图包括以下详细信息:滑坡的类型,发生破坏的边坡的范围和滑坡运动将影响的最大范围。这些图件将能用来预测某地区的滑坡灾害的相对程度。其结果,一个地区可能被分成不同的灾害层次,如低、中、高危险性地段。
B.2 遥感技术和其他监测滑坡活动特征的手段
地图和其他形式的信息有时可以通过GIS(地理信息系统)进行叠加,从而使得人们在一张图上可以看到不同类型的信息。如果没有计算机化的GIS系统,通过使用透明图进行迭加的方法也可以达到这一目的。重要的是,在叠加的时候,各种图件和数据必须使用同样的比例尺。下面列举并描述了多种可以构筑GIS数据层以进行滑坡潜在危险性分析的信息类型。
Topographic Map 地形图标明边坡坡度,地形构成,排水型式。
Terrain Map 地形地体图标明物质,深度,地质过程,地形构成,地表和地下排水,边坡坡度(也被称为表面地质或第四纪地质图)。
Bedrock Map 基岩图确定基岩类型,地表和地下结构,地表覆盖层(上覆层)和岩石的年龄。
Engineering Soil Map 工程土性图确定地表物质类型,排水条件,被限定的工程特征,土质特征和植被条件。
Forest Cover Map 森林覆盖图确定地表植被,地形特征,地表排水型式。在某些情况下,还需确定土体排水特征。
Research Studies 综合研究可以提供以上内容的有关信息,以及控制性因素的定量数据,在可能的情况下,还可以提供相关地区的局部稳定性的风险评价。
Aerial Photography Remote Sensing 航空遥感照片解译通过航空遥感照片的解译,可以确定植被、地形、排水型式、土体的排水特征、基岩地质、滑坡类型以及与其他影响因素之间的关系(实例显示于图B4~图B7中)。通过使用倾斜航空照片(Oblique aerial photographs)和垂直立体照片对(Vertical stereo photograph pairs)对某特定地区地形特征进行细致解译,可以得到有关滑坡类型、发生频度以及管理活动的重要信息。如果可能的话,可以对比同一地区过去的和现在的航空照片。因为在现在的照片上,有些滑坡可能变得不是那么太明显。在航空照片上辨识的特征可以帮助确定滑坡类型,并形成对上覆层的合理评价。因此,这些结果也为场地的灾害危险性评价提供了有用的手段。
InSAR lmaging 远红外干涉图像法
InSAR是Interferometric Synthetic Aperture Radar(远红外干涉开口雷达)的缩写。无论是远红外干涉法InSAR,还是激光探测测距法LIDAR(将在下面详细叙述),都是通过高灵敏度传感器从卫星上发射能量脉冲,然后接收从地面上反射回来的脉冲。大多数远红外干涉设备发射的能量脉冲可以穿透云雾和雨滴,可以用于那些难于徒步进入的地区。通过使用从雷达卫星上发射,并经过地面反射回来的信号,可以建立能够显示地面地形的数字高程模型(DEM)地图。在两个不同的时期的作业,可以得到两个影像。将这两个影像进行合并处理,可以得到称为干涉图样(Interferogram)的图件。如果在取得两张影像的时间间隔内,该地区发生过任何位移的话,则位移发生的位置会在将这两张影像图进行合并处理时显示出来。通过这样的方式,可以确定某个地区是否发生过变形。
由于常规的卫星天线的尺寸多受到限制,典型的环绕地球的卫星雷达仅仅只具有3到4英里(4.8~6.4千米)的地面解像度。合成开口雷达(Synthetic Aperture Radar,SAR)利用航天飞船沿着其轨道轨迹运动的优势,在数学上合成可操作的更大的天线,从而得到高空间分辨率的影像(数百英尺,约为几十米)。
LiDAR lmaging 激光探测测距图像法
LiDAR是Light Detection and Ranging(激光探测测距)的缩写,有时也称为机载激光条带状填图(ALSM,Airborne Laser Swath Mapping)。使用狭窄的激光束穿透地表浓密的覆盖,如森林中的树丛,LiDAR可以产生高精度的地形图,即使在普通照片中根本看不到地面也会裸露出来。这一技术产生了非常高精度的数字标高模型(DEM)图(图B7)。最精确的裸地DEM图可以通过在树叶凋零的季节得到的LiDAR影像中解译后获得。
Li DAR填图系统的基本要素为:一个安装在飞行器上的扫描激光测距器(Rangefinder),一个为飞行器定位的差分全球定位系统(Global Positioning system,GPS),和一个测试飞行器旋转的内部测试单元(internal measurement unit,IMU)。由于以下三个原因,Li DAR成为了有用的地形填图工具。第一是精度。第二是生产效率,其测试是以每秒10,000到80,000次激光脉冲的速度来进行的。第三是因为Li DAR采用单光作业,并提供其自身的照明。这些特征克服了森林地形中照相测量中的主要缺点。Li DAR填图系统提供的图件非常清晰,由于消除了树木的影响,这些图件提供了历史滑坡的详细细节,而这些细节是其他的一般手段根本无法辨别的。Li DAR填图系统价格昂贵、技术复杂,主要被政府机关、大学和一些私人企业使用。一个主要缺点是许多Li DAR填图系统使用近-远红外激光,无法穿透雨和雾。
B.3 滑坡实时观测及滑坡监测仪器安装
能够实时提供滑坡活动观测结果的滑坡实时监测系统对于保护生命财产的安全起着至关重要的作用。传统的现场观测方法,即使采用定期观测的方式,也无法在运动发生的瞬间马上得到变形的结果。而且在正在活动的滑坡上进行调查工作往往是很危险的。再者,滑坡的大规模运动多发生在暴雨等能见度很差的条件下。由设置在远距离处的实时观测系统提供的连续数据可以比较全面地反映滑坡的动力学行为(变形发展阶段),这使得工程技术人员能够在对现场情况进行判断之后,制定出有效的防治滑坡的设计方案。滑坡观测可能会比较昂贵,而且大多数的观测系统需要由滑坡专家来进行安装。使用滑坡实时监测系统的一大优势是可以将滑坡运动的监测系统与预警系统协调使用。滑坡观测如图B8~图B10所示。
滑坡泥石流等地质灾害的不确定性决定了其评估方法采用非确定性分析方法。该类方法是基于地质灾害预测理论的广义系统科学原理,在类比法的基础上发展起来的一类研究方法。随着概率论、数理统计及信息理论、模糊数学理论用于地质灾害预测,目前已形成了多种预测模型,其预测成果可相互对比、检验,从而可使预测成果更具合理性、科学性。目前常用的非确定性分析方法主要有以下几种。
一、参数合成法
参数合成法又称专家经验指数综合评判法。它是最为简单的定量评估方法。该类模型主要是建立在专家丰富的经验基础之上的,通过专家打分法等途径获取专家经验知识,专家选择影响地质灾害的因子并编制成图。根据专家的经验,赋予每个因子一个适当的权重,最后进行加权叠加或合成,形成地质灾害危险性分区图。
它的主要优点是:①可以同时考虑大量的参数;②可以应用于任意比例尺的区域和单体斜坡稳定性评估;③大大降低了隐含规则的使用,定量化程度提高;④整个流程可以在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)
“数字滑坡”技术,就是实现“信息化滑坡”的技术。具体方法为:在滑坡地学理论指导下,以遥感和空间定位(GCPs或GPS)方法为主,结合其他调查手段识别滑坡,获取数字形式的与地理坐标配准的滑坡基本信息(滑坡各要素):位置、形态、土地覆盖、变形和位移、地质构成等;利用GIS技术存贮和管理这些数字信息;在此基础上,对滑坡及其发育环境信息进行空间分析,服务于滑坡调查、监测、研究、灾害评价、危险预测、灾情评估、减灾和防治措施等。
“数字滑坡”技术使滑坡遥感前期的获取、处理、存贮和显示滑坡信息的方式发生根本变化,使我们能更准确地定性、定位、定量地认识滑坡,科学有效地存贮和管理遥感调查结果,方便、快捷地传输及交流滑坡信息,从而为改善滑坡灾害调查、监测、防治、预测及滑坡理论研究工作提供更加准确、丰富的滑坡信息。
“数字滑坡”技术的实现主要依赖于四类信息科学技术的支持:①遥感技术;②数字摄影测量及图像处理技术;③GIS技术;④计算机技术。其他如灵镜(虚拟现实)技术、网络技术等,也在滑坡过程研究、滑坡遥感监测、滑坡信息存贮和传输等过程中使用。
“数字滑坡”技术大致可分为七个部分:信息源、建立解译基础、滑坡遥感解译、时空分析、现场验证、制图和滑坡数据库。