(School of Geoscience and Info-Physics,Central South University,Changsha 410083,China)
geological disaster; analytic hierarchy process; GIS; risk assessment
备注
文中选取灾害点密度、坡度分布、地貌类型及高程、岩土体类型、水系分布、断裂发育程度、降雨量分布、人类工程活动强度8个评价指标因子参与汉阴县地质灾害危险性评价; 根据各评价指标因子对地质灾害发育的影响程度,利用层次分析法计算各评价指标的权重,并运用MAPGIS为每个评价指标因子图层赋属性。基于MAPGIS空间分析功能完成了汉阴县地质灾害危险性区划工作; 区划结果分为高、中、低3种危险性区,结合汉阴县地质灾害危险性区划成果和发展规划,为汉阴县的地质灾害防治及群测群防提供一些理论性建议。
This paper chooses density,gradient distribution,geomorphic type and altitude,rock soil mass type,river distribution,fault developmental degree,precipitation and human engineering activity intensity of geological disasters as assessment indexes for risk assessment; According to the degree of influence of each evaluation factor for the development of geological disasters,completing the geological disaster risk zoning of Hanyin county by MAPGIS using the analytic hierarchy process model; the zoning results are classified as high,medium and low risk areas.Combining the results of geological disaster risk zoning and the development planning of Hanyin county,this paper provides geological disaster prevention and control,mass prediction and disaster preventionof Hanyin county with some theoretical suggestions.
引言
20世纪90年代以来,随着数理知识在各学科中的应用,相互穿插渗透,特别是GIS与数学模型结合的评价方法,使得地质灾害危险性评价有了得心应手的方法,当影响因素多、分布复杂时,层次分析理论使地质灾害评价更加清晰化、层次化、定量化,在大区域以及其他领域中的应用也尤为方便、简介。经过30多年的发展,GIS技术和数学模型的结合应用,在我国地质灾害评价及预测等方面都取得了丰硕成果,但统一的评价理论体系尚未形成,理论与实践的结合有待完善,结合遥感(RS)、全球定位系统(GPS)在地质灾害评价中的应用仍处于探讨和总结过程中。
文中结合现有地质灾害信息,拟在查明汉阴县地质灾害发育概况及其类型的基础上,阐明汉阴县地质灾害的发育特征、分布规律,分析汉阴县地质灾害的形成条件,运用GIS、层次分析法实现汉阴县地质灾害危险性评价。
1 汉阴县地质灾害发育概况
汉阴县发育的地质灾害主要有滑坡、崩塌以及泥石流3类(图1)。截止2012年底,调查的地质灾害点共计267处,其中滑坡244处,崩塌12处,泥石流11处,分别占地质灾害点总数的91.4%,4.5%,4.1%.
2 汉阴县地质灾害危险性评价模型
3 汉阴县地质灾害危险性评价成果
根据汉阴县地质灾害危险性区划结果,统计出汉阴县高危险区面积为206.87 km2,中危险区为456.2 km2,低危险区为702.81 km2,分别占全县总面积的15.15%,33.40%,51.45%.
3.1 地质灾害危险性大区这一类区主要包括月河盆地的双乳镇、城关镇、平梁镇和汉江两岸的低山区。区内地质环境脆弱,地质灾害密集发育,其地质灾害的发育与松散层的分布、厚度及人类工程活动密切相关。根据其地质环境条件的不同,将其划分为2个亚区
1)月河盆地亚区:该区位于汉阴县中部的双乳镇-城关镇-平梁镇及蒲溪镇北部。面积为164.42 km2,占全县总面积的12.40%,地貌以盆地、丘陵为主,地势相对平缓,但该区人口密集,修路、建房开挖坡脚多形成高陡临空面,在强降雨下会诱发堆积层滑坡等地质灾害。
2)漩涡-汉阳亚区:该区位于汉江两岸的低山区,面积为741.49 km2,占全县总面积的3.13%,区内地形多为缓坡洼地,有利于地表水的汇集,坡面垦耕程度较高,植被较差,人口较为稠密,在久雨或暴雨季节,地表水汇集使土体浸水失重,诱发滑坡等地质灾害。
3.2 地质灾害危险性中等区这一类区主要包括月河及北部的中河、青泥河、观音河、东沙河沿岸,漩涡镇和上七镇的低山区。根据地质环境条件及人类工程活动程度差异,将其划分为3个亚区。
1)中河-青泥河-观音河-东沙河亚区:该区为中低山河谷区,面积为242.85 km2,占全县总面积的18.31%,该区的河谷多呈不太对称的“V”字形,一侧有盘山路,该区人口较稠密,工程削坡较为严重,多形成高陡临空面,挡护措施较少,在降雨等外力的作用下,会使边坡失稳。
2)月河盆地亚区:该区位于双乳镇北部和凤凰山北麓的部分地区,面积为83.18 km2,占全县总面积的6.27%,主要为低山沟谷区,两岸边坡较为稳定。但在一些低洼的山坡凹槽地带,堆积着较厚的松散层,在久雨或暴雨的诱发下易产生滑坡等地质灾害。
3)漩涡-上七亚区:该区位于汉江北岸的黄龙洞河、冷水河和南岸的上七镇,面积为111.56 km2,占全县总面积的8.41%,属中低山河谷区,河谷两侧坡面坡度多在20°~35°,区内人口较为稠密,人类工程活动频繁,易诱发滑坡、泥石流等地质灾害。
3.3 地质灾害危险性小区这一类区主要位于汉阴县北部的秦岭中、低山区,中部的月河区和凤凰山区、南部的巴山区,按其所处的地理位置及地质环境条件的不同,将其划分为4个亚区
1)大寨子-红岩沟亚区:该区主要位于县北部的中、低山区,面积为272.95 km2,占全县总面积的20.58%,区内人烟稀少,植被发育,松散层分布范围有限,发生地质灾害的频率较低。
2)月河盆地亚区:该区位于县中部的月河一级阶地,面积为70.90 km2,占全县总面积的5.35%,虽然区内人类工程活动强烈,但该区地势平坦,发生地质灾害的可能性较小。
3)凤凰山亚区:该区位于凤凰山—梨花山一带,属中山区,面积为204.07 km2,占全县总面积的15.39%,区内人烟稀少,植被发育,人类工程活动轻微,一般不易形成地质灾害,即使形成地质灾害,也不会造成大范围的危害。
4)巴山亚区:该区位于汉江以南的双坪镇-汉阳镇,属中低山区,面积为134.75 km2,占全县总面积的10.16%,区内人口密度小,植被发育,发生地质灾害的可能性较小。
4 结 论
文中以汉阴县为研究对象,以汉阴县地质灾害调查与区划报告和地质环境为基础资料,分析总结汉阴县地质灾害的发育特征及分布规律,运用GIS,层次分析法实现汉阴县地质灾害危险性评价,并取得以下结论
1)野外实地调查表明,汉阴县发育的地质灾害有滑坡、崩塌、泥石流3类,截止2012年底,地质灾害点共计267处,其中滑坡244处,崩塌12处,泥石流11处,分别占地质灾害总数的91.4%,4.5%,4.1%.
2)在构建地质灾害危险性评价体系和选取评价指标时,选取地质灾害现状、地质环境条件和诱发因素作为评价体系的一级评价指标,选取地质灾害点密度、地形地貌、坡度分布、水系分布、岩土体类型、地质构造、降雨量和人类工程活动作为二级评价指标,根据各评价指标对地质灾害的影响程度,利用层次分析法计算各评价指标的权重,并运用MAPGIS为每个评价指标图层赋属性。
3)运用MAPGIS的空间分析功能,结合各个评价因子的属性和权重得到汉阴县地质灾害危险性区划图,其中高危险性区为206.87 km2,中危险性区为456.23 km2,低危险性区为702.78 km2,分别占研究区面积的15.15%,33.40%,51.45%.
2.1 评价模型本研究以汉阴县为例,考虑到研究区内地质灾害发育特点和现有地质灾害信息资料情况,决定采用层次分析法对汉阴县地质灾害进行危险性评价。根据属性的不同,将汉阴县地质灾害危险性评价研究分成3个层次,将最上层为目标层,最下层为方案层,中间层即为准则(约束)层,可以有一个或几个层次,同一层的各因素受上一层各因素的影响,同时又影响(支配)下一层的各因素。根据汉阴的情况,建立危险性评价模型,如图2所示。
2.2 危险性评价步骤及评价因子计算2.2.2 汉阴县地质灾害危险性评价指标的选取影响汉阴县地质灾害发育的因素有多种,主要有以下3大类:灾害历史因素、主控因素和诱发因素。灾害历史因素包括灾害点影响范围、灾害点的密度、灾害规模大小的分布,鉴于目前所掌握的地质灾害历史资料不够详细,文中仅选取灾害点密度为灾害历史因素评价因子; 主控因素即灾害的地质环境背景,均可从地质资料中获取,文中主控因素选取的评价因子有:地形地貌因子、岩土体类型因子、水系分布因子、坡度分布因子、断裂构造因子; 诱发因素选取的评价因子为:降雨量因子、人类工程活动强度。
2.2.3 评价指标的量化及权重的计算根据各评价因子对区域地质灾害危险性影响程度的大小,对评价因子量化并计算权重,然后编辑MAPGIS能够识别的图形文件,并对各评价因子添加特征属性。评价因子的量化带有人为主观性,但对评价结果影响较小。表1为本研究的评价因子量化分级表。
地质灾害危险性一般划分为3个级别,本研究的评价因子也量化为高、中、低3个级别,其量化值分别为3,2,1.
岩土体类型划分依据是建造类型和力学强度; 水系将汉江、月河划分为高危险性区,东沙河、观音河、中河、青泥河、冷水河划分为中危险性区,其他河流及其支流划分为低危险性区。
表2 地质灾害危险性评价指标重要度比较矩阵表
Tab.2 Comparative matrix of the importance of geologic hazards risk evaluation indexes1)计算矩阵A,A1,A2,A3的特征向量W及最大特征值λmax:
文中利用Excel计算判断矩阵,使得判断矩阵的构造、层次单排序和层次总排序的计算以及一致性检验和检验之后对判断矩阵的调整变得十分简单[4-6]。
对矩阵A,A1,A2,A3进行一致性检验
A,A2一致性比率分别为CR(A),CR(A2)则有
式中 A1为一阶矩阵; A3为二阶矩阵,总是完全一致的,矩阵A,A1,A2,A3均通过一致性检验,满足一致性要求。
3)计算组合权重。将特征向量W(A),W(A1),W(A2),W(A3)见表4,计算各元素组合权重值。
4)对研究区划分评价单元网格,然后将各因子区特征属性添加到单元网格中,使区域内所有网格带有所有因子的特征属性。网格越小,数量越多,评价精度越高,但计算工作量增大。
5)利用MAPGIS对各评价因子进行属性叠加,并计算各评价单元的危险性指标值,然后根据计算完的危险性指标值,生成地质灾害危险性区划图。
- [1] Forster,Alan.The assessment of landslide hazard potential as a guide to land use Planning in the South Wales Coalfield[J].National Museum of Wales,2005,24:81-87.
- [2] Wooten,Richard M.Landslide hazards and landslide hazard mapping in North Carolina[J].Geological Society of America,2006,38(3):28.
- [3] 彭一民.我国灾害学研究蓬勃发展[J].中国地展报,1988(4):1.
- [4] 石菊松,张永双,董 诚,等.基于GIS技术的巴东新城区滑坡灾害危险性区划[J].地球学报,2005,26(3):275-282.
- [5] 唐 川,齐 信.汶川地震高烈度区暴雨滑坡活动的遥感动态分析[J].地球科学,2010(2):317-323.
- [6] 柳 源.中国地质灾害(以崩、滑、流为主)危险性分析与区划[J].中国地质灾害与防治学报,2003,14(1):95-99.
- [7] 张茂省,雷学武.遥感技术在黄土高原区地质灾害详细调查中的应用[J].西北地质,2007(3):92-97.
- [8] 张宏芳,潘留杰.汉阴县近46 年气候变化特征[J].陕西气象,2009(1):23-26.
- [9]吕远强,林杜军,罗伟强.基于人工神经网络的区域地质灾害危险性预测评价[J].中国地质灾害与防治学报,2007,19(1):44-47.
- [10]孙佳伟.略阳县地质风险评价研究[D].西安:长安大学,2012.
- [11]汉阴县志编篡委员会.汉阴志[M].西安:陕西人民出版社,1991.
- [12]王国良.层次分析法在地质灾害危险性评估中的应用[J].西部探矿工程,2006(9):286-288.
- [13]魏海霞.基于MapGIS的永平县地质灾害易发性评价[J].地矿测绘,2007,23(2):34-35.
- [14]唐亚明,薛 强,毕俊擘,等.陕北黄土崩塌灾害风险评价指标体系构[J].地质通报,2012,31(6):979-988.