A coupling simulation method and device for surface two-dimensional water power and infiltration based on horton method

Abstract

The application provides a surface two-dimensional water power and Horton method-based infiltration coupling simulation method and device, which comprises the following steps: on the basis of a known current surface water storage amount, without considering infiltration, performing surface two-dimensional water power calculation of a next calculation period in advance to determine a water source of the next calculation period containing infiltration; if the water source containing infiltration is not 0, determining an actual infiltration water amount based on the Horton method; if the water source containing infiltration is 0, determining that the actual infiltration water amount is 0; determining a cumulative infiltration water amount after the next calculation period, correcting the surface two-dimensional water power calculation of the next calculation period according to the actual infiltration water amount, obtaining a surface water storage amount after the next calculation period, and repeating the above process until the end of the infiltration simulation process. The application can make up for the problem that the existing surface two-dimensional water power algorithm ignores the infiltration water amount.

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 Horton 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 Horton 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 Horton 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] If the amount of water source containing infiltration is not zero, based on the Horton method, the maximum infiltration amount of the grid node in the next calculation cycle is determined according to the current cumulative infiltration amount, and the actual infiltration amount of the grid node in the next calculation cycle is determined according to the maximum infiltration amount and the water source containing infiltration; if the amount of water source containing infiltration is zero, the actual infiltration amount of the grid node in the next calculation cycle is determined to be zero.

[0009] Based on the actual infiltration volume and the current cumulative infiltration volume, determine the cumulative infiltration volume of the grid nodes after the next calculation cycle;

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

[0011] In one implementation, based on the Horton method, the maximum infiltration volume of the grid node in the next calculation cycle is determined according to the current cumulative infiltration volume, including:

[0012] Based on the Horton method, the time parameter in the Horton method is calculated according to the current cumulative infiltration volume. The sum of the time parameter in the Horton method and the time step of the two-dimensional hydrodynamic algorithm is used as the new time parameter in the Horton method. Then, based on the Horton method, the maximum cumulative infiltration volume of the grid node after the next calculation cycle is determined according to the new time parameter in the Horton method.

[0013] 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 of the grid node in the next calculation cycle.

[0014] In one implementation, based on the Horton method, the maximum infiltration volume of the grid node in the next calculation cycle is determined according to the current cumulative infiltration volume, and the method further includes:

[0015] Based on the improved Horton method, the maximum infiltration rate of the grid node in the next calculation cycle is determined according to the current cumulative infiltration volume and the time parameters of the two-dimensional hydrodynamic algorithm. The product of the maximum infiltration rate and the time step of the two-dimensional hydrodynamic algorithm is taken as the maximum infiltration volume of the grid node in the next calculation cycle.

[0016] In one implementation, determining the maximum infiltration volume of a grid node in the next calculation cycle based on the current cumulative infiltration volume further includes:

[0017] If a specified maximum cumulative infiltration volume is configured, and the sum of the current cumulative infiltration volume and the maximum infiltration volume for the next calculation period determined based on the Horton method is greater than the specified maximum cumulative infiltration volume, then the difference between the specified maximum cumulative infiltration volume and the current cumulative infiltration volume will be used as the maximum infiltration volume for the next calculation period.

[0018] In one implementation, determining the actual infiltration volume of a grid node in the next calculation cycle based on the maximum infiltration volume and the water source containing infiltration includes:

[0019] The maximum infiltration volume in the next calculation cycle and the minimum value among the water sources containing infiltration are used as the actual infiltration volume of the grid nodes in the next calculation cycle.

[0020] In one implementation, determining the cumulative infiltration volume of a grid node after the next calculation cycle based on the actual infiltration volume in the next calculation cycle and the current cumulative infiltration volume includes:

[0021] If the actual infiltration volume in the next calculation cycle is not zero, the sum of the current cumulative infiltration volume and the actual infiltration volume will be used as the cumulative infiltration volume of the grid node in the next calculation cycle.

[0022] In one implementation, determining the cumulative infiltration volume of the grid node after the next calculation cycle based on the actual infiltration volume of the next calculation cycle and the current cumulative infiltration volume further includes:

