«上一篇 下一篇»
  计算机工程  2019, Vol. 45 Issue (10): 314-320  DOI: 10.19678/j.issn.1000-3428.0052813
0

引用本文  

许硕, 唐作其, 王鑫. 基于D-AHP与TOPSIS的突发事件应急管理能力评估[J]. 计算机工程, 2019, 45(10), 314-320. DOI: 10.19678/j.issn.1000-3428.0052813.
XU Shuo, TANG Zuoqi, WANG Xin. Emergency Management Capability Assessment Based on D-AHP and TOPSIS[J]. Computer Engineering, 2019, 45(10), 314-320. DOI: 10.19678/j.issn.1000-3428.0052813.

基金项目

贵州省科技支撑计划(黔科合人才平台[2018]5616);贵州大学青年教师科研基金(贵大青合字(2013)01号)

通信作者

唐作其(通信作者), 副教授

作者简介

许硕(1993-), 男, 硕士研究生, 主研方向为应急管理系统;
王鑫, 硕士研究生

文章历史

收稿日期:2018-10-08
修回日期:2018-11-09
基于D-AHP与TOPSIS的突发事件应急管理能力评估
许硕 , 唐作其 , 王鑫     
贵州大学 计算机科学与技术学院, 贵阳 550025
摘要:在突发事件应急管理评估过程中,参评专家经验及语言描述差异会导致评估信息存在不确定性和模糊性,直接影响评估结果。针对该问题,基于D数偏好关系改进层次分析法(D-AHP)和逼近理想解排序法(TOPSIS),提出一种突发事件应急管理能力评估方法。根据应急管理实际情况建立评估指标体系,通过构建D-AHP层次结构模型求解各评估指标的影响权重,同时结合专家的评估结果,利用TOPSIS方法对评估对象的突发事件应急管理能力水平进行排序。对某市下辖4个区县的评估结果表明,该方法能有效识别应急管理过程中的关键环节,直观展示各个评估对象的实际水平并有助于给出完善应急管控建设的针对性建议。
关键词突发事件应急管理    多准则决策    D数理论    逼近理想解排序法    层次分析法    模糊偏好关系    
Emergency Management Capability Assessment Based on D-AHP and TOPSIS
XU Shuo , TANG Zuoqi , WANG Xin     
College of Computer Science and Technology, Guizhou University, Guiyang 550025, China
Abstract: In the emergency management assessment process, the differences in the experience and language description of the experts can cause uncertainty and ambiguity in the assessment information, directly affecting the assessment results.Aiming at this problem, an assessment method for emergency management capability is proposed based on Analytic Hierarchy Process extended by D numbers preference relation (D-AHP) and Technique for Order Preference by Similarity to Ideal Solution(TOPSIS).According to the actual situation of emergency management, an assessment index system is established.By constructing the D-AHP hierarchical structure model, the impact weights of each assessment index is solved.Combined with the assessment results of the experts, the TOPSIS method is used to rank the emergency management capability level of the assessment objects.The assessment results of four districts and counties of a certain city show that the method can effectively identify the key links in the emergency management process, visually display the actual level of each assessment object and help to give specific suggestions for improving the emergency management and control construction.
Key words: emergency management    Multi-Criterion Decision Making(MCDM)    D number theory    Technique for Order Preference by Similarity to Ideal Solution(TOPSIS)    Analytic Hierarchy Process(AHP)    fuzzy preference relation    
0 概述

在我国经济社会快速发展的同时, 各类突发事件也呈频发态势, 其危害日益凸显, 对社会稳定及人民群众的生命财产安全造成极大威胁。对突发事件应急管理能力进行评估, 是减少突发事件附带损失、完善应急管控措施的重要手段, 一直是国内外研究的热点。

近年来, 一些学者提出各种模型及方法对突发事件应急管理能力进行评估。文献[1]考虑复杂环境中评估者的犹豫性和模糊性对评估结果的影响, 提出基于犹豫模糊Einstein算子的突发事件风险管理能力评估模型。文献[2]提出针对应急管理能力评估的模糊综合评价法, 在一定程度上消除了模糊性, 提高了评估精度。文献[3]提出基于数据包络分析(Data Envelopment Analysis, DEA)和逼近理想解排序方法(Technique for Order Preference by Similarity to Ideal Solution, TOPSIS)的评估方法, 通过建立数据包络函数求解权重, 并结合TOPSIS进行综合排序。文献[4]提出基于灰色层次分析法(Gray Analytic Hierarchy Process, G-AHP)的突发事件应急管理能力评价方法, 建立层次化的评估体系并量化显示评估结果。文献[5]将D数理论与决策实验室法(Decision-Making Trial and Evaluation Laboratory, DEMATEL)相结合, 通过融合专家评估信息, 分析影响应急管理各种因素之间的关联关系, 识别应急管理过程中的关键因素。

对突发事件应急管理能力进行评估, 确定评估对象的应急能力水平, 涉及到大量评估准则和影响因素, 是典型的多准则决策(Multi-Criterion Decision Making, MCDM)问题。TOPSIS是MCDM问题的代表方法, 被用于求解距离正理想解最近、负理想解最远的方案。虽然该方法具有一定的客观性和逼近性, 但当评估指标体系过于复杂时, 较难求解指标权重。而AHP常用于MCDM问题中的指标权重确定, 该方法通过建立层次化结构模型, 系统分析评估对象内部各指标之间的相互关系, 求得的指标权重更加合理客观。但由于评估专家经验差异, 因此产生的不确定性信息无法构建合理成对比较矩阵。通过D数偏好关系对AHP进行改进, 在利用不确定性评估信息求解权重时具有明显优势。目前该方法已应用于桥梁状态评估[6]、帷幕灌浆效率评估[7]、投资决策[8]等方面, 且取得了较好的效果。

结合突发事件应急管理的实际情况, 本文提出一种基于D-AHP和TOPSIS的评估方法。使用D-AHP方法分配指标权重, 识别出应急管理过程中的关键环节, 为应急管控建设提出针对性建议。在此基础上, 结合TOPSIS方法综合考虑应急管理过程中的多种影响因素, 对各评估对象按照应急管理能力水平优劣进行排序, 为应急管理评估提供新的思路。

1 理论分析 1.1 D数理论

突发事件应急管理能力进行评估较为依赖专家评估经验, 而专家评估经验的差异会造成评估信息的不确定性, 同时语言性评估信息存在着一定的模糊性。在对不确定性及模糊性信息的处理上, 常用方法有D-S证据理论、模糊集理论、灰色系统理论等。其中D-S证据理论优势如下[9]:

1) 该理论可处理随机性导致的不确定性及模糊性导致的不确定性, 作用对象相对较广。

2) 通过融合规则, 该理论可将已有证据或信任函数合成为新的证据或信任函数。

但证据理论的适用范围有限, 如要求满足识别框架中命题之间的互斥性假设和基本概率分配(Basic Probability Assignment, BPA)的完整性约束[10]。针对上述问题, 文献[11]提出了D数理论。

D数理论继承了证据理论在处理不确定性信息方面的优势, 尤其可以更好地处理语言评估所产生的不确定性, 扩大了其适用范围。关于D数的定义和性质如下[11]:

定义1  设存在有限非空集Ω, 存在映射D:

$ D:\mathit{\Omega } \to [0,1] $ (1)

满足:

$ \sum\limits_{B \subseteq \mathit{\Omega }} D (B) \le 1,D(\theta ) = 0 $ (2)

则称映射D为D数, 其中BΩ的一个子集, θ为空集。若$\sum\limits_{B \subseteq \mathit{\Omega }} D (B) = 1$, 则说明由D数所表示的信息是完整的; 否则, 说明信息是不完整的。

定义2  设存在D数D和非空有限集Ω, 则D的信息完整度Q可量化表示为:

$ Q = \sum\limits_{B \subseteq \mathit{\Omega }} D (B) $ (3)

定义3  设离散集Ω={b1, b2, …, bn}, D数的特殊表达形式为:

$ \begin{array}{*{20}{c}} {D\left( {\left\{ {{b_1}} \right\}} \right) = {v_1}}\\ {D\left( {\left\{ {{b_2}} \right\}} \right) = {v_2}}\\ \vdots \\ {D\left( {\left\{ {{b_n}} \right\}} \right) = {v_n}} \end{array} $ (4)

或简单表示为:

$ D = \left\{ {\left( {{b_1},{v_1}} \right),\left( {{b_2},{v_2}} \right), \cdots ,\left( {{b_n},{v_n}} \right)} \right\} $ (5)

其中, vi≥0且$\sum\limits_{i = 1}^n {{v_i}} \le 1$

该表示形式可方便高效地描述不确定性问题。

定义4  设存在D数:

D={(b1, v1), (b2, v2), …, (bn, vn)}, 则其集成可表示为:

$ I(D) = \sum\limits_{i = 1}^n {{b_i}} {v_i} $ (6)

可将D数的集成称为I值。

1.2 D数偏好关系

在对评估样本进行比较赋值, 构造决策矩阵时, 较常用的方法为模糊偏好关系。模糊偏好关系可以较好地反应评估专家对具体方案相对重要性的重要程度, 使用“$ \succ $”表示[12]

定义5  设存在评估样本A集为{A1, A2, …, An}, 其模糊偏好关系为:

$ {\mu _R}:A \times A \to [0,1] $ (7)

以矩阵形式表示为R=[rij]n×n:

