CN116611249A - Urban area rainwater infiltration amount calculation method - Google Patents
Urban area rainwater infiltration amount calculation method Download PDFInfo
- Publication number
- CN116611249A CN116611249A CN202310592499.9A CN202310592499A CN116611249A CN 116611249 A CN116611249 A CN 116611249A CN 202310592499 A CN202310592499 A CN 202310592499A CN 116611249 A CN116611249 A CN 116611249A
- Authority
- CN
- China
- Prior art keywords
- permeable
- grid
- permeable surface
- water
- infiltration
- 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
- 230000008595 infiltration Effects 0.000 title claims abstract description 103
- 238000001764 infiltration Methods 0.000 title claims abstract description 103
- 238000004364 calculation method Methods 0.000 title abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 188
- 230000005540 biological transmission Effects 0.000 claims abstract description 43
- 238000000034 method Methods 0.000 claims abstract description 34
- 230000011218 segmentation Effects 0.000 claims abstract description 4
- 239000002689 soil Substances 0.000 claims description 18
- 230000035699 permeability Effects 0.000 claims description 16
- 230000008020 evaporation Effects 0.000 claims description 15
- 238000001704 evaporation Methods 0.000 claims description 15
- 239000012466 permeate Substances 0.000 claims description 6
- 230000035515 penetration Effects 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 description 6
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- 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
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Mathematical Physics (AREA)
- Databases & Information Systems (AREA)
- Computational Mathematics (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Sewage (AREA)
Abstract
The invention discloses a calculation method of urban area rainwater infiltration amount, which comprises the steps of carrying out grid segmentation on an urban area to obtain a permeable grid set and a impermeable grid set; determining the flow direction of water flow of the grid; combining the permeable grids to obtain a permeable surface, and identifying the permeable surface as a storage type permeable surface or a transmission type permeable surface according to the water flow direction of each permeable grid in the permeable surface and the elevation relation between the permeable grids and surrounding grids, and obtaining the storage type permeable surface area and the transmission type permeable surface area; calculating the effective runoff of the stagnation type water permeable surface; calculating the rainwater infiltration amount of the stagnant water permeable surface; calculating the infiltration amount of the rainwater on the transmission type permeable surface; and (5) calculating the infiltration amount of rainwater in the urban area. According to the invention, based on the spatial pattern characteristics of the permeable surface and the impermeable surface of the urban area, the rain infiltration amount of the underlying surfaces of different types is calculated according to the rain infiltration principle and the rain infiltration process of the underlying surfaces of different types, so that the scientificity and the precision of the calculation of the rain infiltration amount of the urban area can be effectively improved.
Description
Technical Field
The invention belongs to the technical field of rainwater infiltration quantity quantification, and particularly relates to a calculation method of urban area rainwater infiltration quantity.
Background
The rainwater infiltration amount of urban areas is generally calculated by using the infiltration coefficient, the area of the permeable surface and the rainfall duration. The influence of the space pattern features of the underlying surface of the urban area on the runoff path and the infiltration process is large, the continuous infiltration of the rainwater runoff quantity accumulated on the low-lying water permeable surface and the sponge facilities after rainfall is considered, and the existing calculation means cannot meet the fine requirements.
Disclosure of Invention
The invention provides a calculation method of urban area rainwater infiltration amount, which aims at solving the problem of how to quantify urban area rainwater infiltration amount.
In order to achieve the above purpose, the technical scheme provided by the invention is as follows: the method for calculating the rainwater infiltration amount of the urban area comprises the following steps:
s1: grid segmentation is carried out on the urban area to obtain a permeable grid set and a impermeable grid set;
s2: determining the flow direction of water flow of the grid according to the elevation data of the grid;
s3: combining the permeable grids to obtain a permeable surface, and identifying the permeable surface as a storage type permeable surface or a transmission type permeable surface according to the water flow direction of each permeable grid in the permeable surface and the elevation relation between the permeable grids and surrounding grids, and obtaining the storage type permeable surface area and the transmission type permeable surface area;
s4: calculating the effective runoff of the stagnation type water permeable surface according to the area of the stagnation type water permeable surface;
s5: according to the effective runoff of the stagnant water permeable surface, calculating the rainwater infiltration quantity of the stagnant water permeable surface;
s6: according to the area of the transmission type permeable surface, calculating the infiltration amount of the rainwater of the transmission type permeable surface;
s7: and calculating the rainwater infiltration amount of the urban area according to the rainwater infiltration amount of the stagnant water permeable surface and the rainwater infiltration amount of the transmission type permeable surface.
Based on the technical scheme, the invention can also be improved as follows:
further, the step S2 specifically includes:
s201: calculating elevation data of all grids;
s202: taking any grid as a central grid, comparing elevation data of all adjacent grids of the central grid with the elevation data of all adjacent grids, and connecting centers of the adjacent grids with the largest elevation data difference as water flow directions of the central grid;
s203: step S202 is repeated until the flow directions of all grids are determined.
The invention adopts the further technical proposal and has the beneficial effects that: based on the basic principle of water flowing downwards, the flow direction of grid water flow can be reasonably determined.
Further, the step S3 specifically includes:
s301: sequentially distinguishing grids in the surrounding direction by taking any permeable grid as a center;
if the grid is a permeable grid, combining the central permeable grid and the grid to form a permeable surface;
otherwise, not merging;
s302: sequentially identifying all grids around the water permeable surface by taking the water permeable surface obtained in the step S301 as a center;
if the grid is a permeable grid, combining the grid with the permeable surface and updating the permeable surface;
otherwise, not merging;
s303: repeating step S302 with the updated permeable surface as the center until all the permeable grids are combined to form the permeable surface;
s304: determining the height Cheng Guanji of any permeable grid in the permeable surface and the grids around the permeable surface;
if the elevation of the permeable grid is higher than or equal to that of the surrounding grids, the permeable grid is divided from the permeable surface, and the permeable surface is updated;
otherwise, not dividing;
s305: repeating the step S304 until all the permeable grids in the permeable surface are too high Cheng Guanji compared with the grids around the permeable surface;
s306: repeating the steps S301 and S302 on the water permeable grid segmented in the step S304 to obtain a new water permeable surface;
s307: judging the type of the permeable surface according to the water flow direction of each permeable grid in the permeable surface and the elevation relation between the permeable grids and surrounding grids, wherein the type of the permeable surface comprises a stagnation type permeable surface and a transmission type permeable surface;
s308: repeating step S307 until all the water permeable surfaces are marked as the storage type water permeable surface or the transmission type water permeable surface;
s309: and counting the identified water-permeable area of the stagnation type and the transmission type.
The invention adopts the further technical proposal and has the beneficial effects that: the water flow direction through single grid can't confirm that the permeable surface grid is stagnating and holds type permeable surface grid or transmission type permeable surface grid, and this step is through combining the permeable surface grid, can be with the permeable surface grid of hydrologic connection relation to merge as a whole, according to the holistic positional relationship of permeable surface and impervious surface, confirms the type of permeable surface. According to the rainwater infiltration principle and process of different types of underlying surfaces, the rainwater infiltration amount can be calculated by adopting different methods.
Further, step S4 specifically includes:
s401: according to the area of the single permeable grid in the stagnant water permeable surface, calculating the effective water storage volume V of the single permeable grid in the stagnant water permeable surface:
V＝H×a
wherein H is the effective water storage height of the permeable grids, and a is the area of a single permeable grid in the stagnant water permeable surface;
s402: statistics of stagnancyCalculating the effective water storage volume V of the stagnated water permeable surface by using the number n of all the water permeable grids in the water permeable surface z ：
V Z ＝nV
S403: determining the total runoff amount R t ：
Wherein P is rainfall, A pi For the area of the grid i in the permeate water,as the runoff coefficient of the grid i in the permeate water,is the area of the grid j in the impermeable water, +.>The runoff coefficient of the impermeable middle grid j;
s404: comparing the total runoff amount R t Effective water storage volume V of water permeable surface z Is of a size of (2);
if the total runoff amount R t Is greater than the effective water storage volume V z The effective runoff of the stagnated water permeable surface is the effective water storage volume V z ；
Otherwise, the effective runoff amount of the stagnation type water permeable surface is the total runoff amount R t 。
The invention adopts the further technical proposal and has the beneficial effects that: the influence of the water collecting range and rainfall factor on the rainwater runoff collected by the stagnation and storage type permeable surface grid is comprehensively considered, and the effective runoff of the stagnation and storage type permeable surface, namely the stagnation and storage amount, can be accurately calculated.
Further, step S5 specifically includes:
s501: when the effective runoff of the stagnant water permeable surface is the effective water storage volume V z And setting the infiltration time to be the same as the evaporation time, and obtaining the infiltration time of the permeable grid:
μ i ×t i +ε i ×t i ＝H i
wherein mu is i The soil permeability, t, of the permeable grid i i For the infiltration time, ε, of the water-permeable grid i i For evaporation rate of water permeable grid i, H i The effective water storage height of the permeable grid i;
according to the area of the stagnant water permeable surface and the infiltration time of the permeable grid, calculating the infiltration quantity I of the stagnant water permeable surface z ：
In the method, in the process of the invention,is the area mu of the stagnation type permeable surface grid i i The soil permeability, t, of the permeable grid i i The infiltration time of the permeable grid i;
s502: when the effective runoff amount of the stagnant water permeable surface is the total runoff amount R of the catchment range t When the water level line is measured, the height h from the nearest permeable grid to the water level line is obtained:
A×[h+(h+E 1 )+(h+E 2 )...+(h+E i )]＝R t
wherein A is the area of a single permeable grid of a stagnant water permeable surface, E 1 、E 2 …E i Is the corresponding grid elevation difference;
according to the height from the nearest permeable grid to the water surface line, setting the infiltration time to be the same as the evaporation time, and solving the infiltration time t of the nearest permeable grid from the water surface line:
μ×t+ε×t＝h
wherein mu is the soil permeability of the nearest permeable grid to the water surface line, epsilon is the evaporation rate of the nearest permeable grid to the water surface line;
sequentially calculating the infiltration time t of other permeable grids i i '：
μ i '×t i '+ε i '×t i '＝h+E 1 +E 2 +...E i '
Wherein mu is the soil permeability rate of the permeable grid i, t i ' is the downpermeation time of the permeable grid i, ε is the evaporation rate of the permeable grid i, E 1 、E 2 …E i ' is the corresponding grid elevation difference;
according to the area of the stagnant water permeable surface and the infiltration time of the permeable grid, calculating the infiltration quantity I of the stagnant water permeable surface z '：
In the method, in the process of the invention,is the area mu of the stagnation type permeable surface grid i i ' soil penetration rate, t, of permeable grid i i ' is the time of infiltration of the water permeable grid i.
The invention adopts the further technical proposal and has the beneficial effects that: based on the relation between the total runoff amount in the catchment range and the effective water storage volume of the stagnation water permeable surface, the stagnation water storage levels under different conditions are respectively determined, the height from each water permeable surface grid to the water surface line is calculated, and then the rainwater infiltration amount of the stagnation water permeable surface is obtained.
Further, according to the area of the transmission type permeable surface, the rainwater infiltration quantity I of the transmission type permeable surface is calculated c ：
In the method, in the process of the invention,is the area mu of the transmission type permeable surface grid i i The soil permeability coefficient, t, of the permeable grid i i Is the infiltration time of the water permeable grid i.
Further, according to the stagnant water permeable surface rainwater infiltration amount and the transmission type permeable surface rainwater infiltration amount, calculating the urban area rainwater infiltration amount:
urban area rainwater infiltration amount=i z +I c
Wherein I is z I is the infiltration amount of rainwater on the stagnant water-permeable surface c Is the infiltration amount of the rainwater on the transmission type permeable surface.
The beneficial effects of the invention are as follows: the method is used for calculating the rainwater infiltration amount of the urban area, so that the calculation accuracy can be effectively improved, and meanwhile, the method is beneficial to quantitatively analyzing the hydrologic mechanism of the water permeability surface of the urban area and promotes the development of related researches.
Drawings
FIG. 1 is a method and a flow for calculating the infiltration amount of rainwater in urban areas;
FIG. 2 is a cross-sectional view of a water permeable surface;
FIG. 3 is a diagram showing an example of water permeable surface type determination;
FIG. 4 is a graph showing an example of calculation of the height of the water stored in the grid when the effective runoff of the water permeable surface is the effective water storage volume;
FIG. 5 is a graph showing an example of calculation of the height of the water stored in the grid when the effective runoff amount of the water permeable surface is the total runoff amount in the catchment area.
Detailed Description
The principles and features of the present invention are described below with reference to the drawings so as to facilitate understanding of the present invention by those skilled in the art, and the examples are provided for explanation of the present invention only and are not intended to limit the scope of the present invention. The specific conditions are not noted in the examples and are carried out according to conventional conditions or conditions recommended by the manufacturer.
Examples:
a calculation method of urban area rainwater infiltration amount is shown in figure 1, and comprises the following steps:
s1: grid segmentation is carried out on the urban area to obtain a permeable grid set and a impermeable grid set;
s2: according to the elevation data of the grids, determining the flow direction of the water flow of the grids, wherein the specific steps are as follows:
s201: calculating elevation data of all grids;
s202: taking any grid as a central grid, comparing elevation data of all adjacent grids of the central grid with the elevation data of all adjacent grids, and connecting centers of the adjacent grids with the largest elevation data difference as water flow directions of the central grid;
s203: repeating the step S202 until the water flow directions of all grids are determined;
s3: combining the permeable grids to obtain a permeable surface, identifying the permeable surface as a stagnation type permeable surface or a transmission type permeable surface according to the water flow direction of each permeable grid in the permeable surface and the elevation relation between the permeable grids and surrounding grids, and obtaining the stagnation type permeable surface area and the transmission type permeable surface area, wherein the specific steps are as follows:
s301: sequentially distinguishing grids in the surrounding direction by taking any permeable grid as a center;
if the grid is a permeable grid, combining the central permeable grid and the grid to form a permeable surface;
otherwise, not merging;
s302: sequentially identifying all grids around the water permeable surface by taking the water permeable surface obtained in the step S301 as a center;
if the grid is a permeable grid, combining the grid with the permeable surface and updating the permeable surface;
otherwise, not merging;
s303: repeating step S302 with the updated permeable surface as the center until all the permeable grids are combined to form the permeable surface;
s304: determining the height Cheng Guanji of any permeable grid in the permeable surface and the grids around the permeable surface;
if the elevation of the permeable grid is higher than or equal to that of the surrounding grids, the permeable grid is divided from the permeable surface, and the permeable surface is updated;
otherwise, not dividing;
s305: repeating the step S304 until all the permeable grids in the permeable surface are too high Cheng Guanji compared with the grids around the permeable surface;
s306: repeating the steps S301 and S302 on the water permeable grid segmented in the step S304 to obtain a new water permeable surface, as shown in fig. 2;
s307: judging the type of the permeable surface according to the water flow direction of each permeable grid in the permeable surface and the elevation relation between the permeable grids and surrounding grids, wherein the type of the permeable surface comprises a stagnation type permeable surface and a transmission type permeable surface;
as shown in fig. 3, if the water flow direction of each permeable grid in the permeable surface is all inward and all permeable grids are lower in elevation than surrounding grids, the permeable surface is a stagnant permeable surface; or the water flow directions of the water permeable grids in the water permeable surface are all inward, but parts with the heights equal to or higher than those of surrounding grids exist, the parts with the heights lower than those of the surrounding grids are identified as storage type water permeable surfaces, and the parts with the heights equal to and higher than those of the surrounding grids are identified as transmission type water permeable surfaces;
if each permeable grid in the permeable surface has an outward water flow direction, taking the height of the impermeable surface with the outward water flow direction as a boundary, identifying the impermeable surface as a stagnation type permeable surface at a height lower than the boundary, and identifying the impermeable surface as a transmission type permeable surface at a height equal to or higher than the boundary;
s308: repeating step S307 until all the water permeable surfaces are marked as the storage type water permeable surface or the transmission type water permeable surface;
s309: counting the identified water-permeable area of the stagnation type and the transmission type;
s4: according to the area of the stagnation type permeable surface, calculating the effective runoff of the stagnation type permeable surface, wherein the effective runoff comprises the following specific steps:
s401: according to the area of the single permeable grid in the stagnant water permeable surface, calculating the effective water storage volume V of the single permeable grid in the stagnant water permeable surface:
V＝H×a
wherein H is the effective water storage height of the permeable grids, and a is the area of a single permeable grid in the stagnant water permeable surface; the effective water storage height H of the permeable grid is the difference between the permeable grid with the lowest elevation and the impermeable grid with the lowest elevation at the periphery, if the impermeable surface is not arranged at the periphery of the permeable grid, the outward search is continued until the impermeable surface with the lowest elevation at the periphery is searched;
s402: statistics of stagnancyCalculating the effective water storage volume V of the stagnated water permeable surface by using the number n of all the water permeable grids in the water permeable surface z ：
V Z ＝nV
S403: determining the total runoff amount R t ：
Wherein P is the rainfall,for the area of the grid i in the permeate water +.>For the runoff coefficient of grid i in the permeate water, +.>Is the area of the grid j in the impermeable water, +.>The runoff coefficient of the impermeable middle grid j;
s404: comparing the total runoff amount R t Effective water storage volume V of water permeable surface z Is of a size of (2);
if the total runoff amount R t Is greater than the effective water storage volume V z The effective runoff of the stagnated water permeable surface is the effective water storage volume V z ；
Otherwise, the effective runoff amount of the stagnation type water permeable surface is the total runoff amount R t ；
S5: according to the effective runoff of the stagnant water permeable surface, calculating the rainwater infiltration quantity of the stagnant water permeable surface, wherein the method comprises the following specific steps of:
s501: when the effective runoff of the stagnant water permeable surface is the effective water storage volume V z In this case, as shown in fig. 4, the infiltration time was set to be the same as the evaporation time, and the infiltration time of the water permeable grid was obtained:
μ i ×t i +ε i ×t i ＝H i
wherein mu is i The soil permeability, t, of the permeable grid i i For the infiltration time, ε, of the water-permeable grid i i For evaporation rate of water permeable grid i, H i The effective water storage height of the permeable grid i;
according to the area of the stagnant water permeable surface and the infiltration time of the permeable grid, calculating the infiltration quantity I of the stagnant water permeable surface z ：
In the method, in the process of the invention,is the area mu of the stagnation type permeable surface grid i i The soil permeability, t, of the permeable grid i i The infiltration time of the permeable grid i;
s502: when the effective runoff amount of the stagnant water permeable surface is the total runoff amount R of the catchment range t At this time, as shown in fig. 5, the height h from the nearest water permeable grid to the water surface line is obtained:
A×[h+(h+E 1 )+(h+E 2 )...+(h+E i )]＝R t
wherein A is the area of a single permeable grid of a stagnant water permeable surface, E 1 、E 2 …E i Is the corresponding grid elevation difference;
according to the height from the nearest permeable grid to the water surface line, setting the infiltration time to be the same as the evaporation time, and solving the infiltration time t of the nearest permeable grid from the water surface line:
μ×t+ε×t＝h
wherein mu is the soil permeability of the nearest permeable grid to the water surface line, epsilon is the evaporation rate of the nearest permeable grid to the water surface line;
sequentially calculating the infiltration time t of other permeable grids i i '：
μ i '×t i '+ε i '×t i '＝h+E 1 +E 2 +...E i '
Wherein mu is the soil permeability rate of the permeable grid i, t i ' is the downpermeation time of the permeable grid i, ε is the evaporation rate of the permeable grid i, E 1 、E 2 …E i ' is the corresponding grid elevation difference;
according to the area of the stagnant water permeable surface and the infiltration time of the permeable grid, calculating the infiltration quantity I of the stagnant water permeable surface z '：
In the method, in the process of the invention,is the area mu of the stagnation type permeable surface grid i i ' soil penetration rate, t, of permeable grid i i ' is the infiltration time of the water permeable grid i;
s6: according to the area of the transmission type permeable surface, calculating the infiltration quantity I of the rainwater of the transmission type permeable surface c ：
In the method, in the process of the invention,is the area mu of the transmission type permeable surface grid i i The soil permeability coefficient, t, of the permeable grid i i For the infiltration time of the permeable grid i, t for the transmission type permeable surface i Taking rainfall duration;
s7: according to the stagnant water permeable surface rainwater infiltration amount and the transmission type permeable surface rainwater infiltration amount, calculating the urban area rainwater infiltration amount:
urban area rainwater infiltration amount=i z +I c
Wherein I is z I is the infiltration amount of rainwater on the stagnant water-permeable surface c Is the infiltration amount of the rainwater on the transmission type permeable surface.
The foregoing description of the preferred embodiments of the invention is not intended to limit the invention to the precise form disclosed, and any such modifications, equivalents, and alternatives falling within the spirit and scope of the invention are intended to be included within the scope of the invention.
Claims (7)
1. The method for calculating the rainwater infiltration amount of the urban area is characterized by comprising the following steps of:
s1: grid segmentation is carried out on the urban area to obtain a permeable grid set and a impermeable grid set;
s2: determining the flow direction of water flow of the grid according to the elevation data of the grid;
s3: combining the permeable grids to obtain a permeable surface, and identifying the permeable surface as a storage type permeable surface or a transmission type permeable surface according to the water flow direction of each permeable grid in the permeable surface and the elevation relation between the permeable grids and surrounding grids, and obtaining the storage type permeable surface area and the transmission type permeable surface area;
s4: calculating the effective runoff of the stagnation type water permeable surface according to the area of the stagnation type water permeable surface;
s5: according to the effective runoff of the stagnant water permeable surface, calculating the rainwater infiltration quantity of the stagnant water permeable surface;
s6: according to the area of the transmission type permeable surface, calculating the infiltration amount of the rainwater of the transmission type permeable surface;
s7: and calculating the rainwater infiltration amount of the urban area according to the rainwater infiltration amount of the stagnant water permeable surface and the rainwater infiltration amount of the transmission type permeable surface.
2. The method for calculating the rainwater infiltration amount of the urban area according to claim 1, wherein the step S2 is specifically:
s201: calculating elevation data of all grids;
s202: taking any grid as a central grid, comparing elevation data of all adjacent grids of the central grid with the elevation data of all adjacent grids, and connecting centers of the adjacent grids with the largest elevation data difference as water flow directions of the central grid;
s203: step S202 is repeated until the flow directions of all grids are determined.
3. The method for calculating the rainwater infiltration capacity of the urban area according to claim 1, wherein the step S3 is specifically:
s301: sequentially distinguishing grids in the surrounding direction by taking any permeable grid as a center;
if the grid is a permeable grid, combining the central permeable grid and the grid to form a permeable surface;
otherwise, not merging;
s302: sequentially identifying all grids around the water permeable surface by taking the water permeable surface obtained in the step S301 as a center;
if the grid is a permeable grid, combining the grid with the permeable surface and updating the permeable surface;
otherwise, not merging;
s303: repeating step S302 with the updated permeable surface as the center until all the permeable grids are combined to form the permeable surface;
s304: determining the height Cheng Guanji of any permeable grid in the permeable surface and the grids around the permeable surface;
if the elevation of the permeable grid is higher than or equal to that of the surrounding grids, the permeable grid is divided from the permeable surface, and the permeable surface is updated;
otherwise, not dividing;
s305: repeating the step S304 until all the permeable grids in the permeable surface are too high Cheng Guanji compared with the grids around the permeable surface;
s306: repeating the steps S301 and S302 on the water permeable grid segmented in the step S304 to obtain a new water permeable surface;
s307: judging the type of the permeable surface according to the water flow direction of each permeable grid in the permeable surface and the elevation relation between the permeable grids and surrounding grids, wherein the type of the permeable surface comprises a stagnation type permeable surface and a transmission type permeable surface;
s308: repeating step S307 until all the water permeable surfaces are marked as the storage type water permeable surface or the transmission type water permeable surface;
s309: and counting the identified water-permeable area of the stagnation type and the transmission type.
4. The method for calculating the rainwater infiltration capacity of the urban area according to claim 1, wherein the step S4 is specifically:
s401: according to the area of the single permeable grid in the stagnant water permeable surface, calculating the effective water storage volume V of the single permeable grid in the stagnant water permeable surface:
V＝H×a
wherein H is the effective water storage height of the permeable grids, and a is the area of a single permeable grid in the stagnant water permeable surface;
s402: counting the number n of all permeable grids in the stagnant water permeable surface, and calculating the effective water storage volume V of the stagnant water permeable surface z ：
V Z ＝nV
S403: determining the total runoff amount R t ：
Wherein P is the rainfall,for the area of the grid i in the permeate water +.>For the runoff coefficient of grid i in the permeate water, +.>Is the area of the grid j in the impermeable water, +.>The runoff coefficient of the impermeable middle grid j;
s404: comparing the total runoff amount R t Effective water storage volume V of water permeable surface z Is of a size of (2);
if the total runoff amount R t Is greater than the effective water storage volume V z The effective runoff of the stagnated water permeable surface is the effective water storage volume V z ；
Otherwise, the effective runoff amount of the stagnation type water permeable surface is the total runoff amount R t 。
5. The method for calculating the rainwater infiltration capacity of the urban area according to claim 1, wherein the step S5 is specifically:
s501: when the effective runoff of the stagnant water permeable surface is the effective water storage volume V z And setting the infiltration time to be the same as the evaporation time, and obtaining the infiltration time of the permeable grid:
μ i ×t i +ε i ×t i ＝H i
wherein mu is i The soil permeability, t, of the permeable grid i i For the infiltration time, ε, of the water-permeable grid i i For evaporation rate of water permeable grid i, H i The effective water storage height of the permeable grid i;
according to the area of the stagnant water permeable surface and the infiltration time of the permeable grid, calculating the infiltration quantity I of the stagnant water permeable surface z ：
In the method, in the process of the invention,is the area mu of the stagnation type permeable surface grid i i The soil permeability, t, of the permeable grid i i The infiltration time of the permeable grid i;
s502: when the effective runoff amount of the stagnant water permeable surface is the total runoff amount R of the catchment range t When the water level line is measured, the height h from the nearest permeable grid to the water level line is obtained:
A×[h+(h+E 1 )+(h+E 2 )...+(h+E i )]＝R t
wherein A is the area of a single permeable grid of a stagnant water permeable surface, E 1 、E 2 …E i Is the corresponding grid elevation difference;
according to the height from the nearest permeable grid to the water surface line, setting the infiltration time to be the same as the evaporation time, and solving the infiltration time t of the nearest permeable grid from the water surface line:
μ×t+ε×t＝h
wherein mu is the soil permeability of the nearest permeable grid to the water surface line, epsilon is the evaporation rate of the nearest permeable grid to the water surface line;
sequentially calculating the infiltration time t of other permeable grids i i '：
μ i '×t i '+ε i '×t i '＝h+E 1 +E 2 +...E i '
Wherein mu is the soil permeability rate of the permeable grid i, t i ' is the downpermeation time of the permeable grid i, ε is the evaporation rate of the permeable grid i, E 1 、E 2 …E i ' is the corresponding grid elevation difference;
according to the area of the stagnant water permeable surface and the infiltration time of the permeable grid, calculating the infiltration quantity I of the stagnant water permeable surface z '：
In the method, in the process of the invention,is the area mu of the stagnation type permeable surface grid i i ' soil penetration rate, t, of permeable grid i i ' is the time of infiltration of the water permeable grid i.
6. The method for calculating the rainwater infiltration amount of the urban area according to claim 1, wherein: according to the area of the transmission type permeable surface, calculating the infiltration quantity I of the rainwater of the transmission type permeable surface c ：
In the method, in the process of the invention,is the area mu of the transmission type permeable surface grid i i The soil permeability coefficient, t, of the permeable grid i i Is the infiltration time of the water permeable grid i.
7. The method for calculating the rainwater infiltration amount of the urban area according to claim 1, wherein: according to the stagnant water permeable surface rainwater infiltration amount and the transmission type permeable surface rainwater infiltration amount, calculating the urban area rainwater infiltration amount:
urban area rainwater infiltration amount=i z +I c
Wherein I is z I is the infiltration amount of rainwater on the stagnant water-permeable surface c Is the infiltration amount of the rainwater on the transmission type permeable surface.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310592499.9A CN116611249B (en) | 2023-05-23 | 2023-05-23 | Urban area rainwater infiltration amount calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310592499.9A CN116611249B (en) | 2023-05-23 | 2023-05-23 | Urban area rainwater infiltration amount calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116611249A true CN116611249A (en) | 2023-08-18 |
| CN116611249B CN116611249B (en) | 2023-11-21 |
Family
ID=87679611
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310592499.9A Active CN116611249B (en) | 2023-05-23 | 2023-05-23 | Urban area rainwater infiltration amount calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116611249B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107220496A (en) * | 2017-05-26 | 2017-09-29 | 上海市气象灾害防御技术中心 | A kind of urban rainstorm waterlogging assesses modeling method |
| CN111275266A (en) * | 2020-02-25 | 2020-06-12 | 华南师范大学 | Rapid optimization method for spatial layout of impervious ground surface for urban rainstorm and waterlogging prevention and control |
| WO2022032872A1 (en) * | 2020-08-14 | 2022-02-17 | 贵州东方世纪科技股份有限公司 | Big data-based hydrologic forecasting method |
| CN114547531A (en) * | 2022-02-25 | 2022-05-27 | 中国水利水电科学研究院 | Urban impervious surface effectiveness quantification method |
| CN114754838A (en) * | 2022-03-31 | 2022-07-15 | 北京城市气象研究院 | Early warning method suitable for natural slope and urban inland inundation water accumulation depth |
-
2023
- 2023-05-23 CN CN202310592499.9A patent/CN116611249B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107220496A (en) * | 2017-05-26 | 2017-09-29 | 上海市气象灾害防御技术中心 | A kind of urban rainstorm waterlogging assesses modeling method |
| CN111275266A (en) * | 2020-02-25 | 2020-06-12 | 华南师范大学 | Rapid optimization method for spatial layout of impervious ground surface for urban rainstorm and waterlogging prevention and control |
| WO2022032872A1 (en) * | 2020-08-14 | 2022-02-17 | 贵州东方世纪科技股份有限公司 | Big data-based hydrologic forecasting method |
| CN114547531A (en) * | 2022-02-25 | 2022-05-27 | 中国水利水电科学研究院 | Urban impervious surface effectiveness quantification method |
| CN114754838A (en) * | 2022-03-31 | 2022-07-15 | 北京城市气象研究院 | Early warning method suitable for natural slope and urban inland inundation water accumulation depth |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116611249B (en) | 2023-11-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111582755B (en) | Mountain torrent disaster comprehensive risk dynamic assessment method based on multi-dimensional set information | |
| CN111507375B (en) | Urban waterlogging risk rapid assessment method and system | |
| CN103940974B (en) | Based on the middling scale watershed soil erosion Dynamic Spatial And Temporal Analysis method of GIS | |
| CN109657841B (en) | Deep extraction method for urban rainstorm waterlogging | |
| CN108416468B (en) | Flood advanced early warning and forecasting method | |
| CN103093114B (en) | A kind of distributed basin water deficit measuring method based on landform and soil characteristic | |
| CN110232737B (en) | Urban catchment area division method | |
| CN114647881B (en) | Urban waterlogging modeling method considering microscopic hydrologic process of building | |
| CN115407049B (en) | Method for measuring and calculating carbon sink potential of soil in water and soil conservation engineering measures | |
| CN113435630B (en) | Basin hydrological forecasting method and system with self-adaptive runoff yield mode | |
| CN112268595A (en) | Method for monitoring regional average groundwater level | |
| CN116611249B (en) | Urban area rainwater infiltration amount calculation method | |
| CN113902259A (en) | Wet avalanche disaster risk assessment method and device | |
| CN112084643A (en) | Drainage basin extraction method based on digital elevation and soil parameters | |
| CN115859040A (en) | Method for calculating conversion coefficient of rainstorm point surface by adopting moving point fixed surface relation | |
| CN106980764B (en) | A kind of general hydrological simulation method directly evaporated using remote sensing | |
| CN110569595A (en) | data-free area rainfall station network site selection method based on numerical simulation | |
| CN111369102B (en) | Method and device for extracting waterlogging risk points | |
| Baldwin | Catchment-scale soil water retention characteristics and delineation of hydropedological functional units in the Shale Hills Catchment | |
| CN109409674B (en) | Drainage basin brown infrastructure current situation regulation and storage capacity assessment method | |
| CN116992783B (en) | Full-effective grid cell flow field simulation method under large-depth drainage of underground water | |
| CN111524054A (en) | Computing method for computing rechecked sponge index by using micro-partition/volume method | |
| Khanbilvardi et al. | Mathematical model of erosion and deposition on a watershed | |
| Manaouch et al. | Integrating WaTEM/SEDEM model and GIS-based FAHP Method for Identifying Ecological Rainwater Harvesting Sites in Ziz upper watershed, SE Morocco | |
| CN110674469B (en) | Hydrological frequency calculation method suitable for arid karst-deficient data area |
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 |