A coupling simulation method for surface two-dimensional water power and infiltration based on SCS curve number method

Abstract

This invention provides a coupled simulation method for surface two-dimensional hydrodynamics and infiltration based on the SCS curve number method. The method includes: for each grid node of the surface two-dimensional hydrodynamic algorithm, based on the known current surface water volume, without considering infiltration, pre-executing the surface two-dimensional hydrodynamic calculation for the next calculation cycle to determine whether there are new water sources and existing water sources including infiltration, and their specific values; calculating the infiltration volume for the new water sources in the next calculation cycle based on the SCS curve number method, and then calculating the infiltration volume for the existing water sources in the next calculation cycle to determine the actual infiltration volume for the next calculation cycle; and correcting the surface two-dimensional hydrodynamic calculation for the next calculation cycle based on the actual infiltration volume to obtain the surface water volume after the next calculation cycle. This invention can overcome the problem of existing surface two-dimensional hydrodynamic algorithms ignoring infiltration volume.

Description

Technical Field

[0001] This invention relates to the field of two-dimensional hydrodynamic algorithms for the Earth's surface, and in particular to a coupled simulation method of two-dimensional hydrodynamics for the Earth's surface and infiltration based on the SCS curve number method. Background Technology

[0002] Existing two-dimensional hydrodynamic algorithms for the Earth's surface discretize a specific surface region into multiple grids. Precipitation is input into the gridded region, or inflow rate is input into a certain boundary of the region. For the grid nodes, numerical methods are used to solve the shallow water equations to simulate the conservation of surface water and the dynamic process of inertial convection, propagation, and diffusion of surface water flow under the influence of gravity and surface resistance, so as to calculate the changes in surface water volume and flow velocity at each grid node.

[0003] However, existing two-dimensional surface hydrodynamic algorithms can only simulate the conservation of surface water volume and the dynamic processes of inertial convection, propagation, and diffusion of surface water flow under the influence of gravity and surface resistance. In reality, surface water flow will inevitably undergo an infiltration process, which existing two-dimensional surface hydrodynamic algorithms ignore. This results in the simulated surface water volume being significantly overestimated and deviating from reality, thus affecting the practical value of the two-dimensional surface hydrodynamic algorithms in relevant business applications. Summary of the Invention

[0004] In view of this, the purpose of this invention is to provide a coupled simulation method of surface two-dimensional hydrodynamics and infiltration based on the SCS curve number method, which can overcome the problem of existing surface two-dimensional hydrodynamic algorithms ignoring infiltration volume and significantly improve the actual effect of surface two-dimensional hydrodynamic algorithms in relevant business applications.

[0005] In a first aspect, the present invention provides a coupled simulation method for surface two-dimensional hydrodynamics and infiltration based on the SCS curve number method, comprising:

[0006] Without considering infiltration, perform the two-dimensional hydrodynamic calculation of the surface in advance for the next calculation cycle;

[0007] Based on the known current surface water volume of the grid node and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation, determine whether the grid node has new water sources containing infiltration and existing water sources containing infiltration in the next calculation cycle.

[0008] In the case of a new water source containing infiltration, based on the SCS curve number method, the first infiltration volume of the grid node in the next calculation cycle is determined according to the new water source containing infiltration. The first infiltration volume is the infiltration volume for the new water source containing infiltration. In the case of an existing water source containing infiltration, the second infiltration volume of the grid node in the next calculation cycle is determined according to the existing water source containing infiltration. The second infiltration volume is the infiltration volume for the existing water source containing infiltration.

[0009] Based on the first infiltration volume, or the second infiltration volume, or the first infiltration volume and the second infiltration volume, determine the actual infiltration volume of the grid node in the next calculation cycle;

[0010] Taking into account the actual infiltration volume, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is modified to obtain the surface water accumulation of the grid nodes after the next calculation cycle.

[0011] In one implementation, based on the SCS curve number method, the first infiltration volume of the grid node in the next calculation cycle is determined according to the newly added water source including infiltration, including:

[0012] The sum of the current cumulative water sources containing infiltration and the newly added water sources containing infiltration in the next calculation cycle is used as the cumulative water sources containing infiltration for the grid node after the next calculation cycle.

[0013] Based on the SCS curve number method, the cumulative infiltration volume of the grid nodes after the next calculation cycle is determined according to the cumulative water source including infiltration after the next calculation cycle and the maximum cumulative infiltration volume set in the SCS curve number method.

[0014] The difference between the cumulative infiltration volume after the next calculation cycle and the current cumulative infiltration volume is used as the first infiltration volume of the grid node in the next calculation cycle.

[0015] In one implementation, determining the second infiltration volume of a grid node in the next calculation cycle based on an existing water source including infiltration includes:

[0016] The product of the preset soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm on the surface is taken as the saturated infiltration volume.

[0017] The minimum value between the existing water source including infiltration and the saturated infiltration volume is used as the second infiltration volume of the grid node in the next calculation cycle.

[0018] In one implementation, determining the actual infiltration volume of a grid node in the next calculation cycle based on a first infiltration volume, or a second infiltration volume, or both the first and second infiltration volumes, includes:

[0019] If there is a new water source with infiltration, and the first infiltration volume is not 0, and there is an existing water source with infiltration, and the second infiltration volume is not 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be the sum of the first infiltration volume and the second infiltration volume.

[0020] If there is a new water source that includes infiltration, and the first infiltration volume is not 0, and there is no existing water source that includes infiltration, and the second infiltration volume is 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be: the first infiltration volume.

[0021] If there is no new water source containing infiltration, the first infiltration volume is 0, there is an existing water source containing infiltration, and the second infiltration volume is not 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be the second infiltration volume.

[0022] In one implementation, the method further includes:

[0023] If there are no new water sources containing infiltration, the first infiltration volume is 0, there are no existing water sources containing infiltration, and the second infiltration volume is 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be 0.

