Simulation deduction method and system for dynamic gradual change of three-dimensional equivalent of underground water flow field

Abstract

The invention discloses a simulation deduction method and system for dynamic gradual change of a three-dimensional equivalent of an underground water flow field, and particularly relates to the field of underground water. Underground water level monitoring equipment is used for collecting underground water level data, preprocessing the collected underground water level data and integrating the data to form a unified data set; the method comprises the following steps: converting irregular underground water level data into regular three-dimensional grid data by adopting an interpolation method, generating a three-dimensional equivalent body of an underground water flow field according to an interpolation result, representing underground water level distribution of different depths and regions, carrying out dynamic simulation on the underground water flow field by utilizing a numerical simulation method, generating dynamic changes in a time sequence, and carrying out dynamic simulation on the underground water flow field. And constructing a prediction model of the underground water flow field change, deducing the change of the underground water flow field in the future by using the prediction model, and evaluating the water resource conditions in different scenes.

Description

Technical Field

[0001] The present invention relates to the field of groundwater, and more specifically, to a method and system for simulating and deducing the dynamic gradual change of three-dimensional isosurfaces of groundwater flow fields. Background Art

[0002] Groundwater is one of the important water resources on the earth, widely affecting multiple fields such as water ecology, water environment, and hydrological cycle. Under the dual influence of global climate change and human activities, the sustainable management of groundwater resources faces increasing challenges. Therefore, studying the dynamic changes of groundwater flow fields, especially the generation and visualization of three-dimensional isosurfaces, is of great significance for the rational utilization of water resources and the protection of water ecology.

[0003] Traditional groundwater monitoring mainly relies on point-like water level monitoring, which has a relatively limited coverage area and is difficult to reflect the changing trend of groundwater flow fields in real time and comprehensively. With the progress of monitoring technologies, especially the development of sensor technologies and data acquisition networks, the acquisition of groundwater level data has become more convenient and efficient. Through the method and system for simulating and deducing the dynamic gradual change of three-dimensional isosurfaces of groundwater flow fields, it serves the scientific management and sustainable utilization of groundwater resources, and has important application value especially in the fields of water resources, water environment protection, and geological exploration. Summary of the Invention

[0004] In order to overcome the above-mentioned defects of the prior art, embodiments of the present invention provide a method and system for simulating and deducing the dynamic gradual change of three-dimensional isosurfaces of groundwater flow fields to solve the problems raised in the above background art.

[0005] To achieve the above object, the present invention provides the following technical solution. The method for simulating and deducing the dynamic gradual change of three-dimensional isosurfaces of groundwater flow fields specifically includes the following steps:

[0006] Step S1: Use groundwater level monitoring equipment to collect groundwater level data, preprocess the collected groundwater level data, and integrate it to form a unified data set;

[0007] Step S2: Use an interpolation method to convert irregular groundwater level data into regular three-dimensional grid data, and generate three-dimensional isosurfaces of the groundwater flow field according to the interpolation results, representing the groundwater level distribution at different depths and regions;

[0008] Step S3: Use a numerical simulation method to dynamically simulate the groundwater flow field and generate dynamic changes in the time series;

[0009] Step S4: Construct a prediction model for the change of the groundwater flow field, use the prediction model to deduce the future change of the groundwater flow field, and evaluate the water resource status under different scenarios.

[0010] In a preferred embodiment, in step S1, groundwater level data is collected using a groundwater level monitoring device, the collected groundwater level data is preprocessed, and then integrated to form a unified data set. The specific steps are as follows:

[0011] Step A1, data collection: Monitoring points are arranged in the study area, and groundwater level data is collected at each monitoring point using a groundwater level monitoring device. Different monitoring depths are set according to different groundwater layers to obtain groundwater level data at different levels. The collection frequency is set to once per hour, and the timestamp and geographical location information of the data are recorded;

[0012] Step A2, data preprocessing: The original collected data is checked to remove outliers and incorrect data. For data obtained from different devices, unified unit conversion and standardization processing are performed to ensure that all data can be compared and analyzed under the same standard;

[0013] Step A3, data integration: A data table is created in the database for data storage. Each monitoring point, timestamp, and corresponding groundwater level data are used as record items in the data table to form a complete data set, and the data is classified according to hydrological regions and groundwater layer partitions.

