Method and device for simulating action process of surface ponding coupled with urban pipe network

Abstract

The invention provides a surface ponding action process simulation method and device coupled with an urban pipe network. The method comprises the steps that underlying surface data, pipe network data, rainfall data and actually-measured hydrological data of a research area are acquired; constructing a one-dimensional drainage pipe network model of the research area according to the underlying surface data, the pipe network data and the rainfall data; according to the underlying surface data and the rainfall data, constructing an earth surface two-dimensional submerging model of the research area; determining at least one simulation session and a simulation total duration corresponding to each simulation session according to the rainfall data and the actually measured hydrological data; and based on a preset time step, performing data interaction between the one-dimensional drainage pipe network model and the earth surface two-dimensional submerging model according to the overflow data and the backflow data of the plurality of pipe network nodes until the total simulation time length of each simulation session is reached, so as to complete a simulation event for the surface ponding overflow and regression process of the research area. According to the method, the comprehensiveness and accuracy of urban rainstorm waterlogging simulation can be improved.

Description

Technical Field

[0001] The present application relates to the field of urban hydrology and pipe network drainage, and in particular to a method and device for simulating the surface water action process coupled to an urban pipe network. Background Art

[0002] Urban rainstorm waterlogging simulation is an important means to formulate urban flood prevention and mitigation measures and forecast and warn of rainstorm waterlogging. Timely and accurate simulation of urban rainstorm waterlogging process is of great significance to improving urban waterlogging risk prevention capabilities and reducing waterlogging disaster losses. In urban waterlogging simulation, the Storm Water Management Model (SWMM) is the most widely used, but this model cannot accurately give the surface water depth and flooding range; in addition, the hydrodynamic model based on the two-dimensional shallow water equation can characterize the surface runoff movement process under complex terrain conditions, but this model does not consider the influence of underground drainage pipe networks.

[0003] At present, the existing simulation method for the surface water action process of coupled urban pipe networks usually connects the SWMM pipe network model and the two-dimensional hydrodynamic model or flooding model in a unidirectional manner, that is, the overflow of the pipe network node simulated by the SWMM model is used as the point source boundary condition to drive the two-dimensional hydrodynamic model or flooding model. Although the implementation process is relatively simple and easy to use, it ignores the backflow of surface water through the pipe network nodes. The simulation results of urban waterlogging lack accuracy and comprehensiveness, and it is difficult to meet the simulation requirements for the overflow and receding process of surface water in coupled urban pipe networks. Summary of the invention

[0004] The present application provides a method and device for simulating the process of ground waterlogging coupled to an urban pipe network, which couples a one-dimensional drainage pipe network model and a two-dimensional surface inundation model in a bidirectional manner, simulates in steps, and drives each other, thereby realizing a two-way water exchange between the drainage pipe network and its overlying surface, accurately simulating the overflow and retreat processes of surface water, thereby improving the accuracy and comprehensiveness of urban waterlogging simulation.

[0005] In a first aspect, the present application provides a method for simulating a ground surface water action process coupled to an urban pipe network, the method comprising:

[0006] Obtain the underlying surface data, pipe network data, rainfall data and measured hydrological data of the study area;

[0007] A one-dimensional drainage pipe network model of the study area is constructed according to the underlying surface data, the pipe network data and the rainfall data; and a two-dimensional surface inundation model of the study area is constructed according to the underlying surface data and the rainfall data, wherein the one-dimensional drainage pipe network model is bidirectionally coupled with the two-dimensional surface inundation model through a plurality of pipe network nodes, and the plurality of pipe network nodes are exchange nodes for surface water flow and underground pipe flow;

[0008] Determine at least one simulation session and a total simulation duration corresponding to each simulation session according to the rainfall data and the measured hydrological data;

[0009] Based on a preset time step, data interaction is performed between the one-dimensional drainage network model and the two-dimensional surface inundation model according to the overflow data and backflow data of the multiple network nodes until the total simulation time of each simulation session is reached, so as to complete the simulation events of at least one simulation session of the overflow and receding process of ground water in the study area.

[0010] In a second aspect, the present application provides a ground surface water action process simulation device coupled to an urban pipe network, comprising:

[0011] The acquisition unit is used to obtain the underlying surface data, pipe network data, rainfall data and measured hydrological data of the study area;

[0012] A processing unit is used to construct a one-dimensional drainage pipe network model of the study area according to the underlying surface data, the pipe network data and the rainfall data; and, to construct a two-dimensional surface inundation model of the study area according to the underlying surface data and the rainfall data, wherein the one-dimensional drainage pipe network model is bidirectionally coupled with the two-dimensional surface inundation model through a plurality of pipe network nodes, and the plurality of pipe network nodes are exchange nodes for surface water flow and underground pipe flow; and, to determine at least one simulation session and a total simulation duration corresponding to each simulation session according to the rainfall data and the measured hydrological data; and, based on a uniform time step, to perform data interaction between the one-dimensional drainage pipe network model and the two-dimensional surface inundation model according to the overflow data and backflow data of the plurality of pipe network nodes until the total simulation duration of each simulation session is reached, so as to complete the simulation event of the at least one simulation session for the overflow and receding process of surface water in the study area.

[0013] In a third aspect, the present application provides an electronic device, including a processor and a memory, wherein the memory stores a computer program, and when the processor calls the computer program in the memory, it executes the method as described in any one of the first aspects.

[0014] It can be seen that in the embodiment of the present application, the processor obtains the underlying surface data, pipe network data, rainfall data and measured hydrological data of the study area; constructs a one-dimensional drainage pipe network model of the study area according to the underlying surface data, pipe network data and rainfall data; and constructs a two-dimensional surface inundation model of the study area according to the underlying surface data and rainfall data, the one-dimensional drainage pipe network model and the two-dimensional surface inundation model are bidirectionally coupled through multiple pipe network nodes, and the multiple pipe network nodes are exchange nodes for surface water flow and underground pipe flow; determines at least one simulation session and the total simulation time corresponding to each simulation session according to the rainfall data and the measured hydrological data; based on a preset time step, data interaction is performed between the one-dimensional drainage pipe network model and the two-dimensional surface inundation model according to the overflow data and backflow data of multiple pipe network nodes until the total simulation time of each simulation session is reached, so as to complete the simulation event of the at least one simulation session of the overflow and receding process of the surface water in the study area. In this way, compared with the existing urban waterlogging simulation scheme that is unidirectionally connected through a one-dimensional drainage network model and a two-dimensional surface inundation model, in the present application, the one-dimensional drainage network model and the two-dimensional surface inundation model are bidirectionally coupled and simulated separately. After completing the simulation of each time step, a data exchange is performed. The surface water flow and the underground pipe flow can be exchanged with each other through the network nodes, realizing the accurate simulation of the surface water overflow and the water retreat process, which is conducive to improving the accuracy and comprehensiveness of urban waterlogging simulation. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0016] Figure 1 is a structural schematic diagram of an electronic device provided in an embodiment of the present application;