$ \begin{array}{*{20}{c}} {\;\;\;\;\;\;\;\;{A_1}\quad {A_2}\quad \cdots \quad {A_n}}\\ {\mathit{\boldsymbol{R}} = \begin{array}{*{20}{c}} {{A_1}}\\ {{A_2}}\\ \vdots \\ {{A_n}} \end{array}\left[ {\begin{array}{*{20}{c}} {{r_{11}}}&{{r_{12}}}& \cdots &{{r_{1n}}}\\ {{r_{21}}}&{{r_{22}}}& \cdots &{{r_{2n}}}\\ \vdots & \vdots &{}& \vdots \\ {{r_{n1}}}&{{r_{n2}}}& \cdots &{{r_{nn}}} \end{array}} \right]} \end{array} $ (8)

其中:

1) rij≥0。

2) rij+rji=1, $\forall $i, j∈{1, 2, …, n}。

3) rii=0.5, $\forall $i∈{1, 2, …, n}。

矩阵元素rij表示评估专家认为Ai相对于Aj的重要程度, 其取值范围及含义如下:

$ {r_{ij}} = {\mu _R}\left( {{A_i},{A_j}} \right) $

1) 如果rij=0, 则Aj明显比Ai重要。

2) 如果0 < rij < 0.5, 则Aj在某种程度上比Ai重要。

3) 如果rij=0.5, 则AiAj无差别。

4) 如果0.5 < rij < 1, 则Ai在某种程度上比Aj重要。

5) 如果rij=1, 则Ai明显比Aj重要。

但该方法存在一定的使用限制, 在评估信息不确定的情况下, 无法构建出合适的偏好关系矩阵。

针对模糊偏好关系这一局限性, 文献[12]使用D数理论对模糊偏好关系进行改进, 将模糊偏好关系的适用范围扩大到不确定信息领域。对应的评判矩阵被称为D数偏好矩阵, 简称D矩阵。

定义6  设一组评估样本A ={A1, A2, …, An}, 其D数偏好关系为:

$ {\mathit{\boldsymbol{R}}_\mathit{\boldsymbol{D}}}:A \times A \to D $ (9)

以矩阵形式表示为RD=[Dij]n×n:

$ \begin{array}{*{20}{c}} {\;\;\;\;\;\;\;\;\;\;\;\;{A_1}\quad {A_2}\quad \cdots \quad {A_n}}\\ {{\mathit{\boldsymbol{R}}_\mathit{\boldsymbol{D}}} = \begin{array}{*{20}{c}} {{A_1}}\\ {{A_2}}\\ \vdots \\ {{A_n}} \end{array}\left[ {\begin{array}{*{20}{c}} {{D_{11}}}&{{D_{12}}}& \cdots &{{D_{1n}}}\\ {{D_{21}}}&{{D_{22}}}& \cdots &{{D_{2n}}}\\ \vdots & \vdots &{}& \vdots \\ {{D_{n1}}}&{{D_{n2}}}& \cdots &{{D_{nn}}} \end{array}} \right]} \end{array} $ (10)

其中:

1) Dij={(bij1, vij1), (bij2, vij2), …, (bmij, vmij)}, Dji={(1-bij1, vij1), (1-bij2, vij2), …, (1-bmij, vmij)}, $\forall $i, j∈{1, 2, …, n}。

2) bkij∈[0, 1], $\forall $k∈{1, 2, …, m}。

3) Dij={(0.5, 1.0)}, $\forall $i∈{1, 2, …, n}。

bkij表示第k位专家认为第i个方案相对于第j个方案的重要程度, vkij表示该专家对该重要度的支持程度。

由D数性质可知, 通过D数偏好关系可以较好地消除评估专家经验背景的差异及评估的模糊性对评估结果的影响。

文献[13]给出了D数偏好关系计算评估对象优先权重的过程, 可分为以下步骤:

步骤1  组织参评专家, 根据评估标度对评估对象进行成对比较, 以偏好关系的形式表示各评估对象的重要程度, 并构建D矩阵RD

步骤2  由式(6), 将矩阵RD集成表示为确定数矩阵RC

步骤3  计算评估对象间的偏好概率, 构建基于确定数矩阵RC的概率矩阵RP

步骤4  对概率矩阵RP进行三角化处理, 得到三角概率矩阵RPT, 并按照重要程度对评估对象进行排序。

步骤5  根据得到的指标排序, 将确定数矩阵RC表示为矩阵RCT, 并计算评估对象权重。

在计算权重时, 由式(11)计算D矩阵RD的不一致度系数。

$ \begin{array}{*{20}{l}} {\sum\limits_{i = 1}^n {\mathit{\boldsymbol{R}}_P^{\rm{T}}} (i,j)}\\ {\mathit{I}{\rm{. }}D. = \frac{{j < i}}{{n(n - 1)/2}}} \end{array} $ (11)

其中, RPT(i, j)表示三角化概率矩阵RPT中的元素, n表示成对比较的评估对象的个数。