[0014] In a preferred embodiment, in step S2, an interpolation method is used to convert irregular groundwater level data into regular three-dimensional grid data, and based on the interpolation results, a three-dimensional isosurface of the groundwater flow field is generated to represent the groundwater level distribution at different depths and regions. The specific steps are as follows:

[0015] Step B1, constructing a regular grid: The coordinates of each monitoring point are represented as (x i , y i , z i ), and the corresponding groundwater level value is represented as h i . Equally spaced coordinate points are selected in the x, y, and z directions respectively to form a three-dimensional regular grid, and the coordinates of each grid point are obtained as (x j , y k , z l ). The grid spacings are set as Δx, Δy, and Δz, and the coordinates of the grid points satisfy:

[0016] x j = x min + (j - 1)Δx, j = 1, 2,..., N x ;

[0017] y k = y min + (k - 1)Δy, k = 1, 2,..., N y ;

[0018] z l =z min +(l-1)Δz,l=1,2,...,N z ;

[0019] Where j and k represent the positions of the grid points in the x and y directions respectively, and N x 、N y 、N z are the total number of points of the grid in the x, y, and z directions respectively;

[0020] Step B2, interpolation calculation: For each regular grid point (x j ,y k ,z l ), and calculate the corresponding groundwater level h by using the inverse distance weighted interpolation method j,k,l , the specific calculation formula is: Among them, h i is the groundwater level at monitoring point i, d i is a regular grid of points (x j ,y k ,z l ) to the monitoring point i, the calculation formula is: (x i ,y i ,z i ) is the coordinate of the monitoring point, α is a weight factor, and the interpolation method is used to calculate the value of each grid point (x j ,y k ,z l ) for groundwater level h j,k,l The groundwater level distribution data on the entire three-dimensional grid is obtained by calculation.

[0021] In a preferred embodiment, in step S3, a numerical simulation method is used to dynamically simulate the groundwater flow field to generate dynamic changes in a time series. The specific steps are as follows:

[0022] Step C1, establish a groundwater flow model: the control equation of groundwater flow is expressed as: Groundwater flow indication Among them, S is the water storage, t is the time, Q is the flow rate, K is the permeability coefficient, is the water level gradient;

[0023] Step C2, solve the numerical equation: divide the model area into a finite number of grid cells, each grid cell represents an area in the groundwater flow field, each grid cell defines a water level and flow, set the grid division and time step to Δt, use the finite difference method to discretize the equation, and update the water level through numerical calculation: Among them, hn+1 is the water level at the next time step, h n is the water level at the current time step, A is the area of ​​the grid cell, and the flow is the gradient of the water level at the current time step n;

[0024] Step C3, dynamic simulation: set the initial condition to the groundwater level h(x, y, z, t=0) at time t=0, the boundary conditions include a fixed water level boundary and an outflow boundary, the outflow boundary is the inflow and outflow water volume of the set boundary, and the water level is updated cyclically according to the set time step Δt until the predetermined simulation time is reached. The water level data is recorded at each time step, a time series is generated, the dynamic changes of the groundwater flow field are analyzed, and the response characteristics of the flow field are identified.

[0025] In a preferred embodiment, in step S4, a prediction model for groundwater flow field changes is constructed, and the prediction model is used to deduce future groundwater flow field changes to evaluate water resource conditions under different scenarios. The specific steps are as follows:

[0026] Step D1: different water resource scenarios are set to evaluate their impact on the groundwater flow field, including a baseline scenario, a climate change scenario, and a pumping enhancement scenario. For each set scenario, a prediction model for groundwater flow field changes is constructed, further comprising the following steps:

[0027] Step D101: The baseline scenario is that groundwater recharge and extraction remain at a stable state without human intervention and climate change. The specific formula is: Where h(x,y,t) is the groundwater level at location (x,y) at time t, T 水 is the water storage capacity of the aquifer, K is the permeability coefficient, Q ext is the external recharge, which is a constant, and A is the area of ​​each grid cell;

[0028] Step D102: The climate change scenario is the change in precipitation pattern, temperature, and evaporation factors in the future. The specific formula is: Among them, Q ext (t) is a function of time t, representing the dynamic recharge under the climate change scenario, h(x,y,t) is the groundwater level at location (x,y) at time t, T 水 is the water storage capacity of the aquifer, K is the permeability coefficient, and A is the area of ​​each grid cell;