[0017] Figure 2 It is a flow chart of a method for simulating the surface water action process of a coupled urban pipe network provided in an embodiment of the present application;

[0018] Figure 3 It is a schematic diagram of a simulation process of ground surface water action coupled to an urban pipe network provided in an embodiment of the present application;

[0019] Figure 4 It is a coupling node diagram of a one-dimensional drainage network model and a two-dimensional surface inundation model provided in an embodiment of the present application;

[0020] Figure 5It is a functional unit block diagram of a ground water surface action process simulation device coupled to an urban pipe network provided in an embodiment of the present application. DETAILED DESCRIPTION

[0021] In order to enable those skilled in the art to better understand the solution of the present application, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0022] The terms "first", "second", etc. in the specification and claims of this application and the above-mentioned drawings are used to distinguish different objects, rather than to describe a specific order. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions. For example, a process, method, system, product or device that includes a series of steps or units is not limited to the listed steps or units, but optionally includes steps or units that are not listed, or optionally includes other steps or units inherent to these processes, methods, products or devices.

[0023] Reference to "embodiments" herein means that a particular feature, structure, or characteristic described in conjunction with the embodiments may be included in at least one embodiment of the present application. The appearance of the phrase in various locations in the specification does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment that is mutually exclusive with other embodiments. It is explicitly and implicitly understood by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0024] In the embodiments of the present application, "and / or" describes the association relationship of the associated objects, indicating that three relationships may exist. For example, A and / or B can represent the following three situations: A exists alone; A and B exist at the same time; B exists alone. Among them, A and B can be singular or plural.

[0025] In the embodiment of the present application, the symbol " / " can indicate that the objects associated with each other are in an "or" relationship. In addition, the symbol " / " can also indicate a division sign, that is, performing a division operation. For example, A / B can indicate A divided by B.

[0026] In the embodiments of the present application, "at least one item" or similar expressions refer to any combination of these items, including any combination of single items or plural items, and refer to one or more, and multiple refers to two or more. For example, at least one item of a, b, or c can represent the following seven situations: a, b, c, a and b, a and c, b and c, a, b, and c. Among them, each of a, b, and c can be an element or a set containing one or more elements.

[0027] In the embodiments of the present application, "equal to" can be used in conjunction with greater than, and is applicable to the technical solution adopted when greater than, and can also be used in conjunction with less than, and is applicable to the technical solution adopted when less than. When equal to is used in conjunction with greater than, it is not used in conjunction with less than; when equal to is used in conjunction with less than, it is not used in conjunction with greater than.

[0028] At present, in the existing simulation scheme for the ground waterlogging process coupled with the urban pipe network, the one-dimensional drainage pipe network model and the two-dimensional surface inundation model are unidirectionally connected, that is, the overflow of the pipe network node simulated by the one-dimensional drainage pipe network model is used as the point source boundary condition to drive the two-dimensional surface inundation model. Although the implementation process is relatively simple and easy to use, the backflow of surface water through the pipe network nodes is ignored, and the simulation results of urban waterlogging are not accurate enough.

[0029] In response to the above problems, an embodiment of the present application provides a method and device for simulating the surface water action process coupled to an urban pipe network. The embodiment of the present application is described in detail below in conjunction with the accompanying drawings.

[0030] See also Figure 1 , Figure 1 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present application. For the sake of convenience, only the parts related to the embodiment of the present application are shown. For specific technical details not disclosed, please refer to the method part of the embodiment of the present application. The electronic device 1 can be any terminal device including a mobile phone, a computer, a PDA, a POS, a car computer, etc. Figure 1 As shown, the electronic device 1 includes a memory 20, a processor 10, a communication bus 40, a communication interface 30 and one or more programs 21, the one or more programs 21 are stored in the memory 20 and are configured to be executed by the processor 10, and the one or more programs 21 include instructions for executing any step in the following method embodiments. In a specific implementation, the processor 10 is used to execute any step in the following method embodiments, and when performing data transmission such as sending, the communication interface 30 can be selectively called to complete the corresponding operation.

[0031] In the embodiment of the present application, the processor 10 included in the computer device may have the function corresponding to any method step in the embodiment.

[0032] Those skilled in the art will understand that Figure 1 The computer device structure shown in the figure does not constitute a limitation on the computer device, and may include more or less components than shown in the figure, or combine certain components, or arrange the components differently.

[0033] See also Figure 2 , Figure 2 A schematic diagram of a flow chart of a method for simulating the surface water action process of a coupled urban pipe network provided in an embodiment of the present application, wherein the method is applied to Figure 1 The processor 10 shown, the method includes:

[0034] Step S201, the processor obtains the underlying surface data, pipe network data, rainfall data and measured hydrological data of the study area.

[0035] In a possible embodiment, the underlying surface data includes terrain data, elevation data and land use type, and the land use type includes construction land, green land, transportation land and public land;

[0036] The pipe network data includes rainwater pipe information, drain outlet information and node information. The rainwater pipe information includes the spatial distribution of the pipe network, pipe diameter, pipe section length, pipe roughness, and upstream and downstream bottom elevations. The drain outlet information includes the coordinates, elevation, and drainage type information of the drain outlet. The node information includes the node bottom elevation, well depth, and coordinate information. The rainwater pipe information, the drain outlet information, and the node information need to be restricted to the same coordinate system.

[0037] The rainfall data includes at least one rainfall amount of a measured event, and there is a corresponding relationship between the measured event and the simulated event;

[0038] The measured hydrological data include measured runoff data of pipelines and node monitoring water levels during measured rainfall events. The measured hydrological data correspond to the measured events of rainfall data, and the test duration of the measured events is the same as the total simulation duration of the corresponding simulation events.

[0039] Among them, the terrain data and elevation data in the underlying surface data can be obtained through the Geographic Information System (GIS), and the land use type can be obtained through remote sensing satellite image analysis, or obtained from urban planning departments, land administration bureaus, etc. In addition, the underlying surface data, pipe network data, rainfall data, and measured hydrological data of the study area in this embodiment can also be obtained through local water affairs departments or municipal engineering departments.

