李正兆^1,2, 傅大放¹, 王君娴^1,3, 闵克东¹, 张骏彧^3,4
1. 东南大学 土木工程学院，南京 210096，中国;
2. 常州市规划设计院，常州 213000，中国;
3. 东南大学-蒙纳士大学未来城市联合研究中心，苏州 215123，中国;
4. 荷兰国际水利环境工程学院 水科学与工程系，代尔夫特 2611AX，荷兰
收稿日期：2021-03-29
基金项目：国家“十三五”重点研发计划(2018YFC0809900)
作者简介：李正兆(1983—), 男, 博士研究生
通讯作者：张骏彧, 博士后, E-mail: j.zhang@un-ihe.org

摘要：传统的防灾减灾理念下, 内涝风险评估通常只考虑基础设施的抵御能力和气象、地形等自然条件。随着系统韧性理论的发展, 政府和市民的防控意识、应急处理和灾后恢复等方面均需要纳入到风险评估和对策之中。该文提出了一套用于评估城市应对内涝灾害的综合能力的指标体系, 包含13个方面共41个基础指标, 并在该指标体系基础上构建了韧性评估模型; 以内涝风险较高的圩区城市江苏省昆山市为例, 进行了韧性评估模型的应用, 结果表明昆山市高新区韧性最高, 应对内涝灾害的能力最强。结合指标体系和评估结果, 对韧性较差区域提出了针对性提升策略。
关键词：防灾减灾内涝灾害城市韧性韧性评估
Urban resilience assessment model for waterlogging disasters and its application
LI Zhengzhao^1,2, FU Dafang¹, WANG Junxian^1,3, MIN Kedong¹, ZHANG Junyu^3,4
1. School of Civil Engineering, Southeast University, Nanjing 210096, China;
2. Changzhou Institute of Planning and Design, Changzhou 213000, China;
3. SEU-Monash Joint Research Center for Future Cities, Suzhou 215123, China;
4. Department of Water Science and Engineering, IHE Delft Institute for Water Education, Delft 2611AX, The Netherlands

