A hydrological simulation calculation method, device and computer equipment

Abstract

The present invention discloses a hydrological simulation calculation method, device and computer equipment, the method comprising: obtaining watershed data of a target area; dividing the watershed data to obtain a plurality of sub-watershed data; creating at least one independent hydrological unit according to the plurality of sub-watershed data, wherein each independent hydrological unit includes at least one sub-watershed data among the plurality of sub-watershed data, and there is no hydraulic connection between any two independent hydrological units at the current moment; taking the current moment as the starting moment, after each preset time, according to the sub-watershed data included in each independent hydrological unit and pre-configured rainfall runoff conditions, obtaining the runoff data corresponding to each independent hydrological unit; if after the first preset time, the runoff corresponding to each independent hydrological unit has no hydraulic connection with other independent hydrological units, integrating the runoff data corresponding to each independent hydrological unit to obtain the hydrological simulation result of the target area.

Description

Technical Field

[0001] The present invention relates to the technical field of watershed hydrological simulation, and in particular to a hydrological simulation calculation method, device and computer equipment. Background Art

[0002] Hydrological models are important technical means to simulate watershed hydrology and reveal the hydrological response process of watersheds. Among them, distributed hydrological models can fully reflect the impact of spatial changes in rainfall and underlying surface elements on runoff in the watershed. In the development of hydrological simulation technology, distributed hydrological models have been widely used and developed with the above-mentioned simulation advantages.

[0003] However, with the development of GIS (Geographic Information System or Geo-Information system, GIS for short) and RS (Remote Sensing, RS for short) technology and its application in hydrological models, the resolution of spatial data such as DEM (Digital Elevation Model, DEM for short), soil, vegetation, and meteorology used in distributed hydrological model calculations has become increasingly fine. In order to optimize model parameters, find model parameters that can accurately reflect the characteristics of the watershed, and obtain accurate simulation results, the amount of simulation calculations has increased dramatically. At the same time, in order to meet the needs of hydrological simulation under a changing environment, the multi-scenario simulation based on distributed hydrological simulation has greatly increased the amount of calculations for distributed hydrological simulation. In addition, with the diversification of simulation needs, the emergence of two-dimensional simulation calculations has further increased the amount of simulation calculations. With the increasing amount of simulation calculations, the high requirements of real-time decision-making systems for simulation timeliness have made the improvement of simulation efficiency a difficult problem that needs to be solved urgently.

[0004] Faced with the huge amount of simulation calculations and the urgent need to improve the simulation efficiency, traditional serial calculations can no longer meet the computing power requirements of hydrological models. Parallel computing has begun to be applied in hydrological models and has gradually developed. At present, in the field of distributed hydrological simulation, the more mature research on the parallel characteristics of the basin itself is to discretize the basin river network space with the help of river network binary tree encoding for parallel calculation. When this method is used for calculation, the downstream sub-basin cannot start the calculation before the upstream sub-basin is calculated. As a result, a large number of computing units are in a suspended calculation state, and the parallelization degree of model calculation will drop rapidly in a short period of time. Based on the river network binary tree encoding method, the parallelism of basin calculation is limited to the fixed river network topology network, and the hydraulic connection between sub-basins that changes with the dynamic evolution of runoff generation and confluence process is not taken into account. Therefore, the parallelism potential of the basin's own runoff generation and confluence calculation is not fully explored. Summary of the invention

[0005] Therefore, the technical problem to be solved by the present invention is to overcome the defect of parallel computing limit efficiency of existing distributed hydrological simulation calculations due to their reliance on a fixed river network topology structure, thereby providing a distributed hydrological simulation calculation method, device and computer equipment based on dynamic independent hydrological unit division.

[0006] According to the first aspect, an embodiment of the present invention discloses a hydrological simulation calculation method, including: acquiring watershed data of a target area; dividing the watershed data to obtain multiple sub-watershed data; creating at least one independent hydrological unit based on the multiple sub-watershed data, wherein each independent hydrological unit includes at least one sub-watershed data among the multiple sub-watershed data, and there is no hydraulic connection between every two independent hydrological units at the current moment; taking the current moment as the starting moment, every preset time, according to the sub-watershed data included in each independent hydrological unit, and pre-configured rainfall and runoff conditions, the runoff data corresponding to each independent hydrological unit is obtained; if after the first preset time, the runoff corresponding to each independent hydrological unit has no hydraulic connection with other independent hydrological units, then the runoff data corresponding to each independent hydrological unit are integrated to obtain the hydrological simulation result of the target area.

[0007] Optionally, when there are multiple independent hydrological units, the method also includes: if after a first preset time, there is at least one independent hydrological unit whose corresponding runoff has a hydraulic connection with other independent hydrological units, then after the first preset time, the independent hydrological units among the multiple independent hydrological units that have a hydraulic connection with each independent hydrological unit in at least one independent hydrological unit are integrated to form a new independent hydrological unit, and the first preset time is any preset time after the start time; within a second preset time, according to the sub-basin data included in the new independent hydrological unit and the pre-configured rainfall runoff conditions, the corresponding runoff data are obtained, and the second preset time is the next time period of the first preset time; if the corresponding runoff of the new independent hydrological unit has no hydraulic connection with other independent hydrological units, then the runoff data corresponding to each new independent hydrological unit are integrated with the runoff data of all independent hydrological units among the multiple independent hydrological units that have no hydraulic connection with other independent hydrological units to obtain the hydrological simulation results of the target area.

