Novel simulation method for representing space-time evolution of underground water level

Abstract

The invention discloses a novel simulation method for representing underground water level space-time evolution, and belongs to the technical field of computer machine learning. A novel simulation method for representing underground water level space-time evolution comprises the following steps: acquiring long-time sequence underground water level data based on an underground water long-term observation hole; according to the underground water level information, multiple hydrogeological parameter information is superposed, a two-dimensional underground water seepage equation is used as a physical mechanism, a neural network model with the physical mechanism in a region is constructed, and then time and space evolution of the underground water level in the region is simulated; reasonable use of groundwater resources can be enhanced, scientific and effective supervision is implemented, and a scientific basis is provided for protecting the safety of regional groundwater and promoting sustainable development of the groundwater resources.

Description

Technical Field

[0001] The invention relates to the technical field of groundwater simulation, and in particular to a novel simulation method for characterizing the temporal and spatial evolution of groundwater levels. Background Art

[0002] As one of the most important freshwater resources on earth, groundwater not only affects the ecological environment, but also controls the development of social economy. Especially in areas lacking surface water resources, groundwater is often the only reliable water source. Globally, more than one-third of the water used by humans comes from groundwater. However, affected by human activities, climate change and other factors, the shortage of groundwater resources is becoming increasingly serious, causing many geological environmental problems such as land subsidence, groundwater pollution and ecosystem degradation. In order to effectively manage and protect groundwater resources, the simulation and prediction of groundwater is imminent.

[0003] Traditional groundwater simulation methods require clear discretization equations and the construction of networks based on the geometric structure and boundary conditions of the problem. For complex geometric shapes or multi-physics field coupling problems, modeling and solving may be very complex and time-consuming. Traditional groundwater simulation methods are mostly based on control equations, so there are obvious limitations in dealing with nonlinear problems. Compared with traditional numerical simulation methods, machine learning and deep learning algorithms have obvious advantages in solving nonlinear problems and improve computational efficiency. However, these methods are all driven by data, and the model lacks physical constraints, which makes the internal mechanism of the model difficult to explain, and the reliability of the model results is questioned. In order to solve the above problems, the present invention proposes a new simulation method for characterizing the spatiotemporal evolution of groundwater levels. Summary of the invention

[0004] The purpose of the present invention is to obtain the changing information of groundwater in time series and to quantitatively describe the evolution law, thereby proposing a new simulation method for characterizing the spatiotemporal evolution of groundwater level and performing efficient groundwater level simulation.

[0005] In order to achieve the above object, the present invention adopts the following technical solutions: A new simulation method for characterizing the temporal and spatial evolution of groundwater levels includes the following steps: S1. Obtain monthly precipitation data covering the study area through long-term satellite remote sensing precipitation products; obtain monthly groundwater level information through long-term groundwater observation holes within the study area.

[0006] S2. Obtain basic hydrogeological data within the study area, combine it with the groundwater control equation, and construct a two-dimensional groundwater seepage equation for the study area.

[0007] S3. Based on the two-dimensional groundwater seepage equation in the study area and combined with the neural network model, a deep learning neural network model integrating physical knowledge - PINNs model is created.

[0008] S4, using the precipitation data and groundwater level data obtained in S1 as input to train the deep learning neural network model integrating physical knowledge - PINNs model, and obtaining the data error term between the network output and the training data; S5, bringing the network output result obtained by the model in S4 into the physical equation, solving the physical equation by automatic differentiation, and obtaining the physical residual term; S6. Superimpose the data error term obtained in S4 and the physical residual term obtained in S5 to obtain the overall loss of the PINNs model, which is back-propagated back to the neural network model. Learning is performed through model iteration to achieve accurate simulation of the groundwater level.

[0009] Preferably, when S1 obtains long-time satellite remote sensing precipitation products, it uses geoscientific related software for processing, and the geoscientific related processing software is ArcGIS, QGIS or GeoScence .

