Solute transport calculation method based on hydrogeological model

Abstract

The invention discloses a solute transport calculation method based on a hydrogeological model. The solute transport calculation method comprises the following steps: performing three-dimensional space superposition on a hydrogeological grid in a grid model obtained in solute transport simulation and a hydrogeological model; each hydrogeological grid automatically adjusts the top and bottom plate height and the layered structure of the corresponding cell according to the top and bottom plate elevation and the layered structure of the corresponding stratum in the superposed hydrogeological model, and performs parameter assignment on the corresponding hydrogeological grid according to the hydrogeological parameters in the hydrogeological model; assigning values to other simulation data of the hydrogeological grids in the grid model according to the actual condition of the simulation area; parameters required by migration simulation calculation are extracted from the grid model, and migration simulation calculation of groundwater flow and solute in groundwater is carried out through a migration simulation calculation module. According to the invention, the working efficiency of solute transport simulation is improved, and the accuracy and consistency of hydrogeological data are ensured.

Description

Technical Field

[0001] The invention relates to the field of geological information system management software development, and in particular to a solute transport calculation method based on a hydrogeological model. Background Art

[0002] A hydrogeological model is a model used to study the interaction between groundwater and geological conditions. It is a three-dimensional model that combines knowledge of geology and hydrology and uses information such as hydrogeological strata and groundwater resources. Through hydrogeological models, researchers can predict the flow direction, speed and water quality changes of groundwater, as well as evaluate the sustainable utilization of groundwater resources. These models play an important role in groundwater management, environmental protection and water resources planning. The hydrogeological model construction technology we are currently using is based on the three-dimensional layered geological body automatic modeling technology of ZGIS.

[0003] The attribute data in the hydrogeological model is used to describe the permeability characteristics of underground rock formations or soils, including parameters such as permeability, porosity, and pore structure, which play an important role in simulating groundwater flow and hydrological cycle processes. Solute transport simulation is a simulation method used to study the migration and transmission process of dissolved substances (solutes) in groundwater. This simulation is usually based on a model of groundwater flow and takes into account processes such as diffusion, adsorption, reaction, and degradation of solutes in groundwater. Through solute transport simulation, researchers can predict the concentration distribution, migration speed, and diffusion range of solutes in groundwater, as well as evaluate the response and impact of groundwater pollution. These simulations have important application value in groundwater pollution assessment, environmental risk assessment, and groundwater resource management.

[0004] At present, the tools used for solute transport simulation (hereinafter referred to as simulation) in China are mainly MODFLOW and the calculation tool software based on MODFLOW secondary development. The data input, data processing operations and steps required for simulation vary according to the development manufacturers, and the existing hydrogeological models cannot be directly used. In the simulation, the corresponding hydrogeological data often need to be re-entered, which limits the efficiency of data analysis and application. Even based on the same hydrogeological data, due to the different abstract understanding of the model by the users, the preprocessing and input of the data required for simulation are also different, which affects the accuracy and consistency of the simulation results. Summary of the invention

[0005] The present invention provides a solute transport calculation method based on a hydrogeological model, which aims to solve the problem that the simulation program data input is cumbersome and cannot be well combined with the existing hydrogeological model. The method uses a three-dimensional spatial analysis method to perform spatial superposition analysis on the grids divided in the simulation and the hydrogeological model, automatically brings the existing hydrological attributes in the model into the simulated grid, and automatically adjusts the simulated grid structure based on the layered structure of the hydrogeological model body, thereby achieving the purpose of automatically constructing most of the data required for the simulation. The specific scheme is as follows:

[0006] A solute transport calculation method based on a hydrogeological model, the method comprising:

[0007] Step 1, superimposing the hydrogeological grid in the grid model obtained in the solute transport simulation with the hydrogeological model in three-dimensional space;

