Groundwater flow calculation method

Abstract

The present invention discloses a groundwater flow calculation method. The method comprises steps: S1. solving a numerical model of a groundwater flow by using a triangular continuous Galerkin finite element method, and forming an original triangular element mesh and obtaining water head field data of the groundwater; S2. carrying out refinement of the original triangular element mesh to construct a partial equilibrium domain; and and S3. calculating a relevant flow value of the partial equilibrium domain. The method disclosed by the invention proposes two flow calculation methods that can ensure conservation of mass in an equilibrium domain by using the continuous Galerkin finite element method to calculate a water head field of the groundwater, so as to improve the accuracy of a flow field of the groundwater based on the triangular continuous Galerkin finite element method and the relevant calculations thereof.

Description

technical field [0001] The invention relates to the technical field of groundwater calculation, in particular to a calculation method for groundwater flow. Background technique [0002] Partial equilibrium is a very important flow field property in many numerical calculations of groundwater flow. Zhang et al pointed out that in the study of two-phase flow, local equilibrium plays an important role. When Sun and Liu studied the groundwater flow field, they found that in the study of streamline drawing and groundwater movement law, local equilibrium is necessary to ensure the accuracy of calculation results. Dogrul and Kadir analyzed the importance of local equilibrium properties in groundwater resource management. When CassMiller summed up the numerical simulation of groundwater resources, he called the local equilibrium, the heterogeneity and uncertainty of aqueous media, and the simulation of unconfined water flow the three major challenges of groundwater flow simulation,...

AI Technical Summary

Problems solved by technology

Cordes and Kinzelbach discussed the local equilibrium domain and the corresponding flow calculation method in the stable situation without unit vertical compensation, no pumping and injection wells, and calculation errors, but because the assumptions were too strict, it was not suitable for in actual engineering projects

Hughes proved the equilibrium of any equilibrium area of ​​the continuous Galerkin finite element method, and proposed a flow calculation method that satisfies the equilibrium of the triangular element, but the flow value of the element boundary calculated by using his algorithm is not unique

Berger and Howington discussed the local equilibrium domain of the one-dimensional continuous Galerkin finite element method and its flow calculation method, and proposed the importance of keeping the flow calculation consistent with the calculation equation of the continuous Galerkin finite element method, but they did not extend their method to The two-dimensional and three-dimensional cases belong to the theoretical discussion

Delineating the local equilibrium domain based on the triangular Galerkin finite element method and giving the corresponding flow calculation method is one of the most urgent issues in the field of hydrogeological numerical calculation, which severely limits the visualization of groundwater flow field, pollutant migration, and groundwater flow. Computational accuracy for studies related to management and multiphase flow

Images