Abstract: In traditional disaster prevention and mitigation analyses, risk assessments for waterlogging disasters usually only consider the infrastructure resilience and natural conditions such as weather and terrain. However, with the development of resilience theory, risk assessment and countermeasures should also include flood prevention and control, emergency responses and post-disaster recovery by the government and residents. This paper presents a set of indicators for evaluating the ability of cities to deal with waterlogging disasters, including 41 basic indicators covering 13 aspects incorporated into a resilience assessment model. The resilience assessment model is then applied to Kunshan, a low lying city in Jiangsu Province, China. The results show that the Kunshan High-Tech Zone has the best resilience and the strongest ability to deal with waterlogging disasters. The index system and model evaluation results are used to develop a targeted improvement strategy for areas with poor resilience to waterlogging disasters.
Key words: disaster preventionwaterlogging disastersurban resilienceresilience assessment
近几十年，在全球气候变暖的背景下，极端天气所引发的自然灾害事件频繁发生，严重影响城市的公共安全。其中内涝灾害居于前列，中国约有50%的人口和70%的财产暴露在内涝的威胁下^[1]。为了更好地应对以内涝为代表的城市灾害，韧性城市的理念被提了出来^[2]，人们对城市内涝灾害风险及其防治也随着该理念的不断发展有了新的认识。洛克菲勒基金会^[3]认为韧性城市强调人、社区、机关、企业、系统在遭受长期压力和突发性的动荡时，能够生存、适应及增长。黄弘等^[4]将韧性城市定义为：能够有效应对来自外部与内部的经济社会、技术系统和基础设施的冲击和压力，能在遭受重大突发事件后维持城市的基本功能、结构和系统，并能迅速恢复、进行适应性调整、可持续发展的城市。因此，在面对内涝灾害时，城市不仅需要单纯地防灾，还应该具备适应内涝和快速恢复的能力。在涉及风险评估时，不能仅考虑基础设施的排涝能力，还需要考虑政府的应急能力、城市各系统的恢复速度^[5]、市民的防灾减灾意识等更多社会层面的指标。
现阶段，已有国内外****探索应对内涝灾害的城市韧性评估模型。陈长坤等^[6]从抵抗力、恢复力和适应力三大属性出发，建立了雨洪情境下的韧性评估模型；郑艳等^[7]构建了评估城市发展、绿色基础设施和灰色基础设施应对内涝灾害能力的韧性指数；程元^[8]研究了暴雨灾害下城市韧性的测度问题，选取与内涝最有直接关系的指标构建了陕西省韧性评价体系；Balica等^[9]、Rosenzweig等^[10]、Mugume等^[11]则侧重于研究雨洪预测模型以考量城市内涝韧性。也有很多****将研究主体缩小到社区，Orencio等^[12]构建的韧性指标评估模型包含7个方面；Qasim等^[13]从社会、经济、制度和物理4个层面确定了社区洪水灾害韧性指标；李彤玥^[14]分析了韧性城市的理论演变及特性，并分别构建了应对洪水的韧性城市评价指标体系和基于社会、经济、制度、基础设施4个维度的社区韧性评价指标体系；姜宇逍^[15]明确了韧性社区与海绵社区的区别，建立了三级指标体系。但是，现有的内涝灾害相关的韧性评估体系大多为灾损评估，过度重视客观因子，却忽视主观因子；已有的内涝灾害韧性评估研究虽然考虑了工程-非工程相结合，但是仍然局限于顶层政策层面，而缺乏对底层策略的研究，即缺乏对抗灾实施者、承灾受害者的提升指导。
在模型算法方面，定量评估是韧性评价的主要方法，即通过构建指标体系、选用相应的数学方法、建立运算模型来达到量化的目的。国内外使用最为广泛的是层次分析法(analytic hierarchy process, AHP)，即对指标进行分层，通过专家咨询确定同一层中各指标的重要性，计算指标所占权重，最后通过层层递进的加权计算得出最终的韧性值。Orencio等^[12]采用AHP计算了沿海社区的韧性指数；Qasim等^[13]采用专家打分法计算了巴基斯坦开伯尔-普赫图赫瓦省3个地区的社区韧性指数；姜宇逍^[15]运用AHP-综合模糊分析法构建了社区内涝灾害韧性评价体系; 邱桐等^[16]采用AHP确定主观权重，采用熵权法(entropy weight method, EWM)确定客观权重，研究了城市地下空间的综合韧性。此外，陈长坤等^[6]利用逼近理想解排序法(technique for order of preference by similarity to ideal solution based on Kullback-Leibler divergence, KL-TOPSIS)对武汉市的韧性进行评估；苏伯尼等^[17]通过二维水力学模型模拟积水的时空分布，建立了一套针对城市暴雨内涝灾害的定量风险评估方法，但与前述基于AHP方法的研究相比，其指标涉及面相对较少。
本文遴选出影响中国城市内涝的致灾、抗灾、承灾三类特征指标，形成一套能完整描述城市各层面应对内涝灾害的韧性指标体系，通过可视化方法，展示各指标在适应性、应对性和恢复性上的性能。结合地理信息系统(geographic information system, GIS)和AHP，建立了城市内涝灾害韧性定量评估模型。最后，以江苏省昆山市为例开展应用研究，结合分析结果提出韧性强化策略。本研究成果为城市应对内涝灾害韧性理论的完善提供了数据支撑，对提高城市抗涝减灾能力具有重要的现实意义。
1 指标体系构建在参考国内外已有指标体系的基础上，结合中国现行标准，如《全国综合减灾示范社区标准》《海绵城市建设指标体系》等，以城市韧性作为研究主体，以“应对灾害、恢复能力、适应变动”作为研究目的，考虑内涝灾害致灾因素复杂、抗灾设施单一、承灾系统多维的特点，进行指标初选。现有韧性城市评价及内涝风险评价体系涉及降雨、地形、地貌、排水设施、次生灾害、避灾设施、灾情传递、灾害应对、医疗条件、生活保障、产业经济、受灾对象、就业恢复这13个方面。在此基础上，依据独立性、全面性、可空间量化性原则，筛选设计了41个基础评价指标，如表 1所示。
表 1 指标选取说明

类别	说明	基础评价指标	参考文献
降雨	内涝形成的最主要原因是因为降雨，内涝灾害的严重程度与降雨量的大小、降雨持续时间、暴雨强度、降雨的分布情况等均有关系	年均降雨量
地形	地势是造成内涝并影响内涝严重程度的最主要因素	绝对高程、相对高程	[9, 18-21]
地貌	地貌与内涝的形成也有较大的关系，地貌一般与土地、河流等有关	径流系数、水面率
排水设施	由于研究对象为内涝灾害，排水设施主要是指排涝方面的指标，区域的排涝能力与径流的流量有很大的关系	排涝标准、排涝泵站使用年数	[22-23]
次生灾害	次生灾害是影响内涝灾害损失的一项不可忽视的指标，很多灾害都会引发次生灾害，因此考察次生灾害的潜在风险也是非常必要的	重污染企业数、河网密度	[10, 24-25]
避灾设施	一座城市的避灾设施在灾害来临时尤为重要，避灾设施建设与维护得越好，灾害所造成的损失就越小	道路密度、万人拥有公共汽车数、每百户拥有私家车数、避难场所可容纳人数、星级宾馆床位数	[10, 21, 23, 26]
灾情传递	灾情传递的速度决定了灾损的严重程度	水文水位站点密度、互联网用户数、手机用户数	[10, 21-22, 24, 26]
灾害应对	灾害应对指的是城市居民及社区在灾害来临前、来临时、来临后所作出的反应，拥有良好的灾害应对能力是减少灾损的有效途径	居民减灾意识、居民减灾知识、居民减灾行为、社区减灾意识	[9, 21, 23]
医疗条件	医疗条件是影响灾后恢复的最直接因素，医疗条件越好，灾后得到救治的人越多，城市恢复得越快	医疗保险覆盖率、人均医疗财政投入、万人拥有的病床数、万人拥有卫生技术人员数、心理援助设施数	[9, 10, 22, 24, 26]
生活保障	生活保障是指在灾害发生后可以暂时满足城市居民基本生活需求的因素	农作物播种面积、大型零售店数	[10, 21-22]
产业经济	一座城市产业的构成、经济的发达程度决定了城市受到灾害创伤的程度以及恢复能力	单位面积GDP、农业总产值占比、第三产业产值占比、非私营单位从业人员数	[20-21, 23-24, 26]
受灾对象	受灾对象主要分为居民和社区，因此人口的构成和社区的抗灾能力是影响城市恢复的关键因素	综合减灾示范社区数、人口密度、小学以下文化程度人口占比、旅游人次、60岁以上人口占比、最低生活保障人口占比	[9, 20-21, 23-24, 26]
就业恢复	就业恢复也是灾后恢复的重要组成部分	人均可支配收入、每个就业者平均负担人数、失业率	[20-21]

