基于斜坡单元的大区域地质灾害风险性评价——以莲花镇为例

朱 镕，柯 波，叶 武，吴丽斌

（1.武汉理工大学安全科学与应急管理学院，武汉 430070；2.武汉理工大学资源与环境工程学院，武汉 430070；3.浙江海川勘察有限公司，杭州 310000）

我国遭受地质灾害的破坏非常严重，特别是山地因为受到多种因素的影响非常容易发生地质灾害，地质灾害的频繁发生对我国的经济高速高效发展和居民的身体健康造成很大的隐患[1-6]。莲花镇是一个多山的区域，因此开展地质灾害风险评价有着非常重要的意义。

对于地质灾害风险评价的研究，国内外许多学者做了很多有益的研究。国外学者如Cotecchia 等[7]把研究对象瞄准了滑坡，对意大利南部斜坡的滑坡危险性进行分区研究。Brugner 等[8]在1972 年以1∶50 000 的地质图为基础，根据实地情况选取了地貌、坡角等影响因子来分析斜坡稳定性，并且形成了所研究区域的稳定性分区图。Remondo 等[9]在2005年把Deva 做为示例，研究了该区域的地质灾害分区工作，在此次研究中提出了新的地质灾害风险性评价方法，这种新方法是基于统计学原理的。我国地质灾害风险评价比国外发展较晚，开始阶段更多的只是对地质灾害的定性分析，目前我国也进入到定量评价阶段，比如周鑫[10]选用AHP-PCA-ICM 模型与PCA-PSO-SVM 模型对研究区金沙江上游进行地质灾害危险性分区，选用了高差、覆盖层、降雨和覆盖植物种类等8 个因子进行评价。朱良峰等[11]在2002年基于GIS 系统提出了新的进行地质灾害风险性分区的评价方法和评价过程。可以看出，随着研究学者对地质灾害的理论和防治技术的逐步研究，现研究学者对区域地质灾害的研究也更趋向于更精确及大比例尺[12]。基于此本研究通过对莲花镇的孕灾地质条件和承灾体等进行调查，来分析研究莲花镇地质灾害发育类型、形成特征和破坏模式，开展大比例尺下地质灾害风险评价研究，对莲花镇范围进行斜坡单元划分，以斜坡单元为基础来逐个评价该地区地质灾害的易发性、危险性及风险性，为研究区地质灾害防治提供支撑。

1 研究区地质环境条件概况

莲花镇总体地势东北高西南低，海拔高度65.28～1 041.55 m，山脉多分布于西北部呈条带状，形成中低山逐渐向中部丘陵区过渡的地貌形态。除了河流沟谷略平坦外，大部分地形起伏较大，山体自然地形坡度在20~50°。

莲花镇属亚热带季风气候，温暖湿润，雨量充沛，四季分明，光照充足，年平均气温17 ℃，莲花镇最大年降雨量1 852.9 mm，每年6 月中下旬—7 月上旬为集中降雨期，受冷暖气流控制连绵阴雨，常暴雨成灾，该时期灾害性天气较多，易引发地质灾害。

莲花镇位于扬子准地台钱塘台褶带，主要出露二叠系、奥陶系等早古生代地层，以唐家坞组分布最为广泛，石炭系及泥盆系等仅局部地段有零星分布。

2 评价体系构建

目前，国内外研究人员对于地质灾害风险评价是由地质灾害危险性和承灾体易损性组成的，根据联合国提出的风险性=危险性×易损性。本文结合以往研究从4 个方面对地质灾害风险进行评价，分别是易发性、危险性、承灾体易损性和风险性评价。

2.1 评价单元确定

目前国内外地质灾害风险评价比较常用的基本单元有网格单元、流域单元、斜坡单元及地形单元等。其中斜坡单元可以提高与实际地貌的吻合程度，并且斜坡单元是地质灾害发生的基本单元。国内外众多学者利用斜坡单元来评价地质灾害[13]。因此，本次地质灾害风险评价采用斜坡单元作为基本评价单元。本次斜坡单元的划分基于DEM 模型数据，利用Arcgis 水文分析模块来提取沟谷线和山脊线形成斜坡单元，再根据微地形来修正斜坡单元，本次将研究区域总共划分为1121个斜坡单元。

2.2 易发性评价体系构建

