CN107194156A - 一种水量动态分配的有源积水扩散算法 - Google Patents
一种水量动态分配的有源积水扩散算法 Download PDFInfo
- Publication number
- CN107194156A CN107194156A CN201710302748.0A CN201710302748A CN107194156A CN 107194156 A CN107194156 A CN 107194156A CN 201710302748 A CN201710302748 A CN 201710302748A CN 107194156 A CN107194156 A CN 107194156A
- Authority
- CN
- China
- Prior art keywords
- diffusion
- grid
- node
- water
- grids
- 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
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
-
- 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
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Alarm Systems (AREA)
Abstract
本发明公开了一种水量动态分配的有源积水扩散算法,在计算出每个扩散源的积水量后,以扩散栅格为中心,根据水量平衡原理,在中心栅格及其周围八个栅格中扩散,每扩散完一圈进行扩散结果的合理性判断,判断已形成的积涝水位是否大于其外圈某个栅格的高程值,如果存在,则将扩散范围增加一圈再次进行扩散,直到积涝水位小于其外圈的任何一个栅格的高程值。在所有节点扩散完成之后,判断是否存在栅格参与了多次节点扩散过程,若有,则对这些节点进行合并扩散处理。在扩散过程中,从地域连通性及地形阻碍作用出发,不断判断积水路径的合理性,动态地进行水量分配,避免了人为选择多流向影响因子的主观性,提高了暴雨积水模拟的准确性。
Description
技术领域
本发明涉及一种水量动态分配的有源积水扩散算法,属于城市内涝监测预警技术领域。
背景技术
城市暴雨积涝模拟过程主要分为降雨过程模拟、产汇流计算、暴雨淹没分析。城市暴雨积水即指降雨下落城市区域形成径流后由于地势低洼、排水不及时等原因造成的一定径流深度的地面积水。城市暴雨淹没分析是指在求得城市暴雨积水总量后,根据城市地形将积水合理分布到区域空间中,得到城市可能的积水深度以及积水的淹没情况。对于这一问题的研究基本从两个角度进行:(1)给定水位下的积水淹没分析,这些算法都是基于递归算法改进的角度进行的,需要给定积涝水位,但是在实际降雨中积涝水位获取存在一定难度;(2)已知洪水体积的积涝扩散:有些算法未突破汇水区边界限制,实际积水扩散时并不存在任何边界,会导致汇水区边界处的扩散结果不合理。有些算法虽然不存在边界限制问题,但其扩散过程可能会使得一个节点形成的扩散区域的积涝水位不一致,这与实际水塘形成的情景不一致。
数字高程模型(Digital Elevation Model,简称DEM),是以数字的形式按一定结构组织在一起,表示实际地形特征空间分布的模型,也是地形形状大小起伏的数字描述。最基本的DEM由一系列地面点(x,y)(x为广义的该点的经度坐标值,y为广义该点的纬度坐标值)位置及其相联系的高程Z所组成,用数学函数式表达是:Z=(x,y)。DEM通常用地表规则网格单元构成的高程矩阵表示,广义的DEM还包括等高线、三角网等所有表达地面高程的数字表示。在本发明中,DEM是指覆盖研究区域的每个规则网格的高程值的集合。
发明内容
本发明的目的在于克服现有技术存在的不足,提供一种计算速度快准确度高的水量动态分配的有源积水扩散算法。
本发明采用的技术方案为:一种水量动态分配的有源积水扩散算法,包括以下步骤:
(1)读取DEM数据,获得区域每个栅格的高程值;
(2)生成与DEM栅格一一对应的网格编号数据,该数据唯一地确定了每个栅格的编号,读取编号数据,并根据经纬度坐标与DEM数据构成一一匹配的关系,并将其作为DEM栅格的一个属性;
(3)根据地理位置,确定扩散节点(即扩散源)的栅格编号及其积涝水量,同样将节点的编号与积水量作为DEM栅格属性。
(4)从第一个节点开始,依次进行积水扩散,当前节点的溢流量为Vi,假设积涝高程为H,每个DEM栅格的面积为area,节点所在栅格及其他周围的栅格与积涝高程的高程依次为X1,X2…X9,则根据水量平衡原理有((X-X1)+(X-X2)+…+(X-X9))*area=Vi,得到X。
(5)依次判断X与相邻8个栅格相邻的16个栅格的高程值差,当所有的差值都小于设定的阈值s时则扩散结束。
(6)当出现一个差值大于阈值s时,将扩散范围扩大到为以节点为中心的5*5栅格,然后重复步骤(4)得到所有节点扩散完成。(当一个栅格参与了两个及以上的节点扩散时,其每次的高程以原始高程参与计算。)
(7)一轮所有节点扩散完成,判断是否有栅格参与了两个及以上的节点扩散过程,如果不存在,则扩散结束,得到最后积涝结果。
(8)如果存在栅格参与了两个及以上的节点扩散过程,则将这些节点的水量相加,在其已形成的积涝范围的外接多边形内进行水量扩散过程,得到新的扩散结果,并更新这些栅格的扩散结果,直到所有栅格判断完毕,得到积涝结果。
本发明的有益效果:本发明方法在扩散过程中,从地域连通性及地形阻碍作用出发,不断判断积水路径的合理性,动态地进行水量分配,避免了人为选择多流向影响因子的主观性,提高了暴雨积水模拟的准确性。
附图说明
图1-图6为本发明具体实施方式中积涝情况示意图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步说明。
现假设研究区域如图1所示,每个栅格大小为10m*10m,图上所示的数字为每个栅格单元的高程值(单位为m),五角星代表的是节点所在位置。
情景1:现假设在一场降雨过后,该节点的积水量为430m3,则得到图2所示的积涝情况,最终的积涝水位为8.7m,小于周围所有相邻栅格的高程值,该节点的扩散到此结束,最终各栅格的积涝情况如图2所示。
情景2:现假设在一场降雨过后,该节点的积水量为730m3,则得到图3所示的积涝情况,围绕节点进行第一圈扩散的结果如图3所示,此时存在积涝水位大于其相邻栅格的情况,则该节点需要重新扩散,扩散范围增大,最终得到如图4所示的情况,最终积涝水位为9.00714,取近似值9.01。
情景3:现假设该区域存在两个位置相近的节点,在一场降雨过后,这两个节点的积水量分别为730m3和110m3,对这两个节点分别进行积涝扩散,得到如图5所示的积涝情况,这两个节点形成的积涝区域存在交叉的现象,所以需要对这两个节点做积涝合并处理,既把这两个节点的水量相加,形成新的节点,新节点的扩散范围为这两个节点积涝区域的外接多边形,最终得到如图6所示的积涝结果,最终积涝水位为9.0533,取近似值9.05。
应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (1)
1.一种水量动态分配的有源积水扩散算法,其特征在于:包括以下步骤:
(1)读取DEM数据,获得区域每个栅格的高程值;
(2)生成与DEM栅格一一对应的网格编号数据,该数据唯一地确定了每个栅格的编号,读取编号数据,并根据经纬度坐标与DEM数据构成一一匹配的关系,并将其作为DEM栅格的一个属性;
(3)根据地理位置,确定扩散节点的栅格编号及其积涝水量,同样将节点的编号与积水量作为DEM栅格属性;
(4)从第一个节点开始,依次进行积水扩散,当前节点的溢流量为Vi,假设积涝高程为H,每个DEM栅格的面积为area,节点所在栅格及其他周围的栅格与积涝高程的高程依次为X1,X2…X9,则根据水量平衡原理有((X-X1)+(X-X2)+…+(X-X9))*area=Vi,得到X;
(5)依次判断X与相邻8个栅格相邻的16个栅格的高程值差,当所有的差值都小于设定的阈值s时则扩散结束;
(6)当出现一个差值大于阈值s时,将扩散范围扩大到为以节点为中心的5*5栅格,然后重复步骤(4)得到所有节点扩散完成;
(7)一轮所有节点扩散完成,判断是否有栅格参与了两个及以上的节点扩散过程,如果不存在,则扩散结束,得到最后积涝结果;
(8)如果存在栅格参与了两个及以上的节点扩散过程,则将这些节点的水量相加,在其已形成的积涝范围的外接多边形内进行水量扩散过程,得到新的扩散结果,并更新这些栅格的扩散结果,直到所有栅格判断完毕,得到积涝结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710302748.0A CN107194156B (zh) | 2017-05-03 | 2017-05-03 | 一种水量动态分配的有源积水扩散方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710302748.0A CN107194156B (zh) | 2017-05-03 | 2017-05-03 | 一种水量动态分配的有源积水扩散方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107194156A true CN107194156A (zh) | 2017-09-22 |
CN107194156B CN107194156B (zh) | 2019-03-08 |
Family
ID=59873968
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710302748.0A Active CN107194156B (zh) | 2017-05-03 | 2017-05-03 | 一种水量动态分配的有源积水扩散方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107194156B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109871621A (zh) * | 2019-02-25 | 2019-06-11 | 中国水利水电科学研究院 | 城市暴雨内涝汇水区分析方法 |
CN110686861A (zh) * | 2019-10-24 | 2020-01-14 | 中国科学院地理科学与资源研究所 | 一种基于四向高程差的净雨量转栅格格式方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103870699A (zh) * | 2014-03-21 | 2014-06-18 | 中国地质大学(武汉) | 基于双层异步迭代策略的水动力学洪水演进模拟方法 |
CN103927389A (zh) * | 2014-04-30 | 2014-07-16 | 北京中有联科技有限公司 | 一种洪涝灾害地理分析评估动态模型的构建方法 |
CN104460343A (zh) * | 2014-11-13 | 2015-03-25 | 长江水利委员会长江科学院 | 一种基于水位监测数据的河道洪水淹没模拟方法 |
CN104898183A (zh) * | 2015-05-29 | 2015-09-09 | 杭州辰青和业科技有限公司 | 暴雨城市内涝建模评估方法 |
-
2017
- 2017-05-03 CN CN201710302748.0A patent/CN107194156B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103870699A (zh) * | 2014-03-21 | 2014-06-18 | 中国地质大学(武汉) | 基于双层异步迭代策略的水动力学洪水演进模拟方法 |
CN103927389A (zh) * | 2014-04-30 | 2014-07-16 | 北京中有联科技有限公司 | 一种洪涝灾害地理分析评估动态模型的构建方法 |
CN104460343A (zh) * | 2014-11-13 | 2015-03-25 | 长江水利委员会长江科学院 | 一种基于水位监测数据的河道洪水淹没模拟方法 |
CN104898183A (zh) * | 2015-05-29 | 2015-09-09 | 杭州辰青和业科技有限公司 | 暴雨城市内涝建模评估方法 |
Non-Patent Citations (3)
Title |
---|
FENGCHANG XUE ET.: "Numerical Simulation of Urban Waterlogging Based on F1oodArea Model", 《ADVANCES IN METEOROLOGY》 * |
石赟赟等: "基于GIS和SWMM的城市暴雨内涝淹没模拟分析", 《水电能源科学》 * |
黄国如等: "基于GIS和SWMM模型的城市暴雨积水模拟", 《水资源与水工程学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109871621A (zh) * | 2019-02-25 | 2019-06-11 | 中国水利水电科学研究院 | 城市暴雨内涝汇水区分析方法 |
CN110686861A (zh) * | 2019-10-24 | 2020-01-14 | 中国科学院地理科学与资源研究所 | 一种基于四向高程差的净雨量转栅格格式方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107194156B (zh) | 2019-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106884405B (zh) | 一种无资料地区溃堤型山洪灾害分析评价方法 | |
CN108874936B (zh) | 一种基于改进新安江模型的适用于山丘区的水文预报方法 | |
CN109543275B (zh) | 一种城区地表径流二维数值模拟方法 | |
CN110543984B (zh) | 一种基于水安全的城市水面率规划方法 | |
CN108388714A (zh) | 流域水系和城市管网耦合的平原河网城市洪水模拟方法 | |
CN109871621B (zh) | 城市暴雨内涝汇水区分析方法 | |
CN102902893B (zh) | 一种基于dem的汇水区降雨积水深度的计算方法 | |
CN108269199B (zh) | 一种面向对象的小水库群时空分布式出流计算方法 | |
CN110570517B (zh) | 一种基于下垫面特征的重配置产流模拟方法 | |
CN105760666A (zh) | 一种基于dem的汇水区临界雨量计算方法 | |
CN104750985B (zh) | 一种蓄水关系近似一致的复杂平原河网概化方法 | |
CN116305933B (zh) | 一种基于dem数据的坡面产汇流简易计算方法及装置 | |
CN107194156B (zh) | 一种水量动态分配的有源积水扩散方法 | |
CN104750936A (zh) | 一种基于gis的城市水文径流计算方法及*** | |
CN105279317A (zh) | 一种基于dem的平地河网水流方向估算方法 | |
CN110543692B (zh) | 一种基于下垫面特征的重配置汇流模拟方法 | |
CN110717232B (zh) | 一种考虑梯田淤地坝影响的子流域汇流模拟方法 | |
Du et al. | Impacts of socio-economic factors on sediment yield in the Upper Yangtze River | |
CN113869804B (zh) | 一种洪涝灾害下的电网设备风险预警方法及*** | |
CN111079065B (zh) | 一种时空动态优化的嵌套式自稳定汇流演算方法 | |
CN115859855B (zh) | 基于静水平衡原理的河网汇水域提取方法 | |
CN114528360B (zh) | 一种基于流速变化的流水分布快速生成方法 | |
CN103208135A (zh) | 一种基于cd-tin的城区地表特征点汇水量计算方法 | |
CN110110910B (zh) | 一种估算无资料地区库塘坝蓄量的方法 | |
CN118410728A (zh) | 一种基于光滑粒子法的人工景观水系水位预测方法、处理器、介质、设备 |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220513 Address after: 210000 building 7, Zhicheng Park, 6 Zhicheng Road, Jiangbei new district, Nanjing, Jiangsu. Patentee after: Jiangsu beilian Guoxin Technology Co.,Ltd. Address before: 210019 no.219, ningliu Road, Nanjing City, Jiangsu Province Patentee before: NANJING University OF INFORMATION SCIENCE & TECHNOLOGY |