A simple method for creating a groundwater cell infiltration matrix

Abstract

The application discloses a simple method for creating a groundwater unit infiltration matrix, comprising the following steps: calculating horizontal seepage flow and vertical seepage flow between nodes based on Darcy formula; calculating gravity water release amount and elastic water release amount of the nodes; and substituting the coefficients of groundwater head in the horizontal seepage flow, the vertical seepage flow, the gravity water release amount and the elastic water release amount into corresponding positions in a pre-constructed groundwater unit infiltration matrix to obtain a final groundwater unit infiltration matrix. According to the Darcy formula, the groundwater unit infiltration matrix can be simply constructed, and complex mathematical derivation is no longer needed, that is, the method adopts a more intuitive physical law instead of abstract mathematical transformation, and the design and development difficulty of a groundwater numerical simulation model is greatly reduced.

Description

Technical Field

[0001] This invention relates to a simple method for creating a groundwater unit permeability matrix, belonging to the technical field of groundwater numerical simulation models. Background Technology

[0002] Numerical simulation of groundwater plays a crucial role in groundwater engineering problems such as groundwater resource management, groundwater pollution prevention and remediation, and is an important method for solving scientific and engineering problems related to groundwater. The construction of the permeability matrix of groundwater units is the core of the design and development of groundwater numerical simulation models. Current methods first construct basis functions, then use these basis functions to transform the differential equations of groundwater flow motion into integral equations, thereby obtaining the permeability matrix of the groundwater units. The basis function construction process involves complex mathematical operations, such as integration and polynomial interpolation, resulting in high computational complexity. This is especially true when dealing with heterogeneous aquifers or multidimensional, multi-field problems, where the derivation of the integral equations becomes significantly more difficult. The mathematical derivation in this process is complex, requiring high mathematical ability and is operationally challenging. Summary of the Invention

[0003] The technical problem to be solved by this invention is that the construction of groundwater unit permeability matrix is ​​complex or even difficult. In order to overcome the defects of the prior art, a simple method for creating groundwater unit permeability matrix is ​​provided. Based on Darcy's formula and basic hydrogeological parameters, the design and development of groundwater numerical simulation model is made simpler, the logical steps are simplified, and the efficiency of building groundwater numerical model is improved.

[0004] Prior to this invention, a simplified method for creating a groundwater unit permeability matrix is ​​provided, comprising: calculating the horizontal seepage flow between nodes based on Darcy's formula; calculating the vertical seepage flow from one node to the next layer node based on Darcy's formula; calculating the gravity release of the nodes; calculating the elastic release of the nodes; adding the coefficients of the same groundwater head among the horizontal seepage flow, vertical seepage flow, gravity release, and elastic release respectively to obtain the coefficient of the final groundwater head; and substituting the coefficient of the final groundwater head into the corresponding position in the pre-constructed groundwater unit permeability matrix to obtain the final groundwater unit permeability matrix.

[0005] Prioritizes calculating the horizontal seepage flow between nodes based on Darcy's formula, including: calculating the horizontal seepage flow from node ② into node ①. : Equation (1), where, It is the average hydraulic conductivity between node ② and node ①; It is the distance from the circumcenter of the triangle to the sides of the triangle containing nodes ① and ②. It is the length of the line segment between node ② and node ①; It is the water head value of node ② at time n+1. It is the head value of node ① at time n+1; The time period is long; calculate the horizontal seepage flow from node ③ into node ①. .

[0006] Equation (2).

[0007] In the formula, It is the average hydraulic conductivity between node ③ and node ①; It is the distance from the circumcenter of the triangle to the sides of the triangle containing nodes ① and ③. It is the length of the line segment between node ③ and node ①; It is the head value of node ③ at time n+1.

[0008] Calculate the horizontal seepage flow from node ③ into node ②. : Equation (3).

[0009] In the formula, It is the average hydraulic conductivity between node ③ and node ②; It is the distance from the circumcenter of the triangle to the sides of the triangle containing nodes ② and ③. It is the length of the line segment between node ③ and node ②.

[0010] The coefficient of groundwater head includes the coefficient generated by formula (1). = The coefficients generated by formula (1) = The coefficients generated by formula (2) = The coefficients generated by formula (2) = The coefficients generated by formula (3) = The coefficients generated by formula (3) = .

[0011] Prioritize the horizontal seepage flow from node ① into node ② as - The horizontal seepage flow from node ① into node ③ is - The horizontal seepage flow from node ② into node ③ is - .

[0012] Firstly, calculate the vertical seepage flow from node ① into the next layer node, including: assuming node ④ is the next layer node corresponding to node ①, node ⑤ is the next layer node corresponding to node ②, and node ⑥ is the next layer node corresponding to node ③, calculate the vertical seepage flow from node ④ into node ①. Vertical seepage flow from node ⑤ into node ② Vertical seepage flow from node 6 into node 3 .

[0013] Equation (4) Equation (5) In equation (6), The vertical permeability coefficient between node ① and node ④ is given. The vertical permeability coefficient between node ② and node ⑤ is given. The vertical permeability coefficient between node ③ and node ⑥; Let be the head value of node ④ at time n+1. Let be the head value of node ⑤ at time n+1. Let the head value of node ⑥ be at time n+1; For the festival ① and nodes The vertical distance between them For nodes and nodes The vertical distance between them For nodes and nodes The vertical distance between them; The area of ​​the sub-equilibrium region at node ①. The area of ​​the sub-equilibrium region of node ②. Let be the area of ​​the sub-equilibrium region of node ③.

[0014] The coefficient of groundwater head includes the coefficient generated by formula (4). = The coefficients generated by formula (4) = The coefficients generated by formula (5) = The coefficients generated by formula (5) = The coefficients generated by formula (6) = The coefficients generated by formula (6) = .

[0015] Prioritize the vertical seepage flow from node ① into node ④ as follows: The vertical seepage flow from node ② into node ⑤ - The vertical seepage flow from node ③ into node ⑥ - .

[0016] Prioritize calculating the gravity-induced water release at the node, including: calculating the gravity-induced water release at node ①. Gravity-induced water release at node ② and the amount of water released by gravity at node ③ : Equation (7) Equation (8) Equation (9), where, It refers to the water supply degree; It is the water head value at time n of node ①. It is the water head value at time n of node ②. This is the water head value at time n of node ③. The coefficient of groundwater head includes the coefficient generated by formula (7). = The coefficients generated by formula (8) = The coefficients generated by formula (9) = .

[0017] Prioritize calculating the elastic water release of the node, including: when there is confined water, calculating the elastic water release of node ①. Elastic water release at node ② And the elastic water release of node ③ : Equation (10) Equation (11). Equation (12), where, It is the elastic water release coefficient; the coefficient of groundwater head includes the generation coefficient of formula (10). = The coefficients generated by formula (11) = Formula (12) generates coefficients = .

[0018] Prioritize adding the coefficients of the same groundwater head among horizontal seepage flow, vertical seepage flow, gravity release, and elastic release to obtain the final groundwater head coefficient, including the coefficients generated by formulas (1), (2), (4), and (7). Add the values ​​together, and use the coefficients generated by formula (1) Add the values ​​together, and use the coefficients generated by formula (2) Add the values ​​together, and use the coefficients generated by formula (4) Add them together to obtain the coefficients generated by formulas (3), (5), and (8). Add the values ​​together, and use the coefficients generated by formula (3) Add the values ​​together, and use the coefficients generated by formula (5) Add the values ​​together, and use the coefficients generated by formula (6) Add the values ​​together, and combine the coefficients generated by formulas (6) and (9). Adding the values ​​together yields the following: , , , , , , , , The coefficient of the final groundwater head.

[0019] Preferably, the present invention provides an electronic device including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps of any of the methods described herein.

[0020] Preferably, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of any of the methods described herein.

[0021] The beneficial effects achieved by this invention are as follows: Existing methods, when constructing the permeability matrix of a groundwater unit, first require the construction of basis functions, and then the transformation of the differential equation of groundwater flow motion into an integral equation using basis functions, residual weighting, or variational methods to obtain the permeability matrix of the groundwater unit. This process involves complex mathematical derivations, requires high mathematical ability, and is difficult to operate. Existing methods also present complex challenges in handling some special problems, such as dry-wet transitions and confined-unconfined transitions.

[0022] The method of this invention takes a different approach. Instead of relying on basis functions and differential equations of groundwater flow, it can simply construct the permeability matrix of groundwater units based mainly on Darcy's formula. It also eliminates the need for complex mathematical derivations. In other words, the method adopts more intuitive physical laws rather than abstract mathematical transformations, which greatly reduces the difficulty of designing and developing groundwater numerical simulation models. In particular, it can be quickly mastered by non-groundwater professionals, enabling them to easily complete model design and development. Attached Figure Description