[0008] Step 2, each hydrogeological grid automatically adjusts the top and bottom plate heights and layered structures of the corresponding cells according to the top and bottom plate elevations and layered structures of the corresponding strata in the superimposed hydrogeological model, and assigns parameters to the corresponding hydrogeological grids using the hydrogeological parameters in the hydrogeological model;

[0009] Step 3, assigning values ​​to other simulation data of the hydrogeological grid in the grid model according to the actual conditions of the simulation area, wherein the other simulation data include simulation initial condition data, historical simulation data and hydrological and meteorological data;

[0010] Step 4: extract the parameters required for the transport simulation operation from the grid model, and perform groundwater flow and solute transport simulation calculations in groundwater through the transport simulation calculation module.

[0011] Furthermore, in step 1, a hydrogeological model is established by collecting hydrogeological data; or an already constructed hydrogeological model file is directly imported.

[0012] Furthermore, step 1 also includes:

[0013] Initialize the solute transport simulation, including: setting the measurement units of the basic simulation data and the default parameters of the solver equations before starting the simulation calculation;

[0014] Generate a grid model based on the grid accuracy and simulation range in the basic data of the simulation.

[0015] Furthermore, the method further comprises: after the three-dimensional space superposition, using the hydrogeological model to automatically input boundary conditions, specifically comprising:

[0016] Obtain a three-dimensional view after three-dimensional spatial superposition, perform spatial superposition analysis based on the hydrogeological model, automatically identify the environmental information in the hydrogeological model, superimpose the environmental information in the hydrogeological model with the grid model in the three-dimensional view, and assign corresponding boundary conditions to the corresponding hydrogeological grid based on the environmental information.

[0017] Furthermore, the environmental information includes information about lakes, rivers, oceans and groundwater bodies.

[0018] Furthermore, the method also includes: after adjusting the top and bottom plate heights and layered structures of the corresponding cells and assigning parameters to the corresponding hydrogeological grids, artificially adjusting the grid parameters and grid structure of the hydrogeological grids through human-computer interaction based on the actual hydrogeological conditions.

[0019] Furthermore, other simulation data include simulation initial condition data, historical simulation data and hydrological and meteorological data. When there is no simulation initial condition data, the historical simulation data is used as the simulation initial condition data for this simulation. The hydrological and meteorological data includes hydrological and meteorological information such as rainfall, evaporation, and seepage.

[0020] Furthermore, in step 4, MODFLOW and MT3DMS are used as transport simulation calculation modules, and the input parameters required by MODFLOW and MT3DMS are extracted from the grid model to perform simulation calculations of groundwater flow and solute transport in groundwater.

[0021] Furthermore, the method also includes: in step 4, the migration simulation calculation results are discrete point head and concentration data with time series, the spatial position information of the discrete points corresponds to the cell space information in the grid model, and after the calculation results are read, the calculation data is associated with the cells in the grid model and stored in the database.

[0022] Furthermore, the method also includes visually displaying the calculation results, specifically including:

[0023] The calculated data is extracted from the database, and the center points of the cells in the grid model are used as the spatial coordinates of the discrete points. The simulation results are converted into three-dimensional models, plane diagrams, and solute transport animations using spatial interpolation, and corresponding statistical charts are generated.

[0024] The present invention has the following beneficial effects:

[0025] The present invention provides a solute transport calculation method based on a hydrogeological model, which can greatly improve the efficiency and accuracy of simulation work. Compared with manual input of hydrogeological data, this method can save researchers time and energy in parsing and constructing hydrogeological models, and facilitate researchers to establish solute transport simulation calculations more quickly. At the same time, the system also provides a visualization tool, allowing users to more intuitively understand and analyze solute transport simulations. Through automated hydrogeological data input, the workflow of solute transport simulations can be reduced, and groundwater flow and hydrological cycle processes can be simulated more accurately, providing important support for groundwater pollution assessment, environmental risk assessment, and groundwater resource management. BRIEF DESCRIPTION OF THE DRAWINGS

[0026] Figure 1 A schematic flow chart of a solute transport calculation method based on a hydrogeological model provided in an embodiment of the present invention;

[0027] Figure 2 A schematic diagram of a hydrogeological model provided by an embodiment of the present invention;

[0028] Figure 3 A schematic diagram of a grid model for adjusting the top and bottom plate heights of a unit cell and a layered structure provided in an embodiment of the present invention.

[0029] Figure 4 The present invention provides an overall flow chart of the embodiment of the present invention. DETAILED DESCRIPTION

[0030] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0031] Figure 1 A flow chart of a solute transport calculation method based on a hydrogeological model provided in an embodiment of the present invention is shown in FIG. Figure 1 As shown, a solute transport calculation method based on a hydrogeological model provided by an embodiment of the present invention includes:

[0032] Step 1, superimposing the hydrogeological grid in the grid model obtained in the solute transport simulation with the hydrogeological model in three-dimensional space;

[0033] Step 2: Each hydrogeological grid automatically adjusts the top and bottom plate heights and layered structures of the corresponding cells according to the top and bottom plate elevations and layered structures of the corresponding strata in the superimposed hydrogeological model, and assigns parameters to the corresponding hydrogeological grids using the hydrogeological parameters in the hydrogeological model. In this way, the assignment and structural adjustment of the hydrogeological grids are completed.

[0034] Step 3, assigning values ​​to other simulation data of the hydrogeological grid in the grid model according to the actual conditions of the simulation area, wherein the other simulation data are parameters required for the migration simulation operation and not in the hydrogeological model;

[0035] Step 4: extract the parameters required for the transport simulation operation from the grid model, and perform groundwater flow and solute transport simulation calculations in groundwater through the transport simulation calculation module.

[0036] The schematic diagram of the hydrogeological model is as follows: Figure 2 As shown, the grid model is composed of multiple hydrogeological grids, and the schematic diagram of the grid model with adjusted top and bottom plate heights of the cells and layered structure is shown in Figure 3 shown.

[0037] Among them, other simulation data include simulation initial condition data, historical simulation data and hydrological and meteorological data. When there is no simulation initial condition data, the historical simulation data will be used as the simulation initial condition data for this simulation. The hydrological and meteorological data includes hydrological and meteorological information such as rainfall, evaporation, and seepage. In this way, data input can be combined with the actual situation of the simulation area.

[0038] The present invention fully considers the problems that many solute transport simulation software have various and complicated data input items during the simulation process, which require high model understanding ability of researchers and the manual processing is cumbersome and non-intuitive. It uses an existing hydrogeological model to automatically generate the grid structure and related boundary data input required for the simulation, which not only improves the work efficiency of solute transport simulation but also ensures the accuracy and consistency of hydrogeological data.

[0039] In some embodiments, a hydrogeological model is established by collecting hydrogeological data; or an already constructed hydrogeological model file is directly imported, wherein the hydrogeological model includes various hydrogeological parameters required for simulation.

[0040] In some embodiments, step 1 further comprises:

[0041] Initialize the solute transport simulation, including: setting the measurement units of the basic simulation data and the default parameters of the solver equations before starting the simulation calculation;

[0042] Generate a grid model based on the grid accuracy and simulation range in the basic data of the simulation.

[0043] This includes basic information such as the spatial area of ​​the simulation, the simulation time, the mesh accuracy of the simulation, the units of measurement for basic data such as length and time, and the default coefficients of the solver equations.

[0044] In some embodiments, the method further includes: after performing the three-dimensional space superposition, using the hydrogeological model to automatically input boundary conditions, specifically including:

[0045] Obtain a three-dimensional view after three-dimensional spatial superposition, perform spatial superposition analysis based on the hydrogeological model, automatically identify the environmental information in the hydrogeological model, superimpose the environmental information in the hydrogeological model with the grid model in the three-dimensional view, and assign corresponding boundary conditions to the corresponding hydrogeological grid based on the environmental information.

[0046] The environmental information includes information about lakes, rivers, oceans and groundwater bodies.

[0047] In the above embodiment, by automatically identifying information such as lakes, rivers, oceans, groundwater bodies, etc. in the hydrogeological model, it is superimposed and displayed with the grid model in a three-dimensional view, and the corresponding simulation unit grids are automatically assigned corresponding boundary conditions.

[0048] Boundary conditions are used to describe the interaction between the groundwater system and the external environment. By setting these boundary conditions, the groundwater flow and solute transport processes can be simulated more accurately.

[0049] In some embodiments, the method further includes: after adjusting the top and bottom plate heights of the corresponding cells, the layered structure and the parameter assignment of the corresponding hydrogeological grid, based on the actual hydrogeological conditions, manually adjusting the grid parameters and grid structure of the hydrogeological grid through human-computer interaction, thereby realizing manual intervention for grid correction.

[0050] In some embodiments, in step 4, MODFLOW and MT3DMS are used as transport simulation calculation modules, and the input parameters required by MODFLOW and MT3DMS are extracted from the grid model to perform simulation calculations of groundwater flow and solute transport in groundwater.

[0051] Among them, MODFLOW is a finite difference model used to simulate groundwater flow. MT3DMS is an extension module of MODFLOW, focusing on the simulation of solute migration. It is suitable for the migration of chemical substances in groundwater. MT3DMS expands the function of solute migration on the basis of MODFLOW and can handle complex solute migration processes.

[0052] In this step, all cells in the grid model, i.e. the hydrogeological grid, are traversed, and information such as the cell elevation, spatial coordinates, hydrogeological parameters, and boundary conditions are organized into an input format recognizable by the calculation module.

[0053] In some embodiments, the method also includes: in step 4, the migration simulation calculation results are discrete point head and concentration data with time series, the spatial position information of the discrete points corresponds to the cell space information in the grid model, and after the calculation results are read, the calculation data is associated with the cells in the grid model and stored in the database.

[0054] After the data required for the simulation is input, the system will perform integrity and feasibility checks. If the checks pass, the MODFLOW and MT3DMS calculation programs will be called to perform simulation calculations. The specific calculation method is not the focus of this article and will not be elaborated in detail.

[0055] In some embodiments, the method further includes visually displaying the calculation results, specifically including:

[0056] The calculated data is extracted from the database, and the center points of the cells in the grid model are used as the spatial coordinates of the discrete points. The simulation results are converted into three-dimensional models, plane diagrams, and solute transport animations using spatial interpolation, and corresponding statistical charts are generated. The export of simulation results is also supported.

[0057] Figure 4 It is a schematic diagram of the overall process of the present invention, such as Figure 4 As shown, the overall process of the present invention is as follows:

[0058] (1) Construct or import a hydrogeological model;

[0059] (2) Initialize the solute transport simulation, set the measurement units of the basic simulation data and the default parameters of the solver equations, and generate a grid model;

[0060] (3) superimposing the hydrogeological grid in the grid model obtained from the solute transport simulation with the hydrogeological model in three dimensions, automatically assigning values ​​to the hydrogeological grid in the grid model and adjusting its structure;

[0061] (4) Automatic input of boundary conditions using hydrogeological models;

[0062] (5) Based on the actual hydrogeological conditions, the grid parameters and grid structure of the hydrogeological grid are manually adjusted through human-computer interaction;

[0063] (6) assigning values ​​to other simulation data of the hydrogeological grid in the grid model according to the actual conditions of the simulation area;

[0064] (7) Using MODFLOW and MT3DMS as transport simulation calculation modules, the input parameters required by MODFLOW and MT3DMS were extracted from the grid model to simulate the groundwater flow and solute transport in groundwater;

[0065] (8) If MODFLOW and MT3DMS cannot be executed normally, return to step (5) and perform manual intervention adjustment again;

[0066] (8) Read and analyze the calculation results;

[0067] (9) Visualization and export of results.

[0068] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims

1. A solute transport calculation method based on a hydrogeological model, characterized in that: The method comprises: Step 1, superimposing the hydrogeological grid in the grid model obtained in the solute transport simulation with the hydrogeological model in three-dimensional space; Step 2, each hydrogeological grid automatically adjusts the top and bottom plate heights and layered structures of the corresponding cells according to the top and bottom plate elevations and layered structures of the corresponding strata in the superimposed hydrogeological model, and assigns parameters to the corresponding hydrogeological grids using the hydrogeological parameters in the hydrogeological model; Step 3, assigning values ​​to other simulation data of the hydrogeological grid in the grid model according to the actual conditions of the simulation area, wherein the other simulation data are parameters required for the migration simulation operation and not in the hydrogeological model; Step 4: extract the parameters required for the transport simulation operation from the grid model, and perform groundwater flow and solute transport simulation calculations in groundwater through the transport simulation calculation module.

2. The solute transport calculation method based on the hydrogeological model according to claim 1 is characterized in that: In step 1, a hydrogeological model is established by collecting hydrogeological data; or an already constructed hydrogeological model file is directly imported.

3. The solute transport calculation method based on the hydrogeological model according to claim 1 is characterized in that: Step 1 also includes: Initialize the solute transport simulation, including: setting the measurement units of the basic simulation data and the default parameters of the solver equations before starting the simulation calculation; Generate a grid model based on the grid accuracy and simulation range in the basic data of the simulation.

4. The solute transport calculation method based on the hydrogeological model according to claim 1 is characterized in that: The method further includes: after performing the three-dimensional space superposition, using the hydrogeological model to automatically input boundary conditions, specifically including: Obtain a three-dimensional view after three-dimensional spatial superposition, perform spatial superposition analysis based on the hydrogeological model, automatically identify the environmental information in the hydrogeological model, superimpose the environmental information in the hydrogeological model with the grid model in the three-dimensional view, and assign corresponding boundary conditions to the corresponding hydrogeological grid based on the environmental information.

5. The solute transport calculation method based on the hydrogeological model according to claim 4 is characterized in that: The environmental information includes information about lakes, rivers, oceans and groundwater bodies.

6. The solute transport calculation method based on the hydrogeological model according to claim 1 is characterized in that: The method further includes: after adjusting the top and bottom plate heights and layered structures of the corresponding cells and assigning parameters to the corresponding hydrogeological grids, artificially adjusting the grid parameters and grid structure of the hydrogeological grids through human-computer interaction based on the actual hydrogeological conditions.

7. The solute transport calculation method based on the hydrogeological model according to claim 1 is characterized in that: Other simulation data include simulation initial condition data, historical simulation data and hydrological and meteorological data. When there is no simulation initial condition data, the historical simulation data is used as the simulation initial condition data for this simulation. The hydrological and meteorological data includes hydrological and meteorological information such as rainfall, evaporation, and seepage.

8. The solute transport calculation method based on the hydrogeological model according to claim 1 is characterized in that: In step 4, MODFLOW and MT3DMS are used as transport simulation calculation modules, and the input parameters required by MODFLOW and MT3DMS are extracted from the grid model to simulate the groundwater flow and the transport of solutes in groundwater.

9. The solute transport calculation method based on the hydrogeological model according to claim 1, characterized in that: The method also includes: in step 4, the migration simulation calculation results are discrete point head and concentration data with time series, the spatial position information of the discrete points corresponds to the cell space information in the grid model, and after the calculation results are read, the calculation data is associated with the cells in the grid model and stored in the database.

10. The solute transport calculation method based on the hydrogeological model according to claim 9, characterized in that: The method further includes visually displaying the calculation results, specifically including: The calculated data is extracted from the database, and the center points of the cells in the grid model are used as the spatial coordinates of the discrete points. The simulation results are converted into three-dimensional models, plane diagrams, and solute transport animations using spatial interpolation, and corresponding statistical charts are generated.

Images