[0029] Step D103: The pumping enhancement scenario is that the groundwater extraction volume increases, which affects the change of the groundwater flow field. The specific formula is: Among them, Q withdrawis the flow rate of artificial pumping, indicating the amount of groundwater extracted, h(x,y,t) is the groundwater level at position (x,y) at time t, T 水 is the water storage capacity of the aquifer, K is the permeability coefficient, A is the area of ​​each grid cell, Q ext is the external supply, which is a constant;

[0030] Step D2: Under each scenario, numerical solutions are used to solve the problem, and the changes in water level and flow are recorded at each time step. Water level change curves and flow change curves under different scenarios are generated. The changes in water level and flow under the baseline scenario, climate change scenario, and pumping enhancement scenario are compared to analyze their impact on groundwater resources.

[0031] The present application also provides a system for simulating and deducing the dynamic gradual change of the three-dimensional equivalent volume of the groundwater flow field, which specifically includes a data acquisition and preprocessing module, a three-dimensional modeling module, a dynamic simulation and numerical simulation module, and a scenario assessment module;

[0032] The data acquisition and preprocessing module uses groundwater level monitoring equipment to collect groundwater level data, preprocesses the collected groundwater level data, and integrates them to form a unified data set;

[0033] The three-dimensional modeling module converts irregular groundwater level data into regular three-dimensional grid data using an interpolation method, and generates a three-dimensional isovoid of the groundwater flow field based on the interpolation result, indicating the distribution of groundwater levels at different depths and regions;

[0034] The dynamic simulation and numerical simulation module uses numerical simulation methods to dynamically simulate the groundwater flow field and generate dynamic changes in the time series;

[0035] The scenario assessment module constructs a prediction model for groundwater flow field changes, uses the prediction model to deduce future groundwater flow field changes, and assesses water resource conditions under different scenarios.

[0036] The beneficial effects of the present invention are as follows: The groundwater level monitoring equipment is used to collect groundwater level data, and the collected groundwater level data is preprocessed and integrated to form a unified data set. The interpolation method is adopted to convert the irregular groundwater level data into regular three-dimensional grid data, improving the spatial accuracy of the model. According to the interpolation result, a three-dimensional isosurface of the groundwater flow field is generated, representing the groundwater level distribution at different depths and regions, making the simulation of the groundwater flow field more realistic. The numerical simulation method is used to dynamically simulate the groundwater flow field, generating the dynamic changes in the time series, constructing a prediction model for the changes in the groundwater flow field, and using the prediction model to deduce the future changes in the groundwater flow field, evaluating the water resource status under different scenarios. The present invention predicts the decline of the groundwater level through a simulation model, makes emergency response measures in advance, reduces disaster losses, can improve the management efficiency and scientific nature of groundwater resources, and also provides a solid foundation for sustainable development and ecological protection. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Figure 1 It is a flowchart of the method of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0038] Next, the technical solutions in the embodiments of the present application will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative efforts fall within the scope of protection of the present application.

[0039] In the description of the present application, the terms "first" and "second" are only used for descriptive purposes and cannot be understood as indicating or implying relative importance or implicitly specifying the quantity of the indicated technical features. Thus, the features defined with "first" and "second" may explicitly or implicitly include one or more of the described features. In the description of the present application, "a plurality" means two or more, unless otherwise specifically defined.

[0040] In the description of the present application, the term "for example" is used to mean "serving as an example, illustration, or explanation". Any embodiment described in the present application as "for example" is not necessarily construed as being more preferred or advantageous than other embodiments. The following description is provided to enable any person skilled in the art to implement and use the present invention. In the following description, details are set forth for purposes of explanation. It should be understood that those of ordinary skill in the art can recognize that the present invention can be implemented without the use of these specific details. In other instances, well-known structures and processes are not elaborated in detail to avoid obscuring the description of the present invention with unnecessary details. Therefore, the present invention is not intended to be limited to the embodiments shown, but rather to be in line with the broadest scope consistent with the principles and features disclosed in the present application.

[0041] Embodiment 1

[0042] This embodiment provides a method for simulating and inferring the dynamic gradual change of the three-dimensional isosurface of the groundwater flow field as shown in Figure 1 The method specifically includes the following steps:

[0043] Step S1: Use groundwater level monitoring equipment to collect groundwater level data, preprocess the collected groundwater level data, and integrate it to form a unified data set;