[0023] To more clearly illustrate the technical solution of this application, the drawings used in the embodiments will be briefly introduced below. Obviously, for those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0024] Figure 1 This is a schematic diagram illustrating the calculation of horizontal seepage flow between nodes in some embodiments of this application.

[0025] Figure 2 This is a structural diagram of the permeability matrix of the groundwater unit in some embodiments of this application.

[0026] Figure 3 These are simulated partition diagrams in some embodiments of this application.

[0027] Figure 4 A comparison of the water level duration curves of observation well O-1 for groundwater flow numerical simulation using FEFLOW and GFModel.

[0028] Figure 5 A comparison of the water level duration curves of observation well O-2 for groundwater flow numerical simulation using FEFLOW and GFModel.

[0029] Figure 6 A comparison of the water level duration curves of observation well O-3 for groundwater flow numerical simulation using FEFLOW and GFModel.

[0030] Figure 7 A comparison of the water level duration curves of observation well O-4 for groundwater flow numerical simulation using FEFLOW and GFModel.

[0031] Figure 8 This is a water level map calculated by the method of the present invention. Detailed Implementation

[0032] Example 1

[0033] See Figure 1 This application discloses a simple method for creating a groundwater unit permeability matrix, which uses triangles to divide the study area. Figure 1 This diagram illustrates the calculation of horizontal seepage flow between nodes within a triangle. The three nodes of the triangle are numbered ①, ②, and ③. It is the circumcenter of the triangle and also the node in the equilibrium region corresponding to the triangle. q is the midpoint of the side of the triangle containing nodes ① and ②, q is the midpoint of the side of the triangle containing nodes ① and ③, and m is the midpoint of the side of the triangle containing nodes ② and ③.

[0034] The specific steps and description of the method of the present invention are as follows: (1) Calculate the horizontal seepage flow between nodes using Darcy's formula. Specifically, calculate the horizontal seepage flow from node ② into node ①. : Equation (1), where, It is the average hydraulic conductivity between node ② and node ①; It is the distance from the circumcenter of the triangle to the sides of the triangle containing nodes ① and ②. It is the length of the line segment between node ② and node ①; It is the water head value of node ② at time n+1. It is the head value of node ① at time n+1; The time period is long. The horizontal seepage flow from node ① into node ② is - .

[0035] Calculate the horizontal seepage flow from node ③ into node ①. : Equation (2), where, It is the average hydraulic conductivity between node ③ and node ①; It is the distance from the circumcenter of the triangle to the sides of the triangle containing nodes ① and ③. It is the length of the line segment between node ③ and node ①; This is the head value of node ③ at time n+1. The horizontal seepage flow from node ① into node ③ is - .

[0036] Calculate the horizontal seepage flow from node ③ into node ②. : Equation (3), where, It is the average hydraulic conductivity between node ③ and node ②; It is the distance from the circumcenter of the triangle to the sides of the triangle containing nodes ② and ③. This is the length of the line segment between node ③ and node ②. The horizontal seepage flow from node ② into node ③ is - .

[0037] Calculate the vertical seepage flow from node ① into the next level node: When establishing a three-dimensional groundwater flow numerical model, assume that the next level node corresponding to node ① is node ④, the next level node corresponding to node ② is node ⑤, and the next level node corresponding to node ③ is node ⑥. Calculate the vertical seepage flow from node ④ into node ①. Vertical seepage flow from node ⑤ into node ② Vertical seepage flow from node 6 into node 3 : Equation (4) Equation (5) Equation (6), where, The vertical permeability coefficient between node ① and node ④ is given. The vertical permeability coefficient between node ② and node ⑤ is given. The vertical permeability coefficient between node ③ and node ⑥; Let be the head value of node ④ at time n+1. Let be the head value of node ⑤ at time n+1. Let the head value of node ⑥ be at time n+1; For the festival ① and nodes The vertical distance between them For nodes and nodes The vertical distance between them For nodes and nodes The vertical distance between them; The area of ​​the sub-equilibrium region at node ①. The area of ​​the sub-equilibrium region of node ②. Let be the area of ​​the sub-equilibrium region of node ③. Specifically... For node ①, The area of ​​the enclosed quadrilateral. For node ②, The area of ​​the enclosed quadrilateral. For node ③, The area of ​​the enclosed quadrilateral. The vertical seepage flow from node ① into node ④ is - The vertical seepage flow from node ② into node ⑤ - The vertical seepage flow from node ③ into node ⑥ - .