[0008] Optionally, multiple sub-basin data are obtained based on the watershed data, specifically including: obtaining the river network water system in the target area based on the watershed data; classifying the river network water system to obtain the classified river network water system; dividing the sub-basin according to the classified river network water system to obtain multiple sub-basin data.

[0009] Optionally, runoff generation includes area runoff generation, the preconfigured rainfall runoff conditions include rainfall, evaporation and infiltration parameters, and the calculation process of area runoff generation includes: obtaining the infiltration amount of the first independent hydrological unit based on the infiltration parameters of the first independent hydrological unit; obtaining the area runoff generation based on rainfall, evaporation and infiltration.

[0010] Optionally, the calculation of runoff also includes the calculation of slope runoff data and the calculation of river runoff data. The preconfigured rainfall runoff conditions also include water flow parameters. The calculation process of slope runoff and river runoff includes: obtaining the flow of the first independent hydrological unit based on the preset flow area and water flow parameters of the first independent hydrological unit, wherein the first independent hydrological unit is any one of at least one independent hydrological unit; obtaining slope runoff data based on the flow and preset slope dynamic parameters; obtaining river runoff data based on the flow and preset river dynamic parameters.

[0011] Optionally, before calculating the runoff of each independent hydrological unit, the method also includes: allocating a calculation kernel according to the area of ​​a first independent hydrological unit, and obtaining data of a first calculation kernel used to calculate the runoff of the first independent hydrological unit, wherein the first independent hydrological unit is any one of at least one independent hydrological unit.

[0012] According to the second aspect, an embodiment of the present invention further discloses a hydrological simulation calculation device, comprising: an acquisition module for acquiring watershed data of a target area; a division module for obtaining a plurality of sub-watershed data by division based on the watershed data; a module for creating an independent hydrological unit for creating at least one independent hydrological unit based on the plurality of sub-watershed data, wherein each independent hydrological unit includes at least one sub-watershed data among the plurality of sub-watershed data, and there is no hydraulic connection between every two independent hydrological units at the current moment; a first runoff calculation module for taking the current moment as the starting moment, after a preset time, obtaining the runoff data corresponding to each independent hydrological unit according to the sub-watershed data included in each independent hydrological unit and the preconfigured rainfall runoff conditions; a first verification module for verifying whether the runoff corresponding to each independent hydrological unit has a hydraulic connection with other independent hydrological units after a first preset time, and if the runoff corresponding to each independent hydrological unit has no hydraulic connection with other independent hydrological units, the runoff data corresponding to each independent hydrological unit are integrated to obtain the hydrological simulation result of the target area.

[0013] Optionally, when there are multiple independent hydrological units, the device further includes: a dynamic integration module, which is used to, if after a first preset time, verify that the runoff and confluence corresponding to at least one independent hydrological unit has a hydraulic connection with other independent hydrological units, then after the first preset time, the independent hydrological units in the multiple independent hydrological units that have a hydraulic connection with each independent hydrological unit in at least one independent hydrological unit are formed into a new independent hydrological unit, and the first preset time is any preset time after the start time; a second runoff and confluence calculation module, which is used to calculate the runoff and confluence calculation according to the sub-units included in the new independent hydrological unit within the second preset time. The watershed data and the preconfigured rainfall-runoff conditions are used to obtain the corresponding runoff data. The second preset time is the next time period of the first preset time. The second verification module is used to verify whether the runoff corresponding to the new independent hydrological unit has a hydraulic connection with other independent hydrological units. If the runoff corresponding to the new independent hydrological unit has no hydraulic connection with other independent hydrological units, the runoff data corresponding to each new independent hydrological unit are integrated with the runoff data of all independent hydrological units in the multiple independent hydrological units that have no hydraulic connection with other independent hydrological units to obtain the hydrological simulation results of the target area.

[0014] According to the third aspect, an embodiment of the present invention also discloses a computer device, comprising: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor performs the steps of the hydrological simulation calculation method as in the first aspect or any optional embodiment of the first aspect.

[0015] According to the fourth aspect, an embodiment of the present invention further discloses a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the hydrological simulation calculation method of the first aspect or any optional embodiment of the first aspect.

[0016] The technical solution of the present invention has the following advantages:

[0017] The hydrological simulation calculation method provided by the present invention includes: obtaining the watershed data of the target area, and obtaining a plurality of sub-watershed data by dividing the watershed data. According to the plurality of sub-watershed data, the plurality of sub-watershed data are divided into at least one independent water unit, and there is no hydraulic connection between each independent hydrological unit. Taking the current moment as the starting moment, after each preset time, the runoff corresponding to each independent hydrological unit is obtained according to the sub-watershed data included in each independent hydrological unit and the pre-configured rainfall runoff conditions. If there is no hydraulic connection between the runoff of each independent hydrological unit and other independent hydrological units, the runoff data of each independent hydrological unit are integrated to output the hydrological simulation result corresponding to the target area.