1.3 D-AHP方法

AHP是一种灵活实用的MCDM方法[14-15], 将问题分解成多个因素, 通过支配关系形成层次结构, 再以两两比较的方式确定待决策方案的相对重要性。

但不完整的评估信息限制了典型AHP方法的使用。而通过D数偏好关系扩展AHP得到的D-AHP方法, 由D数偏好关系求解指标权重, 结合AHP层次结构对各级指标权重进行逐层集成, 可以较好地应用于信息不确定环境下的复杂评估决策问题, 有效降低评估过程中的主观性和不确定性对评估结果的影响。D-AHP层次结构如图 1所示。

Download:
图 1 D-AHP层次结构

在D-AHP层次结构中, 目标层包括评估目标、待决策问题等, 准则层包括各级评估指标及属性, 而方案层包括各个具有上述指标及属性的待评估对象。

D-AHP通过对各层次元素权重进行集成求解权重, 其过程为:设准则层指标Cj(j=1, 2, …, m)相对于评估目标的权重为wcj, wcj≥0且$\sum\limits_{j = 1}^m {{w_{{c_j}}}} = 1$; 设方案层评估对象Ai(i=1, 2, …, m)相对于准测层指标Cj的权重为wAij, wAij≥0且$\sum\limits_{i = 1}^n {{w_{{A_{ij}}}}} = 1$。由式(12)求解评估对象Ai相对于评估目标的综合权重:

$ \begin{array}{*{20}{l}} {{w_i} = \sum\limits_{j = 1}^m {{w_{{A_{ij}}}}} {w_{{c_j}}},i = 1,2, \cdots ,n}\\ {{w_i} \ge 0\;且\;\sum\limits_{i = 1}^n {{w_i}} = 1}。\end{array} $ (12)
1.4 TOPSIS方法

TOPSIS[16]是求解多目标决策问题的常用方法, 其求解思路为:首先假设存在一个正理想解和一个负理想解, 通过计算评估样本与正负理想解的距离及贴进度, 作为各个评估样本的优劣程度的标准, 距离正理想解最近的样本为最优方案, 距离负理想解最近的样本为最劣方案。该方法的计算步骤如下:

步骤1  构建初始评判矩阵。针对研究问题, 设存在评估样本集A={A1, A2, …, An}, 每个样本的指标集C={C1, C2, …, Cm}。构建初始评判矩阵MA=[aij]n×m:

$ {\mathit{\boldsymbol{M}}_A} = \left[ {\begin{array}{*{20}{c}} {{a_{11}}}&{{a_{12}}}& \cdots &{{a_{1m}}}\\ {{a_{21}}}&{{a_{22}}}& \cdots &{{a_{2m}}}\\ \vdots & \vdots &{}& \vdots \\ {{a_{n1}}}&{{a_{n2}}}& \cdots &{{a_{nm}}} \end{array}} \right] $ (13)

其中, 矩阵元素aij表示第i个评估样本的第j个指标的量化评估值。

步骤2  对初始评判矩阵进行规范化处理。评估指标分为极大性指标和极小性指标。对于极大性指标, 其值越大越好; 对于极小性指标, 其值越小越好。根据评估指标的不同极性, 由式(14)和式(15)对上述2种评估指标进行规范化处理。

对于极大性指标:

$ {b_{ij}} = \frac{{{a_{ij}} - \min \left( {{a_{ij}}} \right)}}{{\max \left( {{a_{ij}}} \right) - \min \left( {{a_{ij}}} \right)}} $ (14)

对于极小性指标:

$ {b_{ij}} = \frac{{\max \left( {{a_{ij}}} \right) - {a_{ij}}}}{{\max \left( {{a_{ij}}} \right) - \min \left( {{a_{ij}}} \right)}} $ (15)

bij构成规范化评判矩阵MB=[bij]n×m

步骤3  构建加权评判矩阵。评估指标权重与规范化判断矩阵MB对应相乘, 可得加权矩阵MC=[cij]n×m, cij=wjbij

$ {\mathit{\boldsymbol{M}}_C} = \left[ {\begin{array}{*{20}{c}} {{w_1}{b_{11}}}&{{w_2}{a_{12}}}& \cdots &{{w_m}{a_{1m}}}\\ {{w_1}{b_{21}}}&{{w_2}{b_{22}}}& \cdots &{{w_m}{a_{2m}}}\\ \vdots & \vdots &{}& \vdots \\ {{w_1}{b_{n1}}}&{{w_2}{b_{n2}}}& \cdots &{{w_m}{b_{nm}}} \end{array}} \right] $ (16)

其中, wj表示第j个指标的权重。

步骤4  计算正负理想解。由式(17)和式(18)计各评估样本指标的理想解。

