基于二元Logistic对城市内涝风险区预测研究

陈韬; 钟传胤; 赵大维; 王前朋; 裘娜; Tao CHEN; Chuan-yin ZHONG; Da-wei ZHAO; Qian-peng WANG; Na QIU

中国农村水利水电 ›› 2022 ›› (5) : 57-61.

水环境与水生态

基于二元Logistic对城市内涝风险区预测研究

作者信息 +

Research on the Prediction of Urban Waterlogging Risk Area Based on Binary Logistic

Author information +

稿件信息 +

摘要

在气候变化和快速城市化背景下，城市内涝已经成为一种“城市病”。为了识别内涝风险区及预测研究，从城市空间因素的角度，以北京城六区作为研究区域，收集整理内涝信息和城市空间因素数据。采用二元Logistic回归分析，探究内涝风险区的特点，得出影响内涝的主要城市空间因素依次为汛期降雨量、综合径流系数、高速路网、雨水排放系统及涝点与桥梁的距离，并建立城六区内涝风险区预测模型，使用网络积水地图的积水点进行验证，结果表明预测模型结果准确，在现实中暴雨天气可以根据预测模型在高风险地区或是内涝发生高概率区域地点做好内涝防治措施。

Abstract

Under the background of climate change and rapid urbanization， urban waterlogging has become an “urban disease”. For identifying waterlogging risk area and prediction research， this paper explores the characteristics of waterlogging risk area， the main urban space factors affecting the waterlogging of flood season rainfall， runoff coefficient， highway system， drainage system and flood point distance and bridge from the perspectives of urban space factors， six areas as the study area， north of the capital collected waterlogging data information and the urban space factors， the binary logistic regression analysis， and waterlogging risk area prediction model is established by using the network map of water of water points for validation， results show that the prediction model is correct. In reality， rainstorm weather can be based on the forecast model in high-risk areas or areas with high probability of waterlogging prevention and control measures.

关键词

城市内涝 / 城市空间因素 / 内涝风险区 / 回归分析

Key words

urban waterlogging / urban space factor / waterlogging risk area / regression analysis

基金

北京市社会科学基金项目“北京市政给排水基础设施增量优化和管理机制研究”(18GLB013)

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用

陈韬 , 钟传胤 , 赵大维 , 王前朋 , 裘娜. 基于二元Logistic对城市内涝风险区预测研究[J].中国农村水利水电, 2022(5): 57-61

Tao CHEN , Chuan-yin ZHONG , Da-wei ZHAO , Qian-peng WANG , Na QIU. Research on the Prediction of Urban Waterlogging Risk Area Based on Binary Logistic[J].China Rural Water and Hydropower, 2022(5): 57-61

0 引言

随着城市的快速发展，近些年来城市气候条件也在不断变迁，这在一定程度上导致了城市内涝灾害的发生，并且已然成为各大城市的通病，应予以重视。十多年间，北京已开发区域面积拓展近500 km²，增长约40％。遭遇暴雨事件时，若城市的排水能力不足以将积水点内的雨水及时排出，便会发生内涝灾害，进而造成不必要的损失。其中下垫面硬化是导致城市积水点堆积的重要原因^［1，2］，城市中各种空间因素也与内涝的发生息息相关^［3］，交通网络迅速发展、桥梁分布不均、建筑物集中建造、人口不断增长、绿植面积被挤压等，致使部分马路在暴雨时极易“堆积成河”，对城市内居民的人身及财产安全造成严重威胁。北京地区从2012年开始，分三批开展城区内桥梁的修复优化工作，预计修补七十多座下凹式立交桥泵站，并在2016年整改完成，取得了良好的内涝防治效果。所以如何识别高风险内涝区，进行针对性地改造，具有重要意义。

国内外学者从不同角度开展城市内涝问题研究：Lowe等^［4］利用MIKEFLOOD在不同情景下探究降低内涝风险的方法；Lin等^［5］使用互联网公开数据基于POI风险评估方法对厦门城市内涝风险进行评估；宋晓猛等^［6］利用社会经济统计数据和水文气象资料，探讨了城市化背景下北京的洪涝特征、形成机制和影响因素；G. Tsakiris^［7］使用FLOW进行建模提出了洪涝评估的新模式，重点研究了城市地区；唐磊^［8］等人按照汇水区划分的方法分析北方城市内涝积水成因和优化方案；李芮等^［9］人采用InfoWorks ICM构建上清桥区域精细化洪涝模型，研究4种内涝减控措施的改善情况。