[0018] In this way, the watershed data of the target area is extracted and divided to obtain multiple sub-watershed data of the target area. Furthermore, at least one independent hydrological unit is created based on the hydraulic connection between the sub-watershed data, so that there is no hydraulic connection between each independent hydrological unit. Since there is no hydraulic connection between each independent hydrological unit, when calculating the runoff corresponding to each independent hydrological unit, there is no need to wait for the calculation of the upstream node to be completed before calculating the runoff of the downstream node, thereby solving the problem of the efficiency limit of parallel computing that exists in the traditional distributed hydrological simulation calculation that relies on a fixed river network topology network. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0020] Figure 1 A flowchart of a specific example of a hydrological simulation calculation method in an embodiment of the present invention;

[0021] Figure 2 A flowchart of a specific example of a hydrological simulation calculation method in an embodiment of the present invention;

[0022] Figure 3 A schematic diagram of a specific example of a hydrological simulation calculation method in an embodiment of the present invention;

[0023] Figure 4 A flowchart of a specific example of a hydrological simulation calculation method in an embodiment of the present invention;

[0024] Figure 5 A flowchart of a specific example of a hydrological simulation calculation method in an embodiment of the present invention;

[0025] Figure 6 A flowchart of a specific example of a hydrological simulation calculation method in an embodiment of the present invention;

[0026] Figure 7 A flowchart of a specific example of a hydrological simulation calculation method in an embodiment of the present invention;

[0027] Figure 8 A flowchart of a specific example of a hydrological simulation calculation method in an embodiment of the present invention;

[0028] Fig. 9 is a principle block diagram of a specific example of a hydrological simulation computing device in an embodiment of the present invention;

[0029] Fig.10 FIG. 4 is a specific example diagram of a computer device in an embodiment of the present invention. DETAILED DESCRIPTION

[0030] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0031] In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc., indicating the orientation or positional relationship, are based on the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as limiting the present invention. In addition, the terms "first", "second", and "third" are used for descriptive purposes only, and cannot be understood as indicating or implying relative importance.

[0032] In the description of the present invention, it should be noted that, unless otherwise clearly specified and limited, the terms "installed", "connected", and "connected" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be a direct connection, or it can be indirectly connected through an intermediate medium, it can also be the internal connection of two components, it can be a wireless connection, or it can be a wired connection. For ordinary technicians in this field, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0033] In addition, the technical features involved in the different embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0034] In view of the technical problems mentioned in the background technology, the embodiment of the present invention provides a hydrological simulation calculation method, see Figure 1 As shown, Figure 1 A flow chart of a hydrological simulation calculation method provided by an embodiment of the present invention, the method comprising the following steps:

[0035] Step 101, obtaining watershed data of the target area.

[0036] For example, the target area may be a specific geographical area where hydrological simulation is required, wherein the watershed data of the target area may be digital elevation model (DEM) data, soil type data, land use data, and remote sensing image data, etc. The embodiment of the present invention does not limit the scope of the target area and the type and quantity of watershed data, and those skilled in the art may determine them according to actual needs.

[0037] Step 102: divide the watershed data into multiple sub-watershed data.

[0038] For example, after obtaining the watershed data, due to the spatial differences in different geographical locations, in order to conduct a more accurate analysis of the target area, it is necessary to perform preprocessing operations such as filling depressions, flow direction analysis, flow calculation, and river linking on the DEM data in the watershed data before dividing it into sub-watersheds.

[0039] Optionally, in a specific embodiment, Figure 2 As shown in the figure, the division based on the watershed data to obtain multiple sub-watershed data can be achieved in the following ways:

[0040] Step 1021, analyzing the watershed data to obtain the river network in the target area.

[0041] Step 1022, classify the river network water system to obtain a classified river network water system.

[0042] Step 1023, divide the sub-basin according to the classified river network system to obtain multiple sub-basin data.

[0043] For example, after the DEM data is preprocessed, the river network is extracted based on the preprocessed DEM data, and the extracted river network is classified. The division of sub-basins is based on the classified river network.

[0044] For example, Figure 3As shown in the figure, based on the DEM data of the target area, the target area is divided into 7 sub-basins, and the sub-basin classification adopts the same method as the water system classification. After the sub-basins are divided, the characteristic parameters of the 7 sub-basins can be extracted respectively with the divided sub-basins as units. The characteristic parameters of each sub-basin required for hydrological model simulation are extracted through the basin data (DEM, soil type data, land use data and remote sensing image data, etc.), mainly including area, slope, permeability coefficient, impermeability coefficient and river channel roughness, etc. The extracted parameter characteristics and the corresponding sub-basins can be stored in the database in a mapping relationship for subsequent runoff calculation.

[0045] Step 103: create at least one independent hydrological unit based on the multiple sub-watershed data, wherein each independent hydrological unit includes at least one sub-watershed data among the multiple sub-watershed data, and there is no hydraulic connection between any two independent hydrological units at the current moment.

