丁建闯,钟海仁,许礼林,余 洁
(1.首都师范大学资源环境与旅游学院,北京 100048;2.华北科技学院建筑工程学院,河北 三河 065201;3.国信司南(北京)地理信息技术有限公司,北京 100142)
减灾能力是综合风险管理的重要组成部分,科学评价减灾能力、摸清区域综合减灾能力底数、查明重点区域减灾能力,是提升灾害风险管理能力的有效途径[1]。如何采用科学、合理的评价方法来更好开展区域减灾能力评估,一直是减灾能力研究的关键,国内外学者在该领域相继开展了大量的研究工作。
在国内,单灾种减灾能力研究有借助灰色关联分析法、GIS技术、云模型、模糊综合(聚类)评价法、广义目标函数关系、层次分析法、熵权TOPSIS-PCA耦合优化模型、数据包络分析法等分别针对地震、干旱、火灾、泥石流、洪涝、风暴潮、气象等灾害的减灾能力进行了研究[2-14];综合减灾能力研究方面:借助模糊数学法、层次分析法、可变模糊聚类法、熵权-灰靶模型和GIS技术、综合评价法等分别针对城市综合减灾能力、脆弱性综合评价能力、城市灾害应急能力、区域减灾基础能力、社区减灾能力、区域综合减灾能力等进行了研究[15-24];灾害系统减灾能力研究方面:基于灾害群、灾害链、灾害遭遇等灾害系统情景,分别从理念与框架、灾害链适应对策与适应模式、多级模糊综合评估模型的角度对灾害系统减灾能力的评估进行了研究[25-27]。国外,美国是世界上最早、最成功将防灾减灾能力评价运用到实践的国家,20世纪80年代起就着手建立社区、州、国家三个层级的综合防灾减灾能力建设指标体系,于1997-2000年间应用研发的政府、企业、社区、家庭联动的灾害应急能力评估程序在全国开展了实践评估[28]。英国内阁办公室制订《关于形成社区减灾能力的战略框架》[29],细化了在形成社区减灾能力建设中个人、社区和其他参与者行为的指导原则。CUTTER等[30-31]提出了广泛应用于减灾能力定量评估的基于社会结构、经济条件、制度结构、基础设施和社区资本5个维度36项指标的社区减灾能力基线指标。KUNIHIRO等[32]通过阪神地震震害分析,认为城市抗震能力或城市地震易损性取决于城市自然条件、空间结构、社会特性、救灾资源和风险管理能力等多种因素影响,提出一种综合性、定量的城市抗震能力评价方法。PARSONS等[33]基于应急能力和适应能力构建了澳大利亚自然灾害减灾能力指数。国际上,欧盟FP7项目“城市防洪减灾能力协作研究”[34]通过量化减灾措施的成本效益,针对相关驱动因素量化了洪水管理计划。《2005-2015年兵库行动框架:加强国家和社区的抗灾能力》[35]和《2015-2030年仙台减轻灾害风险框架》[36]更是将加强减灾能力建设推向了新的维度和高度,提倡将减灾能力纳入可持续发展与灾害风险管理的过程中,强调国际协同减灾的理念。
减灾能力影响因素涉及社会经济、自然系统的多个领域,具有复杂性、动态性和相对性,其评估模型以多指标综合评价法最为常见,有定性评价法,如适用于战略层决策分析的专家会议法和德尔菲法;运筹学法,如适用于评价经济学中技术及规模有效性等的数据包络分析法;统计分析法,如适用于主体对象分类评价等的主成分分析、因子分析、聚类分析等;系统工程法,如适用于系统安全性评价、成本效益决策等的层次分析法、评分法、关联矩阵法等;模糊数学法,如适用于专家系统等的模糊综合评价、模糊模式识别、模糊积分等;智能化评价法,如适用于城市综合发展水平评价等的基于神经网络的评价[37-44]。而减灾能力评估实质属于评价指标的加权聚合过程,对其进行定量化评估属多项减灾能力的归一统算问题,应重点关注指标权重的制定与指标聚合过程[24]。前述可知,减灾能力评估模型并不唯一,国内外专家学者在减灾能力定量化评估方法方面并未达成一致。因此,还没有一个公认和普适的方法供大规模使用,有必要引入新的评估方法并基于已有方法来开展基于实例数据评估的深入对比分析研究;乡镇(街道)减灾能力是区域减灾能力的一种呈现形式,国务院第一次全国自然灾害综合风险普查领导小组办公室将其定义为:乡镇(街道)在灾害管理(应急预案、风险评估、资金投入)、防灾备灾(物资储备、医疗保障)、自救互救(专业和志愿者队伍、公众避险、转移安置)方面具备的各种能力,并力荐且已成功使用TOPSIS(Technique for Order Preference by Similarity to an Ideal Solution,优劣解距离法)法进行评估,该法适用于多组可比较样本数据,或与标准数据进行比较获得指标的归一化数值,再将各指标数值进行加权叠加,能最大化保留原始数据信息,指标归一化处理消除了不同量纲带来的影响,评估结果量化准确直观,但缺点是原始数据正负向标准化没有统一的标准,评估结果极易受少数极端值的影响,优劣分档欠佳[45-46]。同时,在全国范围的第一次全国自然灾害综合风险普查与综合减灾能力评估项目的开展、执行与协调中,相关专家结合本地区该法评估的减灾能力在大小、等级,及各层级综合减灾能力评估总结建议方面,对TOPSIS法提出了一些异议;RSR(Rank Sum Ratio,秩和比法)法是田凤调[47]提出的、在医药卫生领域应用较为成熟的一种多指标综合评价法,采用秩次参与评估计算能有效消除异常值干扰,且待评对象间的评估结果差距显著,应用时需将定量的指标值转换为半定量的秩次,即定量信息定性化(等级化),但用秩次代替原始指标值会丢失大量的源数据信息,即较多地损失了原始指标值提供的定量信息,无法体现源数据准确性;而对RSR法改进后的非整RSR法集成了两者的优点,为验证RSR法及非整RSR法的可行性与高效性,以TOPSIS法评估结果数据为参照,选取某市区18个乡镇的基础调查数据为例进行论证。基于此,在对RSR法进行改进与完善的基础上,以实际调查数据为例开展基于三种评估方法的比较研究,可为区域减灾能力定量化评估提供可行的较优方法。
以某市区“第一次全国自然灾害综合风险普查”数据为基础,该调查数据主要包括常住人口数量、本级灾害管理工作人员总数、是否开展乡镇(街道)灾害风险评估、上一年度防灾减灾救灾资金投入总金额、现有储备物资装备折合金额、实有住院床位数、消防员志愿者民兵预备役人数、上一年度组织的应急管理培训和演练参与人次、本级灾害应急避难场所容量等九类数据,减灾能力评估中的评估指标和权重数据即是基于此基础调查数据所派生出的。
乡镇减灾能力评估数据主要通过《乡镇与社区减灾能力调查技术规范》[48]调查获取,如医疗保障能力指标数据通过《公共服务设施调查技术规范》[49]医疗卫生机构表格调查获取,自救互救能力指标数据通过《政府减灾能力调查技术规范》[50]政府专职和企事业专职消防队伍表格调查获取,乡镇(街道)减灾能力的评估指标包括灾害管理能力、灾害备灾能力和自救转移能力三方面(表1)。
表1 乡镇(街道)减灾能力评估指标
每个指标对评估结果的影响不同,权重大小反映了各指标在综合评估中的重要性程度。权重确定方法总体可分为侧重专家经验的主观赋权法、侧重样本各属性指标间相关特征或内在数学关系的客观赋权法与综合考虑专家经验和客观评估数据的组合(综合)赋权法三种,使用中应综合考虑主观赋权法、客观赋权法和组合赋权法的应用条件,结合区域评估要求进行选取。根据国务院普查办印发的《评估指标权重制定技术规范》[51]要求,国务院普查办组织有关行业部门、省级普查办组建指标权重制定专家团队,集中完成减灾能力评估指标权重打分工作,前期制定了“乡镇(街道)减灾能力评估指标权重(参考值)”[46](表2、表3)。
表2 乡镇(街道)减灾能力评估指标权重
表3 乡镇减灾能力评估指标及权重数据
某市区乡镇减灾能力评估是在该区18个乡镇的评估指标和权重数据基础上,借助评估模型开展的。出于数据保密性,实际评估中用别名代替各乡镇名称,每个乡镇(用Ai表示)的灾害管理能力、灾害备灾能力和自救转移能力3个一级指标下队伍管理、风险评估、财政投入、物资储备、医疗保障、自救互救、公众避险、转移安置(用Xi表示)8个二级指标及权重值(用Wi表示)(表3)。
利用评估模型开展减灾能力评估,实质就是评估指标的加权聚合过程,在此借助以下三种模型方法。
秩和比是内涵极为丰富的一个统计量,表示不同计量单位多个指标的综合水平或指数,它是非参数统计向参数统计的“接口”,是行或列秩次的平均值,受分组数和指标数两项因素影响,具有0~1连续变量的特征。RSR法集参数统计与非参数统计于一身,用一个指标来综合反映各计量单位不同、功能类别各异的多个指标的综合水平,综合评价能力较强,可广泛应用于综合评价与鉴别分类、因素与关联分析、统计监控与预测决策等领域。RSR法评价原理是:首先在一个n行m列矩阵中,通过秩转换获得无量纲统计量RSR,再运用参数统计分析的概念与方法来研究RSR的分布,最后依据RSR值对评价对象进行排序或分档排序进而得出综合评价[47]。
1)编秩:将n个评价对象的m个评价指标列成n行m列的原始数据表或矩阵,依次编出各评价对象每个指标的秩。在指标编秩时,高优指标编秩法的最大指标值编以最高的秩次n,次大的编以n-1,……次小的指标值编以2,最小的编以1;低优指标编法则与高优指标相反,最小的指标值编以最高的秩次n,次小的编以n-1,……次大的编以2,最大的编以1;若几个指标值相同,则均编以平均秩次。
2)计算RSR:通常按行或列或加权形式计算秩和比。
若按行计算:
(1)
若按列计算:
(2)
若按加权秩和比计算:
(3)
式中:R为各评价对象每个指标的秩;W为各指标权重。
3)确定RSR的分布或计算概率单位(Probit):根据累计频率,对照“百分数与概率单位对照表”,求其所对应的Probit值。
4)分档排序:0≤RSR≤1,RSR→1表明评价对象越优。
非整RSR法(非整秩次秩和比法)是对RSR法的进一步改进和完善。高优指标编秩中,最小指标值的秩编为1,最大指标值的秩编为n,其余指标值的秩由小到大分别编为1与n之间的线性递增的非整秩次,所编秩次与原指标值之间存在定量的线性对应关系,即原指标值被定量地转换为秩次,而不是生硬的等级化;低优指标的编秩原理与高优指标相同,但大小方向相反,非整RSR法编秩计算见式(4)和式(5)。
对高优指标编秩:
(4)
对低优指标编秩:
(5)
式中:X为各评价指标;Xmax、Xmin分别为指标最大、最小值。
非整RSR法对RSR法的编秩方法进行了优化,避免了秩次化后原指标值定量信息的损失,尤其在待评区域同时具有高优指标、低优指标、不分高低及其组合指标时,非整RSR法最大化保留源数据信息的优势更加明显。除编秩方法外,非整RSR法与RSR法的计算步骤几乎相同,其过程是在将n个评价对象的m个评价指标列成n行m列的原始矩阵后,判定高优、低优、不分指标,首先按式(4)或式(5)依次编出各评价对象每个指标的秩,再按式(1)、式(2)或式(3)计算出非整RSR值。
TOPSIS法是一种常用的组内综合评价法,它能充分利用原始调查数据信息,在与标准数据进行比较获得指标归一化数值后,再将各指标数值进行加权叠加的过程,优点是评估结果量化准确、直观。第一次全国自然灾害综合风险普查推荐利用TOPSIS法开展评估,该法对数据分布及样本含量没有严格限制,数据计算简单易行,着重体现评估单元能力与目标值的差异,评估结果能准确反映各评价区域或单元间的减灾能力差距。TOPSIS法评价原理是:首先基于归一化后的标准化矩阵,采用余弦法找出有限方案中的最优方案和最劣方案,再分别计算各评价对象与最优方案和最劣方案间的距离,进而获得各评价对象与最优方案的相对接近程度,以此作为评价优劣的依据[45-46]。
1)构建标准化矩阵,假设有n个评价对象(起算单元),每个对象都有m个属性指标,则原始数据矩阵为:
(6)
进行属性归一化:
(7)
根据专家知识经验判定法,得到每个指标的权重Wj,进行向量定权:
Zij=Wj×Pij。
(8)
构造加权规范矩阵,求得归一化和定权后的标准化矩阵Z:
(9)
2)确定最优方案和最劣方案,最优方案Z+由Z中每列元素的最大值构成:
(10)
最劣方案Z-由Z中每列元素的最小值构成:
(11)
3)计算各评价对象与最优、最劣方案的距离:
(12)
(13)
4)计算各评价对象与最优方案的贴近程度Ci:
(14)
式中:0≤Ci≤1,Ci→1表明评价对象最优。
5)根据Ci大小进行排序,给出评价结果。
基于RSR模型评估原理及计算过程,某市区乡镇减灾能力评估计算与结果见表4。
表4 RSR法评估RSRW值分布
基于非整RSR模型评估原理及计算过程,某市区各乡镇减灾能力评估计算与结果见表5。
表5 非整RSR法评估RSRW值分布
基于TOPSIS模型评估原理及计算过程,某市区各乡镇减灾能力评估计算与结果见表6。
表6 TOPSIS法评估Ci值分布
将表4、表5与表6中的数据进行汇总,三种模型评估的结果数据及排名对照见表7。
表7 乡镇减灾能力评估结果对照表
将表7中的评估结果数据分别绘制成图1、图2和图3,由图中趋势线的斜率及极值差异可知,基于不同模型单个乡镇或整体区域的评估结果差异显著。但从评估效果来看,三种方法都能对乡镇减灾能力进行定量化评估,采用任一方法结合均值-标准差等分级法均能对各乡镇减灾能力大小进行有效地分等定级。
图1 TOPSIS法评估结果图
图2 RSR法评估结果图
图3 非整RSRW法评估结果图
从评估结果数值来看,基于相同的评估指标和权重,同一乡镇的三种模型评估结果均存在不同程度的差异,区域整体上以RSR法评价值最大、TOPSIS法评价值最小、非整RSR法评价值居中(图4)。
图4 三种方法评估结果数值对比图
对评估结果深入分析可知,TOPSIS法评估结果的平均值为0.285、极值差为0.283、标准差为0.115,均为三者中最小;非整RSR法评估结果的平均值为0.354、极值差为0.382、标准差为0.134,处于居中;RSR法评估结果的平均值为0.528、极值差为0.472、标准差为0.141,均为三者中最大(表8)。表明RSR法下各乡镇间评估结果差异明显、区分度最大,即各乡镇优劣分档的效果最佳。
表8 评估结果数据处理分析
从评估对象排名来看,RSR法与TOPSIS法的评估结果差异最大,18个乡镇中,减灾能力排名差距为2级的为A2、A4、A5、A7、A9、A12、A13、A17,排名差距为3级的为A6、A14、A18,排名差距为4、5、6、10级的分别为A10、A8、A1、A3;其次为非整RSR法与TOPSIS法,减灾能力排名差距为2级的为A5、A9、A18,排名差距为3级的为A4、A6、A10、A12、A13,排名差距为11级的为A3;非整RSR法与RSR法排名差异相对较小,减灾能力排名差距为2级的为A14、A15,排名差距为4级的为A8、A9,排名差距为5级的为A1、A12(图5)。
图5 三种方法评估结果排名对比图
从评估结果整体趋势来看,非整RSR法评估结果整体处于TOPSIS与RSR法之间,该法既吸取了TOPSIS法最大化保留源数据信息的优点,又继承了RSR法较好的优劣分档的优势(图6)。
图6 三种方法评估趋势图
对评估结果数据进行深入对比分析,并结合三种模型的评估原理,可知:
1)实用性方面。非整RSR法与RSR法评估差异相对较小,两种方法与TOPSIS法的差异显著。评估指标及权重数据虽相同,由于三种方法的计算方式不尽相同,评估结果差异明显,但三种方法均可实现区域减灾能力定量化评估,仅对各乡镇减灾能力评估值整体大小进行对比意义不太明显,而基于各乡镇减灾能力的排名差距,探究其主导因素并对应提出提升减灾能力的措施意义重大。
2)数值大小方面。以TOPSIS法评估计算的数值最小、标准差最小;RSR法评估计算的数值最大、标准差最大;非整TOPSIS法评估计算的数值居中、标准差居中。基于三种方法的评估原理及式(1)~(14),表明RSR法降低了异常指标的干扰、优劣分档最佳;TOPSIS法最大化保留了源数据信息,但易受异常指标干扰;非整RSR法继承了两者的优点,是三者中既重视源数据信息又强调排名的较好方法。
3)等级排名方面。非整RSR法与RSR法的排名相差2级以上的乡镇有6个,最大相差5级;非整RSR法与TOPSIS法排名相差2级以上的乡镇有9个,最大相差11级;RSR法与TOPSIS法排名相差2级以上的乡镇有15个,最大相差10级。
4)整体趋势方面。非整RSR法居于TOPSIS法、RSR法两者之中,以非整RSR法为最优。
5)计算工作量方面。由三种模型评估公式及实际计算知,TOPSIS法计算工作量最大,是另两种方法的几倍,非整RSR法与RSR法的计算工作量相当。
经采用RSR法、非整RSR法与TOPSIS法对某市区各乡镇减灾能力进行定量评估与对比分析,主要结论如下:
1)三种方法均适用于区域减灾能力定量化评估,采用任一种方法均可量化出各区域分异的减灾能力,应根据区域评估目的和要求进行选择。
2)RSR法优劣分档最佳、信息损失最大;TOPSIS法信息损失最小、优劣分档最差;非整RSR法信息损失和优劣分档均居中,为三者中的较好方法。
3)从实用性、数值大小、等级排名、整体趋势、计算工作量等五方面综合考量,非整RSR法评估模型最为理想。
4)三种方法评估结果产生差异的内在机理:借助评估模型进行层层计算时,对指标数据的近似取值程度或处理方式,及指标与权重的加权聚合过程不同,如TOPSIS法的式(6)~(14)和RSR法的式(1)~(5)所示;各乡镇单一指标极端值对该指标及18个乡镇减灾能力值的影响与贡献程度不同,如A3、A6、A10、A18在三种方法下评估数值及排名的差异等。
因此,后续研究应对三种模型下减灾能力排名差距较大的主导诱因进行深入分析。