[0023] If the actual infiltration volume in the next calculation cycle is 0, perform the following steps:

[0024] Based on the Horton method, the time parameter in the Horton method is calculated according to the current cumulative infiltration volume. Then, the new time parameter in the Horton method is calculated according to the infiltration capacity recovery rate. Based on the new time parameter in the Horton method, the cumulative infiltration volume of the grid node after the next calculation cycle is determined.

[0025] Alternatively, if the actual infiltration volume in the next calculation cycle is 0, perform the following steps:

[0026] Based on the improved Horton method, the cumulative infiltration volume of the grid nodes after the next calculation cycle is determined according to the current cumulative infiltration volume and the infiltration capacity recovery rate.

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

[0028] The external water source conditions of the surface two-dimensional hydrodynamic algorithm are modified according to the actual infiltration volume in the next calculation cycle.

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

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

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

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

[0033] The actual infiltration volume determination module is used to determine the maximum infiltration volume of the grid node in the next calculation cycle based on the Horton method, when the water source containing infiltration is not zero, and to determine the actual infiltration volume of the grid node in the next calculation cycle based on the current cumulative infiltration volume, and the maximum infiltration volume and the water source containing infiltration. When the water source containing infiltration is zero, the actual infiltration volume of the grid node in the next calculation cycle is determined to be zero.

[0034] The cumulative infiltration volume update module is used to determine the cumulative infiltration volume of grid nodes after the next calculation cycle based on the Horton method, the actual infiltration volume of the next calculation cycle, the current cumulative infiltration volume, and the infiltration capacity recovery rate.

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

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

[0037] The present invention provides a coupled simulation method and apparatus for two-dimensional surface hydrodynamics and infiltration based on the Horton 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. If the water source containing infiltration is not zero, based on the Horton method, the water source containing infiltration for the grid nodes in the next calculation cycle is determined according to the current cumulative infiltration volume. The maximum infiltration volume for each calculation cycle is used to determine the actual infiltration volume of the grid node in the next calculation cycle based on the maximum infiltration volume and the water source containing infiltration. If the water source containing infiltration is 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be 0. Then, based on the actual infiltration volume, the current cumulative infiltration volume, and the infiltration capacity recovery rate, the cumulative infiltration volume of the grid node after the next calculation cycle is determined. Finally, 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.

[0038] 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, making the algorithm more realistic in simulating surface water flow evolution and significantly enhancing its application value in related business applications.

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

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

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

[0042] Figure 1A flowchart illustrating a coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Horton method, provided for an embodiment of the present invention;

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

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

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

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

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

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

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

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

[0051] Step S104: Based on the known current surface water volume of the grid node and the surface water volume of the grid node obtained from the pre-executed two-dimensional hydrodynamic calculation, determine the water source including infiltration for the grid node in the next calculation cycle.

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

[0053] Step S106: If the amount of water source containing infiltration is not zero, based on the Horton method, determine the maximum infiltration amount of the grid node in the next calculation cycle according to the current cumulative infiltration amount, and determine the actual infiltration amount of the grid node in the next calculation cycle according to the maximum infiltration amount and the water source containing infiltration; if the amount of water source containing infiltration is zero, determine that the actual infiltration amount of the grid node in the next calculation cycle is zero.

[0054] In one example, when the water source containing infiltration is not zero, the maximum infiltration volume of the grid node in the next calculation cycle can be determined based on the Horton method or a modified Horton method, according to the current cumulative infiltration volume. Then, the minimum value between the maximum infiltration volume and the water source containing infiltration is taken as the actual infiltration volume of the grid node in the next calculation cycle. When the water source containing infiltration is zero, the actual infiltration volume of the grid node in the next calculation cycle can be directly determined to be zero.

[0055] Step S108: Determine the cumulative infiltration volume of the grid node after the next calculation cycle based on the actual infiltration volume, the current cumulative infiltration volume, and the infiltration capacity recovery rate.