[0038] (2) The structural design of the permeability matrix of the groundwater unit in the three-dimensional groundwater flow numerical model is as follows: Figure 2 As shown. Figure 2 In and The values ​​in the permeability matrix of a groundwater unit are generally non-zero. Figure 2 The spaces in the text represent 0. These values ​​are all formed by the current triangle. This represents the coefficient of the head at node ① in the equation corresponding to node ①. This represents the coefficient of the head at node 2 in the equation corresponding to node ①. This represents the coefficient of the head at node ① in the equation corresponding to node ②, and so on. This represents the coefficient of the head at node ③ in the equation corresponding to node ①. This represents the coefficient of the head at node ④ in the equation corresponding to node ①. This represents the coefficient of the head at node ② in the equation corresponding to node ②. This represents the coefficient of the head at node ③ in the equation corresponding to node ②. This represents the coefficient of the head at node ⑤ in the equation corresponding to node ②. This represents the coefficient of the head at node ① in the equation corresponding to node ③. This represents the coefficient of the head at node ② in the equation corresponding to node ③. This represents the coefficient of the head at node ③ in the equation corresponding to node ③. This represents the head coefficient of node ⑥ in the equation corresponding to node ③. This represents the coefficient of the head at node ① in the equation corresponding to node ④. This represents the coefficient of the head at node ④ in the equation corresponding to node ④. This represents the coefficient of the head at node ② in the equation corresponding to node ⑤. This represents the coefficient of the head at node ⑤ in the equation corresponding to node ⑤. This represents the coefficient of the head at node ③ in the equation corresponding to node ⑥. This represents the coefficient of the head at node ⑥ in the equation corresponding to node ⑥.

[0039] Put all the coefficients in formulas (1) to (6) into Figure 2 Among the various values, formula (1) produces... Value Formula (1) generates Value Formula (2) generates Value Formula (2) generates Value Formula (3) generates Value Formula (3) generates Value Formula (4) produces Value Formula (4) generates Value Formula (5) generates Value Formula (5) generates Value Formula (6) generates Value Formula (6) generates Value .

[0040] The results above show that, Figure 2 middle, =-( + + ),and = , = , = , = , = . , , , , All have The property that all elements in a given row are equal to the negative of the sum of the elements in that row. The permeability matrix of a groundwater unit exhibits significant symmetry and is diagonally dominant. When performing two-dimensional numerical simulations of groundwater flow, Figure 2 The number of elements in is 9, that is Figure 2 The 3 in the top left corner 3-matrix.

[0041] (3) The amount of water released also contributes to the permeability matrix of the groundwater unit. For unconfined groundwater, the gravity release at node ① Gravity-induced water release at node ② and the amount of water released by gravity at node ③ They are respectively: Equation (7) Equation (8) Equation (9), where, It refers to the water supply degree; It is the water head value at time n of node ①. It is the water head value at time n of node ②. It is the head value at time n of node ③; generated by formula (7) for Formula (8) generates for Formula (9) generates The value is .

[0042] When the water is under pressure, the elastic water release at node ① Elastic water release at node ② And the elastic water release of node ③ They are respectively: Equation (10) Equation (11). Equation (12), where, It is the elastic water release coefficient; generated by formula (10) The value is Formula (11) generates The value is Formula (12) generates The value is .

[0043] The above formulas (1) to (12) generate Figure 2 Each of the Sum the corresponding values ​​and then put them into... Figure 2 In this process, a groundwater unit permeability matrix is ​​formed.

[0044] As can be seen from the above derivation process, the method of the present invention does not use complex mathematical expressions such as basis functions, differential equations and integral equations, but only Darcy's formula to complete the construction of the permeability matrix of the groundwater unit. The method is very simple, which greatly reduces the process of constructing the permeability matrix of the groundwater unit and improves efficiency.

[0045] Example 2

[0046] The groundwater unit permeability matrix construction method proposed in this invention has developed a groundwater flow numerical model and compared it with the world-renowned groundwater numerical simulation software FEFLOW to demonstrate the effectiveness of the method.