[0010] Preferably, it is characterized in that said S2 constructing the two-dimensional groundwater seepage equation of the study area specifically comprises the following steps: S2.1. Based on the groundwater control equation, the physical laws used for the neural network constraints are determined. Under heterogeneous isotropic conditions, the basic differential equation for the unstable movement of groundwater can be expressed as:

[0011] Where h represents the groundwater level, represents the water storage rate of the water-containing medium; and are the permeability coefficients of the aquifer in the x and y directions, respectively. Under isotropic conditions, for any point in space, the value of K is the same in all directions ( ), but in space, the value of K will be different - heterogeneity; W is the source and sink term, which represents the influence of external recharge or pumping; S2.2, the equation satisfies the following conditions under the initial conditions:

[0012] In the formula, represents the groundwater level at the initial moment of the model; t represents time; Preferably, the S3 constructs a deep learning neural network model incorporating physical knowledge, and the construction specifically includes the following steps: S3.1. Construct the model input layer to receive groundwater level information including: water level height, water level position; physical parameter information including permeability coefficient; source and sink information. The input layer is connected to the first hidden layer through a fully connected layer. If the input data is recorded as X, the output data is:

[0013] In the formula, represents the output of the first hidden layer, represents the weight matrix of the first layer, represents the hidden layer bias vector of the first layer, Represents the activation function of the network, common ones include ReLU, Sigmoid and Tanh; S3.2 When there are multiple hidden layers, each hidden layer is connected to the previous layer through a fully connected layer:

[0014] In the formula, represents the output of the nth hidden layer, represents the weight matrix of the nth hidden layer, represents the bias vector of the nth hidden layer, represents the activation function of the network, represents the input of the nth hidden layer, which is also the output of the n-1th hidden layer; S3.3. Finally, the hidden layer and the output layer are also connected through a fully connected layer. The output layer can be expressed as:

[0015] In the formula, is the prediction result of the network, is the weight matrix of the output layer, is the bias vector of the output layer, is the output of the last hidden layer; S3.4. Through automatic differentiation technology, the output layer prediction results are connected with the physical equations, and the physical equations are used as global constraints to guide the neural network training process.

[0016] Preferably, the data error term between the PINNs model network output and the training data in S4 is calculated by model MSE:

[0017] In the formula, represents the true value of the groundwater level, represents the groundwater value predicted by the model, and N represents the sample size.

[0018] Preferably, the process of solving the physical equation in S5 to obtain the calculation of the physical residual term specifically includes the following steps: S5.1. Water level values ​​predicted using the network , and the partial derivatives are calculated by automatic differentiation: Compute the first-order derivative: , and ; Compute the second-order derivative: and ; S5.2. Calculate the loss function: Substitute the derivative of the network output into the physical equation to calculate:

[0019] In the formula, N represents the sample size.

[0020] Preferably, the model back propagation and parameter iterative update in S6 specifically include the following steps: S6.1. The total loss function of the PINNs network is composed of the physical residual term of the physical equation and the data error term of the neural network model, which are weighted as follows:

[0021] In the formula, represents the total loss term of the PINNs model, and is the weight coefficient used to balance the data error term and the physical residual, and are the physical residual term and the data error term respectively; S6.2. In the process of back propagation, The gradient of the loss with respect to the network output is calculated:

[0022] In the formula, represents the total loss term of the PINNs model, and is the weight coefficient, and are the physical residual term and the data error term respectively, N is the sample size, Indicates the true groundwater level value, represents the groundwater level value predicted by the model; S6.3. Using the calculated gradient, the weights in the network are updated through an optimization algorithm, and the training is iterated continuously until the model loss function converges; the optimization algorithm is SGD or Adam.