[0044] Step S2: Use the interpolation method to convert the irregular groundwater level data into regular three-dimensional grid data, and generate a three-dimensional isosurface of the groundwater flow field according to the interpolation result, representing the groundwater level distribution at different depths and regions;

[0045] Step S3: Use the numerical simulation method to dynamically simulate the groundwater flow field and generate the dynamic changes in the time series;

[0046] Step S4: Build a prediction model for the change of the groundwater flow field, use the prediction model to infer the future change of the groundwater flow field, and evaluate the water resource status under different scenarios.

[0047] Preferably, in step S1, using the groundwater level monitoring equipment to collect groundwater level data, preprocess the collected groundwater level data, and integrate it to form a unified data set, the specific steps are as follows:

[0048] Step A1: Data collection: Reasonably arrange monitoring points in the study area, use groundwater level monitoring equipment to collect groundwater level data at each monitoring point, set different monitoring depths according to different groundwater layers to obtain groundwater level data at different levels, set the collection frequency to once per hour, and record the time stamp and geographical location information of the data;

[0049] Step A2, Data Preprocessing: Check the original collected data, clear outliers and error data. For data obtained from different devices, perform unified unit conversion and standardization processing to ensure that all data can be compared and analyzed under the same standard;

[0050] Step A3, Data Integration: Create a data table in the database for data storage. Take each monitoring point, timestamp, and the corresponding groundwater level data as record items of the data table, unify them into a complete data set, and classify the data according to hydrological regions and groundwater layer partitions.

[0051] Preferably, in step S2, the interpolation method is used to convert irregular groundwater level data into regular three-dimensional grid data, and according to the interpolation results, generate a three-dimensional isosurface of the groundwater flow field to represent the groundwater level distribution at different depths and regions. The specific steps are as follows:

[0052] Step B1, Construct a Regular Grid: Represent the coordinates of each monitoring point as (x i , y i , z i ) and the corresponding groundwater level value as h i . Select equally spaced coordinate points in the x, y, and z directions respectively to form a three-dimensional regular grid, and obtain the coordinates of each grid point (x j , y k , z l ). Set the grid spacing as Δx, Δy, and Δz. The coordinates of the grid points satisfy:

[0053] x j = x min + (j - 1)Δx, j = 1, 2,..., N x ;

[0054] y k = y min + (k - 1)Δy, k = 1, 2,..., N y ;

[0055] z l = z min + (l - 1)Δz, l = 1, 2,..., N z ;

[0056] Among them, j and k respectively represent the positions of the grid points in the x and y directions. N x , N y , N z are the total number of grid points in the x, y, and z directions respectively;

[0057] Step B2, Interpolation Calculation: For each regular grid point (x j , yk , z l ), calculate the corresponding groundwater level h by using the inverse distance weighted interpolation method j,k,l , and the specific calculation formula is as follows: where h i is the groundwater level value of monitoring point i, and d i is the distance from the regular grid point (x j , y k , z l ) to monitoring point i, and the calculation formula is: (x i , y i , z i ) are the coordinates of the monitoring point, α is a weight factor, and the interpolation method is used to calculate the groundwater level h j , y k , z l ) for each grid point, and the groundwater level distribution data on the entire three-dimensional grid is obtained. j,k,l

[0058] Preferably, in step S3, the numerical simulation method is used to dynamically simulate the groundwater flow field to generate the dynamic changes in the time series, and the specific steps are as follows:

[0059] Step C1: Establish a groundwater flow model: Express the control equation of groundwater flow as: The groundwater flow rate is expressed as where S is the water storage, t is the time, Q is the flow rate, K is the permeability coefficient, is the water level gradient;

[0060] Step C2: Solve the numerical equation: Divide the model area into finite grid cells, each grid cell represents an area in the groundwater flow field, each grid cell defines a water level and a flow rate, set the grid division and the time step as Δt, discretize the equation using the finite difference method, and update the water level through numerical calculation: where h n+1 is the water level at the next time step, h n is the water level at the current time step, A is the area of the grid cell, and the flow rate is the gradient of the water level at the current time step n;

[0061] ​Step C3, dynamic simulation: set the initial condition to the groundwater level h(x, y, z, t=0) at time t=0, the boundary conditions include a fixed water level boundary and an outflow boundary, the outflow boundary is the inflow and outflow water volume of the set boundary, and the water level is updated cyclically according to the set time step Δt until the predetermined simulation time is reached. The water level data is recorded at each time step, a time series is generated, the dynamic changes of the groundwater flow field are analyzed, and the response characteristics of the flow field are identified.

