CN116049963B - Quick simulation method based on urban oversized flood underground space regulation - Google Patents
Quick simulation method based on urban oversized flood underground space regulation Download PDFInfo
- Publication number
- CN116049963B CN116049963B CN202310124926.0A CN202310124926A CN116049963B CN 116049963 B CN116049963 B CN 116049963B CN 202310124926 A CN202310124926 A CN 202310124926A CN 116049963 B CN116049963 B CN 116049963B
- Authority
- CN
- China
- Prior art keywords
- building
- underground space
- flood
- urban
- water
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 76
- 238000004088 simulation Methods 0.000 title claims abstract description 37
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 93
- 230000008569 process Effects 0.000 claims abstract description 40
- 238000004364 calculation method Methods 0.000 claims description 20
- 239000003673 groundwater Substances 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 4
- 239000006185 dispersion Substances 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 4
- 239000012530 fluid Substances 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 11
- 230000009467 reduction Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 11
- 238000004891 communication Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 6
- 238000006424 Flood reaction Methods 0.000 description 5
- 238000009825 accumulation Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 239000002352 surface water Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000011664 signaling Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
-
- 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)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a rapid simulation method based on urban oversized flood underground space regulation, which relates to the technical field of flood control and disaster reduction science and comprises the following steps: s1, collecting building distribution of a city, and setting a characteristic elevation value based on the building distribution of the city; s2, identifying the building according to the characteristic elevation value to obtain the size of an underground space of the building to be used as the water storage capacity of the building; s3, calculating to obtain the inflow of flood into the underground space of the building by utilizing the Saint View path group based on the water storage capacity of the building; s4, calculating to obtain an urban underground space regulation model through a shallow water equation based on water inflow, and using the model to complete quick simulation of urban oversized flood underground space regulation; the invention simulates the influence of the flood on the urban flood process after entering the underground space, thereby realizing the regulation and storage effect of the simulated underground space on the urban oversized flood.
Description
Technical Field
The invention relates to the technical field of flood control and disaster reduction science, in particular to a rapid simulation method based on urban oversized flood underground space regulation.
Background
Under extreme rainstorm conditions, urban inland inundation can often break through urban underground space flood control facilities, so that underground space is flushed, economic loss is caused, and the urban inland inundation evolution process is influenced. Under climate change conditions, the frequency of occurrence of extreme weather increases; and with the development of the urban process, the urban underground space is developed and utilized more and more fully, which aggravates the loss of underground space after water inflow and the influence of the underground space on the urban waterlogging evolution process.
In the past urban flood simulation, modeling work aiming at the underground space water inflow simulation process is complex, and more data are needed; the water quantity entering the underground space is calculated by adopting a weir flow formula, so that the description of the flood process is simpler; the method is only suitable for small-scale fine simulation, and is difficult for large-scale modeling and simulation of city scale. However, in urban oversized flooding scenarios, the impact on urban flood evolution after the underground space is filled with water is not negligible. Therefore, a rapid simulation method based on urban oversized flood underground space regulation is urgently needed to solve the problems.
Disclosure of Invention
The invention aims to provide a rapid simulation method for regulating and accumulating an underground space based on urban oversized floods, which can rapidly, stably and accurately calculate the amount of floods entering the underground space on the basis of not increasing modeling workload and required data, simulate the influence of floods entering the underground space on urban floods, and further realize the regulation and accumulating effect of the simulated underground space on the urban oversized floods.
The technical scheme of the invention is as follows:
in a first aspect, the present application provides a method for rapid simulation of urban oversized flood underground space regulation, comprising the steps of:
s1, collecting building distribution of a city, and setting a characteristic elevation value based on the building distribution of the city;
s2, identifying the building according to the characteristic elevation value to obtain the size of an underground space of the building to be used as the water storage capacity of the building;
s3, calculating to obtain the inflow of flood into the underground space of the building by utilizing the Saint View path group based on the water storage capacity of the building;
s4, calculating to obtain an urban underground space regulation model through a shallow water equation based on water inflow, and using the model to complete quick simulation of urban oversized flood underground space regulation.
Further, step S2 includes the steps of:
s201, identifying building units with adjacent spatial positions as the same building object;
s202, calculating the building area according to the characteristic elevation value and the building object and the unit area contained in the building object, calculating and identifying the size of the underground space of the building object according to the area and the building type of the building object, and giving the water storage capacity of the building object as the characteristic attribute of water storage.
Further, in step S202, the calculation formula of the water storage capacity is as follows:
wherein V is Total (S) For the water storage capacity of the current building, α is the building type coefficient, S is the area of the current building,is the average height of the local underground space.
Further, in step S3, the process of calculating the inflow of flood into the underground space of the building using the san francisco procedure set is:
wherein A is the water flow cross-section area, t is the time, Q is the flow, x is the distance, g is the gravity acceleration, H is the water head difference, S f The head loss per unit length, n is Manning coefficient, R is hydraulic radius, and U is flow velocity.
Further, the above process of calculating the inflow of flood into the underground space of the building by using the san france equation group further includes the application of a finite difference method, and the simplified calculation according to the characteristics of the ideal river model, wherein the calculation process is as follows:
wherein Deltat is the time step, deltaQ is the variation of the flow in Deltat time,for average flow rate +.>A is the average overcurrent area variation in delta t time 2 For the water passing area of the flow position, A 1 For the water passing area of the inflow position, L is idealLength of river course(s)>For average water passing area H 2 Is a ground water head H 1 Is the water head of the underground space, eta is a parameter, n is a Manning coefficient, < ->Is the average hydraulic radius.
Further, in step S4, the formula of the urban underground space regulation model obtained by calculation through the shallow water equation is as follows:
wherein h is water depth, t is time, u and v are flow velocity in x and y directions respectively, sce is fluid source or sink, Z is free surface elevation, v e For effective viscosity coefficient, g is gravitational acceleration, F x And F y Is the friction term in the x and y directions.
The calculation of the water inflow in the step S3 and the simulation of the urban underground space regulation model in the step S4 are alternately performed and are mutually influenced. Specifically, the water inflow in the step S3 influences the surface water accumulation of the urban underground space regulation model in the step S4; in contrast, the ground water depth calculated in step S4 also affects the water inflow of the underground space in step S3, and the two affect each other and alternate to complete the rapid simulation of the regulation of the underground space of the urban oversized flood.
In a second aspect, the present application provides an electronic device, including:
a memory for storing one or more programs;
a processor;
when the one or more programs are executed by the processor, a fast simulation method based on urban oversized flooding underground space regulation as in any of the above first aspects is implemented.
In a third aspect, the present application provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements a fast simulation method based on urban extra-large flood underground space regulation as in any one of the first aspects above.
Compared with the prior art, the invention has at least the following advantages or beneficial effects:
(1) According to the rapid simulation method based on urban oversized flood underground space regulation, the distribution of the building is approximately used as the distribution of the underground space according to the distribution characteristics of the underground space of the modern city, so that the collection work of the distribution data of the underground space is avoided;
(2) The invention is based on modeling work of the traditional urban flood evolution model, can realize rapid identification of building objects and endow the building objects with characteristic properties;
(3) In the flood evolution process, the water quantity of flood entering the underground space is calculated according to the Saint Violet range group, so that the problems of numerical oscillation and inaccurate exchange flow existing in the process of calculating the exchange flow by a weir flow formula are solved;
(4) According to the invention, under the situation of urban oversized flood, the underground space can simulate the regulation and storage effect of the flood process, so that the preparation work can be simply and rapidly completed, the evolution process of the urban oversized flood can be more accurately and thoroughly reflected, the efficient and accurate calculation of the urban oversized flood can be realized, and the improvement of the simulation precision of the urban oversized flood is facilitated.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings that are needed in the embodiments will be briefly described below, it being understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a step diagram of a rapid simulation method based on urban oversized flood underground space regulation;
FIG. 2 is a diagram of a dual grid link;
FIG. 3 is a graph of the water storage change of a rainfall underground space in city A;
FIG. 4 is a flooding pattern of an area of city A6 hours after rainfall occurs;
fig. 5 is a schematic block diagram of an electronic device.
Icon: 101. a memory; 102. a processor; 103. a communication interface.
Detailed Description
For the purposes of making the objects, technical solutions and advantages of the embodiments of the present application more clear, the technical solutions of the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is apparent that the described embodiments are some embodiments of the present application, but not all embodiments. The components of the embodiments of the present application, which are generally described and illustrated in the figures herein, may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present application, as provided in the accompanying drawings, is not intended to limit the scope of the application, as claimed, but is merely representative of selected embodiments of the application. All other embodiments, which can be made by one of ordinary skill in the art based on the embodiments herein without making any inventive effort, are intended to be within the scope of the present application.
It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further definition or explanation thereof is necessary in the following figures.
It should be noted that, in this document, the term "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrase "comprising … …" does not exclude the presence of other like elements in a process, method, article or apparatus that comprises the element.
In the description of the present application, it should also be noted that, unless explicitly specified and limited otherwise, the terms "disposed," "connected," and "connected" are to be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communication between two elements. The specific meaning of the terms in this application will be understood by those of ordinary skill in the art in a specific context.
Some embodiments of the present application are described in detail below with reference to the accompanying drawings. The various embodiments and features of the embodiments described below may be combined with one another without conflict.
Example 1
Referring to fig. 1, fig. 1 is a step diagram of a rapid simulation method for urban oversized flood underground space regulation according to an embodiment of the present application.
The application provides a rapid simulation method based on urban oversized flood underground space regulation, which comprises the following steps:
s1, collecting building distribution of a city, and setting a characteristic elevation value based on the building distribution of the city;
s2, identifying the building according to the characteristic elevation value to obtain the size of an underground space of the building to be used as the water storage capacity of the building;
s3, calculating to obtain the inflow of flood into the underground space of the building by utilizing the Saint View path group based on the water storage capacity of the building;
s4, calculating to obtain an urban underground space regulation model through a shallow water equation based on water inflow, and using the model to complete quick simulation of urban oversized flood underground space regulation.
The beneficial effects of the invention are as follows: according to the underground space distribution characteristics of modern cities, the distribution of buildings is approximately used as the distribution of underground spaces, so that the collection work of underground space distribution data is avoided; based on modeling work of a traditional urban flood evolution model, the building object is rapidly identified, and characteristic attributes are given to the building object; in the flood evolution process, the water quantity of flood entering the underground space is calculated according to the Saint View south equation group, and the problems of numerical oscillation and inaccurate exchange flow existing in the process of calculating the exchange flow by using a weir flow formula are solved. By using the method for simulating the regulation effect of the underground space on the flood process under the condition of the urban oversized flood, the preparation work can be simply and rapidly completed, the evolution process of the urban oversized flood can be more accurately and thoroughly reflected, the efficient and accurate calculation of the urban oversized flood can be realized, and the simulation precision of the urban oversized flood can be improved.
As a preferred embodiment, step S2 comprises the steps of:
s201, identifying building units with adjacent spatial positions as the same building object;
s202, calculating the building area according to the characteristic elevation value and the building object and the unit area contained in the building object, calculating and identifying the size of the underground space of the building object according to the area and the building type of the building object, and giving the water storage capacity of the building object as the characteristic attribute of water storage.
In particular, in step S202, for the unstructured triangular mesh data structure with unobvious regular spatial coordinate distribution, the proximity judgment is difficult, and a structured rectangular mesh with similar mesh size can be established to link with the nodes in the triangular mesh, the linking method is shown in fig. 2, and fig. 2 is a double-mesh linking schematic diagram. And accordingly, performing rapid approach judgment, wherein the judgment steps are as follows:
(1) cyclically identifying a first building node in the triangular mesh;
(2) linking to corresponding units in the rectangular grid, and searching adjacent rectangular units according to the coordinates of the rectangular grid;
(3) the adjacent rectangular units are linked back to the corresponding nodes in the triangular meshes; wherein, due to different mesh types, there may be more than one node linked to;
(4) calculating the distance between nodes to confirm whether the adjacent node is an adjacent node and judging whether the adjacent node is a building unit;
(5) if the building unit is a building unit, the steps (2) - (5) are circulated until all nodes contained in the building are all identified;
(6) after the end of identifying a single building, returning to step (1) again, starting the identification of the next building object until all the buildings are identified.
The beneficial effects of the above steps are: the double-grid link data structure is adopted to establish the link between the rectangular grid and the triangular grid nodes, so that the characteristics of simple rectangular grid structure, obvious grid coordinate rule and quick approach judgment can be brought into play, and the advantages of finer boundary characterization, strong rough difference resistance and high reliability of the triangular grid in flood calculation can be maintained. And the steps are designed based on a computer to store the grid data, so that automatic approach judgment is facilitated, and the building can be identified quickly.
As a preferred embodiment, the water storage capacity in S202 is calculated by the following formula:
wherein V is Total (S) The water storage capacity of the current building; alpha is a building type coefficient; s is the area of the current building;is the average height of the local underground space.
In particular, the heights of the underground spaces of different building types are often different, the building types can be divided into residential buildings, office buildings and business buildings, and different alpha values are given to the different building types, so that the description of different sizes of the underground spaces of different buildings is realized.
The beneficial effects of the above-mentioned further scheme are: a double-grid linked data structure may be employed for the triangular grid cells to enable fast approach determination. And each step is designed based on a computer related data storage format, is suitable for computer execution, and is convenient for automatic processing so as to quickly realize identification and attribute addition of building objects. The method can adapt to the existing data and modeling method, and is suitable for modeling and simulating underground space in a large scale of city scale.
As a preferred embodiment, the inflow of the flood into the underground space of the building in step S3 is calculated mainly by the san francisco group.
The ideal river model is a hydraulic model abstracted from the flow process of the underground space water inlet, and has the characteristics of uniform gradient, fixed section shape, uniform roughness and the like, so that the calculation of the inflow is facilitated. Compared with a common weir flow formula for calculating flow, the Saint View southern range group can describe the water flow exchange process more accurately and obtain more accurate and stable exchange flow.
Specifically, the san velam equation set includes a continuous equation (2) and a momentum equation (3), and the specific form is as follows:
wherein A is the water flow cross-section area, t is the time, Q is the flow, x is the distance, g is the gravity acceleration, H is the water head difference, S f Is the head loss per unit length.
Wherein S is f The calculation is performed according to the following formula:
wherein n is a Manning coefficient, R is a hydraulic radius, and U is a flow velocity.
The combination of the continuity equation (2) and the momentum equation (3) of the san velam equation set can obtain the following target equation:
the target equation (5) is solved by applying a finite difference method and bringing the difference approximation into the equation (5), so that the following equation can be obtained:
wherein Deltat is the time step, deltaQ is the variation of the flow in Deltat time,for average flow rate +.>A is the average overcurrent area variation in delta t time 2 For the water passing area of the flow position, A 1 L is the ideal river length for the water passing area of the inflow position>For average water passing area H 2 Is a ground water head H 1 Is the underground space water head, eta is a parameter, and +.>n is Manning coefficient, <>Is the average hydraulic radius.
Equation (6) can be simplified as follows according to the characteristics of the ideal river model:
to meet the large-scale modeling and computing requirements, existing data is adapted, decentralized corrections are made to the flow calculated by the san View equation, and the water inflow to calculate individual building objects is described as follows:
wherein,for the inflow of the ith building in Δt time, j is the number of the surrounding nodes of the building, (j=1, 2,., n), n is the total number of surrounding nodes of the building, α 1 For dispersion factor (underground space water inlet width/building boundary length),>for the ideal river channel exchange flow calculated by the formula (7) in the t time of the ith building, V Remainder of the process Is the underground residual water storage space.
It should be noted that, by setting the critical water level at which the surface water begins to enter the underground space, the flood control facility of the underground space is truly simulated. And an ideal river channel model is established for describing the flow process of the underground space water inlet, the flow exchange process can be more accurately described by solving the Saint Vignan equation to calculate the underground space water inflow, the more stable exchange flow is obtained, and the dispersion coefficient alpha is introduced, so that the acquisition and modeling work of the underground space inlet distribution data are avoided, and the large-scale simulation considering the underground space regulation and storage effect can be rapidly performed.
As a preferred embodiment, in step S4, the form of the shallow water equation is as follows:
wherein h is water depth, t is time, u and v are flow velocity in x and y directions respectively, sce is fluid source or sink, Z is free surface elevation, v e For effective viscosity coefficient, g is gravitational acceleration, F x And F y Is the friction term in the x and y directions.
The method for describing the flow rate of the underground space and urban flood exchange comprises the following steps: and (3) dynamically assigning the calculated water exchange quantity of the flood and the underground space to a variable Sce, and dynamically simulating the water flow process of the ground flood into the underground space in the urban flood evolution calculation process. And further simulate the regulation and storage effect of the underground space on urban flood evolution.
Example 2
Assuming that the rapid simulation of the underground space regulation effect in the ultra-large flood process of the city A is required, according to the investigation of the underground space of the city, the distribution of the underground space is found to be basically consistent with the building distribution, and the condition that the building distribution is adopted to replace the underground space distribution is satisfied. And the waterproof capacity of the underground space in A market is that the accumulated water around is less than or equal to 0.5m; the average height of the underground space is 2m; the gradient of the underground space inlet is 0.075; the ratio of the entrance width of the underground space to the perimeter of the building was 0.1. According to the conditions, the simulation of the underground space regulation effect in the process of applying extra-large flooding to the city A is carried out.
After the land utilization and elevation data for city a is obtained, conventional modeling work is performed on the area, and in assigning a grid to the elevation data for the building portion is assigned 9999 to identify the distribution of the building. In the calculation process, a unit with the height of 9999 is identified as a building unit and is continuously determinedIs identified as a building object and is given its characteristic properties. The average height of the underground space of the city is 2m, and the height coefficient of the underground space of different types of buildings is alpha 1 ，α 2 ，α 3 ，α 4 .., the water storage capacity calculation formula is:
V total (S) ＝2αS (12)
In the flood evolutionary computation process, the water depth around the building is obtained, when the water depth around the building is all smaller than 0.5m, underground space water inflow computation is not performed, and when the water depth larger than 0.5m exists, computation is started. Firstly, confirming the imaging parameters of an ideal river model, wherein the gradient is 0.075, the section shape is rectangular, the roughness is 0.017, the gravity acceleration is 9.8m/s2, and the calculation step length is 1s, so that the formula (7) can be embodied as follows:
since the ratio of the entrance width of the underground space to the perimeter of the building is 0.1, α=0.1, and Δt=1s, formula (8) can be embodied as:
according to the exchange water flow characteristic at the time t-1 and the water flow characteristic of the ground and the underground space at the time t, carrying out calculation in a formula (10), and carrying the calculation result into a formula (11), the water flow entering the underground space at the time t can be calculated.
In the flood evolution process, the exchanged water quantity of the ground flood and the underground space is assigned to a source-sink item in a shallow water equation so as to achieve the hydraulic process of simulating the ground flood to enter the underground space in the flood evolution process. In the continuous evolution process of the flood simulation, the water quantity entering the underground space is continuously and circularly calculated, and is assigned to a source and sink item, so that the flood process simulation of the extra-large city considering the regulation and storage effect of the underground space can be completed.
As shown in FIG. 3, a diagram of the water storage change of a rainfall underground space in a city A is shown in FIG. 3. As can be seen from the graph, the water inflow of the underground space is basically 0 in 1h before rainfall, the ground water accumulation of the two scenes has no obvious difference, the water storage of the underground space is rapidly increased in the period of 1-3h, the time period is the time for the ground water accumulation to break through the underground space flood control facility, the ground water accumulation difference of the two scenes is gradually increased, the maximum ground water accumulation is obviously smaller than that of the non-considered scenes in consideration of the water inflow of the underground space, the difference is 623444m3, and the water storage of the underground space is reduced by 16.8%. After 4 hours, the water storage capacity of the underground space is slowly increased in the water-withdrawal stage of urban inland inundation, and the difference between the two situations is gradually increased.
Fig. 4 is a view showing the flooding of a city a in 6 hours after the occurrence of rainfall, and it can be seen that the flooding depth of the region is significantly reduced by considering the underground space regulation effect. Therefore, the flood simulation considering the underground space can reduce the severity of urban waterlogging, and has obvious regulation and storage effects on urban oversized flooding.
Example 3
Referring to fig. 5, fig. 5 is a schematic block diagram of an electronic device according to embodiment 3 of the present application.
An electronic device comprises a memory 101, a processor 102 and a communication interface 103, wherein the memory 101, the processor 102 and the communication interface 103 are directly or indirectly electrically connected with each other to realize data transmission or interaction. For example, the components may be electrically connected to each other via one or more communication buses or signal lines. The memory 101 may be used to store software programs and modules that are stored within the memory 101 for execution by the processor 102 to perform various functional applications and data processing. The communication interface 103 may be used for communication of signaling or data with other node devices.
The Memory 101 may be, but is not limited to, a random access Memory (Random Access Memory, RAM), a Read Only Memory (ROM), a programmable Read Only Memory (Programmable Read-Only Memory, PROM), an erasable Read Only Memory (Erasable Programmable Read-Only Memory, EPROM), an electrically erasable Read Only Memory (Electric Erasable Programmable Read-Only Memory, EEPROM), etc.
The processor 102 may be an integrated circuit chip with signal processing capabilities. The processor 102 may be a general purpose processor including a central processing unit (Central Processing Unit, CPU), a network processor (Network Processor, NP), etc.; but also digital signal processors (Digital Signal Processing, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), field programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.
It will be appreciated that the configuration shown in the figures is illustrative only and that a rapid simulation method based on urban oversized flood subterranean space regulation may also include more or fewer components than shown in the figures or have a different configuration than shown in the figures. The components shown in the figures may be implemented in hardware, software, or a combination thereof.
In the embodiments provided in the present application, it should be understood that the disclosed method may be implemented in other manners as well. The above-described embodiments are merely illustrative, for example, of the flowcharts or block diagrams in the figures that illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present application. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
In addition, the functional modules in the embodiments of the present application may be integrated together to form a single part, or each module may exist alone, or two or more modules may be integrated to form a single part.
The functions, if implemented in the form of software functional modules and sold or used as a stand-alone product, may be stored in a computer-readable storage medium. Based on such understanding, the technical solution of the present application may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, including several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
The foregoing description is only of the preferred embodiments of the present application and is not intended to limit the same, but rather, various modifications and variations may be made by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.
It will be evident to those skilled in the art that the present application is not limited to the details of the foregoing illustrative embodiments, and that the present application may be embodied in other specific forms without departing from the spirit or essential characteristics thereof. The present embodiments are, therefore, to be considered in all respects as illustrative and not restrictive, the scope of the application being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.
Claims (5)
1. A rapid simulation method based on urban oversized flood underground space regulation is characterized by comprising the following steps:
s1, collecting building distribution of a city, and setting a characteristic elevation value based on the building distribution of the city;
s2, identifying the building according to the characteristic elevation value to obtain the size of an underground space of the building to be used as the water storage capacity of the building;
step S2 comprises the following sub-steps:
s201, identifying building units with adjacent spatial positions as the same building object;
s202, calculating the building area according to the characteristic elevation value and the building object and the unit area contained in the building object, calculating and identifying the size of the underground space of the building object according to the area and the building type of the building object, and giving the water storage capacity of the building object with the characteristic attribute of water storage;
in step S202, the calculation formula of the water storage capacity is:
wherein V is Total (S) For the water storage capacity of the current building, α is the building type coefficient, S is the area of the current building,is the average height of the local underground space;
s3, calculating to obtain the inflow of flood into the underground space of the building by utilizing the Saint View path group based on the water storage capacity of the building;
s4, calculating to obtain an urban underground space regulation model through a shallow water equation based on water inflow, and using the model to complete quick simulation of urban oversized flood underground space regulation;
in step S3, the process of calculating the inflow of flood into the underground space of the building using the san francisco procedure set is:
wherein A is the water flow cross-section area, t is the time, Q is the flow, x is the distance, g is the gravity acceleration, H is the water head difference, S f The head loss is per unit length, n is Manning coefficient, R is hydraulic radius, and U is flow velocity;
dispersion corrections are made to the flow calculated by the san france program and the water inflow to calculate individual building objects is described as follows:
wherein,for the inflow of the ith building in delta t time, j is the number of the surrounding nodes of the building, j=1, 2, …, n, n is the total number of the surrounding nodes of the building, alpha 1 For dispersion factor, i.e. underground space water inlet width/building boundary length +.>Is the firstideal river channel exchange flow of i buildings in t time, V Remainder of the process Is the underground residual water storage space.
2. The rapid simulation method based on urban oversized flood underground space regulation of claim 1, wherein the process of calculating the inflow of flood into the underground space of the building by using the san france equation group further comprises the steps of applying a finite difference method and performing simplified calculation according to the characteristics of an ideal river model, wherein the calculation process is as follows:
wherein Deltat is the time step, deltaQ is the variation of the flow in Deltat time,for average flow rate +.>A is the average overcurrent area variation in delta t time 2 For the water passing area of the flow position, A 1 L is the ideal river length for the water passing area of the inflow position,for average water passing area H 2 Is a ground water head H 1 Is the water head of the underground space, eta is a parameter, n is a Manning coefficient, < ->Is the average hydraulic radius.
3. The rapid simulation method based on urban oversized flood underground space regulation as claimed in claim 1, wherein in the step S4, the formula for obtaining the urban underground space regulation model through calculation of the shallow water equation is as follows:
wherein h is water depth, t is time, u and v are flow velocity in x and y directions respectively, sce is fluid source or sink, Z is free surface elevation, v e For effective viscosity coefficient, g is gravitational acceleration, F x And F y Is the friction term in the x and y directions.
4. An electronic device, comprising:
a memory for storing one or more programs;
a processor;
a fast simulation method based on urban extra-large flood underground space regulation according to any one of claims 1-3, when said one or more programs are executed by said processor.
5. A computer readable storage medium, on which a computer program is stored, which computer program, when being executed by a processor, implements a fast simulation method based on urban extra-large flood underground space regulation as claimed in any one of claims 1-3.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211591526 | 2022-12-12 | ||
| CN2022115915262 | 2022-12-12 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116049963A CN116049963A (en) | 2023-05-02 |
| CN116049963B true CN116049963B (en) | 2024-01-23 |
Family
ID=86131377
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310124926.0A Active CN116049963B (en) | 2022-12-12 | 2023-02-07 | Quick simulation method based on urban oversized flood underground space regulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116049963B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101446275B1 (en) * | 2014-02-17 | 2014-10-02 | (주) 한길산업 | Rainwater retention and infiltration system preventing road and sidewalk flood damage and lowering non-point source cotaminant and its construction method |
| CN112785053A (en) * | 2021-01-15 | 2021-05-11 | 北京市水科学技术研究院 | Method and system for forecasting urban drainage basin flood |
| CN113191054A (en) * | 2021-04-30 | 2021-07-30 | 西安理工大学 | High-precision urban rainfall flood simulation method based on display card accelerated coupling pipe network |
| CN114647881A (en) * | 2022-03-14 | 2022-06-21 | 南京师范大学 | Urban inland inundation modeling method considering building micro hydrological process |
| CN114970315A (en) * | 2022-04-19 | 2022-08-30 | 河海大学 | Urban accumulated water simulation and rapid prediction method based on spatial dynamic characteristic deep learning |
| CN115391712A (en) * | 2022-05-18 | 2022-11-25 | 河海大学 | Urban flood risk prediction method |
-
2023
- 2023-02-07 CN CN202310124926.0A patent/CN116049963B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101446275B1 (en) * | 2014-02-17 | 2014-10-02 | (주) 한길산업 | Rainwater retention and infiltration system preventing road and sidewalk flood damage and lowering non-point source cotaminant and its construction method |
| CN112785053A (en) * | 2021-01-15 | 2021-05-11 | 北京市水科学技术研究院 | Method and system for forecasting urban drainage basin flood |
| CN113191054A (en) * | 2021-04-30 | 2021-07-30 | 西安理工大学 | High-precision urban rainfall flood simulation method based on display card accelerated coupling pipe network |
| CN114647881A (en) * | 2022-03-14 | 2022-06-21 | 南京师范大学 | Urban inland inundation modeling method considering building micro hydrological process |
| CN114970315A (en) * | 2022-04-19 | 2022-08-30 | 河海大学 | Urban accumulated water simulation and rapid prediction method based on spatial dynamic characteristic deep learning |
| CN115391712A (en) * | 2022-05-18 | 2022-11-25 | 河海大学 | Urban flood risk prediction method |
Non-Patent Citations (4)
| Title |
|---|
| Chang-He Li等.A Flooding Warning System based on RFID Tag Array for Energy Facility.IEEE.2018,全文. * |
| 城市下垫面空间格局对社区尺度内涝过程的影响模拟;刘家宏等;水科学进展;参见第1节 * |
| 城市水文水动力藕合模型及其应用研究;梅超;中国优秀博硕士学位论文全文数据库;参见第二、三章 * |
| 基于区域排水模型的地下空间内涝危险性评估;王渲等;中国农村水利水电(第1期期);参见第1-2节 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116049963A (en) | 2023-05-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Westoby et al. | Numerical modelling of glacial lake outburst floods using physically based dam-breach models | |
| Maskell et al. | Investigating river–surge interaction in idealised estuaries | |
| Zhu et al. | An effective depression filling algorithm for DEM-based 2-D surface flow modelling | |
| Nardi et al. | Quantification of riverbank erosion and application in risk analysis | |
| Lafdani et al. | Research article daily rainfall-runoff prediction and simulation using ANN, ANFIS and conceptual hydrological MIKE11/NAM models | |
| CN115130396A (en) | Distributed hydrological model modeling method for riverway type reservoir area | |
| CN107978138A (en) | A kind of disaster monitoring method for early warning based on mountain torrents dynamical evolution simulation model | |
| Morianou et al. | Assessing hydro-morphological changes in Mediterranean stream using curvilinear grid modeling approach-climate change impacts | |
| Iskender et al. | Evaluation of surface runoff estimation in ungauged watersheds using SWAT and GIUH | |
| CN116049963B (en) | Quick simulation method based on urban oversized flood underground space regulation | |
| Avila-Aceves et al. | Geospatial modelling of floods: A literature review | |
| Bashar | Comparative performance of soil moisture accounting approach in continuous hydrologic simulation of the Blue Nile | |
| Frank et al. | Use of an integrated one dimensional-two dimensional hydraulic modelling approach for flood hazard and risk mapping | |
| Kurniyaningrum et al. | Sensitivity of flow depth inundation based on the micro-scale topography in Krukut River, Jakarta, Indonesia | |
| CN116561476A (en) | Method for automatically realizing small-basin mountain torrent flood flow and inundation forecast | |
| Khan et al. | Assessment of hydro-morphological change of Surma-Kushiyara River system | |
| Gong et al. | A case study of flood modeling with adaptive mesh refinement | |
| Alemseged et al. | Uncertainty issues in hydrodynamic flood modeling | |
| Hraniciuc et al. | Creating Flood Hazard Maps Using 2d Hydraulic Models | |
| Nwe et al. | Sediment transport due to flooding from levee breach: a numerical application of the 2004 Niigata flood | |
| Zulkifli et al. | River water level prediction for flood risk assessment using NARX neural network | |
| Bahtiar et al. | A Numerical Algorithm and Visualization Software for Flood Simulation in Urban Area: A Case Study of West Jakarta, Indonesia | |
| Hadihardaja et al. | Decision support system for predicting flood characteristics based on database modelling development (case study: Upper Citarum, West Java, Indonesia) | |
| Mahmood et al. | Assessment of drainage capacity of Chaktai and Rajakhali Khal in Chittagong City and inundation adjacent of urban areas | |
| CN113094876B (en) | Reservoir flood control risk calculation method and system based on ARIMA model under future situation |
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 |