表选项






将这13个方面分别对应至内涝灾害致灾、抗灾、承灾过程。致灾因子表示的是影响形成径流的自然因素，包含降雨、地形和地貌；抗灾因子表示的是影响灾损的人为因素，包含排水设施、次生灾害、避灾设施、灾情传递、灾害应对；承灾因子表示的是影响灾后恢复的人为因素，包含医疗条件、生活保障、产业经济、受灾对象和就业恢复。再将致灾因子、抗灾因子、承灾因子融合为内涝韧性指数，体现城市应对内涝灾害的综合能力。具体指标和层次结构如图 1所示。

图 1 城市应对内涝灾害的韧性指标体系

图选项

本套指标体系是基于韧性城市的“适应性、应对性、恢复性”三大核心特性构建的。适应性既体现为城市根据外界环境的变化而调整、适应的能力^[27]，也体现在其理解风险并在规划建设过程中进行准备以更好地与未来风险相适应的能力；应对性体现城市对灾害事件进行预防准备、监测预警并在事件发生后应急处置的能力；恢复性体现城市维持业务连续性水平并从灾害事件影响中迅速恢复的能力^[4]。
由于三类因子及其细化指标与三大核心特征并非简单的一一对应关系，传统指标体系较难直观地将每个指标与城市韧性进行联系。为立体展现各指标在应对内涝灾害不同阶段、不同层面的韧性表现，设计了城市应对内涝灾害的韧性概念图，如图 2所示。

图 2 城市应对内涝灾害的韧性概念图

图选项

致灾因子中，降雨类指标是城市内涝灾害风险基准的体现，影响城市在规划建设过程中应具备的基础适应力；地形类指标是城市自身特性和纵向规划建设的主要体现，既影响城市在规划建设过程中所具备的可调整性，也影响灾害事件中的预防准备能力；地貌类指标是城市下垫面规划及建设的主要体现，既影响内涝灾害中形成的径流总量(风险基准)，也影响灾害中径流调蓄及灾害后积水排空能力。
抗灾因子中，排水设施类及灾害应对类指标分别是城市抗灾设施和措施的体现，合理的设施建设和民众应对可以有效提高城市处置内涝事件的能力，针对未来风险而开展的设施规划和应对提升也可以一定程度提高城市的适应性；避灾设施类指标是城市对受灾民众安置能力的体现，既影响城市对受灾群体的保障能力，也影响民众自身从灾害中恢复的速度；次生灾害类及灾情传递类指标分别是城市对灾害预防准备和监测预警能力的体现，影响内涝事件发生时的应急处置能力。
承灾因子中，受灾对象类、产业经济类、生活保障类及医疗条件类指标既在灾害事件中影响城市受灾程度和救援能力，也在灾害发生后与就业恢复类指标共同影响民众生活及城市经济的恢复速度；此外，部分生活保障类指标(如农作物面积)在内涝这一特定灾害中也影响城市下垫面特性，因而对城市适应力也有一定作用。
2 韧性评估模型框架及数据处理方法2.1 数据标准化对于矢量数据，需要在已有行政区的矢量图中自行添加字段进行数据的录入和处理。数据表是数据库中重要的数据管理对象，每一个数据表均包含一个具体主题的信息，所有字段都无空值。
数据库的建立是基于各指标原始的意义和计量单位，通常每个指标原始数据的量纲和数量级各不相同。当各指标间的水平相差很大时，直接用原始数据进行分析会突出数值较高的指标在综合分析中的作用，而削弱数值较低指标的作用。因此，为了保证分析结果的可靠性，需要对各指标的原始数据进行标准化处理。在遴选出的41项指标中，有些指标对城市韧性起积极作用，称之为正向指标；有些指标对城市韧性起消极作用，称之为负向指标。采用极值法(离差标准化)对原始数据进行标准化处理。对于正向指标，

$Y_{i}=\frac{x_{i}-\min x_{i}}{\max x_{i}-\min x_{i}}.$