[0062] Preferably, in step S4, a prediction model for groundwater flow field changes is constructed, and the prediction model is used to deduce future groundwater flow field changes to evaluate water resource conditions under different scenarios. The specific steps are as follows:

[0063] Step D1: different water resource scenarios are set to evaluate their impact on the groundwater flow field, including a baseline scenario, a climate change scenario, and a pumping enhancement scenario. For each set scenario, a prediction model for groundwater flow field changes is constructed, further comprising the following steps:

[0064] Step D101: The baseline scenario is that groundwater recharge and extraction remain at a stable state without human intervention and climate change. The specific formula is: Where h(x,y,t) is the groundwater level at location (x,y) at time t, T 水 is the water storage capacity of the aquifer, K is the permeability coefficient, Q ext is the external recharge, which is a constant, and A is the area of ​​each grid cell;

[0065] Step D102: The climate change scenario is the change in precipitation pattern, temperature, and evaporation factors in the future. The specific formula is: Among them, Q ext (t) is a function of time t, representing the dynamic recharge under the climate change scenario, h(x,y,t) is the groundwater level at location (x,y) at time t, T 水 is the water storage capacity of the aquifer, K is the permeability coefficient, and A is the area of ​​each grid cell;

[0066] Step D103: The pumping enhancement scenario is that the groundwater extraction volume increases, resulting in a drop in the groundwater level and affecting the change of the groundwater flow field. The specific formula is: Among them, Q withdraw is the flow rate of artificial pumping, indicating the amount of groundwater extracted, h(x,y,t) is the groundwater level at position (x,y) at time t, T 水 is the water storage capacity of the aquifer, K is the permeability coefficient, A is the area of ​​each grid cell, Q ext is the external supply, which is a constant;

[0067] Step D2: Under each scenario, numerical solutions are used to solve the problem, and the changes in water level and flow are recorded at each time step. Water level change curves and flow change curves under different scenarios are generated. The changes in water level and flow under the baseline scenario, climate change scenario, and pumping enhancement scenario are compared to analyze their impact on groundwater resources.

[0068] Example 2

[0069] This embodiment provides a system for simulating and deducing the dynamic gradual change of a three-dimensional isovolumetric volume of a groundwater flow field, which specifically includes a data acquisition and preprocessing module, a three-dimensional modeling module, a dynamic simulation and numerical simulation module, and a scenario assessment module;

[0070] The data acquisition and preprocessing module uses groundwater level monitoring equipment to collect groundwater level data, preprocesses the collected groundwater level data, and integrates them to form a unified data set;

[0071] The three-dimensional modeling module converts irregular groundwater level data into regular three-dimensional grid data using an interpolation method, and generates a three-dimensional isovoid of the groundwater flow field based on the interpolation result, indicating the distribution of groundwater levels at different depths and regions;

[0072] The dynamic simulation and numerical simulation module uses numerical simulation methods to dynamically simulate the groundwater flow field and generate dynamic changes in the time series;

[0073] The scenario assessment module constructs a prediction model for groundwater flow field changes, uses the prediction model to deduce future groundwater flow field changes, and assesses water resource conditions under different scenarios.

[0074] It should be noted that in the above embodiments, the description of each embodiment has its own emphasis, and for parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0075] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0076] The present invention will be described with reference to the flowcharts and / or block diagrams of methods, apparatuses (systems), and computer program products according to embodiments of the present invention. It should be understood that each flow and / or block in the flowchart and / or block diagram, and the combination of flows and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded computer, or other programmable data processing devices to generate a machine, such that the instructions executed by the processor of the computer or other programmable data processing devices produce a means for implementing the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or multiple blocks.

[0077] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, such that the instructions stored in the computer-readable memory produce a manufactured article including an instruction means that implements the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or multiple blocks.

[0078] These computer program instructions can also be loaded onto a computer or other programmable data processing device, such that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process, so that the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one flow Figure 1 one flow or multiple flows and / or blocks Figure 1 or multiple blocks.

[0079] Although the preferred embodiments of the present invention have been described, those skilled in the art can make additional changes and modifications once they learn the basic creative concepts. Therefore, the appended claims are intended to be construed to include the preferred embodiments and all changes and modifications that fall within the scope of the present invention.