[0046] For example, after the sub-basins are divided, all sub-basins are divided into at least one independent hydrological unit according to whether there is hydraulic connection between the sub-basins. Hydraulic connection refers to the hydraulic exchange or mutual influence of hydraulic states between hydrological units.

[0047] For example, when dividing independent hydrological units, the division is carried out according to the key node locations in the basin, the distribution of observation sites and the specific simulation requirements. Specifically, for urban basins, the distribution of key nodes (such as outlets, storage tanks, etc.) can be considered on the basis of the known distribution of water systems and the layout of urban rainwater pipe networks, and the independent hydrological units can be divided so that each sub-basin that converges into the key node constitutes an independent hydrological unit. When conducting two-dimensional flood simulation, the sub-basins contained in the frequently flooded area or the waterlogged area can be integrated into independent hydrological units with reference to historical flood data.

[0048] like Figure 3 As shown in the figure, when the division is based on the observation station, when the discharge outlet of sub-basin 2 is equipped with an observation station but there is no observation station in the upper-level sub-basin 1, sub-basin 1 and sub-basin 2 are integrated into independent hydrological unit 1; when the division is based on the simulation requirements, when the runoff area to be calculated includes sub-basins 3 and 4, sub-basins 3 and 4 are integrated into independent hydrological unit 2; when the division is based on the presence or absence of hydraulic connection, when sub-basin 5 has no intersection with any tributary of other sub-basins, it means that sub-basin 5 has no hydraulic connection with other sub-basins, and it is set as independent hydrological unit 3. At the same time, sub-basins 6 and 7 in the floodplain are integrated into independent hydrological unit 4. At this time, a total of 4 independent hydrological units are obtained to complete the subsequent runoff calculation.

[0049] Step 104, starting from the current time, after a preset time, the runoff data corresponding to each independent hydrological unit is obtained according to the sub-basin data included in each independent hydrological unit and the pre-configured rainfall runoff conditions.

[0050] Exemplarily, after obtaining each independent hydrological unit, after each preset time has passed, the runoff data within the corresponding preset time period is calculated based on the sub-basin data in the corresponding independent hydrological unit and the preconfigured rainfall-runoff conditions, wherein the preconfigured rainfall-runoff conditions may be data such as rainfall, evaporation, and infiltration parameters within the preset time, and the preset time may be 1 hour. The runoff calculation is based on the rainfall-runoff conditions within 1 hour.

[0051] The calculation of runoff generation and confluence includes the calculation of runoff generation, slope confluence and river confluence of the corresponding area in the independent hydrological unit. The calculation of runoff generation in the area includes evaporation and infiltration calculation by Horton formula for each hydrological unit. The runoff confluence calculation formed by deducting evaporation and infiltration from rainfall uses Manning formula and dynamic wave method to calculate the slope confluence and river confluence process in the independent hydrological unit. After the calculation is completed, the corresponding runoff generation and confluence data can be stored in the database for subsequent verification of the hydraulic connection between independent hydrological units.

[0052] Optionally, after the division of independent hydrological units is completed, computing kernels can be allocated according to the area size of each independent hydrological unit, where the computing kernel is the specific operator size used to calculate each independent hydrological unit, thereby achieving reasonable allocation of computing kernels and avoiding insufficient computer computing power in the process of calculating runoff.

[0053] For example, when the areas corresponding to independent hydrological units 1, 2, and 3 are 50 km 2 、30km 2 and 20km 2 When , the corresponding operator sizes account for 50%, 30% and 20% of the total operators respectively.

[0054] Optionally, in a specific embodiment, Figure 4 As shown, runoff generation includes runoff generation based on area. The pre-configured rainfall runoff conditions include rainfall, evaporation, and infiltration parameters. The runoff generation calculation for the runoff generation area is realized by the following method:

[0055] Step 1041, obtaining the infiltration amount of the first independent hydrological unit based on the infiltration parameter of the first independent hydrological unit;

[0056] Step 1042, obtaining the runoff generation of the area based on the rainfall, evaporation and infiltration.

[0057] Exemplarily, the infiltration rate is calculated by the following formula,

[0058]

[0059] in,

[0060]

[0061] Where f is the infiltration rate, K s is the saturated hydraulic conductivity, h f Soil water suction at the wetting front, θ soil water content, θ s Saturated soil moisture content and antecedent soil moisture content at θ0.

[0062] After obtaining the infiltration amount, the area runoff is obtained by subtracting the infiltration amount and evaporation amount from the rainfall.

[0063] Alternatively, in another specific embodiment, Figure 5 As shown in the figure, the slope confluence and river confluence that produce runoff are achieved by the following methods:

[0064] Step 1044, based on the preset flow area and water flow parameters of the first independent hydrological unit, obtain the flow of the first independent hydrological unit, wherein the first independent hydrological unit is any one of the at least one independent hydrological unit;

[0065] Step 1045, obtaining slope confluence data based on the flow rate and preset slope dynamic parameters.

[0066] Step 1046, obtaining river confluence data based on the flow and preset river dynamic parameters.

[0067] For example, the calculation of runoff in runoff generation and runoff mainly includes the runoff calculation of slope and river channel. The slope runoff and river channel runoff are calculated by the following formula:

[0068] When calculating runoff, you first need to calculate runoff, where runoff is formed by deducting evaporation and infiltration from rainfall. The runoff calculation formula is as follows:

[0069] S=P-ET-f

[0070] Among them, S is surface runoff, P is rainfall, and ET is evaporation.

[0071] When calculating the confluence, the Manning formula and the dynamic wave method are used to calculate the slope confluence and river confluence processes within independent hydrological units respectively.

[0072] Manning's formula:

[0073]

[0074] Among them, Q is the flow rate, A is the flow area, R is the hydraulic radius, S0 is the channel bottom slope, and η is the Manning roughness coefficient.

[0075] Dynamic wave method:

[0076]

[0077] Among them, Q is the flow rate, H is the water head, U is the flow velocity, A is the flow area, S f When calculating the slope confluence and river confluence, you only need to input the corresponding slope, river roughness, river section area and other parameters.

[0078] Step 105 , if after the first preset time, the runoff corresponding to each independent hydrological unit has no hydraulic connection with other independent hydrological units, the runoff data corresponding to each independent hydrological unit are integrated to obtain the hydrological simulation result of the target area.

[0079] Exemplarily, when the calculated runoff of each independent hydrological unit has no hydraulic connection with other independent hydrological units, the runoff of each independent hydrological unit is integrated to obtain the total runoff data in the target area, wherein the runoff integration method is the same as the calculation method of river runoff, both of which are calculated using the dynamic wave method, which will not be repeated here.

[0080] In a specific embodiment, when there are multiple independent hydrological units, and each runoff data calculated after a preset time has a hydraulic connection with other independent hydrological units, such as Figure 6 As shown, the integration of hydrological units that do not conform to hydraulic connections and the subsequent calculation and verification of runoff generation and confluence are achieved through the following methods:

[0081] Step 106, if after the first preset time, there is at least one independent hydrological unit whose corresponding runoff is hydraulically connected with other independent hydrological units, then after the first preset time, the independent hydrological units in the multiple independent hydrological units that are hydraulically connected with each independent hydrological unit in at least one independent hydrological unit are integrated to form a new independent hydrological unit, and the first preset time is any preset time after the start time.

[0082] Step 107, within the second preset time, according to the sub-basin data included in the new independent hydrological unit and the pre-configured rainfall runoff conditions, the corresponding runoff generation and confluence data are obtained.

[0083] The second preset time is the next time period of the first preset time.

[0084] Step 108, if the runoff corresponding to the new independent hydrological unit has no hydraulic connection with other independent hydrological units, the runoff data corresponding to each new independent hydrological unit is integrated with the runoff data of all independent hydrological units in the multiple independent hydrological units that have no hydraulic connection with other independent hydrological units to obtain the hydrological simulation results of the target area.

[0085] For example, when calculating the runoff of each independent hydrological unit, the corresponding runoff within each preset time is calculated and counted, and then the calculated runoff is verified, and then the runoff at the second preset time is calculated. According to the different preset times, the calculation is divided into two paths, one of which is: Figure 7 As shown, each preset time is a unit time; the second case is: Figure 8 As shown, the first preset time is the total simulation time, and the second preset time is the remaining simulation time after the moment when the independent hydrological unit hypothesis is verified to be invalid.

[0086] For the first case, the runoff of each independent hydrological unit is calculated within a unit time (first preset time). After the calculation is completed, it is immediately determined whether the independent hydrological unit assumption between each independent hydrological unit is established (i.e., whether there is a hydraulic connection). When the independent hydrological unit assumption is established, the runoff calculation of the next unit time (second preset time) is performed until the entire simulation time is completed, and the runoff data of the independent hydrological unit are integrated and output. If it is not established, the unestablished hydrological units in the independent hydrological units are reintegrated to obtain new independent hydrological units, and the runoff calculation of the next unit time (second preset time) is performed again according to the area allocation calculation kernel of the reintegrated independent hydrological units until there is no hydraulic connection between all independent hydrological units and the simulation time is completed, and all the runoff data at the current moment are integrated to output the hydrological simulation results in the target area. This method avoids the problem of repeated calculation of runoff of some independent hydrological units caused by the failure of the independent hydrological unit assumption in the second case, and reduces the amount of calculation, but the integration calculation per unit time increases the communication time to a certain extent.

[0087] For the second case, after all simulation time is finished, the runoff data for each unit time is obtained. During verification, the first set of runoff data for the first unit time is verified first. When all runoff calculation results for the first unit time are verified to have no hydraulic connection, the second set of runoff data for the second unit time is verified. When there is a hydraulic connection with other independent hydrological units, the hydrological unit with hydraulic connection with the hydrological unit is reintegrated to obtain a new independent hydrological unit, and the runoff of the new independent hydrological unit in all simulation times (second preset time) after this moment is calculated. The runoff data of each moment are then verified in turn until there is no hydraulic connection between all independent hydrological units. The runoff data of each independent hydrological unit are integrated and the hydrological simulation results of the basin are output. When calculating the runoff of the new independent hydrological unit as mentioned above, it is not necessary to calculate the runoff of the independent hydrological unit without hydraulic connection again. It is only necessary to calculate the runoff of the reintegrated independent hydrological unit. This method can effectively reduce the communication time, but when the number of repeated calculations is large, the simulation calculation amount will increase.

