CN111507375A - 一种城市内涝风险快速评估方法及系统 - Google Patents

Abstract

本发明公开了一种城市内涝风险快速评估方法及系统，以内涝发生快慢为衡量风险的标准：对于给定内涝点，降雨之后达到危险积水深度所需的时间越短，留给应急响应的时间越有限，则相应的风险越大。以上方案只需要设计暴雨、DEM、土地利用/覆盖和排水系统设计标准即可计算，不涉及复杂的水动力模型，对排水管网数据要求低，基于GIS平台即可完成计算。本发明方法及系统解决了现有内涝风险评估对基础数据及操作人员的建模能力要求高、且其存在计算效率低、实用性差的问题。

Description

一种城市内涝风险快速评估方法及系统

技术领域

本发明涉及内涝灾害风险评估技术领域，尤其涉及一种城市内涝风险快速评估方法及系统。

背景技术

在全球气候变暖和城市化背景下，城市地区的“雨岛效应”、产流增加和汇流加快等水文效应的综合作用使得整体上城市内涝的压力与风险呈增加趋势，需要通过海绵城市建设来减缓内涝风险，其减缓效果则取决于对实际内涝风险的精准理解与评估。已有内涝风险评估多采用基于“降雨-产汇流-管道排水-地表溢流”水动力内涝模拟，这类方法对基础数据要求高，对操作人员的建模能力要求高，且计算效率低，因而在实际应用中难以推广使用。

发明内容

本发明为解决现有的内涝风险评估对基础数据及操作人员的建模能力要求高、且其存在计算效率低、实用性差的问题，提供了一种城市内涝风险快速评估方法及系统。

为实现以上发明目的，而采用的技术手段是：

一种城市内涝风险快速评估方法，包括以下步骤：

S1.对于目标区域，根据DEM识别潜在的内涝发生位置，划定每一个内涝发生位置的汇水范围，并构成相应的小流域；

S2.计算每个所述小流域的属性，包括其汇水面积A、平均坡度S、积水体积V；

S3.设定降雨情景，根据水量平衡模型，考虑包括排水系统设计标准、土壤渗透和蒸散发的影响，计算在地表形成积水的内涝雨量R_i；

S4.以所述小流域为单元，以内涝发生快慢作为衡量标准，计算目标区域的内涝风险Risk：Risk＝ln(Rⁱ×A×S^1/2/V)。

上述方案中，以内涝发生快慢为衡量标准的风险指标：对于给定内涝点，降雨之后其达到危险积水深度所需的时间越短，留给应急响应的时间越有限，则相应的风险越大。以上方案只需要设计暴雨、DEM、土地利用/覆盖和排水系统设计标准即可计算，计算过程不涉及复杂的水动力模型，对排水管网数据要求低，基于GIS平台即可完成计算。

优选的，所述步骤S1包括以下步骤：

S11.洼地识别：对原始DEM(数字高程模型)进行填洼，计算填洼后DEM与原始DEM的差值DEM；在所述差值DEM中进行搜索，一组空间上相邻的非零值栅格即为一个洼地；所述空间相邻基于八邻域定义，即在3×3的栅格范围内，中心栅格与周边8个栅格是相邻的；

S12.阈值分析：对于识别得到的所有洼地，按面积从小到大排序，确定预设的面积阈值；

S13.洼地填充：对于所述面积小于面积阈值的洼地，用填洼后DEM的高程值表示，仅保留面积大于或等于面积阈值的目标洼地；

S14.流域搜索：对于每个目标洼地，采用种子填充算法搜索确定其汇水范围，所述汇水范围与洼地自身构成一个小流域。

优选的，所述步骤S12中，所述面积阈值的大小设置应为：使得所保留的面积大于或等于面积阈值的目标洼地，占识别得到的所有洼地总面积中的比例不小于80％。在本优选方案中，筛选洼地需要在减少洼地数量和保持地形代表性之间达到平衡。洼地数量随其面积的增加而急剧减少，呈帕累托分布。通过设定适当的面积阈值，去除小于该阈值的洼地，并使得保留洼地占洼地总面积中的比例不小于80％，这样既能表达主要地形特征，又能显著减少分析单元的数量。

优选的，步骤S14中所述的采用种子填充算法搜索确定其汇水范围具体为：

以洼地的栅格为种子，按八邻域进行搜索：对于一个种子，将其邻域中高程不小于种子的非洼地栅格标记为汇水栅格，并加入种子集合；不断执行此操作，直至找不到新的汇水栅格；得到的种子集合即为该洼地的汇水范围。

优选的，步骤S1所述的DEM采用由激光雷达技术生产的高精度高分辨率DEM。在本优选方案中，由于该步骤都是基于DEM进行的，DEM空间分辨率和高程精度对于内涝风险评估结果的合理性至关重要，在具备条件时应尽可能使用由激光雷达技术生产的高精度高分辨率DEM。选择DEM空间分辨率时，需要在微地形表达精度(分辨率越高，则精度越高)和计算效率(分辨率越低，则效率越高)之间达到一定的平衡。

优选的，所述步骤S2包括以下步骤：

S21.所述小流域的汇水面积A、平均坡度S计算：所述小流域的汇水面积A为该小流域对应洼地的自身面积与其汇水范围面积之和；平均坡度S采用洼地汇水范围内栅格的平均坡度值；

S22.所述小流域的积水体积V计算：积水体积V为洼地达到危险积水深度所需的水量；先设定积水深度，得到该积水深度对应的水面高程H，统计洼地低于水面高程H的栅格共有m个，且对应的高差为h_i(i＝1,2,...,m)；设定DEM的栅格大小为d；则积水体积V为：

优选的，所述步骤S3包括以下步骤：

S31.计算降雨量R：

设定一个或多个降雨情景，对于每个降雨情景，执行以下计算：

式中q为设定暴雨强度，t为设定降雨历时，P为设定重现期；其中A₁，C，b，n根据统计方法进行确定；计算得到该降雨情景暴雨强度q后，乘以降雨历时即得到其降雨量：R＝q×t；

S32.计算产流雨量R_f：

当目标区域具备实测的蒸散发量和土壤渗透能力数据时：利用已有的不透水面成果，或者采用监督分类方法从高分遥感影像或航空正射影像中提取不透水面的空间分布，并结合渗透模型计算通过土壤的下渗量；降雨量R扣除实测的蒸散发量和计算下渗量后即得到洼地小流域的产流雨量R_f；

当目标区域不具备实测的蒸散发量和土壤渗透能力数据时：根据洼地流域的土地利用现状确定不同土地利用类型的面积比例，再使用《室外排水设计规范》给出的径流系数推荐值，采用面积加权平均的方式计算每个洼地小流域的综合径流系数；降雨量R乘以综合径流系数即得到洼地小流域的产流雨量R_f；

S33.计算管道排水P：管道排水P根据目标区域的城市排水系统的设计资料确定或采用水力模型对目标区域的实际管道排水能力进行校核或根据目标区域的实际降雨事件对应的内涝情况进行经验性的反推；

S34.计算内涝雨量R_i：R_i＝R_f-P，内涝雨量为产流雨量R_f与管道排水P之差；内涝雨量R_i值应为正值，若为非正值则说明无内涝发生。

优选的，步骤S32所述的渗透模型为Horton下渗模型或Green-Ampt下渗模型。

本发明还提供了一种城市内涝风险快速评估系统，包括：

小流域构建模块，用于根据DEM识别目标区域潜在的内涝发生位置，划定每一个内涝发生位置的汇水范围，并构成相应的小流域；

小流域属性计算模块，用于计算每个所述小流域的属性，包括其汇水面积A、平均坡度S、积水体积V；

内涝雨量计算模块，用于设定降雨情景，根据水量平衡模型，考虑包括排水系统设计标准、土壤渗透和蒸散发的影响，计算在地表形成积水的内涝雨量R_i；

内涝风险评估模块，用于以所述小流域构建模块得到的小流域为单元，以内涝发生快慢作为衡量标准，根据所述小流域属性计算模块及内涝雨量计算模块的计算结果计算目标区域的内涝风险Risk：Risk＝ln(R_i×A×S^1/2/V)。

与现有技术相比，本发明技术方案的有益效果是：

本发明的城市内涝风险快速评估方法及系统以内涝发生快慢为衡量风险的标准：对于给定内涝点，降雨之后其达到危险积水深度所需的时间越短，留给应急响应的时间越有限，则相应的风险越大。以上方案只需要设计暴雨、DEM、土地利用/覆盖和排水系统设计标准即可计算，计算过程不涉及复杂的水动力模型，对排水管网数据要求低，基于GIS平台即可完成计算。本发明方法及系统解决了现有内涝风险评估对基础数据及操作人员的建模能力要求高、且其存在计算效率低、实用性差的问题。

附图说明

图1为本发明内涝风险计算的原理示意图。

图2为本发明方法的流程图。

图3为本发明中当目标区域具备实测的蒸散发量和土壤渗透能力数据时内涝雨量的计算流程图。

图4为本发明中当目标区域不具备实测的蒸散发量和土壤渗透能力数据时内涝雨量的计算流程图。

图5为实施例2中对洼地进行流域搜索后的结果示意图。

图6为实施例2中对目标区域进行内涝分析得到的结果示意图。

图7为本发明系统的模块示意图。

具体实施方式

附图仅用于示例性说明，不能理解为对本专利的限制；

为了更好说明本实施例，附图某些部件会有省略、放大或缩小，并不代表实际产品的尺寸；

对于本领域技术人员来说，附图中某些公知结构及其说明可能省略是可以理解的。

下面结合附图和实施例对本发明的技术方案做进一步的说明。

实施例1

本实施例提供了一种城市内涝风险快速评估方法，如图1和2所示，包括以下步骤：

S1.对于目标区域，根据DEM识别潜在的内涝发生位置，划定每一个内涝发生位置的汇水范围，并构成相应的小流域；由于该步骤是基于DEM进行的，DEM空间分辨率和高程精度对于内涝风险评估结果的合理性至关重要，在具备条件时应尽可能使用由激光雷达技术生产的高精度高分辨率DEM。选择DEM空间分辨率时，需要在微地形表达精度(分辨率越高，则精度越高)和计算效率(分辨率越低，则效率越高)之间达到一定的平衡。

该步骤S1的具体步骤如下：

S11.洼地识别：对原始DEM进行填洼，计算填洼后DEM与原始DEM的差值DEM；在所述差值DEM中进行搜索，一组空间上相邻的非零值栅格即为一个洼地；所述空间相邻基于八邻域定义，即在3×3的栅格范围内，中心栅格与周边8个栅格是相邻的；

S12.阈值分析：对于识别得到的所有洼地，按面积从小到大排序，确定预设的面积阈值；

由于筛选洼地需要在减少洼地数量和保持地形代表性之间达到平衡。洼地数量随其面积的增加而急剧减少，呈帕累托分布。通过设定适当的面积阈值，使得保留洼地占洼地总面积中的比例不小于80％，这样既能表达主要地形特征，又能显著减少分析单元的数量。同时，选择洼地面积阈值时，需要考虑洼地小流域空间尺度对汇流时间的影响，即小流域面积有上限。当小流域面积过大时，由于地表径流耗时和地表调蓄作用的影响，部分汇水面积并不会对目标洼地的积水体积有直接贡献。

S13.洼地填充：对于所述面积小于面积阈值的洼地，用填洼后DEM的高程值表示，仅保留面积大于或等于面积阈值的目标洼地；

S14.流域搜索：对于每个目标洼地，采用种子填充算法搜索确定其汇水范围，所述汇水范围与洼地自身构成一个小流域。其中采用种子填充算法搜索确定其汇水范围具体为：

以洼地的栅格为种子，按八邻域进行搜索：对于一个种子，将其邻域中高程不小于种子的非洼地栅格标记为汇水栅格，并加入种子集合；不断执行此操作，直至找不到新的汇水栅格；得到的种子集合即为该洼地的汇水范围。

S2.计算每个所述小流域的属性，包括其汇水面积A、平均坡度S、积水体积V；具体步骤包括：

S21.所述小流域的汇水面积A、平均坡度S计算：所述小流域的汇水面积A为该小流域对应洼地的自身面积与其汇水范围面积之和；平均坡度S采用洼地汇水范围内栅格的平均坡度值；

S22.所述小流域的积水体积V计算：积水体积V为洼地达到危险积水深度所需的水量；先设定积水深度，在本实施例中该积水深度设为50厘米此深度积水会阻断交通。实际应用本方法时，该积水深度亦可根据特定分析目标由本领域技术人员另行设定。确定积水深度后，可得到该积水深度对应的水面高程H(米)，统计洼地低于水面高程H(米)的栅格共有m个，且对应的高差为h_i(i＝1,2,...,m)(米)；设定DEM的栅格大小为d(米)；则积水体积V(立方米)为：

S3.设定降雨情景，根据水量平衡模型，考虑包括排水系统设计标准、土壤渗透和蒸散发的影响，计算在地表形成积水的内涝雨量R_i；具体包括以下步骤：

S31.计算降雨量R：

设定一个或多个降雨情景，其设定基于目标区域的排水现状，包括目标区域的排水系统的设计标准、实际内涝严重程度以及预期应对的暴雨重现期；对于每个降雨情景，执行以下计算：

式中q为设定暴雨强度(升/(秒·公顷))，t为设定降雨历时(分钟)，P为设定重现期(年)；其中A₁，C，b，n根据统计方法进行确定；计算得到该降雨情景暴雨强度q后，乘以降雨历时即得到其降雨量：R＝q×t(毫米)；

S32.计算产流雨量R_f：

如图3所示，当目标区域具备实测的蒸散发量和土壤渗透能力数据时：利用已有的不透水面成果，或者采用监督分类方法从高分遥感影像或航空正射影像中提取不透水面的空间分布，并结合Horton渗透模型计算通过土壤的下渗量；降雨量R扣除实测的蒸散发量和计算下渗量后即得到洼地小流域的产流雨量R_f；

如图4所示，当目标区域不具备实测的蒸散发量和土壤渗透能力数据时：根据洼地流域的土地利用现状确定不同土地利用类型的面积比例，再使用《室外排水设计规范》给出的径流系数推荐值，采用面积加权平均的方式计算每个洼地小流域的综合径流系数；降雨量R乘以综合径流系数即得到洼地小流域的产流雨量R_f；

S33.计算管道排水P：管道排水P根据目标区域的城市排水系统的设计资料确定(如设定重现期为1年，则认为其管道排水即为设定重现期和历时对应的产流量)；必要时亦可采用水力模型对目标区域的实际管道排水能力进行校核或根据目标区域的实际降雨事件对应的内涝情况进行经验性的反推来确定管道排水P；

S34.计算内涝雨量R_i：R_i＝R_f-P，内涝雨量为产流雨量R_f与管道排水P之差；内涝雨量R_i值应为正值，若为非正值则说明无内涝发生。

S4.以所述小流域为单元，以内涝发生快慢作为衡量标准，计算目标区域的内涝风险Risk：Risk＝ln(R_i×A×S^1/2/V)。

本实施例所提供的城市内涝风险快速评估方法，以内涝发生快慢为衡量标准的风险指标：对于给定内涝点，降雨之后其达到危险积水深度所需的时间越短，留给应急响应的时间越有限，则相应的风险越大。以上方案只需要DEM、土地利用/覆盖和排水系统设计标准即可计算，且计算过程不涉及复杂的水动力模型，对排水管网数据要求低，可基于GIS平台即完成计算。

实施例2

为验证上述实施例1城市内涝风险快速评估方法的有效性，进行了如下内涝风险的实例分析：

该实例分析基于广州市某区域，面积约20平方公里，不透水率为73％，设计排水能力可应对1年重现期暴雨。本实施例分析5年重现期降雨情景下该区域的内涝风险，具体步骤和结果说明如下：

(1)采用机载激光雷达技术生产1米分辨率DEM，经洼地识别和阈值分析后设定本例的面积阈值为5000平方米，经筛选后保留有53个洼地；

(2)采用流域搜索方法确定每一个洼地的汇水范围，该汇水范围与洼地自身构成一个小流域。结果如图5所示，图中黑色表示洼地，其他灰度表示洼地对应汇水范围。进一步计算小流域的汇水面积A、平均坡度S和积水体积V，其中计算积水体积V时将积水深度设定为50厘米；

(3)模拟5年重现期、1小时历时的暴雨情景，根据广州市暴雨公式可计算其降雨量为75.5毫米，相应的1年重现期降雨量为54.5毫米；

(4)管道排水即为设计排水能力对应的产流量，根据1年重现期降雨量和综合径流系数计算。本实施例中将该区域用于排水管道设计的综合径流系数统一设为0.5，即管道排水P＝54.5×0.5＝27.25毫米；

(5)采用监督分类方法从高精度航空正射影像提取透水面与不透水面，二者的径流系数分别设为0.3和0.9，基于面积加权计算每一个洼地流域的综合径流系数；

(6)对每一个洼地流域，根据5年重现期降雨量和实际综合径流系数计算产流雨量R_f，然后从产流雨量R_f中扣除管道排水27.25毫米，即得到内涝雨量R_i。若内涝雨量R_i非正值，则认为该洼地无内涝风险。

(7)对于内涝雨量R_i大于零的洼地，根据以下公式计算其在5年重现期、1小时历时暴雨下的内涝风险：

Risk＝ln(R_i×A×S^1/2/V)。

根据自然断点法将计算得到的内涝风险分为低、中、高三类，其空间分布如图6所示，其中A、B、C为典型的内涝黑点。这说明该方法能准确定位内涝发生位置，并评价其风险等级，结果与历史内涝情况一致。

实施例3

本实施例提供了一种城市内涝风险快速评估系统，如图7所示，包括：

小流域构建模块1，用于根据DEM识别目标区域潜在的内涝发生位置，划定每一个内涝发生位置的汇水范围，并构成相应的小流域；

小流域属性计算模块2，用于计算每个所述小流域的属性，包括其汇水面积A、平均坡度S、积水体积V；

内涝雨量计算模块3，用于设定降雨情景，根据水量平衡模型，考虑包括排水系统设计标准、土壤渗透和蒸散发的影响，计算在地表形成积水的内涝雨量R_i；

内涝风险评估模块4，用于以所述小流域构建模块2得到的小流域为单元，根据小流域属性计算模块2及内涝风险评估模块4的计算结果，以内涝发生快慢作为衡量标准，计算目标区域的内涝风险Risk：Risk＝ln(R_i×A×S^1/2/V)。

该城市内涝风险快速评估系统基于实施例1的方法构建，其原理与上述方法一致，此处不再进行赘述。

附图中描述位置关系的用语仅用于示例性说明，不能理解为对本专利的限制；

显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。

Claims (9)

1.一种城市内涝风险快速评估方法，其特征在于，包括以下步骤：

S1.对于目标区域，根据DEM识别潜在的内涝发生位置，划定每一个内涝发生位置的汇水范围，并构成相应的小流域；

S2.计算每个所述小流域的属性，包括其汇水面积A、平均坡度S、积水体积V；

S3.设定降雨情景，根据水量平衡模型，考虑包括排水系统设计标准、土壤渗透和蒸散发的影响，计算在地表形成积水的内涝雨量R_i；

S4.以所述小流域为单元，以内涝发生快慢作为衡量标准，计算目标区域的内涝风险Risk：Risk＝ln(R_i×A×S^1/2/V)。

2.根据权利要求1所述的城市内涝风险快速评估方法，其特征在于，所述步骤S1包括以下步骤：

S11.洼地识别：对原始DEM进行填洼，计算填洼后DEM与原始DEM的差值DEM；在所述差值DEM中进行搜索，一组空间上相邻的非零值栅格即为一个洼地；所述空间相邻基于八邻域定义；

S12.阈值分析：对于识别得到的所有洼地，按面积从小到大排序，确定预设的面积阈值；

S13.洼地填充：对于所述面积小于面积阈值的洼地，用填洼后DEM的高程值表示，仅保留面积大于或等于面积阈值的目标洼地；

S14.流域搜索：对于每个目标洼地，采用种子填充算法搜索确定其汇水范围，所述汇水范围与洼地自身构成一个小流域。

3.根据权利要求2所述的城市内涝风险快速评估方法，其特征在于，所述步骤S12中，所述面积阈值的大小设置应为：使得所保留的面积大于或等于面积阈值的目标洼地，占识别得到的所有洼地总面积中的比例不小于80％。

4.根据权利要求2所述的城市内涝风险快速评估方法，其特征在于，步骤S14中所述的采用种子填充算法搜索确定其汇水范围具体为：

以洼地的栅格为种子，按八邻域进行搜索：对于一个种子，将其邻域中高程不小于种子的非洼地栅格标记为汇水栅格，并加入种子集合；不断执行此操作，直至找不到新的汇水栅格；得到的种子集合即为该洼地的汇水范围。

5.根据权利要求1所述的城市内涝风险快速评估方法，其特征在于，步骤S1所述的DEM采用由激光雷达技术生产的高精度高分辨率DEM。

6.根据权利要求2所述的城市内涝风险快速评估方法，其特征在于，所述步骤S2包括以下步骤：

S21.所述小流域的汇水面积A、平均坡度S计算：所述小流域的汇水面积A为该小流域对应洼地的自身面积与其汇水范围面积之和；平均坡度S采用洼地汇水范围内栅格的平均坡度值；

S22.所述小流域的积水体积V计算：积水体积V为洼地达到危险积水深度所需的水量；先设定积水深度，得到该积水深度对应的水面高程H，统计洼地低于水面高程H的栅格共有m个，且对应的高差为h_i(i＝1,2,...,m)；设定DEM的栅格大小为d；则积水体积V为：

7.根据权利要求6所述的城市内涝风险快速评估方法，其特征在于，所述步骤S3包括以下步骤：

S31.计算降雨量R：

设定一个或多个降雨情景，对于每个降雨情景，执行以下计算：

式中q为设定暴雨强度，t为设定降雨历时，P为设定重现期；其中A₁，C，b，n根据统计方法进行确定；计算得到该降雨情景暴雨强度q后，乘以降雨历时即得到其降雨量：R＝q×t；

S32.计算产流雨量R_f：

当目标区域具备实测的蒸散发量和土壤渗透能力数据时：利用已有的不透水面成果，或者采用监督分类方法从高分遥感影像或航空正射影像中提取不透水面的空间分布，并结合渗透模型计算通过土壤的下渗量；降雨量R扣除实测的蒸散发量和计算下渗量后即得到洼地小流域的产流雨量R_f；

当目标区域不具备实测的蒸散发量和土壤渗透能力数据时：根据洼地流域的土地利用现状确定不同土地利用类型的面积比例，再使用《室外排水设计规范》给出的径流系数推荐值，采用面积加权平均的方式计算每个洼地小流域的综合径流系数；降雨量R乘以综合径流系数即得到洼地小流域的产流雨量R_f；

S33.计算管道排水P：管道排水P根据目标区域的城市排水系统的设计资料确定或采用水力模型对目标区域的实际管道排水能力进行校核或根据目标区域的实际降雨事件对应的内涝情况进行经验性的反推；

S34.计算内涝雨量R_i：R_i＝R_f-P，内涝雨量为产流雨量R_f与管道排水P之差；内涝雨量R_i值应为正值，若为非正值则说明无内涝发生。

8.根据权利要求7所述的城市内涝风险快速评估方法，其特征在于，步骤S32所述的渗透模型为Horton下渗模型或Green-Ampt下渗模型。

9.一种城市内涝风险快速评估系统，其特征在于，包括：

小流域构建模块，用于根据DEM识别目标区域潜在的内涝发生位置，划定每一个内涝发生位置的汇水范围，并构成相应的小流域；

小流域属性计算模块，用于计算每个所述小流域的属性，包括其汇水面积A、平均坡度S、积水体积V；

内涝雨量计算模块，用于设定降雨情景，根据水量平衡模型，考虑包括排水系统设计标准、土壤渗透和蒸散发的影响，计算在地表形成积水的内涝雨量R_i；

内涝风险评估模块，用于以所述小流域构建模块得到的小流域为单元，以内涝发生快慢作为衡量标准，根据所述小流域属性计算模块及内涝雨量计算模块的计算结果计算目标区域的内涝风险Risk：Risk＝ln(R_i×A×S^1/2/V)。

Images