[0080] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalent technologies, the present invention is also intended to include these changes and modifications.

Claims

1. A method for simulating the dynamic gradual change of a three-dimensional isovolumetric volume of a groundwater flow field, characterized in that: The specific steps include: Step S1, using groundwater level monitoring equipment to collect groundwater level data, pre-processing the collected groundwater level data, and integrating them to form a unified data set; Step S2, using an interpolation method to convert irregular groundwater level data into regular three-dimensional grid data, and generating a three-dimensional isovoid of the groundwater flow field based on the interpolation result, representing the groundwater level distribution at different depths and regions; Step S3, using a numerical simulation method to dynamically simulate the groundwater flow field to generate dynamic changes in a time series; Step S4: construct a prediction model for groundwater flow field changes, use the prediction model to deduce future groundwater flow field changes, and evaluate water resource conditions under different scenarios.

2. The method for simulating and deducing the dynamic gradual change of the three-dimensional equivalent volume of the groundwater flow field according to claim 1 is characterized by: In step S1, groundwater level data is collected using groundwater level monitoring equipment, the collected groundwater level data is preprocessed, and the data is integrated to form a unified data set. The specific steps are as follows: Step A1, data collection: monitoring points are set up in the study area, and groundwater level data are collected at each monitoring point using groundwater level monitoring equipment. Different monitoring depths are set according to different groundwater layers to obtain groundwater level data at different levels. The collection frequency is set to every hour, and the timestamp and geographic location information of the data are recorded; Step A2, data preprocessing: check the original collected data, remove abnormal values ​​and erroneous data, and perform unified unit conversion and standardization processing on the data obtained from different devices to ensure that all data are compared and analyzed under the same standard; Step A3, data integration: Create a data table in the database for data storage, take each monitoring point, timestamp and corresponding groundwater level data as record items in the data table, unify them into a complete data set, and classify the data according to hydrological regions and groundwater layer divisions.

3. The method for simulating and deducing the dynamic gradual change of the three-dimensional equivalent volume of the groundwater flow field according to claim 1 is characterized by: In step S2, an interpolation method is used to convert irregular groundwater level data into regular three-dimensional grid data, and a three-dimensional isovoid of the groundwater flow field is generated according to the interpolation result to represent the groundwater level distribution at different depths and regions. The specific steps are as follows: Step B1: Construct a regular grid: The coordinates of each monitoring point are expressed as (x i ,y i ,z i ) and the corresponding groundwater level is expressed as h i , select equally spaced coordinate points in the x, y, and z directions to form a three-dimensional regular grid, and obtain the coordinates of each grid point (x j ,y k ,z l ), set the grid spacing to Δx, Δy and Δz, and the coordinates of the grid points satisfy: x j =x min +(j-1)Δx,j=1,2,...,N x ; y k =y min +(k-1)Δy,k=1,2,...,N y ; z l =z min +(l-1)Δz,l=1,2,...,N z ; Where j and k represent the positions of the grid points in the x and y directions respectively, and N x 、N y 、N z are the total number of points of the grid in the x, y, and z directions respectively; Step B2, interpolation calculation: For each regular grid point (x j ,y k ,z l ), and calculate the corresponding groundwater level h by using the inverse distance weighted interpolation method j,k,l , the specific calculation formula is: Among them, h i is the groundwater level at monitoring point i, d i is a regular grid of points (x j ,y k ,z l ) to the monitoring point i, the calculation formula is: (x i ,y i ,z i ) is the coordinate of the monitoring point, α is a weight factor, and the interpolation method is used for each grid point (x j ,y k ,z l ) for groundwater level h j,k,l The groundwater level distribution data on the entire three-dimensional grid is obtained by calculation.