[0088] In a specific embodiment, the verification of the hydraulic connection can be, Figure 3 When the water level at the outflow node of the hydrological unit 1 is lower than the water level of the key node (such as the outlet, storage tank, etc.), water level support occurs at the node and the independent hydrological unit assumption no longer holds; when a flood area appears in the hydrological unit 3, the hydrological units 3 and 4 involved in the flood area no longer meet the independent hydrological unit assumption.

[0089] The reintegration method may be, for example, Figure 3 When the node flooding outflow occurs in hydrological unit 1 as shown in the figure, the free outflow node is formed after integrating the flooding outflow involving hydrological units 1 and 2, and the newly integrated independent hydrological unit is formed as hydrological unit a. For the flood and waterlogging problem in the two-dimensional simulation, the hydrological units within the flooding range are integrated until the newly formed hydrological unit has no hydraulic connection with the other hydrological units. For example, when the confluence is integrated and it is checked that Figure 1 If flooding occurs in the hydrological units 3 and 4, the hydrological units 3 and 4 are integrated, and the integrated independent hydrological unit is hydrological unit b.

[0090] The hydrological simulation calculation method provided by the present invention obtains the watershed data of the target area, and divides the watershed data to obtain multiple sub-watershed data. According to the multiple sub-watershed data, the multiple sub-watershed data are divided into at least one independent water unit, and there is no hydraulic connection between each independent hydrological unit. Taking the current moment as the starting moment, every preset time, according to the sub-watershed data included in each independent hydrological unit, and the pre-configured rainfall runoff conditions, the runoff data corresponding to each independent hydrological unit is obtained. If there is no hydraulic connection between the runoff of each independent hydrological unit and other independent hydrological units, the runoff data of each independent hydrological unit are integrated to output the hydrological simulation result corresponding to the target area.

[0091] In this way, the watershed data of the target area is extracted and divided to obtain multiple sub-watershed data of the target area. Furthermore, at least one independent hydrological unit is created based on the hydraulic connection between the sub-watershed data, so that there is no hydraulic connection between each independent hydrological unit. Since there is no hydraulic connection between each independent hydrological unit, when calculating the runoff corresponding to each independent hydrological unit, there is no need to wait for the calculation of the upstream node to be completed before calculating the runoff of the downstream node, thereby solving the problem of the efficiency limit of parallel computing that exists in the traditional distributed hydrological simulation calculation that relies on a fixed river network topology network.

[0092] The embodiment of the present invention also discloses a hydrological simulation computing device, such as Fig. 9 As shown, the device comprises:

[0093] An acquisition module 201 is used to acquire watershed data of a target area;

[0094] A division module 202, used for dividing the watershed data to obtain a plurality of sub-watershed data;

[0095] The division module 202 further includes:

[0096] The river network extraction module is used to obtain the river network in the target area based on the watershed data.

[0097] The river network water system classification module is used to classify the river network water system and obtain the classified river network water system.

[0098] The sub-basin division module is used to divide the sub-basin according to the classified river network system to obtain multiple sub-basin data.

[0099] An independent hydrological unit creation module 203 is used to create at least one independent hydrological unit according to the multiple sub-basin data, wherein each independent hydrological unit includes at least one sub-basin data among the multiple sub-basin data, and there is no hydraulic connection between any two independent hydrological units at the current moment;

[0100] The first runoff calculation module 204 is used to obtain the runoff data corresponding to each independent hydrological unit according to the sub-basin data included in each independent hydrological unit and the pre-configured rainfall runoff conditions, starting from the current time and after each preset time.

[0101] Optionally, after the independent hydrological units are divided,

[0102] The first flow generation and confluence calculation module 204 further includes:

[0103] An infiltration calculation module, used for obtaining the infiltration amount of the first independent hydrological unit based on the infiltration parameter of the first independent hydrological unit;

[0104] The module for calculating the runoff generation of an area is used to obtain the runoff generation of an area based on rainfall, evaporation and infiltration.

[0105] A flow calculation module, used for obtaining the flow of the first independent hydrological unit based on a preset flow area and water flow parameters of the first independent hydrological unit, wherein the first independent hydrological unit is any one of the at least one independent hydrological unit;

[0106] The confluence calculation module is used to obtain the slope confluence based on the flow and preset slope dynamic parameters; and to obtain the river confluence based on the flow and preset river dynamic parameters.

[0107] The first verification module 205 is used to verify whether the runoff corresponding to each independent hydrological unit has a hydraulic connection with other independent hydrological units after a first preset time has passed. If the runoff corresponding to each independent hydrological unit has a hydraulic connection with other independent hydrological units, the runoff corresponding to each independent hydrological unit is integrated to obtain the hydrological simulation result of the target area.