[0023] Compared with the prior art, the present invention provides a new simulation method for characterizing the temporal and spatial evolution of groundwater levels, which has the following beneficial effects: The present invention obtains regional precipitation information based on remote sensing image data; further obtains long-term groundwater level data based on long-term groundwater observation wells; based on the groundwater level information, multiple hydrogeological parameter information is superimposed, and a two-dimensional groundwater seepage equation is used as the physical mechanism to construct a regional neural network model with a physical mechanism, thereby efficiently simulating the temporal and spatial evolution of the regional groundwater level, providing a scientific basis for protecting the safety of regional groundwater and promoting the sustainable development of groundwater resources. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] Figure 1 A schematic diagram of a novel simulation method for characterizing the temporal and spatial evolution of groundwater levels proposed by the present invention; Figure 2 The network structure of a new neural network model for simulating the spatiotemporal evolution of groundwater level proposed by the present invention; Figure 3 This is a schematic diagram of the contribution of each characteristic factor to the network. DETAILED DESCRIPTION

[0025] The technical solutions in the embodiments of the present invention will be described clearly and completely below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all the embodiments.

[0026] As one of the most important freshwater resources on earth, groundwater not only affects the ecological environment, but also controls the development of social economy. Traditional groundwater simulation methods mainly rely on numerical solutions, such as finite element method (FEM), finite difference method (FDM) and finite volume method (FVM). Such methods require explicit discretization of equations and construction of networks based on the geometric structure and boundary conditions of the problem. For complex geometries or multi-physics coupling problems, modeling and solving can be very complex and time-consuming. In addition, such methods show obvious limitations in dealing with nonlinear problems. In order to solve the above problems, the present invention proposes a new simulation method for characterizing the spatiotemporal evolution of groundwater levels.

[0027] PINNs is a neural network framework that integrates physical knowledge into deep learning. The core idea of ​​PINNs is to embed physical laws into the loss function of the neural network, so that the model can not only fit the data, but also guide the training process through the constraints of the physical process, thereby enhancing the interpretability of the model. The prediction results of the model can be explained by physical laws, thereby reducing the "black box" characteristics of the neural network. Because the embedded physical mechanism is universal, PINNs can show better generalization ability when dealing with unseen physical scenes or conditions. Although PINNs have been successful in many fields, their application research in groundwater simulation is still relatively limited. The present invention generalizes the groundwater flow in the study area into a two-dimensional unsteady flow problem, integrates the groundwater time series water level, water level geological parameters and precipitation data, and explores the potential of PINNs in groundwater level simulation and prediction by taking the basic differential equations of unsteady groundwater motion as neural network constraints, providing a new method for accurate and efficient groundwater simulation. Specifically include the following content.

[0028] Embodiment 1: See also Figure 1 , a new simulation method for characterizing the spatiotemporal evolution of groundwater level, including the following steps; S1. Obtain monthly precipitation data covering the study area through long-term satellite remote sensing precipitation products; obtain monthly groundwater level information through long-term groundwater observation holes within the study area.

[0029] S1 When acquiring remote sensing image data, it is necessary to use geoscientific related software for processing. The geoscientific related processing software is ArcGIS, QGIS or GeoScence .

[0030] S2. Obtain basic hydrogeological data within the study area, combine it with the groundwater control equation, and construct a two-dimensional groundwater seepage equation for the study area.

[0031] S2.1. Based on the groundwater control equation, the physical laws used for the neural network constraints are determined. Under heterogeneous isotropic conditions, the basic differential equation for the unstable movement of groundwater can be expressed as:

[0032] Where h represents the groundwater level, represents the water storage rate of the water-containing medium; and are the permeability coefficients of the aquifer in the x and y directions, respectively. Under isotropic conditions, for any point in space, the value of K is the same in all directions ( ), but in space, the value of K will be different - heterogeneity; W is the source and sink term, which represents the influence of external recharge or pumping; S2.2, the equation satisfies the following conditions under the initial conditions:

[0033] In the formula, represents the groundwater level at the initial moment of the model; t represents time; S3. Based on the two-dimensional groundwater seepage equation in the study area and combined with the neural network model, a deep learning neural network model integrating physical knowledge - the PINNs model is constructed.