[0040] It can be understood that the rainfall data and measured hydrological data of the study area respectively include test data of at least one measured session. Among them, the rainfall data of each measured session corresponds one-to-one to the measured hydrological data of each measured session, and the start time and test duration of the same measured session are the same. Specifically, the rainfall data and measured hydrological data of the study area respectively include rainfall data and measured hydrological data of 60 sessions from January to October 2024. Then the start time of the first measured session can be 10:00 on January 2, and the test duration is 3 hours. Then the start time of the rainfall data and measured hydrological data of the first session are both 10:00 on January 2, and the test duration is 3 hours.

[0041] In step S202, the processor constructs a one-dimensional drainage network model of the study area according to the underlying surface data, the network data and the rainfall data; and constructs a two-dimensional surface inundation model of the study area according to the underlying surface data and the rainfall data, wherein the one-dimensional drainage network model is bidirectionally coupled with the two-dimensional surface inundation model through a plurality of network nodes, and the plurality of network nodes are exchange nodes for surface water flow and underground pipe flow.

[0042] Among them, the pipe network nodes mainly include rainwater inlets and inspection wells; the two-dimensional surface inundation model is constructed based on the two-dimensional shallow water equations (SWE), which is the control equation of the two-dimensional surface inundation model; and the one-dimensional drainage pipe network model uses the one-dimensional Saint-Venant equations to calculate the pipe network water flow, which is the control equation of the one-dimensional drainage pipe network model.

[0043] Specifically, the two-dimensional shallow water equations include the following formulas:

[0044]

[0045] Where h is the water depth, u and v are the flow velocities in the x and y directions respectively, t is the time, B(x,y) is the bottom slope elevation, τ is the friction term, the subscripts bx and by are the friction components of the friction term in the x and y directions respectively, and g is the gravitational acceleration.

[0046] Specifically, the one-dimensional Saint-Venant equations include the following formulas:

[0047]

[0048] Among them, Q is the flow rate, A is the cross-sectional area of ​​the water flow, t is the time, h is the water depth, g is the gravitational acceleration, S f is the friction slope; x represents the direction.

[0049] Step S203: The processor determines at least one simulation session and a total simulation duration corresponding to each simulation session according to the rainfall data and the measured hydrological data.

[0050] Among them, there is a corresponding relationship between the measured hydrological data and the measured times of the rainfall data, there is a corresponding relationship between the measured times and the simulated times, and the test duration of the measured times is the same as the total simulation duration of the corresponding simulated times. Therefore, the number of simulation times can be determined according to the number of measured times included in the rainfall data, and the total simulation duration of each simulation time can be determined according to the test duration of each measured time in the measured hydrological data.

[0051] For example, when the rainfall data includes 60 measured sessions from January to October 2024, the simulation sessions are also 60 times, and when the test duration of the first measured session is 3 hours, the total simulation duration of the first simulation session is also 3 hours, and when the test duration of the second measured session is 1 hour, the total simulation duration of the second simulation session is also 2 hours.

[0052] In step S204, the processor performs data interaction between the one-dimensional drainage network model and the two-dimensional surface inundation model based on a preset time step and according to the overflow data and backflow data of the multiple network nodes, until the total simulation time of each simulation session is reached, so as to complete the simulation events of at least one simulation session of the overflow and receding process of the ground water in the study area.

[0053] In a possible embodiment, performing data interaction between the one-dimensional drainage pipe network model and the two-dimensional surface inundation model according to the overflow data and the backflow data of the plurality of pipe network nodes includes:

[0054] In a first time step, performing surface water flow simulation by using the surface two-dimensional flooding model to determine backflow data of the plurality of pipe network nodes in the first time step, wherein the backflow data includes backflow volumes of the plurality of pipe network nodes;

[0055] Sending the backflow data to the one-dimensional drainage pipe network model to drive the one-dimensional drainage pipe network model to perform underground pipe flow simulation in a second time step;

[0056] In the second time step, the underground pipe flow simulation is performed by the one-dimensional drainage pipe network model to determine overflow data of the plurality of pipe network nodes in the second time step, wherein the overflow data includes overflow amounts of the plurality of pipe network nodes;

[0057] The input file of the surface two-dimensional flooding model is updated according to the overflow data, and the surface two-dimensional flooding model is driven according to the input file to repeat the surface water flow simulation and determine the backflow data of the multiple pipe network nodes in the third time step within a third time step.

[0058] In a possible embodiment, performing surface water flow simulation according to the surface two-dimensional inundation model to determine the plurality of pipe network nodes before the return flow data in the first time step; the method further comprises:

[0059] In a fourth time step, an underground pipe flow simulation is performed through the one-dimensional drainage pipe network model to determine the coordinates, water depths and flow rates of the pipe network nodes where overflow or backflow occurs in the fourth time step;

[0060] The surface two-dimensional inundation model is driven to perform the surface water flow simulation within the first time step according to the input file of the surface two-dimensional inundation model and the flow of the pipe network node, and the flow of the pipe network node is used as the input boundary condition for the surface two-dimensional inundation model to perform the surface water flow simulation.

[0061] It can be seen that in the embodiment of the present application, the one-dimensional drainage network model first simulates the one-dimensional drainage network water flow in the study area based on the underlying surface data, the network data and the rainfall data. After completing the simulation of one time step, the coordinates, water depth and flow rate of the network nodes where overflow or backflow occurs are recorded; then, the flow rate of the network nodes is used as the input boundary condition of the surface two-dimensional inundation model simulation, and the two-dimensional surface overflow and retreat process of the study area is simulated in combination with the model input file. After completing the simulation of one time step, the backflow data of multiple network nodes are determined; then, the backflow data is fed back to the one-dimensional drainage network model for underground network simulation. After completing the simulation of one time step, the overflow data of the network nodes is determined, and the surface two-dimensional inundation model input file is updated. At this point, a two-way data interaction between the one-dimensional drainage network model and the surface two-dimensional inundation model is completed. The one-dimensional drainage network model and the surface two-dimensional inundation model perform data interaction once each time a time step of simulation is completed until the total simulation time of each simulation session is reached, and the data interaction and simulation process are stopped.

[0062] In a possible embodiment, when the water level of the pipe network node is greater than the surface water level, overflow calculation is performed to obtain the overflow volume of the plurality of pipe network nodes, and, when the water level of the pipe network node is less than the surface water level, backflow calculation is performed to obtain the backflow volume of the plurality of pipe network nodes, and, when the water level of the pipe network node is equal to the surface water level, flow exchange is not performed, wherein,

[0063] The calculation formula for overflow and return flow of each pipe network node is as follows:

[0064]