$ S_j^ + = \mathop {\max }\limits_{1 \le i \le n} \left\{ {{c_{ij}}} \right\},j = 1,2, \cdots ,m $ (17)
$ S_j^ - = \mathop {\min }\limits_{1 \le i \le n} \left\{ {{c_{ij}}} \right\},j = 1,2, \cdots ,m $ (18)

其中, Sj+表示正理想解, Sj表示负理想解。

步骤5  计算评估样本和理想解之间的欧式距离。由式(19)和式(20)计算样本与理想解之间的欧式距离:

$ D_i^ + = \sqrt {\sum\limits_{j = 1}^m {{{\left( {{c_{ij}} - S_j^ + } \right)}^2}} } ,i = 1,2, \cdots ,n $ (19)
$ D_i^ - = \sqrt {\sum\limits_{j = 1}^m {{{\left( {{c_{ij}} - S_j^ - } \right)}^2}} } ,i = 1,2, \cdots ,n $ (20)

其中, Di+Di分别表示评估样本与正、负理想解之间的欧式距离。

步骤6  贴进度分析。由式(21)计算各样本的贴进度, 并根据贴进度大小对各样本进行相对优劣排序。

$ {D_i} = \frac{{D_i^ - }}{{D_i^ + + D_i^ - }},i = 1,2, \cdots ,n $ (21)
2 评估指标体系的构建

突发事件应急管理能力评估过程复杂, 涉及到大量影响因素。为确保评估结果的客观准确, 必须建立科学有效的评估指标体系。在总结前人工作[17-18]并通过专家访谈的基础上, 本文结合应急管理的实际情况, 构建突发事件应急管理能力评估指标体系, 确定应急预防与准备、监测与预警、应急处置与救援和事后恢复与重建4个一级指标以及预案编制、演练周期、保障力度、信息监测等14个二级指标。构建的突发事件应急管理评估指标体系如图 2所示。

Download:
图 2 突发事件应急管理能力评估指标体系
3 算例分析

本文以某市下辖的A县、B县、C县和D区为研究对象, 根据建立的突发事件应急管理能力评估指标体系, 在实地调研的基础上收集相关数据, 进行如下评估工作。

3.1 D-AHP层次结构模型 3.1.1 D-AHP层次结构的构建

结合突发事件应急管理能力评估指标体系, 根据突发事件应急管理实际情况, 构建基于D-AHP的突发事件应急管理能力评估层次结构模型。

确定准则层一级指标集C={C1, C2, C3, C4}={应急预防与准备, 监测与预警, 应急处置与救援, 事后恢复与重建}, 并确定一、二级评估指标间隶属关系:C1={P1, P2, P3}={预案编制, 演练周期, 保障力度}, C2={P4, P5, P6, P7}={信息监测, 预警研判, 报警流程, 预警信息收发}, C3={P8, P9, P10}={险情报送, 响应时间, 救援能力}, C4={P11, P12, P13, P14}={经济补偿, 重建计划, 损失评估, 总结报告}。

在该层次结构模型中, 确定突发事件应急管理能力为目标层元素; 由一、二级评估指标组成D-AHP层次结构的准则层; A县、B县、C县和D区4个评估对象组成方案层。构建的突发事件应急管理能力评估D-AHP层次结构模型如图 3所示。

Download:
图 3 突发事件应急管理能力评估D-AHP层次结构模型
3.1.2 权重计算

为定量表示指标之间的相对重要性, 采用0.1~0.9标度进行成对比较赋值。标度及对应含义如表 1所示。

下载CSV 表 1 评估标度及含义

参考上述评估标准, 首先计算准测层中各一级指标相对于突发事件应急管理能力的权重:构建基于D数偏好关系的D矩阵RD; 由D数的集成表示, 将矩阵RD转化为确定数矩阵RC; 构建基于RC的概率矩阵RP; 对概率矩阵RP进行三角化处理, 得到三角概率矩阵RPT