[0034] S3.1. Construct the model input layer to receive groundwater level information including: water level height, water level position; physical parameter information including permeability coefficient; source and sink information. The input layer is connected to the first hidden layer through a fully connected layer. If the input data is recorded as X, the output data is:

[0035] In the formula, represents the output of the first hidden layer, represents the weight matrix of the first layer, represents the hidden layer bias vector of the first layer, Represents the activation function of the network, common ones include ReLU, Sigmoid and Tanh; S3.2 When there are multiple hidden layers, each hidden layer is connected to the previous layer through a fully connected layer:

[0036] In the formula, represents the output of the nth hidden layer, represents the weight matrix of the nth hidden layer, represents the bias vector of the nth hidden layer, represents the activation function of the network, represents the input of the nth hidden layer, which is also the output of the n-1th hidden layer; S3.3. Finally, the hidden layer and the output layer are also connected through a fully connected layer. The output layer can be expressed as:

[0037] In the formula, is the prediction result of the network, is the weight matrix of the output layer, is the bias vector of the output layer, is the output of the last hidden layer; S3.4. Through automatic differentiation technology, the output layer prediction results are connected with the physical equations, and the physical equations are used as global constraints to guide the neural network training process.

[0038] Preferably, the data error term between the PINNs model network output and the training data in S4 is calculated by model MSE:

[0039] In the formula, represents the true value of the groundwater level, represents the groundwater value predicted by the model, and N represents the sample size.

[0040] S4. The precipitation data and groundwater level data obtained in S1 are used as input to train the deep learning neural network model integrating physical knowledge - PINNs model, and the data error term between the network output and the training data is obtained.

[0041] The data error term between the PINNs model network output and the training data in S4 is calculated by the model MSE:

[0042] In the formula, represents the true value of the groundwater level, represents the groundwater value predicted by the model, and N represents the sample size.

[0043] S5. Substitute the network output result obtained from the model in S4 into the physical equation, solve the physical equation through automatic differentiation, and obtain the physical residual term.

[0044] S5.1. Water level values ​​predicted using the network , and the partial derivatives are calculated by automatic differentiation: Compute the first-order derivative: , and ; Compute the second-order derivative: and ; S5.2. Calculate the loss function: Substitute the derivative of the network output into the physical equation to calculate:

[0045] In the formula, N represents the sample size.

[0046] S6. Superimpose the data error term obtained in S4 and the physical residual term obtained in S5 to obtain the overall loss of the PINNs model, which is back-propagated back to the neural network model. Learning is performed through model iteration to achieve accurate simulation of the groundwater level.

[0047] S6.1. The total loss function of the PINNs network is composed of the physical residual term of the physical equation and the data error term of the neural network model, which are weighted as follows:

[0048] In the formula, represents the total loss term of the PINNs model, and is the weight coefficient used to balance the data error term and the physical residual, and are the physical residual term and the data error term respectively; S6.2. In the process of back propagation, The gradient of the loss with respect to the network output is calculated:

[0049] In the formula, represents the total loss term of the PINNs model, and is the weight coefficient, and are the physical residual term and the data error term respectively, N is the sample size, Indicates the true groundwater level value, represents the groundwater level value predicted by the model; S6.3. Using the calculated gradient, the weights in the network are updated through an optimization algorithm, and the training is iterated continuously until the model loss function converges; the optimization algorithm is SGD or Adam.

[0050] Embodiment 2: Based on Example 1, but different from Example 1, the present invention obtains monthly precipitation data from 1984 to 2023 through ERA5-Land remote sensing images on the GEE platform. Linyi City, Shandong Province is used as the study area, the water level of the existing long-term groundwater observation well is obtained, its time series evolution law is analyzed, and simulation is performed. The main steps are as follows: Step 1: Through long-term satellite remote sensing precipitation products, monthly precipitation data with a resolution of 0.1° covering the study area are obtained; combined with long-term groundwater observation wells within the study area, monthly groundwater level information is obtained.

[0051] Step 2: The basic hydrogeological data within the study area were obtained and combined with the groundwater control equation to construct a two-dimensional groundwater seepage equation for the study area.

[0052] The permeability coefficient describes the water permeability of underground media and directly affects the flow velocity and flow distribution of groundwater. The water storage rate describes the ability to release or store water volume from a column with a unit horizontal area and height equal to the thickness of the aquifer when the water head of the aquifer drops or rises by a unit height.

[0053] Step 3: According to the two-dimensional groundwater seepage equation in the study area and combined with the neural network model, a deep learning neural network model incorporating physical knowledge - the PINNs model - was created.

[0054] Step 4: Input data, train the deep learning neural network model that integrates physical knowledge, and obtain the data error term between the network output and the training data.

[0055] Step 5: The network output results obtained by the model are substituted into the physical equation, and the physical equation is solved by automatic differentiation to obtain the physical residual term.

[0056] Step 6: The data error term of the neural network is superimposed with the physical residual term in the process of solving the physical equation to obtain the overall loss term of the PINNs model, which is backpropagated back to the neural network model. Through model iteration, learning is performed to achieve accurate simulation of the groundwater level.

[0057] Step 7: The constructed groundwater level simulation model was interpreted using SHAP, and the importance results of each characteristic factor in the PINNs network were obtained.

[0058] The model simulation results show that the PINNs model score is 0.87 and the RMSE is 1.98 m. The results show that the model has a good fit and high accuracy in groundwater level simulation.

[0059] Depend on Figure 3 It can be seen that in the PINNs network, the contribution of the permeability coefficients Kx and Ky to the network is significantly higher than that of the traditional MLP network, followed by the spatial position. This shows that for the PINNs network, the model prediction pays more attention to the physical process than the spatial position, which enhances the interpretability of the model prediction results.

[0060] The above are only preferred specific implementation modes of the present invention, but the protection scope of the present invention is not limited thereto. Any technician familiar with the technical field can make equivalent replacements or changes according to the technical solutions and inventive concepts of the present invention within the technical scope disclosed by the present invention, which should be covered by the protection scope of the present invention.

Claims

1. A new simulation method for characterizing the spatiotemporal evolution of groundwater levels, characterized in that: The following steps are included: S1. Obtain monthly precipitation data covering the study area through long-term satellite remote sensing precipitation products; obtain monthly groundwater level information through long-term groundwater observation holes within the study area; S2. Obtain basic hydrogeological data within the study area, combine it with the groundwater control equation, and construct a two-dimensional groundwater seepage equation for the study area; S3. Based on the two-dimensional groundwater seepage equation in the study area and combined with the neural network model, a deep learning neural network model (PINNs model) incorporating physical knowledge is created; S4, using the precipitation data and groundwater level data obtained in S1 as input to train the deep learning neural network model PINNs model integrating physical knowledge, and obtaining the data error term between the network output and the training data; S5, bringing the network output result obtained by the model in S4 into the physical equation, solving the physical equation by automatic differentiation, and obtaining the physical residual term; S6. Superimpose the data error term obtained in S4 and the physical residual term obtained in S5 to obtain the overall loss of the PINNs model, which is back-propagated back to the neural network model. Learning is performed through model iteration to achieve accurate simulation of the groundwater level.

2. A new simulation method for characterizing the temporal and spatial evolution of groundwater levels according to claim 1, characterized in that: When S1 obtains long-term satellite remote sensing precipitation products, it uses geoscientific related software for processing. The geoscientific related processing software is ArcGIS, QGIS or GeoScence .

3. A new simulation method for characterizing the temporal and spatial evolution of groundwater levels according to claim 1, characterized in that: The S2 construction of the two-dimensional groundwater seepage equation in the study area specifically includes the following steps: S2.

