Regional inundation deduction method and device, electronic equipment and storage medium

Abstract

The invention relates to a regional inundation deduction method and device, electronic equipment and a storage medium, and the method comprises the steps: determining an overflow path at a first moment based on a two-dimensional model of a river channel and a pipeline according to a first parameter, corresponding to the first moment, of water flow at a plurality of positions in a pipeline and / or a river channel in a target region; based on a surface water power two-dimensional model, determining the surface water accumulation range and depth at the first moment according to second parameters of water flow at multiple positions of the surface of the target area corresponding to the first moment, and based on the two-dimensional grids and the second parameters of the two-dimensional grids, determining third parameters corresponding to the three-dimensional grids of each sub-area in the sub-area water power three-dimensional model; determining the water accumulation range and depth of the sub-region at the first moment based on the third parameter; and deducing the submerging process of the target area based on the overflow paths, the surface ponding ranges and depths, and the subarea ponding ranges and depths at the plurality of first moments. According to the method disclosed by the invention, the efficiency and precision of large-range regional inundation deduction can be improved at the same time.

Description

TECHNICAL FIELD

[0001] The present disclosure relates to the technical field of disaster simulation, and particularly relates to a regional inundation deduction method and device, electronic equipment and storage medium. BACKGROUND

[0002] Under the background of climate change intensification and urbanization process superposition, water disasters such as floods and rainstorms show high frequency and high intensity of spatiotemporal evolution characteristics.

[0003] The current water disaster simulation technology system has significant technical bottlenecks: although the two-dimensional water power model has the advantage of calculation efficiency, it cannot represent the refined water flow dynamics of three-dimensional key disaster-bearing bodies such as urban underground space, important infrastructure and immovable cultural relics due to the limitation of plane simplification assumption; and the three-dimensional water power model can realize refined reconstruction of physical field, but it faces problems such as complex grid division and huge consumption of computing resources, and the simulation of the inundation process of the target region for a single time can take tens of hours.

[0004] Therefore, the efficiency and accuracy of large-scale regional inundation deduction cannot be considered at the same time. SUMMARY

[0005] Therefore, the efficiency and accuracy of large-scale regional inundation deduction cannot be considered at the same time.

[0006] According to an aspect of the present disclosure, a regional inundation deduction method is provided, which comprises:

[0007] Based on the two-dimensional model of the river and the pipeline, the first parameter of the water flow at a plurality of positions in the pipeline and / or the river in the target region corresponding to the first time is determined to determine the overflow path corresponding to the first time; based on the two-dimensional model of the surface water power, the second parameter of the water flow at a plurality of positions on the surface in the target region corresponding to the first time is determined to determine the surface water accumulation range and depth corresponding to the first time, the two-dimensional grid of the surface is included in the two-dimensional model of the surface water power, and each two-dimensional grid corresponds to at least one second parameter; based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, the third parameter corresponding to each sub-area three-dimensional grid in the sub-area water power three-dimensional model is determined; based on the third parameter, the sub-area water accumulation range and depth corresponding to the first time are determined; based on the overflow path, the surface water accumulation range and depth, and the sub-area water accumulation range and depth corresponding to each of the plurality of first times, the inundation process of the target region is deduced.

[0008] In a possible implementation, the method further includes: determining a coupling region based on the two-dimensional grid and the three-dimensional grid of the sub-region, the coupling region representing a connection region of the ground surface and the sub-region; obtaining an intermediate three-dimensional grid representing the coupling region, a grid precision of the intermediate three-dimensional grid being between a grid precision of the two-dimensional grid and a grid precision of the three-dimensional grid of the sub-region; determining the third parameter corresponding to each three-dimensional grid of the sub-region based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid includes: determining a transition parameter corresponding to the intermediate three-dimensional grid based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid; determining the third parameter corresponding to at least part of the three-dimensional grid of the sub-region based on the transition parameter; and performing hydrodynamic simulation based on the third parameter corresponding to at least part of the three-dimensional grid of the sub-region and the three-dimensional hydrodynamic model of the sub-region to determine the third parameter corresponding to each three-dimensional grid of the sub-region.

[0009] In a possible implementation, the method of determining the transition parameter corresponding to the intermediate three-dimensional grid based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid includes: determining an initial transition parameter corresponding to the intermediate three-dimensional grid intersecting with the two-dimensional grid based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, and taking the intermediate three-dimensional grid intersecting with the two-dimensional grid as a first intersecting three-dimensional grid; and recursively assigning a value to each intermediate three-dimensional grid in a direction from near to far of the first intersecting three-dimensional grid based on the initial transition parameter to obtain each transition parameter.

[0010] In a possible implementation, the method of determining the third parameter corresponding to the three-dimensional grid of the sub-region intersecting with the intermediate three-dimensional grid based on the transition parameter includes: determining the third parameter corresponding to the three-dimensional grid of the sub-region intersecting with the intermediate three-dimensional grid based on the transition parameter; and the method further includes: taking the three-dimensional grid of the sub-region intersecting with the intermediate three-dimensional grid as a second intersecting three-dimensional grid; and taking the third parameter corresponding to the second intersecting three-dimensional grid as a hydrodynamic initial parameter of the three-dimensional hydrodynamic model of the sub-region.

[0011] In a possible implementation, the method of determining the third parameter corresponding to the three-dimensional grid of the sub-region intersecting with the intermediate three-dimensional grid based on the transition parameter includes: determining at least one intermediate three-dimensional grid closest to a single second intersecting three-dimensional grid for the single second intersecting three-dimensional grid; and determining an initial parameter corresponding to the single second intersecting three-dimensional grid based on the transition parameter corresponding to the at least one intermediate three-dimensional grid.

[0012] In a possible implementation, the initial transition parameter includes a first flow rate, and the determining, based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, of the initial transition parameter corresponding to an intermediate three-dimensional grid intersecting with the two-dimensional grid includes: vertically assigning the first flow rate corresponding to the first intersecting three-dimensional grid to a plurality of intermediate three-dimensional grids under the first intersecting three-dimensional grid to obtain a second flow rate corresponding to each of the plurality of intermediate three-dimensional grids under the first intersecting three-dimensional grid; and the assigning, based on the initial transition parameter, of each of the intermediate three-dimensional grids in a direction from near to far of the first intersecting three-dimensional grid to obtain each of the transition parameters includes: interpolating each of the first flow rate and the second flow rate to determine the transition parameter, and the transition parameter includes a flow rate.

[0013] In a possible implementation, the initial transition parameter further includes a first water level, and the transition parameter further includes a pressure and a water phase volume fraction, and the method further includes: determining a liquid surface height based on the first water level, determining a water phase integral corresponding to the intermediate three-dimensional grid based on the liquid surface height and an elevation corresponding to the intermediate three-dimensional grid, and determining a pressure corresponding to the intermediate three-dimensional grid based on the liquid surface height, an initial pressure corresponding to the liquid surface height, and the elevation corresponding to the intermediate three-dimensional grid.

[0014] According to another aspect of the present disclosure, there is provided a regional inundation deduction device, which comprises:

[0015] An overflow path determination unit is configured to determine, based on a two-dimensional model of a river channel and a pipeline, an overflow path corresponding to a first time based on first parameters of water flow at a plurality of positions in the pipeline and / or the river channel in a target region at the first time.

[0016] A surface water accumulation range and depth determination unit is configured to determine, based on a two-dimensional surface hydrodynamic model, a surface water accumulation range and depth corresponding to the first time based on second parameters of water flow at a plurality of positions on the surface in the target region at the first time, the two-dimensional surface hydrodynamic model including a two-dimensional grid representing the surface, each of the two-dimensional grids corresponding to at least one second parameter.

[0017] A third parameter determination unit is configured to determine, based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, third parameters corresponding to three-dimensional grids in a sub-region hydrodynamic three-dimensional model.

[0018] A sub-region water accumulation range and depth determination unit is configured to determine, based on the third parameters, a sub-region water accumulation range and depth corresponding to the first time.

[0019] The target area submergence deduction unit is configured to deduce a submergence process of the target area based on the overflow path, the surface water accumulation range and depth, and the sub-area water accumulation range and depth corresponding to each of the first time points.

[0020] In a possible implementation, the apparatus further includes:

[0021] The coupling area determination unit is configured to determine a coupling area based on the two-dimensional grid and the sub-area three-dimensional grid, the coupling area representing a connection area between the surface and the sub-area.

[0022] An intermediate three-dimensional grid representing the coupling area is obtained, the grid precision of the intermediate three-dimensional grid being between the grid precision of the two-dimensional grid and the grid precision of the sub-area three-dimensional grid.

[0023] The third parameter determination unit is further configured to:

[0024] determine a transition parameter corresponding to the intermediate three-dimensional grid based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid.

[0025] determine a third parameter corresponding to at least part of the sub-area three-dimensional grid based on the transition parameter.

[0026] perform hydrodynamic simulation based on the third parameter corresponding to at least part of the sub-area three-dimensional grid and the sub-area hydrodynamic three-dimensional model to determine the third parameter corresponding to each of the sub-area three-dimensional grids.

[0027] In a possible implementation, the third parameter determination unit is further configured to:

[0028] determine an initial transition parameter corresponding to an intermediate three-dimensional grid intersecting with the two-dimensional grid based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, and take the intermediate three-dimensional grid intersecting with the two-dimensional grid as a first intersection three-dimensional grid.

[0029] recursively assign a value to each of the intermediate three-dimensional grids in a direction from near to far of the first intersection three-dimensional grid based on the initial transition parameter to obtain each of the transition parameters.

[0030] In a possible implementation, the third parameter determination unit is further configured to:

[0031] determine a third parameter corresponding to a sub-area three-dimensional grid intersecting with the intermediate three-dimensional grid based on the transition parameter.

[0032] The apparatus further includes:

[0033] a second interface three-dimensional grid determination unit, configured to determine a sub-region three-dimensional grid that interfaces with the intermediate three-dimensional grid as a second interface three-dimensional grid;

[0034] a hydrodynamics initial parameter determination unit, configured to determine a third parameter corresponding to the second interface three-dimensional grid as a hydrodynamics initial parameter of the sub-region hydrodynamic three-dimensional model.

[0035] In a possible implementation, the third parameter determination unit is further configured to:

[0036] for a single second interface three-dimensional grid, determine at least one intermediate three-dimensional grid closest to the single second interface three-dimensional grid;

[0037] determine an initial parameter corresponding to the single second interface three-dimensional grid based on a transition parameter corresponding to the at least one intermediate three-dimensional grid.

[0038] In a possible implementation, the initial transition parameter includes a first flow rate, and the third parameter determination unit is further configured to:

[0039] vertically assign a first flow rate corresponding to the first interface three-dimensional grid to a plurality of intermediate three-dimensional grids under the first interface three-dimensional grid to obtain a second flow rate corresponding to each of the plurality of intermediate three-dimensional grids under the first interface three-dimensional grid;

[0040] interpolate each of the first flow rate and the second flow rate to determine the transition parameter, the transition parameter including a flow rate.

[0041] In a possible implementation, the initial transition parameter further includes a first water level, and the transition parameter further includes a pressure and a water phase volume fraction, and the third parameter determination unit is further configured to:

[0042] determine a liquid surface height based on the first water level, and determine a water phase integral corresponding to the intermediate three-dimensional grid based on the liquid surface height and an elevation corresponding to the intermediate three-dimensional grid;

[0043] determine a pressure corresponding to the intermediate three-dimensional grid based on the liquid surface height, an initial pressure corresponding to the liquid surface height, and the elevation corresponding to the intermediate three-dimensional grid.

[0044] According to another aspect of the present disclosure, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory, the processor executes the computer program to implement the steps of the above method.

[0045] According to another aspect of the present disclosure, there is provided a non-transitory computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the above method.

[0046] According to another aspect of the present disclosure, there is provided a computer program product comprising a computer program, or a non-transitory computer readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above method.

[0047] According to the embodiments of the present disclosure, for each first time, the overflow path can be determined based on the first parameter, the river channel and the two-dimensional model of the pipeline; the surface water accumulation range and depth can be determined based on the second parameter and the two-dimensional model of the surface water dynamics; and the third parameter corresponding to the sub-area three-dimensional grid can be determined based on the second parameter corresponding to the two-dimensional grid. The third parameter is used to determine the sub-area water accumulation range and depth corresponding to each time. In this way, the two-dimensional model of water dynamics and the three-dimensional model of water dynamics can be coupled for each first time, so as to complete the dynamic coupling of the two-dimensional model of water dynamics and the three-dimensional model of water dynamics. Not only the efficiency of the simulation of the regional inundation process is improved, but also the water disaster evolution process of the urban scale is retained, and the micro-scale sub-area can be finely simulated, so that the accuracy of disaster deduction is significantly enhanced. Therefore, the method of the present disclosure can improve the efficiency and accuracy of the large-scale regional inundation deduction.

[0048] Other features and aspects of the present disclosure will become apparent from the following detailed description of exemplary embodiments, taken in conjunction with the accompanying drawings. BRIEF DESCRIPTION OF DRAWINGS

[0049] The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate exemplary embodiments, features, and aspects of the present disclosure and serve to explain the principles of the present disclosure.

[0050] Figure 1 A flowchart of a regional inundation deduction method provided by the embodiments of the present disclosure.

[0051] Figure 2 A structure diagram of a regional inundation deduction device provided by the embodiments of the present disclosure.

[0052] Figure 3 A structure diagram of an electronic device for regional inundation deduction provided by the embodiments of the present disclosure. DETAILED DESCRIPTION

[0053] Various exemplary embodiments, features, and aspects of the present disclosure will be described in detail below with reference to the accompanying drawings. The same reference signs in the drawings represent functionally identical or similar elements. Although various aspects of the embodiments are shown in the drawings, the drawings are not necessarily drawn to scale unless specifically indicated.

[0054] As used herein, the terms "comprise", "comprising", "have", "having", "include", "including", "contain", "containing", or variants thereof, are open-ended, and include one or more stated features, integers, elements, steps, components or functions but do not preclude the presence or addition of one or more other features, integers, elements, steps, components, functions or groups thereof.

[0055] When an element is referred to as being "connected", "coupled", "responsive", or "in communication" with, to or with one or more other elements, it can be directly connected, coupled, responsive, or in communication with the one or more other elements or can be connected, coupled, responsive, or in communication with the one or more other elements via one or more intervening elements.

[0056] Although the terms first, second, third, etc. can be used herein to describe various elements / operations, these elements / operations should not be limited by these terms. These terms are only used to distinguish one element / operation from another. Thus, a first element / operation in some embodiments could be termed a second element / operation in other embodiments without departing from the teachings of the present inventive concept.

[0057] The word "exemplary" is used herein to mean "serving as an example, instance, or illustration." Any implementation described herein as "exemplary" is not necessarily to be construed as preferred or advantageous over other implementations.

[0058] In addition, for the purpose of convenience and brevity, detailed descriptions of well-known devices, methods, procedures, components and circuits will not be described in detail since they would be apparent to one skilled in the art. The skilled person will understand that the present disclosure can be practiced without certain specific details, which are provided in the following detailed description for purposes of illustrating the present disclosure.

[0059] Under the dual influence of climate change and urbanization, the frequency of extreme rainfall and urban flood disasters is rising, posing significant challenges to the protection of immovable cultural relics, the safety of urban underground space, and the operation of critical infrastructure. Urban flood simulation mainly relies on two-dimensional hydrodynamic models for water dynamics and prediction, which have high computational efficiency but have the following technical shortcomings:

[0060] Two-dimensional hydrodynamic models (e.g., LISFLOOD-FP model) lack the ability to simulate the water dynamics of complex three-dimensional features, making it difficult to meet the needs of detailed simulation. Existing two-dimensional hydrodynamic models are based on elevation grids and simplified boundary settings, which cannot accurately simulate the impact of complex three-dimensional structures in urban underground space, immovable cultural relics, or critical infrastructure on water flow paths, water depth, and velocity distribution, thereby limiting the ability to identify risks in complex three-dimensional urban features, especially in the microscopic scale of subway facilities, tunnel sections, and cultural heritage units.

[0061] Although the hydrodynamic three-dimensional model (for example, OpenFOAM) has complex three-dimensional modeling capability, the running efficiency is low. The hydrodynamic three-dimensional model adopts a method based on computational fluid dynamics (CFD) to simulate the flow state of water. The CFD method is sensitive to grid quality and boundary conditions, which leads to the fact that the hydrodynamic three-dimensional model often needs to independently set the initial condition, and it is difficult to reproduce the real response process of the overall water system of the city. In view of this, the present disclosure proposes a regional inundation process deduction method, which can balance the efficiency and accuracy of large-scale inundation deduction.

[0062] Figure 1 The flowchart of the regional inundation process deduction method provided by the embodiments of the present disclosure is shown. The method comprises:

[0063] S11, based on the river and pipeline two-dimensional model, according to the first parameter of the water flow at a plurality of positions in the pipeline and / or river in the target region corresponding to the first time, determining the overflow path corresponding to the first time.

[0064] The river and pipeline two-dimensional model can be a hydrodynamic two-dimensional model for simulating the flow state of the river and / or pipeline. The first time can be any time, or a specified time. The target region can be a region to be simulated for the inundation process. The first parameter can include the river depth corresponding to a plurality of positions in the river, the river flow rate. The first parameter can include the pipeline water level corresponding to a plurality of positions in the pipeline, the pipeline flow rate.

[0065] In the embodiments of the present disclosure, the first parameter can be input to the river and pipeline two-dimensional model. The river and pipeline two-dimensional model can simulate the flow state in the river and / or pipeline based on the first parameter to obtain the overflow path corresponding to the first time. The overflow path can show the overflow position of the river and / or pipeline, and the flow path after the water flow flows out of the river and / or pipeline, etc.

[0066] S12, based on the surface hydrodynamic two-dimensional model, according to the second parameter of the water flow at a plurality of positions on the surface in the target region corresponding to the first time, determining the surface water accumulation range and depth corresponding to the first time, the surface hydrodynamic two-dimensional model comprising a two-dimensional grid representing the surface, each two-dimensional grid corresponding to at least one second parameter.

[0067] The surface hydrodynamic two-dimensional model can be a hydrodynamic two-dimensional model for simulating the flow state of the surface. The second parameter can include the surface water depth corresponding to a plurality of positions on the surface, the surface water flow rate.

[0068] The two-dimensional model of surface water can include a plurality of two-dimensional grids representing the surface. The two-dimensional grids can represent a two-dimensional structure of the surface. A single node of the two-dimensional grids can correspond to at least one second parameter representing a water depth and a water flow rate of the surface at a position corresponding to the node. Alternatively, a center of the two-dimensional grids can correspond to at least one second parameter representing a water depth and a water flow rate of the surface at a position corresponding to the center.

[0069] In the embodiments of the present disclosure, the second parameters can be input into a two-dimensional model of surface water dynamics. The two-dimensional model of surface water dynamics can simulate a flow state of the surface water based on the second parameters to obtain a range and a depth of accumulated water of the surface at the first time.

[0070] S13, based on the two-dimensional grids and the second parameters corresponding to the two-dimensional grids, determine third parameters corresponding to three-dimensional grids of each sub-region in a three-dimensional model of water dynamics of a sub-region.

[0071] The three-dimensional model of water dynamics of a sub-region can be a three-dimensional model of water dynamics for simulating a flow state around and / or inside a key geographical object. The sub-region can include at least one of the following geographical objects: immovable cultural relics, infrastructures, and urban three-dimensional spaces. The immovable cultural relics can include ancient architectural complexes, grotto sites, etc. The infrastructures can include subway tunnels and underground comprehensive pipe galleries. The urban three-dimensional spaces can include above-ground architectural complexes and underground transportation networks. The three-dimensional model of water dynamics of a sub-region can include a plurality of three-dimensional grids of sub-regions representing the sub-region. The three-dimensional grids of sub-regions can represent a three-dimensional structure of a surface or an interior of the sub-region.

[0072] The third parameters can include sub-region water depths and sub-region water flow rates at a plurality of positions around and / or inside the region. A single node of the three-dimensional grids of sub-regions can correspond to at least one third parameter representing a sub-region water depth and a sub-region water flow rate at a position around and / or inside the sub-region corresponding to the node. Alternatively, a center of the three-dimensional grids of sub-regions can correspond to at least one third parameter representing a sub-region water depth and a sub-region water flow rate at a position around and / or inside the sub-region corresponding to the center.

[0073] In the embodiments of the present disclosure, the third parameters can be calculated based on the second parameters.

[0074] For example, the third parameters corresponding to the three-dimensional grids of sub-regions can be determined based on distances between the three-dimensional grids of sub-regions and the two-dimensional grids, and the second parameters corresponding to the two-dimensional grids.

[0075] For example, the second parameters can be mapped to the three-dimensional grids of sub-regions based on the first mapping relationship and a first mapping loss condition to obtain the third parameters. The first mapping loss condition can represent a magnitude of increase or decrease of the second parameters with changes in distances and heights from the second grids.

[0076] The above is only an example, and the method of calculating the third parameter based on the second parameter is not limited in the embodiments of the present disclosure.

[0077] In S14, the sub-region waterlogging range and depth corresponding to the first time are determined based on the third parameter.

[0078] In the embodiments of the present disclosure, the third parameter can be input into the sub-region hydrodynamic three-dimensional model. The sub-region hydrodynamic three-dimensional model can simulate the water flow state around and / or inside the sub-region based on the third parameter to obtain the sub-region waterlogging range and depth corresponding to the first time.

[0079] In S15, the submergence process of the target region is deduced based on the overflow path, the surface waterlogging range and depth, and the sub-region waterlogging range and depth corresponding to each of the first time.

[0080] In the embodiments of the present disclosure, the sub-region hydrodynamic three-dimensional model can deduce the submergence process of the target region based on the overflow path, the surface waterlogging range and depth, and the sub-region waterlogging range and depth corresponding to each of the first time. Alternatively, the sub-region hydrodynamic three-dimensional model can deduce the submergence process of the target region based on the third parameter corresponding to each of the first time. The sub-region hydrodynamic three-dimensional model can visually present the submergence process of the target region.

[0081] In the embodiments of the present disclosure, for each first time, the overflow path can be determined based on the first parameter and the river and pipeline two-dimensional model; the surface waterlogging range and depth can be determined based on the second parameter and the surface hydrodynamic two-dimensional model; and the third parameter corresponding to the two-dimensional grid can be determined based on the second parameter corresponding to the two-dimensional grid; and the sub-region waterlogging range and depth corresponding to each time can be determined using the third parameter. In this way, the hydrodynamic two-dimensional model and the hydrodynamic three-dimensional model can be coupled for each first time to complete the dynamic coupling of the hydrodynamic two-dimensional model and the hydrodynamic three-dimensional model. Not only the efficiency of the simulation of the regional submergence process is improved, but also the water disaster evolution process at the urban scale is retained, and the fine simulation of the sub-region at the microscopic scale is performed, which significantly enhances the accuracy of disaster deduction. Therefore, the method of the present disclosure can improve the efficiency and accuracy of the submergence deduction of a large range of regions.

[0082] In a possible implementation, the method further includes: determining a coupling region based on the two-dimensional grid and the three-dimensional grid of the sub-region, the coupling region representing a connection region of the ground surface and the sub-region; obtaining an intermediate three-dimensional grid representing the coupling region, a grid precision of the intermediate three-dimensional grid being between a grid precision of the two-dimensional grid and a grid precision of the three-dimensional grid of the sub-region; and determining the third parameter corresponding to each three-dimensional grid of the sub-region based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, including: determining a transition parameter corresponding to the intermediate three-dimensional grid based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid; determining the third parameter corresponding to at least part of the three-dimensional grid of the sub-region based on the transition parameter; and performing hydrodynamic simulation based on the third parameter corresponding to at least part of the three-dimensional grid of the sub-region and the three-dimensional hydrodynamic model of the sub-region to determine the third parameter corresponding to each three-dimensional grid of the sub-region.

[0083] The coupling region represents a connection region of the ground surface and the sub-region. The coupling region can include a part of the ground surface region. The two-dimensional grid is a planar grid. The three-dimensional grid of the sub-region is a three-dimensional grid, and the three-dimensional grid of the sub-region is distributed on multiple grid layers. The coupling region can include a region in the real world corresponding to the three-dimensional grid of the sub-region and the two-dimensional grid that are adjacent to each other.

[0084] A plurality of intermediate three-dimensional grids can be constructed for the coupling region. The plurality of intermediate three-dimensional grids can represent a three-dimensional structure of the coupling region. The grid precision of the intermediate three-dimensional grid can be between the grid precision of the two-dimensional grid and the grid precision of the three-dimensional grid of the sub-region.

[0085] In an example, in a direction from the two-dimensional grid to the three-dimensional grid of the sub-region, the grid precision of the intermediate three-dimensional grid can gradually change from the grid precision of the two-dimensional grid to the grid precision of the three-dimensional grid of the sub-region.

[0086] A single intermediate three-dimensional grid can correspond to at least one transition parameter. The transition parameter can represent a water flow state in the coupling region. The transition parameter can be a hydrodynamic parameter. A single node of the intermediate three-dimensional grid can correspond to at least one transition parameter to represent a water flow state of a position in the coupling region corresponding to the node. Alternatively, a center of the intermediate three-dimensional grid can correspond to at least one transition parameter to represent a water flow state of a position in the coupling region corresponding to the center.

[0087] In the embodiments of the present disclosure, the transition parameter corresponding to the intermediate grid can be calculated based on the second parameter.

[0088] For example, the transition parameter corresponding to each intermediate grid can be determined based on a distance between the intermediate three-dimensional grid and the two-dimensional grid and the second parameter corresponding to the two-dimensional grid intersecting with the intermediate grid.

[0089] Exemplarily, the second parameter can be mapped to the intermediate three-dimensional grid based on the second mapping relationship and a second mapping loss condition to obtain a transition parameter. The second mapping loss condition can represent the amplitude of increase or decrease of the second parameter with the change of the distance and height from the second grid.

[0090] The above is only an example, and the method of calculating the transition parameter corresponding to the intermediate grid based on the second parameter is not limited in the embodiments of the present disclosure. The process of determining the transition parameter will be exemplarily introduced below.

[0091] In the embodiments of the present disclosure, the intermediate three-dimensional grid and the sub-area three-dimensional grid located in the same layer can perform parameter equivalent transfer. In this way, at least part of the sub-area three-dimensional grid can obtain the corresponding third parameter. Here, the part of the sub-area three-dimensional grid can be the sub-area three-dimensional grid adjacent to the intermediate three-dimensional grid. In another example, the third parameter corresponding to a single sub-area three-dimensional grid can be determined based on the transition parameter corresponding to the intermediate three-dimensional grid around the single sub-area three-dimensional grid. The above is only an example, and the method of determining the third parameter corresponding to at least part of the sub-area three-dimensional grid based on the transition parameter is not limited in the embodiments of the present disclosure.

[0092] In the embodiments of the present disclosure, the sub-area hydrodynamic three-dimensional model can simulate the water flow state around and / or inside the sub-area based on the third parameter corresponding to at least part of the sub-area three-dimensional grid to obtain the third parameter corresponding to each of the sub-area three-dimensional grid.

[0093] A single sub-area three-dimensional grid can correspond to at least one third parameter. The third parameter can represent the water flow state around and / or inside the sub-area. The third parameter can be a hydrodynamic parameter. A single node of the sub-area three-dimensional grid can correspond to at least one third parameter to represent the water flow state of the position around or inside the sub-area corresponding to the node. Or the center of the sub-area three-dimensional grid can correspond to at least one third parameter to represent the water flow state of the position around or inside the sub-area corresponding to the center.

[0094] In the embodiments of the present disclosure, an intermediate three-dimensional grid is established. The intermediate three-dimensional grid and the sub-area three-dimensional grid are both three-dimensional grids. Moreover, the grid precision of the intermediate three-dimensional grid is between that of the two-dimensional grid and the sub-area three-dimensional grid. In this way, in the process of determining the three-dimensional parameter, the intermediate three-dimensional grid plays a transition and progressive role, reduces the error introduced in the process of determining the third parameter, and makes the connection between the hydrodynamic two-dimensional model and the hydrodynamic three-dimensional model smoother.

[0095] In a possible implementation, the determining the transition parameters corresponding to the intermediate three-dimensional grids based on the two-dimensional grid and the second parameters corresponding to the two-dimensional grid comprises: determining initial transition parameters corresponding to the intermediate three-dimensional grids intersecting with the two-dimensional grid based on the two-dimensional grid and the second parameters corresponding to the two-dimensional grid, and taking the intermediate three-dimensional grids intersecting with the two-dimensional grid as first intersecting three-dimensional grids; and recursively assigning values to each of the intermediate three-dimensional grids in a direction from near to far of the first intersecting three-dimensional grids based on the initial transition parameters to obtain the transition parameters.

[0096] The coupling region and the ground surface can be adjacent and partially overlapped. Thus, part of the two-dimensional grid can be adjacent to or intersect with part of the three-dimensional intermediate grid. For ease of description, the intermediate three-dimensional grids intersecting with the two-dimensional grid are named as first intersecting three-dimensional grids. The second parameters of the two-dimensional grid can be assigned to the first intersecting three-dimensional grids intersecting with the second grid as initial transition parameters corresponding to the first intersecting three-dimensional grids.

[0097] Then, the initial transition parameters can be used to recursively assign values to other intermediate three-dimensional grids. For example, the transition parameters corresponding to the intermediate three-dimensional grids adjacent to the first intersecting three-dimensional grids can be determined based on the first intersecting three-dimensional grids, the distances between the centers or the nodes of the first intersecting three-dimensional grids and the intermediate three-dimensional grids, and the initial transition parameters corresponding to the first intersecting three-dimensional grids. Next, the transition parameters corresponding to the other intermediate three-dimensional grids adjacent to the intermediate three-dimensional grids can be determined based on the intermediate three-dimensional grids, the distances between the centers or the nodes of the intermediate three-dimensional grids and the other intermediate three-dimensional grids, and the transition parameters corresponding to the intermediate three-dimensional grids. The same method can be used to determine the transition parameters corresponding to all the intermediate three-dimensional grids.

[0098] Using the method of the embodiments of the present disclosure, the second parameters can be converted from the two-dimensional plane to the three-dimensional grids in multiple levels recursively. This makes the transition parameters smoothly advance to the three-dimensional intermediate grids. This lays a foundation for subsequently determining the third parameters corresponding to the sub-area three-dimensional grids, and improves the stability of the hydrodynamic three-dimensional model in which the sub-area three-dimensional grids are located.

[0099] In a possible implementation, the determining the third parameters corresponding to at least part of the sub-area three-dimensional grids based on the transition parameters comprises: determining the third parameters corresponding to the sub-area three-dimensional grids intersecting with the intermediate three-dimensional grids based on the transition parameters; and the method further comprises: taking the sub-area three-dimensional grids intersecting with the intermediate three-dimensional grids as second intersecting three-dimensional grids; and taking the third parameters corresponding to the second intersecting three-dimensional grids as the initial hydraulic parameters of the sub-area hydrodynamic three-dimensional model.

[0100] The partial intermediate three-dimensional mesh can be adjacent to, overlap with, or intersect with the partial subregion three-dimensional mesh. For ease of description, the subregion three-dimensional mesh adjacent to, overlapping with, or intersecting with the partial intermediate three-dimensional mesh can be named as a second interface three-dimensional mesh.

[0101] For example, the average of the transition parameters corresponding to one or more intermediate three-dimensional meshes interfacing with the second interface three-dimensional mesh can be taken as the third parameter of the second interface three-dimensional mesh. For example only, the embodiments of the present disclosure do not limit this.

[0102] The intermediate three-dimensional mesh and the second interface three-dimensional mesh can each include a plurality of layers, and the number of layers can be the same or different. The second interface three-dimensional mesh can be a three-dimensional mesh representing the edge of the subregion. The third parameter corresponding to the second interface three-dimensional mesh is taken as the initial hydraulic parameter of the subregion hydrodynamic three-dimensional model, and the hydrodynamic three-dimensional model obtains the hydraulic parameter of the edge of the subregion. Thus, initial conditions are provided for the simulation of water flow around and inside the subregion.

[0103] In the embodiments of the present disclosure, the third parameter of the edge of the subregion is determined based on the transition parameter, and the second parameter is successfully passed to the partial subregion three-dimensional mesh, thereby improving the accuracy and stability of the simulation process of the regional hydrodynamic three-dimensional model.

[0104] In a possible implementation, the determining, based on the transition parameter, of the third parameter corresponding to the subregion three-dimensional mesh interfacing with the intermediate three-dimensional mesh includes: determining, for a single second interface three-dimensional mesh, at least one intermediate three-dimensional mesh closest to the single second interface three-dimensional mesh; and determining, based on the transition parameter corresponding to the at least one intermediate three-dimensional mesh, an initial parameter corresponding to the single second interface three-dimensional mesh.

[0105] In the embodiments of the present disclosure, the at least one intermediate three-dimensional mesh closest to the single second interface three-dimensional mesh can be at least one intermediate three-dimensional mesh within a preset radius range centered on the single second interface three-dimensional mesh. The transition parameter corresponding to the at least one intermediate three-dimensional mesh can be a transition parameter corresponding to a node of the at least one intermediate three-dimensional mesh falling within the preset radius range. The initial parameter corresponding to the single second interface three-dimensional mesh can be determined based on the transition parameter corresponding to the at least one intermediate three-dimensional mesh through nearest-neighbor interpolation.

[0106] In this way, the transition parameter most relevant to the second interface three-dimensional mesh can be selected, and the accuracy of determining the initial parameter is further improved.

[0107] In a possible implementation, the initial transition parameter includes a first flow rate, and the method of determining the initial transition parameter of the intermediate three-dimensional grid intersecting with the two-dimensional grid based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid includes: vertically assigning the first flow rate corresponding to the first intersecting three-dimensional grid to a plurality of intermediate three-dimensional grids under the first intersecting three-dimensional grid, to obtain a second flow rate corresponding to each of the plurality of intermediate three-dimensional grids under the first intersecting three-dimensional grid; and recursively assigning each of the intermediate three-dimensional grids in a direction from near to far of the first intersecting three-dimensional grid based on the initial transition parameter, to obtain each of the transition parameters, which includes a flow rate.

[0108] As described above, the second parameter can include a ground water flow rate corresponding to a plurality of positions on the ground. In a case where the second parameter is transmitted to the first intersecting three-dimensional grid, the first intersecting three-dimensional grid corresponds to the initial transition parameter. The initial transition parameter can include a first flow rate.

[0109] There can be a plurality of intermediate three-dimensional grids under the first intersecting three-dimensional grid. The second flow rate corresponding to each of the plurality of intermediate three-dimensional grids under the first intersecting three-dimensional grid can be equal to the first flow rate. This is because the accumulated water in the city belongs to shallow water, and there is no obvious vertical flow.

[0110] In the embodiment of the present disclosure, the first flow rate and the second flow rate can be transmitted to each of the intermediate three-dimensional grids by using an interpolation method and taking the distance between the intermediate three-dimensional grids as a weight, so that the intermediate three-dimensional grid can correspond to a flow rate. The interpolation method here is preferably bilinear interpolation, so as to improve the determination efficiency of the transition parameter.

[0111] In this way, the accuracy, smoothness and efficiency of determining the flow rate in the transition parameter can be improved.

[0112] In a possible implementation, the initial transition parameter further includes a first water level, and the transition parameter further includes a pressure and a water phase volume fraction. The method further includes: determining a liquid surface height based on the first water level, determining a water phase integral corresponding to the intermediate three-dimensional grid based on the liquid surface height and an elevation corresponding to the intermediate three-dimensional grid, and determining a pressure corresponding to the intermediate three-dimensional grid based on the liquid surface height, an initial pressure corresponding to the liquid surface height and the elevation corresponding to the intermediate three-dimensional grid.

[0113] As described above, the second parameter can include a ground water level corresponding to a plurality of positions on the ground. In a case where the second parameter is transmitted to the first intersecting three-dimensional grid, the first intersecting three-dimensional grid corresponds to the initial transition parameter. The initial transition parameter can include a first water level.

[0114] In the embodiments of the present disclosure, the liquid level can be determined based on the first water level or the ground water level in the second parameter. The height of the intermediate three-dimensional grid can be compared with the liquid level, and the water phase integral of the intermediate grid can be determined based on the comparison result.

[0115] When the height of the intermediate three-dimensional grid is greater than the liquid level, the water phase integral of the intermediate three-dimensional grid can be a first value. When the height of the intermediate three-dimensional grid is not greater than the liquid level, the water phase integral of the intermediate three-dimensional grid can be a second value. The first value is not equal to the second value. For example, the first value is 0, and the second value is 1.

[0116] In the embodiments of the present disclosure, the initial pressure corresponding to the liquid level can be set, and the initial pressure can be equal to the atmospheric pressure. For example, the initial pressure can be 0. The pressure of the intermediate three-dimensional grid can be determined based on the initial pressure and the height of the intermediate three-dimensional grid according to the hydrostatic equilibrium gradient.

[0117] In the embodiments of the present disclosure, more types of transition parameters including pressure and water phase integral can be further determined based on the initial transition parameter. In this way, the three-dimensional water simulation model can be provided with multiple types of accurate and smooth parameters, and the accuracy of the sub-region water flow simulation can be improved.

[0118] In a possible implementation, the method further includes: determining the river and pipe two-dimensional model based on the rainstorm disaster data, the target region terrain data, the river pipe network data, the land use data, the river flood data, and the river downstream water level data; determining the surface water power two-dimensional model based on the target region terrain data, the overflow water quantity data, the land use data, the soil data, and the parameter file; determining the sub-region water power three-dimensional model based on the geometry data and the fluid parameter data of the sub-region; and determining the sub-region waterlogging range and depth corresponding to the first time based on the third parameter, including: inputting the third parameter into the sub-region water power three-dimensional model to obtain the sub-region waterlogging range and depth.

[0119] The storm disaster data can include data representing rainfall intensity, such as a rainfall intensity time series. The urban terrain data can be data including surface elevation information, such as a digital elevation model (DEM) of the target region. The river channel network data can include information such as the locations, sizes, slopes, and the like of rivers, rainwater pipes, inspection wells, water outlets, and the like. Because the probability and speed of water infiltration are different due to different land use types, land use data also needs to be provided. The land use data can include the permeability, roughness, and surface coverage type of different surface types. The river flood data can include upstream flow of the river. The downstream water level data of the river can include astronomical tide data, storm surge data, sea level rise data, and the like. The overflow water quantity data can be the location of overflow of the river or pipe, and the time series of the overflow water quantity. The soil data can be parameters representing soil type, permeability, and water holding capacity. The model parameter file can include parameters defining the model operation.

[0120] The geometry data of the sub-region can represent the structure of the sub-region. It can include point cloud data, high-precision digital surface model (DSM), or CAD file data corresponding to the sub-region. The geometry data can be data that removes redundant structures, repairs grid holes, and normalizes model scales. The three-dimensional water simulation model can generate an initial Cartesian grid representing the structure of the sub-region based on the collected data. For complex features, grid refinement is performed on the initial Cartesian grid to improve the accuracy level. High-precision geometry data is obtained. The fluid parameter data can be the density and viscosity of water in the flow, and can also include turbulence model parameters such as dissipation rate and turbulent kinetic energy.

[0121] The third parameter at a single time is input into the sub-region water simulation three-dimensional model as input data, and the sub-region waterlogging range and depth corresponding to the time can be obtained. If the third parameters corresponding to multiple times are input into the sub-region water simulation three-dimensional model, the sub-region surrounding or internal flooding inference process can be obtained. In one example, the sub-region three-dimensional grid can also be divided to obtain multiple partitions, and the number of partitions can be determined based on the number of CPU cores. In this way, multiple partitions can be simulated in parallel to improve simulation efficiency.

[0122] Using the method in the embodiments of the present disclosure, a two-dimensional model of the river and pipe, a two-dimensional model of surface water simulation, and a three-dimensional model of sub-region water simulation corresponding to the city in the scenario of disastrous rainstorm, flood, and other water disasters can be constructed. Not only can the surface, river, and pipe be efficiently overflowed and flooded, but the water flow around or inside the sub-region can also be accurately simulated. Both the accuracy and efficiency of the entire target region flooding process are taken into account.

[0123] Figure 2A structural schematic diagram of a regional inundation deduction device is provided for an embodiment of the present disclosure. The device 20 comprises:

[0124] An overflow path determination unit 21 is configured to determine an overflow path corresponding to a first time based on a first parameter of water flow at a plurality of positions in a pipeline and / or a river channel in a target region at the first time according to a two-dimensional model of the pipeline and the river channel.

[0125] A ground surface water accumulation range and depth determination unit 22 is configured to determine a ground surface water accumulation range and depth corresponding to the first time based on a second parameter of water flow at a plurality of positions on a ground surface in the target region at the first time according to a two-dimensional hydrodynamic model of the ground surface, the two-dimensional hydrodynamic model of the ground surface comprising a two-dimensional grid representing the ground surface, each two-dimensional grid corresponding to at least one second parameter.

[0126] A third parameter determination unit 23 is configured to determine a third parameter corresponding to a three-dimensional grid of each sub-region in a sub-region hydrodynamic three-dimensional model based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid.

[0127] A sub-region water accumulation range and depth determination unit 24 is configured to determine a sub-region water accumulation range and depth corresponding to the first time based on the third parameter.

[0128] A target region inundation deduction unit 25 is configured to deduce an inundation process of the target region based on a plurality of overflow paths corresponding to the first time respectively, the ground surface water accumulation range and depth, and the sub-region water accumulation range and depth.

[0129] In a possible implementation, the device 20 further comprises:

[0130] A coupling region determination unit is configured to determine a coupling region based on the two-dimensional grid and the three-dimensional grid of the sub-region, the coupling region representing a connecting region of the ground surface and the sub-region.

[0131] An intermediate three-dimensional grid representing the coupling region is obtained, a grid precision of the intermediate three-dimensional grid being between a grid precision of the two-dimensional grid and a grid precision of the three-dimensional grid of the sub-region.

[0132] The third parameter determination unit 23 is further configured to:

[0133] determine a transition parameter corresponding to the intermediate three-dimensional grid based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid;

[0134] determine the third parameter corresponding to at least part of the three-dimensional grid of the sub-region based on the transition parameter;

[0135] The third parameter corresponding to each of the sub-area three-dimensional grids is determined based on the third parameter corresponding to the at least part of the sub-area three-dimensional grids and the sub-area water power three-dimensional model.

[0136] In a possible implementation, the third parameter determination unit 23 is further configured to:

[0137] Based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, an initial transition parameter corresponding to an intermediate three-dimensional grid intersecting with the two-dimensional grid is determined, and the intermediate three-dimensional grid intersecting with the two-dimensional grid is taken as a first intersecting three-dimensional grid.

[0138] Based on the initial transition parameter, each of the intermediate three-dimensional grids is recursively assigned a value in a direction from near to far of the first intersecting three-dimensional grid, to obtain the transition parameter.

[0139] In a possible implementation, the third parameter determination unit 23 is further configured to:

[0140] Based on the transition parameter, a third parameter corresponding to a sub-area three-dimensional grid intersecting with the intermediate three-dimensional grid is determined.

[0141] The apparatus 20 further includes:

[0142] A second intersecting three-dimensional grid determination unit is configured to take a sub-area three-dimensional grid intersecting with the intermediate three-dimensional grid as a second intersecting three-dimensional grid.

[0143] A hydraulics initial parameter determination unit is configured to take the third parameter corresponding to the second intersecting three-dimensional grid as a hydraulics initial parameter of the sub-area water power three-dimensional model.

[0144] In a possible implementation, the third parameter determination unit 23 is further configured to:

[0145] For a single second intersecting three-dimensional grid, at least one intermediate three-dimensional grid closest to the single second intersecting three-dimensional grid is determined.

[0146] Based on the transition parameter corresponding to the at least one intermediate three-dimensional grid, an initial parameter corresponding to the single second intersecting three-dimensional grid is determined.

[0147] In a possible implementation, the initial transition parameter includes a first flow rate, and the third parameter determination unit 23 is further configured to:

[0148] A first flow rate corresponding to the first intersecting three-dimensional grid is vertically assigned to a plurality of intermediate three-dimensional grids under the first intersecting three-dimensional grid, to obtain a second flow rate corresponding to each of the plurality of intermediate three-dimensional grids under the first intersecting three-dimensional grid.

[0149] interpolating each of the first flow rate and the second flow rate to determine the transition parameter, the transition parameter comprising a flow rate.

[0150] In a possible implementation, the initial transition parameter further comprises a first water level, and the transition parameter further comprises a pressure and a water phase volume fraction, and the third parameter determination unit 23 is further configured to:

[0151] determine a liquid surface height based on the first water level, and determine a water phase integral corresponding to the intermediate three-dimensional grid based on the liquid surface height and an elevation corresponding to the intermediate three-dimensional grid.

[0152] determine a pressure corresponding to the intermediate three-dimensional grid based on the liquid surface height, an initial pressure corresponding to the liquid surface height, and the elevation corresponding to the intermediate three-dimensional grid.

[0153] In some embodiments, the apparatus provided by the embodiments of the present disclosure has functions or includes modules for performing the methods described in the above method embodiments, and the specific implementation can refer to the description of the above method embodiments. For brevity, it will not be described here.

[0154] The embodiments of the present disclosure also provide an electronic device, including a memory, a processor and a computer program stored in the memory, and the processor executes the computer program to implement the steps of the above method.

[0155] The embodiments of the present disclosure also provide a non-volatile computer readable storage medium having a computer program stored thereon, and the computer program is executed by a processor to implement the steps of the above method.

[0156] The embodiments of the present disclosure also provide a computer program product, including a computer program or a non-volatile computer readable storage medium carrying a computer program, and the computer program is executed by a processor to implement the steps of the above method.

[0157] Figure 3 The structure schematic diagram of the electronic device for area flooding deduction provided by the embodiments of the present disclosure is provided. For example, the electronic device 1900 can be provided as a server or a terminal device. Referring to Figure 3 , the electronic device 1900 includes a processing component 1922, which further includes one or more processors, and a memory resource represented by a memory 1932, for storing instructions executable by the processing component 1922, such as an application program. The application program stored in the memory 1932 can include one or more than one module each corresponding to a set of instructions. In addition, the processing component 1922 is configured to execute the instructions to perform the above method.

[0158] The electronic device 1900 can further include a power supply component 1926 configured to perform power management of the electronic device 1900, a wired or wireless network interface 1950 configured to connect the electronic device 1900 to a network, and an input / output interface 1958 (I / O interface). The electronic device 1900 can operate based on an operating system stored in the memory 1932, such as Windows Server TM , Mac OS X TM , Unix TM , Linux TM , FreeBSD TM , or the like.

[0159] In exemplary embodiments, there is also provided a non-transitory computer readable storage medium, such as the memory 1932 including computer program instructions, which can be executed by the processing component 1922 of the electronic device 1900 to perform the above method.

[0160] The computer readable storage medium can be a tangible device that can retain and store instructions for execution by a processor. The computer readable storage medium can be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. More specific examples (a non-exhaustive list) of the computer readable storage medium include the following: a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanically encoded device such as punch-cards or punched tape, and any suitable combination of the foregoing. A computer readable storage medium, as used herein, is not to be construed as being transitory signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through a waveguide or other transmission media (e.g., light pulses passing through a fiber-optic cable), or electrical signals transmitted through a wire.

[0161] The computer program (or computer readable program instructions) described herein can be downloaded from a computer readable storage medium to various computing / processing devices by way of a network, such as the Internet, a local area network, a wide area network, and / or a wireless network. The network can include copper transmission cables, optical transmission fibers, wireless transmission, routers, firewalls, switches, gateway computers and / or edge servers. A network adapter card or network interface in each computing / processing device receives computer readable program instructions from the network and forwards the computer readable program instructions for storage in a computer readable storage medium within the respective computing / processing device.

[0162] Computer readable program instructions for carrying out operations of the present disclosure can be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, or either source code or object code written in any combination of one or more programming languages, including an object oriented programming language such as Smalltalk, C++ or the like, and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The computer readable program instructions can execute entirely on the user's computing / processing device, partly on the user's computing / processing device, as a stand-alone software package, partly on the user's computing / processing device and partly on a remote computing / processing device or entirely on the remote computing / processing device or server. In the latter scenario, the remote computing / processing device can be connected to the user's computing / processing device through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection can be made to an external computing / processing device, for example, through the Internet using an Internet Service Provider. In some embodiments, electronic circuitry including, for example, programmable logic circuitry, field-programmable gate array (FPGA), or programmable logic array (PLA) can execute the computer readable program instructions by utilizing state information of the computer readable program instructions to personalize the electronic circuitry, in order to perform aspects of the present disclosure.

[0163] The computer readable program instructions can also be loaded onto a computing / processing device, other programmable data processing apparatus, or other device to cause a series of operations to be performed on the computing / processing device, other programmable apparatus or other device to produce a computer implemented process, such that the instructions which execute on the computing / processing device, other programmable apparatus, or other device implement the operations specified in the flow diagrams and / or block diagrams.

[0164] These computer readable program instructions can be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions / acts specified in the flowchart and / or block diagram block or blocks. These computer readable program instructions can also be stored in a computer readable storage medium that can include a non-transitory computer readable storage medium that can be a computer- readable storage medium having no data storage cycles that change state. The instructions can be executed by one or more processors of a computer, to cause a series of operational elements or steps to be performed on the computer to produce a computer implemented process. Such instructions can also be stored and / or executed by other computer-readable media. Computer-readable media storing the computer readable instructions can include computers, processors, or other programmable data processing apparatuses capable of receiving, storing, and / or executing instructions.

[0165] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational elements or steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer, other programmable elements, or other

[0166] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational elements or steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer, other programmable elements, or other

[0167] Embodiments of the present disclosure have been described above, and the description is intended to be illustrative of the embodiments and not restrictive of the disclosure. Many modifications and variations of the described embodiments are possible in light of this disclosure without departing from the scope and spirit of the described embodiments. The choice of words in this document is intended to best explain the principles of the embodiments, the practical application, or technical improvement over prior art, or to enable others skilled in the art to understand the embodiments disclosed herein.

Claims

1. A method of area flooding deduction, characterized by, The method comprises: determining, based on a two-dimensional model of a river channel and a pipeline, an overflow path corresponding to a first time according to a first parameter of water flow at a plurality of positions in the pipeline and / or the river channel in a target region corresponding to the first time; determining, based on a two-dimensional surface water dynamic model, a surface water accumulation range and depth corresponding to the first time according to a second parameter of water flow at a plurality of positions on the surface in the target region corresponding to the first time, the two-dimensional surface water dynamic model comprising a two-dimensional grid representing the surface, each two-dimensional grid corresponding to at least one second parameter; determining, based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, a third parameter corresponding to each three-dimensional grid in a sub-region water dynamic three-dimensional model; determining, based on the third parameter, a sub-region water accumulation range and depth corresponding to the first time; deducing a flooding process of the target region based on a plurality of overflow paths, surface water accumulation ranges and depths, and sub-region water accumulation ranges and depths corresponding to the first time respectively.

2. The method of claim 1, wherein, The method further comprises: determining a coupling region based on the two-dimensional grid and the three-dimensional grid of the sub-region, the coupling region representing a connection region between the surface and the sub-region; obtaining an intermediate three-dimensional grid representing the coupling region, the grid precision of the intermediate three-dimensional grid being between the grid precision of the two-dimensional grid and the grid precision of the three-dimensional grid of the sub-region; the determining, based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, a third parameter corresponding to each three-dimensional grid in a sub-region water dynamic three-dimensional model, comprises: determining, based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, a transition parameter corresponding to the intermediate three-dimensional grid; determining, based on the transition parameter, a third parameter corresponding to at least part of the three-dimensional grid of the sub-region; performing a hydrodynamic simulation based on the third parameter corresponding to at least part of the three-dimensional grid of the sub-region and the sub-region water dynamic three-dimensional model to determine the third parameter corresponding to each three-dimensional grid of the sub-region.

3. The method of claim 2, wherein, the determining, based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, a transition parameter corresponding to the intermediate three-dimensional grid, comprises: determining, based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid, an initial transition parameter corresponding to the intermediate three-dimensional grid intersecting with the two-dimensional grid, and taking the intermediate three-dimensional grid intersecting with the two-dimensional grid as a first intersection three-dimensional grid; based on the initial transition parameter, recursively assigning a value to each intermediate three-dimensional grid in a direction from near to far of the first intersection three-dimensional grid to obtain each transition parameter.

4. The method of claim 2, wherein, the determining, based on the transition parameter, a third parameter corresponding to at least part of the three-dimensional grid of the sub-region, comprises: determining, based on the transition parameter, a third parameter corresponding to the three-dimensional grid of the sub-region intersecting with the intermediate three-dimensional grid; the method further comprises: taking the three-dimensional grid of the sub-region intersecting with the intermediate three-dimensional grid as a second intersection three-dimensional grid; taking the third parameter corresponding to the second intersection three-dimensional grid as a hydrodynamics initial parameter of the sub-region water dynamic three-dimensional model.

5. The method of claim 4, wherein, The method further comprises: For a single second interface three-dimensional grid, determining at least one intermediate three-dimensional grid closest to the single second interface three-dimensional grid; Based on the transition parameters corresponding to the at least one intermediate three-dimensional grid, determining the initial parameters corresponding to the single second interface three-dimensional grid.

6. The method of claim 3, wherein, The initial transition parameters include a first flow rate, and the method further comprises: The first flow rate corresponding to the first interface three-dimensional grid is vertically assigned to a plurality of intermediate three-dimensional grids under the first interface three-dimensional grid, to obtain a second flow rate corresponding to each of the plurality of intermediate three-dimensional grids under the first interface three-dimensional grid; The method further comprises: The initial transition parameters include a first water level, and the transition parameters further include a pressure and a water phase volume fraction, and the method further comprises:

7. The method of claim 6, wherein, Based on the first water level, determining a liquid surface height, and based on the liquid surface height and an elevation corresponding to the intermediate three-dimensional grid, determining a water phase integral corresponding to the intermediate three-dimensional grid; Based on the liquid surface height, an initial pressure corresponding to the liquid surface height, and the elevation corresponding to the intermediate three-dimensional grid, determining a pressure corresponding to the intermediate three-dimensional grid. The method further comprises:

8. An area flooding deduction apparatus, characterized by comprising: An overflow path determination unit configured to determine an overflow path corresponding to a first time based on a first parameter of water flow at a plurality of positions in a pipe and / or a river in a target region at the first time according to a two-dimensional model of the river and the pipe; A surface water accumulation range and depth determination unit configured to determine a surface water accumulation range and depth corresponding to the first time based on a second parameter of water flow at a plurality of positions on the surface in the target region at the first time according to a two-dimensional surface water dynamic model, the two-dimensional surface water dynamic model including a two-dimensional grid representing the surface, each of the two-dimensional grids corresponding to at least one second parameter; A third parameter determination unit configured to determine a third parameter corresponding to each sub-region three-dimensional grid in a sub-region water dynamic three-dimensional model based on the two-dimensional grid and the second parameter corresponding to the two-dimensional grid; A sub-region water accumulation range and depth determination unit configured to determine a sub-region water accumulation range and depth corresponding to the first time based on the third parameter; A target region inundation deduction unit configured to deduce an inundation process of the target region based on a plurality of overflow paths, surface water accumulation ranges and depths, and sub-region water accumulation ranges and depths corresponding to the first time respectively. The processor executes the computer program to implement the steps of the method of any one of claims 1 to 7.

9. An electronic device comprising a memory, a processor, and a computer program stored on the memory, wherein the computer program, when executed by the processor, is arranged to perform the method of any one of claims 1 to 8. The computer program is executed by the processor to implement the steps of the method of any one of claims 1 to 7.

10. A non-transitory computer readable storage medium having stored thereon a computer program, characterized in that, ​

Images