(1)

对于负向指标，

$Y_{i}=\frac{\max x_{i}-x_{i}}{\max x_{i}-\min x_{i}} .$

(2)

式中：x_i为指标原始数据，Y_i为标准化处理后的数据^[28]。
2.2 模型架构图 3展示了城市应对内涝灾害的韧性评估模型架构。

图 3 模型流程图

图选项

对城市韧性的综合评估基于加权平均数计算，

$R=\omega_{1} Y_{1}+\omega_{2} Y_{2}+\cdots+\omega_{n} Y_{n}.$

(3)

式中：R为每层次计算得到的分数，ω_i为权重。
韧性评估模型共分为3个模块：指标录入模块、指标处理模块和韧性评估模块。指标录入模块为模型的基本输入模块，基于前期调研分析，判断各指标的分辨率即指标统计的口径是以格网为单位还是以行政区为单位。例如，河网密度、水面率等可以从各类数据库里下载获取分辨率较高的栅格数据；而交通、医疗等数据的统计口径多为行政区甚至是全市，就只能获得矢量数据。通过ArcGIS的转换功能将矢量数据转化成栅格数据，实现指标基本单元的统一。指标处理模块对无法直接反映指标的数据进行空间化处理，并对所有数据进行标准化处理得到Y_i。之后，将数据输入至韧性评估模块，使用AHP方法计算每个指标的权重ω_i，通过对每一层加权平均值的计算，得到栅格层面城市应对内涝的韧性值，再转化为基于城市行政区划的韧性值。
2.3 权重确定采用层次分析法确定各指标的权重。邀请防灾领域、城市规划领域、土建领域、社会学领域的15位****、从业人员和管理人员对指标体系中同一层次的指标的重要程度进行打分，取平均值后选择最接近的标度。根据影响内涝灾害的主要因子及其相互间的属性关系，构造层次结构模型，由上到下分为目标层、准则层和指标层。本套指标体系中，目标层为一级指标，准则层为二级指标，指标层为三级和四级指标。对同一层次要素的重要性进行两两比较，构造判断矩阵，利用数学方法计算矩阵的特征向量值，最大特征向量值对应的列矩阵即为每一个指标的权重。应用MATLAB 2011b对所有指标进行权重计算，结果如表 2所示。表 2中，权重表示本指标对上级指标的贡献度，组合权重表示本指标对一级指标(城市韧性)的贡献度。
表 2 应对内涝灾害的韧性评估指标权重汇总

准则层				指标层1				指标层2
二级指标	权重	组合权重		三级指标	权重	组合权重		四级指标	权重	组合权重
致灾因子	0.493 4	0.493 4		降雨	0.200 0	0.098 7		年均降雨量	1.000 0	0.098 7
				地形	0.200 0	0.098 7		绝对高程 相对高程	0.500 0 0.500 0	0.049 3 0.049 3
				地貌	0.600 0	0.296 0		径流系数 水面率	0.750 0 0.250 0	0.222 0 0.074 0
抗灾因子	0.310 8	0.310 8		排水设施	0.409 4	0.152 4		排涝标准 排涝泵站使用年数	0.833 3 0.166 7	0.127 0 0.025 4
				次生灾害	0.053 4	0.016 6		重污染企业数 河网密度	0.166 7 0.833 3	0.002 8 0.013 8
				避灾设施	0.097 4	0.030 3		道路密度 万人拥有公共汽车数 每百户拥有私家车数 避难场所可容纳人数 星级宾馆床位数	0.424 9 0.273 3 0.175 3 0.087 6 0.038 9	0.012 9 0.008 3 0.005 3 0.002 7 0.001 2
				灾情传递	0.183 4	0.057 0		水文水位站点密度 互联网用户数 手机用户数	0.400 0 0.200 0 0.400 0	0.022 8 0.011 4 0.022 8
				灾害应对	0.256 4	0.079 7		居民减灾知识 居民减灾意识 居民减灾行为 社区减灾意识	0.166 6 0.166 6 0.333 4 0.333 4	0.013 3 0.013 3 0.026 6 0.026 6
承灾因子	0.195 8	0.195 8		医疗条件	0.160 0	0.031 3		人均医疗投入 医疗保险覆盖率 万人拥有病床数 万人拥有卫生技术人员数 心理援助设施数	0.120 4 0.316 8 0.316 7 0.182 7 0.063 4	0.003 8 0.009 9 0.009 6 0.005 7 0.002 0
				生活保障	0.262 5	0.051 4		农业播种面积 大型零售店数	0.333 3 0.666 7	0.017 1 0.034 3
				产业经济	0.097 2	0.019 0		单位面积GDP 农业总产值占比 第三产业值占比 非私营从业人员	0.563 7 0.109 5 0.069 2 0.257 6	0.010 7 0.002 1 0.001 3 0.004 9
				受灾对象	0.418 5	0.081 9		综合减灾示范社区占比 人口密度 小学以下文化人口占比 旅游人次 60岁以上老年人口占比 最低生活保障人口占比	0.231 1 0.365 7 0.051 5 0.080 3 0.135 7 0.135 7	0.018 9 0.030 0 0.004 2 0.006 6 0.011 1 0.011 1
				就业恢复	0.061 8	0.012 1		人均可支配收入 每个就业者平均负担人数 失业率	0.297 0 0.539 6 0.163 4	0.003 6 0.006 5 0.002 0

