A coupling simulation method and device for surface two-dimensional water power and Green-Ampt method-based infiltration

Abstract

The application provides a coupling simulation method and device for surface two-dimensional water dynamics and infiltration based on a Green-Ampt method, which comprises the following steps: in the absence of consideration of infiltration, performing surface two-dimensional water dynamics calculation of a next calculation period in advance; combining a known current surface water storage to determine a water source of the next calculation period containing infiltration; determining actual infiltration water volume and change of a subsurface water saturation layer of the next calculation period based on the Green-Ampt method, updating a cumulative infiltration water volume and an initial subsurface water saturation rate difference, correcting the surface two-dimensional water dynamics calculation of the next calculation period, and obtaining the surface water storage after the next calculation period. The application can significantly improve the actual effect of the surface two-dimensional water dynamics algorithm in related business applications.

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 and apparatus for two-dimensional hydrodynamics of the Earth's surface and infiltration based on the Green-Ampt 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 the present invention is to provide a coupled simulation method and apparatus for two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method, which can overcome the problem of existing two-dimensional surface hydrodynamic algorithms ignoring infiltration volume and significantly improve the actual effect of two-dimensional surface hydrodynamic algorithms in relevant business applications.

[0005] In a first aspect, the present invention provides a coupled simulation method for two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt 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 of the grid node obtained from the pre-executed two-dimensional hydrodynamic calculation, determine the water source including infiltration for the grid node in the next calculation cycle.

[0008] Based on the Green-Ampt method, the infiltration status of the grid nodes in the next calculation cycle is determined according to whether there is a groundwater saturation layer, the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference between the initial groundwater content and the actual infiltration volume. The infiltration status includes the actual infiltration volume and the change of the groundwater saturation layer.

[0009] Based on whether there is a water-saturated underground layer, the actual infiltration volume in the next calculation cycle, and the current cumulative infiltration volume, determine the difference between the cumulative infiltration volume and the initial groundwater content of the grid nodes after the next calculation cycle;

[0010] Based on the actual infiltration volume in the next calculation cycle, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is revised and executed to obtain the surface water volume after the next calculation cycle.

[0011] In one implementation, based on the Green-Ampt method, the infiltration status of grid nodes in the next calculation cycle is determined according to the presence of a current groundwater saturation layer, the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference between the initial groundwater content and the current infiltration volume. This includes:

[0012] In the absence of a groundwater saturation layer, based on the Green-Ampt method, the remaining time to restore to the initial state is updated according to the current remaining time and the time step of the two-dimensional hydrodynamic algorithm. Then, based on the water source including infiltration, the current cumulative infiltration volume and the difference between the initial groundwater content and the current groundwater content, the actual infiltration volume of the grid node in the next calculation cycle and whether a groundwater saturation layer will be formed in the next calculation cycle are determined.

[0013] If a groundwater saturation layer exists, the remaining time to restore to the initial state is reset to the time specified by the Green-Ampt method. Based on the water source including infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference in initial groundwater content, the actual infiltration volume of the grid node in the next calculation cycle and whether the groundwater saturation layer is lost in the next calculation cycle are determined.

[0014] In one implementation, based on the water source including infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference in initial groundwater content, the actual infiltration volume of the grid node in the next calculation cycle and whether a groundwater-saturated layer will form in the next calculation cycle are determined, including:

[0015] In the absence of a currently saturated underground aquifer;

[0016] If the water source including infiltration is 0 in the next calculation cycle, the actual infiltration volume of the grid node in the next calculation cycle is determined to be 0, and no groundwater saturation layer is formed in the next calculation cycle.

[0017] If the water source including infiltration in the next calculation cycle is not 0 and does not exceed the product of the soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm on the surface, the actual infiltration volume of the grid node in the next calculation cycle is determined to be the water source including infiltration in the next calculation cycle, and no groundwater saturation layer is formed in the next calculation cycle.

[0018] If the water source containing infiltration in the next calculation cycle is not zero and exceeds the product of the soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm, the remaining time to recover to the initial state is reset to the recovery time to the initial state specified by the Green-Ampt method. Based on the Green-Ampt method, the actual infiltration volume of the grid node in the next calculation cycle and whether a groundwater saturation layer is formed in the next calculation cycle are determined according to the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference between the initial groundwater content and the current groundwater content.

[0019] In one implementation, based on the Green-Ampt method, the actual infiltration volume of the grid nodes in the next calculation cycle and whether a groundwater-saturated layer will form in the next calculation cycle are determined according to the water source including infiltration, the current cumulative infiltration volume, and the difference in initial groundwater content. This includes:

[0020] In the absence of a currently saturated underground aquifer;

[0021] If the water source including infiltration in the next calculation cycle is not zero and exceeds the product of the soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm on the surface;

[0022] Based on the Green-Ampt method, the cumulative infiltration water volume that will form a groundwater-saturated layer in the next calculation period is determined according to the infiltrated water source and the current initial groundwater content difference. Then, the following steps are performed in conjunction with the current cumulative infiltration water volume and the initial groundwater content difference:

[0023] If the sum of the current cumulative infiltration volume and the infiltration-inclusive water source in the next calculation period is less than the cumulative infiltration volume that forms a groundwater saturation layer in the next calculation period, the actual infiltration volume in the next calculation period is determined as the infiltration-inclusive water source in the next calculation period, and a groundwater saturation layer is not formed in the next calculation period.

[0024] If the current cumulative infiltration volume is less than the cumulative infiltration volume that will form a groundwater saturation layer in the next calculation period, and the sum of the current cumulative infiltration volume and the infiltrated water source in the next calculation period is not less than the cumulative infiltration volume that will form a groundwater saturation layer in the next calculation period, then it is determined that a groundwater saturation layer will form in the next calculation period. Based on the Green-Ampt method, the updated cumulative infiltration volume is determined according to the difference between the cumulative infiltration volume that will form a groundwater saturation layer in the next calculation period and the current initial groundwater content. The difference between the updated cumulative infiltration volume and the current cumulative infiltration volume is taken as the actual infiltration volume in the next calculation period.

[0025] If the current cumulative infiltration volume is not less than the cumulative infiltration volume that will form a groundwater saturation layer in the next calculation cycle, it is determined that a groundwater saturation layer will form immediately. Based on the Green-Ampt method, the maximum cumulative infiltration volume after the next calculation cycle is determined according to the difference between the current cumulative infiltration volume and the initial groundwater content. Combined with the water source that includes infiltration in the next calculation cycle, the actual infiltration volume in the next calculation cycle and whether the groundwater saturation layer will be lost in the next calculation cycle are determined.

[0026] In one implementation, determining the actual infiltration volume for the next calculation cycle and whether the groundwater saturation layer is lost in the next calculation cycle, based on the infiltrated water source for the next calculation cycle, includes:

[0027] In the absence of a currently saturated underground aquifer;

[0028] If the water source including infiltration in the next calculation cycle is not zero and exceeds the product of the soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm on the surface;

[0029] Provided that the current cumulative infiltration volume is not less than the cumulative infiltration volume that will form the groundwater saturation layer in the next calculation period;

[0030] The difference between the maximum cumulative infiltration volume after the next calculation period determined based on the Green-Ampt method and the current cumulative infiltration volume is taken as the maximum infiltration volume of the next calculation period. This value is then compared with the water source containing infiltration in the next calculation period, and the minimum value is taken as the actual infiltration volume of the next calculation period.

[0031] If the maximum infiltration volume in the next calculation cycle is greater than the infiltration water source in the next calculation cycle, it is determined that the groundwater saturation layer will be lost in the next calculation cycle.

[0032] In one implementation, based on the Green-Ampt method, the actual infiltration volume of the grid node in the next calculation cycle and whether it loses its groundwater saturation layer in the next calculation cycle are determined according to the difference between the infiltrated water source, the current cumulative infiltration volume, and the initial groundwater content in the next calculation cycle. This includes:

[0033] Given the current existence of a water-saturated underground aquifer;

[0034] If the water source containing infiltration in the next calculation cycle is not zero, based on the Green-Ampt method, the maximum cumulative infiltration volume after the next calculation cycle is determined according to the difference between the current cumulative infiltration volume and the initial groundwater content. Combined with the water source containing infiltration in the next calculation cycle, the actual infiltration volume in the next calculation cycle and whether the groundwater saturation layer is lost in the next calculation cycle are determined.

[0035] If the water source including infiltration is 0 in the next calculation cycle, the actual infiltration volume in the next calculation cycle is determined to be 0, and the groundwater saturation layer is lost in the next calculation cycle.

[0036] In one implementation, the difference between the cumulative infiltration volume and the initial groundwater content of the grid node after the next calculation cycle is determined based on whether a groundwater saturation layer exists, the actual infiltration volume in the next calculation cycle, and the current cumulative infiltration volume. This includes:

[0037] If there is no saturated aquifer in the ground and the actual infiltration rate in the next calculation period is 0, the cumulative infiltration rate after the next calculation period is determined based on the current cumulative infiltration rate using the Green-Ampt method. The updated groundwater surface moisture content difference is then determined based on the current groundwater surface moisture content difference using the Green-Ampt method. If the remaining time to return to the initial state is not greater than 0, the cumulative infiltration rate after the next calculation period is set to 0, and the initial groundwater moisture content difference in the Green-Ampt method is set as the updated groundwater surface moisture content difference.

[0038] If there is no saturated groundwater layer at present and the actual infiltration volume in the next calculation period is not zero, or if there is a saturated groundwater layer at present, the sum of the current cumulative infiltration volume and the actual infiltration volume in the next calculation period is taken as the cumulative infiltration volume after the next calculation period. Based on the Green-Ampt method, the updated groundwater surface moisture content difference is determined according to the current groundwater surface moisture content difference.

[0039] In one implementation, based on the known current surface water volume of the grid node and the surface water volume of the grid node obtained from a pre-performed two-dimensional hydrodynamic calculation, the water source including infiltration for the grid node in the next calculation cycle is determined, including:

[0040] If the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is greater than the current surface water volume, then the water source including infiltration for this grid node in the next calculation cycle is determined to be: the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation.

[0041] If the current surface water volume is not 0, and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is less than the current surface water volume, then the water source including infiltration for the grid node in the next calculation cycle is determined to be the arithmetic mean of the current surface water volume and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation.

[0042] 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, then the water source including infiltration for this grid node in the next calculation cycle is determined to be: the current surface water volume.

[0043] If the current surface water volume is 0, and the surface water volume calculated in the pre-executed two-dimensional hydrodynamic calculation is also 0, then the water source including infiltration for this grid node in the next calculation cycle is determined to be 0.

[0044] Secondly, the present invention also provides a coupled simulation device for two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method, comprising:

[0045] The pre-execution module pre-executes the two-dimensional hydrodynamic calculations of the surface for the next calculation cycle without considering infiltration.

[0046] The water source determination module is used to determine the water source, including infiltration, of the grid node 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 from the pre-executed two-dimensional hydrodynamic calculation.

[0047] The infiltration status determination module is used to determine the infiltration status of grid nodes in the next calculation cycle based on the Green-Ampt method, according to whether there is a groundwater saturation layer, the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference in initial groundwater content. The infiltration status includes the actual infiltration volume and the change of the groundwater saturation layer.

[0048] The parameter update module is used to determine the difference between the cumulative infiltration volume and the initial groundwater content of the grid nodes after the next calculation cycle, based on whether there is a groundwater saturation layer, the actual infiltration volume in the next calculation cycle, and the current cumulative infiltration volume.

[0049] The correction execution module is used to correct the surface two-dimensional hydrodynamic calculation for the next calculation cycle based on the actual infiltration volume of the next calculation cycle, and to obtain the surface water volume after the next calculation cycle.

[0050] 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.

[0051] This invention provides a coupled simulation method and apparatus for two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method. First, without considering infiltration, the two-dimensional surface hydrodynamic calculation for the next calculation cycle is performed in advance. Then, based on the known current surface water volume of the grid nodes and the surface water volume of the grid nodes obtained from the pre-performed two-dimensional surface hydrodynamic calculation, the water source containing infiltration for the grid nodes in the next calculation cycle is determined. Next, based on the Green-Ampt method, the simulation is performed according to whether a groundwater-saturated aquifer exists, the water source containing infiltration in the next calculation cycle, and... The difference between the current cumulative infiltration volume and the initial groundwater content determines the infiltration status of the grid node in the next calculation cycle. The infiltration status includes the actual infiltration volume and the change in the groundwater saturation layer. Then, based on whether the groundwater saturation layer exists, the actual infiltration volume in the next calculation cycle, and the current cumulative infiltration volume, the cumulative infiltration volume and the initial groundwater content difference of the grid node after the next calculation cycle are determined. Finally, based on the actual infiltration volume in the next calculation cycle, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is corrected to obtain the surface water volume after the next calculation cycle.

[0052] 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.

[0053] 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.

[0054] 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

[0055] 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.

[0056] Figure 1 A schematic diagram of the process for a coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method, provided for an embodiment of the present invention;

[0057] Figure 2 A technical framework diagram of a coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method provided for embodiments of the present invention;

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

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

[0060] 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.

[0061] Currently, surface water flow inevitably undergoes an infiltration process. However, existing two-dimensional surface hydrodynamic algorithms neglect this infiltration process, resulting in an overestimation of surface water volume in simulations that significantly deviates from reality. This affects the practical value of two-dimensional surface hydrodynamic algorithms in relevant business applications. Based on this, this invention provides a coupled simulation method and apparatus for two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method. This method can overcome the problem of existing two-dimensional surface hydrodynamic algorithms neglecting infiltration volume and significantly improve the actual effect of two-dimensional surface hydrodynamic algorithms in relevant business applications.

[0062] To facilitate understanding of this embodiment, a detailed description of the coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt 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 Green-Ampt method. The method mainly includes the following steps S102 to S110:

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

[0064] 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.

[0065] Step 104: Based on the known 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, determine the water source including infiltration for the grid node in the next calculation cycle.

[0066] In one example, the water source containing infiltration for the grid node in the next calculation cycle can be determined by comparing the known current surface water volume of the grid node with the surface water volume of the grid node obtained from a pre-executed two-dimensional hydrodynamic calculation, and by checking whether both are zero. Specifically, the following cases are considered: the water source containing infiltration is 0 or not 0.

[0067] Step S106: Based on the Green-Ampt method, determine the infiltration status of the grid nodes in the next calculation cycle according to whether there is a groundwater saturation layer, the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume, and the initial groundwater content difference.

[0068] The infiltration status includes the actual infiltration volume and changes in the groundwater saturation layer. In one example, if there is currently no groundwater saturation layer, the remaining time to restore the initial state is updated based on the Green-Ampt method. Then, based on the Green-Ampt method, the remaining time to restore the initial state, the actual infiltration volume for the next calculation cycle, and whether a groundwater saturation layer will form in the next calculation cycle are determined based on the difference between the infiltrated water source, the current cumulative infiltration volume, and the initial groundwater content. Conversely, if there is no groundwater saturation layer, the remaining time to restore the initial state is reset, and the actual infiltration volume for the next calculation cycle is determined based on the difference between the current cumulative infiltration volume and the initial groundwater content, as well as whether a groundwater saturation layer will be lost in the next calculation cycle.

[0069] Step S108: Based on the actual infiltration volume and the current cumulative infiltration volume, determine the difference between the cumulative infiltration volume and the initial groundwater content of the grid node in the next calculation cycle.

[0070] In one example, after the aforementioned step S106, the following scenarios can be identified: 1) There is currently no saturated groundwater layer, and the actual infiltration volume in the next calculation cycle is 0; 2) There is currently no saturated groundwater layer, and the actual infiltration volume in the next calculation cycle is not 0; 3) There is currently a saturated groundwater layer. Based on these scenarios, the Green-Ampt method can be used to determine the cumulative infiltration volume and initial groundwater content difference of the grid nodes after the next calculation cycle in different ways for each of these different scenarios.

[0071] Step S110: Based on the actual infiltration volume of the next calculation cycle, correct the two-dimensional hydrodynamic calculation of the surface for the next calculation cycle to obtain the surface water volume after the next calculation cycle.

[0072] In one example, the external water source conditions of the surface two-dimensional hydrodynamic algorithm are modified based on the actual infiltration volume in the next calculation cycle. 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.

[0073] The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method provided in this invention, compared with existing two-dimensional surface hydrodynamic algorithms, further improves the simulation of 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 method provided by this invention makes existing two-dimensional surface 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.

[0074] To facilitate understanding, this invention provides a core technical approach for a coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method: For each grid node, in each calculation cycle, a two-dimensional surface hydrodynamic calculation is pre-executed without considering infiltration to obtain a water source containing infiltration. Then, the actual infiltration volume is calculated according to the Green-Ampt method. Finally, the actual infiltration volume is set as an external water source condition, and the two-dimensional surface hydrodynamic calculation is corrected, thereby realizing the coupling method of two-dimensional surface hydrodynamics and the Green-Ampt method.

[0075] For specific implementation details, please refer to [link / reference]. Figure 2 The diagram shown illustrates a technical framework for a coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method, including:

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

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

[0078] (III) Based on the known current surface water volume of the grid nodes and the surface water volume of the grid nodes obtained from the pre-executed two-dimensional hydrodynamic calculation, determine the water source, including infiltration, for the grid nodes in the next calculation cycle. Please continue to see... Figure 2 This includes four scenarios: Scenario 1-A, Scenario 1-B, Scenario 1-C, and Scenario 1-D. Specifically:

[0079] Scenario 1-A: If the surface water volume calculated in the pre-executed 2D hydrodynamic calculation is greater than the current surface water volume, then the water source including infiltration for this grid node in the next calculation cycle is determined to be: the surface water volume calculated in the pre-executed 2D hydrodynamic calculation. The expression for the water source including infiltration in the next calculation cycle is:

[0080] ;

[0081] In the formula, For water sources that include 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.

[0082] Scenario 1-B: If the current surface water volume is not zero, and the surface water volume calculated in the pre-executed two-dimensional hydrodynamic calculation is less than the current surface water volume, then the water source including infiltration for this grid node in the next calculation cycle is determined to be the arithmetic mean of the current surface water volume and the surface water volume calculated in the pre-executed two-dimensional hydrodynamic calculation. The expression for the water source including infiltration in the next calculation cycle is:

[0083] ;

[0084] In the formula, For water sources that include 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.

[0085] Scenario 1-C: If the current surface water volume is not zero, and the surface water volume calculated in the pre-executed two-dimensional hydrodynamic calculation is equal to the current surface water volume, then the water source including infiltration for this grid node in the next calculation cycle is determined to be: the current surface water volume. The expression for the water source including infiltration in the next calculation cycle is:

[0086] ;

[0087] In the formula, For water sources that include 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.

[0088] Scenario 1-D: If the current surface water volume is 0, and the surface water volume calculated by the pre-executed two-dimensional hydrodynamic calculation is also 0, then the water source including infiltration for this grid node in the next calculation cycle is determined to be 0.

[0089] (iv) Based on the Green-Ampt method, the actual infiltration volume of the grid nodes and the change in the groundwater saturation layer in the next calculation cycle are determined according to whether there is a saturated groundwater layer, the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference between the initial groundwater content and the current infiltration volume. Please continue to see... Figure 2 This includes two scenarios: Scenario 2-A and Scenario 2-B. Specifically:

[0090] Scenario 2-A: In the absence of a groundwater-saturated layer, based on the Green-Ampt method, the remaining time to restore the initial state is updated according to the current time remaining until the initial state is reached and the time step of the 2D surface hydrodynamic algorithm. Then, based on the water source including infiltration, the current cumulative infiltration volume, and the difference between the initial groundwater content and the current infiltration volume, the actual infiltration volume of the grid nodes in the next calculation cycle and whether a groundwater-saturated layer has formed in the next calculation cycle are determined. Specifically:

[0091] (1) Based on the Green-Ampt method, update (specifically, "reduce") the remaining time to recover to the initial state. Specifically, update the remaining time to recover to the initial state according to the following formula:

[0092] ;

[0093] In the formula, This refers to the remaining time variable from restoring the initial state, which is updated in the Green-Ampt method. In the Green-Ampt method, the remaining time variable is the distance from the current state to the initial state. This represents the time step in the two-dimensional hydrodynamic algorithm for the Earth's surface. In practice, the remaining time to recover to the initial state is substituted into the above formula. The remaining time to restore the initial state is calculated.

[0094] (2) Based on the water source including infiltration, the current cumulative infiltration volume, and the initial groundwater content difference for the next calculation cycle, determine the actual infiltration volume of the grid nodes in the next calculation cycle and whether a groundwater-saturated layer will form in the next calculation cycle. Please continue to see Figure 2 This includes three scenarios: Scenario 3-A, Scenario 3-B, and Scenario 3-C. Specifically:

[0095] Scenario 3-A: If the water source including infiltration is 0 in the next calculation cycle, the actual infiltration volume of the grid node in the next calculation cycle is determined to be 0, and no groundwater saturation layer will form in the next calculation cycle. Furthermore, since the actual infiltration volume in the next calculation cycle is 0, there is no need to modify the external water source conditions subsequently.

[0096] In Scenario 3-B, if the water source including infiltration in the next calculation cycle is not zero and does not exceed the product of the soil saturated hydraulic conductivity and the time step of the 2D surface hydrodynamic algorithm, the actual infiltration volume of the grid node in the next calculation cycle is determined as the water source including infiltration in the next calculation cycle, and no groundwater saturation layer is formed in the next calculation cycle. Furthermore, since the actual infiltration volume in the next calculation cycle is not zero, the external water source conditions need to be modified subsequently.

[0097] In Scenario 3-C, if the water source including infiltration in the next calculation cycle is not zero and exceeds the product of the soil saturated hydraulic conductivity and the time step of the 2D surface hydrodynamic algorithm, the remaining time to restore to the initial state is reset to the time specified by the Green-Ampt method. Based on the Green-Ampt method, the actual infiltration volume of the grid node in the next calculation cycle and whether a groundwater-saturated layer will form in the next calculation cycle are determined according to the water source including infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference between the initial groundwater content and the actual infiltration volume. Furthermore, since the actual infiltration volume in the next calculation cycle is not zero, the external water source conditions need to be modified subsequently.

[0098] Specifically:

[0099] (I) Reset the remaining time to restore to the initial state to the duration of restoration to the initial state specified by the Green-Ampt method.

[0100] (II) Based on the Green-Ampt method, the cumulative infiltration volume forming a groundwater saturation layer in the next calculation period is determined according to the difference between the infiltrated water source and the current initial groundwater content in the next calculation period. Specifically, the cumulative infiltration volume forming a groundwater saturation layer in the next calculation period is determined according to the following formula:

[0101] ;

[0102] In the formula, The cumulative infiltration volume of the underground water-saturated layer is calculated for the next calculation cycle. For saturated hydraulic conductivity, For underground moist frontal suction head. This represents the initial groundwater difference in the Green-Ampt method. For the next calculation cycle, including the infiltrated water source, This represents the time step in the two-dimensional hydrodynamic algorithm for the Earth's surface. In practice, the water source including infiltration from the next calculation cycle is substituted into the above formula. The cumulative infiltration volume of the underground water-saturated layer that will form in the next calculation cycle is calculated.

[0103] (III) For any grid node, based on the Green-Ampt method, the actual infiltration volume of the grid node in the next calculation cycle and whether a groundwater saturation layer will form in the next calculation cycle are determined according to the water source including infiltration in the next calculation cycle, the current cumulative infiltration volume and the initial groundwater content difference, and the cumulative infiltration volume that will form a groundwater saturation layer in the next calculation cycle. Please continue to see... Figure 2 This includes three scenarios: Scenario 4-A, Scenario 4-B, and Scenario 4-C. Specifically:

[0104] Scenario 4-A: If the sum of the current cumulative infiltration volume and the infiltration-including water source in the next calculation period is less than the cumulative infiltration volume that forms a groundwater saturation layer in the next calculation period, then the actual infiltration volume in the next calculation period is determined as the infiltration-including water source in the next calculation period, and a groundwater saturation layer does not form in the next calculation period.

[0105] Scenario 4-B: If the current cumulative infiltration volume is less than the cumulative infiltration volume that will form a groundwater-saturated layer in the next calculation period, and the sum of the current cumulative infiltration volume and the infiltration-including water source in the next calculation period is not less than the cumulative infiltration volume that will form a groundwater-saturated layer in the next calculation period, then it is determined that a groundwater-saturated layer will form in the next calculation period. Based on the Green-Ampt method, the updated cumulative infiltration volume is determined according to the difference between the cumulative infiltration volume that will form a groundwater-saturated layer in the next calculation period and the current initial groundwater content. The difference between the updated cumulative infiltration volume and the current cumulative infiltration volume is taken as the actual infiltration volume in the next calculation period.

[0106] Specifically, the updated cumulative infiltration volume is determined according to the following formula:

[0107] ;

