洪水灾害是当今世界最为频发的自然灾害之一,而洪水灾害风险评价是防洪减灾的重要支撑。本文建立浑太胡同防洪保护区溃堤洪水灾害风险评价的方法框架,并对研究区因典型溃堤洪水导致的洪灾风险进行评价。首先,通过构建浑太胡同溃堤洪水演进二维数值模型,模拟分析了研究区内溃堤洪水的时空演变过程;其次,以栅格为评价对象,综合考虑研究区的实际情况,建立了溃堤洪水风险评价指标体系,并通过级差最大化合成法组合了指标的层次分析法权重和熵权法权重,得到了各评价指标的组合权重;最后,利用K-Means聚类算法对评价对象综合得分进行了聚类分析,最终得到了研究区发生典型溃堤洪水时的风险情况。结果表明,级差最大化合成法能够对研究区洪水风险进行合理区分,且能避免乘法合成法可能出现的“倍增效应”;当迟拓子与马家堡子组合溃口发生100年一遇的流量过程时,中等和较高洪灾风险区域占研究区总面积的89.64%。 Floods are considered as one of the most frequent natural disasters and flood risk assessment is an important support for flood prevention and disaster reduction. This study establishes a methodological framework for evaluating the flood risk of Hun-Tai flood protection zone due to the typical floods. Firstly, by building a two-dimensional numerical model of flooding, the spatial and temporal evolution process of dike-break flood in the study area was simulated and analyzed. Secondly, taking grids as research objects and considering the actual situation, a risk evaluation index system was established, and the subjective weights of the indices calculated by the AHP and the objective weights of the indices calculated by the entropy weight method were combined through the level difference maximization method (LDMM). Finally, the K-Means clustering algorithm was used to perform a cluster analysis on the comprehensive scores of the evaluation objects, and the risk situation of a typical dike-breaking flooding was obtained. The results of the study show that, the LDMM can reasonably evaluate the flood risk, and can avoid the “multiplication effect” that may occur in the multiplicative synthesis method; when a 100-year flow process occurs at the combination of Chituozi and Majiapuzi dike-break, the moderate and higher flood risk areas cover 89.64% of total area.
穆喜平1,王运涛1*,刘海星1,王睿1,张弛1,彭勇1,杨坤2
1大连理工大学水利工程学院,辽宁 大连
2上海城市水资源开发利用国家工程中心有限公司,上海
收稿日期:2020年4月10日;录用日期:2020年5月13日;发布日期:2020年6月2日
洪水灾害是当今世界最为频发的自然灾害之一,而洪水灾害风险评价是防洪减灾的重要支撑。本文建立浑太胡同防洪保护区溃堤洪水灾害风险评价的方法框架,并对研究区因典型溃堤洪水导致的洪灾风险进行评价。首先,通过构建浑太胡同溃堤洪水演进二维数值模型,模拟分析了研究区内溃堤洪水的时空演变过程;其次,以栅格为评价对象,综合考虑研究区的实际情况,建立了溃堤洪水风险评价指标体系,并通过级差最大化合成法组合了指标的层次分析法权重和熵权法权重,得到了各评价指标的组合权重;最后,利用K-Means聚类算法对评价对象综合得分进行了聚类分析,最终得到了研究区发生典型溃堤洪水时的风险情况。结果表明,级差最大化合成法能够对研究区洪水风险进行合理区分,且能避免乘法合成法可能出现的“倍增效应”;当迟拓子与马家堡子组合溃口发生100年一遇的流量过程时,中等和较高洪灾风险区域占研究区总面积的89.64%。
关键词 :溃堤洪水,风险评价,HEC-RAS,指标体系,级差最大化
Copyright © 2020 by author(s) and Wuhan University.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
《减轻灾害风险全球评估报告》指出 [
近年来,国内众多学者在洪灾风险评价方面做了大量的研究 [
本文在现有研究的基础上,旨在从致灾因子、孕灾环境和承灾体三方面来构建综合考虑情景模拟结果的基于指标体系的溃堤洪水风险评价方法,并研究评价指标权重的计算方法,探索K-Means聚类算法在洪水风险等级区划中的应用,最后以此作为方法框架对浑太胡同进行溃堤洪水风险评价研究。结果表明,该方法方能够对洪灾风险进行有效的评价。
浑太胡同防洪保护区位于辽宁省,“浑太”指的是浑河和太子河,如图1所示,自小北河起到下游的三岔河,浑河和太子河大体呈平行状向西南流去,在两河的堤防间形成了一块封闭的类似胡同状的区域(以下简称:浑太胡同),该区域的面积约为300 km2。
图1. 浑太胡同防洪保护区及河流示意图
浑太胡同涉及到的乡镇级行政单位包括辽阳县境内的小北河镇、唐马寨镇及海城市境内的高坨镇、温香镇。根据2012年的统计数据 [
区域名称 | 所属县区 | 区域面积(km2) | 常住人口(人) | 耕地面积(公顷) | 农业产值(万元) | 工业产值(万元) | 服务业产值(万元) |
---|---|---|---|---|---|---|---|
小北河镇 | 辽阳县 | 134 | 45,937 | 8000 | 56,000 | 115,900 | 284,630.07 |
唐马寨镇 | 139 | 41,287 | 7936 | 48,171.52 | 45,259 | 143,027.80 | |
高坨镇 | 海城市 | 59.33 | 24,532 | 3466.67 | 24,595 | 134,935 | 288,196.66 |
温香镇 | 140.27 | 23,854 | 13,779.42 | 97,761 | 0 | 280,231.66 |
表1. 浑太胡同防洪保护区社会经济数据
模型输入的数据主要有:DEM数据、糙率数据和入流边界条件。研究区的地形数据来源于松辽委,如图1所示,DEM的分辨率为30 m,本文对研究区域进行了网格剖分,共划分了大小平均约为10,000 m2的网格31,273个(下文所指的评价对象);参考对研究区洪水模拟的已有成果及研究区的实际情况,将居民地、水田、旱田、树丛和空地5种不同土地利用类型的糙率值分别设置为0.070、0.050、0.060、0.065和0.035;对于入流边界条件,根据研究区堤坝和洪水特点,本次风险评价选取迟拓子与马家堡子组合溃口作为典型溃口,溃口位置如图1所示,并将组合溃口100年一遇的流量过程作为模型的入流边界条件。各溃口的流量过程如图2所示:
图2. 组合溃口100年一遇流量过程
本文基于指标体系和情景模拟的风险评价方法,首先,利用HEC-RAS软件构建研究区溃堤洪水演进模型来获取致灾因子指标,通过灾损评估得到承灾体洪灾经济脆弱性指标,并利用ArcGIS处理研究区地形数据得到孕灾环境指标;其次,由级差最大化合成法对层次分析法权重和熵权法权重进行优化组合,得到各评价指标的组合权重,并计算各评价对象的综合得分;最后,通过K-Means聚类算法对各评价对象的综合得分进行聚类分析,可得到研究区的洪水风险图。洪水灾害风险评价的思路如图3所示:
图3. 洪水灾害风险评价
蓄滞洪区的洪水要素分析方法可分为水文学法和水动力学法 [
连续方程:
∂ H ∂ t + ∇ ⋅ h V + q = 0 (1)
动量方程:
∂ V ∂ t + V ⋅ ∇ V = − g ⋅ ∇ H + υ 1 ∇ 2 V − c f V + f k × V (2)
式中,H为水面高程(m);t为时间(s);h为水深(m);V为水流速度(m/s);q为旁侧入流量(m2/s);g为重力加速度(m/s2); υ 1 为运动粘度(m2/s); c f 为摩擦系数;f为科里奥利系数;k为单位矢量。
HEC-RAS作为一款水动力学软件,其2D模块便是基于浅水方程原理开发的。该软件的主要功能包含以下四个方面:① 恒定流水面计算;② 1D和2D非恒定流模拟;③ 可变边界的泥沙输移计算;④ 水温和水质传输模拟。利用以上功能,HEC-RAS既可以模拟桥梁、溢洪道、堰、闸门、滞洪池等水工建筑物对水流的影响,也可以对洪水进行二维演进数值模拟。
1) 指标体系
根据研究区的特点,参考国内外关于洪水灾害风险评价的研究成果,本文在选择评估指标时,综合考虑了洪水对自然、经济、社会等方面的影响,从致灾因子、孕灾环境和承灾体脆弱性三个方面出发,确定了用以评价研究区溃堤洪水风险的指标体系。
指标体系可以划分为目标层、准则层和指标层,其中,目标层是以浑太胡同溃堤洪水风险评价为总目标,准则层包含致灾因子、孕灾环境和承灾体。致灾因子指标包括:最大淹没水深(MD)、洪水到达时间(AT)、洪水淹没时间(DT);孕灾环境指标包括:地表高程(DEM)、地形坡度(SL);承灾体指标:洪灾GDP损失(LOP)。根据指标与洪水灾害风险的相互关系,可以将指标划分为正向指标或负向指标,风险随正向指标的增大而增加,随负向指标的增大而减小。指标体系如表2所示:
目标层 | 准则层 | 指标层 | 趋向性 |
---|---|---|---|
浑太胡同溃堤 洪水风险评价 | 致灾因子 | 最大淹没水深(MD) | + |
洪水淹没时间(DT) | + | ||
洪水到达时间(AT) | - | ||
孕灾环境 | 地表高程(DEM) | - | |
地形坡度(SL) | + | ||
承灾体 | GDP损失(LOP) | + |
表2. 风险评价指标体系
2) 指标权重
指标权重的计算方法可以分为主观赋权法、客观赋权法和主客观组合赋权法 [
级差最大化组合赋权模型以各评价对象的区分度最大为目标函数、以各指标权重取值区间与指标权重和等于1为约束条件来对各指标的权重进行优选,其中各评价对象的区分度可以用其综合得分的方差s2来度量。记k个评价对象的综合得分为 x = ( x 1 , x 2 , ⋯ , x k ) ,k个评价对象n个评价指标的标准化矩阵为X,n个评价指标的权重为 θ = ( θ 1 , θ 2 , ⋯ , θ n ) , x ¯ 为各评价对象得分的平均值,则有:
X = ( x 11 x 12 ⋯ x 1 k x 21 x 22 ⋯ x 2 k ⋮ ⋮ ⋱ ⋮ x n 1 x n 2 ⋯ x n k ) (3)
第i个指标的权重 θ i ( i = 1 , 2 , ⋯ , n ) 的取值范围 [ θ i − , θ i + ) 由层次分析法计算所得的权重和熵权法计算所得的权重共同确定,进而,组合权重优化模型的数学表达式如(4)所示:
s 2 = max 1 k − 1 ∑ i = 1 k ( x i − x ¯ ) 2 s . t . { ∑ i = 1 n θ i = 1 θ i − ≤ θ i < θ i + (4)
K均值聚类算法 [
S S E = ∑ i = 1 k ∑ p ∈ C i | p − m i | 2 (5)
式中,k为聚类数目, C i 为第i个样本簇,p为 C i 的样本点, m i 为 C i 的质心。SSE随k值变化的规律是:k取值越大时,样本数据便被划分的越细,每个样本簇的聚合程度就越高,进而误差平方和SSE变小。根据k与SSE的这种关系便可确定K-Means的聚类数目。因k与SSE的关系曲线形似手肘,故而该方法被称为手肘法。
2016年,中国水利水电科学研究院等单位对浑太地区的洪水风险进行了研究。本文根据水科院等单位的研究成果报告 [
通过HEC-RAS对研究区的溃堤洪水演进进行模拟,可直接获得致灾因子指标MD、AT和DT。基于研究区洪水淹没计算结果,结合不同财产的灾前价值,根据分类财产损失率与水深的关系 [
图4. 迟拓子与马家堡子组合溃口发生100年一遇流量过程时研究区各指标空间分布
为了计算层次分析法权重,本文通过参考其它相关研究成果和专家的意见,对比了各评价指标的相对重要程度,形成了用以计算指标层次分析法权重的判断矩阵,判断矩阵如表3所示。
指标 | DEM | SL | MD | AT | DT | LOP |
---|---|---|---|---|---|---|
DEM | 1 | 4 | 1/3 | 2 | 1/3 | 1/3 |
SL | 1 | 1/6 | 1/2 | 1/6 | 1/5 | |
MD | 1 | 5 | 3 | 1/2 |
AT | 1 | 1/3 | 1/5 | |||
---|---|---|---|---|---|---|
DT | 1 | 1/5 | ||||
LOP | 1 |
表3. 层次分析法判断矩阵
由各种权重计算方法求得的评价指标权重结果见表4。从表中的计算结果可以看出,由层次分析法和熵权法计算所得的LOP指标权重均排序第一,熵权法计算得到的SL指标权重远高于由层次分析法计算所得的权重;级差最大化合成法综合考虑了主、客观权重两方面的因素,所得权重中,LOP、MD指标权重分别排序第一、第二。
(1)准则 | (2)指标 | 单一权重 | 组合权重 | |||||||
---|---|---|---|---|---|---|---|---|---|---|
AHP权重 | 熵权重 | 级差最大化组合权重 | ||||||||
(3)权重 | (4)排名 | (5)权重 | (6)排名 | (7)组合权重区间 | (8)权重 | (9)排名 | ||||
致灾因子 | MD | 0.2651 | 2 | 0.1401 | 3 | (0.1401,0.2651) | 0.2017 | 2 | ||
AT | 0.0594 | 5 | 0.1324 | 4 | (0.0594,0.1324) | 0.1165 | 4 | |||
DT | 0.1629 | 3 | 0.1324 | 5 | (0.1324,0.1629) | 0.1508 | 3 | |||
孕灾环境 | DEM | 0.1057 | 4 | 0.1297 | 6 | (0.1057,0.1297) | 0.1103 | 5 | ||
SL | 0.0393 | 6 | 0.1752 | 2 | (0.0393,0.1752) | 0.0987 | 6 | |||
承灾体 | LOP | 0.3676 | 1 | 0.2902 | 1 | (0.2902,0.3676) | 0.3220 | 1 | ||
表4. 风险评价指标权重计算结果
图5为各评价对象的综合得分经手肘法分析后的结果。从图中可以看出,对于级差最大化合成法来说,当k大于4时,随着k值的增加,SSE的下降幅度越来越小。为了方便区划研究区洪水风险等级,可将K-Means聚类数目k定为5,也即研究区溃堤洪水的风险等级可划分为5级:I最高、II较高、III中等、IV较低、V最低。将k和存有各评价对象综合得分的数据文件输入到K-Means计算程序中进行聚类分析,由分析结果可知,当迟拓子与马家堡子组合溃口发生100年一遇的溃口流量时,各级洪水风险的取值范围分别为:I (0.5012, 1)、II (0.3812, 0.5012)、III (0.3369, 0.3812)、IV (0.2591, 0.3369)、V (0, 0.2591)。当迟拓子与马家堡子组合溃口发生100年一遇溃口流量过程时,图6中的(a)为对应求得的研究区洪水风险图。
图5. SSE与k的关系曲线
图6. 由不同权重计算法求得的研究区洪水风险图
1) 评价结果合理性
图6中(b)和(c)分别为由层次分析法和熵权法求得的研究区洪水风险图。对比图6中的(a)和(b),可以看出由级差最大化组合权重求得的风险图和由层次分析法所求得的风险图整体趋势基本相同,但(b)图在研究区上游的中等、较高风险区域面积小于(a)图的;对比图6中的(a)和(c),可以看出(a)和(c)的上半部分基本相同,但下半部分差别较大;而图6中的(b)和(c)整体差别较大。由此可见,利用级差最大化合成法组合层次分析法权重和熵权法权重,求各评价指标的组合权重,能够改善单一权重法的不足,使得评价结果能够兼顾主、客观两方面的因素。
2) 洪水风险分析
从图6中的(a)可以看出,研究区溃堤洪水风险整体呈现由上游到下游增大的趋势,这与溃堤洪水水深由上游到下游呈增大的趋势一致;研究区洪水风险等级最高的区域多为居民居住地,较高风险区域多为下游受淹没的水田和旱田,而风险最低的区域多处于研究区上游,这主要与研究区上游地势较高、水深较浅有一定的关系。表5总结了研究区各级洪水风险区域的面积和占比情况,从表可以看出,中等风险区域和较高风险区域的面积和占到了研究区总面积的89.64%。
风险等级 | 最低 | 较低 | 中等 | 较高 | 最高 |
---|---|---|---|---|---|
所占面积 | 6.83 | 20.83 | 93.43 | 189.58 | 5.04 |
所占比例 | 2.16% | 6.60% | 29.59% | 60.05% | 1.60% |
表5. 研究区各风险等级区域统计(面积单位:km2)
3) 评价指标分析
图7为不同风险级别下各指标的平均取值,图8为图6中(a)所示区域各指标的取值情况,各区域均由程序随机选取得到,其中,1和2为最低风险区域,3和4为较低风险区域,5和6为中等风险区域,7和8为较高风险区域,9和10为最高风险区域。从图7可以看出,风险最高区域的洪灾经济损失最高,从前边分析知,这些区域多是人们生产和生活的聚集地,在研究区内呈现点状分布。从图7和图8可以看出,研究区整体的洪水风险等级、不同区域的洪水风险等级不依某一单一指标取值的增减而呈现增大或减小的趋势,体现了不同区域的风险等级是各评价指标综合作用的结果,比如,虽然区域4、5的DEM取值相同,但二者的风险级别不同。
图7. 不同风险等级下各指标的平均取值
图8. 代表区域各指标取值
4) 权重组合方法对比
通过乘法合成法 [
本文利用HEC-RAS模型对研究区当迟拓子与马家堡子组合溃口发生100年一遇溃堤流量时进行了洪水演进模拟,构建了综合考虑情景模拟结果的基于指标体系的溃堤洪水风险评价方法,利用级差最大化合成法,求得了用以评价研究区溃堤洪水风险的各评价指标的组合权重,并利用K-Means聚类算法对研究区洪水风险进行了聚类分析。得出以下结论:
1) 当迟拓子与马家堡子组合溃口发生100年一遇的溃口流量过程时,研究区的洪水风险整体呈现由上游到下游增大的趋势,风险最高的地方多为居民聚居地,中等风险和较高风险区域占到了研究区总面积的89.64%。
2) 利用乘法合成法组合主、客观权重时,可能会出现倍增效应,使得某一指标的组合权重远高于由层次分析法或熵权法所得的该指标的权重,因此弱化了其他指标在评价体系中的作用。
3) 由级差最大化合成法得到的组合权重能够对评价对象进行有效区分,且各区域的风险高低不依某一单一指标的增减而降低或升高,而是各指标综合作用的结果。
该研究是在中央高校基本科研业务费(DUT18RC(3)072)的资助下完成的。
穆喜平,王运涛,刘海星,王 睿,张 弛,彭 勇,杨 坤. 基于HEC-RAS的浑太胡同溃堤洪水灾害风险评价Flood Risk Assessment of Hun-Tai Area Due to Dike Break Based on HEC-RAS[J]. 水资源研究, 2020, 09(03): 289-300. https://doi.org/10.12677/JWRR.2020.93030