表选项






2.4 研究区域及数据类型选取江苏省昆山市作为研究区域。该区域地处太湖流域下游的阳澄淀泖区，河网密布，地势平坦，自然坡度较小，内涝灾害频繁。昆山市内涝灾害主要由梅雨和台风暴雨两种雨型所致。梅雨历时长、分布广、雨量大，多发生在6~7月间，灾情较为严重的1954、1991、1999年洪涝灾害都是由梅雨造成的。台风是昆山市夏秋季节主要的灾害性天气之一，当台风登陆过境并伴有强暴雨发生时，对暴雨集中区的雨水管网及某些圩区来说，极易造成较严重的灾情。本文通过各指标数据的采集及综合评估，分析昆山市应对内涝灾害的韧性状况，识别薄弱环节，并针对性地设计韧性提升策略。
根据指标统计口径的不同，原始数据可分为栅格数据和矢量数据。栅格数据大多来自自然属性的指标，如降雨、地形、地貌等，矢量数据大多来自社会属性的指标，如灾害应对、灾情传递等。图 4展示了昆山市的行政区域划分。行政区域是大部分矢量数据的基本评估单元。

图 4 昆山市行政区域划分

图选项

1) 栅格数据。
栅格数据的处理方式依具体的指标而定，大多会应用到ArcGIS的裁剪、转换、栅格计算器等功能。以径流系数为例，从中国科学院资源环境科学数据中心下载江苏省的土地利用图，通过裁剪获得昆山市的土地利用栅格图，涉及11类土地利用类型：水田、旱地、其他林地、高覆盖草地、河渠、湖泊、水库、城镇用地、农村居民地、沼泽地和裸土地。通过栅格计算器功能，利用条件函数，将11类土地利用类型转化为各地块的径流系数。图 5a以昆山市各地块径流系数为例，展示了栅格数据的评估结果。其中，径流系数的确定参考《昆山市防洪排涝规划》。

图 5 指标层处理结果栅格图

图选项

2) 矢量数据。
矢量数据一般从政府官网、城市统计年鉴、社会调研等途径获得，处理方式基本相同，即在建立的数据表中输入每个行政区的原始数值，通过简单的四则运算获得指标值，再通过要素转栅格功能得到栅格数据。
社会调研数据以居民减灾知识为例。首先设计包含内涝灾害、预警预防、应急处置等相关知识的问卷(含单选题和多选题)。在昆山市发放线上和线下问卷共600份，了解各区域居民对内涝灾害应急措施的掌握程度。剔除无效问卷(8份)后，进行计分统计(单选题正确得1分，错误得0分；多选题漏选按比例得分，错选得0分)。以行政区为单位，计算正确率(平均得分/总分)，取max x_i=1、min x_i=0进行标准化，得出各区域居民减灾知识掌握程度，录入矢量数据表中，最后转化为栅格数据。图 5b以居民减灾知识为例，展示了矢量数据转化为栅格数据后的评估结果。
3) 准则层和目标层数据。
应用ArcGIS的栅格计算器工具对上述数据进行标准化处理，将所有数据进行量纲归一化，之后将所有指标的权重代入栅格计算器中，计算准则层和目标层的数值。
3 结果及讨论图 6反映了昆山市致灾、抗灾、承灾三大因子的综合评估结果。

图 6 准则层评估结果栅格图

图选项

致灾性能反映了城市发生内涝灾害的难易程度，是一个反向指标。从图 6a可以看出，昆山市中心城区的致灾性能最高，其韧性最低，这主要是由于中心城区的地势相对比较平坦，高程变化幅度小；就地貌而言，中心城区多为建设用地和居住用地，且建筑密度大，其径流系数也大，降雨无法就地入渗，形成内涝灾害的概率也相对较大。此外，水体附近的致灾性能低，这主要是由于水面率大，有足够的水体对强降雨进行调蓄。
抗灾性能反映的是区域抵抗内涝灾害的能力，基础设施的建设和居民应对灾害的能力这两方面尤为重要。图 6b中颜色较深且可以明显看出行政边界的是高新区。高新区的抗灾性能较强，从城市化进程角度来看，高新区一直是昆山市发展的重要区域，其基础设施建设完善，没有河流、重污染企业等次生灾害风险；从居民应对的角度而言，高新区的居民在减灾知识、减灾意识和减灾行为上均是各行政区中最好的。
承灾性能主要体现区域在内涝灾害来临后的恢复能力，经济实力越强的地区承受灾害的能力也越强，灾害所带来的损失越能够被接受。承灾性能涉及的指标大多以行政区划作为统计口径，从图 6c中可以明显看出行政边界。对于高新区、开发区、花桥镇等经济实力较好的行政区，医疗条件充足、多以二三产业为主、生活方面的保障设施众多、人均可支配收入高、文化层次高，其韧性自然就高。其他行政区大多是乡镇，以农村人口为主，虽然其人口密度小，但是人口的文化层次低、老龄化严重，且经济发展较差；内涝发生时，可能会出现预警信息接收不及时、转移困难等问题，故其韧性低。
由致灾性能、抗灾性能和承灾性能乘以各自的权重得到昆山市应对内涝灾害韧性分布图，如图 7a所示。将该栅格数据转为要素点，并和昆山市行政区划进行空间连接，则这些代表韧性的要素点就分布在昆山市的各行政区内。经过数据汇总，得到每个行政区内的韧性平均值，将各区的值等比扩大为0~1之间的数值，则可以得出每一个行政区综合相对韧性值(韧性值最高的行政区相对韧性值为1)，如图 7b所示。