地质灾害易发性是在某特定范围内发生地质灾害概率的反映，是由地质环境条件所决定的，是一种自然属性。目前国内外对易发性评价常用方法主要有聚类分析法、人工神经网络法、综合指数法、信息量法和支持向量机等[14-16]。本次研究选取综合指数法进行易发性评价。综合指数法是把会导致地质灾害发生的有关因子进行无量纲化，是将各个因子确定不同的权重得分来分析因子对地质灾害发生的不同的影响程度。再进行加和计算来确定各个斜坡单元的易发性指数，其评价模型公式为

式中：Yi为第i个斜坡单元易发程度综合指数；Fj为第i个斜坡单元第j类指标权重；Sj为第i个斜坡单元第j类指标赋值；n为项数。

对于影响因子权重的确定，本次研究对研究区域地质灾害特征和规律进行分析后确定评价指标因子为坡度、坡向、高差、坡形、覆盖层厚度、岩性与岩土结构、斜坡结构、与构造间距和切坡高度共9 项。采用层次分析法，比较各个影响因子的相对重要性，通过判断矩阵确定9 个影响因子的权重，并且对各个指标进行赋值。最后确定适合本研究区域的地质灾害易发性评价指标体系及量化分值表，具体见表1。

表1 以斜坡为单元的地质灾害易发程度评价指标体系及量化分值表

本次研究根据研究区域各个斜坡单元的评价结果，将地质灾害易发程度等级分为高易发、中易发、低易发和不易发4 个等级。分级标准见表2。

表2 地质灾害易发程度划分标准表

2.3 地质灾害危险性评价方法

易发性评价的局限性在于只研究了影响地质灾害发生的静态因子，对降雨这种动态因子并没有考虑进去。因此危险性评价是在易发性评价的基础上叠加降雨工况进行评价。本次研究采用危险性指数法在易发性评价的基础上对研究区域的危险性进行评价。危险性综合指数法计算方式如下

式中：Hi为某种工况下第i个评价单元危险性指数（危险性概率）；Yi为第i个评价斜坡单元易发性指数；Ymax为最大易发性指数，Ymax=5；Pi为某种工况下第i个评价单元的给定时间段内的失稳概率。

其中对于Pi的计算方法采用极值降雨假设来确定。极值降雨假设是基于研究区历史上有地质灾害发生的事实，假设有监测纪录以来，24 h 最大降雨量Lmax/day为灾害发生的触发因素，不同降雨工况下失稳概率Pi则可表达为Pi=L/Lmax/day。其中L对照4 种工况，分别取24 h 计降雨量35、75、175 和250 mm，Lmax/day的取值则根据搜集研究区及周边区域近几十年以来的最大24 h 降雨量来确定，本次研究取值为190 mm。

根据此评价方法，计算得到每个评价单元的危险性指数值，根据地质灾害危险程度划分标准（表3）对各评价单元地质灾害危险程度等级进行分级。

表3 地质灾害危险程度划分标准表

2.4 地质灾害风险性评价方法

地质灾害风险评价是在危险性评价的基础上叠加上承灾体易损性评价来评价的。承灾体是被地质灾害伤害的如人口、建筑物、线性工程和环境等对象。易损性是指在可能被灾害威胁范围内的承灾体在地质灾害发生时，遭到的伤害程度的预估。本次研究的承灾体分为人员和经济2 种类型。本次风险性评价计算公式如下

式中：Ri为某工况下第i个评价单元风险值；Hj为某工况下第i个评价单元内j号承灾体地质灾害危险性指数；Ej为某工况下第i个评价单元内j号承灾体价值；Vj为某工况下第i个评价单元内j号承灾体易损性。

本次风险性评价中承灾体经济价值是根据建筑结构进行划分的，其价值见表4。

表4 承灾体经济价值

本次承灾体易损性评价分为人员易损性和经济易损性。其中人员易损性按最大风险原则考虑在影响范围之内的人员其易损性值都为1，影响范围之外的人员易损性都为0。经济易损性根据构筑物本身结构确定易损性，其数值区间为0~1。

本次研究将风险等级划分为极高风险、高风险、中风险和低风险4 个等级，按照人员伤亡风险和经济损失风险，综合考虑以就高原则确定其最终风险等级。划分标准见表5。

表5 地质灾害风险等级划分标准表

3 莲花镇地质灾害风险评价

3.1 野外调查