[0065]

[0066] Among them, Q0 is the overflow of the pipe network node, Qs is the return flow of the network node, C0 is the orifice flow coefficient, and its value range is [0,1]; A m is the water flow area of ​​the network node, g is the gravitational acceleration, h 1d is the water level at the network node, h 2d is the surface water level, h z is the ground elevation, c w is the weir flow coefficient, ranging from [0,1]; w is the circumference of the network node or the width of the stormwater inlet.

[0067] Among them, restrictive conditions are usually added to the overflow and backflow of the pipeline network nodes to improve the stability of the data.

[0068] The calculation formula for adding restrictive conditions to the overflow volume of the pipe network node is as follows:

[0069]

[0070] Among them, V0 is the water storage capacity of the network node, Δt e is the time step of water flow exchange.

[0071] The calculation formula for adding restrictive conditions to the return flow of the pipe network node is as follows:

[0072]

[0073] Among them, V c is the current water volume in the surface grid cell corresponding to the pipe network node, Δt e is the time step of water flow exchange.

[0074] In a possible embodiment, the study area in the two-dimensional surface inundation model includes a plurality of closely arranged grid cells, the backflow data also includes the water depth of each grid cell, and the surface water flow simulation is performed by the two-dimensional surface inundation model to determine the backflow data of the plurality of pipe network nodes within the first time step, and the method further includes:

[0075] Get the initial water depth of each grid cell;

[0076] Determining the change in water accumulation in each grid cell after the first time step;

[0077] The calculation formula for the water flow change of each grid unit is as follows:

[0078] Q grid =Q up +Q down +Q left +Q right (3)

[0079]

[0080] Among them, Q grid is the change in water accumulation in the grid unit, Q up , Q down , Q left , Q right are the flow rates from the upstream, downstream, left and right grid cells respectively, which are obtained by formula (4); Q ij is the flow at the junction of adjacent grid cells i and j, A ij is the cross-sectional area at the junction of adjacent grid cells i and j, R ij is the hydraulic radius at the junction of adjacent grid cells i and j, S ij is the water surface slope between adjacent grid cells i and j, and n is the Manning coefficient;

[0081] The water accumulation depth of each grid unit is determined according to the water accumulation change, the initial water accumulation depth and the area of ​​each grid unit.

[0082] Among them, the surface two-dimensional inundation model mainly refers to the Storm Water Management Model, namely the SWMM model. The SWMM model can divide the study area into multiple subset waters, each subset water can be further divided into multiple grid cells, these grid cells can have different sizes, their size depends on the accuracy requirements of the model and the detail of the available data, the user can define the size of the grid cell according to the actual situation and the needs of the model. In addition, the SWMM model can also be used in conjunction with other geographic information systems, such as GIS and terrain data, to further improve the simulation accuracy of surface water depth and inundation range.

[0083] In a possible embodiment, determining the water accumulation depth of each grid cell according to the water accumulation change, the initial water depth and the grid cell area includes:

[0084] Determine the amount of change in ponding depth according to the amount of change in ponding of each grid unit and the area of ​​the grid unit;

[0085] The water accumulation depth of each grid unit is determined according to the sum of the water accumulation depth change and the initial water depth.

[0086] In a possible embodiment, after data interaction is performed between the one-dimensional drainage pipe network model and the two-dimensional surface inundation model according to the overflow data and the backflow data of the plurality of pipe network nodes, the method further includes:

[0087] Acquire at least one simulated hydrological data according to at least one simulated event;

[0088] Calculate the simulation accuracy of the one-dimensional drainage network model and the two-dimensional surface inundation model in the at least one simulation event according to the at least one simulated hydrological data and the at least one measured hydrological data, wherein the simulation accuracy includes a Nash efficiency coefficient and a determination coefficient;

[0089] The calculation formulas of the Nash efficiency coefficient and the determination coefficient are as follows:

[0090]

[0091]

[0092] Where NSE is the Nash efficiency coefficient, R 2 is the coefficient of determination, X m,i is the measured hydrological data of the ith measured session, X s,i is the simulated hydrological data of the ith measured event, X m is the average value of the measured hydrological data of all measured sessions, X s is the average value of the simulated hydrological data of all measured events, n is the total number of measured events, and the value range of NSE is (-∞,1], where the closer NSE is to 1, the higher the accuracy.

[0093] The validity of the at least one simulated hydrological data is determined according to the Nash efficiency coefficient and the determination coefficient.

[0094] In a possible embodiment, determining the validity of the at least one simulated hydrological data according to the Nash efficiency coefficient and the determination coefficient includes:

[0095] If it is detected that the Nash efficiency coefficient and the determination coefficient are both greater than a preset threshold, and the relative error between the measured hydrological data and the simulated hydrological data is less than a preset error, then the at least one simulated hydrological data is determined to be valid data; and,

[0096] If it is detected that the Nash efficiency coefficient and the determination coefficient are not both greater than the preset threshold, and / or the relative error between the measured hydrological data and the simulated hydrological data is not less than the preset error, then it is determined that at least one simulated hydrological data is invalid data, and the step of "performing data interaction between the one-dimensional drainage network model and the two-dimensional surface inundation model based on the overflow data and backflow data of the multiple network nodes" is re-executed.

[0097] Among them, the relative error (ε) refers to the ratio of the absolute error (E) to the measured hydrological data. The absolute error (E) refers to the absolute value of the difference between the measured hydrological data and the simulated hydrological data. The calculation formula of the relative error (ε) is ε = E / measured hydrological data = |measured hydrological data-simulated hydrological data| / measured hydrological data.

[0098] The preset threshold is 0.60 and the preset error is 15%. Specifically, when the Nash efficiency coefficient and the determination coefficient are both greater than 0.60 and the relative error is less than 15%, the model accuracy is determined to be high and the simulated hydrological data is valid data. Otherwise, the simulated hydrological data is determined to be invalid data, and the simulation operation is repeated or the model parameters are adjusted until the Nash efficiency coefficient and the determination coefficient are both greater than 0.60 and the relative error is less than 15%.

[0099] The measured hydrological data may include measured water depth and measured flooding range, and the simulated hydrological data may include simulated water depth and simulated flooding range.

[0100] It can be seen that in the embodiments of the present application, based on the measured hydrological data and the simulated hydrological data, the model is calibrated and verified through the Nash efficiency coefficient, the determination coefficient and the relative error. When the preset conditions are met, it means that the model accuracy is high and the obtained simulated hydrological data has a high accuracy. When the preset conditions are not met, it means that the model accuracy is low and the obtained simulated hydrological data has a low accuracy. It is necessary to re-simulate or adjust the model parameters and repeat the simulation process until the calculated Nash efficiency coefficient, the determination coefficient and the relative error meet the preset conditions.