[0047] The simulation area is rectangular, the aquifer is homogeneous and isotropic, and the aquifer floor is used as the reference surface. Calculations are performed over 30 days per month. The planar boundary conditions are set as impermeable boundaries. Table 1 shows the specific parameters for the groundwater flow numerical model.

[0048] Table 1. Specific parameters for setting up the numerical model of groundwater flow

[0049] The simulation region is partitioned using triangular elements, with 7504 elements and 3850 nodes per layer. The partitioning result is as follows: Figure 3 As shown in Table 2. In addition, four observation wells were set up in the groundwater flow numerical model, and the specific information of the observation wells is shown in Table 2.

[0050] Table 2. Observation Well Information

[0051] Observation well number X-coordinate (m) Y coordinate (m) O-1 200 500 O-2 700 500 O-3 1200 500 O-4 1000 750

[0052] Numerical simulations of groundwater flow were performed using FEFLOW and GFModel respectively. The results of the water level duration curves are shown below. Figures 4-7 As shown, the horizontal axis is in days, and the vertical axis is in meters. The groundwater level contour lines at the end of the second year are shown below. Figure 8 As shown, the units for both the horizontal and vertical axes are meters. Table 3 shows the comparison results of water balance simulations obtained by FEFLOW and the method of this invention.

[0053] Table 3. Comparison of water balance results between FEFLOW and the method of the present invention.

[0054] time <![CDATA[Change in storage volume in FEFLOW (m 3 )]]> <![CDATA[The change in storage volume (m 3 )]]> Relative error (%) First year 156197.00 158192.50 1.26 The second year 513399.00 518341.93 0.95 Third year 873364.00 878394.67 0.57

[0055] Table 3 shows that the method of the present invention is closest to the simulation results of FEFLOW, which also demonstrates that the method of the present invention is effective.

[0056] In this embodiment of the application, the present invention provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps of any of the methods described above.

[0057] In this application embodiment, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of any of the methods described above.

[0058] The various embodiments in this specification are described in a progressive manner. The same or similar parts between the various embodiments can be referred to each other. Each embodiment focuses on describing the differences from other embodiments.

[0059] Other embodiments of the invention will readily occur to those skilled in the art upon consideration of the specification and practice of the invention described herein. This application is intended to cover any variations, uses, or adaptations of the invention that follow the general principles of the invention and include common knowledge or customary techniques in the art not invented herein. The specification and embodiments are to be considered exemplary only.

[0060] The above specific embodiments further illustrate the purpose, technical solution and beneficial effects of this application. It should be understood that the above are only specific embodiments of this application and are not intended to limit the scope of protection of this application. Any modifications, equivalent substitutions, improvements, etc., made on the basis of the technical solution of this application should be included within the scope of protection of this application.

Claims

1. A simplified method for creating a groundwater unit permeability matrix, characterized in that, include: Based on Darcy's formula, calculate the horizontal seepage flow between nodes; Based on Darcy's formula, calculate the vertical seepage flow from one node to the next level node; Calculate the amount of water released by gravity at the nodes; Calculate the elastic water release volume of the node; The coefficients of the same groundwater head among horizontal seepage flow, vertical seepage flow, gravity release, and elastic release are added together to obtain the final groundwater head coefficient. Substitute the coefficient of the final groundwater head into the corresponding position in the pre-constructed groundwater unit permeability matrix to obtain the final groundwater unit permeability matrix. Based on Darcy's formula, the horizontal seepage flow between nodes is calculated, including: Calculate the horizontal seepage flow from node ② into node ①. : Equation (1) In the formula, It is the average hydraulic conductivity between node ② and node ①; It is the distance from the circumcenter of the triangle to the sides of the triangle containing nodes ① and ②. It is the length of the line segment between node ② and node ①; It is the water head value of node ② at time n+1. It is the head value of node ① at time n+1; It is a long period of time; Calculate the horizontal seepage flow from node ③ into node ①. : Equation (2) In the formula, It is the average hydraulic conductivity between node ③ and node ①; It is the distance from the circumcenter of the triangle to the sides of the triangle containing nodes ① and ③. It is the length of the line segment between node ③ and node ①; It is the head value of node ③ at time n+1; Calculate the horizontal seepage flow from node ③ into node ②. : Equation (3) In the formula, It is the average hydraulic conductivity between node ③ and node ②; It is the distance from the circumcenter of the triangle to the sides of the triangle containing nodes ② and ③. It is the length of the line segment between node ③ and node ②; The coefficient of groundwater head includes the coefficient generated by formula (1). = The coefficients generated by formula (1) = The coefficients generated by formula (2) = The coefficients generated by formula (2) = The coefficients generated by formula (3) = and the coefficients generated by formula (3) = .

