CN107194156B - Active water accumulation diffusion method for dynamically distributing water amount - Google Patents
Active water accumulation diffusion method for dynamically distributing water amount Download PDFInfo
- Publication number
- CN107194156B CN107194156B CN201710302748.0A CN201710302748A CN107194156B CN 107194156 B CN107194156 B CN 107194156B CN 201710302748 A CN201710302748 A CN 201710302748A CN 107194156 B CN107194156 B CN 107194156B
- Authority
- CN
- China
- Prior art keywords
- diffusion
- grid
- node
- water
- product
- 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.)
- Active
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
The invention discloses an active water accumulation diffusion algorithm for dynamically distributing water quantity, which is characterized in that after the water accumulation quantity of each diffusion source is calculated, diffusion is carried out in a central grid and eight grids around the central grid by taking a diffusion grid as a center according to a water quantity balance principle, the rationality judgment of a diffusion result is carried out after one circle of diffusion is finished, whether the formed accumulated water level is larger than the elevation value of a certain grid on the outer ring of the central grid is judged, if the accumulated water level is larger than the elevation value of any grid on the outer ring of the central grid, the diffusion is carried out again by adding one circle of diffusion range until the accumulated water level is smaller than the elevation value. After all the nodes are diffused, judging whether grids participate in the node diffusion process for many times, and if yes, carrying out merging diffusion processing on the nodes. In the diffusion process, the reasonability of the ponding path is continuously judged from the aspects of regional connectivity and terrain obstruction, water distribution is dynamically carried out, the subjectivity of manually selecting multi-flow-direction influence factors is avoided, and the accuracy of rainstorm ponding simulation is improved.
Description
Technical field
The present invention relates to the active ponding broadcast algorithms that a kind of water dynamically distributes, and belong to urban waterlogging monitoring and warning technology
Field.
Background technique
The flooded simulation process of urban storm product be broadly divided into rainfall simulation, produce runoff concentration calculation, heavy rain floods analysis.City
Storm water is certain caused by referring to after rainfall whereabouts urban area forms runoff due to hypsography low-lying, draining not in time etc.
The excess surface water of height of run-off.Urban storm floods analysis and refers to after acquiring urban storm ponding total amount, according to City Terrain
By ponding reasonable layout into regional space, obtain the possible depth of accumulated water in city and ponding floods situation.For this
The research of problem is carried out from two angles substantially: (1) giving the ponding under water level and flood analysis, these algorithms are all based on recurrence
What the angle of algorithm improvement carried out, it needs to give constant volume excess water position, but product excess water position is obtained in the presence of certain difficult in practical rainfall
Degree；(2) the product flood diffusion of flood volume known to: some algorithms do not break through the limitation of water catchment area boundary, when practical ponding is spread not
There are any boundary, the diffusion result that will lead to water catchment area boundary is unreasonable.Although there is no boundary limitations to ask for some algorithms
Topic, but its diffusion process may be such that the product excess water position for the diffusion zone that a node is formed is inconsistent, this and practical pool
The scene of formation is inconsistent.
Digital elevation model (Digital Elevation Model, vehicle economy M) is in digital form by certain knot
Structure organizes together, and indicates the number description that the model of actual landform feature space distribution and shape of mountain size rise and fall.Most
By millet cake (x, y) in a series, (x is the latitude coordinate value of the point of broad sense to basic DEM, and y is the latitude coordinate of the broad sense point
Value) position and its elevation Z that is associated formed, be with mathematical functional expression expression: Z=(x, y).DEM usually uses earth's surface regular
The elevation matrix that grid cell is constituted indicates that the DEM of broad sense further includes the number of all expression ground elevations such as contour, triangulation network
Word indicates.In the present invention, DEM refers to the set of the height value of each regular grid of covering survey region.
Summary of the invention
It is an object of the invention to overcome the shortcomings of the prior art, a kind of water that the fast accuracy of calculating speed is high is provided
Measure the active ponding broadcast algorithm dynamically distributed.
The technical solution adopted by the present invention are as follows: a kind of active ponding broadcast algorithm that water dynamically distributes, including following step
It is rapid:
(1) dem data is read, the height value of each grid in region is obtained；
(2) it generates and has uniquely determined the volume of each grid with the one-to-one grid number data of DEM grid, the data
Number, number data are read, and constitute matched relationship one by one according to latitude and longitude coordinates and dem data, and as DEM grid
One attribute of lattice；
(3) according to geographical location, the grid number and its product excess water amount of diffusion node (i.e. diffusion source) are determined, it equally will section
The number and water accumulating volume of point are as DEM grid attribute.
(4) since first node, ponding diffusion is successively carried out, the spillway discharge of present node is Vi, it is assumed that product flood is high
Journey is H, and the area of each DEM grid is area, the elevation of grid where node and the grid around other and the flooded elevation of product according to
Secondary is X1,X2…X9, then ((X-X is had according to principle of water balance1)+(X-X2)+…+(X-X9)) * area=Vi, obtain X.
(5) the elevation value difference for successively judging X Yu 16 adjacent grids of adjacent 8 grids, when all differences are both less than
Then diffusion terminates when the threshold value s of setting.
(6) when difference occur and being greater than threshold value s, range of scatter is expanded to as the 5*5 grid centered on node,
Then step (4) are repeated and obtain all node diffusion completions.(when a grid takes part in the diffusion of two or more nodes,
Its each elevation participates in calculating with original elevation.)
All node diffusions of (7) one wheels are completed, and are judged whether there is grid and are taken part in two or more node diffusion processes,
If it does not exist, then diffusion terminates, the flooded result of product to the end is obtained.
(8) two or more node diffusion processes are taken part in if there is grid, are then added the water of these nodes,
Water diffusion process is carried out in the external contact zone of its established flooded range of product, obtains new diffusion as a result, and updating this
The diffusion of a little grids obtains the flooded result of product as a result, until the judgement of all grids finishes.
Beneficial effects of the present invention: the method for the present invention goes out in diffusion process from ground domain connectivity and landform inhibition
Hair constantly judges the reasonability in ponding path, dynamically carries out water operation, avoids artificial selection multithread to impact factor
Subjectivity improves the accuracy of storm water simulation.
Detailed description of the invention
Fig. 1-Fig. 6 is the flooded situation schematic diagram of product in the specific embodiment of the invention.
Specific embodiment
Present invention will be further explained below with reference to the attached drawings and specific embodiments.
It is now assumed that survey region is as shown in Figure 1, each grid size is 10m*10m, and the number shown on figure is each grid
(unit m), what five-pointed star represented is node position to the height value of lattice unit.
Scene 1: it is now assumed that after a rainfall, the water accumulating volume of the node is 430m3, then product flood shown in Fig. 2 is obtained
Situation, final product excess water position are 8.7m, and less than the height value of all adjacent cells of surrounding, the diffusion of the node leaves it at that,
The flooded situation of the product of final each grid is as shown in Figure 2.
Scene 2: it is now assumed that after a rainfall, the water accumulating volume of the node is 730m3, then product flood shown in Fig. 3 is obtained
Situation carries out the result of first lap diffusion around node as shown in figure 3, there are the feelings that product excess water position is greater than its adjacent cells at this time
Condition, then the node needs to spread again, and range of scatter increases, and finally obtains situation as shown in Figure 4, and final product excess water position is
9.00714 taking approximation 9.01.
Scene 3: it is now assumed that there are nodes similar in two positions in the region, after a rainfall, the two nodes
Water accumulating volume is respectively 730m3And 110m3, the two nodes are carried out with product flood diffusion respectively, obtains the flooded situation of product as shown in Figure 5,
There is the phenomenon that intersecting in the flooded region of the product that the two nodes are formed, so needing to do the two nodes the flooded merging treatment of product, both
The water of the two nodes is added, new node is formed, the range of scatter of new node is the outer of the flooded region of the two nodes product
Connect polygon, finally obtain it is as shown in FIG. 6 product flood as a result, finally product excess water position be 9.0533, take approximation 9.05.
It should be pointed out that for those skilled in the art, without departing from the principle of the present invention,
Several improvements and modifications can also be made, these modifications and embellishments should also be considered as the scope of protection of the present invention.In the present embodiment not
The available prior art of specific each component part is realized.
Claims (1)
1. the active ponding method of diffusion that a kind of water dynamically distributes, it is characterised in that: the following steps are included:
(1) dem data is read, the height value of each grid in region is obtained；
(2) it generates and has uniquely determined the number of each grid with the one-to-one grid number data of DEM grid, the data,
Number data are read, constitute matched relationship one by one according to latitude and longitude coordinates and dem data, and as the one of DEM grid
A attribute；
(3) according to geographical location, the grid number and its product excess water amount of diffusion node are determined, equally by the number of node and product flood
Water is as DEM grid attribute；
(4) since first node, ponding diffusion is successively carried out, the spillway discharge of present node is Vi, it is assumed that product flood elevation is X,
The area of each DEM grid is area, and the elevation of grid and surrounding grid where node and the flooded elevation of product is followed successively by X1,
X2…X9, then ((X-X is had according to principle of water balance1)+(X-X2)+…+(X-X9)) * area=Vi, obtain X；
(5) the elevation value difference for successively judging X Yu 16 adjacent grids of adjacent 8 grids, when all differences are both less than set
Threshold value s when then diffusion terminate；
(6) when difference occur and being greater than threshold value s, range of scatter is expanded to as the 5*5 grid centered on node, then
Step (4) are repeated to complete until all nodes are spread；
All node diffusions of (7) one wheels are completed, and are judged whether there is grid and are taken part in two or more node diffusion processes, if
It is not present, then diffusion terminates, and obtains the flooded result of product to the end；
(8) two or more node diffusion processes are taken part in if there is grid, be then added the water of these nodes, at it
Water diffusion process is carried out in the external contact zone of the established flooded range of product, obtains new diffusion as a result, and updating these grid
The diffusion of lattice obtains the flooded result of product as a result, until the judgement of all grids finishes.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710302748.0A CN107194156B (en) | 2017-05-03 | 2017-05-03 | Active water accumulation diffusion method for dynamically distributing water amount |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710302748.0A CN107194156B (en) | 2017-05-03 | 2017-05-03 | Active water accumulation diffusion method for dynamically distributing water amount |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107194156A CN107194156A (en) | 2017-09-22 |
| CN107194156B true CN107194156B (en) | 2019-03-08 |
Family
ID=59873968
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710302748.0A Active CN107194156B (en) | 2017-05-03 | 2017-05-03 | Active water accumulation diffusion method for dynamically distributing water amount |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107194156B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109871621B (en) * | 2019-02-25 | 2021-06-01 | 中国水利水电科学研究院 | Urban rainstorm waterlogging catchment area analysis method |
| CN110686861B (en) * | 2019-10-24 | 2020-09-01 | 中国科学院地理科学与资源研究所 | Net rainfall grid format conversion method based on four-way elevation difference |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103870699A (en) * | 2014-03-21 | 2014-06-18 | 中国地质大学（武汉） | Hydrodynamics flood routing analogy method based on double-deck asynchronous iteration strategy |
| CN103927389A (en) * | 2014-04-30 | 2014-07-16 | 北京中有联科技有限公司 | Method for establishing flood disaster geographical analysis and evaluation dynamic model |
| CN104460343A (en) * | 2014-11-13 | 2015-03-25 | 长江水利委员会长江科学院 | Riverway flood inundation simulation method based on water level monitoring data |
| CN104898183A (en) * | 2015-05-29 | 2015-09-09 | 杭州辰青和业科技有限公司 | Modeling evaluation method for urban heavy rain inundation |
-
2017
- 2017-05-03 CN CN201710302748.0A patent/CN107194156B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103870699A (en) * | 2014-03-21 | 2014-06-18 | 中国地质大学（武汉） | Hydrodynamics flood routing analogy method based on double-deck asynchronous iteration strategy |
| CN103927389A (en) * | 2014-04-30 | 2014-07-16 | 北京中有联科技有限公司 | Method for establishing flood disaster geographical analysis and evaluation dynamic model |
| CN104460343A (en) * | 2014-11-13 | 2015-03-25 | 长江水利委员会长江科学院 | Riverway flood inundation simulation method based on water level monitoring data |
| CN104898183A (en) * | 2015-05-29 | 2015-09-09 | 杭州辰青和业科技有限公司 | Modeling evaluation method for urban heavy rain inundation |
Non-Patent Citations (3)
| Title |
|---|
| Numerical Simulation of Urban Waterlogging Based on F1oodArea Model;Fengchang Xue et.;《Advances in Meteorology》;20161231;第1-9页 * |
| 基于GIS和SWMM模型的城市暴雨积水模拟;黄国如等;《水资源与水工程学报》;20150831;第26卷(第4期);第1-6页 * |
| 基于GIS和SWMM的城市暴雨内涝淹没模拟分析;石赟赟等;《水电能源科学》;20140630;第32卷(第6期);第57-60、12页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107194156A (en) | 2017-09-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106884405B (en) | Inrush type mountain flood assay method for a kind of Cross Some Region Without Data | |
| CN109543275B (en) | A kind of city rainwash Two-dimensional numerical simulation method | |
| CN109871621B (en) | Urban rainstorm waterlogging catchment area analysis method | |
| CN102902893B (en) | Method for calculating rainfall ponding depth of catchment area based on DEM (digital elevation model) | |
| Rao et al. | A distributed model for real-time flood forecasting in the Godavari Basin using space inputs | |
| CN111538798B (en) | Urban catchment area refined extraction method considering DSM and DLG | |
| Nayak et al. | GIS based mapping of groundwater fluctuations in Bina basin | |
| CN103093114A (en) | Distributed-type river basin water deficit calculating method based on terrain and soil characteristics | |
| CN107194156B (en) | Active water accumulation diffusion method for dynamically distributing water amount | |
| Basarudin et al. | Event-based rainfall-runoff modelling of the Kelantan River Basin | |
| CN105760666A (en) | Catchment area critical rainfall calculating method based on DEM | |
| Topno et al. | SCS-CN Runoff estimation for Vindhyachal region using remote sensing and GIS | |
| CN108269199B (en) | Object-oriented small reservoir group time-space distribution type outflow calculation method | |
| CN115471078A (en) | Flood risk point assessment method and device based on urban water affair system | |
| CN104750936A (en) | GIS (Geographic Information System)-based urban hydrological run-off calculating method and system | |
| CN106096129B (en) | A kind of foot of the hill water surface scale analysis method calculated based on mountainous region charge for remittance | |
| Sabri et al. | Estimation of soil losses and reservoir sedimentation: a case study in tillouguite sub-basin (high atlas-Morocco). | |
| CN105279317A (en) | DEM based flat ground river network water flow direction estimation method | |
| Wang et al. | Spatiotemporal characteristics of cultural sites and their driving forces in the Ili River Valley during historical periods | |
| CN116305933B (en) | Simple slope yield confluence calculation method and device based on DEM data | |
| Mondal et al. | Evaluation of morphometric parameters of drainage networks derived from topographic map and digital elevation model using remote sensing and GIS | |
| CN110543692B (en) | Reconfiguration convergence simulation method based on underlying surface characteristics | |
| CN114385959A (en) | Method and device for dividing sub-basin units in dam-near region and storage medium | |
| Mokhtari et al. | Long term simulation of Shazand Plain Aquifer under changing resources and applications | |
| CN105740812A (en) | Digital surface model based city catchment area extraction method |
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 |