[0101] It can be seen that in the embodiment of the present application, by bidirectionally coupling the one-dimensional drainage network model and the two-dimensional surface flooding model, bidirectional flow exchange can be performed through overflow and backflow of the network nodes. The two-dimensional surface flooding model feeds back the backflow of the network nodes to the one-dimensional drainage network model and performs a one-dimensional network simulation, outputs the overflow data of the network nodes, and updates the input file of the two-dimensional surface flooding model. In this way, compared with the existing one-way coupling method in which the water flow between models can only overflow from the node to the surface, and the surface water flow cannot return to the surface through the node, the surface water flow and the underground pipe flow in the bidirectional coupling method provided in the embodiment of the present application can be exchanged with each other through the node, which can accurately simulate the overflow and retreat process of surface water accumulation, which is conducive to improving the accuracy and comprehensiveness of urban waterlogging simulation.

[0102] See also Figure 3 , Figure 3 The present application provides a schematic diagram of a simulation process of ground surface water effect coupled to an urban pipe network, such as Figure 3 As shown, the ground water accumulation process includes ground water overflowing and ground water receding.

[0103] Among them, the process of surface water overflow refers to the process from the beginning of surface water accumulation to the maximum water depth, and from the beginning of surface flooding to the maximum flooding range; and the process of surface water receding refers to the process from the maximum water depth of surface water accumulation to the end of ...

[0104] in, Figure 3 3-1 in the figure shows the overflow process of ground water, i.e. overflow occurs, where the arrow direction indicates the direction of water flow. Specifically, the surface water level is lower than the water level of the pipe network node, and the pipe network water flows into the ground through the pipe network node, i.e., the rainwater inlet or the inspection well mouth, and merges into the ground water flow.

[0105] When ground water overflows, the calculation formula for the overflow volume of the pipe network node is as follows:

[0106]

[0107] Among them, Q0 is the overflow of the pipe network node, C0 is the orifice flow coefficient, and its value range is [0,1]; A m is the water flow area of ​​the network node, g is the gravitational acceleration, h 1d is the water level at the network node, h 2d The surface water level.

[0108] in, Figure 3 3-2 in Figure 3 shows a situation in the process of ground water receding, that is, backflow occurs, where the arrow direction indicates the direction of water flow. Specifically, the relationship between the surface water level, the ground elevation and the water level of the pipe network node is: the water level of the pipe network node < ground elevation < surface water level. The surface water flows through the pipe network node, that is, the rainwater outlet or the inspection well mouth, into the underground pipe and merges into the underground pipe network water flow.

[0109] When the ground water recedes, the calculation formula of the backflow volume of the pipe network node is as follows:

[0110]

[0111] Among them, Q s is the return flow of the network node, c w is the weir flow coefficient, ranging from [0,1]; w is the circumference of the network node or the width of the stormwater inlet.

[0112] in, Figure 3Figure 3-3 shows another situation of the ground water receding process, that is, backflow occurs, where the arrow direction indicates the direction of water flow. Specifically, the relationship between the surface water level, the ground elevation and the water level of the pipe network node is: ground elevation < water level of the pipe network node < surface water level. The surface water flows through the pipe network node, that is, the rainwater outlet or inspection well mouth, into the underground pipe and merges into the underground pipe network water flow.

[0113] When another situation occurs where the ground water recedes, the calculation formula for the return flow of the pipe network node is as follows:

[0114]

[0115] It can be seen that in this embodiment, when the surface water level is lower than the network node water level, overflow calculation is performed, and water flows from the network to the surface; when the surface water level is higher than the network node water level, backflow calculation is performed, and water flows from the surface to the network; when the surface water level is equal to the network node water level, no flow exchange is performed.

[0116] See also Figure 4 , Figure 4 The embodiment of the present application provides a one-dimensional drainage network model and a two-dimensional surface flooding model coupling node diagram, such as Figure 4 As shown, the coupling node diagram is a coupling node diagram of a one-dimensional drainage network model and a two-dimensional surface inundation model in the study area.

[0117] Among them, the one-dimensional drainage network model is used to simulate the water flow of the underground network, and the two-dimensional surface inundation model is used to simulate the surface runoff. The one-dimensional drainage network model and the two-dimensional surface inundation model are bidirectionally coupled through multiple network nodes, such as Figure 4 As shown, the pipe network node can be a rainwater wellhead, and multiple pipe network nodes are exchange nodes for surface water flow and underground pipe flow.

[0118] The study area includes drainage pipes. This embodiment only shows part of the drainage pipes in the study area, and the drainage pipes include multiple pipe network nodes.

[0119] Furthermore, the two-dimensional surface inundation model divides the study area into multiple grid units. The user can define the shape and size of the grid units according to the actual situation and the needs of the model. This embodiment does not limit the shape and size of the grid units, but only provides a case where multiple grid units are the same square grid. Each pipe network node corresponds to a grid unit, and the pipe network node and the corresponding grid unit can be defined as a coupling node and a coupling unit, respectively. The one-dimensional drainage pipe network model and the two-dimensional surface inundation model are specifically bidirectionally coupled through coupling nodes and coupling units.

[0120] It can be seen that in the embodiment of the present application, by bidirectionally coupling the one-dimensional drainage network model and the two-dimensional surface hydrodynamic model, the surface water flow and the underground pipe flow can be exchanged with each other through nodes, which can accurately simulate the overflow and receding process of surface water, which is beneficial to improving the accuracy of the simulation of the surface water action process.

[0121] See also Figure 5 , Figure 5 This is a functional unit block diagram of a ground surface water action process simulation device coupled to a city pipe network provided in an embodiment of the present application. The ground surface water action process simulation device 5 coupled to a city pipe network includes:

[0122] An acquisition unit 501 is used to acquire underlying surface data, pipe network data, rainfall data and measured hydrological data of a research area;

[0123] Processing unit 502 is used to construct a one-dimensional drainage pipe network model of the study area according to the underlying surface data, the pipe network data and the rainfall data; and, to construct a two-dimensional surface inundation model of the study area according to the underlying surface data and the rainfall data, wherein the one-dimensional drainage pipe network model is bidirectionally coupled with the two-dimensional surface inundation model through multiple pipe network nodes, and the multiple pipe network nodes are exchange nodes for surface water flow and underground pipe flow; and, to determine at least one simulation session and the total simulation time corresponding to each simulation session according to the rainfall data and the measured hydrological data; and, based on a unified time step, to perform data interaction between the one-dimensional drainage pipe network model and the two-dimensional surface inundation model according to the overflow data and backflow data of the multiple pipe network nodes until the total simulation time of each simulation session is reached, so as to complete the simulation event of the at least one simulation session of the overflow and receding process of the surface water in the study area.

