CN115329607A - System and method for evaluating underground water pollution - Google Patents
System and method for evaluating underground water pollution Download PDFInfo
- Publication number
- CN115329607A CN115329607A CN202211256594.3A CN202211256594A CN115329607A CN 115329607 A CN115329607 A CN 115329607A CN 202211256594 A CN202211256594 A CN 202211256594A CN 115329607 A CN115329607 A CN 115329607A
- Authority
- CN
- China
- Prior art keywords
- groundwater
- boundary
- numerical
- concentration
- calculation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000003911 water pollution Methods 0.000 title claims abstract description 22
- 238000004088 simulation Methods 0.000 claims abstract description 80
- 238000005070 sampling Methods 0.000 claims abstract description 49
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 31
- 238000011156 evaluation Methods 0.000 claims abstract description 21
- 238000009792 diffusion process Methods 0.000 claims abstract description 18
- 239000003245 coal Substances 0.000 claims abstract description 14
- 238000005516 engineering process Methods 0.000 claims abstract description 10
- 238000003895 groundwater pollution Methods 0.000 claims abstract description 10
- 239000003673 groundwater Substances 0.000 claims description 70
- 238000004364 calculation method Methods 0.000 claims description 68
- 238000005065 mining Methods 0.000 claims description 23
- 238000012544 monitoring process Methods 0.000 claims description 16
- 238000011160 research Methods 0.000 claims description 13
- 230000009189 diving Effects 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 239000000470 constituent Substances 0.000 claims description 8
- 238000011109 contamination Methods 0.000 claims description 7
- 239000003814 drug Substances 0.000 claims description 7
- 238000005315 distribution function Methods 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 6
- 239000011148 porous material Substances 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 239000004480 active ingredient Substances 0.000 claims description 3
- 239000008346 aqueous phase Substances 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000013480 data collection Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000008595 infiltration Effects 0.000 claims description 3
- 238000001764 infiltration Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 238000004080 punching Methods 0.000 claims description 3
- 239000002349 well water Substances 0.000 claims description 3
- 235000020681 well water Nutrition 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 2
- 239000012528 membrane Substances 0.000 claims description 2
- 238000005325 percolation Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 238000001223 reverse osmosis Methods 0.000 abstract description 2
- 230000007547 defect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/18—Water
- G01N33/182—Specific anions in water
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Geophysics And Detection Of Objects (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
The invention relates to the field of data processing of prediction evaluation, and particularly discloses an evaluation system and method for underground water pollution. According to the invention, the geological condition information and hydrogeological condition information of the closed pit coal mine area can be acquired through the arrangement of the data acquisition module, and SO is used 4 2‑ The diffusion of (2) represents the diffusion of groundwater pollution caused by the reverse osmosis of accumulated water in a mine pit, and the mass conservation and SO are utilized 4 2‑ Establishing a numerical model according to diffusion rules, arranging grids in a three-dimensional region to be researched, solving the determined position coordinates of a pollution source by combining known boundary conditions and initial conditions, sampling at a calibrated target, analyzing the credibility of a numerical simulation result, and solving the problem that the existing underground water pollution evaluation system does not combine numerical valuesThe simulation technology is used for inversing and tracking the position of the coal mine series layer polluted well, and the accuracy of treatment of the coal mine series layer polluted well is improved.
Description
Technical Field
The invention relates to the field of data processing of prediction evaluation, in particular to an evaluation system and method for groundwater pollution.
Background
Because the groundwater is below the earth's surface and exists in the gap of rock, therefore groundwater is difficult to discover after being polluted, the groundwater can not be excavated all the earth's surface to look for pollution sources after being polluted, only a plurality of sampling wells can be excavated in the earth's surface, groundwater pollution and pollution sources can be evaluated through the groundwater data of the sampling wells, the identification of groundwater pollution sources through the numerical simulation method combined with groundwater sampling monitoring and inversion becomes very important, and especially, the assessment of groundwater pollution through the numerical simulation means in coal mines is a technology to be developed urgently.
In the literature, "research progress and development direction of underground water pollution prevention and control technology of coal mine site", it is mentioned that a pit is closed and production is stopped after long-term large-scale mining, underground water pollution is caused by stopping drainage of pit water after pit closing and production stopping, but complex mining can be performed according to the position of a coal seam in the excavation process of a coal mine, the trend and distribution of the pit are complicated, so that the position of pit water accumulation of polluted underground water is difficult to find, a string layer polluted well causing underground water pollution is found, effective plugging and feeding of physicochemical materials or medicaments are important means for preventing and controlling the underground water pollution of the mine, but the existing underground water pollution assessment system does not use a numerical simulation means to utilize underground water sampling monitoring data to reversely identify the string layer polluted well position of the polluted underground water, and the string layer polluted well position causing underground water pollution can be effectively plugged or the string layer polluted well causing underground water pollution can be prevented from being aggravated after the string layer polluted well position causing underground water pollution is determined by utilizing the numerical simulation means and the underground water sampling monitoring data to reversely identify.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the invention provides an underground water pollution evaluation system and method, and aims to overcome the defect that the positions of series-layer polluted wells polluting underground water are inversely identified by means of underground water sampling monitoring data combined with a numerical simulation method in the conventional underground water pollution evaluation system.
The invention solves the technical problems through the following technical scheme, and the device comprises a data acquisition module, a numerical inversion tracking module and a pollution evaluation module, wherein the data acquisition module is used for transmitting geological and hydrogeological conditions of a research mining area and underground water monitoring data of a sampling well to the numerical inversion tracking module;
the method for tracking the underground water pollution source of the closed-pit coal mine area by the underground water pollution evaluation system through numerical simulation technology inversion comprises the following steps:
the method comprises the following steps: the data acquisition module acquires the underground water quality and water level information of a sampling well, the geological conditions and hydrogeological conditions of a mining area and the range of a cross-layer polluted well, and transmits data to the numerical inversion tracking module;
step two: the numerical inversion tracking module carries out grid division on a researched area, locally encrypts the area in the range of the series layer contaminated well, and is based on SO according to data provided by the data acquisition module 4 2- Diffusion of (2) to establish multiple phasesThe flow numerical simulation model sets numerical simulation initial conditions and boundary conditions according to data provided by the data acquisition module, defines the range of the series layer polluted well and solves the three-dimensional coordinates of the pollution source;
step three: obtaining three-dimensional coordinates according to a numerical simulation technology, punching and sampling at the coordinate position, and verifying SO in the groundwater sample at the position of the pollution source 4 2- Whether the concentration of (b) is close to the calculated concentration obtained by numerical simulation calculation;
step four: performing numerical simulation on different string layer polluted wells by the same method, and collecting numerical simulation results and sampling results;
step five: analyzing the reliability of the calculation result of the numerical simulation, and calculating to obtain SO by using the numerical simulation 4 2- Concentration minus actual SO in the sampled groundwater 4 2- Absolute value of concentration, divided by SO in the actual sampled groundwater 4 2- The obtained percentage is used as the basis of credibility evaluation, the obtained percentage value which is more than 80 percent is recorded as a credible sample group, and the stratum polluted well corresponding to the credible sample group is effectively plugged or a medicament is put in for treatment.
As a further scheme of the present invention, the real-time monitoring data of the mining area geological condition, the hydrogeological condition and the groundwater sampling well collected by the data collection module includes:
three-dimensional zone boundaries of the actual study site: the method is used for defining a three-dimensional calculation area of a numerical simulation calculation flow field, setting a calculated space step length L according to the actual field size, and dividing the actual research field size into
In the formula, aL, bL and cL are recorded as the length, width and depth of an actual research site, and a coordinate set in a calculation area is recorded as
；
And (3) the stratum properties and classification of the mining area to be evaluated: the device is used for setting numerical values to calculate the type of the diving aquifer and setting the property and the parameter of the diving aquifer according to most of the types of the diving aquifers;
surface type: setting the upper boundary of the numerical simulation calculation as a water exchange boundary according to the field openness and rainwater infiltration characteristics
；
Boundary formation properties of the underlying groundwater: the bottom boundary property and parameter used for setting the numerical simulation calculation area are set as the zero flux boundary according to the characteristics of the bottom micro-penetration
；
And (3) researching underground water flow and water pressure information around the mining area: setting the types and parameters of the upper boundary, the lower boundary, the left boundary and the right boundary of the calculation area, setting the upper boundary and the lower boundary of the simulation area to be watertight zero-flux boundaries, and respectively recording the upper boundary and the lower boundary as water-impervious zero-flux boundaries
And
the left and right boundaries are set as boundaries of known pressure and are noted as the left and right boundaries, respectively
And
；
setting the sampling time interval of the monitoring data: as a basis for dividing monitoring time of numerical simulation calculation;
sampling well information: the numerical simulation system is used for numbering the sampling wells and carrying out numerical calculation in the numerical simulation calculation;
sampling SO in well water quality information 4 2- Concentration: from the desired position SO of the pollution source by means of sampling 4 2- Extracting 300 groups of data from the concentration range, substituting the data into the multiphase flow numerical simulation model, inverting the model, and taking 100 groups of close sampling wells SO 4 2- And (4) checking the concentration data, and performing statistical calculation on the data with the accuracy reaching the preset requirement.
As a further aspect of the present invention, the numerical modeling information of the numerical tracking inversion module includes:
initial conditions were as follows: initial pressure distribution, initial SO including simulation zone 4 2- And an aqueous phase and an initial SO 4 2- The initial conditions are as follows:
in the formula:
boundary conditions: setting known flow boundary and pressure boundary according to hydrogeological conditions, setting source and sink terms and SO according to the rainstorm intensity of mining area and the information of peripheral rivers 4 2- Concentration boundary, boundary conditions are as follows:
in the formula:
record as SO in groundwater 4 2- The concentration of (c);
is recorded as the actual average flow velocity of the groundwater,
is marked as SO 4 2- The degree of saturation of (a) is,
record as SO in groundwater 4 2- N is expressed as the normal direction of the boundary,
noted as a known pressure distribution function on the upper boundary,
the known pressure function noted as the left boundary,
the known pressure function noted as the right boundary,
is denoted as upper boundary SO 4 2- Is determined by the known concentration function of (a),
SO marked as left boundary 4 2- With the knowledge of the function of the concentration,
SO marked as right boundary 4 2- A known concentration function;
time step length: setting a numerical simulation calculation time step length T according to the sampling time interval of the sampling well, and selecting calculation moments T =10T,11T,12T,13T and 14T which need to be subjected to simulation calculation;
the multiphase flow numerical simulation calculates partial differential equation: establishing a multiphase flow partial differential equation according to the mass conservation law, wherein the multiphase flow partial differential equation is as follows:
in the formula: s is the kind of the constituent, s =1 is the water as the constituent, s =2 is the SO as the constituent 4 2- ，
Is recorded as studying the porosity of the formation in the simulated zone,
record as SO in groundwater 4 2- (ii) the total concentration of (a) in (b),
as a result of the density of the constituent materials,
record as SO in groundwater 4 2- The concentration of (a) in (b),
record as SO in groundwater 4 2- The degree of saturation of (a) is,
record as SO in groundwater 4 2- Diffusion coefficient tensor of (a);
in the formula: d is denoted as SO 4 2- The coefficient of molecular diffusion in groundwater,
in the form of the distortion coefficient of the formation,
is marked as the function of the kronecker,
and
are respectively marked as SO 4 2- Transverse and longitudinal dispersivity in groundwater,
is marked as SO 4 2- The speed of the seepage of (a) is,
and
respectively represent SO 4 2- The seepage partial velocity in the i direction and in the j direction;
in the formula:
is marked as SO 4 2- The relative permeability of the porous material to the porous material,
is marked as SO 4 2- The intrinsic permeability of the membrane is such that,
in the notation of the viscosity of the groundwater,
recording as the pressure of the groundwater, g as the acceleration of gravity,
is recorded as the density of the groundwater.
As a further scheme of the invention, the numerical inversion tracking module is connected with a display screen, the display screen is used for displaying numerical simulation calculation results, displaying the calculation results on a three-dimensional map of a mining area, calibrating three-dimensional coordinates of a stratum contamination well with highest credibility, and displaying SO at the position 4 2- The concentration of (c).
As a further scheme of the present invention, the contamination evaluation module evaluates the reliability of the calculation result of the numerical inversion tracking module, and the reliability calculation formula is:
in the formula:
the data is recorded as the reliability of the data,
is recorded as SO obtained by numerical calculation 4 2- The concentration of the active ingredients in the mixture is,
record as the actual SO sampled 4 2- And (4) concentration.
Compared with the prior art, the invention has the following advantages: the system and the method for evaluating the groundwater pollution can acquire geological condition information and hydrogeological condition information of a closed pit coal mine area through the arrangement of a data acquisition module, and use SO 4 2- The diffusion of (2) represents the diffusion of groundwater pollution caused by the reverse osmosis of accumulated water in a mine pit, and the mass conservation and SO are utilized 4 2- The method comprises the steps of establishing a numerical model according to a diffusion rule, arranging grids for a researched three-dimensional area, solving a determined position coordinate of a pollution source by combining known boundary conditions and initial conditions, sampling at a calibrated target after a calculation result is obtained, analyzing the reliability of a numerical simulation result, and effectively plugging or feeding a medicament for a series layer polluted well of an area with the reliability meeting requirements, so that the defect that the position of the coal mine series layer polluted well is not inversely tracked by combining a numerical simulation technology in the existing underground water pollution assessment system is overcome, the manpower input of a large number of well drilling and sampling is reduced, and the accuracy of treatment of the coal mine series layer polluted well is improved.
Drawings
Fig. 1 is an overall structural view of the present invention.
Detailed Description
The following examples are given for the detailed implementation and the specific operation procedures, but the scope of the present invention is not limited to the following examples.
As shown in fig. 1, the present embodiment provides a technical solution: an assessment system for groundwater pollution comprises a data acquisition module, a numerical inversion tracking module and a pollution assessment module, wherein the data acquisition module is used for transmitting geological and hydrogeological conditions of a research mining area and sampling well groundwater monitoring data to the numerical inversion tracking module, the numerical inversion tracking module sets the hydrogeological conditions, calculates a flow field boundary and known boundary conditions according to acquired data information, and establishes a multiphase flow numerical simulation model, and the numerical inversion tracking module is connected with the pollution assessment module;
the method for tracking the underground water pollution source of the closed-pit coal mine area by the underground water pollution evaluation system through numerical simulation technology inversion comprises the following steps:
the method comprises the following steps: the data acquisition module acquires the underground water quality and water level information of the sampling well, the geological conditions and hydrogeological conditions of the mining area and the range of the series layer polluted well, and transmits the data to the numerical inversion tracking module;
step two: the numerical inversion tracking module carries out grid division on a researched area, locally encrypts the area in the range of the series layer contaminated well, and is based on SO according to data provided by the data acquisition module 4 2- The method comprises the steps of (1) establishing a multiphase flow numerical simulation model by diffusion, setting numerical simulation initial conditions and boundary conditions according to data provided by a data acquisition module, defining the range of a series layer polluted well, and solving the three-dimensional coordinates of a pollution source;
step three: obtaining three-dimensional coordinates according to a numerical simulation technology, punching and sampling at the coordinate position, and verifying SO in the groundwater sample at the position of the pollution source 4 2- Whether the concentration of (a) is close to the calculated concentration obtained by numerical simulation calculation;
step four: performing numerical simulation on different string layer polluted wells by the same method, and collecting numerical simulation results and sampling results;
step five: analyzing the reliability of the calculation result of the numerical simulation, and calculating to obtain SO by using the numerical simulation 4 2- Concentration minus actual SO in the sampled groundwater 4 2- Absolute value of concentration, divided by SO in the actual sampled groundwater 4 2- The obtained percentage is used as the basis of credibility evaluation, the obtained percentage value which is more than 80 percent is recorded as a credible sample group, and the stratum polluted well corresponding to the credible sample group is effectively plugged or a medicament is put in for treatment.
Geological condition information and hydrogeological condition information of a closed pit coal mine area can be acquired through the arrangement of the data acquisition module, and SO is used 4 2- The diffusion represents the diffusion of groundwater pollution caused by the reverse seepage of the accumulated water in the pit, and the quality is utilizedConservation and SO 4 2- The method comprises the steps of establishing a numerical model according to a diffusion rule, arranging grids for a three-dimensional region to be researched, solving a determined position coordinate of a pollution source by combining known boundary conditions and initial conditions, sampling and analyzing the credibility of a numerical simulation result at a calibrated target after the calculation result is obtained, and effectively plugging or feeding a medicament into a stratum-series polluted well in a region with the credibility meeting requirements.
The mining area geological condition, hydrogeological condition and groundwater sampling well real-time monitoring data collected by the data collection module comprise:
three-dimensional zone boundaries of the actual research site: the method is used for defining a three-dimensional calculation area of a numerical simulation calculation flow field, setting a calculated space step length L according to the actual field size, and dividing the actual research field size into
In the formula, aL, bL and cL are recorded as the length, width and depth of an actual research site, and a coordinate set in a calculation area is recorded as
；
And (3) the stratum properties and classification of the mining area to be evaluated: the device is used for setting numerical values to calculate the type of the diving aquifer and setting the property and the parameter of the diving aquifer according to most of the types of the diving aquifers;
surface type: setting the upper boundary of the numerical simulation calculation, setting the upper boundary as a water exchange boundary according to the field openness and the rainwater infiltration characteristic, and recording the boundary as a water exchange boundary
；
Boundary formation properties of bottom groundwater: bottom for setting numerical simulation calculation areaBoundary properties and parameters, the bottom boundary of the investigation region was set to be a zero-flux boundary, recorded as
；
And (3) researching underground water flow and water pressure information around the mining area: setting the types and parameters of the upper boundary, the lower boundary, the left boundary and the right boundary of the calculation area, setting the upper boundary and the lower boundary of the simulation area to be watertight zero-flux boundaries, and respectively recording the upper boundary and the lower boundary as water-impervious zero-flux boundaries
And
the left and right boundaries are set as boundaries of known pressure and are noted as the left and right boundaries, respectively
And
；
setting a sampling time interval of monitoring data: as a basis for dividing monitoring time of numerical simulation calculation;
sampling well information: the numerical simulation system is used for numbering the sampling wells in numerical simulation calculation and carrying out numerical calculation;
SO in sampling well water quality information 4 2- Concentration: from the desired position SO of the pollution source by means of sampling 4 2- Extracting 300 groups of data from the concentration range, substituting the data into the multiphase flow numerical simulation model, inverting the model, and taking 100 groups of close sampling wells SO 4 2- And (4) checking the concentration data, and performing statistical calculation on the data with the accuracy reaching the preset requirement.
The numerical modeling information of the numerical tracking inversion module includes:
initial conditions: initial pressure distribution, initial SO including simulation zone 4 2- And an aqueous phase and an initial SO 4 2- The initial conditions are as follows:
in the formula:
boundary conditions: setting known flow boundary and pressure boundary according to hydrogeological conditions, setting source and sink terms and S according to the rainstorm intensity of mining area and the information of peripheral riversO 4 2- Concentration boundary, the boundary conditions are:
in the formula:
record as SO in groundwater 4 2- The concentration of (d);
is recorded as the actual average flow velocity of the groundwater,
is marked as SO 4 2- Is saturated withThe degree of the magnetic field is measured,
record as SO in groundwater 4 2- N is expressed as the normal direction of the boundary,
noted as a known pressure distribution function on the upper boundary,
the known pressure function noted as the left boundary,
the known pressure function noted as the right boundary,
is denoted as upper boundary SO 4 2- Is determined by the known concentration function of (a),
SO marked as left boundary 4 2- With the knowledge of the function of the concentration,
SO marked as right boundary 4 2- A known concentration function;
step length of time: setting a numerical simulation calculation time step length T according to the sampling time interval of the sampling well, and selecting calculation moments T =10T,11T,12T,13T and 14T which need to be subjected to simulation calculation;
and (3) performing numerical simulation calculation of partial differential equation of multiphase flow: establishing a multiphase flow partial differential equation according to a mass conservation law, wherein the multiphase flow partial differential equation is as follows:
in the formula: s is the type of the constituent, s =1 is the water, and s =2 is the SO 4 2- ，
Is recorded as studying the porosity of the formation in the simulated zone,
record as SO in groundwater 4 2- (ii) the total concentration of (a) in (b),
as a result of the density of the constituent materials,
record as SO in groundwater 4 2- The concentration of (a) in (b),
record as SO in groundwater 4 2- The degree of saturation of (a) is,
record as SO in groundwater 4 2- Diffusion coefficient tensor of (a);
in the formula: d is denoted as SO 4 2- The coefficient of molecular diffusion in groundwater,
in the form of the distortion coefficient of the formation,
is marked as the function of the kronecker,
and
are respectively marked as SO 4 2- Transverse and longitudinal dispersivity in groundwater,
is marked as SO 4 2- The speed of the seepage of (a) is,
and
respectively represent SO 4 2- The percolation partial velocity in the i direction and in the j direction;
in the formula:
is marked as SO 4 2- The relative permeability of the porous material to the porous material,
is marked as SO 4 2- The inherent permeability of (a) to (b),
and is recorded as the viscosity of the groundwater,
is recorded as the pressure of the groundwaterThe force, g, is recorded as the acceleration of gravity,
is recorded as the density of the groundwater.
The numerical inversion tracking module is connected with a display screen, the display screen is used for displaying numerical simulation calculation results, displaying the calculation results on a three-dimensional map of a mining area, calibrating three-dimensional coordinates of a string layer polluted well with highest credibility and displaying SO at the position 4 2- The concentration of (c).
The result of numerical value calculation can be displayed on the display screen through the setting of the display screen, and the position of the string layer polluted well is specifically displayed by combining a three-dimensional map on the display screen, so that researchers can conveniently analyze the position of the string layer polluted well and discuss the plugging mode and the dosing metering, and the underground water pollution can be effectively prevented and treated.
The pollution evaluation module evaluates the reliability of the calculation result of the numerical inversion tracking module, and the reliability calculation formula is as follows:
in the formula:
is recorded as the reliability of the data,
is recorded as SO obtained by numerical calculation 4 2- The concentration of the active ingredients in the mixture is,
record as the actual SO sampled 4 2- And (4) concentration.
The result of numerical calculation can be judged through the setting of the credibility formula, and then the stratum contamination well treatment can be implemented to the calculation coordinate with the credibility meeting the requirement, so that the waste of plugging materials and the thrown medicament caused by inaccurate calculation result of numerical simulation is avoided.
Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention, and that variations, modifications, substitutions and alterations can be made to the above embodiments by those of ordinary skill in the art within the scope of the present invention.
Claims (5)
1. An assessment system for groundwater pollution comprises a data acquisition module, a numerical inversion tracking module and a pollution assessment module, and is characterized in that the data acquisition module is used for transmitting geological and hydrogeological conditions of a research mining area and groundwater monitoring data of a sampling well to the numerical inversion tracking module, the numerical inversion tracking module sets the hydrogeological conditions, calculates flow field boundaries and known boundary conditions according to acquired data information, and establishes a multiphase flow numerical simulation model, and the numerical inversion tracking module is connected with the pollution assessment module; the numerical modeling information of the numerical inversion tracking module comprises:
initial conditions: initial pressure distribution, initial SO including simulation zone 4 2- And an aqueous phase and an initial SO 4 2- The initial conditions are as follows:
in the formula:
boundary conditions: setting known flow boundary and pressure boundary according to hydrogeological conditions, setting source and sink terms and SO according to the rainstorm intensity of mining area and the information of peripheral rivers 4 2- Concentration boundary, boundary conditions are as follows:
in the formula:
record as SO in groundwater 4 2- The concentration of (c);
is recorded as the actual average flow velocity of the groundwater,
is marked as SO 4 2- The degree of saturation of (a) is,
record as SO in groundwater 4 2- N is expressed as the normal direction of the boundary,
noted as a known pressure distribution function on the upper boundary,
the known pressure function noted as the left boundary,
noted as the known function of pressure at the right boundary,
is denoted as upper boundary SO 4 2- Is determined by the known concentration function of (a),
SO marked as left boundary 4 2- With the knowledge of the function of the concentration,
SO marked as right boundary 4 2- A known concentration function;
time step length: setting a numerical simulation calculation time step length T according to the sampling time interval of the sampling well, and selecting calculation moments T =10T,11T,12T,13T and 14T which need to be subjected to simulation calculation;
and (3) performing numerical simulation calculation of partial differential equation of multiphase flow: establishing a multiphase flow partial differential equation according to the mass conservation law, wherein the multiphase flow partial differential equation is as follows:
in the formula: s is the type of the constituent, s =1 is the water, and s =2 is the SO 4 2- ，
Is recorded as the porosity of the stratum in the research simulation area,
record as SO in groundwater 4 2- (ii) the total concentration of (a) in (b),
as a result of the density of the constituent materials,
record as SO in groundwater 4 2- The concentration of (a) in (b),
record as SO in groundwater 4 2- The degree of saturation of (a) is,
record as SO in groundwater 4 2- Diffusion coefficient tensor of (a);
in the formula: d is denoted as SO 4 2- The coefficient of molecular diffusion in groundwater,
in the form of the distortion coefficient of the formation,
is marked as the function of the kronecker,
and
are respectively marked as SO 4 2- Transverse and longitudinal dispersivity in groundwater,
is marked as SO 4 2- The speed of the seepage of (a) is,
and
respectively represent SO 4 2- The percolation partial velocity in the i direction and in the j direction;
in the formula:
is marked as SO 4 2- The relative permeability of the porous material to the porous material,
is recorded as SO 4 2- The intrinsic permeability of the membrane is such that,
in the notation of the viscosity of the groundwater,
recording as the pressure of the groundwater, g as the acceleration of gravity,
and is recorded as the density of the groundwater.
2. An evaluation system for groundwater contamination according to claim 1, wherein: the mining area geological condition, hydrogeological condition and groundwater sampling well real-time monitoring data collected by the data collection module comprise:
three-dimensional zone boundaries of the actual research site: the method is used for defining a three-dimensional calculation area of a numerical simulation calculation flow field, setting a calculated space step length L according to the actual field size, and dividing the actual research field size into
In the formula, aL, bL and cL are recorded as the length, width and depth of an actual research site, and a coordinate set in a calculation area is recorded as
；
And (3) the stratum properties and classification of the mining area to be evaluated: the device is used for setting numerical values to calculate the type of the diving aquifer and setting the property and the parameter of the diving aquifer according to most of the types of the diving aquifers;
surface type: setting the upper boundary of the numerical simulation calculation as a water exchange boundary according to the field openness and rainwater infiltration characteristics
；
Boundary formation properties of bottom groundwater: the bottom boundary property and parameter used for setting the numerical simulation calculation area are set as the zero flux boundary according to the characteristics of the bottom micro-penetration
；
And (3) researching underground water flow and water pressure information around the mining area: setting the types and parameters of the upper boundary, the lower boundary, the left boundary and the right boundary of the calculation area, setting the upper boundary and the lower boundary of the simulation area to be watertight zero-flux boundaries, and respectively recording the upper boundary and the lower boundary as water-impervious zero-flux boundaries
And
the left and right boundaries are set as boundaries of known pressure and are noted as the left and right boundaries, respectively
And
；
setting the sampling time interval of the monitoring data: as a basis for dividing monitoring time of numerical simulation calculation;
sampling well information: the numerical simulation system is used for numbering the sampling wells in numerical simulation calculation and carrying out numerical calculation;
SO in sampling well water quality information 4 2- Concentration: from the desired position SO of the pollution source by means of sampling 4 2- Extracting 300 groups of data from the concentration range, substituting the data into the multiphase flow numerical simulation model, inverting the model, and taking 100 groups of close sampling wells SO 4 2- And (4) checking the data of the concentration, and performing statistical calculation on the data with the accuracy reaching the preset requirement.
3. An evaluation system for groundwater contamination according to claim 1, wherein: the numerical inversion tracking module is provided with a display screen, the display screen is used for displaying numerical simulation calculation results and displaying the calculation results on a three-dimensional map of a mining area, calibrating three-dimensional coordinates of a string layer polluted well with highest credibility and displaying SO at the position 4 2- The concentration of (c).
4. An evaluation system for groundwater contamination according to claim 1, wherein: the pollution evaluation module evaluates the reliability of the calculation result of the numerical inversion tracking module, and the reliability calculation formula is as follows:
in the formula:
the data is recorded as the reliability of the data,
is recorded as SO obtained by numerical calculation 4 2- The concentration of the active ingredients in the mixture is,
record as the actual SO sampled 4 2- And (4) concentration.
5. An evaluation system for groundwater contamination according to claim 1, wherein: the method for tracking the underground water pollution source of the closed-pit coal mine area by the underground water pollution evaluation system through numerical simulation technology is as follows:
the method comprises the following steps: the data acquisition module acquires the underground water quality and water level information of the sampling well, the geological conditions and hydrogeological conditions of the mining area and the range of the series layer polluted well, and transmits the data to the numerical inversion tracking module;
step two: the numerical inversion tracking module carries out grid division on a researched area, locally encrypts the area in the range of the series layer contaminated well, and is based on SO according to data provided by the data acquisition module 4 2- The method comprises the steps of (1) establishing a multiphase flow numerical simulation model by diffusion, setting numerical simulation initial conditions and boundary conditions according to data provided by a data acquisition module, defining the range of a series layer polluted well, and solving the three-dimensional coordinates of a pollution source;
step three: obtaining three-dimensional coordinates according to a numerical simulation technology, punching and sampling at the coordinate position, and verifying SO in the groundwater sample at the position of the pollution source 4 2- Whether the concentration of (a) is close to the calculated concentration obtained by numerical simulation calculation;
step four: carrying out numerical simulation on different string layer polluted wells by adopting the same method, and collecting numerical simulation results and sampling results;
step five: analyzing the reliability of the calculation result of the numerical simulation, and calculating to obtain SO by using the numerical simulation 4 2- Concentration minus actual SO in the sampled groundwater 4 2- Absolute value of concentration, divided by SO in the actual sampled groundwater 4 2- The obtained percentage is used as the basis of credibility evaluation, the obtained percentage value which is more than 80 percent is recorded as a credible sample group, and the stratum polluted well corresponding to the credible sample group is effectively plugged or a medicament is put in for treatment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211256594.3A CN115329607B (en) | 2022-10-14 | 2022-10-14 | System and method for evaluating underground water pollution |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211256594.3A CN115329607B (en) | 2022-10-14 | 2022-10-14 | System and method for evaluating underground water pollution |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115329607A true CN115329607A (en) | 2022-11-11 |
| CN115329607B CN115329607B (en) | 2023-02-03 |
Family
ID=83914649
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211256594.3A Active CN115329607B (en) | 2022-10-14 | 2022-10-14 | System and method for evaluating underground water pollution |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115329607B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116611274A (en) * | 2023-07-21 | 2023-08-18 | 中南大学 | Visual numerical simulation method for groundwater pollution migration |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109190280A (en) * | 2018-09-18 | 2019-01-11 | 东北农业大学 | A kind of pollution source of groundwater inverting recognition methods based on core extreme learning machine alternative model |
| US20200252283A1 (en) * | 2017-10-16 | 2020-08-06 | Chinese Research Academy Of Environmental Sciences | Multi-objective optimization method for groundwater pollution monitoring network |
| CN112149871A (en) * | 2020-08-25 | 2020-12-29 | 河海大学 | Pollutant point source analysis method based on combination of GIS (geographic information System) space statistics and random simulation |
| CN112697849A (en) * | 2020-11-20 | 2021-04-23 | 中科盛联(北京)新材料有限公司 | Dynamic pollution source positioning method |
| CN113239598A (en) * | 2021-06-08 | 2021-08-10 | 中国环境科学研究院 | Underground water pollution source space comprehensive identification method applying numerical simulation |
| WO2021208393A1 (en) * | 2020-04-15 | 2021-10-21 | 北京工业大学 | Inversion estimation method for air pollutant emission inventory |
| CN114254572A (en) * | 2021-12-16 | 2022-03-29 | 西北工业大学太仓长三角研究院 | Aero-engine compressor flow field performance prediction method and system considering pollutant deposition |
| CN114386329A (en) * | 2022-01-11 | 2022-04-22 | 南京大学 | Multi-objective optimization method for DNAPL (deoxyribonucleic acid-Polypropylene) polluted site remediation under uncertain conditions |
| CN114492164A (en) * | 2021-12-24 | 2022-05-13 | 吉林大学 | Organic pollutant migration numerical model substitution method based on multi-core extreme learning machine |
| CN114707395A (en) * | 2021-12-24 | 2022-07-05 | 吉林大学 | Groundwater organic pollution source inversion method based on heuristic homotopy algorithm |
| CN114861502A (en) * | 2022-05-27 | 2022-08-05 | 北京林业大学 | Safe drinking water area determination method based on three-dimensional dynamic groundwater pollution simulation of Modflow model |
-
2022
- 2022-10-14 CN CN202211256594.3A patent/CN115329607B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20200252283A1 (en) * | 2017-10-16 | 2020-08-06 | Chinese Research Academy Of Environmental Sciences | Multi-objective optimization method for groundwater pollution monitoring network |
| CN109190280A (en) * | 2018-09-18 | 2019-01-11 | 东北农业大学 | A kind of pollution source of groundwater inverting recognition methods based on core extreme learning machine alternative model |
| WO2021208393A1 (en) * | 2020-04-15 | 2021-10-21 | 北京工业大学 | Inversion estimation method for air pollutant emission inventory |
| CN112149871A (en) * | 2020-08-25 | 2020-12-29 | 河海大学 | Pollutant point source analysis method based on combination of GIS (geographic information System) space statistics and random simulation |
| CN112697849A (en) * | 2020-11-20 | 2021-04-23 | 中科盛联(北京)新材料有限公司 | Dynamic pollution source positioning method |
| CN113239598A (en) * | 2021-06-08 | 2021-08-10 | 中国环境科学研究院 | Underground water pollution source space comprehensive identification method applying numerical simulation |
| CN114254572A (en) * | 2021-12-16 | 2022-03-29 | 西北工业大学太仓长三角研究院 | Aero-engine compressor flow field performance prediction method and system considering pollutant deposition |
| CN114492164A (en) * | 2021-12-24 | 2022-05-13 | 吉林大学 | Organic pollutant migration numerical model substitution method based on multi-core extreme learning machine |
| CN114707395A (en) * | 2021-12-24 | 2022-07-05 | 吉林大学 | Groundwater organic pollution source inversion method based on heuristic homotopy algorithm |
| CN114386329A (en) * | 2022-01-11 | 2022-04-22 | 南京大学 | Multi-objective optimization method for DNAPL (deoxyribonucleic acid-Polypropylene) polluted site remediation under uncertain conditions |
| CN114861502A (en) * | 2022-05-27 | 2022-08-05 | 北京林业大学 | Safe drinking water area determination method based on three-dimensional dynamic groundwater pollution simulation of Modflow model |
Non-Patent Citations (5)
| Title |
|---|
| GAO QI等: "Source Inversion of a Sudden Pollution Accident Occurred in Plain River Network Based on Pattern Search Method", 《IEEE XPLORE》 * |
| 侯泽宇: "地下水DNAPLs污染源反演识别及不确定性分析研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
| 杨之珍: "三维非稳定流数学模型在模拟废弃煤矿区地下水串层污染中的应用研究", 《能源与环境》 * |
| 杨佳立: "地下水水源保护区边界污染源防控的V-M数值模拟探析", 《科技广场》 * |
| 闫雪嫚: "基于贝叶斯理论的地下水DNAPLs污染源反演识别研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116611274A (en) * | 2023-07-21 | 2023-08-18 | 中南大学 | Visual numerical simulation method for groundwater pollution migration |
| CN116611274B (en) * | 2023-07-21 | 2023-09-29 | 中南大学 | Visual numerical simulation method for groundwater pollution migration |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115329607B (en) | 2023-02-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Desaulniers et al. | Origin, age and movement of pore water in argillaceous Quaternary deposits at four sites in southwestern Ontario | |
| Banks et al. | Assessing spatial and temporal connectivity between surface water and groundwater in a regional catchment: Implications for regional scale water quantity and quality | |
| Sauter | Quantification and forecasting of regional groundwater flow and transport in a karst aquifer (Gallusquelle, Malm, SW. Germany) | |
| Berhane et al. | Water leakage investigation of micro-dam reservoirs in Mesozoic sedimentary sequences in Northern Ethiopia | |
| Scheiber et al. | Hydrochemical apportioning of irrigation groundwater sources in an alluvial aquifer | |
| CN115329607B (en) | System and method for evaluating underground water pollution | |
| Bonanno et al. | Flow directions of stream‐groundwater exchange in a headwater catchment during the hydrologic year | |
| Ronchetti et al. | Tracer test to assess flow and transport parameters of an earth slide: The Montecagno landslide case study (Italy) | |
| Li et al. | Groundwater circulation patterns in bedrock aquifers from a pre-selected area of high-level radioactive waste repository based on two-dimensional numerical simulation | |
| Rushton et al. | Significance of hydraulic head gradients within horizontal wells in unconfined aquifers of limited saturated thickness | |
| Bourke et al. | A hydrological framework for persistent pools along non-perennial rivers | |
| Vorlicek et al. | Quantitative hydrogeological studies of the Treviso alluvial plain, NE Italy | |
| Michalski et al. | Characterization of transmissive fractures by simple tracing of in‐well flow | |
| Hauwert | Groundwater flow and recharge within the Barton Springs segment of the Edwards Aquifer, southern Travis and northern Hays counties, Texas | |
| Chiu et al. | The effect of hydrological conditions and bioactivities on the spatial and temporal variations of streambed hydraulic characteristics at the subtropical alpine catchment | |
| Rickel et al. | Seasonal shifts in surface water-groundwater connections in a ferricrete-impacted stream estimated from electrical resistivity | |
| CN109978273B (en) | Mine water inflow prediction method based on landform water control | |
| Belan | Characterizating a fractured rock aquifer with hydraulic testing at a contaminated municipal well using flexible liner methods and depth discrete monitoring | |
| YOKOYAMA et al. | A review of groundwater observation methods for slow-moving landslide | |
| Wilkes et al. | Characterisation of the hydrogeology of the Augustus River catchment, Western Australia | |
| Middleton | Aquifer–stream connectivity at various scales: Application of sediment–water interface temperature and vulnerability assessments of groundwater dependent streams | |
| Schalla et al. | Interim characterization report for the 300 Area process trenches | |
| Lee et al. | Investigating the origin and dynamics of salinity in a confined aquifer system in southeast Australia (Western Port Basin) | |
| Mahindawansha et al. | Evaluation of Hydraulic Exchange in the Hyporheic Zone: A Depth-Oriented Analysis Focusing on Upstream and Downstream Conditions | |
| Alberti et al. | Characterization of salt-water intrusion in the lower Esino Valley, Italy using a three-dimensional numerical model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |