CN109308560A - 一种城市暴雨内涝及积涝的风险预警方法 - Google Patents
一种城市暴雨内涝及积涝的风险预警方法 Download PDFInfo
- Publication number
- CN109308560A CN109308560A CN201810790661.7A CN201810790661A CN109308560A CN 109308560 A CN109308560 A CN 109308560A CN 201810790661 A CN201810790661 A CN 201810790661A CN 109308560 A CN109308560 A CN 109308560A
- Authority
- CN
- China
- Prior art keywords
- leakage rate
- net leakage
- pipe net
- grid
- water
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 67
- 230000003993 interaction Effects 0.000 claims abstract description 11
- 230000005540 biological transmission Effects 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000011144 upstream manufacturing Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 9
- 241000196324 Embryophyta Species 0.000 claims description 7
- 238000005086 pumping Methods 0.000 claims description 7
- 241001672694 Citrus reticulata Species 0.000 claims description 3
- 239000003673 groundwater Substances 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 239000002352 surface water Substances 0.000 claims description 3
- 238000013316 zoning Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 7
- 238000013461 design Methods 0.000 description 14
- 230000008569 process Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 3
- 230000006378 damage Effects 0.000 description 3
- 238000006424 Flood reaction Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000008595 infiltration Effects 0.000 description 2
- 238000001764 infiltration Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 101100366940 Mus musculus Stom gene Proteins 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 241000124033 Salix Species 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000005304 joining Methods 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000004078 waterproofing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
- G06Q50/265—Personal security, identity or safety
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Computer Security & Cryptography (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Alarm Systems (AREA)
Abstract
本发明提出了一种城市暴雨内涝及积涝的风险预警方法，包括：建立产流模型、一维管网模型和二维管网模型，产流模型根据降雨数据计算每一时间步的径流量，并将径流量数据传送给一维管网模型，在没有管网的区域，将径流量传送给二维管网模型，即假定径流量在有管网的地方先流入管网，没有管网的地方，直接在地表上流动；进行水量交互；产流模型、一维管网模型和二维管网模型分别各自推算至下一时间步，依次类推，从而完成对城市暴雨内涝及积涝进行的预警分析。
Description
技术领域
本发明涉及洪水灾害防控领域，特别涉及一种城市暴雨内涝及积涝的风险预警方法。
背景技术
近年来，随着城市建设的快速发展和强对流天气的增多，城市内涝及积涝呈现出一些新的特点，其造成的交通堵塞、河道蔓延、地质灾害等现象明显增加，城市建设及行业防灾工作面临新的课题。为此，需要根据城市积涝的新的灾害特征开展深入细致的风险评估研究，以评估结果指导各级生产管理部门进行有针对性的进行防灾、建设、运行和改造工作。然而，目前对暴雨灾害引起的城市内涝及积涝模型主要基于地表径流、水文动力学来实现，而对于气象部门而言，暴雨引起的内涝着重关注于对积涝灾害的预警及风险判断, 随着评估要求的提高，对评估的精细程度、空间的表现方式、时效性等提出了更高的要求, 传统的方法与技术不能满足需要，且现有的方法只能对降雨量和水位进行简单的估算，无法对暴雨带来的风险、损失及危害面积进行计算，进而导致暴雨前期人们的准备工作不充分，造成大量的人身事故和财产损失。
改革开放以来，由于城市化的快速发展，强降雨常造成城市内涝及积涝，致使交通堵塞、甚至造成人员伤亡，城市内涝和积涝已成为城市特有的气象灾害，迫切需要开展城市内涝监测预报预警服务。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此，本发明的目的在于提出一种城市暴雨内涝及积涝的风险预警方法。
2、为了实现上述目的，本发明的实施例提供一种城市暴雨内涝及积涝的风险预警方法，其中，包括如下步骤：
步骤S1，建立产流模型、一维管网模型和二维管网模型，产流模型根据降雨数据计算每一时间步的径流量，并将径流量数据传送给一维管网模型，在没有管网的区域，将径流量传送给二维管网模型，即假定径流量在有管网的地方先流入管网，没有管网的地方，直接在地表上流动；
步骤S2，进行水量交互，在每一时间步上比较河道断面水位或者管网节点水位与地表二维管网模型的水位的大小，如果河道或者管道的水位高于地表二维管网模型的水位，水流发生漫溢，水流从一维管网模型进入二维管网模型，如果地面水位高于河道或者管网节点水位，地表水回流至管网或者河道，则水流从二维管网模型进入一维管网模型；
步骤S3，进行水量交互之后，产流模型、一维管网模型和二维管网模型分别各自推算至下一时间步，依次类推。
进一步，所述建立产流模型：第一步，进行子汇水区划分工作，即将整个计算区域划分成小的子汇水区分别进行产流计算；第二步，确定产流参数；第三步，以子汇水区为单位进行产流计算。
进一步，所述建立一维管网模型又包括如下步骤：
步骤1，排除部分不在研究区域内、没有明确出水口以及重要属性信息缺失的管道和节点，对于部分重要管道，如果属性信息缺少，则根据上下游管道及周边管道信息进行人工补齐；
步骤2，排除部分完全重合的管道和节点；
步骤3，拓扑关系检查与重建；
步骤4，短管道合并；
步骤5，外层支管删除。
进一步，所述建立二维管网模型包括网格剖分，网格剖分又包括网格划分、属性赋值、排水泵站的处理和边界条件的确定。
进一步，所述网格的属性赋值主要包括：对网格逐一进行编号、类型、高程和糙率的赋值
进一步，所述网格的属性赋值包括如下步骤：首先需要对每个网格进行编号，赋予每个网格一个ID；其次，根据基础底图将网格划分为水面型网格、河道型网格、农田网格、非农田陆地网格等不同类型网格，分别赋予网格相应的糙率。
进一步，所述上边界条件的确定包括上游入流边界条件、下游出流边界条件和降雨边界条件。
根据本发明实施例的一种城市暴雨内涝及积涝的风险预警方法，通过建立产流模型、一维管网模型和二维管网模型，进行水量交互，并根据进行水量交互将产流模型、一维管网模型和二维管网模型分别各自推算至下一时间步，依次类推，从而完成对城市暴雨内涝及积涝进行的预警分析。
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1为本发明一种城市暴雨内涝及积涝的风险预警方法的模型构建示意图；
图2为本发明中产流模型的产流区划分示意图；
图3为本发明中产流模型的坡度计算结构的示意图；
图4为本发明中一维管网模型的河道断面输入和排水管网处理概化的分布图；
图5为本发明实施例中成都市3小时暴雨雨型的条形图；
图6为本发明实施例中成都市6小时暴雨雨型的条形图；
图7为本发明实施例中成都市望江楼降雨过程线示意图；
图8为本发明实施例中成都市望江楼站洪水过程计算值曲线图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
如图1所示，本发明实施例的一种城市暴雨内涝及积涝的风险预警方法，包括如下步骤：
步骤S1，建立产流模型、一维管网模型和二维管网模型，产流模型根据降雨数据计算每一时间步的径流量，并将径流量数据传送给一维管网模型，在没有管网的区域，将径流量传送给二维管网模型，即假定径流量在有管网的地方先流入管网，没有管网的地方，直接在地表上流动；
步骤S2，进行水量交互，在每一时间步上比较河道断面水位或者管网节点水位与地表二维管网模型的水位的大小，如果河道或者管道的水位高于地表二维管网模型的水位，水流发生漫溢，水流从一维管网模型进入二维管网模型，如果地面水位高于河道或者管网节点水位，地表水回流至管网或者河道，则水流从二维管网模型进入一维管网模型；
步骤S3，进行水量交互之后，产流模型、一维管网模型和二维管网模型分别各自推算至下一时间步，依次类推。
模型的建立：又包括产流模型的建立、一维管网模型的建立和二维管网模型的建立。
建立产流模型：
第一步，进行子汇水区划分工作，即将整个计算区域划分成小的子汇水区分别进行产流计算。子汇水区的划分方法采用根据管网节点的位置，采用泰森多边形方法批量划分子汇水区。这种方法在节点数目比较多的情况下，误差在可接受范围内，对于大型区域的建模计算比较适用。如图2所示，为产流区划分示意图，从图2中可以看出基本上每个节点均会对应一个产流区，部分没有节点的地方，则产流区与二维管网模型单元一致，产流直接流入二维管网模型单元格。
第二步，确定产流参数。产流计算实际包含两个部分：一是由降雨到净雨的过程，二是净雨汇入雨水口的过程。上述两个过程都放在了产流模型里计算，需要确定的参数主要有子汇水区不透水率、坡度和下渗率。参数确定主要采用三种方法：1)采用收集到的基础数据直接提取或者计算，这类参数选取受人工因素的影响比较小，通过这种方法确定的参数主要包括子汇水区的坡度和不透水率。结合遥感数据和DEM数据，利用ArcGIS和ENVI 计算和提取每个产流区的不透水率和坡度等重要产汇流参数。图3所示，为该区域坡度计算结果。2)根据实测数据进行率定得到，这类参数主要包括下渗率参数。3)根据经验方法确定。部分参数属于概念性参数，很难通过物理测量来获取，又属于对结果影响较小的参数，在参数率定时，比较小可能会调整的，因而此类参数主要根据经验来选定。
第三步，以子汇水区为单位进行产流计算。
建立一维管网模型：一维管网模型包括河道断面的输入和排水管网的处理概化。
1)排除部分不在研究区域内、没有明确出水口以及重要属性信息缺失的管道和节点，对于部分重要管道，如果属性信息缺少，则根据上下游管道及周边管道信息进行人工补齐。
2)排除部分完全重合的管道和节点。对于较大的城市，排水管道信息可能来自不同年代或不同的测绘公司，在不同的数据来源的边界处，可能会存在管道重复测量、重复绘制等问题，需要进行识别和处理。
3)拓扑关系检查与重建。管道与管道、管道与节点之间的相互连接关系称之为管道的拓扑关系，这些拓扑关系在管道测量时就应该明确给出，但有时由于测量或者数据输入出现偏差，导致部分管道拓扑关系出错，对于这类管道需要根据地理位置重建拓扑关系，部分无法重建拓扑关系的管线进行删除处理。
4)短管道合并。管道数目过多对模型稳定，需要对于一些特别短的管道进行合并处理，合并处理时，只合并管径相同的管道。
5)外层支管删除。对于管径小于某个值的最外层管道进行删除简化处理。
如图4为一维管网模型的河道断面输入和排水管网处理概化的分布图。
建立二维管网模型：包括网格剖分，网格的属性赋值主要包括：对网格逐一进行编号、类型、高程、糙率等的赋值。首先需要对每个网格进行编号，赋予每个网格一个ID；其次，根据基础底图将网格划分为水面型网格、河道型网格、农田网格、非农田陆地网格等不同类型网格，分别赋予网格相应的糙率。网格剖分又包括网格划分、属性赋值、排水泵站的处理和边界条件的确定，其中属性赋值主要是确定网格单元的高程、糙率以及修正率等参数。
1)网格属性赋值
网格高程：由DEM数据提取，将DEM与网格关联，取每个网格内所有DEM栅格高程的平均值作为网格形心点的高程；
糙率：按照草地、居民区和街道等不同类型，赋予网格相应的糙率；
面积修正率：提取每个网格内的居民地面积，与网格面积的比值作为该网格的面积修正率。
2)通道属性赋值
将道路、堤沿程的通道挑选出来；
道路通道赋值：将道路通道与DEM空间关联，提取通道高程，并与高程点图层做比较，修改不合理属性。
堤防通道赋值：由于数据比通道稀疏，利用收集到的堤防高程数据，根据距离做插值，赋值给每条堤防通道。
河道型特殊通道：较小的河道概化为特殊性通道，利用已有的河道断面数据，计算各个断面的等效断面，并将河底高程、左堤高程、右堤高程、上宽、下宽分别进行插值，赋值给每条河道型特殊通道。
阻水型特殊通道与其它特殊通道的交口处，根据实际情况设置涵洞，有缺口的特殊通道则采用宽顶堰的方式进行计算。
3)排水泵站的处理
水深大于0.3米时，泵站开启，在泵站的抽排能力范围内，具体每时刻抽排量根据泵站所在位置的网格单元水量来确定，直到水深小于一定值才停止。
4)边界条件的确定
边界类型的确定：模型边界条件主要依据成都市河流分布及洪水来源确定。现以成都市城区洪水分离计算案例为实施例。成都市的主要洪水来源为暴雨内涝，年降水量为800～1400mm，按地域分布由西北向东南递减。降雨过程中的短时强度大，历时长。中心城区小范围、高强度的局部暴雨时有发生，很容易导致内涝。据统计，五十年来，区间暴雨洪水发生了13次，约平均四年发生一次。成都市另一洪水来源为岷江洪水，虽然经过紫平铺水库和都江堰水利工程的调控，成都市的外洪风险大大降低，但在上游暴雨强度很大的情况下，来流量仍会带来内部河流包括府河、清水河、江安河、沙河等河流发生漫堤的风险。
根据成都市洪水来源以及河流分布情况，模型边界条件类型主要为3种：1)上游入流边界条件；2)下游出流边界条件；3)降雨边界条件。
成都市城区的主要边界条件
上游入流边界
成都市外环路以内的十八条中、小河流控制断面的设计洪水过程线按以下两种情况分析计算。
(1)实时测量洪水过程线
(2)利用暴雨推求设计洪水过程线
利用暴雨资料和洪水资料对照分析，成都市小河流的设计洪水过程线有三种类型：1) 涨水慢，退水快；2)涨水快，退水慢；3)前两种类型的平均情况。
对成都市岷、沱江水系的下垫面条件及实测洪水资料分析，属于府河水系的红光七斗、马河、沱江河、杨柳河、磨底河、苏坡排洪渠、黄堰河按第二种类型(涨水快，退水慢)计算各控制断面500年一遇、200年一遇、100年一遇以及50年一遇的设计洪水采用的过程线。
城区主要河道200年一遇、100年一遇和50年一遇设计洪峰流量和采用的洪水过程线如下表：
城区主要河道上游设计洪水过程(P＝0.5％)
城区主要河道上游设计洪水过程(P＝1％)
城区主要河道上游设计洪水过程(P＝2％)
(3)降雨量推求
成都市24小时，6小时和3小时设计暴雨如下表所示。
成都市24小时暴雨量计算成果表
注：任意设计频率的雨量按照Cs＝3.5Cv的皮尔逊III型频率计算。
成都市6小时暴雨量计算成果表
成都市3小时暴雨量计算成果表
(4)暴雨雨型推求
采用Pilgrim&Cordery法推求设计暴雨雨型，依据Pilgrim&Cordery法的原理，将峰值放置在级序最大的位置上，结合级序和比例，则可分析得到与360分钟、180分钟和60分钟的雨型分配比例相对应的雨峰发生位置r，其分别为0.13、0.28、0.33。
1)180分钟设计暴雨雨型：(Pilgrim&Cordery法)180分钟设计暴雨雨型，如图5所示，峰型为单峰，峰时发生在第10位，单位时段为5分钟。
2)360分钟设计暴雨雨型：(Pilgrim&Cordery法)360分钟设计暴雨雨型，如图6所示，峰型为单峰，峰时发生在第9位，单位时段为5分钟。
3)24小时暴雨雨型：24小时设计暴雨的雨型采用典型降雨过程按等比例分配得到，典型降雨选用成都市2013年7月9日的降雨，具体过程如下：
24小时设计暴雨雨型
(5)模型验证
采用成都市2013年7月10-11日长历时暴雨对模型的合理性进行检验，降雨数据采用成都市城区范围内望江楼、永宁等在内的15个站点的实测降雨数据，这些站点包括成都市气象部门和水文部门设立的雨量站和水文站，基本能够反映成都城区范围内降雨分布情况。如图7所示，给出了望江楼站实测暴雨过程线。望江楼站洪水过程计算值如图8所示，望江楼站实测和模型计算的洪峰流量值、洪峰到达时间的对比情况见下表。从表中的结果来看，模型总体上跟实测值较为接近，表明模型产汇流计算参数选取是合理的，模型具有较好的精度。与此同时，由于模型上边界的流量过程并没有实测值，对模型结果精度，特别是退水段产生了一定的影响。
模型计算值与实测值比较
(6)洪水计算方案汇总
成都城区洪水分析计算方案见下表。
1)暴雨内涝
以暴雨为来源的计算方案列表
2)河道洪水
以河道洪水为来源的计算方案列表
| 序号 | 方案情况 |
| 16 | 50年一遇外洪 |
| 17 | 100年一遇外洪 |
| 18 | 200年一遇外洪 |
3)洪水遭遇
以暴雨和河道洪水组合为来源的计算方案列表
| 序号 | 方案情况 |
| 19 | 50年一遇洪水遭遇20年一遇24h暴雨 |
4)典型洪水
典型洪水的计算方案列表
| 序号 | 方案情况 |
| 20 | 遭遇与2011年7月同量级的暴雨洪水(20110703) |
| 21 | 遭遇与2013年7月同量级的暴雨洪水(20130710) |
(2)洪水计算成果
1)淹没过程分析
当计算开始6小时，苏坡排洪渠发生漫溢，但淹没水深相对较小，当计算开始12小时后，苏坡排洪渠及下游地区淹没范围有所扩大，在18小时后，清水河发生比较严重的河道洪水漫溢，府河上游段也有一部分水漫至地表，当计算至24小时后，局部水深有所增加，大部分地区水深开始回退。由于清水河的防洪标准相对较低，从结果上看，总体上符合实际。
2)方案比较
以3小时不同重现期的降雨淹没情况为例进行分析，成都城区范围内遭遇3小时5年一遇、20年一遇、50年一遇暴雨时的淹没情况。从这些图中可以看出，随着降雨由5年一遇增大到50年一遇，淹没范围逐渐扩大。根据统计数据，淹没面积在5年一遇时为15.2km2， 20年一遇暴雨时为42.3km2，当遭遇50年一遇暴雨时，淹没面积达到了75.2km2，占到整个计算范围的10％以上。
根据本发明实施例的一种城市暴雨内涝及积涝的风险预警方法，通过建立产流模型、一维管网模型和二维管网模型，进行水量交互，并根据进行水量交互将产流模型、一维管网模型和二维管网模型分别各自推算至下一时间步，依次类推，从而完成对城市暴雨内涝及积涝进行的预警分析。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。
Claims (7)
1.一种城市暴雨内涝及积涝的风险预警方法，其特征在于，包括如下步骤：
步骤S1，建立产流模型、一维管网模型和二维管网模型，产流模型根据降雨数据计算每一时间步的径流量，并将径流量数据传送给一维管网模型，在没有管网的区域，将径流量传送给二维管网模型，即假定径流量在有管网的地方先流入管网，没有管网的地方，直接在地表上流动；
步骤S2，进行水量交互，在每一时间步上比较河道断面水位或者管网节点水位与地表二维管网模型的水位的大小，如果河道或者管道的水位高于地表二维管网模型的水位，水流发生漫溢，水流从一维管网模型进入二维管网模型，如果地面水位高于河道或者管网节点水位，地表水回流至管网或者河道，则水流从二维管网模型进入一维管网模型；
步骤S3，进行水量交互之后，产流模型、一维管网模型和二维管网模型分别各自推算至下一时间步，依次类推。
2.如权利要求1所述的一种城市暴雨内涝及积涝的风险预警方法，其特征在于，所述建立产流模型：第一步，进行子汇水区划分工作，即将整个计算区域划分成小的子汇水区分别进行产流计算；第二步，确定产流参数；第三步，以子汇水区为单位进行产流计算。
3.如权利要求1所述的一种城市暴雨内涝及积涝的风险预警方法，其特征在于，所述建立一维管网模型又包括如下步骤：
步骤1，排除部分不在研究区域内、没有明确出水口以及重要属性信息缺失的管道和节点，对于部分重要管道，如果属性信息缺少，则根据上下游管道及周边管道信息进行人工补齐；
步骤2，排除部分完全重合的管道和节点；
步骤3，拓扑关系检查与重建；
步骤4，短管道合并；
步骤5，外层支管删除。
4.如权利要求1所述的一种城市暴雨内涝及积涝的风险预警方法，其特征在于，所述建立二维管网模型包括网格剖分，网格剖分又包括网格划分、属性赋值、排水泵站的处理和边界条件的确定。
5.如权利要求4所述的一种城市暴雨内涝及积涝的风险预警方法，其特征在于，所述网格的属性赋值主要包括：对网格逐一进行编号、类型、高程和糙率的赋值。
6.如权利要求4所述的一种城市暴雨内涝及积涝的风险预警方法，其特征在于，所述网格的属性赋值包括如下步骤：首先需要对每个网格进行编号，赋予每个网格一个ID；其次，根据基础底图将网格划分为水面型网格、河道型网格、农田网格、非农田陆地网格等不同类型网格，分别赋予网格相应的糙率。
7.如权利要求1所述的一种城市暴雨内涝及积涝的风险预警方法，其特征在于，所述上边界条件的确定包括上游入流边界条件、下游出流边界条件和降雨边界条件。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810790661.7A CN109308560B (zh) | 2018-07-18 | 2018-07-18 | 一种城市暴雨内涝及积涝的风险预警方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810790661.7A CN109308560B (zh) | 2018-07-18 | 2018-07-18 | 一种城市暴雨内涝及积涝的风险预警方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109308560A true CN109308560A (zh) | 2019-02-05 |
| CN109308560B CN109308560B (zh) | 2019-11-15 |
Family
ID=65225884
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810790661.7A Active CN109308560B (zh) | 2018-07-18 | 2018-07-18 | 一种城市暴雨内涝及积涝的风险预警方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109308560B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110232477A (zh) * | 2019-06-05 | 2019-09-13 | 中国水利水电科学研究院 | 一种反应下垫面变化影响的洪水预警方法和装置 |
| CN110274656A (zh) * | 2019-06-20 | 2019-09-24 | 福州市规划设计研究院 | 一种城市内河水位预报预警方法 |
| CN112287285A (zh) * | 2020-10-29 | 2021-01-29 | 南昌工程学院 | 基于城市内涝基础信息的蓄涝水面率分区处理方法 |
| CN114547922A (zh) * | 2022-04-28 | 2022-05-27 | 南京信息工程大学 | 基于多边形网格的城市内涝分析方法、装置和存储介质 |
| CN115796453A (zh) * | 2022-11-28 | 2023-03-14 | 中国水利水电科学研究院 | 一种市政下水管道人孔盖排气效果评估方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140156232A1 (en) * | 2012-11-30 | 2014-06-05 | International Business Machines Corporation | Methods, systems and computer program storage devices for generating a flooding forecast |
| CN104898183A (zh) * | 2015-05-29 | 2015-09-09 | 杭州辰青和业科技有限公司 | 暴雨城市内涝建模评估方法 |
| CN106373070A (zh) * | 2016-08-29 | 2017-02-01 | 西安理工大学 | 一种应对城市暴雨内涝的四预方法 |
| CN107220496A (zh) * | 2017-05-26 | 2017-09-29 | 上海市气象灾害防御技术中心 | 一种城市暴雨内涝评估建模方法 |
| CN107239607A (zh) * | 2017-05-31 | 2017-10-10 | 中国水利水电科学研究院 | 一种模型要素与计算方案管理方法 |
| CN108182543A (zh) * | 2018-01-17 | 2018-06-19 | 福建四创软件有限公司 | 一种精细化网格内涝水淹预报方法 |
-
2018
- 2018-07-18 CN CN201810790661.7A patent/CN109308560B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140156232A1 (en) * | 2012-11-30 | 2014-06-05 | International Business Machines Corporation | Methods, systems and computer program storage devices for generating a flooding forecast |
| CN104898183A (zh) * | 2015-05-29 | 2015-09-09 | 杭州辰青和业科技有限公司 | 暴雨城市内涝建模评估方法 |
| CN106373070A (zh) * | 2016-08-29 | 2017-02-01 | 西安理工大学 | 一种应对城市暴雨内涝的四预方法 |
| CN107220496A (zh) * | 2017-05-26 | 2017-09-29 | 上海市气象灾害防御技术中心 | 一种城市暴雨内涝评估建模方法 |
| CN107239607A (zh) * | 2017-05-31 | 2017-10-10 | 中国水利水电科学研究院 | 一种模型要素与计算方案管理方法 |
| CN108182543A (zh) * | 2018-01-17 | 2018-06-19 | 福建四创软件有限公司 | 一种精细化网格内涝水淹预报方法 |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110232477A (zh) * | 2019-06-05 | 2019-09-13 | 中国水利水电科学研究院 | 一种反应下垫面变化影响的洪水预警方法和装置 |
| CN110274656A (zh) * | 2019-06-20 | 2019-09-24 | 福州市规划设计研究院 | 一种城市内河水位预报预警方法 |
| CN110274656B (zh) * | 2019-06-20 | 2021-05-04 | 福州市规划设计研究院集团有限公司 | 一种城市内河水位预报预警方法 |
| CN112287285A (zh) * | 2020-10-29 | 2021-01-29 | 南昌工程学院 | 基于城市内涝基础信息的蓄涝水面率分区处理方法 |
| CN112287285B (zh) * | 2020-10-29 | 2023-04-04 | 南昌工程学院 | 基于城市内涝基础信息的蓄涝水面率分区处理方法 |
| CN114547922A (zh) * | 2022-04-28 | 2022-05-27 | 南京信息工程大学 | 基于多边形网格的城市内涝分析方法、装置和存储介质 |
| CN114547922B (zh) * | 2022-04-28 | 2022-07-29 | 南京信息工程大学 | 基于多边形网格的城市内涝分析方法、装置和存储介质 |
| CN115796453A (zh) * | 2022-11-28 | 2023-03-14 | 中国水利水电科学研究院 | 一种市政下水管道人孔盖排气效果评估方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109308560B (zh) | 2019-11-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109308560B (zh) | 一种城市暴雨内涝及积涝的风险预警方法 | |
| CN110298076A (zh) | 一种基于gis和swmm的城市内涝智能建模及分析方法 | |
| CN107085658B (zh) | 一种山洪灾害成灾时间确定方法 | |
| CN111581767B (zh) | 一种管网-河流耦合模型校核特征参数率定方法 | |
| CN104462774A (zh) | 基于水箱模型的城市道路及低洼地区积水预报方法 | |
| Mehta et al. | Water productivity enhancement through controlling the flood inundation of the surrounding region of Navsari Purna river, India | |
| CN116306340A (zh) | 一种模拟不同工况下的城市内涝风险分布的方法 | |
| Hsieh et al. | An assessment of structural measures for flood-prone lowlands with high population density along the Keelung River in Taiwan | |
| CN116561957B (zh) | 一种基于数值模拟的地下空间内涝评估方法 | |
| Tunas | The application of ITS-2 model for flood hydrograph simulation in large-size rainforest watershed, Indonesia | |
| CN109472102B (zh) | 一种基于蚂蚁河流域的桥涵构造推算方法 | |
| CN110889185A (zh) | 一种小流域的洪峰流量分析方法及其应用 | |
| Afrin et al. | Assessment of future flow of Hatirjheel-Begunbari drainage system due to climate change | |
| Szydłowski et al. | Stormwater and snowmelt runoff storage control and flash flood hazard forecasting in the urbanized coastal basin | |
| Ghazavi et al. | Modelling and assessment of urban flood hazards based on rainfall intensity-duration-frequency curves reformation | |
| Narayan et al. | GIS supported geomorphologic instantaneous unit hydrograph (GIUH) of Varuna river basin using geomorphological characteristics | |
| CN113887053A (zh) | 一种面向管网水流计算的市政排水资料质量评估方法及系统 | |
| CN112884279B (zh) | 一种沟道梯级土地整治工程暴雨洪水安全风险评估方法 | |
| Pratik et al. | Development of Flood Hydrograph for Un-gauged Catchment using SUH Approach–A Case Study of Patalganga River Basin | |
| Ekaningtyas | Flood Inundation Prediction of Logung River due to the Break of Logung Dam | |
| Hraniciuc et al. | Coupled Hydrological and Hydraulic Modeling of Dangerous Meteorological Phenomena | |
| Adib et al. | Diagnosis of erodible locations in river bends using a combined method (GIS, RS and the CCHE2D Model)(case study: the Karkheh River in Iran) | |
| Neupane et al. | Modeling of urban storm water drainage using HYKAS for Lamachaur, Pokhara, Nepal | |
| TIAN et al. | Hydrologic Correlation Analysis of Highway Bridge Community Along River | |
| Bohorquez et al. | Parametric study of trends in flood stages over time in the regulated Guadalquivir River (years 1910–2016) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |