摘 要:针对城市洪涝灾害韧性评价过程主观性强、未考虑决策者有限理性和后悔规避的心理行为导致评价结果失真的难题,提出基于有序加权平均算子和后悔理论的城市洪涝灾害韧性评价方法。首先,结合韧性城市概念,从抵抗能力、预警能力、适应能力、恢复能力4 个维度构建指标体系,利用有序加权平均算子将专家极端决策值分布到权重较低位置,通过组合权重的形式求得指标权重,弱化主观偏好极端评价值对权重的影响;然后,考虑专家决策时后悔心理行为特征,构建效用值矩阵、后悔-欣喜值矩阵和感知效用矩阵,得出城市洪涝灾害韧性评价值;最后,将其运用在河南省洪涝灾害韧性评价中,并将计算结果与其他评价方法得出的结果对比分析。结果表明:4 种方法评价结果基本一致,均认为河南省洪涝灾害韧性等级为Ⅲ级,提出的方法评价结果精度更高、更贴近实际情况。此外,各地级市韧性等级差异较大,分布不均,整体呈现西南部城市韧性等级高于东北部城市的特征。
关键词:城市洪涝灾害;韧性评价;有序加权平均算子;后悔理论;河南省
中图分类号:TU998.4 文献标志码:A doi:10.3969/ j.issn.1000-1379.2024.11.009
引用格式:李亚敏,张楠.基于有序加权平均算子和后悔理论的城市洪涝灾害韧性评价[J].人民黄河,2024,46(11):56-62.
在全球气候变暖和极端天气频繁发生的背景下,城市公共安全经常被各类自然灾害威胁。据统计,城市洪涝灾害在自然灾害中占比高达56.8%[1] 。如近年来的济南“7·18”大暴雨、北京“7·21”特大暴雨、郑州“7·20”特大暴雨事件,均对城市公共安全造成巨大冲击,带来巨大经济损失和人员伤亡[2] 。为更好地应对城市洪涝灾害,城市韧性概念随之被提出,使得人们对城市洪涝灾害有更加深刻的认知。相对于传统防灾减灾评价,城市韧性视角下的洪涝灾害评价综合考虑城市洪涝防灾、适应和灾后恢复能力[3] 。
关于城市洪涝灾害韧性评价,国内外学者已有诸多研究成果。Miguez 等[4] 采用洪涝风险指数评价城市的防洪方案优越性;Bertilsson 等[5] 利用空间洪水韧性指数衡量不同城市抗洪涝能力;黄晶等[6] 从自然条件、基础设施、社会经济3 个方面阐述影响城市洪涝灾害韧性指标,采用系统动力学方法预测城市洪涝灾害韧性等级;周铭毅等[7] 采用熵权法和VIKOR 方法评价广东省21 个地级市洪涝灾害韧性等级;蒋硕亮等[8] 基于洪涝背景定性分析影响河南省18 个地级市城市韧性的灾害指标;孟晓静等[9] 采取组合赋权方式确定指标权重,利用TOPSIS 模型评价城市洪涝灾害韧性等级。上述研究在一定程度上推动了城市洪涝灾害韧性评价,但均忽视了专家主观偏好可能造成指标权重失真的问题,且评价结果基于专家在完全理性状态下给予的合理判断。为解决专家主观偏好造成权重失真问题,李志荣等[10] 、董兴芝等[11] 利用有序加权平均算子将专家极端决策值分布在权重较低位置,合理优化极值对权重的干扰。此外,已有研究[12-13] 表明,受外部条件和信息不完整以及认知局限性等因素限制,专家决策结果难以体现专家完全理性,往往呈现规避损失、敏感性递减、概率判断扭曲等心理特征,所做出的评价结果基于主观认知,未能充分体现客观性。因此,为提升决策结果的客观、真实性,如实反映城市洪涝灾害韧性等级,需将富有主观性的评价依据转化为客观依据。后悔理论考虑到专家评分过程中主客观界面的不确定性,通过选择理想评价值作为主客观界面划分的依据,利用效用化后的评分值与理想值之间的差异阐述专家因认知有限理性和心理作用等因素造成的主客观差异,通过修正该差异值实现主观依据客观化的目的,使评价结果如实反映现实情况[14] 。
本文提出基于有序加权平均算子和后悔理论的城市洪涝灾害韧性评价方法,首先利用有序加权平均算子将主观性过强的评价值分布到权重较低位置,使得指标权重仅与所处位置相关,以此弱化专家主观偏好过强对指标权重的影响;然后,通过后悔理论进一步克服专家决策的主客观差异性,使得决策结果更贴近实际情况,提升决策质量;最后,将该方法运用于河南省城市洪涝灾害韧性评价,并与其他方法评价结果进行对比。
1 后悔理论
1982 年,Bell 在前景理论基础上提出后悔理论,以决策者有限行为为前提刻画决策者评价目标时的心理偏好,充分考虑了决策者的心理因素和后悔规避行为,使得决策者对决策结果的后悔程度降到最低[15] 。受决策者心理因素影响,在多方案评价时不由自主地将自身所选方案与剩余方案对比,如果前者优于后者,决策者内心充满欣喜;反之,内心充满后悔。鉴于此,后悔理论可帮助决策者规避后悔心理,提升所选择方案的优越性,使之更加贴近理想方案。换言之,在城市洪涝灾害韧性评价中,由于部分指标无法用具体数值衡量,因此需借助专家经验给予评价,而专家决策结果具有一定主观性,如不加以处理则导致评估结果与实际情况存在差异甚至背道而驰。这种差异在后悔理论中表现为后悔程度,差异性越大,专家后悔心理越强,两者正相关。
后悔理论可用感知效用函数表示,由选择目标方案的效用函数和选择此方案的后悔-欣喜函数两部分组成。设x 和y 分别是选择方案A 和方案B 的结果,选择方案A 的感知效用函数为
q(x,y)= v(x)+R[v(x)-v(y)] (1)
式中:v(x)和v(y)分别为专家选择方案A 和B 得到的效用函数;R[v(x)-v(y)]为后悔-欣喜函数,当R[v(x)-v(y)]>0 时,表示专家为选择方案A 和不选方案B 内心感到欣喜,反之表示内心感到后悔。
对于m 个不同方案,专家的感知效用函数为
qi =v(xi )+R[v(xi )-v(x∗)] (2)
式中: x∗ =max{xi },i =1,2,…,m。
由此可知,R[v(x)-v(x∗)]<0,即后悔-欣喜函数值小于0,此时决策者内心最为后悔,使得规避后悔后的选择结果最能反映实际情况。
2 基于有序加权平均算子和后悔理论的城市洪涝灾害韧性评价模型构建
2.1 构建指标体系
构建城市洪涝灾害韧性评价指标体系的前提是明确其定义,文献[7]中定义如下:城市洪涝灾害韧性是指整个城市系统在被洪涝灾害侵害时,通过原有系统的抵抗能力、预警能力、适应能力降低洪涝灾害损失,同时考虑城市的灾后恢复重建能力,使城市居民生活尽快恢复到灾前水平的能力。从定义中不难发现,城市洪涝灾害韧性评价可从抵抗能力、预警能力、适应能力和恢复能力4 个维度展开。抵抗能力是指城市原有系统能够抵御洪涝灾害的能力,指标主要包括大型商业综合体密度、第三产业产值占比、防汛规划与应急预案等[6] 。预警能力是指提前发现洪涝灾害发生并及时告知广大市民的能力,与气象台站覆盖率、水文站覆盖率、广播电视覆盖率、移动电话用户覆盖率、互联网用户覆盖率关联[7] 。适应能力是指城市遭遇洪涝灾害后自我调节维持稳定的能力,与排水管网密度、人均粮食产量、弱势人口比重、救援响应时间、人均水果产量等有关[8] 。恢复能力则指城市在洪涝灾害后采取相关措施恢复原有城市面貌的能力,与洪涝灾害专项资金、植被覆盖率、登记失业率等相关。结合城市洪涝灾害韧性定义、文献[6-9] 成果、专家意见及统计年鉴相关资料,建立表1 所列的城市洪涝灾害韧性评价指标体系。在城市洪涝灾害韧性等级划分方面,通过查阅《城市内涝风险普查技术规范》(GB/ T 39195—2020)、《气象灾害风险调查技术规范》(DB11/ T 1589—2023)等资料,参考文献[7-8]城市洪涝灾害韧性等级划分原则,结合河南省洪涝灾害韧性特点,将洪涝灾害韧性划分为5 个等级,即Ⅰ级、Ⅱ级、Ⅲ级、Ⅳ级、Ⅴ级,分别对应低韧性、较低韧性、中等韧性、较高韧性、高韧性,具体见表1。
2.2 计算指标权重
有序加权平均算子的核心思想是通过将专家主观偏好极端值(极大值和极小值)调整到权重较低位置,按照降序规则对决策数据二次排序,将指标决策数据与所处位置相关联,加权处理得到指标权重。该方法不仅能充分体现专家决策意志,而且能够合理优化指标权重,有效降低主观偏好带来的权重失真,使得指标权重更具说服力和客观性[16] 。考虑到城市洪涝灾害韧性评价指标多数不易量化的特征,需借助专家主观经验给出合理评判,本文选择有序加权平均算子计算指标权重,具有一定的可行性和契合性。具体步骤如下。
1) 构建新数据bi 。设城市洪涝灾害韧性评价指标为a,ai(i = 1,2,…,n) 为第i 个指标的得分,得到指标原始得分决策数据(a1,a2,…,an )。对n 个指标的决策数据从1 开始重新编号,遵循降序原则,得到新数据b1 > b2 > … > bn [17] 。
2)构建加权向量η 。其元素为
式中: Cin 为第i 个专家的组合数。
2.3 构建基于后悔理论的城市洪涝灾害韧性评价模型
1) 构建评价指标矩阵Q。设有m 个方案、n 个指标,根据指标得分情况构建评价指标矩阵Q:
式中:qij 为方案i 中指标j 的得分。
2) 构建理想点矩阵P。设P 为理想点矩阵,P =(p1,p2,…,pn ),pn 为指标n 的理想值。一般而言,为降低决策时专家后悔心理,取矩阵Q 中各方案的最小值作为理想值[18] 。
3)构建效用值矩阵H。选择合理的效用函数关乎各方案效用值的科学性,参考文献[12]和多轮测试结果,选择幂函数作为效用函数能取得最佳效果且满足后悔理论的假设要求。
h(qij ) = (qij )α (7)
式中: α 为效用函数h 的参数, 结合文献[13],其取值范围为(0,1)。
设H 为效用值矩阵,则
H = (qij )α (8)
4)构建感知效用矩阵D。设R 为后悔-欣喜函数值矩阵,则
R = (amp )ij = 1 - e-β(amp)ij (9)
式中:β 为后悔-欣喜函数的参数,结合文献[14],取值大于0,值越小表明专家心理后悔概率越低;(amp)ij为方案m 和理想点矩阵P 的效用差。
根据效用值矩阵和后悔-欣喜函数值矩阵,求得感知效用矩阵D:
D = (dij ) = H + R = (qij )α + 1 - e-β(amp)ij (10)
5)计算综合评价值。集合有序加权平均算子求得的指标权重和感知效用矩阵,得到城市洪涝灾害韧性综合评价值,即
式中:S 为城市洪涝灾害韧性综合评价值。
6)城市洪涝灾害韧性评价流程。如上所述,城市洪涝灾害韧性评价流程见图1。
3 案例分析
河南省位于我国中东部、黄河中下游地区,整体地势西高东低,属于北亚热带向暖温带过渡的大陆性季风气候区,年均降水量407.7 ~1 295.8 mm,主汛期为6—8 月。2022 年河南省生产总值为6 万亿元,常住人口接近1 亿。主要自然灾害为干旱、洪涝、冰雹、滑坡等,以旱、涝两种灾害危害最大[19] 。
3.1 数据来源
本次评价数据源自2017—2022 年河南省及各地级市统计年鉴、水利年鉴、环境公报以及水务局、气象局等提供的相关材料。
3.2 指标权重
考虑城市洪涝灾害韧性评价指标特点,采取问卷调查形式,结合专家打分法评价指标重要性。为提升专家打分结果的可信性,引入“专家可信度”指标对打分结果初步处理[20] 。在专家选择上主要与其职称、学历、工作年限3 个指标关联。利用SPSS 23.0 软件检测打分数据的信度和效度。首先检测打分结果的相关性,其KMO 值为0.865,显著水平为0,说明打分数据强关联,可进行下一步的因子分析[21] ;然后验证各因子的Cronbach’s α、CR 值(组合信度)和AVE 值(平均方差萃取),其中Cronbach’s α>0.70,CR>0.71,AVE>0.52,说明各因子具有较强的信度和收敛度,专家打分结果可作为赋权的依据[22] 。
邀请拥有丰富经验的城市洪涝灾害韧性评价研究专家、学者对各指标按照10 分制评判指标的重要性,得分越高表明指标越重要。以指标A4 为例,专家打分结果、专家基本情况和可信度分别见表2 和表3。
对A41指标得分按照降序规则重新组合,得到新数据序列b =(9.0,8.5,8.0,7.5,6.5)。根据式(3)、式(4)得到绝对权重w41 =(0.156,0.313,0.313,0.156,0.013) ×(9.0,8.5,8.0,7.5,6.5)T =7.94,则A42 ~A45指标绝对权重: w42 = 7.06, w43 = 5.75, w44 = 6.69, w45 = 7.92。
利用式(5) 计算相对权重: w∗41 = 0. 225, w∗42 =0.200,w∗43 =0.163,w∗44 = 0.189,w∗45 = 0.224。同理,可得到一级指标及A1 ~A3 对应的二级指标相对权重:w∗ =(0.325,0. 208,0. 289,0. 178),w∗1 = (0. 169,0.188,0.179,0.229,0.235),w∗2 =(0.179,0.241,0.159,0.197,0.224),w∗3 =(0.217,0.197,0.143,0.163,0.176,0.104)。
3.3 河南省洪涝灾害韧性综合评价
1) 确定评价指标初始矩阵。以A4 恢复能力指标为例,构建评价矩阵QA4 :
2)构建理想点矩阵。结合理想点矩阵确定原则,根据初始矩阵QA4 ,得到A4 恢复能力的理想点矩阵PA4 :
PA4 = [4 5 5 4 4]T
3)构建效用矩阵。由式(8)得到A4 恢复能力的效用值矩阵HA4 :
由式(11)得到指标A1、A2、A3、A4 的综合评价值:SA 1 = 5.84,SA 2 = 5.16,SA 3 = 4.28,SA 4 = 5.78。通过集合指标权重得到河南省洪涝灾害韧性综合评价值S(河南省) = 5.14。
同理,可以得到河南省各地级市的洪涝灾害韧性综合评价值和韧性等级(见图2)。
依据洪涝灾害韧性等级可把河南省各地级市划分为4 类,具体如下。
第一类为较高韧性城市,韧性等级为Ⅳ级:S(郑州市)= 6.73。
第二类为中等韧性城市,韧性等级为Ⅲ级:S(商丘市)= 5.34,S(驻马店市)= 5.17,S(洛阳市)= 4.76,S(南阳市)= 4.65。
第三类为较低韧性城市,韧性等级为Ⅱ级:S(开封市)= 3. 98,S(周口市) = 3. 87,S (信阳市) = 3. 52,S(新乡市) = 3. 29,S (焦作市) = 3. 27,S (安阳市) =3.13,S(平顶山市)= 2.89,S(许昌市)= 2.74。
第四类为低韧性城市,韧性等级为Ⅰ级:S(漯河市)= 1.97,S(三门峡市) = 1. 89,S (鹤壁市) = 1. 85,S(濮阳市)= 1.63,S(济源市)= 1.56。
3.4 计算结果分析
从整体上看,河南省洪涝灾害韧性等级为Ⅲ级(中等韧性),但各地级市洪涝灾害韧性等级差异较大,在空间上呈现西南部城市韧性等级高于东北部城市韧性等级的特点(见图2)。仅郑州市洪涝灾害韧性等级为Ⅳ级,其余城市洪涝灾害韧性等级为Ⅲ级、Ⅱ级甚至I 级,主要原因是郑州市相对于其他城市在第三产业产值占比、政府救助系统、救援响应时间等方面具有一定的优越性,因此在洪涝灾害面前具有较强的韧性。上述分析结果与河南省防汛抗旱指挥部印发的《河南省防汛应急预案》(2022 年修订版)中的各地级市洪涝灾害韧性等级基本一致。
为进一步验证本文方法评价结果的合理性和可靠性,对已有河南省洪涝灾害韧性研究成果进行分析,其中蒋硕亮等[8] 对河南省18 个地级市洪涝灾害韧性进行评估,认为河南省整体洪涝灾害韧性等级为中等,各城市洪涝灾害韧性等级差异较大,西南部城市韧性等级高于东北部城市的。冯凌彤[23] 利用GIS 技术评价河南省洪涝灾害,认为漯河市、平顶山市、许昌市城市洪涝灾害风险较高,抗风险能力较差,洪涝灾害韧性等级较低。综上所述,本文评价结果与相关研究成果基本一致,具有一定的合理性和可靠性。
在主要影响指标维度上,结合有序加权平均算子得出的指标权重可知抵抗能力和适应能力是影响河南省洪涝灾害韧性的主要因素。其中,抵抗能力中主要影响因素是第三产业产值占比、汛期单日最大降雨量和常住人口密度,预警能力中主要影响因素是水文站覆盖率、互联网用户覆盖率和移动电话用户覆盖率,适应能力中主要影响因素是主城区道路易积水点数量、排水管网密度和救援响应时间,恢复能力中主要影响因素是洪涝灾害专项资金、登记失业率、政府救助系统。
3.5 不同方法对比分析
为检验本文提出的城市洪涝灾害韧性评价方法的合理性,在相同韧性等级下,采用经典的模糊层次分析法、层次分析法和后悔理论、有序加权平均算子和灰色理论3 种方法评价河南省洪涝灾害韧性,与本文评价结果比较见表4。
由表4 可知,4 种评价方法均认为河南省洪涝灾害韧性等级为Ⅲ级,结论基本一致,表明本文所建立的洪涝灾害韧性评价方法具有可行性和可靠性,可为城市洪涝灾害韧性评价提供新思路。此外,方法4 较方法2 评价值增大0.13,主要原因是有序加权平均算子较层次分析法能更好地处理专家极端决策值带来的权重偏差。方法4 较方法3 评价值增大0.25,表明后悔理论能有效规避专家决策时的后悔心理,使之作出更加客观的评判。综上,本文提出的基于有序加权平均算子和后悔理论的城市洪涝灾害韧性评价方法计算结果合理、可靠,适用于城市洪涝灾害韧性评价。
4 结论
1)河南省洪涝灾害韧性等级为Ⅲ级(中等韧性),从空间分布看,各地级市洪涝灾害韧性呈现西南部城市韧性等级高于东北部城市韧性等级的特征,仅郑州市洪涝灾害韧性等级为Ⅳ级、达到较高韧性水平,其余城市洪涝灾害韧性等级为Ⅲ级、Ⅱ级甚至Ⅰ级。
2)将本文评价结果与模糊层次分析法、层次分析法和后悔理论、有序加权平均算子和灰色理论3 种方法评价结果进行对比,表明4 种评价方法的结果基本一致。相较于其他评价方法,本方法评价值更为准确、灵活性更强,适用于城市洪涝灾害韧性评价。
3)将有序加权平均算子和后悔理论评价方法运用于城市洪涝灾害韧性评价,并指出抵抗能力和适应能力构成洪涝灾害韧性主要影响指标。该方法不仅能降低专家决策数据主观性过强带来的权重失真,而且进一步考虑专家决策过程的心理特点,充分保留专家决策信息完整性,提升了评价结果的合理性和可靠性。
参考文献:
[1] 黄国如,杨格,曾博威,等.基于绿灰蓝基础设施融合的城市洪涝灾害调控[J].郑州大学学报(工学版),2023,44(2):14-21,74.
[2] 殷杰,许世远,经雅梦,等.基于洪涝情景模拟的城市公共服务灾害应急响应空间可达性评价:以医疗急救为例[J].地理学报,2018,73(9):1737-1747.
[3] 李正兆,傅大放,王君娴,等.应对内涝灾害的城市韧性评估模型及应用[J].清华大学学报(自然科学版),2022,62(2):266-276.
[4] MIGUEZ M G,VER L,ALINE P.A Catchment Scale Integrated Flood Resilience Index to Support Decision Making in Urban Flood Control Design[J].Environment and Planning B:Urban Analytics and City Science,2017,44(5):925-946.
[5] BERTILSSON L,WIKLUND K,MOUDA T I,et al. Urban Flood Resilience:A Multi⁃Criteria Index to Integrate Flood Resilience into Urban Planning[J]. Journal of Hydrology,2019,573:970-982.
[6] 黄晶,佘靖雯,袁晓梅,等.基于系统动力学的城市洪涝韧性仿真研究:以南京市为例[J].长江流域资源与环境,2020,29(11):2519-2529.
[7] 周铭毅,尚志海,蔡灼芬,等.基于VIKOR 方法的广东省城市洪涝灾害韧性评估[J].灾害学,2023,38(1):206-212.
[8] 蒋硕亮,陈贤胜.洪涝情景下城市韧性评估及障碍因素分析[J].统计与决策,2022,38(24):63-67.
[9] 孟晓静,陈鑫,陈佳静,等.组合赋权-TOPSIS 在洪涝灾害下城市区域韧性评估中的应用[J].安全与环境学报,2023,23(5):1465-1473.
[10] 李志荣,朱金善,宋健,等.基于C-OWA 算子和证据理论的船舶夜航光环境安全评价[J].中国安全生产科学技术,2019,15(12):182-187.
[11] 董兴芝,李平,栾中,等.基于OWA-差异驱动组合赋权的高速铁路Wi-Fi 车地通信评价研究[J].铁道运输与经济,2022,44(3):127-133.
[12] 刘同超,李振福,陈卓.基于后悔理论的北极航线安全评价[J].安全与环境学报,2018,18(6):2069-2074.
[13] 郭隆鑫,李希建,刘柱.基于后悔理论-集对势分析的矿山紧急避险系统安全评价[J].矿业研究与开发,2020,40(4):85-90.
[14] 罗景峰,李登峰.考虑后悔行为的变权可变模糊旅游安全预警方法[J]. 中国安全科学学报,2018,28 (1):161-166.
[15] 张鹏,范潮海,王璟.基于前景-后悔理论的油气管道风险模态分析[J].安全与环境学报,2021,21(6):2365-2370.
[16] 孙少楠,潘传旭,张卫东.基于C-OWA 算子的水利设计企业BIM 应用能力评价[J].水利水电技术(中英文),2021,52(9):118-125.
[17] 李华,张旭旭,高红,等.智慧工地本质安全度评价方法研究[J].中国安全生产科学技术,2022,18(1):139-145.
[18] 周英豪,王文杰,卢西洲,等.岩爆灾害博弈论组合赋权预测模型及应用[J].中国安全科学学报,2022,32(7):105-112.
[19] 李梦杰,刘德林.河南省洪涝灾害的灾后恢复力研究[J].水土保持通报,2020,40(6):200-204.
[20] 王欲敏,郑梦琴,马文光,等.基于改进SPA-IAHP 赋权与集对分析法的高海拔公路施工成本风险评价[J].工程管理学报,2021,35(3):129-134.
[21] 马晓雪,刘阳,张靖雯,等.基于模糊AHP 的北极航线安全保障体系结构分析[J].中国航海,2020,43(2):1-6,13.
[22] 艾春安,奉非东,李剑,等.基于区间数判断矩阵与模糊聚类分析的主观群决策AHP 方法[J].统计与决策,2019,35(2):39-43.
[23] 冯凌彤.基于GIS 技术的河南省洪灾风险评估与分析[J].人民黄河,2022,44(5):67-70,74.
【责任编辑 许立新】
基金项目:国家重点研发计划项目(2017YFC0403505)