[0108] A dynamic integration module is used for integrating the independent hydrological units that have hydraulic connections with each of the at least one independent hydrological unit to form a new independent hydrological unit after the first preset time has passed if the verification result after the first preset time is that the runoff corresponding to at least one independent hydrological unit has hydraulic connections with other independent hydrological units, wherein the first preset time is any preset time after the start time;

[0109] A second runoff calculation module is used to obtain the corresponding runoff according to the sub-basin data included in the new independent hydrological unit and the pre-configured rainfall runoff conditions within a second preset time, and the second preset time is the next time period of the first preset time;

[0110] The second verification consistency module is used to integrate the runoff data corresponding to each new independent hydrological unit with the runoff data of all independent hydrological units in multiple independent hydrological units that have no hydraulic connection with other independent hydrological units, so as to obtain the hydrological simulation results of the target area.

[0111] The functions performed by the various components in the hydrological simulation calculation device provided in the embodiment of the present invention have been described in detail in any of the above method embodiments, and therefore will not be repeated here.

[0112] The hydrological simulation calculation device provided by the present invention extracts and divides the watershed data of the target area to obtain multiple sub-watershed data of the target area. Furthermore, at least one independent hydrological unit is created based on the hydraulic connection between the sub-watershed data, so that there is no hydraulic connection between each independent hydrological unit. Since there is no hydraulic connection between each independent hydrological unit, when calculating the runoff corresponding to each independent hydrological unit, it is not necessary to wait for the calculation of the upstream node to be completed before calculating the runoff of the downstream node, thereby solving the defect problem of the parallel computing limit efficiency that the traditional distributed hydrological simulation calculation relies on the fixed river network topology network.

[0113] The embodiment of the present invention also provides a computer device, such as Fig.10 As shown, the computer device may include a processor 301 and a memory 302, wherein the processor 301 and the memory 302 may be connected via a bus or other means. Fig.10 The example of connecting through bus is taken in the following.

[0114] The processor 301 may be a central processing unit (CPU). The processor 301 may also be other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or a combination of the above chips.

[0115] The memory 302 is a non-transient computer-readable storage medium that can be used to store non-transient software programs, non-transient computer executable programs and modules, such as program instructions / modules corresponding to the hydrological simulation calculation method in the embodiment of the present invention. The processor 301 executes various functional applications and data processing of the processor by running the non-transient software programs, instructions and modules stored in the memory 302, that is, the hydrological simulation calculation method in the above method embodiment is implemented.

[0116] The memory 302 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application required by at least one function; the data storage area may store data created by the processor 301, etc. In addition, the memory 302 may include a high-speed random access memory, and may also include a non-volatile memory, such as at least one disk storage device, a flash memory device, or other non-volatile solid-state storage device. In some embodiments, the memory 302 may optionally include a memory remotely arranged relative to the processor 301, and these remote memories may be connected to the processor 301 via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0117] One or more modules are stored in the memory 302, and when executed by the processor 301, the execution is as follows: Figure 1 The hydrological simulation calculation method in the illustrated embodiment.

[0118] For details of the above computer equipment, please refer to Figure 1 The corresponding related descriptions and effects in the illustrated embodiments can be understood and will not be repeated here.

[0119] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program, and the program can be stored in a computer-readable storage medium. When the program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, the storage medium can be a disk, an optical disk, a read-only memory (ROM), a random access memory (RAM), a flash memory, a hard disk drive (HDD) or a solid-state drive (SSD), etc.; the storage medium can also include a combination of the above-mentioned types of memory.

[0120] Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the present invention, and such modifications and variations are all within the scope defined by the appended claims.

Claims

1. A hydrological simulation calculation method, characterized in that: include: Obtain watershed data for the target area; Dividing the watershed data into multiple sub-watershed data; Creating at least one independent hydrological unit according to the plurality of sub-basin data, wherein each of the independent hydrological units includes at least one of the plurality of sub-basin data, and there is no hydraulic connection between any two independent hydrological units at the current moment; Taking the current moment as the starting moment, after each preset time, according to the sub-basin data included in each independent hydrological unit and the pre-configured rainfall runoff conditions, the runoff corresponding to each independent hydrological unit is obtained; If after a first preset time, the runoff data corresponding to each of the independent hydrological units has no hydraulic connection with other independent hydrological units, the runoff data corresponding to each of the independent hydrological units are integrated to obtain a hydrological simulation result of the target area; When the independent hydrological unit includes a plurality of independent hydrological units, the method further includes: If after a first preset time, there is at least one independent hydrological unit whose corresponding runoff is hydraulically connected with other independent hydrological units, then after the first preset time, the independent hydrological units in the plurality of independent hydrological units that are hydraulically connected with each of the at least one independent hydrological unit are integrated to form a new independent hydrological unit, and the first preset time is any preset time after the start time; Within a second preset time, according to the sub-basin data included in the new independent hydrological unit and the pre-configured rainfall runoff conditions, the corresponding runoff data are obtained, and the second preset time is the next time period of the first preset time; If the runoff corresponding to the new independent hydrological unit has no hydraulic connection with other independent hydrological units, the runoff data corresponding to each new independent hydrological unit is integrated with the runoff data of all independent hydrological units in the multiple independent hydrological units that have no hydraulic connection with other independent hydrological units to obtain the hydrological simulation results of the target area.