目前的研究主要以不同的影响因素对城市内涝进行了分析：Zhou ML等^［10］构建城市洪涝灾害风险评价模型探究自然灾害对内涝发生的影响；田子阳等^［11］提出了一种立体性的宏观、中观、微观、多尺度城市内涝风险评价技术；曹梦然等^［12］基于SWMM模型构建南京秦淮区城市排水管网，研究城市暴雨内涝风险的空间分布，但是对于城市哪些因素会对内涝产生影响，同时通过这些因素找出内涝的高风险区域的相关研究还较少。

1 研究方法与数据

1.1 研究区域和方法

研究区域选择在北京市中心城六区，即东城区、西城区、朝阳区、海淀区、丰台区和石景山区，如图1所示。区域内地势西北高、东南低，西北部为山地，东南部为平原，山地有东西走向的燕山和南北走向的太行山。区域内气候为典型的暖温带半湿润大陆性季风气候，夏季高温多雨，且全年的降水绝大部分来自汛期降雨，汛期雨强大，暴雨事件频发。研究区域内诱发内涝灾害的影响因素较多，且人口基数大、车辆往来密切，导致城区上方热量聚集，诱发“城市热岛效应”^［13］，产生局部暴雨事件。

图1 研究区域位置

Fig.1 Study area location

Full size|PPT slide

基于ArcGIS软件，首先对城市内各种空间因素进行处理加工，并利用二元Logistic回归分析方法探究城市内涝区域附近的状态，从而构建出内涝模型，对城六区内可能发生内涝灾害的位置进行预测，最后利用真实内涝点评估模型的准确度。并且可以依据预测结果，在城六区内部署内涝防治工作。

1.2 数据来源

内涝点信息是评价内涝影响因素的重点。获取方式主要是官方平台发布的内涝信息，比如建设部门公布的常年易内涝点，该途径的优点是能够准确把握内涝点的重要性，缺点是缺少内涝点的详细信息。其次可以利用网络检索，登录中新网搜集和整理近十年有关新闻报道，新闻报道内容较为详细，对内涝点的位置以及灾情信息提供准确，但缺点是内涝点信息覆盖不全、存在遗漏。综合以上统计数据，通过整合与筛选，最终得到218个城六区有效内涝点。对于其他城市空间因素，主要通过网络途径进行搜集，包括卫星数据以及一些软件资源，整合所需地图信息。

1.3 数据处理

根据城六区内涝点统计结果，具体位置分别位于东城区左安门桥下、东直门桥下、景泰桥下等；西城区西直门桥下、正阳桥下等；朝阳区东土城路北口、顾家庄桥北、朝阳区广顺南大街东口等；海淀区五路居桥下、颐和园东宫门、西苑桥下等；丰台区：京港澳高速南岗洼、菜户营桥下、玉泉营桥西等；石景山八大处路、金安桥下、石景山游乐园西门等。在ArcGIS中输入内涝点的基本信息，如位置、时刻等并将其数字化，展示出城六区内涝点的空间分布情况，并将其作为风险源，如图2（a）。

图2 城六区城市空间因素数据处理结果

Fig.2 Data processing results of urban spatial factors in the six districts of the city

Full size|PPT slide

高程是影响城市内涝的关键因素；当城市中某条道路的设计高程为相邻道路的最低点，道路两旁均为建设用地，此时周边地势均高于道路，如发生超标降雨，雨水管道不能及时排水，该位置将无疑会成为一处内涝积水点。搜集高程精确到五米的城区卫星地图，并将高程信息输入至ArcGIS中，如图2（b），探究其与内涝发生的相关性。

目前城市内具有多种土地类型，而不同类型间不透水面积与径流系数的大小有所差别，内涝的发生与这些数据信息密切相关^［14］。对北京地区分辨率为30 m的卫星地图处理分析，大致可以将地表覆盖种类分成湿地，人造覆盖（多数由混凝土或沥青等构成），种植地，森林等，如图2（c）。查找相应设计规范以及龚诗涵等人的研究^［15］，各种地表覆盖的径流系数分别取值为：湿地取1，人造覆盖取0.85，裸露地表取0.5，种植地取0.2，绿化用地取0.15，森林取0.1。将不同覆盖类型的径流系数值输入至ArcGIS中，通过加权计算，便可得到研究区域的综合径流系数。

有研究发现住宅用地与内涝发生有着密不可分的关系^［16］，区域内的人口数量反映着居住面积的大小。人口普查信息如表1所示，将人口密度信息细分至各乡镇街道，并在ArcGIS中进行数据提取，得到如图2（d）所示的人口分布图。

表1 城六区人口普查信息

Tab.1 Census information of the six districts of the city

