CN112800631B - Urban waterlogging depth calculation method - Google Patents

Urban waterlogging depth calculation method Download PDF

Info

Publication number
CN112800631B
CN112800631B CN202110330233.8A CN202110330233A CN112800631B CN 112800631 B CN112800631 B CN 112800631B CN 202110330233 A CN202110330233 A CN 202110330233A CN 112800631 B CN112800631 B CN 112800631B
Authority
CN
China
Prior art keywords
value
depth
water
grid unit
target city
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
Application number
CN202110330233.8A
Other languages
Chinese (zh)
Other versions
CN112800631A (en
Inventor
薛丰昌
韩宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Chenxiang Space Information Technology Co ltd
Nanjing University of Information Science and Technology
Original Assignee
Nanjing Chenxiang Space Information Technology Co ltd
Nanjing University of Information Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Chenxiang Space Information Technology Co ltd, Nanjing University of Information Science and Technology filed Critical Nanjing Chenxiang Space Information Technology Co ltd
Priority to CN202110330233.8A priority Critical patent/CN112800631B/en
Publication of CN112800631A publication Critical patent/CN112800631A/en
Application granted granted Critical
Publication of CN112800631B publication Critical patent/CN112800631B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05—Geographic models

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Computer Graphics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method for calculating urban inland inundation water depth, which comprises the steps of initializing a lower limit value of target urban inland water depth, an upper limit value of ponding depth and a ponding depth value of a grid unit corresponding to a minimum elevation value in a target urban by using preset values; and calculating the estimated total ponding amount and the actual total ponding amount of the target area based on the ponding depth value of the grid unit corresponding to the minimum elevation value in the area, further acquiring a ponding amount difference value, and enabling the ponding amount difference value to be smaller than the ponding amount difference threshold value by optimizing the lower limit value of the ponding depth in the target city, the upper limit value of the ponding depth and the ponding depth value of the grid unit corresponding to the minimum elevation value in the target city, and taking the ponding depth value of the grid unit corresponding to the minimum elevation value in the target city meeting the condition as the ponding depth of the target city. The method overcomes the defects that the calculation of the depth of the accumulated water starts from zero, the step length of the calculation is fixed and the like in the traditional calculation method, which influence the calculation efficiency, and improves the calculation efficiency.

Description

Urban waterlogging depth calculation method
Technical Field
The invention relates to the technical field of urban waterlogging monitoring and early warning, in particular to a method for calculating the depth of urban waterlogging.
Background
The existing method for calculating the urban waterlogging depth is to obtain the maximum elevation value of a catchment area, calculate the possible ponding amount of the area based on the maximum elevation value, and approach the real ponding depth from bottom to top according to a fixed step length in the calculation process.
Disclosure of Invention
The purpose of the invention is as follows: the urban waterlogging water depth calculation method is high in calculation efficiency.
The technical scheme is as follows: the invention provides a method for calculating urban waterlogging water depth, which is characterized by being used for calculating the water depth in a target city, and comprising the following steps:
step 1: dividing the target city into I grid units based on the digital elevation data of the target city, and acquiring the elevation value and the side length r of each grid unit;
respectively using preset lower limit initial value of ponding depth, upper limit initial value of ponding depth and initial value of ponding depth of grid unit with minimum elevation value to initialize lower limit value d of ponding depth in target city0Upper limit value d of water accumulation depth1The accumulated water depth value d of the grid unit corresponding to the minimum elevation value in the target city;
then entering step 2;
step 2: acquiring the actual accumulated water total amount W of a target city; then entering step 3;
step 3; calculating the estimated total accumulated water amount W of the target city based on the accumulated water depth value d of the grid unit corresponding to the minimum elevation value in the target city1According to the formula C ═ W1-W | obtaining a total accumulated water difference C; then entering step 4;
and 4, step 4: judging whether the accumulated water total amount difference value C is smaller than the accumulated water total amount difference threshold value y, if so, entering the step 10, otherwise, entering the step 5;
and 5: judgment of d1Whether the water accumulation depth is equal to the initial value of the upper limit of the water accumulation depth preset in the step 1 or not is judged, and if yes, the step 6 is executed; otherwise, entering a step 8;
step 6: judgment of W1If the value is less than W, the value corresponding to d is assigned to d0Then, updating the value corresponding to d by 2 times of the value corresponding to d, and returning to the step 3; otherwise, entering step 7;
and 7: assigning the value corresponding to d1Returning to the step 3;
and 8: judgment of W1If the value is less than W, the value corresponding to d is assigned to d0Then apply
Figure GDA0003082199300000021
Updating the value corresponding to d, and returning to the step 3; otherwise, entering step 9;
and step 9: assigning the value corresponding to d1Then apply
Figure GDA0003082199300000022
Updating the value corresponding to d, and returning to the step 3;
step 10: and taking the ponding depth value d of the grid unit with the minimum elevation value in the target city as the ponding depth in the target city.
As a preferred aspect of the present invention, after step 10, the method further comprises:
step 11: according to the formula hmin+d-hiAcquiring ponding depth d 'of each grid unit i'iWherein h isiIs the elevation value, h, of grid cell iminThe minimum elevation value in the elevation values of each grid unit.
As a preferable scheme of the invention, in step 3, the estimated accumulated water total amount W of the target city is obtained according to the following method1：
Step 3.1: respectively obtaining the water accumulation depth of each grid unit, and further obtaining the sum Q of the water accumulation depths of each grid unit; then entering step 3.2;
step 3.2: according to the formula W1＝r2Obtaining estimated accumulated water total amount W of target city by multiplying Q1(ii) a Where r is the side length of a single grid cell.
As a preferred aspect of the present invention, in step 3.1, the depth of water accumulation of each grid cell is obtained according to the following method:
sequentially performing steps 3.1.1 to 3.1.2 to obtain a water accumulation depth d 'of each grid unit I for each grid unit I, I ═ 1,2.. I'i；
Step 3.1.1: calculate hmin+d-hiA value of (d); wherein h isiIs the elevation value, h, of grid cell iminThe minimum elevation value in the elevation values of each grid unit;
step 3.1.2: judgment hmin+d-hiIf the value of (b) is greater than or equal to 0, if so, h ismin+d-hiValue of (d) as Water depth d 'of grid Unit i'i(ii) a Otherwise, taking accumulated water depth d 'of grid unit i'i＝0。
As a preferable aspect of the present invention, in step 2, according to the following formula:
W＝r2×I×m
acquiring the actual accumulated water total amount W of a target city;
wherein r is the side length of a single grid unit, and m is the precipitation of the target city in a preset time period.
As a preferable scheme of the present invention, in step 1, the preset lower limit of the water accumulation depth is 0mm, the preset upper limit of the water accumulation depth is 0mm, and the water accumulation depth of the grid unit with the minimum preset elevation value is 100 mm.
As a preferable scheme of the present invention, in step 4, a value of the total ponding amount difference threshold y is one hundred thousand times of an actual total ponding amount W of the target city.
Has the advantages that: compared with the prior art, the method provided by the invention initializes the ponding depth value d of the grid unit corresponding to the minimum elevation value in the preset area with the ponding depth initial value of the grid unit with the minimum elevation value, calculates from the preset initial value upwards, and adjusts the calculation step length according to the calculation result in the calculation process.
Drawings
FIG. 1 is a block diagram of a flow chart of a method for calculating urban waterlogging water depth according to an embodiment of the present invention;
FIG. 2 is a schematic diagram illustrating an elevation distribution of grid cells according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of water accumulation of each grid unit when the water accumulation depth of the high-range lowest point grid is 100mm according to the embodiment of the invention;
fig. 4 is a schematic diagram of water accumulation depth results for each grid cell provided in accordance with an embodiment of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
Referring to fig. 1, the method provided by the present invention comprises the following steps:
step 1: dividing the target city into I grid units based on the digital elevation data of the target city, and acquiring the elevation value of each grid unit; respectively using preset lower limit initial value of ponding depth, upper limit initial value of ponding depth and initial value of ponding depth of grid unit with minimum elevation value to initialize lower limit value d of ponding depth in target city0Upper limit value d of water accumulation depth1The accumulated water depth value d of the grid unit corresponding to the minimum elevation value in the target city; then entering step 2;
step 2: according to the following formula:
W＝r2×I×m
and acquiring the actual accumulated water total W of the target city.
Wherein r is the side length of a single grid unit, and m is the precipitation of the target city in a preset time period.
Then step 3 is entered.
Step 3; calculating the estimated total accumulated water amount W of the target city based on the accumulated water depth value d of the grid unit corresponding to the minimum elevation value in the target city1(ii) a According to the formula C ═ W1-W | obtaining a total accumulated water difference C; then step 4 is entered.
Specifically, the estimated total accumulated water amount W of the target city is obtained according to the following method1：
Step 3.1: respectively obtaining the water accumulation depth of each grid unit, and further obtaining the sum Q of the water accumulation depths of each grid unit; then entering step 3.2;
specifically, the water accumulation depth of each grid unit is obtained according to the following method:
sequentially performing steps 3.1.1 to 3.1.2 to obtain a water accumulation depth d 'of each grid unit I for each grid unit I, I ═ 1,2.. I'i；
Step 3.1.1: calculate hmin+d-hiA value of (d); wherein h isiIs the elevation value, h, of grid cell iminThe minimum elevation value in the elevation values of each grid unit;
step 3.1.2: judgment hmin+d-hiIf the value of (b) is greater than or equal to 0, if so, h ismin+d-hiValue of (d) as Water depth d 'of grid Unit i'i(ii) a Otherwise, taking accumulated water depth d 'of grid unit i'i＝0。
And after the water accumulation depth of each grid unit is obtained according to the method, summing the water accumulation depths of the grid units to further obtain the sum Q of the water accumulation depths of the grid units.
Step 3.2: according to the formula W1＝r2Obtaining estimated accumulated water total amount W of target city by multiplying Q1(ii) a Where r is the side length of a single grid element.
And 4, step 4: and (5) judging whether the accumulated water total amount difference value C is smaller than the accumulated water total amount difference threshold value y, if so, entering the step 10, and otherwise, entering the step 5.
And 5: judgment of d1Whether the water accumulation depth is equal to the initial value of the upper limit of the water accumulation depth preset in the step 1 or not is judged, and if yes, the step 6 is executed; otherwise, go to step 8.
Step 6: judgment of W1If the value is less than W, the value corresponding to d is assigned to d0Then, updating the value corresponding to d by 2 times of the value corresponding to d, and returning to the step 3; otherwise step 7 is entered.
And 7: assigning the value corresponding to d1And returning to the step 3.
And 8: judgment of W1If the value is less than W, the value corresponding to d is assigned to d0Then apply
Figure GDA0003082199300000041
Updating the value corresponding to d, and returning to the step 3; otherwise step 9 is entered.
And step 9: assigning the value corresponding to d1Then apply
Figure GDA0003082199300000051
And updating the value corresponding to the d, and returning to the step 3.
Step 10: and taking the water accumulation depth value d of the grid unit with the minimum elevation value as the water accumulation depth in the target city.
In one embodiment, the target city is divided into 5 rows and 5 columns of 25 square grid cells, each grid cell has a side length of 100 meters, and the elevation value of each grid cell is shown in fig. 2. In the area corresponding to the grid cell, the rainfall lasted for 3 hours, and the total rainfall was 40.86mm within the 3 hours.
For the target city, according to W ═ r2Multiplying by I multiplied by m, and calculating to obtain the actual accumulated water total amount W of the target city, which is 10215m3。
Determining a total accumulated water quantity difference threshold value y: taking the accumulated water total quantity difference threshold value y as one hundred thousand of the actual accumulated water total quantity W of the target city, namely that y is approximately equal to 0.1022m3。
Traversing the elevation values of each grid unit to obtain the minimum value h of the elevation valuesmin＝7.0m。
Initializing lower limit value d of depth of water in target city0Upper limit value d of water accumulation depth1The initial values are respectively: d0＝0mm、d1And (4) initializing the accumulated water depth value d of the grid unit corresponding to the minimum elevation value in the target city, wherein d is 100 mm.
When the water accumulation depth d of the grid unit corresponding to the minimum elevation value in the target city is 100mm, the obtained water accumulation depth result of each grid unit is shown in fig. 3.
For the target city, the data of the relevant parameters in the loop process of step 1 to step 10 are shown in table 1.
TABLE 1
Figure GDA0003082199300000052
Figure GDA0003082199300000061
For the target city, the water accumulation depth of the grid unit with the smallest elevation value in the target city obtained in steps 1 to 10 is 280.375mm, and the water accumulation depth of each grid unit is shown in fig. 4.
The method provided by the invention initializes the ponding depth value d of the grid unit corresponding to the minimum elevation value in the preset area with the ponding depth initial value of the grid unit with the minimum elevation value, and calculates upwards from the preset initial value, and adjusts the calculation step length according to the calculation result in the calculation process, thereby eliminating the defects of influence on the calculation efficiency caused by the approximation of the ponding depth from zero to top, the fixation of the calculated step length and the like in the traditional calculation method, and greatly improving the calculation efficiency.
The above description is only a preferred embodiment of the present invention, and it will be apparent to those skilled in the art that various modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be considered as the protection scope of the present invention.

Claims (5)

1. A method for calculating urban waterlogging water depth is used for calculating the water depth in a target city, and comprises the following steps:
step 1: dividing the target city into I grid units based on the digital elevation data of the target city, and acquiring the elevation value and the side length r of each grid unit;
respectively using the products of the preset water accumulation depth lower limit initial value, water accumulation depth upper limit initial value and grid unit with the minimum elevation valueInitializing lower limit value d of depth of water in target city by using initial value of depth of water0Upper limit value d of water accumulation depth1The accumulated water depth value d of the grid unit corresponding to the minimum elevation value in the target city;
then entering step 2;
step 2: acquiring the actual accumulated water total amount W of a target city; then entering step 3;
step 3; calculating the estimated total water volume W of the target city based on the water depth value d of the grid unit corresponding to the minimum elevation value in the target city according to the following steps 3.1 to 3.21According to the formula C ═ W1-W | obtaining a total accumulated water difference C; then entering step 4;
step 3.1: sequentially performing steps 3.1.1 to 3.1.2 to obtain a water accumulation depth d 'of each grid unit I for each grid unit I, I ═ 1,2.. I'iFurther acquiring the sum Q of the water accumulation depths of each grid unit; then entering step 3.2;
step 3.1.1: calculate hmin+d-hiA value of (d); wherein h isiIs the elevation value, h, of grid cell iminThe minimum elevation value in the elevation values of each grid unit;
step 3.1.2: judgment hmin+d-hiIf the value of (b) is greater than or equal to 0, if so, h ismin+d-hiValue of (d) as Water depth d 'of grid Unit i'i(ii) a Otherwise, taking accumulated water depth d 'of grid unit i'i＝0；
Step 3.2: according to the formula W1＝r2Obtaining estimated accumulated water total amount W of target city by multiplying Q1(ii) a Wherein r is the side length of each grid unit;
and 4, step 4: judging whether the accumulated water total amount difference value C is smaller than the accumulated water total amount difference threshold value y, if so, entering the step 10, otherwise, entering the step 5;
and 5: judgment of d1Whether the water accumulation depth is equal to the initial value of the upper limit of the water accumulation depth preset in the step 1 or not is judged, and if yes, the step 6 is executed; otherwise, entering a step 8;
step 6: judgment of W1If the value is less than W, the value corresponding to d is assigned to d0Then using d pairsUpdating the value corresponding to d by 2 times of the corresponding value, and returning to the step 3; otherwise, entering step 7;
and 7: assigning the value corresponding to d1Returning to the step 3;
and 8: judgment of W1If the value is less than W, the value corresponding to d is assigned to d0Then apply
Figure FDA0003082199290000021
Updating the value corresponding to d, and returning to the step 3; otherwise, entering step 9;
and step 9: assigning the value corresponding to d1Then apply
Figure FDA0003082199290000022
Updating the value corresponding to d, and returning to the step 3;
step 10: and taking the ponding depth value d of the grid unit with the minimum elevation value in the target city as the ponding depth in the target city.
2. The method for urban waterlogging depth calculation according to claim 1, wherein after step 10, the method further comprises:
step 11: according to the formula hmin+d-hiAcquiring ponding depth d 'of each grid unit i'iWherein h isiIs the elevation value, h, of grid cell iminThe minimum elevation value in the elevation values of each grid unit.
3. The method for urban waterlogging depth calculation according to claim 1, wherein in step 2, according to the following formula:
W＝r2×I×m
acquiring the actual accumulated water total amount W of a target city;
wherein r is the side length of a single grid unit, and m is the precipitation of the target city in a preset time period.
4. The method for calculating urban waterlogging water depth according to claim 1, wherein in step 1, the preset lower water depth limit is initially 0mm, the preset upper water depth limit is initially 0mm, and the preset grid unit with the smallest elevation value has a water depth of 100 mm.
5. The method for calculating urban waterlogging water depth according to claim 1, wherein in step 4, the value of the total water volume difference threshold y is one hundred thousand of the actual total water volume W of the target city.
CN202110330233.8A 2021-03-29 2021-03-29 Urban waterlogging depth calculation method Active CN112800631B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110330233.8A CN112800631B (en) 2021-03-29 2021-03-29 Urban waterlogging depth calculation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110330233.8A CN112800631B (en) 2021-03-29 2021-03-29 Urban waterlogging depth calculation method

Publications (2)

Publication Number Publication Date
CN112800631A CN112800631A (en) 2021-05-14
CN112800631B true CN112800631B (en) 2021-07-13

Family

ID=75815873

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110330233.8A Active CN112800631B (en) 2021-03-29 2021-03-29 Urban waterlogging depth calculation method

Country Status (1)

Country Link
CN (1) CN112800631B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102902893A (en) * 2012-10-29 2013-01-30 南京信息工程大学 Method for calculating rainfall ponding depth of catchment area based on DEM (digital elevation model)
CN105760666A (en) * 2016-02-04 2016-07-13 南京信息工程大学 Catchment area critical rainfall calculating method based on DEM
CN106295123A (en) * 2016-07-27 2017-01-04 泰华智慧产业集团股份有限公司 The computational methods of multi-source non-single depression surface pond state

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111369059B (en) * 2020-03-05 2023-10-20 山东大学 Urban waterlogging rapid prediction method and system based on rain and flood simulation coupling model

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102902893A (en) * 2012-10-29 2013-01-30 南京信息工程大学 Method for calculating rainfall ponding depth of catchment area based on DEM (digital elevation model)
CN105760666A (en) * 2016-02-04 2016-07-13 南京信息工程大学 Catchment area critical rainfall calculating method based on DEM
CN106295123A (en) * 2016-07-27 2017-01-04 泰华智慧产业集团股份有限公司 The computational methods of multi-source non-single depression surface pond state

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Numerical Simulation of Urban Waterlogging Based on FloodArea Model;Fengchang Xue等;《Advances in Meteorology》;20161231(第2016期);第1-10页 *
基于SWMM模型的沧州市区暴雨积水深度模拟;贾艳梅等;《河北工程技术高等专科学校学报》;20160630(第2期);第15-20页 *
视频监控的城市内涝监测预警;薛丰昌等;《测绘科学》;20180831;第43卷(第8期);第50-55,61页 *

Also Published As

Publication number Publication date
CN112800631A (en) 2021-05-14

Similar Documents

Publication Publication Date Title
CN104537650B (en) The polar orientation decision method and device of a kind of aluminium electrolutic capacitor
CN113435630B (en) Basin hydrological forecasting method and system with self-adaptive runoff yield mode
CN113808193B (en) Light spot centroid positioning method based on blocking threshold
CN110132200B (en) Tailing pond dry beach dynamic monitoring method and system based on Beidou and video identification
CN111046773A (en) Method for judging water retention in pavement based on image technology
CN110442988B (en) Urban surface runoff flow direction calculation method and device based on cellular automaton
CN104134195B (en) Image restoration method based on geometric sparsity of blocks
CN112800631B (en) Urban waterlogging depth calculation method
CN110717900A (en) Pantograph abrasion detection method based on improved Canny edge detection algorithm
CN105389793A (en) Method for automatically identifying trend and width of fracture in object surface
CN115049163B (en) Photovoltaic region prediction method and device based on tabu search and deep stable learning
CN114219252B (en) Basin unit line analysis method based on SCE-UA algorithm
CN104091040A (en) Soil infiltrability calculation method
CN115660196A (en) Production management method and system for pharmaceutical workshop
CN110212592B (en) Thermal power generating unit load regulation maximum rate estimation method and system based on piecewise linear expression
CN109118721B (en) Method and device for calculating composite early warning index based on critical rainfall and early warning method
CN107871315A (en) A kind of video image motion detection method and device
CN113486556B (en) Improved efficient automatic history fitting method for oil and gas reservoir
CN116578931A (en) Urban waterlogging and ponding prediction method based on Bagging-XGB model
CN101923716B (en) Method for improving particle filter tracking effect
CN116628402A (en) Tailing pond flood regulating calculation method and device, electronic equipment and storage medium
CN112819234B (en) Flood forecasting method and system considering initial value correction
CN114155700A (en) Road surface slippery state monitoring method and system
CN111414676B (en) Method for evaluating performance of LID facility by storage tank model
CN111626564B (en) Method for rapidly generating open-drain sand-discharging scheduling scheme of reservoir camera

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