图 7 昆山市应对内涝灾害韧性分布图

图选项

图 7结果表明，昆山市中部的韧性最好，其次是北部和东部，韧性最差的是南部。其中，韧性最高的是高新区，然后依次是开发区、花桥镇、陆家镇、周市镇、千灯镇、周庄镇、张浦镇、巴城镇、淀山湖镇、锦溪镇。
本研究通过从综合到细化的四级指标的分层评估，可以高效定位城市应对内涝韧性的薄弱区域和薄弱环节，从而有针对性地设计提升策略。在设计提升策略时，既可以通过进一步强化优势指标来带动相关联的劣势指标提升，也可以通过多个劣势指标协同改进，实现互补；同时，兼顾各四级指标的组合权重(对城市韧性的贡献度)，可以优先考虑高权重指标的改善。
以锦溪镇锦南圩为例。其行政区韧性为0.64，圩区韧性为0.51，属于韧性较差地区；但其组合权重最高的年均降雨量(0.098 7)、径流系数(0.222 0)、排涝标准(0.127 0) 3个指标与其他地区差异不大。虽然政府已在长期规划中通过常规的城市整体布局、排涝能力提升等方式进行改善，但由于周期长、投资大、提升有上限，因而不能持续有效地提升该区域的韧性。
考虑四级指标组合权重平均值(0.025 0)与中值(0.011 1)，结合该圩区指标评估情况，从居民、城市、政策三方面挖掘区域韧性提升潜力，涉及以下韧性提升策略：
1) 该区域社区减灾意识(组合权重0.026 6)较好，但居民减灾意识(组合权重0.013 3)较差，减灾知识(组合权重0.013 3)和行为(组合权重0.013 3)的认知度也比昆山市其他区域差，体现出社区对居民的相关宣传有效率不足，需进一步落实社区对所辖居民的宣传教育，充分发挥社区在提升城市韧性方面的重要作用。
2) 该区域水面率(组合权重0.074 0)、河网密度(组合权重0.013 8)均较高，说明其蓝色基础设施的工程条件好，但存在大量断头浜，如果对其进行沟通改造，既能够降低发生洪涝灾害的风险，又能够通过水系的连通增加水面率，提高该圩区对内涝的调蓄能力。
3) 该区域人均医疗财政投入(组合权重0.003 8)较低，万人拥有卫生技术人员数(组合权重0.005 7)也比其他区域少。因此，适当增加医疗财政投入、增加卫生技术人员，有助于提高该圩区的恢复性。
由此可见，通过指标体系的评估结果和指标间的关联性，可以有效地实现优势指标和劣势指标之间的联动，使两个或多个韧性相关指标相互影响、相互作用、同步提升，从而有效地提高区域的韧性。
4 结论本文以韧性城市“适应性、应对性、恢复性”三大核心特征提升作为目标，构建了应对内涝灾害的多层次指标体系，分析了各级指标与韧性城市三大特征的一对多映射关系，并通过韧性概念图的构建，实现了该映射关系的直观展现。
基于GIS技术和AHP方法建立了韧性评估模型，并以江苏省昆山市为例进行了评估模型的应用。综合考虑韧性评估结果，可以高效定位城市应对内涝灾害的薄弱区域和环节，结合各四级指标的组合权重(对城市韧性的贡献度)，可以确定各指标的韧性提升潜力和效率。基于韧性评估模型的分析结果，结合现有的城市规划，考虑指标的互利互补，可实现针对特定区域的非常规韧性提升策略的挖掘与设计。
本研究虽然涉及指标众多，考虑因素较为全面，还针对具体城市进行了实地应用，但是仍有待完善优化之处，具体如下：
1) 指标处理。大部分指标是以行政区划作为数据搜集的单元，而在后续的评估中，不能够详细地呈现内涝灾害的空间分布规律，模糊了同一个行政区内不同空间的韧性评估结果，后续应该通过更为细致的数据搜集工作实现数据分布单元的细化。
2) 韧性评估。应用层次分析法进行韧性评估往往基于专家的个人主观认知，因此指标权重的计算存在一定的主观性，后续应该重点考虑如何削减指标权重确定的主观性，使韧性评估更加客观。
3) 提升策略。由于难以量化的指标多以行政区为单位(如社区减灾、居民减灾等)，而可量化的指标则以栅格为单位，因而在评估过程中对同一个行政区内不同空间位置的韧性进行了综合计算，使其失去了空间上的韧性差异。这也导致了城市韧性提升在空间、力度及体量上的精准施策存在困难。
4) 适用范围。本研究的评估只是对一个城市或区域内的相对韧性的评估，对于全国跨城市的评估和比较，还需要考虑城市间的差异，进行模型调整和指标补充。