[0024] In one implementation, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is modified based on the actual infiltration volume to obtain the surface water accumulation at the grid nodes after the next calculation cycle, including:

[0025] The external water source conditions of the two-dimensional hydrodynamic algorithm for the surface are modified according to the actual infiltration volume.

[0026] Based on the modified external water source conditions, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is revised to obtain the surface water volume of the grid nodes after the next calculation cycle.

[0027] In one implementation, based on the known current surface water volume of the grid node and the surface water volume obtained from a pre-executed two-dimensional hydrodynamic calculation, it is determined whether the grid node has any new water sources containing infiltration and existing water sources containing infiltration in the next calculation cycle, including:

[0028] If the current surface water volume is not zero, and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is greater than the current surface water volume, it is determined that the grid node has both newly added water sources containing infiltration and existing water sources containing infiltration in the next calculation cycle. The value of the newly added water source containing infiltration is the difference between the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation and the current surface water volume. The value of the existing water source containing infiltration is the current surface water volume.

[0029] If the current surface water volume is 0, and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is greater than the current surface water volume, it is determined that the grid node has a new water source with infiltration in the next calculation cycle and no existing water source with infiltration. The value of the new water source with infiltration is: the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation.

[0030] If the current surface water volume is not zero, and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is less than the current surface water volume, it is determined that the grid node does not have any new water sources containing infiltration in the next calculation cycle, but has existing water sources containing infiltration. The value of the existing water sources containing infiltration is: the arithmetic mean of the current surface water volume and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation.

[0031] If the current surface water volume is not 0, and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is equal to the current surface water volume, it is determined that the grid node does not have any new water sources containing infiltration in the next calculation cycle, but has existing water sources containing infiltration. The value of the existing water sources containing infiltration is: the current surface water volume.

[0032] Given that the current surface water volume and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation are both 0, it is determined that the grid nodes do not have any new water sources containing infiltration or existing water sources containing infiltration in the next calculation cycle.

[0033] Secondly, the present invention also provides a coupled simulation device for surface two-dimensional hydrodynamics and infiltration based on the SCS curve number method, comprising:

[0034] The pre-execution module is used to pre-execute the two-dimensional hydrodynamic calculations of the surface for the next calculation cycle without considering infiltration.

[0035] The water source assessment module is used to determine whether there are new water sources with infiltration and existing water sources with infiltration in the next calculation cycle based on the known current surface water volume of the grid node and the surface water volume of the grid node obtained by the two-dimensional hydrodynamic calculation of the surface of the grid node in the next calculation cycle.

[0036] The infiltration volume determination module is used to determine the first infiltration volume of the grid node in the next calculation cycle based on the SCS curve number method when there is a new water source containing infiltration. The first infiltration volume is the infiltration volume for the new water source containing infiltration. Furthermore, when there is an existing water source containing infiltration, the module determines the second infiltration volume of the grid node in the next calculation cycle based on the existing water source containing infiltration. The second infiltration volume is the infiltration volume for the existing water source containing infiltration.

[0037] The actual infiltration quantity determination module is used to determine the actual infiltration quantity of the grid node in the next calculation cycle based on the first infiltration quantity, or the second infiltration quantity, or the first infiltration quantity and the second infiltration quantity.

[0038] The correction execution module is used to take into account the actual infiltration volume and correct the surface two-dimensional hydrodynamic calculation for the next calculation cycle to obtain the surface water accumulation of the grid nodes after the next calculation cycle.

[0039] Thirdly, the present invention also provides an electronic device including a processor and a memory, the memory storing computer-executable instructions executable by the processor, the processor executing the computer-executable instructions to implement any of the methods provided in the first aspect.

[0040] Fourthly, the present invention also provides a computer-readable storage medium storing computer-executable instructions, which, when invoked and executed by a processor, cause the processor to implement any of the methods provided in the first aspect.

[0041] The present invention provides a coupled simulation method for surface two-dimensional hydrodynamics and infiltration based on the SCS curve number method. For each grid node of the surface two-dimensional hydrodynamic algorithm, in each calculation cycle of the algorithm, firstly, without considering infiltration, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is pre-executed. Then, based on the known current surface water volume of the grid node and the surface water volume obtained from the pre-executed surface two-dimensional hydrodynamic calculation, it is determined whether the grid node has new water sources containing infiltration and existing water sources containing infiltration in the next calculation cycle. Next, if new water sources containing infiltration exist, based on the SCS curve number method, the grid node is determined according to the new water sources containing infiltration. The first infiltration volume of the grid node in the next calculation cycle is the infiltration volume for newly added water sources containing infiltration. Furthermore, if existing water sources containing infiltration exist, the second infiltration volume of the grid node in the next calculation cycle is determined based on these existing water sources. The second infiltration volume is then determined based on either the first or second infiltration volume, or both. Finally, based on the actual infiltration volume, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is adjusted to obtain the surface water accumulation of the grid node after the next calculation cycle. This achieves a coupling method between surface two-dimensional hydrodynamics and the SCS curve number method. Compared to existing two-dimensional hydrodynamic algorithms for the Earth's surface, this invention further improves upon the existing ability to simulate the conservation of surface water volume and the dynamic processes of inertial convection, propagation, and diffusion of surface water flow under the influence of gravity and surface resistance. It also addresses the significant impact of infiltration on surface water accumulation in reality. The method provided by this invention makes existing two-dimensional hydrodynamic algorithms more realistic in simulating the evolution of surface water flow, significantly enhancing its application value in related business applications.

[0042] Other features and advantages of the invention will be set forth in the description which follows, and will be apparent in part from the description, or may be learned by practicing the invention. The objects and other advantages of the invention are realized and obtained in accordance with the structures particularly pointed out in the description, claims and drawings.

[0043] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, preferred embodiments are described below in detail with reference to the accompanying drawings. Attached Figure Description

[0044] To more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the specific embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained from these drawings without creative effort.

[0045] Figure 1 A schematic flowchart of a coupled simulation method for two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method provided in an embodiment of the present invention;

[0046] Figure 2 A technical framework diagram of a coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method provided in an embodiment of the present invention;

[0047] Figure 3 A schematic diagram of the structure of a coupled simulation device for two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method provided in an embodiment of the present invention;

[0048] Figure 4 This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention. Detailed Implementation

[0049] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below in conjunction with the embodiments. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0050] Currently, existing two-dimensional hydrodynamic algorithms for the Earth's surface neglect the infiltration process, resulting in an overestimation of the surface water volume in the simulation calculations, which significantly deviates from the actual situation and affects the practical value of the two-dimensional hydrodynamic algorithms in relevant business applications. Based on this, this invention provides a coupled simulation method of two-dimensional hydrodynamics for the Earth's surface and infiltration based on the SCS curve number method. This method can overcome the problem of existing two-dimensional hydrodynamic algorithms neglecting the infiltration volume and significantly improve the actual effect of the two-dimensional hydrodynamic algorithms in relevant business applications.

[0051] To facilitate understanding of this embodiment, a detailed description of the coupled simulation method for two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method disclosed in this embodiment of the invention will be provided first. (See [link to relevant documentation]). Figure 1 The diagram shows a flowchart of a coupled simulation method for two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method. The method mainly includes the following steps S102 to S110:

[0052] Step S102: Without considering infiltration, perform the two-dimensional hydrodynamic calculation of the surface for the next calculation cycle in advance.

[0053] In one example, the study area contains multiple grid nodes, and the two-dimensional hydrodynamic calculation of the surface for the next calculation cycle is performed in advance without considering infiltration.

[0054] Step S104: Based on the known current surface water volume of the grid node and the surface water volume of the grid node obtained from the two-dimensional hydrodynamic calculation of the surface in the next calculation cycle, determine whether the grid node has new water sources containing infiltration and existing water sources containing infiltration in the next calculation cycle.

[0055] In one example, the comparison between the known current surface water volume and the surface water volume calculated using a pre-executed two-dimensional hydrodynamic method, and whether both are zero, can be used to determine whether there are any new water sources containing infiltration and their specific values ​​in the next calculation cycle, or whether there are any existing water sources containing infiltration and their specific values ​​in the next calculation cycle. Specifically, this can be categorized into the following cases: the existence of both new and existing water sources containing infiltration; the existence of new water sources containing infiltration but no existing water sources containing infiltration; the absence of new water sources containing infiltration but the existence of existing water sources containing infiltration (furthermore, based on the comparison between the known current surface water volume and the pre-executed two-dimensional hydrodynamic method, this can be further subdivided into two cases); and the absence of both new and existing water sources containing infiltration.

[0056] Step S106: If there is a new water source containing infiltration, based on the SCS curve number method, determine the first infiltration volume of the grid node in the next calculation cycle according to the new water source containing infiltration. The first infiltration volume is the infiltration volume for the new water source containing infiltration. Furthermore, if there is an existing water source containing infiltration, determine the second infiltration volume of the grid node in the next calculation cycle according to the existing water source containing infiltration. The second infiltration volume is the infiltration volume for the existing water source containing infiltration.

[0057] In one example, for cases where there is a new water source containing infiltration (including cases where there is a new water source containing infiltration and an existing water source containing infiltration, and cases where there is a new water source containing infiltration but no existing water source containing infiltration), the infiltration volume of the grid node for the new water source containing infiltration in the next calculation cycle (i.e., the first infiltration volume) can be determined based on the SCS (Soil Conservation Service-Curve Number Method). For cases where there is an existing water source containing infiltration (including cases where there is a new water source containing infiltration and an existing water source containing infiltration, and cases where there is no new water source containing infiltration but an existing water source containing infiltration), the infiltration volume of the grid node for the existing water source containing infiltration in the next calculation cycle (i.e., the second infiltration volume) can be determined directly based on the existing water source containing infiltration. For other cases (including cases where there is no new water source containing infiltration and no existing water source containing infiltration), it is not necessary to calculate the first or second infiltration volume.

[0058] Step S108: Determine the actual infiltration volume of the grid node in the next calculation cycle based on the first infiltration volume, or the second infiltration volume, or the first infiltration volume and the second infiltration volume.

[0059] In one instance, for the case where there is a new water source containing infiltration but no existing water source containing infiltration, the first infiltration volume is taken as the actual infiltration volume; for the case where there is no new water source containing infiltration but an existing water source containing infiltration, the second infiltration volume is taken as the actual infiltration volume; for the case where there is a new water source containing infiltration and an existing water source containing infiltration, the sum of the first and second infiltration volumes is taken as the actual infiltration volume; for other cases, the actual infiltration volume is determined to be 0.

[0060] Step S110: Considering the actual infiltration volume, correct the surface two-dimensional hydrodynamic calculation for the next calculation cycle to obtain the surface water accumulation of the grid nodes after the next calculation cycle.

[0061] In one example, the external water source conditions of the two-dimensional hydrodynamic algorithm for the surface are modified according to the actual infiltration volume. Based on the modified external water source conditions, the two-dimensional hydrodynamic calculation of the surface is performed in the next calculation cycle to obtain the surface water accumulation of the grid nodes after the next calculation cycle.

[0062] Repeating steps S102 to S110 constitutes the complete coupled calculation process of surface two-dimensional hydrodynamics and SCS curve number method infiltration.

[0063] The coupled simulation method of surface two-dimensional hydrodynamics and infiltration based on the SCS curve number method provided in this invention, compared with existing surface two-dimensional hydrodynamic algorithms, further improves upon the existing ability to simulate the dynamic processes of surface water conservation and inertial convection, propagation, and diffusion under the influence of gravity and surface resistance. This further addresses the significant impact of infiltration on surface water volume in reality. The method provided by this invention makes existing surface two-dimensional hydrodynamic algorithms more realistic in simulating surface water flow evolution, significantly enhancing its application value in related business applications. For ease of understanding, this invention provides a specific implementation method of the coupled simulation method of surface two-dimensional hydrodynamics and infiltration based on the SCS curve number method, see [link to implementation details]. Figure 2 The diagram shows a technical framework for a coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method.

[0064] (a) For each grid node, the current surface water volume is known.

[0065] (ii) Without considering infiltration, perform the two-dimensional hydrodynamic calculation of the surface for the next calculation cycle in advance.

[0066] (III) Based on the current surface water volume of the grid node and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation, determine whether the grid node has any new water sources containing infiltration or existing water sources containing infiltration in the next calculation cycle. Please continue to the next section. Figure 2 It includes five scenarios: Scenario A, Scenario B, Scenario C, Scenario D, and Scenario E. Specifically:

[0067] In scenario A, if the current surface water volume is not zero, and the surface water volume calculated by the pre-executed two-dimensional hydrodynamic calculation is greater than the current surface water volume, it is determined that the grid node has both a new water source containing infiltration and an existing water source containing infiltration in the next calculation cycle. The value of the new water source containing infiltration is the difference between the surface water volume calculated by the pre-executed two-dimensional hydrodynamic calculation and the current surface water volume. The value of the existing water source containing infiltration is the current surface water volume.

[0068] The expression for the newly added water source that includes infiltration is: ;

[0069] The expression for existing water sources that include infiltration is: ;

[0070] In the formula, For the addition of water sources including infiltration, 'The amount of surface water obtained from a pre-executed two-dimensional hydrodynamic calculation of the surface.' This represents the current surface water volume. This includes existing water sources that have infiltrated.

[0071] In scenario B, if the current surface water volume is 0 and the surface water volume calculated by the pre-executed two-dimensional hydrodynamic calculation is greater than the current surface water volume, it is determined that the grid node has a new water source with infiltration in the next calculation cycle and no existing water source with infiltration. The value of the new water source with infiltration is: the surface water volume calculated by the pre-executed two-dimensional hydrodynamic calculation.

[0072] The expression for the newly added water source that includes infiltration is: ;

[0073] In the formula, For the addition of water sources including infiltration, 'The amount of surface water accumulated is obtained from a pre-executed two-dimensional hydrodynamic calculation of the surface.'

[0074] In scenario C, if the current surface water volume is not zero and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is less than the current surface water volume, it is determined that the grid node does not have any new water sources containing infiltration in the next calculation cycle, but has existing water sources containing infiltration. The value of the existing water sources containing infiltration is the arithmetic mean of the current surface water volume and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation.

[0075] The expression for existing water sources that include infiltration is: ;

[0076] In the formula, 'The amount of surface water obtained from a pre-executed two-dimensional hydrodynamic calculation of the surface.' This represents the current surface water volume. This includes existing water sources that have infiltrated.

[0077] In scenario D, if the current surface water volume is not 0 and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is equal to the current surface water volume, it is determined that the grid node does not have any new water sources containing infiltration in the next calculation cycle, but has existing water sources containing infiltration. The value of the existing water sources containing infiltration is: the current surface water volume.

[0078] The expression for existing water sources that include infiltration is: ;

[0079] In the formula, This represents the current surface water volume. This includes existing water sources that have infiltrated.

[0080] In scenario E, if the current surface water volume and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation are both 0, it is determined that the grid node will not have any new water sources containing infiltration or existing water sources containing infiltration in the next calculation cycle.

[0081] (iv) In the case of a new water source containing infiltration (i.e., scenario A and scenario B), based on the SCS curve number method, the first infiltration volume of the grid node in the next calculation cycle is determined according to the new water source containing infiltration. The first infiltration volume is the infiltration volume for the new water source containing infiltration. The specific process is as follows:

[0082] Based on the SCS curve number method, the cumulative infiltration volume of the grid nodes after the next calculation cycle is determined according to the cumulative water source including infiltration after the next calculation cycle and the maximum cumulative infiltration volume set in the SCS curve number method.

[0083] The difference between the cumulative infiltration volume after the next calculation cycle and the current cumulative infiltration volume is used as the first infiltration volume of the grid node in the next calculation cycle.

[0084] (4.1) The sum of the current cumulative water sources containing infiltration and the newly added water sources containing infiltration in the next calculation cycle is taken as the cumulative water sources containing infiltration of the grid node after the next calculation cycle.

[0085] Specifically, according to the SCS curve number method, the cumulative water source including infiltration after the next calculation cycle is determined as the current cumulative water source including infiltration plus the newly added water source including infiltration in the next calculation cycle.

[0086] (4.2) Based on the SCS curve number method, the cumulative infiltration water source after the next calculation cycle is determined according to the cumulative infiltration water source after the next calculation cycle and the maximum cumulative infiltration water volume set in the SCS curve number method.

[0087] Specifically, according to the SCS curve number method:

[0088] ;

[0089] In the formula, To accumulate the infiltration volume, For the cumulative water source including infiltration, This is the maximum cumulative infiltration volume specified by the SCS curve number method. Substitute the cumulative infiltration water source from the next calculation period into the above formula. The cumulative infiltration volume of the grid node after the next calculation cycle is calculated.

[0090] (4.3) The difference between the cumulative infiltration volume after the next calculation cycle and the current cumulative infiltration volume is taken as the first infiltration volume of the grid node in the next calculation cycle.

[0091] Specifically, the infiltration volume (i.e., the first infiltration volume) for the next calculation cycle for the newly added water source including infiltration is the cumulative infiltration volume after the next calculation cycle in the SCS curve number method, minus the current cumulative infiltration volume in the SCS curve number method.

[0092] (v) In the case of existing water sources containing infiltration (i.e., scenario A, scenario C, and scenario D), the second infiltration volume of the grid node in the next calculation cycle is determined based on the existing water sources containing infiltration. The second infiltration volume is the infiltration volume for existing water sources containing infiltration.

[0093] Specifically, the product of the preset soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm is used as the saturated infiltration volume, and the minimum value between the existing water source including infiltration and the saturated infiltration volume is used as the second infiltration volume of the grid node in the next calculation cycle.

[0094] (vi) Determine the actual infiltration volume of the grid nodes in the next calculation cycle. This includes:

[0095] (6.1) In the case that there is a new water source with infiltration, and the first infiltration volume is not 0, and there is an existing water source with infiltration, and the second infiltration volume is not 0 (i.e. scenario A), the actual infiltration volume of the grid node in the next calculation cycle is determined to be the sum of the first infiltration volume and the second infiltration volume.

[0096] Specifically, if there is a new water source with infiltration in the next calculation period, and the infiltration volume for the new water source with infiltration is not 0, and there is an existing water source with infiltration, and the infiltration volume for the existing water source with infiltration is not 0, then the actual infiltration volume in the next calculation period is: the sum of the infiltration volumes for the new water source with infiltration and the existing water source in the next calculation period.

[0097] (6.2) In the case that there is a new water source with infiltration, and the first infiltration volume is not 0, and there is no existing water source with infiltration, and the second infiltration volume is 0 (i.e. scenario B), the actual infiltration volume of the grid node in the next calculation cycle is determined as: the first infiltration volume.

[0098] Specifically, if there is a new water source with infiltration in the next calculation period, and the infiltration volume for the new water source with infiltration is not 0, and there is no existing water source with infiltration, and the infiltration volume for the existing water source with infiltration is 0, then the actual infiltration volume in the next calculation period is: the infiltration volume for the new water source with infiltration in the next calculation period.

[0099] (6.3) In the case that there is no new water source containing infiltration, the first infiltration volume is 0, there is an existing water source containing infiltration, and the second infiltration volume is not 0 (i.e. scenario C and scenario D), the actual infiltration volume of the grid node in the next calculation cycle is determined to be the second infiltration volume.

[0100] Specifically, if there are no new water sources containing infiltration in the next calculation period, and the infiltration volume for new water sources containing infiltration is 0, and there are existing water sources containing infiltration, and the infiltration volume for existing water sources containing infiltration is not 0, then the actual infiltration volume in the next calculation period is: the infiltration volume for existing water sources containing infiltration in the next calculation period.

[0101] (6.4) In the case that there is no new water source containing infiltration, the first infiltration volume is 0, there is no existing water source containing infiltration, and the second infiltration volume is 0 (i.e. scenario E), the actual infiltration volume of the grid node in the next calculation cycle is determined to be 0.

[0102] (vii) Based on the actual infiltration volume, revise the two-dimensional hydrodynamic calculation of the surface for the next calculation cycle to obtain the surface water volume of the grid nodes after the next calculation cycle.

[0103] Specifically, the external water source conditions of the surface two-dimensional hydrodynamic algorithm are modified according to the actual infiltration volume; combined with the modified external water source conditions, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is corrected to obtain the surface water accumulation of the grid nodes after the next calculation cycle.

[0104] The actual infiltration volume in the next calculation cycle is treated as an external water source and added to other external water sources in the next calculation cycle. The external water source conditions for the surface two-dimensional hydrodynamic algorithm of this grid node in the next calculation cycle are then modified.

[0105] ;

[0106] In the formula, 'For the modified external water source conditions for the next calculation cycle.' The external water source conditions for the next calculation cycle before the modification. This represents the actual infiltration volume for the next calculation cycle. In the two-dimensional hydrodynamic algorithm for the Earth's surface, external water source conditions can be set for any grid node to simulate the phenomenon of underground overflow and surface drainage. The infiltration process causes the loss of surface water, so it can be considered as an external water source. Although the time step of the two-dimensional hydrodynamic algorithm for the Earth's surface can be very small (e.g., 1 second), it is still discrete, not completely continuous. Therefore, the algorithm ignores the specific changes in the external water source, including actual infiltration, in each calculation cycle, and also ignores the specific changes in the surface water volume in each calculation cycle. However, the algorithm can set the changes in the external water source at the beginning and end of each calculation cycle and output the surface water volume at the beginning and end of each calculation cycle.

[0107] If the actual infiltration volume in the next calculation cycle is not zero, then the two-dimensional hydrodynamic calculation of the surface in the next calculation cycle is modified in combination with the modified external water source conditions. The dynamic process of surface water conservation and inertial convection, propagation and diffusion of surface water flow under the action of gravity and surface resistance is simulated to obtain the surface water volume of the grid node after the next calculation cycle.

[0108] If the actual infiltration volume in the next calculation cycle is 0, the two-dimensional hydrodynamic calculation of the surface is repeated for the next calculation cycle to simulate the conservation of surface water and the dynamic process of inertial convection, propagation and diffusion of surface water under the action of gravity and surface resistance, so as to obtain the surface water volume of the grid node after the next calculation cycle.

[0109] Repeatedly executing the above algorithm process constitutes the complete coupled calculation process of surface two-dimensional hydrodynamics and SCS curve number method infiltration.

[0110] In summary, the coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method provided in this invention includes the following key technical points:

[0111] 1. Coupled with a two-dimensional hydrodynamic algorithm for the surface, the infiltration volume is calculated at every time and at all locations.

[0112] 2. Coupled with the SCS curve number method infiltration algorithm, at each time and at all locations, the surface water volume calculated by the two-dimensional hydrodynamics is corrected based on the calculated infiltration volume.

[0113] 3. In the SCS curve number method, in addition to precipitation, the water source factors (water conservation, gravity, surface resistance, inertial convection, propagation, diffusion, underground overflow, surface drainage, etc.) of the surface two-dimensional hydrodynamic algorithm are added to adapt to the coupling of the SCS curve number method and the surface two-dimensional hydrodynamic algorithm.

[0114] Based on this, the method provided by the embodiments of the present invention can calculate the actual amount of infiltration at every moment and at all locations, and correct the surface water volume calculated by the two-dimensional hydrodynamics of the surface based on the infiltration amount. Compared with existing two-dimensional hydrodynamic algorithms, this invention further improves upon the existing ability to simulate the dynamic processes of surface water conservation and the inertial convection, propagation, and diffusion of surface water flow under the influence of gravity and surface resistance, while also addressing the significant impact of infiltration on surface water volume in reality. The method provided by the embodiments of the present invention makes existing two-dimensional hydrodynamic algorithms more consistent with real-world conditions in simulating surface water flow evolution, significantly enhancing its application value in related business applications.

[0115] Based on the foregoing embodiments, this invention provides a coupled simulation device for surface two-dimensional hydrodynamics and infiltration based on the SCS curve number method. See [link to relevant documentation]. Figure 3 The diagram shows a structural schematic of a coupled simulation device for two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method. The device mainly includes the following parts:

[0116] The pre-execution module 302 is used to pre-estimate the surface water volume of the grid nodes in the study area after the next calculation cycle using a two-dimensional hydrodynamic algorithm without considering infiltration.

[0117] The water source assessment module 304 is used to determine whether there are new water sources containing infiltration and existing water sources containing infiltration in the next calculation cycle based on the known current surface water volume and the pre-estimated surface water volume of the grid node.

[0118] The infiltration volume determination module 306 is used to determine the first infiltration volume of the grid node in the next calculation cycle based on the SCS curve number method when there is a new water source containing infiltration. The first infiltration volume is the infiltration volume for the new water source containing infiltration. And, when there is an existing water source containing infiltration, it determines the second infiltration volume of the grid node in the next calculation cycle based on the existing water source containing infiltration. The second infiltration volume is the infiltration volume for the existing water source containing infiltration.

[0119] The actual infiltration quantity determination module 308 is used to determine the actual infiltration quantity of the grid node in the next calculation cycle based on the first infiltration quantity, or the second infiltration quantity, or the first infiltration quantity and the second infiltration quantity.

[0120] The modified execution module 310 is used to take into account the actual infiltration volume and modify the two-dimensional hydrodynamic calculation of the surface in the next calculation cycle to obtain the surface water accumulation of the grid nodes after the next calculation cycle.

[0121] The coupled simulation device for surface two-dimensional hydrodynamics and infiltration based on the SCS curve number method provided in this invention, compared with existing surface two-dimensional hydrodynamic algorithms, further improves the ability to simulate the dynamic processes of surface water conservation and inertial convection, propagation, and diffusion under the influence of gravity and surface resistance. This is based on the existing ability to simulate the dynamic processes of surface water flow under gravity and surface resistance. The algorithm provided by this invention makes existing surface two-dimensional hydrodynamic algorithms more consistent with the real situation in simulating the evolution of surface water flow, and significantly enhances its application value in related business.

[0122] In one embodiment, the infiltration volume determination module 306 is specifically used for:

[0123] The sum of the current cumulative water sources containing infiltration and the newly added water sources containing infiltration in the next calculation cycle is used as the cumulative water sources containing infiltration for the grid node after the next calculation cycle.

[0124] Based on the SCS curve number method, the cumulative infiltration volume of the grid nodes after the next calculation cycle is determined according to the cumulative water source including infiltration after the next calculation cycle and the maximum cumulative infiltration volume set in the SCS curve number method.

[0125] The difference between the cumulative infiltration volume after the next calculation cycle and the current cumulative infiltration volume is used as the first infiltration volume of the grid node in the next calculation cycle.

[0126] In one embodiment, the infiltration volume determination module 306 is specifically used for:

[0127] The product of the preset soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm on the surface is taken as the saturated infiltration volume.

[0128] The minimum value between the existing water source including infiltration and the saturated infiltration volume is used as the second infiltration volume of the grid node in the next calculation cycle.

[0129] In one implementation, the actual infiltration quantity determination module 308 is specifically used for:

[0130] If there is a new water source with infiltration, and the first infiltration volume is not 0, and there is an existing water source with infiltration, and the second infiltration volume is not 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be the sum of the first infiltration volume and the second infiltration volume.

[0131] If there is a new water source that includes infiltration, and the first infiltration volume is not 0, and there is no existing water source that includes infiltration, and the second infiltration volume is 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be: the first infiltration volume.

[0132] If there is no new water source containing infiltration, the first infiltration volume is 0, there is an existing water source containing infiltration, and the second infiltration volume is not 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be the second infiltration volume.

[0133] In one implementation, the actual infiltration quantity determination module 308 is specifically used for:

[0134] If there are no new water sources with infiltration and no existing water sources with infiltration, the actual infiltration volume of the grid node in the next calculation cycle is determined to be 0.

[0135] In one implementation, the correction execution module 310 is specifically used for:

[0136] The external water source conditions of the two-dimensional hydrodynamic algorithm for the surface are modified according to the actual infiltration volume.

[0137] The surface two-dimensional hydrodynamic calculation is revised for the next calculation cycle to obtain the surface water volume of the grid nodes after the next calculation cycle.

[0138] In one implementation, the water source assessment module 304 is specifically used for:

[0139] If the current surface water volume is known to be non-zero, and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is greater than the current surface water volume, then it is determined that the grid node has both a new water source containing infiltration and an existing water source containing infiltration in the next calculation cycle. The value of the new water source containing infiltration is the difference between the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation and the current surface water volume. The value of the existing water source containing infiltration is the current surface water volume.