[0124] It can be seen that in the embodiment of the present application, the processor obtains the underlying surface data, pipe network data, rainfall data and measured hydrological data of the study area; constructs a one-dimensional drainage pipe network model of the study area according to the underlying surface data, pipe network data and rainfall data; and constructs a two-dimensional surface inundation model of the study area according to the underlying surface data and rainfall data, the one-dimensional drainage pipe network model and the two-dimensional surface inundation model are bidirectionally coupled through multiple pipe network nodes, and the multiple pipe network nodes are exchange nodes for surface water flow and underground pipe flow; determines at least one simulation session and the total simulation time corresponding to each simulation session according to the rainfall data and the measured hydrological data; based on a preset time step, data interaction is performed between the one-dimensional drainage pipe network model and the two-dimensional surface inundation model according to the overflow data and backflow data of multiple pipe network nodes until the total simulation time of each simulation session is reached, so as to complete the simulation event of the at least one simulation session of the overflow and receding process of the surface water in the study area. In this way, compared with the existing urban waterlogging simulation scheme that is unidirectionally connected through a one-dimensional drainage network model and a two-dimensional surface inundation model, in the present application, the one-dimensional drainage network model and the two-dimensional surface inundation model are bidirectionally coupled, and bidirectional flow exchange can be performed through overflow and backflow of the network nodes, thereby achieving accurate simulation of surface water overflow and water retreat process, which is conducive to improving the accuracy and comprehensiveness of urban waterlogging simulation.

[0125] In some embodiments, the underlying surface data includes terrain data, elevation data and land use type, and the land use type includes construction land, green land, transportation land and public land;

[0126] The pipe network data includes rainwater pipe information, drain outlet information and node information. The rainwater pipe information includes the spatial distribution of the pipe network, pipe diameter, pipe section length, pipe roughness, and upstream and downstream bottom elevations. The drain outlet information includes the coordinates, elevation, and drainage type information of the drain outlet. The node information includes the node bottom elevation, well depth, and coordinate information. The rainwater pipe information, the drain outlet information, and the node information need to be restricted to the same coordinate system.

[0127] The rainfall data includes at least one rainfall amount of a measured event, and there is a corresponding relationship between the measured event and the simulated event;

[0128] The measured hydrological data include measured runoff data of pipelines and node monitoring water levels during measured rainfall events. The measured hydrological data correspond to the measured events of rainfall data, and the test duration of the measured events is the same as the total simulation duration of the corresponding simulation events.

[0129] In some embodiments, the data interaction between the one-dimensional drainage pipe network model and the two-dimensional surface flooding model according to the overflow data and the backflow data of the plurality of pipe network nodes includes:

[0130] In a first time step, performing surface water flow simulation by using the surface two-dimensional flooding model to determine backflow data of the plurality of pipe network nodes in the first time step, wherein the backflow data includes backflow volumes of the plurality of pipe network nodes;

[0131] Sending the backflow data to the one-dimensional drainage pipe network model to drive the one-dimensional drainage pipe network model to perform underground pipe flow simulation in a second time step;

[0132] In the second time step, the underground pipe flow simulation is performed by the one-dimensional drainage pipe network model to determine overflow data of the plurality of pipe network nodes in the second time step, wherein the overflow data includes overflow amounts of the plurality of pipe network nodes;

[0133] The input file of the surface two-dimensional flooding model is updated according to the overflow data, and the surface two-dimensional flooding model is driven according to the input file to repeat the surface water flow simulation and determine the backflow data of the multiple pipe network nodes in the third time step within a third time step.

[0134] In some embodiments, the method further comprises: performing surface water flow simulation according to the surface two-dimensional flooding model to determine the return flow data of the plurality of pipe network nodes within the first time step;

[0135] In a fourth time step, an underground pipe flow simulation is performed through the one-dimensional drainage pipe network model to determine the coordinates, water depths and flow rates of the pipe network nodes where overflow or backflow occurs in the fourth time step;

[0136] The surface two-dimensional inundation model is driven to perform the surface water flow simulation within the first time step according to the input file of the surface two-dimensional inundation model and the flow of the pipe network node, and the flow of the pipe network node is used as the input boundary condition for the surface two-dimensional inundation model to perform the surface water flow simulation.

[0137] In some embodiments, when the water level of the pipe network node is greater than the surface water level, overflow calculation is performed to obtain the overflow volume of the plurality of pipe network nodes, and, when the water level of the pipe network node is less than the surface water level, backflow calculation is performed to obtain the backflow volume of the plurality of pipe network nodes, and, when the water level of the pipe network node is equal to the surface water level, flow exchange is not performed, wherein,

[0138] The calculation formula for overflow and return flow of each pipe network node is as follows:

[0139]

[0140]

[0141] Among them, Q0 is the overflow of the pipe network node, Q s is the return flow of the network node, C0 is the orifice flow coefficient, and its value range is [0,1]; A m is the water flow area of ​​the network node, g is the gravitational acceleration, h 1d is the water level at the network node, h 2d is the surface water level, h z is the ground elevation, c w is the weir flow coefficient, ranging from [0,1]; w is the circumference of the network node or the width of the stormwater inlet.

[0142] In some embodiments, the study area in the two-dimensional surface inundation model includes a plurality of closely arranged grid cells, the backflow data also includes the water depth of each grid cell, and the surface water flow simulation is performed by the two-dimensional surface inundation model to determine the backflow data of the plurality of pipe network nodes within the first time step, and the method further includes:

[0143] Get the initial water depth of each grid cell;

[0144] Determining the change in water accumulation in each grid cell after the first time step;

[0145] The calculation formula for the water flow change of each grid unit is as follows:

[0146] Q grid =Q up +Q down +Q left +Q right (3)

[0147]

[0148] Among them, Q grid is the change in water accumulation in the grid unit, Q up , Q down , Q left , Q right are the flow rates from the upstream, downstream, left and right grid cells respectively, which are obtained by formula (4); Q ij is the flow at the junction of adjacent grid cells i and j, A ij is the cross-sectional area at the junction of adjacent grid cells i and j, R ij is the hydraulic radius at the junction of adjacent grid cells i and j, S ij is the water surface slope between adjacent grid cells i and j, and n is the Manning coefficient;

[0149] The water accumulation depth of each grid unit is determined according to the water accumulation change, the initial water accumulation depth and the area of ​​each grid unit.

[0150] In some embodiments, after data is exchanged between the one-dimensional drainage pipe network model and the two-dimensional surface flooding model according to the overflow data and the backflow data of the plurality of pipe network nodes, the method further includes:

[0151] Acquire at least one simulated hydrological data according to at least one simulated event;

[0152] Calculate the simulation accuracy of the one-dimensional drainage network model and the two-dimensional surface inundation model in the at least one simulation event according to the at least one simulated hydrological data and the at least one measured hydrological data, wherein the simulation accuracy includes a Nash efficiency coefficient and a determination coefficient;

[0153] The calculation formulas of the Nash efficiency coefficient and the determination coefficient are as follows:

[0154]

[0155]

[0156] Where NSE is the Nash efficiency coefficient, R 2 is the coefficient of determination, X m,i is the measured hydrological data of the ith measured session, X s,i is the simulated hydrological data of the ith measured event, X m is the average value of the measured hydrological data of all measured sessions, X s is the average value of the simulated hydrological data of all measured events, n is the total number of measured events, and the value range of NSE is (-∞,1];

[0157] The validity of the at least one simulated hydrological data is determined according to the Nash efficiency coefficient and the determination coefficient.

[0158] In some embodiments, determining the validity of the at least one simulated hydrological data according to the Nash efficiency coefficient and the determination coefficient comprises:

[0159] If it is detected that the Nash efficiency coefficient and the determination coefficient are both greater than a preset threshold, and the relative error between the measured hydrological data and the simulated hydrological data is less than a preset error, then the at least one simulated hydrological data is determined to be valid data; and,

[0160] If it is detected that the Nash efficiency coefficient and the determination coefficient are not both greater than the preset threshold, and / or the relative error between the measured hydrological data and the simulated hydrological data is not less than the preset error, then it is determined that at least one simulated hydrological data is invalid data, and the step of "performing data interaction between the one-dimensional drainage network model and the two-dimensional surface inundation model based on the overflow data and backflow data of the multiple network nodes" is re-executed.

[0161] An embodiment of the present application provides a computer-readable storage medium having a computer program / instruction stored thereon. When the computer program / instruction is executed by a processor, the steps of the method described in any possible embodiment are implemented.

[0162] It should be noted that, for the aforementioned method embodiments, for the sake of simplicity, they are all expressed as a series of action combinations, but those skilled in the art should be aware that the present application is not limited by the described order of actions, because according to the present application, certain steps can be performed in other orders or simultaneously. Secondly, those skilled in the art should also be aware that the embodiments described in the specification are all preferred embodiments, and the actions and modules involved are not necessarily required by the present application.

[0163] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0164] In the several embodiments provided in the present application, it should be understood that the disclosed device can be implemented in other ways. For example, the device embodiments described above are only schematic, such as the division of the above-mentioned units, which is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, and the indirect coupling or communication connection of devices or units can be electrical or other forms.

[0165] The units described above as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0166] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0167] If the above-mentioned integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable memory. Based on this understanding, the technical solution of the present application is essentially or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a memory, including a number of instructions to enable a computer device (which can be a personal computer, server or network device, etc.) to perform all or part of the steps of the above-mentioned methods of each embodiment of the present application. The aforementioned memory includes: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, disk or CD-ROM and other media that can store program codes.

[0168] A person skilled in the art can understand that all or part of the steps in the various methods of the above embodiments can be completed by instructing related hardware through a program, and the program can be stored in a computer-readable memory, and the memory can include: a flash drive, a read-only memory (English: Read-Only Memory, abbreviated as: ROM), a random access memory (English: Random Access Memory, abbreviated as: RAM), a magnetic disk or an optical disk, etc.

[0169] The embodiments of the present application are introduced in detail above. Specific examples are used in this article to illustrate the principles and implementation methods of the present application. The description of the above embodiments is only used to help understand the method and core idea of ​​the present application. At the same time, for general technical personnel in this field, according to the idea of ​​the present application, there will be changes in the specific implementation method and application scope. In summary, the content of this specification should not be understood as a limitation on the present application.

Claims

1. A method for simulating the surface water action process coupled to an urban pipe network, characterized in that: The method comprises: Obtain the underlying surface data, pipe network data, rainfall data and measured hydrological data of the study area; A one-dimensional drainage pipe network model of the study area is constructed according to the underlying surface data, the pipe network data and the rainfall data; and a two-dimensional surface inundation model of the study area is constructed according to the underlying surface data and the rainfall data, wherein the one-dimensional drainage pipe network model is bidirectionally coupled with the two-dimensional surface inundation model through a plurality of pipe network nodes, and the plurality of pipe network nodes are exchange nodes for surface water flow and underground pipe flow; Determine at least one simulation session and a total simulation duration corresponding to each simulation session according to the rainfall data and the measured hydrological data; Based on a preset time step, data interaction is performed between the one-dimensional drainage network model and the two-dimensional surface inundation model according to the overflow data and backflow data of the multiple network nodes until the total simulation time of each simulation session is reached, so as to complete the simulation events of at least one simulation session of the overflow and receding process of ground water in the study area.

2. The method according to claim 1, characterized in that The underlying surface data includes terrain data, elevation data and land use types, and the land use types include construction land, green land, transportation land and public land; The pipe network data includes rainwater pipe information, drain outlet information and node information. The rainwater pipe information includes the spatial distribution of the pipe network, pipe diameter, pipe section length, pipe roughness, and upstream and downstream bottom elevations. The drain outlet information includes the coordinates, elevation, and drainage type information of the drain outlet. The node information includes the node bottom elevation, well depth, and coordinate information. The rainwater pipe information, the drain outlet information, and the node information need to be restricted to the same coordinate system. The rainfall data includes at least one rainfall amount of a measured event, and there is a corresponding relationship between the measured event and the simulated event; The measured hydrological data include measured runoff data of pipelines and node monitoring water levels during measured rainfall events. The measured hydrological data correspond to the measured events of rainfall data, and the test duration of the measured events is the same as the total simulation duration of the corresponding simulation events.