4. The method for simulating and deducing the dynamic gradual change of the three-dimensional equivalent volume of the groundwater flow field according to claim 1 is characterized in that: In step S3, a numerical simulation method is used to dynamically simulate the groundwater flow field to generate dynamic changes in the time series. The specific steps are as follows: Step C1, establish a groundwater flow model: the control equation of groundwater flow is expressed as: Groundwater flow is expressed as Q = -K·▽h, where S is the water storage, t is the time, Q is the flow rate, K is the permeability coefficient, and ▽h is the water level gradient; Step C2, solve the numerical equation: divide the model area into a finite number of grid cells, each grid cell represents an area in the groundwater flow field, each grid cell defines a water level and flow, set the grid division and time step to Δt, use the finite difference method to discretize the equation, and update the water level through numerical calculation: Among them, h n+1 is the water level at the next time step, h n is the water level at the current time step, A is the area of ​​the grid cell, and Q is the flow rate. n =-K·▽h n , ▽h n is the gradient of the water level at the current time step n; Step C3, dynamic simulation: set the initial condition to the groundwater level h(x, y, z, t=0) at time t=0, the boundary conditions include a fixed water level boundary and an outflow boundary, the outflow boundary is the inflow and outflow water volume of the set boundary, and the water level is updated cyclically according to the set time step Δt until the predetermined simulation time is reached. The water level data is recorded at each time step, a time series is generated, the dynamic changes of the groundwater flow field are analyzed, and the response characteristics of the flow field are identified.

5. The method for simulating and deducing the dynamic gradual change of the three-dimensional equivalent volume of the groundwater flow field according to claim 1 is characterized in that: In step S4, a prediction model for groundwater flow field changes is constructed, and the prediction model is used to deduce future groundwater flow field changes and evaluate water resource conditions under different scenarios. The specific steps are as follows: Step D1: Set different water resource scenarios to evaluate their impact on the groundwater flow field, including baseline scenario, climate change scenario and pumping enhancement scenario. For each scenario, construct a prediction model for groundwater flow field changes; Step D2: Under each scenario, numerical solutions are used to solve the problem, and the changes in water level and flow are recorded at each time step. Water level change curves and flow change curves under different scenarios are generated. The changes in water level and flow under the baseline scenario, climate change scenario, and pumping enhancement scenario are compared to analyze their impact on groundwater resources.

6. The method for simulating and deducing the dynamic gradual change of the three-dimensional equivalent volume of the groundwater flow field according to claim 5 is characterized by: The baseline scenario is that groundwater recharge and extraction remain at a stable state without human intervention and climate change. The specific formula is: Where h(x,y,t) is the groundwater level at location (x,y) at time t, T 水 is the water storage capacity of the aquifer, K is the permeability coefficient, Q ext is the external recharge, which is a constant, and A is the area of ​​each grid cell.

7. The method for simulating and deducing the dynamic gradual change of the three-dimensional isovolumetric volume of groundwater flow field according to claim 5, characterized in that: The climate change scenario is that the precipitation pattern, temperature and evaporation factors will change in the future. The specific formula is: Among them, Q ext (t) is a function of time t, representing the dynamic recharge under the climate change scenario, h(x,y,t) is the groundwater level at location (x,y) at time t, T 水 is the water storage capacity of the aquifer, K is the hydraulic conductivity, and A is the area of ​​each grid cell.

8. The method for simulating and deducing the dynamic gradual change of the three-dimensional equivalent volume of the groundwater flow field according to claim 5 is characterized by: The pumping enhancement scenario is that the groundwater extraction volume increases, affecting the change of groundwater flow field. The specific formula is: Among them, Q withdraw is the flow rate of artificial pumping, indicating the amount of groundwater extracted, h(x,y,t) is the groundwater level at position (x,y) at time t, T 水 is the water storage capacity of the aquifer, K is the permeability coefficient, A is the area of ​​each grid cell, Q ext is the external supply, which is a constant.

9. The system for simulating and deducing the dynamic gradual change of the three-dimensional equivalent volume of the groundwater flow field is applied to the method for simulating and deducing the dynamic gradual change of the three-dimensional equivalent volume of the groundwater flow field as claimed in any one of claims 1 to 8, characterized in that: It includes data acquisition and preprocessing module, 3D modeling module, dynamic simulation and numerical simulation module, and scenario assessment module; The data acquisition and preprocessing module uses groundwater level monitoring equipment to collect groundwater level data, preprocesses the collected groundwater level data, and integrates them to form a unified data set; The three-dimensional modeling module converts irregular groundwater level data into regular three-dimensional grid data using an interpolation method, and generates a three-dimensional isovoid of the groundwater flow field based on the interpolation result, indicating the distribution of groundwater levels at different depths and regions; The dynamic simulation and numerical simulation module uses numerical simulation methods to dynamically simulate the groundwater flow field and generate dynamic changes in the time series; The scenario assessment module constructs a prediction model for groundwater flow field changes, uses the prediction model to deduce future groundwater flow field changes, and assesses water resource conditions under different scenarios.

Images