[0140] Given that the current surface water volume is 0 and the surface water volume calculated by the pre-executed two-dimensional hydrodynamic calculation is greater than the current surface water volume, it is determined that the grid node has a new water source with infiltration in the next calculation cycle and no existing water source with infiltration. The value of the new water source with infiltration is: the surface water volume calculated by the pre-executed two-dimensional hydrodynamic calculation.

[0141] If the current surface water volume is known to be non-zero, and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is less than the current surface water volume, then it is determined that the grid node does not have any new water sources containing infiltration in the next calculation cycle, but has existing water sources containing infiltration. The value of the existing water sources containing infiltration is: the arithmetic mean of the current surface water volume and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation.

[0142] Given that the current surface water volume is not zero, and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is equal to the current surface water volume, it is determined that the grid node does not have any new water sources containing infiltration in the next calculation cycle, but has existing water sources containing infiltration. The value of the existing water sources containing infiltration is: the current surface water volume.

[0143] Given that the current surface water volume and the surface water volume calculated in the pre-executed two-dimensional hydrodynamic calculation are both 0, it is determined that the grid nodes do not have any new water sources containing infiltration or existing water sources containing infiltration in the next calculation cycle.

[0144] The device provided in this embodiment of the invention has the same implementation principle and technical effect as the aforementioned method embodiment. For the sake of brevity, any parts not mentioned in the device embodiment can be referred to the corresponding content in the aforementioned method embodiment.

[0145] This invention provides an electronic device, specifically, the electronic device includes a processor and a storage device; the storage device stores a computer program, and the computer program, when run by the processor, executes the method described in any of the above embodiments.

[0146] Figure 4 The present invention provides a schematic diagram of the structure of an electronic device 100, which includes a processor 40, a memory 41, a bus 42 and a communication interface 43. The processor 40, the communication interface 43 and the memory 41 are connected through the bus 42. The processor 40 is used to execute executable modules, such as computer programs, stored in the memory 41.

[0147] The memory 41 may include high-speed random access memory (RAM) or non-volatile memory, such as at least one disk storage device. Communication between this system network element and at least one other network element is achieved through at least one communication interface 43 (which can be wired or wireless), such as the Internet, wide area network, local area network, metropolitan area network, etc.

[0148] Bus 42 can be an ISA bus, PCI bus, or EISA bus, etc. The bus can be divided into address bus, data bus, control bus, etc. For ease of representation, Figure 4 The symbol is represented by a single double-headed arrow, but this does not mean that there is only one bus or one type of bus.

[0149] The memory 41 is used to store programs. After receiving an execution instruction, the processor 40 executes the program. The method executed by the device for defining the flow process disclosed in any of the foregoing embodiments of the present invention can be applied to the processor 40 or implemented by the processor 40.

[0150] Processor 40 may be an integrated circuit chip with signal processing capabilities. In implementation, each step of the above method can be completed by the integrated logic circuitry in the hardware of processor 40 or by instructions in software form. Processor 40 can be a general-purpose processor, including a Central Processing Unit (CPU), a Network Processor (NP), etc.; it can also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field-Programmable Gate Array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components. It can implement or execute the methods, steps, and logic block diagrams disclosed in the embodiments of this invention. The general-purpose processor can be a microprocessor or any conventional processor. The steps of the methods disclosed in the embodiments of this invention can be directly embodied in the execution of a hardware decoding processor, or executed by a combination of hardware and software modules in the decoding processor. The software modules can reside in random access memory, flash memory, read-only memory, programmable read-only memory, electrically erasable programmable memory, registers, or other mature storage media in the art. The storage medium is located in memory 41. The processor 40 reads the information in memory 41 and, in conjunction with its hardware, completes the steps of the above method.

[0151] The computer program product of the readable storage medium provided in the embodiments of the present invention includes a computer-readable storage medium storing program code. The instructions included in the program code can be used to execute the methods described in the foregoing method embodiments. For specific implementation, please refer to the foregoing method embodiments, which will not be repeated here.

[0152] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of this invention, essentially, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0153] Finally, it should be noted that the above-described embodiments are merely specific implementations of the present invention, used to illustrate the technical solutions of the present invention, and not to limit it. The scope of protection of the present invention is not limited thereto. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments within the technical scope disclosed in the present invention, or make equivalent substitutions for some of the technical features; and these modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention, and should all be covered within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A coupled simulation method for surface two-dimensional hydrodynamics and infiltration based on the SCS curve number method, characterized in that, include: Without considering infiltration, perform the two-dimensional hydrodynamic calculation of the surface in advance for the next calculation cycle; Based on the known current surface water volume of the grid node and the surface water volume of the grid node obtained from the two-dimensional hydrodynamic calculation of the surface in the next calculation cycle, it is determined whether the grid node has new water sources with infiltration and existing water sources with infiltration in the next calculation cycle. In the case of a new water source containing infiltration, based on the SCS curve number method, the first infiltration volume of the grid node in the next calculation cycle is determined according to the new water source containing infiltration, and the first infiltration volume is the infiltration volume for the new water source containing infiltration; and, in the case of an existing water source containing infiltration, the second infiltration volume of the grid node in the next calculation cycle is determined according to the existing water source containing infiltration, and the second infiltration volume is the infiltration volume for the existing water source containing infiltration. The actual infiltration volume of the grid node in the next calculation cycle is determined based on the first infiltration volume, or the second infiltration volume, or the first infiltration volume and the second infiltration volume. Taking into account the actual infiltration volume, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is modified to obtain the surface water volume after the next calculation cycle.

2. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method according to claim 1, characterized in that, Based on the SCS curve number method, and according to the newly added water source including infiltration, the first infiltration volume of the grid node in the next calculation cycle is determined, including: The sum of the current cumulative water sources containing infiltration and the newly added water sources containing infiltration in the next calculation cycle is taken as the cumulative water sources containing infiltration of the grid node after the next calculation cycle. Based on the SCS curve number method, the cumulative infiltration water volume of the grid node after the next calculation cycle is determined according to the cumulative water source including infiltration after the next calculation cycle and the maximum cumulative infiltration water volume set in the SCS curve number method. The difference between the cumulative infiltration volume after the next calculation cycle and the current cumulative infiltration volume is taken as the first infiltration volume of the grid node in the next calculation cycle.

3. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method according to claim 1, characterized in that, Determining the second infiltration volume of the grid node in the next calculation cycle based on the existing water source including infiltration includes: The product of the preset soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm on the surface is taken as the saturated infiltration volume. The minimum value between the existing water source containing infiltration and the saturated infiltration volume is taken as the second infiltration volume of the grid node in the next calculation cycle.

4. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method according to claim 1, characterized in that, Determining the actual infiltration volume of the grid node in the next calculation cycle based on the first infiltration volume, or based on the second infiltration volume, or both the first and second infiltration volumes, includes: If there is a new water source containing infiltration, and the first infiltration volume is not 0, and there is an existing water source containing infiltration, and the second infiltration volume is not 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be the sum of the first infiltration volume and the second infiltration volume. If there is a new water source containing infiltration, and the first infiltration volume is not 0, and there is no existing water source containing infiltration, and the second infiltration volume is 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be: the first infiltration volume. If there is no new water source containing infiltration, and the first infiltration volume is 0, and there is an existing water source containing infiltration, and the second infiltration volume is not 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be: the second infiltration volume.

5. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method according to claim 1, characterized in that, The method further includes: If there is no new water source containing infiltration, and the first infiltration volume is 0, and there is no existing water source containing infiltration, and the second infiltration volume is 0, then the actual infiltration volume of the grid node in the next calculation cycle is determined to be 0.

6. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method according to claim 1, characterized in that, Taking into account the actual infiltration volume, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is modified to obtain the surface water volume after the next calculation cycle, including: The external water source conditions of the two-dimensional hydrodynamic algorithm for the surface are modified according to the actual infiltration volume. Based on the modified external water source conditions, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is corrected to obtain the surface water accumulation of the grid node after the next calculation cycle.

7. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method according to claim 1, characterized in that, Based on the current surface water volume of the grid node and the surface water volume of the grid node obtained from the pre-executed two-dimensional hydrodynamic calculation of the next calculation cycle, it is determined whether the grid node has any new water sources containing infiltration and existing water sources containing infiltration in the next calculation cycle, including: If the current surface water volume is not 0, and the surface water volume calculated in the pre-executed calculation is greater than the current surface water volume, it is determined that the grid node has both a newly added water source containing infiltration and an existing water source containing infiltration in the next calculation cycle. The value of the newly added water source containing infiltration is the difference between the surface water volume calculated in the pre-executed calculation and the current surface water volume, and the value of the existing water source containing infiltration is the current surface water volume. If the current surface water volume is 0 and the surface water volume calculated in the pre-execution is greater than the current surface water volume, it is determined that the grid node has a new water source containing infiltration and no existing water source containing infiltration in the next calculation cycle. The value of the new water source containing infiltration is: the surface water volume calculated in the pre-execution. If the current surface water volume is not 0, and the surface water volume calculated in the pre-execution is less than the current surface water volume, it is determined that the grid node does not have any new water sources containing infiltration in the next calculation cycle, but has existing water sources containing infiltration. The value of the existing water sources containing infiltration is the arithmetic mean of the surface water volume calculated in the pre-execution and the current surface water volume. If the current surface water volume is not 0, and the surface water volume calculated in the pre-executed calculation is equal to the current surface water volume, it is determined that the grid node does not have any new water sources containing infiltration in the next calculation cycle, but has existing water sources containing infiltration. The value of the existing water sources containing infiltration is: the current surface water volume. If the surface water volume calculated in the pre-executed calculation and the current surface water volume are both 0, it is determined that the grid node does not have any new water sources containing infiltration or existing water sources containing infiltration in the next calculation cycle.

8. A coupled simulation device for two-dimensional surface hydrodynamics and infiltration based on the SCS curve number method, characterized in that, include: The pre-execution module is used to pre-execute the two-dimensional hydrodynamic calculations of the surface for the next calculation cycle without considering infiltration. The water source assessment module is used to determine whether the grid node has a new water source with infiltration and an existing water source with infiltration in the next calculation cycle, based on the known current surface water volume of the grid node and the surface water volume of the grid node calculated by the two-dimensional hydrodynamics of the surface in the next calculation cycle. The infiltration volume determination module is used to determine, based on the SCS curve number method, the first infiltration volume of the grid node in the next calculation cycle when there is a new water source containing infiltration, and the first infiltration volume is the infiltration volume for the new water source containing infiltration; and, when there is an existing water source containing infiltration, to determine the second infiltration volume of the grid node in the next calculation cycle, and the second infiltration volume is the infiltration volume for the existing water source containing infiltration. The actual infiltration quantity determination module is used to determine the actual infiltration quantity of the grid node in the next calculation cycle based on the first infiltration quantity, or the second infiltration quantity, or the first infiltration quantity and the second infiltration quantity. The correction execution module is used to take into account the actual infiltration volume, correct the surface two-dimensional hydrodynamic calculation of the next calculation cycle, and obtain the surface water accumulation of the grid node after the next calculation cycle.

9. An electronic device, characterized in that, The method includes a processor and a memory, the memory storing computer-executable instructions executable by the processor, the processor executing the computer-executable instructions to implement the method of any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer-executable instructions that, when invoked and executed by a processor, cause the processor to perform the method according to any one of claims 1 to 7.

Images