3. The method according to claim 2, characterized in that The data interaction between the one-dimensional drainage pipe network model and the two-dimensional surface flooding model according to the overflow data and the backflow data of the plurality of pipe network nodes includes: In a first time step, performing surface water flow simulation by using the surface two-dimensional flooding model to determine backflow data of the plurality of pipe network nodes in the first time step, wherein the backflow data includes backflow volumes of the plurality of pipe network nodes; Sending the backflow data to the one-dimensional drainage pipe network model to drive the one-dimensional drainage pipe network model to perform underground pipe flow simulation in a second time step; In the second time step, the underground pipe flow simulation is performed by the one-dimensional drainage pipe network model to determine overflow data of the plurality of pipe network nodes in the second time step, wherein the overflow data includes overflow amounts of the plurality of pipe network nodes; The input file of the surface two-dimensional flooding model is updated according to the overflow data, and the surface two-dimensional flooding model is driven according to the input file to repeat the surface water flow simulation and determine the backflow data of the multiple pipe network nodes in the third time step within a third time step.

4. The method according to claim 3, characterized in that said performing surface water flow simulation according to said surface two-dimensional inundation model to determine the backflow data of said plurality of pipe network nodes in said first time step before; The method further comprises: In a fourth time step, an underground pipe flow simulation is performed through the one-dimensional drainage pipe network model to determine the coordinates, water depths and flow rates of the pipe network nodes where overflow or backflow occurs in the fourth time step; The surface two-dimensional inundation model is driven to perform the surface water flow simulation within the first time step according to the input file of the surface two-dimensional inundation model and the flow of the pipe network node, and the flow of the pipe network node is used as the input boundary condition for the surface two-dimensional inundation model to perform the surface water flow simulation.

5. The method according to claim 4, characterized in that When the water level of the pipe network node is greater than the surface water level, overflow calculation is performed to obtain the overflow volume of the plurality of pipe network nodes, and when the water level of the pipe network node is less than the surface water level, backflow calculation is performed to obtain the backflow volume of the plurality of pipe network nodes, and when the water level of the pipe network node is equal to the surface water level, flow exchange is not performed, wherein, The calculation formula for overflow and return flow of each pipe network node is as follows: Among them, Q0 is the overflow of the pipe network node, Q s is the return flow of the network node, C0 is the orifice flow coefficient, and its value range is [0,1]; A m is the water flow area of ​​the network node, g is the gravitational acceleration, h 1d is the water level at the network node, h 2d is the surface water level, h z is the ground elevation, c w is the weir flow coefficient, ranging from [0,1]; w is the circumference of the network node or the width of the stormwater inlet.

6. The method according to any one of claims 3 to 5, characterized in that: The study area in the two-dimensional surface inundation model includes a plurality of closely arranged grid cells, the backflow data also includes the water depth of each grid cell, the surface water flow simulation is performed by the two-dimensional surface inundation model to determine the backflow data of the plurality of pipe network nodes within the first time step, and the method further includes: Get the initial water depth of each grid cell; Determining the change in water accumulation in each grid cell after the first time step; The calculation formula for the water flow change of each grid unit is as follows: Q grid =Q up +Q down +Q left +Q right (3) Among them, Q grid is the change in water accumulation in the grid unit, Q up , Q down , Q left , Q right are the flow rates from the upstream, downstream, left and right grid cells respectively, which are obtained by formula (4); Q ij is the flow at the junction of adjacent grid cells i and j, A ij is the cross-sectional area at the junction of adjacent grid cells i and j, R ij is the hydraulic radius at the junction of adjacent grid cells i and j, S ij is the water surface slope between adjacent grid cells i and j, and n is the Manning coefficient; The water accumulation depth of each grid unit is determined according to the water accumulation change, the initial water accumulation depth and the area of ​​each grid unit.

7. The method according to claim 6, characterized in that After data is exchanged between the one-dimensional drainage pipe network model and the two-dimensional surface flooding model according to the overflow data and the backflow data of the plurality of pipe network nodes, the method further includes: Acquire at least one simulated hydrological data according to at least one simulated event; Calculate the simulation accuracy of the one-dimensional drainage network model and the two-dimensional surface inundation model in the at least one simulation event according to the at least one simulated hydrological data and the at least one measured hydrological data, wherein the simulation accuracy includes a Nash efficiency coefficient and a determination coefficient; The calculation formulas of the Nash efficiency coefficient and the determination coefficient are as follows: Where NSE is the Nash efficiency coefficient, R 2 is the coefficient of determination, X m,i is the measured hydrological data of the ith measured session, X s,i is the simulated hydrological data of the ith measured event, X m is the average value of the measured hydrological data of all measured sessions, X s is the average value of the simulated hydrological data of all measured events, n is the total number of measured events, and the value range of NSE is (-∞,1]; The validity of the at least one simulated hydrological data is determined according to the Nash efficiency coefficient and the determination coefficient.

8. The method according to claim 7, characterized in that The determining the validity of the at least one simulated hydrological data according to the Nash efficiency coefficient and the determination coefficient comprises: If it is detected that the Nash efficiency coefficient and the determination coefficient are both greater than a preset threshold, and the relative error between the measured hydrological data and the simulated hydrological data is less than a preset error, then the at least one simulated hydrological data is determined to be valid data; and, If it is detected that the Nash efficiency coefficient and the determination coefficient are not both greater than the preset threshold value, and / or the relative error between the measured hydrological data and the simulated hydrological data is not less than the preset error, it is determined that at least one simulated hydrological data is invalid data, and the step of "performing data interaction between the one-dimensional drainage network model and the two-dimensional surface inundation model based on the overflow data and backflow data of the multiple network nodes" is re-executed.

9. A device for simulating the surface water action process coupled to an urban pipe network, characterized in that: include: The acquisition unit is used to obtain the underlying surface data, pipe network data, rainfall data and measured hydrological data of the study area; A processing unit, configured to construct a one-dimensional drainage pipe network model of the study area according to the underlying surface data, the pipe network data and the rainfall data; and to construct a two-dimensional surface inundation model of the study area according to the underlying surface data and the rainfall data, wherein the one-dimensional drainage pipe network model is bidirectionally coupled with the two-dimensional surface inundation model through a plurality of pipe network nodes, and the plurality of pipe network nodes are exchange nodes for surface water flow and underground pipe flow; And, determining at least one simulation session and a total simulation duration corresponding to each simulation session according to the rainfall data and the measured hydrological data; Furthermore, based on a uniform time step, data interaction is performed between the one-dimensional drainage network model and the two-dimensional surface inundation model according to the overflow data and backflow data of the multiple network nodes until the total simulation time of each simulation session is reached, so as to complete the simulation events of at least one simulation session of the ground water overflow and receding process in the study area.

10. An electronic device, characterized in that: The method comprises a processor and a memory, wherein the memory stores a computer program, and the processor executes the method according to any one of claims 1 to 8 when calling the computer program in the memory.

Images