When the groundwater movement is simulated by numerical method, the mathematical model established should objectively reflect the hydrogeological conditions of the actual aquifer system and the basic characteristics of groundwater movement. That is to say, when natural or man-made influences are exerted, the reflection of the mathematical model should be consistent with or very close to that of the actual aquifer. Only in this way can the model be used to predict the state of aquifer system and carry out reasonable management and control. Therefore, it is necessary to test the established mathematical model, that is, to simulate the known actual process (natural groundwater level dynamics or pumping test process) by using the established differential equation and definite solution conditions, and compare the calculated value with the measured value. If the difference is too big, it is necessary to modify the differential equation or boundary conditions or parameters, and re-simulate until they are basically the same.
In fact, the core of the inversion problem is the identification of hydrogeological parameters. Because of the selection of differential equation types and the determination of definite solution conditions, it is relatively easy to solve it through hydrogeological exploration and experiment. In the inversion process, even if these two items are modified, the problem can be easily solved because of the limited types available. However, for the forward problem of numerical calculation, we can't directly use the hydrogeological parameters obtained from the analytical formula based on field test data. This is because the mathematical model of each analytical parameter formula in groundwater dynamics is quite different from the actual situation, such as the assumptions of homogeneous, isotropic, equal thickness and infinite extension of aquifer. Strictly speaking, the hydrogeological parameters obtained from the analytical formula are only applicable to the forward problem of the mathematical model corresponding to the parameter formula. In fact, the parameters correspond to the model, which is the so-called concept of model parameters. On the other hand, the hydrogeological parameters obtained from the field pumping test only represent a small area near the test point. Therefore, we must use the corresponding model to invert the hydrogeological parameters.
The concept of model parameters of 1. 1.2
The groundwater flow system in nature is very complicated in time and space. In order to solve practical problems, it is necessary to generalize the aquifer system, that is, ignore some factors that have nothing to do with the current problems or have little to do with them, and establish corresponding mathematical models. Generally speaking, it includes the following aspects:
1) regional geometric generalization of groundwater flow system;
2) Generalization of boundary properties and simulation of initial flow field;
3) Generalization of parameter properties;
4) Generalization of groundwater dynamics.
Due to the limitation of mathematical tools, hydrogeologic models that are often generalized to be very reasonable cannot be described by existing mathematical tools, or can be described but their solutions cannot be obtained. Therefore, the expedient measure is to "make do", even if the generalization is not lost, the mathematical model is simple and solvable. From the analysis here, we can see that the mathematical model established by generalizing the complex groundwater flow system is only an approximate expression of the actual groundwater flow system. As a part of the mathematical model, parameters are also the approximation and synthesis of the inherent parameters of the aquifer itself, but not the inherent parameters of the aquifer itself. In order to distinguish, we define the parameters obtained from the mathematical model as "model parameters". Obviously, when establishing a hydrogeological model for the same groundwater flow system through exploration, due to the limitations of exploration and the differences of workers' experience, different workers can establish different hydrogeological models and have different mathematical models and "model parameters" accordingly. Therefore, in practical work, it will be difficult to judge whether the parameters are correct or not, and it is impossible to generalize. For example, for the phreatic aquifer system, we can choose Bolton model or Newman model for simulation. Obviously, the number and types of parameters of the two models are different; Even parameters of the same type will have different values due to different models.
1. 1.3 well-posedness of inverse problem
For the mathematical model describing groundwater flow, if the following three items are satisfied:
1) must find its solution (existence of solution);
2) The obtained solution is unique (uniqueness of solution);
3) The solution is continuously dependent on the original data, that is, when the parameters or conditions of the solution change slightly, the change of the solution is also small (the stability of the solution), so the mathematical model is said to be well-posed.
Generally speaking, positive questions are well-posed, while inverse questions are often ill-posed. From the hydrogeology itself, there is no doubt about the existence of the solution of the inverse problem, so the following only discusses the uniqueness and stability of the solution.
1. 1.3. 1 uniqueness of the solution to the rhetorical question
Because when the hydrogeological parameters are different, it is still possible to produce the same head distribution, so it is often impossible to determine the hydrogeological parameters only by observing the head value, such as the one-dimensional steady flow model with pressure:
Identification method of aquifer parameters
Where h is the groundwater head, T(x) is the permeability coefficient, x is the coordinate, and H 1 and H2 are the given head values.
When the aquifer is homogeneous, that is, T(x)=C, its head distribution has nothing to do with the permeability coefficient t and the water storage coefficient s, but only depends on the boundary conditions. It can be seen that although the hydrogeological parameters are different, the same head distribution can still be produced.
The forward problem of the above model is to find h(x) when T(x) and boundary conditions are known. Now let's consider the inverse problem, which is to find T(x) with the known h(x). Integrate the formula (1- 1) and get:
Identification method of aquifer parameters
Namely:
Identification method of aquifer parameters
Although it is known and the constant c is arbitrary, the solution T(x) is not unique. However, if the flow of a certain section is known, for example, when x=L, not only h(x) is known, but also H (x) is known.
Identification method of aquifer parameters
Then C=-q in equation (1-3), so that T(x) is uniquely determined at x = L.
Then consider the two-dimensional steady flow under pressure, and its equation is
Identification method of aquifer parameters
Where h is the groundwater head, T(x) is the permeability coefficient, w is the source-sink term, and x and y are coordinate variables.
Equation (1-5) is rewritten as
Identification method of aquifer parameters
Binling
Identification method of aquifer parameters
Substituting the formula (1-7) into the formula (1-6), we get:
Identification method of aquifer parameters
Since every point h(x, y) and source-sink term in the region ω are known, the equation (1-8) is a first-order linear partial differential equation about t, and its general solution must contain an arbitrary function. For example, t is the solution of equation (1-8) and t' is an equation.
Identification method of aquifer parameters
T+T' must also be the solution of equation (1-8). because
Identification method of aquifer parameters
In other words, the solution of equation (1-5) is not unique.
1. 1.3.2 Instability of the solution of the inverse problem
To illustrate this problem, let's take one-dimensional steady pressure flow as an example, and its mathematical model is as follows.
Identification method of aquifer parameters
According to the formula (1-4):
Identification method of aquifer parameters
Since the regular header h*(x) always has a certain error E, we can write
Identification method of aquifer parameters
Where h(x) is the actual head value. From the formulas (1- 10) and (1-1), we get:
Identification method of aquifer parameters
Therefore, the absolute error between t and T* is
Identification method of aquifer parameters
So the size of T-T* is related to the size of sum. E is small, but it may still be big. if
Identification method of aquifer parameters
therefore
Identification method of aquifer parameters
Obviously, this is not practical. Therefore, under the above circumstances, the solution of the inverse problem is unstable, that is, a small error of the water head h(x) may lead to a large error of the required parameters.
Take the confined two-dimensional steady flow as an example
Identification method of aquifer parameters
Assuming that the head function h(x, y) and the source-sink term w in this area are known, the hydraulic conductivity T(x, y) can be obtained.
According to the formula (1-5):
Identification method of aquifer parameters
Because h(x, y) can only be approximated by actual measurement, there must be errors. Let the error be e(x, y), then the measured water level value h'=h+e(H is the true value of head), and the calculated hydraulic conductivity is
Identification method of aquifer parameters
Even if e(x, y) is small, it may still be large, so the absolute error T-T' may be large. The above discussion shows that the solution of the inverse problem is unstable, that is, a small error of water head may bring a big error to the solution T(x, y) of the inverse problem. Due to the ill-posed nature of inverse problems, constraints are usually added to solve them to avoid unreasonable situations. For example, according to common sense of hydrogeology, both permeability coefficient and hydraulic conductivity coefficient are non-negative, that is, K≥0 and T≥0. The water storage coefficient and specific yield cannot be greater than 1, that is, 0≤μ≤ 1, 0≤S≤ 1, and so on.
Neuman [2] systematically studied the non-uniqueness and instability of the solution of rhetorical questions, and then Chinese scholars Xue Yuqun, Sun Naizheng, Chen Chongxi and others [3 ~ 10] also discussed these problems. Interested readers can refer to related works.