CN107194156A - Active water accumulation diffusion algorithm for water amount dynamic distribution - Google Patents
Active water accumulation diffusion algorithm for water amount dynamic distribution 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
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 a kind of active ponding broadcast algorithm of water dynamically distributes, belong to urban waterlogging monitoring and warning technology
Field.
Background technology
The flooded simulation process of urban storm product is broadly divided into rainfall simulation, production runoff concentration calculation, heavy rain and floods analysis.City
Storm water be refer to rainfall whereabouts urban area formed after runoff due to hypsography low-lying, draining not in time etc. reason cause it is certain
The excess surface water of height of run-off.Urban storm floods analysis and referred to after urban storm ponding total amount is tried to achieve, 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 basic to be carried out from two angles：(1) ponding given under water level floods analysis, and these algorithms are all based on recurrence
It is that the angle of algorithm improvement is carried out, it is necessary to give constant volume excess water position, but in actual rainfall product excess water position obtain exist it is certain difficult
Degree；(2) the product flood diffusion of flood volume known to：Some algorithms do not break through the limitation of water catchment area border, during actual ponding diffusion not
There is any border, the diffusion result of water catchment area boundary can be caused unreasonable.Although border is not present in some algorithms, limitation is asked
Topic, but its diffusion process may be such that the product excess water position of the diffusion zone of node formation is inconsistent, this and actual pool
The scene of formation is inconsistent.
Digital elevation model (Digital Elevation Model, vehicle economy M), is in digital form by certain knot
Structure is organized together, and represents the model of actual landform feature space distribution, is also the numeral description that shape of mountain size rises and falls.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 constituted, be with mathematical functional expression expression：Z=(x, y).DEM is generally with earth's surface rule
The elevation matrix that grid cell is constituted represents that the DEM of broad sense also includes the number of all expression ground elevations such as contour, the triangulation network
Word is represented.In the present invention, DEM refers to the set for covering the height value of each regular grid of survey region.
The content of the invention
It is an object of the invention to overcome the shortcomings of prior art presence, there is provided a kind of high water of the fast degree of accuracy of calculating speed
Measure the active ponding broadcast algorithm of dynamically distributes.
The technical solution adopted by the present invention is：A kind of active ponding broadcast algorithm of water dynamically distributes, including following step
Suddenly：
(1) dem data is read, the height value of each grid in region is obtained；
(2) generation and the one-to-one grid numbering data of DEM grids, the data uniquely determine the volume of each grid
Number, numbering data are read, and the relation matched one by one is constituted with dem data according to latitude and longitude coordinates, and as DEM grids
An attribute；
(3) according to geographical position, the grid numbering and its product excess water amount of diffusion node (spreading source) are determined, equally will section
The numbering of point is used as DEM grid attributes with water accumulating volume.
(4) since first node, ponding diffusion is carried out successively, the spillway discharge of present node is Vi, it is assumed that product flood is high
Journey is H, the areas of each DEM grids be grid where area, node and the grid around other with the elevation of 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 of X 16 grids adjacent with adjacent 8 grids is judged successively, when all differences are both less than
Then spread and terminate during the threshold value s of setting.
(6) when difference occur and being more than threshold value s, range of scatter is expanded to as the 5*5 grids centered on node,
Then repeat step (4) obtains all node diffusion completions.(when a grid take part in the node diffusion of two and the above,
Its each elevation participates in calculating with original elevation.)
All node diffusions of (7) one wheels are completed, and determine whether that grid take part in the node diffusion process of two and the above,
If it does not exist, then diffusion terminates, the flooded result of last product is obtained.
(8) the node diffusion process of two and the above is take part in if there is grid, then is added the water of these nodes,
Water diffusion process is carried out in the external contact zone of its established flooded scope of product, new diffusion result is obtained, and update this
The diffusion result of a little grids, until all grids judge to finish, obtains the flooded result of product.
Beneficial effects of the present invention：The inventive method goes out in diffusion process from ground domain connectivity and landform inhibition
Hair, constantly judge the reasonability in ponding path, dynamically carry out water operation, it is to avoid flow to factor of influence artificial selection more
Subjectivity, improves the accuracy of storm water simulation.
Brief description of the drawings
Fig. 1-Fig. 6 is the flooded situation schematic diagram of product in the specific embodiment of the invention.
Embodiment
The invention will be further described with reference to the accompanying drawings and detailed description.
It is now assumed that survey region is as shown in Figure 1, each grid size is 10m*10m, and the numeral shown on figure is each grid
The height value (unit is m) of lattice unit, what five-pointed star was represented is node position.
Scene 1：It is now assumed that after a rainfall, the water accumulating volume of the node is 430m3, then the product flood shown in Fig. 2 is obtained
Situation, final product excess water position is 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 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 the product flood shown in Fig. 3 is obtained
Situation, carries out the result of first lap diffusion as shown in figure 3, now there are the feelings that product excess water position is more than its adjacent cells around node
Condition, then the node need to spread again, range of scatter increase finally gives situation as shown in Figure 4, and final product excess water position is
9.00714, take approximation 9.01.
Scene 3：It is now assumed that there is the close node in two positions, after a rainfall, the two nodes in the region
Water accumulating volume is respectively 730m3And 110m3, the two nodes are carried out respectively to accumulate flooded diffusion, the flooded situation of product as shown in Figure 5 is obtained,
There is the phenomenon intersected in the flooded region of product of the two nodes formation, 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 accumulates the outer of flooded region for the two nodes
Polygon is connect, the flooded result of product as shown in Figure 6 is finally given, final product excess water position is 9.0533, takes approximation 9.05.
It should be pointed out that for those skilled in the art, under the premise without departing from the principles of the invention,
Some improvements and modifications can also be made, these improvements and modifications also should be regarded as protection scope of the present invention.In the present embodiment not
Clear and definite each part can use prior art to be realized.
Claims (1)
1. a kind of active ponding broadcast algorithm of water dynamically distributes, it is characterised in that：Comprise the following steps：
(1) dem data is read, the height value of each grid in region is obtained；
(2) generation and the one-to-one grid numbering data of DEM grids, the data uniquely determine the numbering of each grid,
Numbering data are read, and the relation matched one by one is constituted with dem data according to latitude and longitude coordinates, and as DEM grids
One attribute；
(3) according to geographical position, the grid numbering and its product excess water amount of diffusion node are determined, equally by the numbering and ponding of node
Amount is used as DEM grid attributes；
(4) since first node, ponding diffusion is carried out successively, the spillway discharge of present node is Vi, it is assumed that product flood elevation is H,
The area of each DEM grids is that the elevation of grid where area, node and the grid around other 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 of X 16 grids adjacent with adjacent 8 grids is judged successively, when all differences are both less than setting
Threshold value s when then spread and terminate；
(6) when difference occur and being more than threshold value s, range of scatter is expanded to as the 5*5 grids centered on node, then
Repeat step (4) obtains all node diffusions and completed；
All node diffusions of (7) one wheels are completed, and determine whether that grid take part in the node diffusion process of two and the above, if
It is not present, then diffusion terminates, obtains the flooded result of last product；
(8) the node diffusion process of two and the above is take part in if there is grid, then is added the water of these nodes, at it
Water diffusion process is carried out in the external contact zone of the established flooded scope of product, new diffusion result is obtained, and update these grid
The diffusion result of lattice, until all grids judge to finish, obtains the flooded result of product.
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 true CN107194156A (en) | 2017-09-22 |
| CN107194156B 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) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109871621A (en) * | 2019-02-25 | 2019-06-11 | 中国水利水电科学研究院 | Urban rainstorm waterlogging water catchment area analysis method |
| CN110686861A (en) * | 2019-10-24 | 2020-01-14 | 中国科学院地理科学与资源研究所 | 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 |
|---|
| 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 (en) * | 2019-02-25 | 2019-06-11 | 中国水利水电科学研究院 | Urban rainstorm waterlogging water catchment area analysis method |
| CN110686861A (en) * | 2019-10-24 | 2020-01-14 | 中国科学院地理科学与资源研究所 | Net rainfall grid format conversion method based on four-way elevation difference |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107194156B (en) | 2019-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106884405B (en) | Inrush type mountain flood assay method for a kind of Cross Some Region Without Data | |
| CN108874936B (en) | Hydrological forecasting method applicable to hilly area based on improved Xinanjiang model | |
| CN109543275B (en) | A kind of city rainwash Two-dimensional numerical simulation method | |
| CN110543984B (en) | Urban water surface rate planning method based on water safety | |
| CN108388714A (en) | The plain river network city flood simulation method of basin water system and urban pipe network coupling | |
| CN109871621B (en) | Urban rainstorm waterlogging catchment area analysis method | |
| CN108269199B (en) | Object-oriented small reservoir group time-space distribution type outflow calculation method | |
| CN102902893A (en) | Method for calculating rainfall ponding depth of catchment area based on DEM (digital elevation model) | |
| CN110570517B (en) | Reconfiguration runoff yield simulation method based on underlying surface characteristics | |
| CN105760666A (en) | Catchment area critical rainfall calculating method based on DEM | |
| CN104750985B (en) | A kind of approximate consistent complicated plain river network of retaining relation generally changes method | |
| CN116305933B (en) | Simple slope yield confluence calculation method and device based on DEM data | |
| CN107194156B (en) | Active water accumulation diffusion method for dynamically distributing water amount | |
| CN104750936A (en) | GIS (Geographic Information System)-based urban hydrological run-off calculating method and system | |
| CN113869804B (en) | Power grid equipment risk early warning method and system under flood disaster | |
| CN105279317A (en) | DEM based flat ground river network water flow direction estimation method | |
| CN110543692B (en) | Reconfiguration convergence simulation method based on underlying surface characteristics | |
| CN117789416A (en) | Regional storm disaster risk early warning method and device | |
| CN110717232B (en) | Sub-basin confluence simulation method considering influence of terrace check dam | |
| Du et al. | Impacts of socio-economic factors on sediment yield in the Upper Yangtze River | |
| CN111079065B (en) | Nested self-stabilization confluence calculation method for space-time dynamic optimization | |
| CN115859855B (en) | River network sink water area extraction method based on still water balance principle | |
| CN114528360B (en) | Rapid flow distribution generation method based on flow velocity change | |
| DONGARDIVE | WATER RESOURCE PLANNING OF MICRO WATERSHED USING GEOSPATIAL TECHNIQUE. | |
| Garba et al. | Effect of slope and runoff trends on the hydrological response of River Kaduna |
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 |