$ {\mathit{\boldsymbol{R}}_D} = \begin{array}{*{20}{c}} {}&{\begin{array}{*{20}{c}} {\;\;\;\;\;\;\;\;\;\;{C_1}\;\;}&{\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{C_2}\;\;\;\;\;}&{\;\;\;\;\;\;\;\;\;\;{C_3}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;}&{\;\;\;\;{C_4}\;\;\;\;\;\;\;\;\;\;} \end{array}\;}\\ {\begin{array}{*{20}{c}} {{C_1}}\\ {{C_2}}\\ {{C_3}}\\ {{C_4}} \end{array}}&{\left[ {\begin{array}{*{20}{c}} {\{ (0.50,1.00)\} }&{\{ (0.20,1.00)\} }&{\{ (0.70,1.00)\} }&{\left\{ {(0.70,0.60),(0.60,0.40)} \right\}}\\ {\{ (0.80,1.00)\} }&{\{ (0.50,1.00)\} }&{\{ (0.60,1.00)\} }&{\{ (0.60,0.80)\} }\\ {\{ (0.30,1.00)\} }&{\{ (0.40,1.00)\} }&{\{ (0.50,1.00)\} }&{\{ (0.80,1.00)\} }\\ {\left\{ {(0.30,0.60),(0.40,0.40)} \right\}}&{\{ (0.40,0.80)\} }&{\{ (0.20,1.00)\} }&{\{ (0.50,1.00)\} } \end{array}} \right]} \end{array} $
$ {\mathit{\boldsymbol{R}}_C} = I\left( {{\mathit{\boldsymbol{R}}_D}} \right) = \begin{array}{*{20}{c}} {}&{\begin{array}{*{20}{c}} {{C_1}\;\;\;}&{{C_2}\;\;}&{\;{C_3}}&{\;\;\;{C_4}} \end{array}}\\ {\begin{array}{*{20}{c}} {{C_1}}\\ {{C_2}}\\ {{C_3}}\\ {{C_4}} \end{array}}&{\left[ {\begin{array}{*{20}{c}} {0.50}&{0.20}&{0.70}&{0.66}\\ {0.80}&{0.50}&{0.60}&{0.48}\\ {0.30}&{0.40}&{0.50}&{0.80}\\ {0.34}&{0.32}&{0.20}&{0.50} \end{array}} \right]} \end{array} $
$ {\mathit{\boldsymbol{R}}_P} = \begin{array}{*{20}{c}} {}&{\begin{array}{*{20}{c}} {{C_1}}&{{C_2}}&{{C_3}}&{{C_4}} \end{array}}\\ {\begin{array}{*{20}{c}} {{C_1}}\\ {{C_2}}\\ {{C_3}}\\ {{C_4}} \end{array}}&{\left[ {\begin{array}{*{20}{c}} 0&0&1&{0.9}\\ 1&0&1&0\\ 0&0&0&1\\ 0&{0.1}&0&0 \end{array}} \right]} \end{array} $
$ \mathit{\boldsymbol{R}}_P^{\rm{T}} = \begin{array}{*{20}{c}} {}&{\begin{array}{*{20}{c}} {{C_2}}&{{C_1}}&{{C_3}}&{{C_4}} \end{array}}\\ {\begin{array}{*{20}{c}} {{C_2}}\\ {{C_1}}\\ {{C_3}}\\ {{C_4}} \end{array}}&{\left[ {\begin{array}{*{20}{c}} 0&1&1&0\\ 0&0&1&{0.9}\\ 0&0&0&1\\ {0.1}&0&0&0 \end{array}} \right]} \end{array} $

由式(11)计算D矩阵的不一致系数, 得I.D.=0.02。经评估专家组讨论, 认为该系数在可接受范围之内。

对各指标进行排序:C2$ \succ $C1$ \succ $C3$ \succ $C4。即对于突发事件应急管理能力评估结果影响程度由高到低依次为监测与预警(C2)、应急预防与准备(C1)、应急处置与救援(C3)和事后恢复与重建(C4)。

由准则层一级指标排序结果, 将矩阵RC表示为RCT:

$ \mathit{\boldsymbol{R}}_C^{\rm{T}} = \begin{array}{*{20}{c}} {}&{\begin{array}{*{20}{c}} {{C_2}\;\;\;}&{{C_1}\;\;}&{{C_3}\;}&{\;{C_4}} \end{array}}\\ {\begin{array}{*{20}{c}} {{C_2}}\\ {{C_1}}\\ {{C_3}}\\ {{C_4}} \end{array}}&{\left[ {\begin{array}{*{20}{c}} {0.50}&{0.80}&{0.60}&{0.58}\\ {0.20}&{0.50}&{0.70}&{0.66}\\ {0.40}&{0.30}&{0.50}&{0.80}\\ {0.42}&{0.34}&{0.20}&{0.50} \end{array}} \right]} \end{array} $

根据该矩阵解方程组:

$ \left\{ {\begin{array}{*{20}{l}} {\lambda \left( {{w_{{c_2}}} - {w_{{c_1}}}} \right) = 0.8 - 0.5}\\ {\lambda \left( {{w_{{c_1}}} - {w_{{c_3}}}} \right) = 0.7 - 0.5}\\ {\lambda \left( {{w_{{c_3}}} - {w_{{c_4}}}} \right) = 0.8 - 0.5}\\ {{w_{{c_1}}} + {w_{{c_2}}} + {w_{{c_3}}} + {w_{{c_4}}} = 1} \end{array}} \right. $

解得:

$ \left\{ {\begin{array}{*{20}{l}} {{w_{{c_1}}} = 1/4 + 0.1/\lambda }\\ {{w_{{c_2}}} = 1/4 + 0.4/\lambda }\\ {{w_{{c_3}}} = 1/4 - 0.1/\lambda }\\ {{w_{{c_4}}} = 1/4 - 0.4/\lambda } \end{array}} \right. $

λ∈[1.6, +∞), 可得wci的取值范围:

$ \left\{ {\begin{array}{*{20}{l}} {{w_{{c_1}}} \in [0.2500,0.3125]}\\ {{w_{{c_2}}} \in [0.2500,0.5000]}\\ {{w_{{c_3}}} \in [0.1875,0.2500]}\\ {{w_{{c_4}}} \in [0.0000,0.2500]} \end{array}} \right. $

其中, wci表示第i个准测层一级指标权重, λ表示评估信息的可信度, 其取值与专家经验和对问题背景的认知水平有关。文献[12]给出了关于λ的取值及说明:

$ \lambda = \left\{ \begin{array}{l} \left\lceil {\underline \lambda } \right\rceil ,信息高度可信\\ n,信息中度可信\\ \frac{{{n^2}}}{2},信息低度可信 \end{array} \right. $

由于参评专家评估经验丰富, 评估信息可信度较高, 因此λ=$\lceil$λ$\rceil$=2。解得准则层各一级指标权重:wc1=0.300 0, wc2=0.450 0, wc3=0.200 0, wc4=0.050 0, 及一级指标权重向量WC=(0.300 0, 0.450 0, 0.200 0, 0.050 0)。

同理, 可求得各二级指标相对于准则层指标的权重wpij, 以及相对于评估目标的综合权重wj, 如表 2所示。

下载CSV 表 2 指标权重计算结果

求得二级指标综合权重向量W=(w1, w2, …, w14)=(0.190 0, 0.070 0, 0.040 0, 0.104 4, 0.101 6, 0.084 4, 0.159 7, 0.072 5, 0.054 7, 0.072 8, 0.015 4, 0.007 4, 0.015 7, 0.011 5)。

3.2 TOPSIS综合排序

结合确定的评估指标体系, 组织10名评估专家对方案层4个评估对象的各项指标进行按照1分~10分的评估标准进行赋值。对评估指标体系进行分析, 经过专家组讨论, 确定极大性指标有10个, 包括预案编制(P1)、保障力度(P3)、信息监测(P4)、预警研判(P5)、预警信息收发(P7)、救援能力(P10)、经济补偿(P11)、重建计划(P12)、损失评估(P13)和总结报告(P14), 其评估赋值越大越好; 极小性指标有4个, 包括演练周期(P2)、报警流程(P6)、险情报送(P8)和响应时间(P9), 其评估赋值越小越好。

将专家赋值取算数平均作为最终结果, 构建初始评判矩阵, 如表 3所示。

下载CSV 表 3 初始评判矩阵

根据式(14)和式(15), 对初始评判矩阵中不同极性的指标赋值进行标准化处理, 得到标准化评判矩阵。结合D-AHP求得的二级指标权重向量W对其进行加权处理, 构建的加权评判矩阵, 如表 4所示。

下载CSV 表 4 加权评判矩阵

取各项指标的最大值和最小值作为正、负理想解, 由式(19)~式(21)求得各评估对象与理想解的欧式距离及贴进度, 并进行排序, 如表 5所示。

下载CSV 表 5 各评估对象的欧式距离及贴进度
3.3 结果分析

表 2可知, 各一级指标对突发事件应急管理能力的影响由高到低依次为监测与预警(C2)、应急预防与准备(C1)、应急处置与救援(C3)和事后恢复与重建(C4), 因此, 故作为评估对象的各区县, 应加强对突发事件的监测与预警工作, 完善各项应急准备措施真正做到防患于未然, 在最大程度上减少突发事件所造成的财产损失和人员伤亡。

图 4直观显示了各二级指标权重分布情况。结合各二级指标权重计算结果, 预案编制(P1)、预警信息收发(P7)、信息监测(P4)、预警研判(P5)和报警流程(P6)这5个二级指标的权重系数较高, 是突发事件应急管理过程中应重点关注的环节, 会对应急管理能力产生较大影响。重建计划(P12)、总结报告(P14)和经济补偿(P11)等指标对应急管理能力的影响程度较低, 但在日常管理过程中仍不可忽视这些方面的建设工作。

Download:
图 4 二级指标权重分布

表 5可知, 根据突发事件应急管理能力水平对各评估对象进行排序为D区、A县、C县、B县, 其贴进度如图 5所示。D区应继续保持其在应急管理方面的优势, 进一步完善应急防护各方面的工作, 真正做到防患于未然。而B县作为该市应急管控建设的薄弱地区, 应引起各方重视。该市在今后的应急管理建设工作中, 应加大对B县的政策及资金倾斜力度, 支持B县进一步完善在应急管理关键环节的建设, 优化应急预案, 提高预警信息收发速度, 简化报警流程等; 同时B县自身应加强面向群众的宣传工作, 提高社会的整体应急意识, 从而全方位提高B县的应急管理能力水平。

Download:
图 5 贴进度比较
4 结束语

本文提出基于D-AHP和TOPSIS的突发事件应急管理能力评估方法。通过引入D-AHP方法处理评估过程中存在的评估信息不完整、不确定等问题, 避免人为主观性对评估结果的影响, 同时识别出应急管理过程中的关键环节, 为完善应急管控建设提供针对性建议; 根据求得的各级指标权重, 结合TOPSIS方法对评估对象的应急能力进行排序, 直观展示各个评估对象的实际水平, 识别出应急管理建设较为薄弱的评估对象, 全面考虑各种因素对应急管理的影响。由于本文构建的评估指标体系未考虑指标间的联系, 因此下一步将对此改进, 并针对不完整信息进行智能评估。

参考文献
[1]
JIANG Zongcai. Evaluation of emergency risk management capability based on hesitant fuzzy Einstein operator[J]. Journal of Intelligent and Fuzzy Systems, 2016, 31(4): 2307-2311. DOI:10.3233/JIFS-169071 (0)
[2]
WU Minghua, ZHAO Jing, ZHANG Taowei. Campus emergency management capability based on fuzzy comprehensive evalua-tion[J]. Journal of Discrete Mathematical Sciences and Cryptography, 2016, 19(4): 863-874. DOI:10.1080/09720529.2016.1182797 (0)
[3]
CHEN Ye, LI K W, XU Haiyan, et al. A DEA-TOPSIS method for multiple criteria decision analysis in emergency management[J]. Journal of Systems Science and Systems Engineering, 2009, 18(4): 489-507. DOI:10.1007/s11518-009-5120-3 (0)
[4]
葛悦, 汪彤, 徐亚博. 基于灰色层次分析法的突发事件应急管理能力评价[J]. 中国安全生产科学技术, 2014, 10(12): 80-86. (0)
[5]
ZHOU Xinyi, SHI Yangqiuyan, DENG Xinyang, et al. D-DEMATEL:a new method to identify critical success factors in emergency management[J]. Safety Science, 2017, 91: 93-104. DOI:10.1016/j.ssci.2016.06.014 (0)
[6]
DENG Xinyang, HU Yong, DENG Yong.Bridge condition assessment using D numbers[EB/OL].[2018-09-15].http://dx.doi.org/10.1155/2014/358057. (0)
[7]
FAN Guichao, ZHONG Denghua, YAN Fugen, et al. A hybrid fuzzy evaluation method for curtain grouting efficiency assessment based on an AHP method extended by D numbers[J]. Expert Systems with Applications, 2016, 44: 289-303. DOI:10.1016/j.eswa.2015.09.006 (0)
[8]
HUANG Xiaozhong, WANG Ningkui, WEI Daijun.Investment decision using D numbers[C]//Proceedings of 2016 Chinese Control and Decision Conference.Washington D.C., USA: IEEE Press, 2016: 4164-4167. (0)
[9]
DEMPSTER A. Upper and lower probabilities induced by a multivalued mapping[J]. Annals of Mathematical Statistics, 1967, 38: 325-329. DOI:10.1214/aoms/1177698950 (0)
[10]
王坚强, 黄思旷. 基于模糊熵和证据推理的语言D数多准则决策方法[J]. 控制与决策, 2016, 31(4): 673-677. (0)
[11]
DENG Yong. D numbers:theory and applications[J]. Journal of Information and Computational Science, 2012, 9(9): 2421-2428. (0)
[12]
DENG Xinyang, HU Yong, DENG Yong. Supplier selection using AHP methodology extended by D numbers[J]. Expert Systems with Applications, 2014, 41(1): 156-167. DOI:10.1016/j.eswa.2013.07.018 (0)
[13]
DENG Xinyang, DENG Yong. D-AHP method with different credibility of information[J]. Soft Computing, 2019, 23(2): 683-691. DOI:10.1007/s00500-017-2993-9 (0)
[14]
王莲芬. 层次分析法中排序权数的计算方法[J]. 系统工程理论与实践, 1987(2): 31-37. (0)
[15]
吴诗辉, 刘晓东, 贾月岭, 等. 一种调整AHP不一致判断矩阵的优化方法[J]. 控制与决策, 2016, 31(11): 2106-2112. (0)
[16]
HWANG C L, YOON K. Multiple attribute decision making:methods and applications[M]. Berlin, Germany: Springer, 1981. (0)
[17]
郑双忠, 邓云峰. 城市突发公共事件应急能力评估体系及其应用[J]. 辽宁工程技术大学学报, 2006, 25(6): 943-946. DOI:10.3969/j.issn.1008-0562.2006.06.044 (0)
[18]
沈惠璋. 突发危急事件应急序贯群决策与支持系统[M]. 北京: 科学出版社, 2011. (0)