野外调查主要是调查研究区孕灾地质环境条件，识别地质灾害隐患或风险，查明研究区内地形地貌、植被、地质构造、岩（土）体工程地质特征、斜坡结构和残坡积层厚度等孕灾条件和气候、水文、人类工程活动等诱发因素，利用调查的数据对斜坡单元进行评价分析。主要包括对斜坡单元的地面调查和承灾体调查。其中地面调查主要是对易发性评价的影响因子的数据进行收集，根据收集的数据对研究区进行易发性评价。

3.2 易发性评价

根据研究区野外调查搜集到的数据和上述评价体系进行地质灾害易发性评价。据统计，莲花镇易发性评价中，中易发斜坡单元88 个，面积5.92 km2，占全镇总斜坡单元的7.85%，占全镇总面积的6.87%；低易发斜坡单元1 002 个，面积77.76 km2，占全镇总斜坡单元的89.38%，占全镇总面积的90.29%；不易发斜坡单元31个，面积2.44 km2，占全镇总斜坡单元的2.77%，占全镇总面积的2.83%。本研究区地质灾害易发性较高区域都是人类工程活动严重、人工切坡高、岩石较破碎且土层较厚的区域。低易发区域主要是人类工程活动弱、土层覆盖薄、植被茂盛和岩石较完整的区域。不易发区域主要是平原和河流区域。

3.3 危险性评价

根据上述危险性评价体系及方法，计算各个斜坡单元不同工况下的地质灾害危险性指数，对其分级后得到莲花镇不同工况下地质灾害危险程度评价结果。评价结果显示，在大雨、暴雨工况下，莲花镇全镇范围均为低危险性。在大暴雨工况下，莲花镇地质灾害危险程度划分为高危险、中危险、低危险3 级，其中高危险斜坡单元51 个，面积3.99 km2，占全镇总斜坡单元的4.55%，占全镇总面积的4.63%；中危险斜坡单元1 013个，面积78.11 km2，占全镇总斜坡单元的90.37%，占全镇总面积的90.70%；低危险斜坡单元57 个，面积4.02 km2，占全镇总斜坡单元的5.08%，占全镇总面积的4.67%。在特大暴雨工况下，莲花镇地质灾害危险程度划分为高危险、中危险、低危险3 级，其中高危险斜坡单元89 个，面积5.95 km2，占全镇总斜坡单元的7.94%，占全镇总面积的6.91%；中危险斜坡单元995个，面积77.43 km2，占全镇总斜坡单元的88.76%，占全镇总面积的89.91%；低危险斜坡单元37 个，面积2.74 km2，占全镇总斜坡单元的3.30%，占全镇总面积的3.18%。

3.4 风险性评价

根据地质灾害风险评价计算方法，在确定不同降雨工况条件下区内地质灾害危险性评价的基础上，通过分析地质灾害危险性评价指数，耦合承灾体人口和经济价值，采用定量计算的方法进行地质灾害风险评价。

综合不同工况下人员和经济伤亡风险在就高原则下的风险性评价结果为：在大雨工况下，区内综合风险2 处中风险，其余为低风险；在暴雨工况下，区内综合风险1 处高风险，6 处中风险，其余均为低风险；在大暴雨工况下，区内分布3 处高风险，20 处中风险，其余为低风险；在特大暴雨工况下，区内分布4 处高风险，21 处中风险，其余为低风险。

4 结束语

1）根据对莲花镇地质灾害特征及分布发育规律的分析研究，结合实际情况选取9 个评价指标构成了莲花镇地质灾害易发性评价指标体系。使用层次分析法确立评价指标权重，并利用综合分析法评价区域易发性，将莲花镇分为了中易发、低易发、不易发3 个等级，其中中易发区面积5.92 km2，占全镇总面积的6.87%；低易发区面积77.76 km2，占全镇总面积的90.29%；不易发区面积2.44 km2，占全镇总面积的2.83%。

2）根据易发性，在不同的降雨工况下，评价了各斜坡的危险性。在特大暴雨工况下，莲花镇地质灾害危险程度划分为高危险、中危险、低危险3 级，其中高危险区面积5.95 km2，占全镇总面积的6.91%；中危险区面积77.43 km2，占全镇总面积的89.91%；低危险区面积2.74 km2，占全镇总面积的3.18 km2。

3）根据危险性数据再叠加影响范围内的承灾体信息，综合经济和人口损失在最不利损失情况下得出，莲花镇地质灾害风险程度划分为高风险、中风险、低风险3 级，其中高风险4 处，中风险21 处，其余为低风险。