行政区	人口数/人	人口密度/（人·km^-2）	面积/km²
西城	1 106 214	21 690	51
东城	708 829	16 877	42
朝阳	3 452 460	7 330	471
海淀	3 133 469	7 270	431
丰台	2 019 764	6 644	304
石景山	567 851	6 603	86

北京全年降雨量分布不均，其中汛期（第三季度）雨量约占七成以上。搜集城六区内气象站近8年汛期降雨量数据，整理后如表2所示，使用克里金插值方法计算出剩余地区的雨量数据，并在ArcGIS中对研究区域内汛期降雨量数值信息进行可视化，如图2（e）所示。

表2 气象站汛期降雨量数据 (mm)

Tab.2 Rainfall data of weather station during flood season

区站号	年份
区站号	2013	2014	2015	2016	2017	2018	2019	2020
北京54511	511.6	388.7	319.6	553.0	451.2	479.4	254.2	525.8
朝阳54433	376.1	383.8	489.1	565.9	528.4	569.6	435.1	739.9
丰台54514	432.0	343.1	520.2	547.9	526.2	354.1	331.9	467.3
海淀54399	486.7	413.8	483.7	583.2	524.2	463.6	354.5	501.3
石景山54513	474.2	251.0	566.5	555.1	545.8	457.2	319.5	494.5

依照北京市市政管网的走向与分布，将雨水管网视为线要素，基于核密度估计（KDE）方法，综合评估现有数据模式，构建平滑密度表面，以此体现出研究区域内的管线密度，并将线要素矢量化分析处理，得到如图2（f）所示的城市管网分布图，评判其与内涝发生的相关性。

许多研究^［17］探究了不透水下垫面与内涝发生的关系，不透水下垫面种类多样，而道桥系统不透水面积最大，各种道路对内涝发生的影响不同，高速路的排水沟通常为中央隔离带，排水管为横向布置^［18］；城市内道路采用市政雨水管网排水以及渗滞系统辅助作用，城市内桥梁下方区域内涝灾害频发，大多数桥梁连接着相邻道路，地表径流会沿不同方向从高速路等道路汇入桥梁下方，若此时超出桥梁下方排水负荷，便会在桥下产生积水^{［19，20］}。如图2（g）所示为道桥分布图，本文将交通系统进行细化分类，包括普通道路、高速路（包含高架）、桥梁，作为内涝发生的不同影响因素加以分析，制作道路密度分布图。根据欧氏距离法计算各像素点与桥梁的最短间距，分析内涝的发生是否还与桥梁坐标有关。

北京历史上曾多次因暴雨导致河流泛滥进而引起洪涝事件^［21-23］。利用卫星地图提取出北京市内的河流水系信息并输入至软件中，加入回归分析。包括潮白河、北运河、永定河、大清河、蓟运河五大水系，以及城六区的河流湖泊如八一湖、青年湖、金海湖、雁栖湖、北海、白龙潭等，如图2（h），分析河流水系是否影响着内涝灾害的发生。

2 北京城六区内涝影响因素分析

2.1 内涝点的空间分布特征

分析研究区域内的内涝点分布情况，可以看出其空间分布特征未表现出规律性，其中二环附近内涝点分布最少；随着环数的增加，内涝点数量逐渐增加，直到扩散至五环外时呈现衰减趋势。对统计出的共218个城六区有效内涝点，使用热点分析（HSA）Getis-Ord Gi*工具分析其与内涝的关联性，通过全部个体的要素均值来评析整体特性，查看临近要素环境中每个要素的自身相关度，将统计出的内涝点与随机生成点位输入至软件中进行热点分析模拟，得到城六区各处内涝概率，如图3所示，其中蓝色表明此处内涝风险概率较低，红色则相反。

图3 城六区热点分析

Fig.3 Six Area Hot Point Analysis

Full size|PPT slide

结果显示高风险点主要集中分布于东城区与朝阳区接壤处附近，高风险点还散布于丰台区与海淀区，而低风险点主要位于城市中心区域。

2.2 局部范围内的城市空间因素

为探究内涝形成原因，需综合内涝点附近全部城市空间因素，并生成概率预测模型^［24］。将内涝灾害辐射范围设为R=500 m，提取范围内全部空间因素。在研究区域内生成816个随机点，利用同样方法提取数据，探究各城市空间因素与内涝灾害的相关性。

利用二元Logistic回归分析方法，而影响内涝发生的自变量设为人口密度、地面高程、雨水排放系统、河流湖泊、高速路网、城区道路、涝点与桥梁的距离、综合径流系数和汛期降雨量。该方法基于最大似然估计的向前逐步回归法（LR），对于模型的比选具有指导意义，逐步迭代补充完整各项解释因子，直至拟合出最优模型结果，并筛除了无关紧要的解释因子。如表3所示为历经五次迭代处理后，得到的主要影响因子数据。结果显示似然比检验结果p=0（＜0.05），表示自变量取值合理，具有分析意义。其次通过霍斯默-莱梅肖检验（Hosmer and Lemeshow Test）p=0.529（＞0.05），可以认为综合利用了数据信息，模型应用具有可行性。

表3 二元Logistic回归结果（LR法）

Tab.3 Binary Logistic Regression Results

自变量	解释系数	标准误差	Wald统计	自由度	显著性	优势比OR
高速路网（a）	0.021	0.007	8.314	1	0.004	1.021
桥梁距离（b）	0.042	0.010	18.441	1	0	1.043
径流系数（c）	-0.004	0.002	4.780	1	0.029	0.996
汛期雨量（d）	0.005	0.001	72.036	1	0	1.005
雨水系统（e）	-0.875	0.181	23.358	1	0	0.417
常量	-14.230	3.621	15.447	1	0	0

根据结果得出模型公式：

l n (\frac{p}{1 - p}) = - 14.228 + 0.005 a - 0.875 b + 0.042 c + 0.021 d - 0.004 e

式中：p为出现内涝概率；1-p为未出现内涝概率，字母表示各种解释因子。

从公式中的系数情况可以看出，呈现正相关的影响因素为汛期降雨量、地表径流系数、高速路网，呈负相关的为雨水排放系统和涝点与桥梁的距离。并且内涝灾害在汛期降雨量大、地表径流系数高的城区发生概率更高；而若城区内雨水排水管网较为完备，或附近桥梁数量较少，则不易发生洪涝灾害。

3 预测模型验证

根据上述公式计算出各区域内涝发生的概率，反映内涝的倾向，并将内涝分布情况在ArcGIS中可视化，如图4所示，各颜色所代表的内涝概率如图例所示，其中红色表示发生概率较高。

图4 内涝概率预测与验证

Fig.4 Prediction and verification of waterlogging probability

Full size|PPT slide

搜集网络信息，统计出2020年内共96个现实内涝点，并将内涝点位置坐标输入至软件中可视化处理，与内涝模型图进行对比，对其预测有效性予以实际验证，拟合情况如表4所示。

表4 现实积水点在模型中的风险分布

Tab.4 Realistic waterlogging point in the risk area in the model

内涝风险概率	内涝点数量
高度风险区（p≥90%）	31
中高度风险区（90%>p≥70%）	26
中度风险区（70%>p≥50%）	24
低度风险区（p<50%）	15

从表4中可以看出，内涝灾害发生概率在50%以上的积水点共有81个，包含内涝灾害发生概率在90%以上的积水点共31个。将内涝概率在50%以上的区域视为易发生内涝处，则预测模型具有84%的准确度，说明模型构建较为合理，内涝风险预测模型在现实生活中具有实际应用意义，能够指导防范内涝风险。

4 结语

（1）本文首先在内涝点附近范围内进行数据分析，使用热点分析方法探究其空间分布特征。其次利用二元Logistic回归分析方法逐步迭代计算，得到与内涝发生呈正负相关性的各种城市空间因素，据此构建城六区内涝预测模型，最后利用实际内涝点分析模型的有效性。结果显示预测模型的准确率很高，模型具有现实应用价值。

（2）不透水下垫面被认为是内涝灾害的罪魁祸首，在交通系统中其具有较大的面积。为了细致研究，本文将系统分为三类路面进行具体分析，得出其中的高速路网和城区桥梁与内涝灾害的关系最为密切。并且结果也表明城区桥梁与径流系数为内涝的主要影响因素。

（3）在未来分析影响内涝发生的城市空间因素时，可对其进一步细化分类研究，如用汛期内典型场次降雨数据替代平均雨量、考虑雨水系统中管长管径等因素、不同桥梁下方排水能力的大小、更为详尽的地区地形特点、城市建筑物参数信息等，分析其与内涝的相关性，不断优化城市内涝预测模型，对于内涝预测概率图中高度风险区域，可在未来防治工作加以重点布置，排除内涝隐患，切实保障人身以及财产安全。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序

1	YANG Y， LIU Y， LI Y， et al. Quantifying spatio-temporal patterns of urban expansion in Beijing during 1985–2013 with rural-urban development transformation［J］. Land Use Policy， 2018，74：220-230. 本文引用 [1]

2	CHE W， YANG Z， ZHAO Y， et al. Analysis of urban flooding control and major and minor drainage systems in China［J］. China Water Wastewater， 2013，29：13-19. 本文引用 [1]

3	SAJIKUMAR N， REMYA R. Impact of land cover and land use change on runoff characteristics［J］. Journal of Environmental Management， 2015，161：460-468. 本文引用 [1]

4	LÖWE R， URICH C， DOMING N S， et al. Assessment of urban pluvial flood risk and efficiency of adaptation options through simulations–A new generation of urban planning tools［J］. Journal of Hydrology， 2017，550：355-367. 本文引用 [1]

5	LIN T， LIU X， SONG J， et al. Urban waterlogging risk assessment based on internet open data： A case study in China［J］. Habitat International， 2018，71：88-96.. 本文引用 [1]

6	宋晓猛，张建云，贺瑞敏，等. 北京城市洪涝问题与成因分析［J］. 水科学进展， 2019，30（2）：153-165. 本文引用 [1]

7	TSAKIRIS G. Flood risk assessment： concepts， modelling， applications［J］. Natural Hazards Earth System Sciences， 2014，14（5）：1 361. 本文引用 [1]

8	唐磊，周飞祥，王巍巍，等. 北方城市典型内涝积水问题的系统化解决方案［J］. 中国给水排水， 2020，36（13）：139-144. 本文引用 [1]

9	李芮，潘兴瑶，邸苏闯，等. 城市内涝应对措施改善效果分析：以清河上清桥为例［J］. 给水排水， 2019，55（）：128-135. 本文引用 [1] 摘要增刊1

10	ZHOU M L. An alternative risk assessment model of urban waterlogging： a case study of Ningbo City［J］. Sustainability， 2021，13（2）：826-826. 本文引用 [1]

11	田子阳，褚俊英，林永寿，等.多尺度城市内涝风险评价技术及应用［J/OL］.水资源保护：1-15［2021-09-18］. https：//kns-cnki-net.door.bucea.edu.cn/kcms/detail/32.1356.TV.20210910.1312.002. html 本文引用 [1]

12	曹梦然，叶亚平，张其成，等.基于SWMM模型的南京秦淮区暴雨内涝风险分析［J］.中国农村水利水电，2019（11）：110-114. 本文引用 [1]

13	YANG J， WANG Y， XIU C， et al. Optimizing local climate zones to mitigate urban heat island effect in human settlements［J］. Journal of Cleaner Production， 2020，275. 本文引用 [1]

14	LINDBERG J. Locating potential flood areas in an urban environment using remote sensing and GIS， case study Lund， Sweden［J］. Student thesis series INES， 2015. 本文引用 [1]

15	龚诗涵，肖洋，方瑜，等. 中国森林生态系统地表径流调节特征［J］. 生态学报， 2016，36（22）：7 472-7 478. 本文引用 [1]

16	吴健生，张朴华. 城市景观格局对城市内涝的影响研究：以深圳市为例［J］. 地理学报， 2017，72（3）：444-456. 本文引用 [1]

17	YU H， ZHAO Y， FU Y， et al. Spatiotemporal variance assessment of urban rainstorm waterlogging affected by impervious surface expansion： a case study of Guangzhou， China［J］. Sustainability， 2018，10（10）：3 761. 本文引用 [1]

18	刘广. 环城高速路快速施工软基处治及施工关键技术［D］. 西安：长安大学， 2019. 本文引用 [1]

19	叶超凡，张一驰，程维明，等. 北京市区快速城市化进程中的内涝现状及成因分析［J］. 中国防汛抗旱， 2018，28（2）：19-25. 本文引用 [1]

20	刘兴哲. 城市下凹桥区内涝积水综合治理设计探讨［J］. 中国给水排水， 2019， 35（4）： 47-53. 本文引用 [1]

21	杨佩国，靳京，赵东升，等. 基于历史暴雨洪涝灾情数据的城市脆弱性定量研究：以北京市为例［J］. 地理科学， 2016，36（5）：733-741. 本文引用 [1]

22	叶瑜，徐雨帆，梁珂，等. 1801年永定河水灾救灾响应复原与分析［J］. 中国历史地理论丛， 2014，29（4）：13-19.

23	段天顺. 谈谈北京历史上的水患［J］. 中国水利， 1982（3）：16-17. 本文引用 [1]

24	赵大维. 北京市政给排水基础设施增量优化研究［D］. 北京：北京建筑大学，2021. 本文引用 [1]

PDF(3751 KB)

521

访问

引用

详细情况

段落导航

收稿日期	出版日期
2021-08-18	2022-05-15
发布日期
2022-06-28

选择文件类型/文献管理软件名称

选择包含的内容