2. The simplified method for creating a groundwater unit permeability matrix according to claim 1, characterized in that, The horizontal seepage flow from node ① into node ② is - ; The horizontal seepage flow from node ① into node ③ is - ; The horizontal seepage flow from node ② into node ③ is - .

3. The simplified method for creating a groundwater unit permeability matrix according to claim 1, characterized in that, Calculate the vertical seepage flow from the node to the next layer node, including: Let node ④ be the next node in the layer below node ①, node ⑤ be the next node in the layer below node ②, and node ⑥ be the next node in the layer below node ③. Calculate the vertical seepage flow from node ④ into node ①. Vertical seepage flow from node ⑤ into node ② Vertical seepage flow from node 6 into node 3 : Equation (4) Equation (5) Equation (6) In the formula, The vertical permeability coefficient between node ① and node ④ is given. The vertical permeability coefficient between node ② and node ⑤. The vertical permeability coefficient between node ③ and node ⑥; Let be the head value of node ④ at time n+1. Let be the head value of node ⑤ at time n+1. Let the head value of node ⑥ be at time n+1; For the festival ① and nodes The vertical distance between them For nodes and nodes The vertical distance between them For nodes and nodes The vertical distance between them; Let be the area of ​​the sub-equilibrium region of node ①. The area of ​​the sub-equilibrium region of node ②. Let be the area of ​​the sub-equilibrium region of node ③; The coefficient of groundwater head includes the coefficient generated by formula (4). = The coefficients generated by formula (4) = The coefficients generated by formula (5) = The coefficients generated by formula (5) = The coefficients generated by formula (6) = The coefficients generated by formula (6) = .

4. A simplified method for creating a groundwater unit permeability matrix according to claim 3, characterized in that, The vertical seepage flow from node ① into node ④ is - The vertical seepage flow from node ② into node ⑤ - The vertical seepage flow from node ③ into node ⑥ - .

5. A simplified method for creating a groundwater unit permeability matrix according to claim 3, characterized in that, The calculation of the gravity-induced water release at the node includes: Calculate the gravity-induced water release at node ① Gravity-induced water release at node ② and the amount of water released by gravity at node ③ : Equation (7) Equation (8) Equation (9) In the formula, It refers to the water supply degree; It is the water head value at time n of node ①. It is the water head value at time n of node ②. It is the head value at time n of node ③; The coefficient of groundwater head includes the coefficient generated by formula (7). = The coefficients generated by formula (8) = The coefficients generated by formula (9) = .

6. A simplified method for creating a groundwater unit permeability matrix according to claim 3, characterized in that, The elastic water release of the calculation node includes: When the water is under pressure, calculate the elastic water release at node ①. Elastic water release at node ② And the elastic water release of node ③ : Equation (10) Equation (11). Equation (12). In the formula, It is the elastic water release coefficient; The coefficient of groundwater head includes the generation coefficient of formula (10). = The coefficients generated by formula (11) = Formula (12) generates coefficients = .

7. A simplified method for creating a groundwater unit permeability matrix according to claim 5, characterized in that, The coefficients for the same groundwater head among horizontal seepage flow, vertical seepage flow, gravity release, and elastic release are added together to obtain the final groundwater head coefficient, which includes: The coefficients generated by formulas (1), (2), (4), and (7) Add the values ​​together, and use the coefficients generated by formula (1) Add the values ​​together, and use the coefficients generated by formula (2) Add the values ​​together, and use the coefficients generated by formula (4) Add them together to obtain the coefficients generated by formulas (3), (5), and (8). Add the values ​​together, and use the coefficients generated by formula (3) Add the values ​​together, and use the coefficients generated by formula (5) Add the values ​​together, and use the coefficients generated by formula (6) Add the values ​​together, and combine the coefficients generated by formulas (6) and (9). Adding the values ​​together yields the following: , , , , , , , , The coefficient of the final groundwater head.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that, When the processor executes the program, it implements the steps of the method according to any one of claims 1 to 7.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the steps of the method according to any one of claims 1 to 7.

Images