2. The method according to claim 1, characterized in that The dividing based on the watershed data to obtain a plurality of sub-watershed data specifically includes: Analyze the watershed data to obtain the river network in the target area; Classifying the river network water system to obtain a classified river network water system; Sub-basin division is performed according to the classified river network system to obtain a plurality of sub-basin data.

3. The method according to any one of claims 1 to 2, characterized in that: The runoff generation and runoff include runoff generation by area, the preconfigured rainfall runoff conditions include rainfall, evaporation and infiltration parameters, and the calculation process of runoff generation by area includes: Based on the infiltration parameter of the first independent hydrological unit, the infiltration amount of the first independent hydrological unit is obtained, wherein the first independent hydrological unit is any one of the at least one independent hydrological unit; The runoff generation of the area is obtained based on the rainfall, the evaporation and the infiltration.

4. The method according to any one of claims 1 to 2, characterized in that: The calculation of runoff generation and runoff also includes the calculation of slope runoff data and the calculation of river runoff data. The preconfigured rainfall runoff condition also includes water flow parameters. The calculation process of the slope runoff data and the river runoff data includes: Based on a preset flow area and water flow parameters of a first independent hydrological unit, the flow of the first independent hydrological unit is obtained, wherein the first independent hydrological unit is any one of the at least one independent hydrological unit; Based on the flow rate and preset slope dynamic parameters, the slope confluence data is obtained; Based on the flow rate and preset river dynamic parameters, the river confluence data is obtained.

5. The method according to claim 3, characterized in that: The calculation of runoff generation and runoff also includes the calculation of slope runoff data and the calculation of river runoff data. The preconfigured rainfall runoff condition also includes water flow parameters. The calculation process of the slope runoff data and the river runoff data includes: Based on a preset flow area and water flow parameters of a first independent hydrological unit, the flow of the first independent hydrological unit is obtained, wherein the first independent hydrological unit is any one of the at least one independent hydrological unit; Based on the flow rate and preset slope dynamic parameters, the slope confluence data is obtained; Based on the flow rate and preset river dynamic parameters, the river confluence data is obtained.

6. The method according to claim 1, characterized in that Before calculating the runoff of each of the independent hydrological units, the method further includes: The calculation kernel is allocated according to the area of ​​the first independent hydrological unit to obtain data of the first calculation kernel used to calculate the runoff generation and runoff of the first independent hydrological unit, wherein the first independent hydrological unit is any one of the at least one independent hydrological unit.

7. A hydrological simulation computing device, characterized in that: include: An acquisition module is used to obtain watershed data of a target area; A division module, used for dividing the watershed data to obtain a plurality of sub-watershed data; A module for creating an independent hydrological unit is used to create at least one independent hydrological unit according to a plurality of sub-basin data, wherein each of the independent hydrological units includes at least one of the plurality of sub-basin data, and there is no hydraulic connection between any two independent hydrological units at the current moment; A first runoff calculation module is used to take the current moment as the starting moment, and obtain the runoff corresponding to each of the independent hydrological units according to the sub-basin data included in each of the independent hydrological units and the pre-configured rainfall runoff conditions every preset time; A first verification module is used to verify whether the runoff corresponding to each of the independent hydrological units has a hydraulic connection with other independent hydrological units after a first preset time has passed. If the runoff corresponding to each of the independent hydrological units has no hydraulic connection with other independent hydrological units, the runoff data corresponding to each of the independent hydrological units are integrated to obtain a hydrological simulation result of the target area; When the independent hydrological unit includes a plurality of independent hydrological units, the device further includes: A dynamic integration module, for forming a new independent hydrological unit from among the multiple independent hydrological units the independent hydrological units that have a hydraulic connection with each of the at least one independent hydrological unit, if after a first preset time, the verification result is that there is a hydraulic connection between the runoff and confluence corresponding to at least one independent hydrological unit and other independent hydrological units, wherein the first preset time is any preset time after the start time; A second runoff calculation module is used to obtain corresponding runoff data according to the sub-basin data included in the new independent hydrological unit and the pre-configured rainfall runoff conditions within a second preset time, wherein the second preset time is a time period next to the first preset time; The second verification module is used to verify that the runoff corresponding to the new independent hydrological unit has no hydraulic connection with other independent hydrological units. The runoff data corresponding to each new independent hydrological unit are integrated with the runoff data of all independent hydrological units in the multiple independent hydrological units that have no hydraulic connection with other independent hydrological units to obtain the hydrological simulation results of the target area.

8. A computer device, characterized in that: include: at least one processor; And a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor so that the at least one processor performs the steps of the hydrological simulation calculation method as described in any one of claims 1-6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the hydrological simulation calculation method according to any one of claims 1 to 6 are implemented.

Images