参考文献

[1]	陈华丽, 陈刚, 丁国平. 基于GIS的区域洪水灾害风险评价[J]. 人民长江, 2003, 34(6): 49-51. CHEN H L, CHEN G, DING G P. Risk evaluation of regional flood disaster based on GIS[J]. Yangtze River, 2003, 34(6): 49-51. DOI:10.3969/j.issn.1001-4179.2003.06.018 (in Chinese)
[2]	MOTESHARREI S, RIVAS J, KALNAY E, et al. Modeling sustainability: Population, inequality, consumption, and bidirectional coupling of the earth and human systems[J]. National Science Review, 2016, 3(4): 470-494.
[3]	The Rockefeller Foundation. City resilience framework: The Rockefeller Foundation[R/OL]. [2021-03-15]. https://www.rockefellerfoundation.org/report/city-resilience-framework/.
[4]	黄弘, 李瑞奇, 于富才, 等. 安全韧性城市构建的若干问题探讨[J]. 武汉理工大学学报(信息与管理工程版), 2020, 42(2): 93-97. HUANG H, LI R Q, YU F C, et al. Discussion on several issues in the construction of safety and resilience city[J]. Journal of Wuhan University of Technology (Information & Management Engineering), 2020, 42(2): 93-97. (in Chinese)
[5]	李瑞奇, 黄弘, 周睿. 基于韧性曲线的城市安全韧性建模[J]. 清华大学学报(自然科学版), 2020, 60(1): 1-8. LI R Q, HUANG H, ZHOU R. Resilience curve modelling of urban safety resilience[J]. Journal of Tsinghua University (Science and Technology), 2020, 60(1): 1-8. (in Chinese)
[6]	陈长坤, 陈以琴, 施波, 等. 雨洪灾害情境下城市韧性评估模型[J]. 中国安全科学学报, 2018, 28(4): 1-6. CHEN C K, CHEN Y Q, SHI B, et al. A model for evaluating urban resilience to rainstorm flood disasters[J]. China Safety Science Journal, 2018, 28(4): 1-6. (in Chinese)
[7]	郑艳, 翟建青, 武占云, 等. 基于适应性周期的韧性城市分类评价: 以我国海绵城市与气候适应型城市试点为例[J]. 中国人口·资源与环境, 2018, 28(3): 31-38. ZHENG Y, ZHAI J Q, WU Z Y, et al. A typology analysis on resilient cities based on adaptive cycle: Taking cases of Chinese sponge cities and climate resilient cites pilot projects[J]. China Population, Resources and Environment, 2018, 28(3): 31-38. (in Chinese)
[8]	程元. 适应暴雨气候的城市韧性评价及规划路径探索: 以陕西省为例[D]. 西安: 西北大学, 2019. CHENG Y. Urban resilience assessment and planning path exploration to adapt to rainstorm climate: By the example of Shaanxi[D]. Xi'an: Northwestern University, 2019. (in Chinese)
[9]	BALICA S F, WRIGHT N G, VAN DER MEULEN F. A flood vulnerability index for coastal cities and its use in assessing climate change impacts[J]. Natural Hazards, 2012, 64(1): 73-105. DOI:10.1007/s11069-012-0234-1
[10]	ROSENZWEIG C, SOLECKI W. Hurricane Sandy and adaptation pathways in New York: Lessons from a first-responder city[J]. Global Environmental Change, 2014, 28: 395-408. DOI:10.1016/j.gloenvcha.2014.05.003
[11]	MUGUME S N, GOMEZ D E, FU G T, et al. A global analysis approach for investigating structural resilience in urban drainage systems[J]. Water Research, 2015, 81: 15-26. DOI:10.1016/j.watres.2015.05.030
[12]	ORENCIO P M, FUJII M. A localized disaster-resilience index to assess coastal communities based on an analytic hierarchy process (AHP)[J]. International Journal of Disaster Risk Reduction, 2013, 3: 62-75. DOI:10.1016/j.ijdrr.2012.11.006
[13]	QASIM S, QASIM M, SHRESTHA R P, et al. Community resilience to flood hazards in Khyber Pukhthunkhwa province of Pakistan[J]. International Journal of Disaster Risk Reduction, 2016, 18: 100-106. DOI:10.1016/j.ijdrr.2016.03.009
[14]	李彤玥. 韧性城市研究新进展[J]. 国际城市规划, 2017, 32(5): 15-25. LI T Y. New progress in study on resilient cities[J]. Urban Planning International, 2017, 32(5): 15-25. (in Chinese)
[15]	姜宇逍. 雨洪防涝视角下韧性社区评价体系及优化策略研究[D]. 天津: 天津大学, 2018. JIANG Y X. Studies on resilience community evaluation system and improvement strategies on an approach of flood disasters[D]. Tianjin: Tianjin University, 2018. (in Chinese)
[16]	邱桐, 陈湘生, 苏栋. 城市地下空间综合韧性防灾抗疫建设框架[J]. 清华大学学报(自然科学版), 2021, 61(2): 117-127. QIU T, CHEN X S, SU D. Comprehensive, resilient disaster and pandemic prevention construction framework for urban underground spaces[J]. Journal of Tsinghua University (Science and Technology), 2021, 61(2): 117-127. (in Chinese)
[17]	苏伯尼, 黄弘, 张楠. 基于情景模拟的城市内涝动态风险评估方法[J]. 清华大学学报(自然科学版), 2015, 55(6): 684-690. SU B N, HUANG H, ZHANG N. Dynamic urban waterlogging risk assessment method based on scenario simulations[J]. Journal of Tsinghua University (Science and Technology), 2015, 55(6): 684-690. (in Chinese)
[18]	杨莹, 林琳, 钟志平, 等. 基于应对公共健康危害的广州社区恢复力评价及空间分异[J]. 地理学报, 2019, 74(2): 266-284. YANG Y, LIN L, ZHONG Z P, et al. Comprehensive evaluation and spatial differentiation of community resilience in Guangzhou based on response of the city to public health hazards[J]. Acta Geographica Sinica, 2019, 74(2): 266-284. (in Chinese)
[19]	JOERIN J, SHAW R, TAKEUCHI Y, et al. The adoption of a climate disaster resilience index in Chennai, India[J]. Disasters, 2014, 38(3): 540-561. DOI:10.1111/disa.12058
[20]	张明斗, 冯晓青. 中国城市韧性度综合评价[J]. 城市问题, 2018(10): 27-36. ZHANG M D, FENG X Q. Comprehensive evaluation of urban resilience in China[J]. Urban Problems, 2018(10): 27-36. (in Chinese)
[21]	TOBIN G A, WHITEFORD L M. Community resilience and volcano hazard: The eruption of Tungurahua and evacuation of the Faldas in Ecuador[J]. Disasters, 2002, 26(1): 28-48. DOI:10.1111/1467-7717.00189
[22]	NORRIS F H, STEVENS S P, PFEFFERBAUM B, et al. Community resilience as a metaphor, theory, set of capacities, and strategy for disaster readiness[J]. American Journal of Community Psychology, 2008, 41(1-2): 127-150. DOI:10.1007/s10464-007-9156-6
[23]	CUTTER S L, BARNES L, BERRY M, et al. A place-based model for understanding community resilience to natural disasters[J]. Global Environmental Change, 2008, 18(4): 598-606. DOI:10.1016/j.gloenvcha.2008.07.013
[24]	白立敏, 修春亮, 冯兴华, 等. 中国城市韧性综合评估及其时空分异特征[J]. 世界地理研究, 2019, 28(6): 77-87. BAI L M, XIU C L, FENG X H, et al. A comprehensive assessment of urban resilience and its spatial differentiation in China[J]. World Regional Studies, 2019, 28(6): 77-87. DOI:10.3969/j.issn.1004-9479.2019.06.2018403 (in Chinese)
[25]	KULIG J C, EDGE D S, JOYCE B. Understanding community resiliency in rural communities through multimethod research[J]. Journal of Rural and Community Development, 2008, 3(3): 76-94.
[26]	HUDEC O, REGGIANI A, ?ISEROVá M. Resilience capacity and vulnerability: A joint analysis with reference to Slovak urban districts[J]. Cities, 2018, 73: 24-35. DOI:10.1016/j.cities.2017.10.004
[27]	黄弘, 李瑞奇, 范维澄, 等. 安全韧性城市特征分析及对雄安新区安全发展的启示[J]. 中国安全生产科学技术, 2018, 14(7): 5-11. HUANG H, LI R Q, FAN W C, et al. Analysis on characteristics of safety resilient city and enlightenments for safe development of Xiongan New Area[J]. Journal of Safety Science and Technology, 2018, 14(7): 5-11. (in Chinese)
[28]	李亚, 翟国方, 顾福妹. 城市基础设施韧性的定量评估方法研究综述[J]. 城市发展研究, 2016, 23(6): 113-122. LI Y, ZHAI G F, GU F M. Review on methods of quantification of urban infrastructure resilience[J]. Urban Development Studies, 2016, 23(6): 113-122. DOI:10.3969/j.issn.1006-3862.2016.06.016 (in Chinese)