[0108] In the formula, For the updated cumulative infiltration volume, For saturated hydraulic conductivity, For time step variable, This represents the current cumulative infiltration volume. For underground moist frontal suction head. This represents the initial groundwater content difference in the Green-Ampt method. In practice, the cumulative infiltration volume of the groundwater-saturated layer in the next calculation cycle, calculated above, is substituted into the above formula. ,Will Substitute the time step variable into the above formula ( The time step of the two-dimensional hydrodynamic algorithm for the Earth's surface. The cumulative infiltration volume of the underground water-saturated layer is calculated for the next calculation cycle. This represents the current cumulative infiltration volume. The updated cumulative infiltration volume is calculated for the next calculation cycle (including the infiltrated water source). The actual infiltration volume for the next calculation cycle is the difference between the updated cumulative infiltration volume and the current cumulative infiltration volume.

[0109] Scenario 4-C: If the current cumulative infiltration volume is not less than the cumulative infiltration volume that will form a groundwater saturation layer in the next calculation cycle, it is determined that a groundwater saturation layer will form immediately. Based on the Green-Ampt method, the maximum cumulative infiltration volume after the next calculation cycle is determined according to the difference between the current cumulative infiltration volume and the initial groundwater content. Combined with the water source that includes infiltration in the next calculation cycle, the actual infiltration volume in the next calculation cycle and whether the groundwater saturation layer will be lost in the next calculation cycle are determined.

[0110] In one example, the process of determining the actual infiltration volume and whether the groundwater saturation layer is lost in the next calculation cycle, based on the infiltrated water source in the next calculation cycle, is as follows: Using the Green-Ampt method, the maximum cumulative infiltration volume after the next calculation cycle is calculated. The difference between the maximum cumulative infiltration volume after the next calculation cycle and the current cumulative infiltration volume is taken as the maximum infiltration volume for the next calculation cycle. This difference is then compared with the infiltrated water source in the next calculation cycle, and the minimum value is taken as the actual infiltration volume for the next calculation cycle. If the maximum infiltration volume for the next calculation cycle calculated using the Green-Ampt method is greater than the infiltrated water source in the next calculation cycle, it is determined that the groundwater saturation layer is lost in the next calculation cycle.

[0111] Specifically, based on the Green-Ampt method, the difference between the current cumulative infiltration volume and the initial groundwater content is substituted into the formula in scenario 4-B, and the time step of the two-dimensional surface hydrodynamic algorithm is adjusted. Substituting the time step variable into the formula in scenario 4-B The maximum cumulative infiltration volume after the next calculation cycle is calculated. Based on the Green-Ampt method, the maximum infiltration volume for the next calculation cycle is the difference between the maximum cumulative infiltration volume after the next calculation cycle and the current cumulative infiltration volume. This difference is then compared with the infiltrated water sources in the next calculation cycle, and the minimum value is taken as the actual infiltration volume for the next calculation cycle. If the maximum infiltration volume for the next calculation cycle calculated based on the Green-Ampt method is not greater than the infiltrated water sources in the next calculation cycle, the groundwater saturation layer is not lost in the next calculation cycle; if the maximum infiltration volume for the next calculation cycle calculated based on the Green-Ampt method is greater than the infiltrated water sources in the next calculation cycle, the groundwater saturation layer is lost in the next calculation cycle.

[0112] Scenario 2-B: Given the current presence of a saturated groundwater layer, the remaining time to restore to the initial state is reset to the duration specified by the Green-Ampt method. Based on the water source including infiltration, the current cumulative infiltration volume, and the initial groundwater content difference in the next calculation cycle, the actual infiltration volume of the grid nodes in the next calculation cycle and whether the saturated groundwater layer is lost in the next calculation cycle are determined. Specifically:

[0113] (1) Reset the remaining time to restore to the initial state to the time specified by the Green-Ampt method to restore to the initial state.

[0114] (2) Based on the Green-Ampt method, the actual infiltration volume of the grid nodes in the next calculation cycle is determined according to the water source including infiltration, the current cumulative infiltration volume, and the initial groundwater content difference. It is also determined whether the groundwater saturation layer is lost in the next calculation cycle. Please continue to see... Figure 2 This includes two scenarios: Scenario 5-A and Scenario 5-B. Specifically:

[0115] Scenario 5-A: If the water source including infiltration in the next calculation cycle is not zero, based on the Green-Ampt method, the maximum cumulative infiltration volume after the next calculation cycle is determined according to the difference between the current cumulative infiltration volume and the initial groundwater content. Combined with the water source including infiltration in the next calculation cycle, the actual infiltration volume in the next calculation cycle and whether the groundwater saturation layer is lost in the next calculation cycle are determined. Furthermore, since the actual infiltration volume in the next calculation cycle is not zero, the external water source conditions need to be modified subsequently.

[0116] Specifically:

[0117] Based on the Green-Ampt method, the difference between the current cumulative infiltration volume and the initial groundwater content is substituted into the formula in scenario 4-B, and the time step of the two-dimensional surface hydrodynamic algorithm is adjusted. Substituting the time step variable into the formula in scenario 4-B The maximum cumulative infiltration volume after the next calculation cycle is calculated. Based on the Green-Ampt method, the maximum infiltration volume after the next calculation cycle is the maximum cumulative infiltration volume after the next calculation cycle minus the current cumulative infiltration volume. This is then compared with the infiltrated water sources in the next calculation cycle, and the minimum value is taken as the actual infiltration volume for the next calculation cycle. If the maximum infiltration volume after the next calculation cycle calculated based on the Green-Ampt method is greater than the infiltrated water sources in the next calculation cycle, then the groundwater saturation layer is lost.

[0118] In Scenario 5-B, if the water source including infiltration is 0 in the next calculation cycle, the actual infiltration volume in the next calculation cycle is determined to be 0, and the groundwater saturation layer is lost in the next calculation cycle. Furthermore, since the actual infiltration volume in the next calculation cycle is 0, there is no need to modify the external water source conditions subsequently.

[0119] (v) Modify the external water source conditions of the surface two-dimensional hydrodynamic algorithm based on the actual infiltration volume. In specific implementation, the actual infiltration volume of the next calculation cycle is regarded as the external water source, and it is accumulated with other external water sources in the next calculation cycle. The external water source conditions of the surface two-dimensional hydrodynamic algorithm for this grid node in the next calculation cycle are modified as follows:

[0120] ;

[0121] 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, within each calculation cycle, and also ignores the specific changes in the surface water volume within 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 current surface water volume at the beginning and end of each calculation cycle.

[0122] (vi) Based on whether there is a groundwater saturation layer, the actual infiltration volume in the next calculation cycle, and the current cumulative infiltration volume, determine the cumulative infiltration volume and the difference in initial groundwater content of the grid nodes after the next calculation cycle.

[0123] If there is currently no saturated aquifer and the actual infiltration rate for the next calculation cycle is 0, based on the Green-Ampt method, the cumulative infiltration rate after the next calculation cycle is determined according to the current cumulative infiltration rate, and the surface water content difference is updated based on the Green-Ampt method. If the remaining time to restore the initial state is not greater than 0, the cumulative infiltration rate after the next calculation cycle is set to 0, and the initial surface water content difference in the Green-Ampt method is set to the updated surface water content difference. Specifically:

[0124] In one example, based on the Green-Ampt method, the cumulative infiltration volume after the next calculation period is determined according to the following formula:

[0125] ;

[0126] In the formula, This represents the cumulative infiltration volume after the next calculation period. This represents the current cumulative infiltration volume. This is the infiltration capacity recovery coefficient in the Green-Ampt method. This represents the maximum groundwater difference in the Green-Ampt method. The time step of the two-dimensional hydrodynamic algorithm for the Earth's surface. This represents the thickness of the underground surface layer. In practice, the current cumulative infiltration volume is substituted into the above formula. The cumulative infiltration volume after the next calculation cycle is calculated.

[0127] In one example, based on the Green-Ampt method, the difference in subsurface surface water content is updated according to the following formula:

[0128] ;

[0129] In the formula, For the updated difference in surface water content, This represents the current difference in surface water content. This is the infiltration capacity recovery coefficient in the Green-Ampt method. This represents the maximum groundwater difference in the Green-Ampt method. This represents the time step in the two-dimensional hydrodynamic algorithm for the Earth's surface. In practice, the current difference in surface water content is substituted into the above formula. The updated groundwater content difference is calculated. If there is no saturated groundwater layer, the actual infiltration volume in the next calculation cycle is 0, and the remaining time to return to the initial state is no greater than 0, the cumulative infiltration volume after the next calculation cycle is determined to be 0, and the initial groundwater content difference in the Green-Ampt method is set as the updated groundwater content difference.

[0130] If there is no saturated aquifer in the current calculation period and the actual infiltration volume in the next calculation period is not zero, or if there is a saturated aquifer in the current calculation period, the sum of the current cumulative infiltration volume and the actual infiltration volume in the next calculation period is taken as the cumulative infiltration volume after the next calculation period. Based on the Green-Ampt method, the updated groundwater surface moisture content difference is determined according to the current difference. In one example, based on the Green-Ampt method, the updated groundwater surface moisture content difference is determined according to the following formula:

[0131] ;

[0132] In the formula, For the updated difference in surface water content, This represents the current difference in surface water content. This represents the actual infiltration volume for the next calculation cycle. This represents the thickness of the underground surface layer. In practice, the current difference in the underground surface water content is substituted into the above formula. The updated difference in surface water content was calculated.

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

[0134] 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.

[0135] 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.

[0136] Repeatedly executing the above algorithm process constitutes the complete coupled calculation process of surface two-dimensional hydrodynamics and Green-Ampt infiltration.

[0137] In summary, the coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method provided in this embodiment of the invention includes at least the following key technical points:

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

[0139] 2. Coupled with the Green-Ampt infiltration algorithm, at each time point and at all locations, the surface water volume calculated by the two-dimensional hydrodynamics is corrected based on the calculated infiltration volume.

[0140] 3. In the Green-Ampt infiltration 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 Green-Ampt method and the surface two-dimensional hydrodynamic algorithm.

[0141] Based on the aforementioned key technologies, this invention enables the calculation of actual infiltration volume at every moment and at all locations, and corrects the surface water accumulation at grid nodes obtained from pre-executed two-dimensional hydrodynamic calculations based on the infiltration volume. Compared to existing two-dimensional hydrodynamic algorithms, 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, thereby addressing the significant impact of infiltration on surface water accumulation. The method provided by this invention makes existing two-dimensional hydrodynamic algorithms more realistic in simulating surface water flow evolution, significantly enhancing its application value in related business applications.

[0142] Based on the foregoing embodiments, this invention provides a coupled simulation device for two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt 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 Green-Ampt method. The device mainly includes the following parts:

[0143] Pre-execution module 302 pre-executes the two-dimensional hydrodynamic calculation of the surface for the next calculation cycle without considering infiltration;

[0144] The water source determination module 304 is used to determine the water source including infiltration of the grid node 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 pre-executed two-dimensional hydrodynamic calculation of the surface.

[0145] The infiltration status determination module 306 is used to determine the infiltration status of the grid nodes in the next calculation cycle based on the Green-Ampt method, according to whether there is a groundwater saturation layer, the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume and the difference in initial groundwater content. The infiltration status includes the actual infiltration volume and the change of the groundwater saturation layer.

[0146] The parameter update module 308 is used to determine the difference between the cumulative infiltration volume and the initial groundwater content of the grid node after the next calculation cycle based on whether there is a groundwater saturation layer, the actual infiltration volume in the next calculation cycle, and the current cumulative infiltration volume.

[0147] The correction execution module 310 is used to correct the two-dimensional hydrodynamic calculation of the surface in the next calculation cycle based on the actual infiltration volume of the next calculation cycle, so as to obtain the surface water volume after the next calculation cycle.

[0148] The coupled simulation device for two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method provided in this invention, compared with existing two-dimensional surface hydrodynamic algorithms, further improves the simulation of 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 method provided by this invention makes existing two-dimensional surface hydrodynamic algorithms more consistent with real-world conditions in simulating surface water flow evolution, significantly enhancing its application value in related business applications.

[0149] In one implementation, the infiltration state determination module 306 is specifically used for:

[0150] In the absence of a groundwater saturation layer, based on the Green-Ampt method, the remaining time to restore to the initial state is updated according to the current remaining time and the time step of the two-dimensional hydrodynamic algorithm. Then, based on the water source including infiltration, the current cumulative infiltration volume and the difference between the initial groundwater content and the current groundwater content, the actual infiltration volume of the grid node in the next calculation cycle and whether a groundwater saturation layer will be formed in the next calculation cycle are determined.

[0151] If a groundwater saturation layer exists, the remaining time to restore to the initial state is reset to the time specified by the Green-Ampt method. Based on the water source including infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference in initial groundwater content, the actual infiltration volume of the grid node in the next calculation cycle and whether the groundwater saturation layer is lost in the next calculation cycle are determined.

[0152] In one implementation, the infiltration state determination module 306 is specifically used for:

[0153] In the absence of a currently saturated underground aquifer;

[0154] If the water source including infiltration is 0 in the next calculation cycle, the actual infiltration volume of the grid node in the next calculation cycle is determined to be 0, and no groundwater saturation layer is formed in the next calculation cycle.

[0155] If the water source including infiltration in the next calculation cycle is not 0 and does not exceed the product of the soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm on the surface, the actual infiltration volume of the grid node in the next calculation cycle is determined to be the water source including infiltration in the next calculation cycle, and no groundwater saturation layer is formed in the next calculation cycle.

[0156] If the water source containing infiltration in the next calculation cycle is not zero and exceeds the product of the soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm, the remaining time to recover to the initial state is reset to the recovery time to the initial state specified by the Green-Ampt method. Based on the Green-Ampt method, the actual infiltration volume of the grid node in the next calculation cycle and whether a groundwater saturation layer is formed in the next calculation cycle are determined according to the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference between the initial groundwater content and the current groundwater content.

[0157] In one implementation, the infiltration state determination module 306 is specifically used for:

[0158] In the absence of a currently saturated underground aquifer;

[0159] If the water source including infiltration in the next calculation cycle is not zero and exceeds the product of the soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm on the surface;

[0160] Based on the Green-Ampt method, the cumulative infiltration water volume that will form a groundwater-saturated layer in the next calculation period is determined according to the infiltrated water source and the current initial groundwater content difference. Then, the following steps are performed in conjunction with the current cumulative infiltration water volume and the initial groundwater content difference:

[0161] If the sum of the current cumulative infiltration volume and the infiltration-inclusive water source in the next calculation period is less than the cumulative infiltration volume that forms a groundwater saturation layer in the next calculation period, the actual infiltration volume in the next calculation period is determined as the infiltration-inclusive water source in the next calculation period, and a groundwater saturation layer is not formed in the next calculation period.

[0162] If the current cumulative infiltration volume is less than the cumulative infiltration volume that will form a groundwater saturation layer in the next calculation period, and the sum of the current cumulative infiltration volume and the infiltrated water source in the next calculation period is not less than the cumulative infiltration volume that will form a groundwater saturation layer in the next calculation period, then it is determined that a groundwater saturation layer will form in the next calculation period. Based on the Green-Ampt method, the updated cumulative infiltration volume is determined according to the difference between the cumulative infiltration volume that will form a groundwater saturation layer in the next calculation period and the current initial groundwater content. The difference between the updated cumulative infiltration volume and the current cumulative infiltration volume is taken as the actual infiltration volume in the next calculation period.

[0163] If the current cumulative infiltration volume is not less than the cumulative infiltration volume that will form a groundwater saturation layer in the next calculation cycle, it is determined that a groundwater saturation layer will form immediately. Based on the Green-Ampt method, the maximum cumulative infiltration volume after the next calculation cycle is determined according to the difference between the current cumulative infiltration volume and the initial groundwater content. Combined with the water source that includes infiltration in the next calculation cycle, the actual infiltration volume in the next calculation cycle and whether the groundwater saturation layer will be lost in the next calculation cycle are determined.

[0164] In one implementation, the infiltration state determination module 306 is specifically used for:

[0165] In the absence of a currently saturated underground aquifer;

[0166] If the water source including infiltration in the next calculation cycle is not zero and exceeds the product of the soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm on the surface;

[0167] Provided that the current cumulative infiltration volume is not less than the cumulative infiltration volume that will form the groundwater saturation layer in the next calculation period;

[0168] The difference between the maximum cumulative infiltration volume after the next calculation period determined based on the Green-Ampt method and the current cumulative infiltration volume is taken as the maximum infiltration volume of the next calculation period. This is compared with the water source containing infiltration in the next calculation period, and the minimum value is taken as the actual infiltration volume of the next calculation period.

[0169] If the maximum infiltration volume in the next calculation cycle is greater than the infiltration water source in the next calculation cycle, it is determined that the groundwater saturation layer will be lost in the next calculation cycle.

[0170] In one implementation, the infiltration state determination module 306 is specifically used for:

[0171] Given the current existence of a water-saturated underground aquifer;

[0172] If the water source containing infiltration in the next calculation cycle is not zero, based on the Green-Ampt method, the maximum cumulative infiltration volume after the next calculation cycle is determined according to the difference between the current cumulative infiltration volume and the initial groundwater content. Combined with the water source containing infiltration in the next calculation cycle, the actual infiltration volume in the next calculation cycle and whether the groundwater saturation layer is lost in the next calculation cycle are determined.

[0173] If the water source including infiltration is 0 in the next calculation cycle, the actual infiltration volume in the next calculation cycle is determined to be 0, and the groundwater saturation layer is lost in the next calculation cycle.

[0174] In one implementation, the parameter update module 308 is specifically used for:

[0175] If there is no saturated aquifer in the ground and the actual infiltration rate in the next calculation period is 0, the cumulative infiltration rate after the next calculation period is determined based on the current cumulative infiltration rate using the Green-Ampt method. The updated groundwater surface moisture content difference is then determined based on the current groundwater surface moisture content difference using the Green-Ampt method. If the remaining time to return to the initial state is not greater than 0, the cumulative infiltration rate after the next calculation period is set to 0, and the initial groundwater moisture content difference in the Green-Ampt method is set as the updated groundwater surface moisture content difference.

[0176] If there is no saturated groundwater layer at present and the actual infiltration volume in the next calculation period is not zero, or if there is a saturated groundwater layer at present, the sum of the current cumulative infiltration volume and the actual infiltration volume in the next calculation period is taken as the cumulative infiltration volume after the next calculation period. Based on the Green-Ampt method, the updated groundwater surface moisture content difference is determined according to the current groundwater surface moisture content difference.

[0177] In one embodiment, the water source determination module 304 is specifically used for:

[0178] If the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is greater than the current surface water volume, then the water source including infiltration for this grid node in the next calculation cycle is determined to be: the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation.

[0179] If the current surface water volume is not 0, and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation is less than the current surface water volume, then the water source including infiltration for the grid node in the next calculation cycle is determined to be the arithmetic mean of the current surface water volume and the surface water volume obtained from the pre-executed two-dimensional hydrodynamic calculation.

[0180] 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, then the water source including infiltration for this grid node in the next calculation cycle is determined to be: the current surface water volume.

[0181] If the current surface water volume is 0, and the surface water volume calculated in the pre-executed two-dimensional hydrodynamic calculation is also 0, then the water source including infiltration for this grid node in the next calculation cycle is determined to be 0.

[0182] 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.

[0183] 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.

[0184] Figure 4This is a schematic diagram of the structure of an electronic device provided in an embodiment of the present invention. The electronic device 100 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.

[0185] 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.

[0186] 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.

[0187] 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.

[0188] 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.

[0189] 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.

[0190] 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.

[0191] 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 two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt 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 pre-executed two-dimensional hydrodynamic calculation of the next calculation cycle, the water source including infiltration of the grid node in the next calculation cycle is determined. Based on the Green-Ampt method, the infiltration status of the grid node in the next calculation cycle is determined according to whether there is a groundwater saturation layer, the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume and the difference in initial groundwater content. The infiltration status includes the actual infiltration volume and the change of the groundwater saturation layer. Based on the actual infiltration volume and the current cumulative infiltration volume, determine the difference between the cumulative infiltration volume and the initial groundwater content of the grid node after the next calculation cycle; Based on the actual infiltration volume, the surface two-dimensional hydrodynamic calculation for the next calculation cycle is corrected to obtain the surface water volume after the next calculation cycle. Based on the Green-Ampt method, the infiltration status of the grid node in the next calculation cycle is determined according to the presence of a water-saturated aquifer, the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference between the initial groundwater content and the current infiltration volume. This includes: In the absence of the aforementioned groundwater saturation layer, based on the Green-Ampt method, the remaining time to restore to the initial state is updated according to the current remaining time and the time step of the two-dimensional hydrodynamic algorithm on the surface. Then, based on the water source including infiltration, the current cumulative infiltration volume, and the difference between the initial groundwater content and the current infiltration volume in the next calculation cycle, the actual infiltration volume of the grid node in the next calculation cycle and whether the aforementioned groundwater saturation layer has been formed in the next calculation cycle are determined. If the underground water-saturated layer is currently present, the remaining time to restore to the initial state is reset to the time specified by the Green-Ampt method for restoring to the initial state. Based on the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference between the initial groundwater content, the actual infiltration volume of the grid node in the next calculation cycle and whether the underground water-saturated layer is lost in the next calculation cycle are determined.

2. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method according to claim 1, characterized in that, Based on the difference between the infiltrated water source, the current cumulative infiltrated water volume, and the initial groundwater content in the next calculation cycle, determine the actual infiltrated water volume of the grid node in the next calculation cycle and whether the groundwater saturation layer has formed in the next calculation cycle, including: In the absence of the aforementioned underground water-saturated layer; If the water source containing infiltration is 0 in the next calculation cycle, it is determined that the actual infiltration volume of the grid node in the next calculation cycle is 0, and the groundwater saturation layer does not form in the next calculation cycle. If the water source including infiltration in the next calculation cycle is not 0 and does not exceed the product of the soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm on the surface, the actual infiltration volume of the grid node in the next calculation cycle is determined to be the water source including infiltration in the next calculation cycle, and the groundwater saturation layer is not formed in the next calculation cycle. If the water source containing infiltration in the next calculation cycle is not zero and exceeds the product of the soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm, the remaining time to recover to the initial state is reset to the recovery time to the initial state specified by the Green-Ampt method. Based on the Green-Ampt method, according to the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume, and the difference between the initial groundwater content, the actual infiltration volume of the grid node in the next calculation cycle and whether the groundwater saturation layer is formed in the next calculation cycle are determined.

3. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method according to claim 2, characterized in that, Based on the Green-Ampt method, the actual infiltration volume of the grid node in the next calculation cycle and whether the groundwater saturation layer has formed in the next calculation cycle are determined according to the water source including infiltration, the current cumulative infiltration volume, and the difference in initial groundwater content. This includes: If the water source including infiltration in the next calculation cycle is not 0 and exceeds the product of the soil saturated hydraulic conductivity and the time step of the two-dimensional hydrodynamic algorithm of the surface; Based on the Green-Ampt method, the cumulative infiltration water volume forming the groundwater saturation layer in the next calculation cycle is determined according to the infiltrated water source and the current initial groundwater content difference in the next calculation cycle. Then, the following steps are performed in conjunction with the current cumulative infiltration water volume and the initial groundwater content difference: If the sum of the current cumulative infiltration volume and the infiltration-including water source in the next calculation cycle is less than the cumulative infiltration volume that forms the groundwater saturation layer in the next calculation cycle, then the actual infiltration volume in the next calculation cycle is determined to be the infiltration-including water source in the next calculation cycle, and the groundwater saturation layer does not form in the next calculation cycle. If the current cumulative infiltration volume is less than the cumulative infiltration volume that will form the groundwater saturation layer in the next calculation period, and the sum of the current cumulative infiltration volume and the water source including infiltration in the next calculation period is not less than the cumulative infiltration volume that will form the groundwater saturation layer in the next calculation period, then it is determined that the groundwater saturation layer will form in the next calculation period. Based on the Green-Ampt method, the updated cumulative infiltration volume is determined according to the difference between the cumulative infiltration volume that will form the groundwater saturation layer in the next calculation period and the current initial groundwater content. The difference between the updated cumulative infiltration volume and the current cumulative infiltration volume is taken as the actual infiltration volume in the next calculation period. If the current cumulative infiltration volume is not less than the cumulative infiltration volume that will form the groundwater saturation layer in the next calculation cycle, it is determined that the groundwater saturation layer will form immediately. Based on the Green-Ampt method, the maximum cumulative infiltration volume after the next calculation cycle is determined according to the difference between the current cumulative infiltration volume and the initial groundwater content. Combined with the water source containing infiltration in the next calculation cycle, the actual infiltration volume in the next calculation cycle and whether the groundwater saturation layer is lost in the next calculation cycle are determined.

4. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method according to claim 3, characterized in that, Based on the infiltrated water source in the next calculation cycle, determine the actual infiltrated water volume in the next calculation cycle and whether the groundwater saturation layer is lost in the next calculation cycle, including: Provided that the current cumulative infiltration volume is not less than the cumulative infiltration volume that will form the underground water-saturated layer in the next calculation cycle; The difference between the maximum cumulative infiltration volume after the next calculation cycle determined based on the Green-Ampt method and the current cumulative infiltration volume is taken as the maximum infiltration volume of the next calculation cycle. This difference is then compared with the water source containing infiltration in the next calculation cycle, and the minimum value is taken as the actual infiltration volume of the next calculation cycle. If the maximum infiltration volume in the next calculation cycle is greater than the water source containing infiltration in the next calculation cycle, it is determined that the groundwater saturation layer is lost in the next calculation cycle.

5. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method according to claim 2, characterized in that, Based on the Green-Ampt method, the actual infiltration volume of the grid node in the next calculation cycle and whether it loses the groundwater saturation layer in the next calculation cycle are determined according to the water source including infiltration, the current cumulative infiltration volume, and the difference in initial groundwater content. This includes: Given the presence of the aforementioned underground water-saturated layer; If the water source containing infiltration in the next calculation cycle is not zero, based on the Green-Ampt method, the maximum cumulative infiltration volume after the next calculation cycle is determined according to the difference between the current cumulative infiltration volume and the initial groundwater content. Combined with the water source containing infiltration in the next calculation cycle, the actual infiltration volume of the next calculation cycle and whether the groundwater saturation layer is lost in the next calculation cycle are determined. If the water source containing infiltration in the next calculation cycle is 0, the actual infiltration volume in the next calculation cycle is determined to be 0, and the underground water-saturated layer is lost in the next calculation cycle.

6. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method according to claim 1, characterized in that, Based on whether the aforementioned saturated groundwater layer currently exists, the actual infiltration volume in the next calculation cycle, and the current cumulative infiltration volume, the difference between the cumulative infiltration volume and the initial groundwater content of the grid node after the next calculation cycle is determined, including: If the current groundwater saturation layer does not exist and the actual infiltration volume in the next calculation cycle is 0, based on the Green-Ampt method, the cumulative infiltration volume after the next calculation cycle is determined according to the current cumulative infiltration volume, and based on the Green-Ampt method, the updated groundwater surface moisture content difference is determined according to the current groundwater surface moisture content difference; if the remaining time to restore the initial state is not greater than 0, the cumulative infiltration volume after the next calculation cycle is set to 0, and the initial groundwater content difference in the Green-Ampt method is set to the updated groundwater surface moisture content difference; If the current groundwater saturation layer does not exist and the actual infiltration volume in the next calculation cycle is not zero, or if the current groundwater saturation layer exists, the sum of the current cumulative infiltration volume and the actual infiltration volume in the next calculation cycle is taken as the cumulative infiltration volume after the next calculation cycle. Based on the Green-Ampt method, the updated groundwater surface moisture content difference is determined according to the current groundwater surface moisture content difference.

7. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method according to claim 1, characterized in that, Based on the known 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, the water source including infiltration for the grid node in the next calculation cycle is determined, including: If the surface water volume of the grid node obtained from the pre-executed two-dimensional hydrodynamic calculation of the next calculation cycle is greater than the current surface water volume, then the water source including infiltration for the grid node in the next calculation cycle is determined to be: the surface water volume of the grid node obtained from the pre-executed two-dimensional hydrodynamic calculation of the next calculation cycle. If the current surface water volume is not 0, and the surface water volume of the grid node obtained by the pre-executed two-dimensional hydrodynamic calculation of the next calculation cycle is less than the current surface water volume, then the water source including infiltration for the grid node in the next calculation cycle is determined to be: the arithmetic mean of the current surface water volume and the surface water volume of the grid node obtained by the pre-executed two-dimensional hydrodynamic calculation of the next calculation cycle; If the current surface water volume is not 0, and the surface water volume of the grid node obtained by the pre-executed two-dimensional hydrodynamic calculation of the next calculation cycle is equal to the current surface water volume, then the water source including infiltration of the grid node in the next calculation cycle is determined to be: the current surface water volume. If the current surface water volume is 0, and the surface water volume of the grid node obtained from the pre-executed two-dimensional hydrodynamic calculation of the next calculation cycle is 0, then the water source including infiltration for the grid node in the next calculation cycle is determined to be 0.

8. An apparatus for implementing the coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Green-Ampt method as described in any one of claims 1-7, characterized in that, include: The pre-execution module pre-executes the two-dimensional hydrodynamic calculations of the surface for the next calculation cycle without considering infiltration. The water source determination module is used to determine the water source, including infiltration, of the grid node 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. The infiltration status determination module is used to determine the infiltration status of the grid node in the next calculation cycle based on the Green-Ampt method, according to whether there is a groundwater saturation layer, the water source containing infiltration in the next calculation cycle, the current cumulative infiltration volume, and the initial groundwater content difference. The infiltration status includes the actual infiltration volume and the change of the groundwater saturation layer. The parameter update module is used to determine the difference between the cumulative infiltration volume and the initial groundwater content of the grid node after the next calculation cycle, based on whether the groundwater saturated layer exists, the actual infiltration volume, and the current cumulative infiltration volume. The correction execution module is used to correct the two-dimensional hydrodynamic calculation of the surface in the next calculation cycle based on the actual infiltration volume, so as to obtain the surface water volume 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.

Images