1. Based on the groundwater control equation, the physical laws used for the neural network constraints are determined. Under heterogeneous isotropic conditions, the basic differential equation for the unstable movement of groundwater can be expressed as: Where h represents the groundwater level, represents the water storage rate of the water-containing medium; and are the permeability coefficients of the aquifer in the x and y directions, respectively. Under isotropic conditions, for any point in space, the value of K is the same in all directions ( ), but in space, the value of K will be heterogeneous; W is the source-sink term, indicating the influence of external recharge or pumping; S2.2, the equation satisfies the following conditions under the initial conditions: In the formula, represents the groundwater level at the initial moment of the model; t represents time.

4. A new simulation method for characterizing the temporal and spatial evolution of groundwater levels according to claim 1, characterized in that: The S3 constructs a deep learning neural network model incorporating physical knowledge, specifically including the following steps: S3.

1. Construct the model input layer to receive groundwater level information including: water level height, water point location; physical parameter information including permeability coefficient; source and sink information. The input layer is connected to the first hidden layer through a fully connected layer. If the input data is recorded as X, the output data is: In the formula, represents the output of the first hidden layer, represents the weight matrix of the first layer, represents the hidden layer bias vector of the first layer, Represents the activation function of the network, common ones include ReLU, Sigmoid and Tanh; S3.2 When there are multiple hidden layers, each hidden layer is connected to the previous layer through a fully connected layer: In the formula, represents the output of the nth hidden layer, represents the weight matrix of the nth hidden layer, represents the bias vector of the nth hidden layer, represents the activation function of the network, represents the input of the nth hidden layer, which is also the output of the n-1th hidden layer; S3.

3. Finally, the hidden layer and the output layer are connected through a fully connected layer; the output layer can be expressed as: In the formula, is the prediction result of the network, is the weight matrix of the output layer, is the bias vector of the output layer, is the output of the last hidden layer; S3.

4. Through automatic differentiation technology, the output layer prediction results are connected with the physical equations, and the physical equations are used as global constraints to guide the neural network training process.

5. A new simulation method for characterizing the temporal and spatial evolution of groundwater level according to claim 1, characterized in that: The data error term between the PINNs model network output and the training data in S4 is calculated by the model MSE: In the formula, represents the true value of the groundwater level, represents the groundwater value predicted by the model, and N represents the sample size.

6. A new simulation method for characterizing the temporal and spatial evolution of groundwater levels according to claim 1, characterized in that: The process of solving the physical equation in S5 to obtain the calculation of the physical residual term specifically includes the following steps: S5.

1. Water level values ​​predicted using the network , and the partial derivatives are calculated by automatic differentiation: Compute the first-order derivative: , and ; Compute the second-order derivative: and ; S5.

2. Calculate the loss function: Substitute the derivative of the network output into the physical equation to calculate: In the formula, N represents the sample size.

7. A new simulation method for characterizing the temporal and spatial evolution of groundwater levels according to claim 1, characterized in that: The overall loss of the PINNs model is obtained in S6, which is back-propagated back to the neural network model and parameter iteration for learning, specifically including the following steps: S6.

1. The total loss function of the PINNs network is composed of the physical residual term of the physical equation and the data error term of the neural network model, which are weighted as follows: In the formula, represents the total loss term of the PINNs model, and is the weight coefficient used to balance the data error term and the physical residual, and are the physical residual term and the data error term respectively; S6.

2. In the process of back propagation, The gradient of the loss with respect to the network output is calculated: In the formula, represents the total loss term of the PINNs model, and is the weight coefficient, and are the physical residual term and the data error term respectively, N is the sample size, Indicates the true groundwater level value, represents the groundwater level value predicted by the model; S6.

3. Using the calculated gradient, the weights in the network are updated through an optimization algorithm, and the training is iterated continuously until the model loss function converges; the optimization algorithm is SGD or Adam.

8. A new simulation method for characterizing the spatiotemporal evolution of groundwater levels according to any one of claims 2 to 6, characterized in that: The calculations are obtained by Python or Matlab.

Images