[0056] In one example, if the actual infiltration volume is not zero, the sum of the current cumulative infiltration volume and the actual infiltration volume can be directly used as the cumulative infiltration volume after the next calculation cycle. If the actual infiltration volume is zero, it is necessary to determine the cumulative infiltration volume after the next calculation cycle based on the current cumulative infiltration volume and the infiltration capacity recovery rate of the grid node, using the Horton method or an improved Horton method.

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

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

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

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

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

[0062] (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 It includes four scenarios: Scenario A, Scenario B, Scenario C, and Scenario D. Specifically:

[0063] Scenario A: If the surface water volume of the grid node obtained from the pre-executed two-dimensional hydrodynamic calculation is greater than the known current surface water volume of the grid node, then the water source including infiltration in the next calculation cycle will be the surface water volume of the grid node obtained from the pre-executed two-dimensional hydrodynamic calculation, that is:

[0064] ;

[0065] In the formula, For water sources that include infiltration, 'The surface water volume of the grid nodes is obtained from the pre-executed two-dimensional hydrodynamic calculation of the surface.' The known current surface water volume for each grid node.

[0066] Scenario B: If the known current surface water volume of a grid node is not zero, and the surface water volume of the grid node obtained from the pre-executed two-dimensional hydrodynamic calculation is less than the known current surface water volume of the grid node, then the water source including infiltration in the next calculation cycle will be the arithmetic mean of 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, that is:

[0067] ;

[0068] In the formula, For water sources that include infiltration, 'The surface water volume of the grid nodes is obtained from the pre-executed two-dimensional hydrodynamic calculation of the surface.' The known current surface water volume for each grid node.

[0069] Scenario C: If the known current surface water volume of a grid node is not zero, and the surface water volume of the grid node obtained from the pre-executed two-dimensional hydrodynamic calculation is equal to the known current surface water volume of the grid node, then the water source including infiltration in the next calculation cycle will be the known current surface water volume of the grid node, that is:

[0070] ;

[0071] In the formula, For water sources that include infiltration, 'The surface water volume of the grid nodes is obtained from the pre-executed two-dimensional hydrodynamic calculation of the surface.' The known current surface water volume for each grid node.

[0072] In scenario D, if the known current surface water volume of a grid node is 0, and the surface water volume of the grid node obtained from the pre-executed two-dimensional hydrodynamic calculation is 0, then the water source including infiltration in the next calculation cycle will be 0, and the actual infiltration volume in the next calculation cycle will be 0.

[0073] (iv) If the amount of water source including infiltration is not zero, perform the following steps to determine the actual infiltration volume for the next calculation period, and then perform (vi), (vii), and (viii):

[0074] (4.1) Based on the Horton method, determine the maximum infiltration volume of the grid node in the next calculation cycle according to the current cumulative infiltration volume.

[0075] Method 1: The maximum infiltration volume of the grid node in the next calculation cycle can be determined based on the Horton method. The process is as follows: Based on the Horton method, the time parameter in the Horton method is calculated according to the current cumulative infiltration volume. The sum of the time parameter in the Horton method and the time step of the two-dimensional hydrodynamic algorithm is used as the new time parameter in the Horton method. Then, based on the Horton method, the maximum cumulative infiltration volume of the grid node after the next calculation cycle is determined according to the new time parameter in the Horton method.

[0076] Specifically, based on the Horton method:

[0077] ;

[0078] In the formula, This represents the cumulative infiltration volume in the Horton method. This represents the final infiltration rate in the Horton method. This is the time variable in the Horton method. This represents the initial infiltration rate in the Horton method. This is the infiltration capacity attenuation coefficient in the Horton method.

[0079] Substitute the current cumulative infiltration volume into the cumulative infiltration volume in the above formula. Calculation ,Should This is the time parameter in the Horton method. Add this time to the time step of the two-dimensional hydrodynamic algorithm for the Earth's surface, and then substitute it into the time variable in the Horton method of the above formula. The maximum cumulative infiltration volume after the next calculation cycle is calculated, and 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 of the grid node in the next calculation cycle.

[0080] Method 2: The maximum infiltration volume of the grid node in the next calculation cycle can be determined based on the improved Horton method. The process is as follows: Based on the improved Horton method, the maximum infiltration rate of the grid node in the next calculation cycle is determined according to the current cumulative infiltration volume and the time parameters of the surface two-dimensional hydrodynamic algorithm. The product of the maximum infiltration rate and the time step of the surface two-dimensional hydrodynamic algorithm is taken as the maximum infiltration volume of the grid node in the next calculation cycle.

[0081] Specifically, based on the improved Horton method:

[0082] ;

[0083] In the formula, This represents the maximum infiltration rate in the Horton method. This represents the final infiltration rate in the Horton method. This represents the initial infiltration rate in the Horton method. This is the infiltration capacity attenuation coefficient in the Horton method. This represents the cumulative infiltration volume in the Horton method. It is a time variable.

[0084] In practical implementation, the current cumulative infiltration volume is substituted into the cumulative infiltration volume variable in the above formula. Substitute the time parameters of the two-dimensional hydrodynamic algorithm for the surface into the time variable in the above formula. The maximum infiltration rate for the next calculation cycle is calculated, and the product of the maximum infiltration rate and the time step of the two-dimensional hydrodynamic algorithm is used as the maximum infiltration volume of the grid node in the next calculation cycle.

[0085] Furthermore, if a specified maximum cumulative infiltration volume is configured, and the sum of the current cumulative infiltration volume and the maximum infiltration volume for the next calculation period determined based on the Horton method or an improved Horton method is greater than the specified maximum cumulative infiltration volume, then the difference between the specified maximum cumulative infiltration volume and the current cumulative infiltration volume will be used as the maximum infiltration volume for the next calculation period. In specific implementation, if a specified maximum cumulative infiltration volume is configured, and the sum of the current cumulative infiltration volume and the maximum infiltration volume for the next calculation period determined based on the Horton method or an improved Horton method exceeds the specified maximum cumulative infiltration volume, then the maximum infiltration volume for the next calculation period will be modified to the specified maximum cumulative infiltration volume minus the current cumulative infiltration volume.

[0086] (4.2) Determine the actual infiltration volume of the grid node in the next calculation cycle based on the maximum infiltration volume and the water source containing infiltration. Specifically, the minimum value between the maximum infiltration volume and the water source containing infiltration in the next calculation cycle is taken as the actual infiltration volume of the grid node in the next calculation cycle.

[0087] (v) When the water source including infiltration is 0, determine that the actual infiltration volume of the grid node in the next calculation cycle is 0, and execute (vii) and (viii).

[0088] (vi) 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 then modified.

[0089] ;

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

[0091] (vii) Determine the cumulative infiltration volume of the grid nodes after the next calculation cycle based on the actual infiltration volume of the next calculation cycle and the current cumulative infiltration volume.

[0092] Scenario 1: If the actual infiltration volume in the next calculation cycle is not zero, the sum of the current cumulative infiltration volume and the actual infiltration volume will be used as the cumulative infiltration volume of the grid node in the next calculation cycle.

[0093] Scenario 2: If the actual infiltration volume in the next calculation cycle is 0, the groundwater drainage during the period without infiltration must be calculated, reflecting a reduction in the cumulative infiltration volume. Based on the Horton method or a modified Horton method, determine the cumulative infiltration volume of the grid nodes after the next calculation cycle.

[0094] Method 1: Determine the cumulative infiltration volume of the grid node after the next calculation cycle based on the Horton method. The specific process is as follows: Based on the Horton method, calculate the time parameter in the Horton method according to the current cumulative infiltration volume, and then calculate the new time parameter in the Horton method according to the infiltration capacity recovery rate. Based on the new time parameter in the Horton method, determine the cumulative infiltration volume of the grid node after the next calculation cycle.

[0095] Specifically, based on the Horton method:

[0096] ;

[0097] In the formula, This is a new time variable in the Horton method. This is the infiltration capacity attenuation coefficient in the Horton method. This is the infiltration recovery coefficient in the Horton method. The time step of the two-dimensional hydrodynamic algorithm for the Earth's surface. This is the time variable in the Horton method.

[0098] In practical implementation, the current cumulative infiltration volume is substituted into the cumulative infiltration volume variable in the Horton method of the formula shown in Method 1 of (4.1). The time parameters in the Horton method are calculated, and then these time parameters are substituted into the time variables in the Horton method in the above formula. The new time parameters in the Horton method are calculated, and then the new time parameters in the Horton method are substituted into the time variables in the Horton method in the formula shown in Method 1 of (4.1). The cumulative infiltration volume after the next calculation cycle is calculated.

[0099] Method 2: Based on the improved Horton method, the cumulative infiltration volume of the grid nodes after the next calculation cycle is determined. The specific process is as follows: Based on the improved Horton method, the cumulative infiltration volume of the grid nodes after the next calculation cycle is determined according to the current cumulative infiltration volume and the infiltration capacity recovery rate.

[0100] Specifically, based on the improved Horton method:

[0101] ;

[0102] In the formula, This is the current cumulative infiltration volume for the next calculation period (also referred to as the corrected cumulative infiltration volume). This represents the current cumulative infiltration volume. This represents the final infiltration rate in the Horton method. For time variables, This is the infiltration recovery coefficient in the Horton method. This represents the time step of the two-dimensional hydrodynamic algorithm for the Earth's surface.

[0103] In practice, the current cumulative infiltration volume, the time parameters of the two-dimensional hydrodynamic algorithm, and the time step of the two-dimensional hydrodynamic algorithm are substituted into the above formula to calculate the cumulative infiltration volume after the next calculation cycle.

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

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

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

[0107] Repeatedly executing the above algorithm process constitutes the complete coupled calculation process of surface two-dimensional hydrodynamics and Horton method infiltration.

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

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

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

[0111] 3. In the Horton 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 the Horton method to the coupling of the surface two-dimensional hydrodynamic algorithm.

[0112] The algorithm provided in this invention can calculate the actual infiltration volume at every moment and at all locations, and correct the surface water volume calculated using two-dimensional hydrodynamics 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, while also addressing the significant impact of infiltration on surface water volume in reality. The algorithm provided in 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.

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

[0114] The pre-execution module 302 is used to pre-execute the two-dimensional hydrodynamic calculation of the surface in the next calculation cycle without considering infiltration;

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

[0116] The actual infiltration volume determination module 306 is used to determine the maximum infiltration volume of the grid node in the next calculation cycle based on the Horton method and the current cumulative infiltration volume when the water source containing infiltration is not zero, and to determine the actual infiltration volume of the grid node in the next calculation cycle based on the maximum infiltration volume and the water source containing infiltration; when the water source containing infiltration is zero, the actual infiltration volume of the grid node in the next calculation cycle is determined to be zero.

[0117] The cumulative infiltration volume update module 308 is used to determine the cumulative infiltration volume of the grid node after the next calculation cycle based on the actual infiltration volume of the next calculation cycle, the current cumulative infiltration volume, and the infiltration capacity recovery rate.

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

[0119] The coupled simulation device of two-dimensional surface hydrodynamics and infiltration based on the Horton method provided in this invention, compared with existing two-dimensional surface hydrodynamic algorithms, further improves the ability to simulate the dynamic processes of surface water conservation and inertial convection, propagation, and diffusion under the influence of gravity and surface resistance. This is based on the existing ability to simulate the dynamic processes of surface water flow under gravity and surface resistance. The algorithm provided by this invention makes existing 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.

[0120] In one implementation, the actual infiltration volume determination module 306 is specifically used for:

[0121] Based on the Horton method, the time parameter in the Horton method is calculated according to the current cumulative infiltration volume. The sum of the time parameter in the Horton method and the time step of the two-dimensional hydrodynamic algorithm is used as the new time parameter in the Horton method. Based on the new time parameter in the Horton method, the maximum cumulative infiltration volume of the grid node after the next calculation cycle is determined.

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

[0123] In one implementation, the actual infiltration volume determination module 306 is specifically used for:

[0124] Based on the improved Horton method, the maximum infiltration rate of the grid node in the next calculation cycle is determined according to the current cumulative infiltration volume and the time parameters of the two-dimensional hydrodynamic algorithm. The product of the maximum infiltration rate and the time step of the two-dimensional hydrodynamic algorithm is taken as the maximum infiltration volume of the grid node in the next calculation cycle.

[0125] In one implementation, the actual infiltration volume determination module 306 is specifically used for:

[0126] If a specified maximum cumulative infiltration volume is configured, and the sum of the current cumulative infiltration volume and the maximum infiltration volume determined based on the Horton method is greater than the specified maximum cumulative infiltration volume, then the difference between the specified maximum cumulative infiltration volume and the current cumulative infiltration volume will be used as the maximum infiltration volume for the next calculation cycle.

[0127] In one implementation, the actual infiltration volume determination module 306 is specifically used for:

[0128] The maximum infiltration volume and the minimum value among the infiltrating water sources are used as the actual infiltration volume of the grid node in the next calculation cycle.

[0129] In one implementation, the cumulative infiltration volume update module 308 is specifically used for:

[0130] If the actual infiltration volume is not zero, 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 of the grid node after the next calculation cycle.

[0131] In one implementation, the cumulative infiltration volume update module 308 is specifically used for:

[0132] If the actual infiltration rate is 0, perform the following steps:

[0133] Based on the Horton method, the time parameter in the Horton method is calculated according to the current cumulative infiltration volume. Then, the new time parameter in the Horton method is calculated according to the infiltration capacity recovery rate. Based on the new time parameter in the Horton method, the cumulative infiltration volume of the grid node after the next calculation cycle is determined.

[0134] Alternatively, if the actual infiltration rate is 0, perform the following steps:

[0135] Based on the improved Horton method, the cumulative infiltration volume of the grid nodes after the next calculation cycle is determined according to the current cumulative infiltration volume and the infiltration capacity recovery rate.

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

[0137] The external water source conditions of the surface two-dimensional hydrodynamic algorithm are modified according to the actual infiltration volume in the next calculation cycle.

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

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

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

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

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

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

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

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

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

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

[0148] 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 Horton method, characterized in that, include: Without considering infiltration, perform the two-dimensional hydrodynamic calculation of the surface in advance for the next calculation cycle; Based on the known current surface water volume of the grid node and the surface water volume of the grid node obtained from the two-dimensional hydrodynamic calculation of the surface in the previously executed calculation cycle, the water source including infiltration of the grid node in the next calculation cycle is determined. If the water source containing infiltration is not zero, based on the Horton method, the maximum infiltration volume of the grid node in the next calculation cycle is determined according to the current cumulative infiltration volume, and the actual infiltration volume of the grid node in the next calculation cycle is determined according to the maximum infiltration volume and the water source containing infiltration. If the water source containing infiltration is 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be 0. Based on the actual infiltration volume, the current cumulative infiltration volume, and the infiltration capacity recovery rate, the cumulative infiltration volume of the grid node after the next calculation cycle is determined. Based on the actual infiltration volume, the external water source conditions for the two-dimensional hydrodynamic calculation of the surface are modified, and combined with the modified external water source conditions, the two-dimensional hydrodynamic calculation of the surface is corrected for the next calculation cycle to obtain the surface water volume after the next calculation cycle. The modification of the external water source conditions for the two-dimensional hydrodynamic calculation of the Earth's surface is achieved using the following formula: （1） In the formula, For the modified external water source conditions of 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.

2. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Horton method according to claim 1, characterized in that, Based on the Horton method, the maximum infiltration volume of the grid node in the next calculation cycle is determined according to the current cumulative infiltration volume, including: Based on the Horton method, the time parameter in the Horton method is calculated according to the current cumulative infiltration volume. The sum of the time parameter in the Horton method and the time step of the two-dimensional hydrodynamic calculation of the surface is used as the new time parameter in the Horton method. Based on the new time parameter in the Horton method, the maximum cumulative infiltration volume of the grid node after the next calculation cycle is determined. 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 of the grid node in the next calculation cycle.

3. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Horton method according to claim 1, characterized in that, Based on the Horton method, the maximum infiltration volume of the grid node in the next calculation cycle is determined according to the current cumulative infiltration volume, and the method further includes: Based on the improved Horton method, the maximum infiltration rate of the grid node in the next calculation cycle is determined according to the current cumulative infiltration volume and the time parameters of the two-dimensional hydrodynamic calculation of the surface. The product of the maximum infiltration rate and the time step of the two-dimensional hydrodynamic calculation of the surface is taken as the maximum infiltration volume of the grid node in the next calculation cycle.

4. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Horton method according to claim 3, characterized in that, Determining the maximum infiltration volume of the grid node in the calculation cycle based on the current cumulative infiltration volume also includes: If a specified maximum cumulative infiltration volume is configured, and the sum of the current cumulative infiltration volume and the maximum infiltration volume determined based on the Horton method is greater than the specified maximum cumulative infiltration volume, then the difference between the specified maximum cumulative infiltration volume and the current cumulative infiltration volume is used as the maximum infiltration volume for the next calculation cycle.

5. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Horton method according to claim 1, characterized in that, Determining the actual infiltration volume of the grid node in the next calculation cycle based on the maximum infiltration volume and the water source containing infiltration includes: The minimum value among the maximum infiltration volume and the infiltrated water source is taken as the actual infiltration volume of the grid node in the next calculation cycle.

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

7. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Horton method according to claim 1, characterized in that, The method for determining the cumulative infiltration volume of a grid node after the next calculation cycle, based on the actual infiltration volume, the current cumulative infiltration volume, and the infiltration capacity recovery rate, further includes: If the actual infiltration volume is 0, perform the following steps: Based on the Horton method, the time parameter in the Horton method is calculated according to the current cumulative infiltration volume, and then the new time parameter in the Horton method is calculated according to the infiltration capacity recovery rate. Based on the new time parameter in the Horton method, the cumulative infiltration volume of the grid node after the next calculation cycle is determined. Alternatively, if the actual infiltration volume is 0, perform the following steps: Based on the improved Horton method, the cumulative infiltration volume of the grid node after the next calculation cycle is determined according to the current cumulative infiltration volume and the infiltration capacity recovery rate.

8. The coupled simulation method of two-dimensional surface hydrodynamics and infiltration based on the Horton method according to claim 1, characterized in that, Based on the actual infiltration volume and the current cumulative infiltration volume, the current cumulative infiltration volume of the grid node in the next calculation cycle is determined, including: If the actual infiltration volume is not zero, the sum of the current cumulative infiltration volume and the actual infiltration volume is taken as the cumulative infiltration volume of the grid node in the next calculation cycle.

9. A coupled simulation device for two-dimensional surface hydrodynamics and infiltration based on the Horton method, characterized in that, include: The pre-execution module is used to pre-execute the two-dimensional hydrodynamic calculations of the surface for the next calculation cycle without considering infiltration. The water source 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 pre-executed two-dimensional hydrodynamic calculation. The actual infiltration volume determination module is used to determine the maximum infiltration volume of the grid node in the next calculation cycle based on the Horton method and the current cumulative infiltration volume when the water source containing infiltration is not zero, and to determine the actual infiltration volume of the grid node in the next calculation cycle based on the maximum infiltration volume and the water source containing infiltration. If the water source containing infiltration is 0, the actual infiltration volume of the grid node in the next calculation cycle is determined to be 0. The cumulative infiltration volume update module is used to determine the cumulative infiltration volume of the grid node after the next calculation cycle based on the actual infiltration volume, the current cumulative infiltration volume, and the infiltration capacity recovery rate. Based on the actual infiltration volume, the external water source conditions for the two-dimensional hydrodynamic calculation of the surface are modified, and the execution module is corrected in combination with the modified external water source conditions. This is used to correct the execution of 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 accumulation after the next calculation cycle. The modification of the external water source conditions for the two-dimensional hydrodynamic calculation of the Earth's surface is achieved using the following formula: （1） In the formula, For the modified external water source conditions of 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.

10. 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 8.

Images