CN117059185A - Simulation method, device, equipment and medium for migration process of target substance - Google Patents
Simulation method, device, equipment and medium for migration process of target substance Download PDFInfo
- Publication number
- CN117059185A CN117059185A CN202311309905.2A CN202311309905A CN117059185A CN 117059185 A CN117059185 A CN 117059185A CN 202311309905 A CN202311309905 A CN 202311309905A CN 117059185 A CN117059185 A CN 117059185A
- Authority
- CN
- China
- Prior art keywords
- heavy metal
- metal salt
- migration
- initial
- target substance
- 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
- 230000005012 migration Effects 0.000 title claims abstract description 469
- 238000013508 migration Methods 0.000 title claims abstract description 469
- 239000013076 target substance Substances 0.000 title claims abstract description 311
- 238000000034 method Methods 0.000 title claims abstract description 237
- 238000004088 simulation Methods 0.000 title claims abstract description 144
- 230000008569 process Effects 0.000 title claims abstract description 139
- 229910001385 heavy metal Inorganic materials 0.000 claims description 454
- 150000003839 salts Chemical class 0.000 claims description 447
- 238000009792 diffusion process Methods 0.000 claims description 111
- 239000011159 matrix material Substances 0.000 claims description 50
- 239000000149 chemical water pollutant Substances 0.000 claims description 48
- 230000007613 environmental effect Effects 0.000 claims description 43
- 239000010410 layer Substances 0.000 claims description 40
- 239000011148 porous material Substances 0.000 claims description 39
- 230000002285 radioactive effect Effects 0.000 claims description 37
- 239000011248 coating agent Substances 0.000 claims description 28
- 238000000576 coating method Methods 0.000 claims description 28
- 239000013077 target material Substances 0.000 claims description 25
- 230000015654 memory Effects 0.000 claims description 24
- 238000001179 sorption measurement Methods 0.000 claims description 24
- 238000002474 experimental method Methods 0.000 claims description 22
- 238000005253 cladding Methods 0.000 claims description 21
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 17
- 239000011707 mineral Substances 0.000 claims description 17
- 238000010276 construction Methods 0.000 claims description 13
- 238000011156 evaluation Methods 0.000 claims description 13
- 230000014509 gene expression Effects 0.000 claims description 12
- 239000011247 coating layer Substances 0.000 claims description 8
- 238000003933 environmental pollution control Methods 0.000 claims description 4
- 238000004891 communication Methods 0.000 claims description 3
- 239000007789 gas Substances 0.000 description 87
- 239000003344 environmental pollutant Substances 0.000 description 26
- 231100000719 pollutant Toxicity 0.000 description 26
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 21
- 230000000694 effects Effects 0.000 description 17
- 238000002386 leaching Methods 0.000 description 16
- 239000000700 radioactive tracer Substances 0.000 description 16
- 230000014759 maintenance of location Effects 0.000 description 14
- 238000011217 control strategy Methods 0.000 description 12
- 239000004927 clay Substances 0.000 description 11
- 241000282414 Homo sapiens Species 0.000 description 10
- 239000000758 substrate Substances 0.000 description 8
- 230000036541 health Effects 0.000 description 7
- 244000005700 microbiome Species 0.000 description 7
- 239000000126 substance Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- VEXZGXHMUGYJMC-UHFFFAOYSA-M Chloride anion Chemical compound [Cl-] VEXZGXHMUGYJMC-UHFFFAOYSA-M 0.000 description 5
- 239000002689 soil Substances 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 239000007788 liquid Substances 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 230000006978 adaptation Effects 0.000 description 3
- 238000003912 environmental pollution Methods 0.000 description 3
- 206010021143 Hypoxia Diseases 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 239000000356 contaminant Substances 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000004090 dissolution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000000945 filler Substances 0.000 description 2
- 230000007954 hypoxia Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000000813 microbial effect Effects 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 239000000941 radioactive substance Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 229910052770 Uranium Inorganic materials 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000031018 biological processes and functions Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Evolutionary Computation (AREA)
- Health & Medical Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to the technical field of numerical simulation and discloses a simulation method, a device, equipment and a medium of a migration process of a target substance, wherein the method comprises the steps of obtaining attribute data and environment data of an underground target substance corresponding to a target area, and constructing an underground target substance migration control equation based on the attribute data and the environment data of the underground target substance corresponding to the target area; performing parameter inversion on an underground target substance migration control equation to generate an underground target substance migration simulation model; and simulating the migration process of the underground target substance by using the underground target substance migration simulation model to generate target substance migration process simulation data. The invention realizes the accurate simulation of the migration process of the target substance in the underground medium.
Description
Technical Field
The invention relates to the technical field of numerical simulation, in particular to a simulation method, a device, equipment and a medium for a migration process of a target substance.
Background
The leaching solution of the landfill and the pollutants such as underground heavy metal salt are diffused to the surface through migration movement, and the pollutants are converged into rivers and lakes along with surface runoff, so that a plurality of environmental pollution problems are caused. Research on migration rules of underground pollutants in underground media can provide powerful technical support for scientific prevention, treatment and repair of environmental pollution problems. However, due to the strong non-uniformity and anisotropy of the subsurface medium, considering the natural decay process of contaminants, the associated convective diffusion theory has difficulty describing complex transport processes in which contaminants are present.
Therefore, how to accurately simulate the migration process of the target substance in the underground medium is a technical problem to be solved.
Disclosure of Invention
In view of the above, the present invention provides a method, apparatus, device and medium for simulating migration process of a target substance, so as to solve the technical problem of how to accurately simulate migration process of a target substance in an underground medium.
In a first aspect, the present invention provides a method for simulating a migration process of a target substance, including: acquiring attribute data and environment data of an underground target substance corresponding to a target area, and constructing an underground target substance migration control equation based on the attribute data and the environment data of the underground target substance corresponding to the target area; performing parameter inversion on an underground target substance migration control equation to generate an underground target substance migration simulation model; and simulating the migration process of the underground target substance by using the underground target substance migration simulation model to generate target substance migration process simulation data.
According to the simulation method for the migration process of the target substance, provided by the embodiment, the migration control equation of the underground target substance is constructed based on the attribute data of the underground target substance and the environmental data of the target area, inversion is performed on the migration control equation of the underground target substance by using experimental data, an underground target substance migration simulation model with high adaptation to the target environment is generated, and the migration process of the underground target substance of the target area is simulated by using the underground target substance migration simulation model, so that the migration process of the target substance in an underground medium is accurately simulated.
In an alternative embodiment, performing parametric inversion on a subsurface target material migration control equation to generate a subsurface target material migration simulation model, comprising: carrying out underground target substance migration experiments by utilizing environment data of a target area to generate a plurality of groups of underground target substance migration experiment data; and inverting parameters in the underground target substance migration control equation by utilizing a plurality of groups of underground target substance migration experimental data to generate an underground target substance migration simulation model.
According to the simulation method for the migration process of the target substance, the parameters in the migration control equation of the underground target substance are inverted by utilizing the environment data of the target area to perform the migration experiment of the underground target substance, so that the migration simulation model of the underground target substance is generated.
In an alternative embodiment, constructing a subsurface target substance migration control equation based on attribute data and environmental data of a subsurface target substance corresponding to a target region, comprises: when the target area is a first area, determining a vertical space position parameter of the first area, the concentration of heavy metal salt in the direction from the ore layer to the ground surface and a flow parameter of the heavy metal salt in the first area based on the environmental data; determining the attenuation period of the radioactive element in the heavy metal salt in the first area and the adsorption efficiency of the target plant on the heavy metal salt based on the attribute data of the underground target substance; determining a first attenuation coefficient based on the attenuation period of the radioactive element in the heavy metal salt in the first region and the adsorption efficiency of the target plant on the heavy metal salt; and acquiring an initial time derivative order, an initial spatial derivative order and an initial super-diffusion coefficient, and constructing an underground target substance migration control equation corresponding to the first region based on the first attenuation coefficient, the vertical spatial position parameter of the first region, the concentration of heavy metal salt in the direction from the ore layer to the earth surface, the flow parameter of the heavy metal salt in the first region, the initial time derivative order, the initial spatial derivative order and the initial super-diffusion coefficient.
According to the simulation method for the migration process of the target substance, provided by the embodiment, the migration control equation of the underground target substance in the target region of the heavy metal in the mineral layer is established, the influence of plants planted on the surface and having the characteristic of adsorbing radioactive substances is considered, so that the migration process of the heavy metal salt in the underground medium is simulated accurately, the diffusion and adsorption speed of the heavy metal salt in the mineral layer are effectively evaluated, and powerful technical support can be provided for the pollution control strategy of the target region.
In an alternative embodiment, an underground target substance migration control equation corresponding to the first region is constructed based on the first attenuation coefficient, the vertical spatial position of the first region, the concentration of heavy metal salt in the direction from the ore layer to the earth surface, the flow parameter of heavy metal salt in the first region, the initial time derivative order, the initial spatial derivative order and the super diffusion coefficient, and the underground target substance migration control equation corresponding to the first region is shown in the following relational expression:
wherein,first attenuation coefficient, < >>Vertical spatial position parameter representing the first area, < ->Representing a time parameter->Indicating the concentration of heavy metal salts in the direction from the seam to the surface,/->Representing the time hausdorff derivative, +.>Representing the spatial hausdorff derivative,/ >Initial time derivative order, ++>Initial spatial derivative order, ++>Flow parameters of heavy metal salts in the first region, < >>Representing the initial superdiffusion coefficient.
In an alternative embodiment, constructing a subsurface target substance migration control equation based on attribute data and environmental data of a subsurface target substance corresponding to a target region, comprises: when the target area is a second area, determining a vertical space position parameter of the second area, a seam position parameter of the second area, a pressure parameter of the second area, a concentration of heavy metal salt in a direction from a disposal warehouse to the ground surface, a concentration of heavy metal salt in a direction from the seam to the ground surface, a flow parameter of heavy metal salt gas migration in the second area and a flow parameter of heavy metal salt solute migration in the second area based on environmental data; determining a second attenuation coefficient based on attribute data of the subsurface target material; acquiring a time fractional order of initial heavy metal salt gas migration, a space fractional order of initial heavy metal salt gas migration and a diffusion coefficient of initial heavy metal salt gas migration, and constructing a control equation of heavy metal salt gas migration corresponding to a second area based on a second attenuation coefficient, a vertical spatial position parameter of the second area, the concentration of gas heavy metal salt in the direction from a disposal warehouse to the ground surface, a flow parameter of heavy metal salt gas migration in the second area, the time fractional order of initial heavy metal salt gas migration, the space fractional order of initial heavy metal salt gas migration and the flow parameter of heavy metal salt gas migration; acquiring a time fractional order of initial heavy metal salt solute migration, a space fractional order of initial heavy metal salt solute migration and a diffusion coefficient of initial heavy metal salt solute migration, and constructing a control equation of heavy metal salt solute migration corresponding to a second region based on a second attenuation coefficient, a vertical spatial position parameter of the second region, a seam position parameter of the second region, a pressure parameter of the second region, a concentration of heavy metal salt radioactive heavy metal salt of the heavy metal salt solute in a direction from a seam to the earth surface, a flow parameter of heavy metal salt solute migration in the second region, the time fractional order of initial heavy metal salt solute migration, the spatial fractional order of initial heavy metal salt solute migration and the flow parameter of heavy metal salt solute migration; and constructing an underground target substance migration control equation corresponding to the second region based on the heavy metal salt gas migration control equation and the heavy metal salt solute migration control equation.
According to the simulation method for the migration process of the target substance, the accurate simulation of the migration process of the heavy metal salt gas and the solute is realized by constructing the migration control equation of the heavy metal salt of the gas in the target area and the migration control equation of the radioactive heavy metal salt of the heavy metal salt solute, and the accurate and effective evaluation of the diffusion speed of the heavy metal salt in the target area is realized, so that powerful technical support can be provided for the pollution control strategy of the target area.
In an alternative embodiment, based on the control equation of heavy metal salt gas migration and the control equation of heavy metal salt solute migration, an underground target substance migration control equation corresponding to the second region is constructed, where the underground target substance migration control equation corresponding to the second region is represented by the following relational expression:
wherein,vertical spatial position parameter representing the second area, < >>Representing a time parameter->Concentration of gaseous heavy metal salt in the direction from the disposal reservoir to the surface,/->Concentration of heavy metal salt solute radioactive heavy metal salt in direction from ore layer to ground surface, < + >>Representing the pressure parameter of the second region, +.>，/>，/>Representing the time fractional derivative, whereinTime fractional order of initial heavy metal salt gas migration, +. >Time fractional order of initial heavy metal salt solute migration, +.>Spatial fractional order indicative of initial heavy metal salt gas migration, +.>Spatial fractional derivative representing Riemann-LiuVil type, ++>The derivative representing this type is the derivative of the Riemann-LiuVil class score order,/I>A seam location parameter indicative of a second zone, < + >>Spatial fractional order of initial heavy metal salt solute migration, +.>Flow parameter indicative of heavy metal salt gas migration in the second zone->A flow parameter indicative of heavy metal salt solute transport in the second zone,/->Diffusion coefficient indicating initial heavy metal salt gas migration, +.>Diffusion coefficient indicating initial heavy metal salt solute migration, +.>Representing a second attenuation coefficient.
In an alternative embodiment, constructing a subsurface target substance migration control equation based on attribute data and environmental data of a subsurface target substance corresponding to a target region, comprises: when the target area is a third area, determining a vertical spatial position parameter of the third area, the concentration of heavy metal salt in the direction from the mineral layer to the ground surface, a flow parameter of heavy metal salt under the upper coating layer and a flow parameter of heavy metal salt under the pore matrix based on the environmental data; determining a third attenuation coefficient based on attribute data of the subsurface target material; acquiring an initial time cut-off fractional order, an initial cut-off coefficient and an initial diffusion coefficient containing heavy metal salt under the upper coating, and constructing a heavy metal salt migration control equation under the upper coating corresponding to a third area based on a third attenuation coefficient, a vertical space position parameter of the third area, the concentration of the heavy metal salt in the direction from the mineral layer to the ground surface, a flow parameter containing the heavy metal salt under the upper coating, the initial time cut-off fractional order, the initial cut-off coefficient and the initial diffusion coefficient containing the heavy metal salt under the upper coating; acquiring an initial time fractional order and a diffusion coefficient of heavy metal salt under an initial pore matrix, and constructing a pore matrix heavy metal salt migration control equation corresponding to a third region based on a third attenuation coefficient, a vertical spatial position parameter of the third region, the concentration of heavy metal salt in the direction from a mineral layer to the earth surface, a flow parameter of heavy metal salt under the pore matrix, the initial time fractional order and the diffusion coefficient of heavy metal salt under the initial pore matrix; and constructing an underground target substance migration control equation corresponding to the third region based on the heavy metal salt migration control equation under the upper coating and the heavy metal salt migration control equation under the pore matrix.
According to the simulation method for the migration process of the target substance, through the migration control equation of the heavy metal salt under the upper coating and the migration control equation of the heavy metal salt under the pore matrix, the strong retention influence of the upper coating on the nuclide migration process is considered, so that the accurate simulation of the nuclide migration process of the heavy metal salt and the accurate and effective evaluation of the influence degree of the nuclide transport on the environment in a mining area are realized, and further powerful technical support is provided for the pollution control strategy of the target area.
In an alternative embodiment, based on the migration control equation of the heavy metal salt under the upper cladding layer and the migration control equation of the heavy metal salt under the pore matrix, an underground target substance migration control equation corresponding to the third region is constructed, and the underground target substance migration control equation corresponding to the third region is shown in the following relational expression:
wherein,vertical spatial position parameter representing third area, < ->Representing a time parameter->Concentration of heavy metal salts in the direction from the seam to the surface,/->Indicating seam location,/->) Indicating the thickness of the upper cladding>Representing the boundary of the earth's surface,time cut-off fractional derivative representing the carport type, wherein +.>Representing that the derivative of this type is the truncated fractional derivative of the carport class, 0 represents the initial moment,/- >Initial time cut-off fractional order, +.>Representing the initial truncated coefficient of the block,time division representing carport typeA derivative of the order>Representing that the derivative of this type is the fractional derivative of the carport class, 0 represents the initial moment, ++>Initial time fraction order, +.>Indicating the flow parameters of the heavy metal salt under the upper coating,indicating the flow parameters of the heavy metal salt under the pore matrix, < + >>Indicating the diffusion coefficient of the heavy metal salt under the initial coating,representing the diffusion coefficient of the heavy metal salt under the original pore matrix, < >>Representing a third attenuation coefficient.
In an alternative embodiment, constructing a subsurface target substance migration control equation based on attribute data and environmental data of a subsurface target substance corresponding to a target region, comprises: when the target area is a landfill, determining a spatial position parameter of the landfill, a concentration of landfill leachate and a flow parameter of the landfill leachate based on the environmental data; determining the attenuation coefficient of landfill leachate in an anaerobic environment based on attribute data of underground target substances; obtaining an initial fractional capacity coefficient, an initial cut-off coefficient, an initial fractional order and an initial diffusion coefficient, and constructing an underground target substance migration control equation in a target area based on an attenuation coefficient of landfill leachate in an anaerobic environment, a spatial position parameter of the landfill, the concentration of landfill leachate, a flow parameter of the landfill leachate, the initial fractional capacity coefficient, the initial cut-off coefficient, the initial fractional order and the initial diffusion coefficient.
According to the simulation method for the migration process of the target substance, the underground target substance migration control equation in the target area is constructed, so that accurate simulation of the leaching liquid migration process of the landfill is realized, the environmental influence degree caused by the transportation of the pollutant leaching liquid of the landfill is accurately and effectively estimated, and further powerful technical support is provided for the pollution control strategy of the target area.
In an alternative embodiment, the underground target substance migration control equation in the target area is constructed based on the attenuation coefficient of the landfill leachate in the anaerobic environment, the spatial location parameter of the landfill, the concentration of the landfill leachate, the flow parameter of the landfill leachate, the initial fractional capacity coefficient, the initial cutoff coefficient, the initial fractional order and the initial diffusion coefficient, and the underground target substance migration control equation in the target area is shown in the following relation:
wherein,spatial location parameters representing the target area, +.>Representing a time parameter->Indicating the concentration of landfill leachate, < + >>Representing the total duration of the simulation, +.>Initial fractional capacity coefficient, ++>Representing an initial sectionBreaking coefficient (F)>Initial fractional order, ++ >Flow parameters of landfill leachate, < ->Initial diffusion coefficient, ++>Attenuation coefficient of landfill leachate in anaerobic environment, < ->Representing the sign of the fractional derivative>Short for the truncated fractional derivative +.>Representing the truncated coefficient>Indicating the moment at which the fractional derivative starts memorization.
In an alternative embodiment, a migration process of a subsurface target substance is simulated using a subsurface target substance migration simulation model, generating simulation data, comprising: the method comprises the steps of obtaining initial concentration of underground target substances, solving an underground target substance migration simulation model by using a finite difference method based on the initial concentration of the underground target substances, and generating time-space distribution of the underground target substances in a target area.
According to the simulation method for the migration process of the target substance, provided by the embodiment, based on the initial concentration of the target substance monitored in real time, the underground target substance migration simulation model is solved by using a finite difference method, objective and reliable time-space distribution data of the underground target substance in the target area are generated, and powerful technical support can be provided for a pollution control strategy of the target area.
In an alternative embodiment, the method further comprises: and evaluating the migration process of the underground target substance by using the simulation data, and determining the environmental pollution control strategy of the target area based on the evaluation result.
The simulation method of the migration process of the target substance provided by the embodiment can provide powerful technical support for the pollution control strategy of the target area by utilizing objective and reliable time-space distribution data of the underground target substance in the target area.
In a second aspect, the present invention provides a simulation apparatus for a migration process of a target substance, comprising: the construction module is used for acquiring attribute data and environment data of the underground target substances corresponding to the target areas and constructing an underground target substance migration control equation based on the attribute data and the environment data of the underground target substances corresponding to the target areas; the inversion module is used for carrying out parameter inversion on the underground target substance migration control equation to generate an underground target substance migration simulation model; and the simulation module is used for simulating the migration process of the underground target substance by using the underground target substance migration simulation model to generate target substance migration process simulation data.
In a third aspect, the present invention provides a computer device comprising: the processor executes the computer instructions, thereby executing the simulation method of the migration process of the target substance according to the first aspect or any implementation manner corresponding to the first aspect.
In a fourth aspect, the present invention provides a computer-readable storage medium having stored thereon computer instructions for causing a computer to execute the simulation method of the target substance migration process of the first aspect or any one of the embodiments corresponding thereto.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the description of the embodiments or the prior art will be briefly described, and it is obvious that the drawings in the description below are some embodiments of the present invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for simulating a migration process of a target substance according to an embodiment of the present invention;
FIG. 2 is a flow chart of a simulation method of a migration process of a target substance according to an embodiment of the present invention;
FIG. 3 is a schematic representation of a heavy metal salt migration process under a fracture matrix according to an embodiment of the present invention;
FIG. 4 is a flow chart of a simulation method of a migration process of a target substance according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of a migration process of heavy metal salts in medium depth ores according to an embodiment of the present invention;
FIG. 6 is a flow chart of a simulation method of a migration process of a target substance according to an embodiment of the present invention;
FIG. 7 is a schematic diagram of a deep geological heavy metal solute migration process according to an embodiment of the present invention;
FIG. 8 is a flow chart (V) of a simulation method of a migration process of a target substance according to an embodiment of the present invention;
FIG. 9 is a schematic diagram of a landfill leachate migration process according to an embodiment of the present invention;
FIG. 10 is a flow chart of a simulation method of a migration process of a target substance according to an embodiment of the present invention;
FIG. 11 is a block diagram of a simulation apparatus of a target substance migration process according to an embodiment of the present invention;
fig. 12 is a schematic diagram of a hardware structure of a computer device according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The simulation method of the migration process of the target substance can be applied to electronic equipment for simulating the migration process or the diffusion process of the underground pollutant; the electronic device may include, but is not limited to, a notebook, desktop, mobile terminal, such as a cell phone, tablet, etc.; of course, the simulation method of the migration process of the target substance provided in the present specification may also be applied to an application program running in the above electronic device.
In recent years, along with the improvement of living standard, the requirements of human beings on the quality of living environment are gradually improved, and the quantity of living garbage is greatly increased, but because the technology and capability of the current garbage treatment cannot keep up with the increasing requirement of the treatment capacity, the leaching solution of the garbage landfill enters rivers and lakes and underground water along with the water body, meanwhile, heavy metal salt generated by the dissolution of underground minerals in water diffuses to the surface, the concentration of the heavy metal salt in the surface soil is finally increased, and the heavy metal salt can be converged into rivers and lakes along with the surface runoff, so that a plurality of environmental pollution problems are brought, and meanwhile, the human health is endangered; because of the complexity of underground media and different depth geological structures and solute migration characteristics, the near-earth side geological fracture matrixes are more, the water migration speed is high, the harm of migration and diffusion of heavy metal salt to the environment is higher, pollution to the water caused by heavy metal salt migration is considered, meanwhile, household garbage is mostly not treated by any technology, the household garbage is randomly stacked or simply buried, garbage leachate containing a large amount of pollutants is generated after the garbage is leached by rainwater, groundwater is continuously polluted, and the garbage leachate enters rivers and lakes through water circulation, and the domestic water is influenced; the heavy metal salt in the medium-depth ore can be dissolved, diffused and migrated, and is influenced by conditions such as hypoxia and microbial activity, gaseous substances can be generated, a crack channel is easy to form at a release position to the ground surface under the driving of gas, so that the heavy metal salt is quickly transported, and meanwhile, the heavy metal salt can be dissolved into a water body and quickly diffused along cracks, so that the pollution problem of the heavy metal salt in the medium-depth ore is also the key and difficult problem of river and lake safety; for deep geological minerals, because the radionuclides in the minerals have the characteristics of long half-life period, strong radioactivity and the like, the influence of the radionuclides can be as long as hundreds of years, various radionuclides diffuse and migrate to biospheres from a mineral layer more or less along with underground water flow or rock cracks through an upper coating, wherein the upper coating has stronger adsorptivity, and can effectively isolate the migration process of the radionuclides to the earth surface; because of strong non-uniformity and anisotropy of underground media, classical convection diffusion theory is difficult to describe the complex migration process of pollutants in the underground media, a multiple-media model needs to establish a plurality of mobile area and non-mobile area coupling models, the number of mobile areas and non-mobile areas needs to be determined according to actual conditions, the migration process of gaseous substances is extremely complex, meanwhile, due to the degradation effect of anaerobic microorganisms in the underground, the natural attenuation process exists in leaching solution pollutants, the model is complex in construction, the parameters are large, and the application is difficult.
In view of the above problems, an embodiment of the present invention provides a method for simulating a migration process of a target substance, by constructing an underground target substance migration control equation based on attribute data of an underground target substance and environmental data of a target area, inverting the underground target substance migration control equation by using experimental data, generating an underground target substance migration simulation model with high adaptation to a target environment, and simulating the migration process of the underground target substance of the target area by using the underground target substance migration simulation model, so as to accurately simulate the migration process of the target substance in an underground medium.
According to an embodiment of the present invention, there is provided an embodiment of a method of simulating a migration process of a target substance, it being noted that the steps shown in the flowcharts of the figures may be performed in a computer system such as a set of computer-executable instructions, and that although a logical order is shown in the flowcharts, in some cases the steps shown or described may be performed in an order different from that herein.
In this embodiment, a method for simulating a migration process of a target substance is provided, which may be used in the above-mentioned notebook, desktop computer, mobile terminal, such as a mobile phone, tablet computer, etc., fig. 1 is a flowchart of a method for simulating a migration process of a target substance according to an embodiment of the present invention, as shown in fig. 1, where the flowchart includes the following steps:
Step S101, obtaining attribute data and environment data of underground target substances corresponding to a target area, and constructing an underground target substance migration control equation based on the attribute data and the environment data of the underground target substances corresponding to the target area.
Specifically, the subsurface target material migrates or diffuses within the target area; underground target materials may include, but are not limited to, landfill leachate and heavy metal salts at different depths in the subsurface; attribute data for a subsurface target substance may include, but is not limited to, initial concentration of the target substance, flow parameters, etc.; the environmental data of the target area may include the depth at which the subsurface target material is located.
And S102, performing parameter inversion on an underground target substance migration control equation to generate an underground target substance migration simulation model.
Step S103, simulating the migration process of the underground target substance by using the underground target substance migration simulation model, and generating target substance migration process simulation data.
Specifically, the initial concentration of the underground target substance is obtained, and the underground target substance migration simulation model is solved by using a finite difference method based on the initial concentration of the underground target substance, so that the time-space distribution of the underground target substance in the target area is generated.
Further, the underground target substance migration process is evaluated by using the simulation data, and the environmental pollution control strategy of the target area is determined based on the evaluation result.
According to the simulation method for the migration process of the target substance, provided by the embodiment, the migration control equation of the underground target substance is constructed based on the attribute data of the underground target substance and the environmental data of the target area, inversion is performed on the migration control equation of the underground target substance by using experimental data, an underground target substance migration simulation model with high adaptation to the target environment is generated, and the migration process of the underground target substance of the target area is simulated by using the underground target substance migration simulation model, so that the migration process of the target substance in an underground medium is accurately simulated.
In this embodiment, a method for simulating a migration process of a target substance is provided, which may be used in the above-mentioned notebook, desktop computer, mobile terminal, such as a mobile phone, tablet computer, etc., fig. 2 is a flowchart of a method for simulating a migration process of a target substance according to an embodiment of the present invention, as shown in fig. 2, where the flowchart includes the following steps:
step S201, obtaining attribute data and environment data of an underground target substance corresponding to a target area, and constructing an underground target substance migration control equation based on the attribute data and the environment data of the underground target substance corresponding to the target area; please refer to step S101 in the embodiment shown in fig. 1 in detail, which is not described herein.
And step S202, performing parameter inversion on the underground target substance migration control equation to generate an underground target substance migration simulation model.
Specifically, the step S202 includes:
in step S2021, the underground target substance transfer experiment is performed using the environmental data of the target area, so as to generate multiple sets of underground target substance transfer experiment data.
Specifically, an experimental model is constructed by utilizing environmental data of a target area, driving conditions of the experimental model, such as initial concentration of underground target substances, depth of the underground target substances in the target area and the like, are set, and based on the driving conditions of the experimental model, the experimental model is driven to perform underground target substance migration experiments, so that multiple groups of underground target substance migration experimental data are generated.
And S2022, inverting parameters in the underground target substance migration control equation by utilizing the plurality of groups of underground target substance migration experimental data to generate an underground target substance migration simulation model.
Step S203, simulating the migration process of the underground target substance by using the underground target substance migration simulation model to generate target substance migration process simulation data; please refer to step S103 in the embodiment shown in fig. 1 in detail, which is not described herein.
According to the simulation method for the migration process of the target substance, the parameters in the migration control equation of the underground target substance are inverted by utilizing the environment data of the target area to perform the migration experiment of the underground target substance, so that the migration simulation model of the underground target substance is generated.
In this embodiment, a method for simulating a migration process of a target substance is provided, which may be used in the above-mentioned notebook, desktop computer, mobile terminal, such as mobile phone, tablet computer, etc., where the attribute data of the underground target substance includes attribute data of heavy metal salt, and the method is applied to simulate a migration process of heavy metal salt in a fracture matrix as shown in fig. 3, and fig. 4 is a flowchart of a simulation method of a migration process of a target substance according to an embodiment of the present invention, as shown in fig. 4, where the flowchart includes the following steps:
step S401, obtaining attribute data and environment data of the underground target substance corresponding to the target area, and constructing an underground target substance migration control equation based on the attribute data and the environment data of the underground target substance corresponding to the target area.
Specifically, it should be noted that, when the underground target substance is a heavy metal salt, that is, when the simulation method of the migration process of the target substance is used to simulate the migration process of the underground heavy metal salt, the target area may be divided into a first area, a second area and a third area according to the depth of the ore deposit where the underground target substance is located; the simulated depth ordering is: the third region is larger than the second region, and the second region is larger than the first region; the first zone is a fracture matrix, such as within 50 meters of the depth of the seam, the second zone is a medium depth seam, such as within 200 meters of the depth of the seam, and the third zone is a deep geological layer, such as greater than 200 meters; when the target area is the first area, the step S401 includes:
step S4011, determining a vertical spatial position parameter of the first region, a concentration of heavy metal salt in a direction from the seam to the earth surface, and a flow parameter of heavy metal salt in the first region based on the environmental data.
Specifically, the environmental data includes a vertical spatial location parameter of the first region, a concentration of heavy metal salt in a seam-to-surface direction, and a flow parameter of heavy metal salt in the first region.
Step S4012, determining the decay period of the radioactive element in the heavy metal salt in the first region and the adsorption efficiency of the target plant to the heavy metal salt based on the attribute data of the underground target substance.
Specifically, the adsorption efficiency of the target plant to the heavy metal saltHeavy metal salt adsorption experiments can be carried out according to a denser planting mode, and the adsorption efficiency of target plants on heavy metal salt is measured>Is a value of (2).
Step S4013, determining a first attenuation coefficient based on the attenuation period of the radioactive element in the heavy metal salt in the first region and the adsorption efficiency of the target plant to the heavy metal salt.
Specifically, the data of the ore deposit to be evaluated, such as the depth h of the ore deposit, the type of heavy metal salt, whether radioactive elements exist in the heavy metal salt or not, and the decay period of the radioactive elements in the heavy metal salt, are obtainedSimultaneously obtaining the adsorption efficiency of the target plant on the heavy metal salt>The method comprises the steps of carrying out a first treatment on the surface of the And simultaneously, parameters such as geological structure, PH value, temperature, conductivity and the like of the fracture matrix are investigated, and the anisotropy and the non-uniformity of the matrix are evaluated and are used as the basis for parameter values in an underground target substance migration control equation in the first area.
Further, a first attenuation coefficientCan be determined by the following relation (1):
（1）
wherein,represents the first attenuation coefficient (i.e. the attenuation coefficient of the heavy metal salt),>represents the adsorption efficiency of the target plant to the heavy metal salt, < ->Represents the decay period of the radioactive element in the heavy metal salt, if there is no radioactivity,/i >。
Step S4014, obtaining an initial time derivative order, an initial spatial derivative order and an initial super-diffusion coefficient, and constructing an underground target substance migration control equation corresponding to the first region based on the first attenuation coefficient, the vertical spatial position parameter of the first region, the concentration of heavy metal salt in the direction from the mineral seam to the earth surface, the flow parameter of heavy metal salt in the first region, the initial time derivative order, the initial spatial derivative order and the initial super-diffusion coefficient.
Specifically, the migration control equation of the underground target substance corresponding to the first region is represented by the following relational expression (2):
（2）
wherein,representing a first attenuation coefficient (i.e., the attenuation coefficient of the heavy metal salt), characterizing the attenuation rate of the heavy metal salt,representation ofVertical spatial position parameter of the first zone (mineral seam to ground direction positive),>representing a time parameter->Indicating the concentration of heavy metal salts in the direction from the seam to the surface,/->Representing the derivative of the time Hausdroff ++>Representing the spatial hausdorff derivative,/>Representing an initial time derivative order, wherein the initial time derivative order adopts a time Hastethodor derivative order to represent the retention effect of local non-connected fracture points on heavy metal salt, and the method is characterized in that ∈10 is used for the retention effect of the local non-connected fracture points on heavy metal salt>Representing the initial spatial derivative order, wherein the initial spatial derivative order adopts spatial Hastethodor derivative order to represent the super diffusion process of heavy metal salt along the direction of the communicated fracture, and the +. >The lower the value is, the>The stronger the retention of the dead spots, the +.>The lower the value, the stronger the superdiffusion effect, here +.>，/>Representing the flow parameters of the heavy metal salt in the first zone, which is valued by the underground fracture matrixQualitative water conductivity determination,/->Representing the initial superdiffusion coefficient.
And step S402, performing parameter inversion on the underground target substance migration control equation to generate an underground target substance migration simulation model.
Specifically, a field crack matrix chloride ion tracer diffusion experiment (the experimental model structure is similar to that shown in fig. 3) is carried out pertinently, the crack matrix is similar to that investigated in step S4011, the position of the tracer chloride ion is located at a distance h from the earth surface, the method is consistent with step S4011, the concentration of the chloride ion tracer in the earth surface is monitored from the release position of the chloride ion tracer until the chloride ion tracer reaches the earth surface concentration which is a preset release concentration, such as 2%, and the process is finished; inversion of experimental data to obtain the order of the time hausdorff derivative in the migration control equation of the subsurface target material in the first zoneSpatial hausdorff derivative->And generating an underground target substance migration simulation model by using the super-diffusion coefficient D, wherein the underground target substance migration simulation model in the embodiment is an underground target substance migration simulation model in the first area and is a time-space Hasteorf fractal convection diffusion adsorption attenuation model.
Step S403, simulating the migration process of the underground target substance by using the underground target substance migration simulation model, and generating target substance migration process simulation data.
Specifically, let theI.e. no plants are on the earth surface, solving a time-space Haoskov fractal convection diffusion adsorption attenuation model in the step S402 by using a finite difference method to obtain the time-space distribution of heavy metal salt in a period of time after the heavy metal salt in the ore layer; wherein the time is generally determined by the time of arrival of the heavy metal salt at the surface, e.g. in the second half of the arrival timeAnd (5) calculating an over-simulation value.
If the heavy metal salt has radioactivity, the adsorption efficiency is measured through a heavy metal salt adsorption experimentAnd (3) solving a time-space Haosduff fractal convection diffusion adsorption attenuation model in the step S402 by using a finite difference method to obtain the time-space distribution of the heavy metal salt in a period of time after the treatment library leaks.
Further, evaluating the migration process of the underground target substances by using the simulation data, and determining a pollution control strategy; if the simulation data indicate that the heavy metal salt does not reach the ground surface within a period of time, the heavy metal salt in the ore deposit is considered to not influence the water quality of the river and the lake; if the simulation data show that the heavy metal salt reaches the concentration which does not exceed the safety standard (meets the ecological environment safety and human health standard) within a period of time, the simulation time is increased by more than the original simulation time, for example, the simulation time is increased by 1 year, and finally, if the heavy metal salt on the surface of the earth is gradually reduced, the heavy metal salt is suggested to be treated by adopting a method of densifying plants; if the simulation data characterize that the heavy metal salt exceeds the safety standard in a period of time, selecting plants with better adsorptivity, continuing simulation, and if the simulation data still do not meet the safety standard, considering unsafe.
Further, if the ore layer contains various heavy metal salts, the steps S401 to S403 are repeated, the pollution is considered for a plurality of times, and the environmental assessment report is submitted to note the conclusion and suggestion.
According to the simulation method for the migration process of the target substance, the migration speed of the water body of the crack substrate is high, the anisotropy and nonlinearity of the crack substrate at the near side are relatively strong, the migration and transformation process of the heavy metal salt is generally expressed as a super-diffusion process along the direction of the communicated cracks and a retention effect along the direction of the dead cavern of the non-communicated cracks, so that the two abnormal diffusion phenomena need to be considered simultaneously in the establishment of a model, the adsorption process of heavy metal salt plants needs to be considered in the established model, if the heavy metal salt has radioactivity, the natural attenuation of the heavy metal salt needs to be considered, the time space Haos-Duofu fractal convection diffusion adsorption attenuation model describing the migration and transformation of the underground heavy metal salt is established in the embodiment, the plant influence with the characteristic of adsorbing radioactive substances is considered in the surface planting, the migration process of the heavy metal salt in the underground medium is accurately simulated, the diffusion and the adsorption speed of the heavy metal salt in the mineral layer can be effectively estimated, and powerful technical support can be provided for pollution prevention and control strategies of target areas.
In this embodiment, a method for simulating a migration process of a target substance is provided, which may be used in the above-mentioned notebook, desktop computer, mobile terminal, such as mobile phone, tablet computer, etc., where the attribute data of the underground target substance includes attribute data of heavy metal salt, and the method is applied to simulate a migration process of heavy metal salt in a medium-depth mine as shown in fig. 5, and fig. 6 is a flowchart of a method for simulating a migration process of a target substance according to an embodiment of the present invention, as shown in fig. 6, where the flowchart includes the following steps:
step S601, obtaining attribute data and environment data of an underground target substance corresponding to a target area, and constructing an underground target substance migration control equation based on the attribute data and the environment data of the underground target substance corresponding to the target area.
Specifically, data of various heavy metal salts to be evaluated, such as water solubility, gas generation proportion, whether radioactive elements exist or not, attenuation period S of the radioactive elements and depth h (within 200 meters) of a mineral deposit are obtained, parameters of geological structure, pH value, temperature, conductivity and the like of an underground substrate are investigated, and substrate anisotropy and non-uniformity are evaluated and used as the basis for taking values of parameters of an underground target substance migration control equation in a second area.
Specifically, when the target area is the second area, the step S601 includes:
step S6011, determining, based on the environmental data, a vertical spatial position parameter of the second area, a seam position parameter of the second area, a pressure parameter of the second area, a concentration of gaseous heavy metal salt in a direction from the disposal reservoir to the earth surface, a concentration of heavy metal salt solute radioactive heavy metal salt in a direction from the seam to the earth surface, a flow parameter of heavy metal salt gas migration in the second area, and a flow parameter of heavy metal salt solute migration in the second area.
Specifically, the environmental data includes a vertical spatial location parameter of the second region, a seam location parameter of the second region, a pressure parameter of the second region, a concentration of gaseous heavy metal salt in a direction from the disposal reservoir to the surface, a concentration of heavy metal salt solute radioactive heavy metal salt in a direction from the seam to the surface, a flow parameter of heavy metal salt gas migration in the second region, and a flow parameter of heavy metal salt solute migration in the second region.
Step S6012 of determining a second attenuation coefficient based on the attribute data of the underground target substance;
specifically, the attribute data of the underground object substance includes the attenuation period of the radioactive element in the heavy metal salt, the second attenuation coefficient Can be determined by the following relation (3):
（3）
wherein,represents the second attenuation coefficient (i.e. the attenuation coefficient of the heavy metal salt),>represents the decay period of the radioactive element in the heavy metal salt.
And step S6013, acquiring the time fractional order of the initial heavy metal salt gas migration, the space fractional order of the initial heavy metal salt gas migration and the diffusion coefficient of the initial heavy metal salt gas migration.
And step S6014, constructing a control equation of heavy metal salt gas migration corresponding to the second region based on the second attenuation coefficient, the vertical spatial position parameter of the second region, the concentration of the heavy metal salt gas in the direction from the treatment warehouse to the ground surface, the flow parameter of heavy metal salt gas migration in the second region, the time fractional order of initial heavy metal salt gas migration, the spatial fractional order of initial heavy metal salt gas migration and the flow parameter of heavy metal salt gas migration.
Specifically, the control equation of the migration of the heavy metal salt gas corresponding to the second region is shown in the following relational expression (4):
（4）
wherein,vertical spatial position parameter representing the second area, < >>Representing a time parameter->Indicating the concentration of gaseous heavy metal salts in the direction from the treatment reservoir to the surface,/- >Representing the pressure parameter of the second region, +.>，/>Representing the time fractional derivative, wherein +.>Time fractional order representing initial heavy metal salt gas migration, +.>Spatial fractional order indicative of initial heavy metal salt gas migration, +.>Flow parameter indicative of heavy metal salt gas migration in the second zone->Diffusion coefficient indicating initial heavy metal salt gas migration, +.>The second attenuation coefficient (the attenuation coefficient of the heavy metal salt) is represented.
Step S6015, obtaining a time fractional order of initial heavy metal salt solute migration, a space fractional order of initial heavy metal salt solute migration and a diffusion coefficient of initial heavy metal salt solute migration.
Step S6016, constructing a control equation of heavy metal salt solute migration corresponding to the second region based on the second attenuation coefficient, the vertical space position parameter of the second region, the seam position parameter of the second region, the pressure parameter of the second region, the concentration of heavy metal salt solute radioactive heavy metal salt in the direction from the seam to the ground surface, the flow parameter of heavy metal salt solute migration in the second region, the time fractional order of initial heavy metal salt solute migration, the space fractional order of initial heavy metal salt solute migration and the flow parameter of heavy metal salt solute migration.
Specifically, the control equation of the heavy metal salt solute migration corresponding to the second region is shown in the following relational expression (5):
（5）
wherein,vertical spatial position parameter representing the second area, < >>Representing a time parameter->Representing the concentration of heavy metal salt solute radioactive heavy metal salt in the direction from the ore layer to the ground surface, ++>Representing the time fractional derivative, wherein +.>Time fractional order representing initial heavy metal salt solute migration, +.>Spatial fractional derivative representing Riemann-LiuVil type, ++>The derivative representing this type is the derivative of the Riemann-LiuVil class score order,/I>A seam location parameter indicative of a second zone, < + >>Representing the initial spatial fraction order, +.>A flow parameter indicative of heavy metal salt solute transport in the second zone,/->Diffusion coefficient indicating initial heavy metal salt solute migration, +.>Representing a second attenuation coefficient (i.e., the attenuation coefficient of the heavy metal salt).
And step S6017, constructing an underground target substance migration control equation corresponding to the second area based on the heavy metal salt gas migration control equation and the heavy metal salt solute migration control equation.
Specifically, the migration control equation of the underground target substance corresponding to the second region is shown in the following relational expressions (4) to (5):
（4）
（5）
Wherein,vertical spatial position parameter representing the second area, < >>Representing a time parameter->Concentration of gaseous heavy metal salt in the direction from the disposal reservoir to the surface,/->Concentration of heavy metal salt solute radioactive heavy metal salt in direction from ore layer to ground surface, < + >>Representing the pressure parameter of the second region, +.>，/>，/>Representing the time fractional derivative, whereinTime fractional order of initial heavy metal salt gas migration, +.>Time fractional order of initial heavy metal salt solute migration, +.>Spatial fractional order indicative of initial heavy metal salt gas migration, +.>Spatial fractional derivative representing Riemann-LiuVil type, ++>The derivative representing this type is the derivative of the Riemann-LiuVil class score order,/I>A seam location parameter indicative of a second zone, < + >>Spatial fractional order of initial heavy metal salt solute migration, +.>Flow parameter indicative of heavy metal salt gas migration in the second zone->A flow parameter indicative of heavy metal salt solute transport in the second zone,/->Diffusion coefficient indicating initial heavy metal salt gas migration, +.>Diffusion coefficient indicating initial heavy metal salt solute migration, +.>Representing a second attenuation coefficient.
And step S602, performing parameter inversion on the underground target substance migration control equation to generate an underground target substance migration simulation model.
Specifically, a soil column experiment (an experimental model structure is similar to that shown in fig. 5) for forming and diffusing indoor gas is performed pertinently, the filling material is similar to the underground medium investigated in the step 601, the filling material is compact, the gas tracer is buried as deeply as possible, and meanwhile, the fluorescent agent is buried as a solute tracer; and monitoring the concentration of the tracer in the surface from the release position of the gas tracer to the surface, recording the time for the tracer to reach the surface, and finishing the gas migration experiment and the fluorescent solute experiment data after the fluorescent solute tracer reaches the surface concentration which is 2% of the preset release concentration.
Order the=0, inversion of the gas migration experimental data to obtain the gas heavy metal salt migration control equation in the second region, i.e. the time fractional order in relation (4)>Spatial fractional order of gas migration +.>Gas diffusion coefficient->The method comprises the steps of carrying out a first treatment on the surface of the And then the obtained time fraction order +.>Spatial fractional order of gas migration +.>Gas diffusion coefficient->Substituting the relation (4) to generate a gas heavy metal salt migration simulation model in the second area, wherein the simulation model is a gas fractional order Darcy-convection-diffusion-attenuation model.
Order the =0, and inverting the experimental data of the fluorescent agent solute transport to obtain a heavy metal salt transport control equation of the heavy metal salt solute radioactivity heavy metal salt in the second region, namely the time fractional order +.>Spatial fraction order +.>Diffusion coefficient->The method comprises the steps of carrying out a first treatment on the surface of the And then the obtained time fraction order +.>Spatial fraction order +.>Diffusion coefficient->Substituting the heavy metal salt into the relational expression (5) to generate a heavy metal salt solute radioactivity heavy metal salt migration simulation model in the second area, wherein the simulation model is a space fractional order convection diffusion attenuation model for rapid diffusion of heavy metal salt.
Step S603, simulating the migration process of the underground target substance by using the underground target substance migration simulation model, and generating target substance migration process simulation data.
Specifically, let the=0, i.e. no radioactivity of the heavy metal salt, solving the gas fractional order darcy-convection-diffusion-attenuation model in step S602 by using a finite difference method to obtain gas heavy metal salt time-space distribution within a period of time after heavy metal salt transport; the time is generally determined by the time of arrival of the gas at the surface, for example, in the next year of arrival, the time of arrival is calculated by an analog value.
Order the=0, i.e. no radioactivity of heavy metal salt, solving a spatial fractional order convective diffusion attenuation model of rapid diffusion of one heavy metal salt in step S602 by using a finite difference method, so as to obtain the time-space distribution of medium-radioactivity heavy metal salt solute in a period of time after heavy metal salt in medium-depth ore is transported; the time is generally determined by the time of arrival of the heavy metal salt of the radioactive gas at the surface, and is calculated by analog values in the following year, for example, in the following year.
Order the=/>Namely, considering the attenuation period of the radioactive heavy metal salt, solving a gas fractional order Darcy-convection-diffusion-attenuation model in the step S602 by using a finite difference method to obtain the time-space distribution of the heavy metal salt in the medium-depth mine in a period of time after the treatment library leaks; the time is generally determined by the time of arrival of the heavy metal salt at the surface, for example, the time of arrival is calculated by analog values in two years after arrival.
Order the=1/S, namely, considering the attenuation period of the radioactive heavy metal salt, solving a spatial fractional order convection diffusion attenuation model of rapid diffusion of the heavy metal salt in the step S602 by using a finite difference method, and obtaining the time-space distribution of heavy metal salt solutes in the medium-depth mine within a period of time after the treatment library leaks; the time is generally determined by the time of arrival of the heavy metal salt at the surface, for example, the time of arrival is calculated by analog values in two years after arrival.
Further, evaluating the migration process of the underground target substances by using the simulation data, and determining a pollution control strategy; if the simulation data indicates that the release amount of the gas heavy metal salt after reaching the gas heavy metal salt is harmful to human health within a period of time (for example, equal to half a month after reaching the gas heavy metal salt), the water body near the ore deposit is considered to need to be monitored and treated in real time.
If the simulation data indicates that the release amount of the heavy metal salt solute after reaching the heavy metal salt solute does not exceed the safety standard (meeting the ecological environment safety and human health standards) within a period of time (for example, the time is equal to the year after the heavy metal salt is in the gas radioactivity), otherwise, the simulation data are considered to be unsafe, and the treatment methods such as a chemical fixing method and the like are recommended to reduce the content of the heavy metal salt solute.
If the simulation data indicates that the release of the gaseous heavy metal salt after the arrival of the gaseous heavy metal salt is harmful to human health within a period of time (e.g., equal to one month after the arrival of the gaseous heavy metal salt), then the vicinity of the seam is considered not to be susceptible to human activity.
If the simulation data indicates that the release amount of the radioactive heavy metal salt solute after reaching the gas radioactive heavy metal salt solute does not exceed the safety standard (meeting the ecological environment safety and human health standards) within a period of time (for example, two years after reaching the gas radioactive heavy metal salt), otherwise, the simulation data are considered to be unsafe, and a treatment method such as a chemical fixing method is suggested to reduce the content of the heavy metal salt solute.
Further, if the medium-depth ore contains multiple heavy metal salts, the steps S601 to S603 are repeated, the pollution of the multiple heavy metal salts is considered for multiple times, the environmental evaluation report is submitted, the conclusion and the suggestion are noted, and if the water solubility of the medium-activity heavy metal salts is poor, the consideration is not needed.
According to the simulation method for the migration process of the target substance, due to the fact that the dissolution, diffusion and migration processes of the heavy metal salt are influenced by the conditions of hypoxia, microbial activity and the like, gaseous substances can be generated, a crack channel is easy to form at a release position to the ground surface under the driving of the gas, so that the rapid migration of the heavy metal salt is caused, meanwhile, the heavy metal salt can be dissolved into a water body and rapidly diffuses along cracks, meanwhile, an underground matrix is relatively strong in anisotropy and nonlinearity, the migration and conversion process of the heavy metal salt shows two stages of gas and solute, wherein the migration and conversion process of the heavy metal salt is relatively high in speed, compared with the migration rate of the gas, the solute is relatively slow, but the main heavy metal salt migrates through the solute transportation stage, the first stage is before the migration of the gas to the ground surface, the communication channel directly reaching the ground surface is not formed at the moment, the solute transportation is generally slow to be influenced by the matrix, and after the migration of the gas to the ground surface, the communication channel is formed at the moment, the solute transportation is hindered by the matrix and the rapid migration and simultaneously influenced, and therefore, a model of the diffusion of the heavy metal salt in the space of the three-dimensional diffusion-order of the heavy metal is required to be established, and the diffusion-phase of the diffusion of the heavy metal in the gas level in the radiation state is relatively high in the three-dimensional diffusion-phase transition modes of the diffusion of the heavy metal; by constructing a target area gas heavy metal salt migration control equation and a heavy metal salt solute radioactivity heavy metal salt migration control equation, accurate simulation of heavy metal salt gas and solute migration processes and accurate and effective evaluation of heavy metal salt diffusion speed in medium-depth ores are realized, and powerful technical support can be provided for pollution control strategies of the target area, such as river and lake water protection.
In this embodiment, a method for simulating a migration process of a target substance is provided, which may be used in the above-mentioned notebook, desktop computer, mobile terminal, such as mobile phone, tablet computer, etc., where the attribute data of the underground target substance includes attribute data of heavy metal salt, and the method is applied to simulate a migration process of deep geological heavy metal solute as shown in fig. 7, and fig. 8 is a flowchart of a method for simulating a migration process of a target substance according to an embodiment of the present invention, as shown in fig. 8, where the flowchart includes the following steps:
step S801, attribute data and environment data of an underground target substance corresponding to a target area are acquired, and an underground target substance migration control equation is constructed based on the attribute data and environment data of the underground target substance corresponding to the target area.
Specifically, acquiring heavy metal salt data of a plurality of nuclides to be evaluated, such as attenuation periods of radioactive elements in heavy metal salt, acquiring nuclide migration process data of natural nuclear ores with similar geological structures, and parameters such as mineral layer depth, thickness and permeability of an upper coating, geological structure, pH value, temperature, conductivity and the like of an underground substrate, evaluating substrate anisotropy and non-uniformity, and taking the parameters as the basis of the values of underground target substance migration control equation parameters; wherein, as shown in FIG. 7, the depth of the seam is equal to ，/>Indicating the depth of burial of ore body->Indicating the thickness of the overlying layer.
Specifically, when the target area is the third area, the step S801 includes:
and step S8011, determining the vertical space position parameter of the third area, the concentration of the heavy metal salt in the direction from the ore layer to the ground surface, the flow parameter of the heavy metal salt under the upper coating layer and the flow parameter of the heavy metal salt under the pore matrix based on the environmental data.
Step S8012, determining a third attenuation coefficient based on attribute data of the subsurface target material.
And step S8013, obtaining an initial time cut-off fractional order, an initial cut-off coefficient and a diffusion coefficient of heavy metal salt under the initial upper coating.
And step S8014, constructing an upper-cladding-containing lower heavy metal salt migration control equation corresponding to the third region based on the third attenuation coefficient, the vertical space position parameter of the third region, the concentration of the heavy metal salt in the direction from the ore layer to the ground surface, the flow parameter of the upper-cladding-containing lower heavy metal salt, the initial time cutoff fractional order, the initial cutoff coefficient and the initial diffusion coefficient of the upper-cladding-containing lower heavy metal salt.
Specifically, the migration control equation of the heavy metal salt under the upper cladding layer in the third region is shown in the following relation (6):
（6）
wherein, Vertical spatial position parameter representing third area, < ->Representing a time parameter->Indicating the concentration of heavy metal salts in the direction from the seam to the surface,/->Indicating seam location,/->) Indicating the thickness of the upper cladding>Representing the surface boundary, typically +.>，/>Time-truncated fractional derivative representing the Caplato (Caputo) type, wherein +.>Representing that the derivative of the type is the truncated fractional derivative of the carport class, 0 represents the initial moment,representing the order of the initial time cut-off fraction, < >>Characterization of the retention of the heavy metal salt by the upper coating, < ->The lower the value, the stronger the retention effect, < ->Representing the initial truncated coefficient, ++>The higher the value, the more the heavy metal salt in the upper coating is near the saturation stage, because the upper coating is thicker and has strong adsorptivity, the more the upper coating is>，/>Indicating the flow parameters of the heavy metal salt under the upper coating, < >>Representing the initial presence of heavy metal salts under the upper coatingDiffusion coefficient, < >>Representing a third attenuation coefficient (i.e., the attenuation coefficient of the heavy metal salt).
And step S8015, obtaining an initial time fractional order and a diffusion coefficient of the heavy metal salt under the initial pore matrix.
And step S8016, constructing a heavy metal salt migration control equation under the pore matrix corresponding to the third region based on the third attenuation coefficient, the vertical space position parameter of the third region, the concentration of the heavy metal salt in the direction from the ore layer to the earth surface, the flow parameter of the heavy metal salt under the pore matrix, the initial time fractional order and the diffusion coefficient of the heavy metal salt under the initial pore matrix.
Specifically, the migration control equation of the heavy metal salt under the pore matrix in the third region is shown in the following relational expression (7):
（7）
wherein,vertical spatial position parameter representing third area, < ->Representing a time parameter->Indicating the concentration of heavy metal salts in the direction from the seam to the surface,/->Representing the boundary of the earth's surface->Time fractional derivative representing the type of carport, wherein +.>Representing that the derivative of this type is the fractional derivative of the carport class, 0 represents the initial moment, ++>Representing the time fractional order, representing the retention degree of the underground medium on the heavy metal salt, +.>，/>The lower the value is, the stronger the retention effect is, and the retention effect of the upper coating layer on the heavy metal salt is larger than that of the common underground medium, so +.>，/>Representing the flow parameters of the heavy metal salt under the pore matrix, determined by the water conductivity of the medium, +.>Representing the diffusion coefficient of the heavy metal salt under the initial pore matrix, the same +.>，/>Representing a third attenuation coefficient (attenuation coefficient of heavy metal salt) characterizing the rate of attenuation of the nuclide in the heavy metal salt,/->，/>Represents the decay period of the radioactive element in the heavy metal salt.
And step S8017, constructing an underground target substance migration control equation corresponding to the third area based on the heavy metal salt migration control equation under the upper cladding and the heavy metal salt migration control equation under the pore matrix.
Specifically, the migration control equations of the underground target substance corresponding to the third region are shown in the following relational expressions (6) to (7).
And step S802, performing parameter inversion on the underground target substance migration control equation to generate an underground target substance migration simulation model.
Specifically, a soil column experiment (an experimental model structure is similar to that shown in fig. 7) of the diffusion of the indoor clay is performed pertinently, the filler is clay, the filler is similar to the structure function of the upper cladding layer, the thickness can be less than half of the actual thickness, meanwhile, the solute tracer is buried, the concentration of the tracer in the clay is monitored, the time of the tracer penetrating through the clay is recorded, and the obtained migration experiment tracer solute experimental data is collated.
Inverting parameters in equations shown in a relation (6) and a relation (7) respectively by using the acquired migration experiment tracer solute experimental data and the acquired nuclide migration process data of the natural uranium ores with similar geological structures, and time-truncated fractional orderCut-off coefficient->Diffusion coefficient->Time fraction order +.>Diffusion coefficient->And generating an underground target substance migration simulation model in a third area, wherein the underground target substance migration simulation model in the third area comprises a time-cut fractional order convection diffusion attenuation model containing heavy metal salt migration under an upper coating layer and a time-cut fractional order convection diffusion attenuation model containing heavy metal salt migration under a pore matrix.
Step S803, simulating the migration process of the underground target substance by using the underground target substance migration simulation model, and generating target substance migration process simulation data.
Specifically, a finite difference method is used to solve the time-cut fractional order convective diffusion attenuation model containing heavy metal salt migration under the upper cladding layer and the time-fractional order convective diffusion attenuation model containing heavy metal salt migration under the pore matrix in the step S802, so as to obtain the time-space distribution of the heavy metal salt containing nuclein.
If the simulation data represent that the radionuclide does not reach the ground surface after twenty years, and the radionuclide does not exceed the safety standard (meets the ecological environment safety and human health standard), the environmental risk is considered to be low, otherwise, the radionuclide is considered to be unsafe, and the treatment methods such as a chemical fixing method and the like are suggested to reduce the content of the heavy metal salt containing the radionuclide.
Further, if the ore layer contains a plurality of highly-radioactive nuclides, the steps S801 to S803 are repeated, the safety of the plurality of nuclides is considered a plurality of times, and environmental impact assessment reports are submitted to note conclusions and suggestions.
According to the simulation method of the migration process of the target substance, which is provided by the embodiment, due to the characteristics of long half life, strong radioactivity and the like of the mine nuclide, the influence of the mine nuclide on the surrounding environment can be hundreds of years, so that various radionuclides migrate into the biosphere in the form of heavy metal salts through the upper coating and the underground medium; the heterogeneity of the underground medium causes the difference of the flow velocity and the path of underground water, and meanwhile, the upper coating has strong retention effect on the nuclide migration process, and the strong retention effect can be attenuated with time, so that the migration process of the underground medium under the two substrate conditions of the upper coating and the underground medium is simulated and needs to be described, and the time-cut fractional order convection diffusion attenuation equation containing the migration of the heavy metal salt under the upper coating and the time-cut fractional order convection diffusion attenuation equation containing the migration of the heavy metal salt under the pore substrate are constructed in the embodiment; the strong retention influence of the upper coating on the nuclide migration process is considered through the heavy metal salt migration control equation under the upper coating and the heavy metal salt migration control equation under the pore matrix, so that the accurate simulation of the heavy metal salt migration process in the deep geological layer and the accurate and effective evaluation of the environmental influence degree of the nuclide transport in the mining area are realized, and further powerful technical support is provided for pollution control strategies of a target area, such as the safety evaluation of mine soil and peripheral rivers and lakes.
In this embodiment, a method for simulating a migration process of a target substance is provided, which may be used in the above-mentioned notebook, desktop computer, mobile terminal such as mobile phone, tablet computer, etc., where the attribute data of the underground target substance includes landfill leachate attribute data, and the method is applied to simulate a migration process of a landfill leachate as shown in fig. 9, and fig. 10 is a flowchart of a method for simulating a migration process of a target substance according to an embodiment of the present invention, as shown in fig. 10, where the flowchart includes the following steps:
in step S1001, attribute data and environment data of an underground target substance corresponding to the target area are obtained, and an underground target substance migration control equation is constructed based on the attribute data and environment data of the underground target substance corresponding to the target area.
Specifically, the position and geological structure of the refuse landfill to be evaluated, the proportion of clay, the depth of the refuse landfill and the oxygen content of the landfill are obtained, and the anisotropy and the non-uniformity of the matrix are evaluated and used as the basis for the value of model parameters in an underground target substance migration control equation; acquiring main soluble pollutants and main anaerobic microorganisms of a landfill, wherein the attenuation coefficient of the pollutants in an anaerobic environment; and continuously monitoring the leaching liquid migration process of the landfill to obtain leaching liquid migration data of the landfill under the similar geological conditions.
Specifically, when the target area is a landfill, the step S1001 includes:
step S10011, determining a spatial position parameter of the landfill, a concentration of landfill leachate and a flow parameter of the landfill leachate based on the environmental data.
In particular, the environmental data may include, but is not limited to, spatial location parameters of the landfill, concentration of landfill leachate, and flow parameters of the landfill leachate.
Step S10012, determining an attenuation coefficient of landfill leachate in an anaerobic environment based on the attribute data of the underground target substance.
Specifically, the attenuation coefficient of the leaching solution of the landfill in an anaerobic environment can be measured by carrying out experiments on the pollutant rate of the leaching solution decomposed by anaerobic microorganisms.
In step S10013, an initial fractional capacity coefficient, an initial truncated coefficient, an initial fractional order and an initial diffusion coefficient are obtained.
Step S10014, constructing an underground target substance migration control equation in the target area based on the attenuation coefficient of the landfill leachate in the anaerobic environment, the spatial position parameter of the landfill, the concentration of the landfill leachate, the flow parameter of the landfill leachate, the initial fractional capacity coefficient, the initial cutoff coefficient, the initial fractional order and the initial diffusion coefficient.
Specifically, the migration control equation of the underground target substance in the target area is shown in the following relational expression (8):
（8）
wherein,spatial location parameters representing the target area, +.>Representing a time parameter->Indicating the concentration of landfill leachate, < + >>Representing the total duration of the simulation, +.>Representing the initial fractional capacity coefficient, representing the clay content in the underground medium, and the clay content is higher than the initial fractional capacity coefficient>High value, and->Representing the initialCut-off coefficient, representing the speed of the clay adsorption leaching liquor pollutant reaching saturation,/for>Represents the initial fractional order, reflects the adsorption retention degree of clay to pollutants, < >>Indicating the flow parameters of the landfill leachate, <' > for the landfill leachate>Representing the initial diffusion coefficient, +.>The method is characterized in that the method shows the attenuation coefficient of leaching liquor in an anaerobic environment of a landfill, shows the attenuation process of leaching liquor pollutants under the anaerobic condition, which is influenced by anaerobic microorganisms, the parameter can be influenced by the control work of human pollution, the influence of the leaching liquor pollutants on surrounding water bodies can be reduced by improving the content and the activity of microorganisms near the landfill, the influence is usually obtained through experiments, and the method is characterized in that>Representing the sign of the fractional derivative>Short for the truncated fractional derivative +. >Representing the truncated coefficient>Indicating the moment at which the fractional derivative starts memorization.
And step S1002, performing parameter inversion on the underground target substance migration control equation to generate an underground target substance migration simulation model.
Specifically, the parameters of equation (8), fractional order capacity coefficients, are respectively inverted using the landfill experimental data obtained in step 1001Cut-off coefficient->Fractional order->Diffusion coefficient->And generating an underground target substance migration simulation model which is a time cut-off fractional order convection diffusion attenuation equation.
In step S1003, the migration process of the underground target substance is simulated by using the underground target substance migration simulation model, and target substance migration process simulation data is generated.
Specifically, a finite difference method is utilized to solve a time cut fractional order convection diffusion attenuation equation describing a pollutant migration attenuation process of the landfill in the step S1002, so as to obtain concentration distribution of the leachate pollutant after diffusion attenuation; analyzing the influence of the leachate pollutants on the surrounding environment by combining the surrounding river and lake underground water, and if the pollutants in the leachate are naturally decomposed for a certain time and have low-concentration residues, the leachate is considered to be harmless; if the result is high concentration residue, the content and activity of microorganisms in the landfill need to be manually interfered, simulation is performed again, and critical attenuation coefficient value is determined The concentration of the final leaching solution pollutant is low; if a critical attenuation coefficient value is not available +.>The leaching solution pollutants are considered to have great influence on the environment; and forming a research result report and sending the report to an environmental protection department.
According to the simulation method for the migration process of the target substances, due to the non-uniformity and the anisotropy of the underground medium, the migration of pollutants in the underground medium shows an abnormal diffusion phenomenon, the landfill is close to the ground surface, the clay matrixes are more, a strong detention effect is generated on the pollutants, but the detention effect is weakened along with the gradual saturation of the clay, meanwhile, the degradation effect of underground anaerobic microorganisms is realized, the natural attenuation process exists on the leaching solution pollutants, and the physical biological process is comprehensively considered, so that the time cut-off fractional order convection diffusion attenuation model describing the pollutant migration attenuation process of the landfill is constructed; by constructing an underground target substance migration control equation in a target area, the accurate simulation of the landfill leachate transfer process and the accurate and effective evaluation of the environmental influence degree caused by landfill pollutant leachate transfer are realized, so that powerful technical support is provided for pollution control strategies of the target area, such as the safety evaluation and treatment of surrounding soil and surrounding rivers and lakes.
In this embodiment, a device for simulating a migration process of a target substance is further provided, and the device is used to implement the foregoing embodiments and preferred embodiments, which are not described herein. As used below, the term "module" may be a combination of software and/or hardware that implements a predetermined function. While the means described in the following embodiments are preferably implemented in software, implementation in hardware, or a combination of software and hardware, is also possible and contemplated.
The present embodiment provides a simulation apparatus for migration process of target substance, as shown in fig. 11, including:
the construction module 1101 is configured to obtain attribute data and environment data of an underground target substance corresponding to the target area, and construct an underground target substance migration control equation based on the attribute data and environment data of the underground target substance corresponding to the target area.
And the inversion module 1102 is used for carrying out parameter inversion on the underground target substance migration control equation to generate an underground target substance migration simulation model.
The simulation module 1103 is configured to simulate a migration process of the underground target substance by using the underground target substance migration simulation model, and generate target substance migration process simulation data.
In some alternative embodiments, the inversion module 1102 includes:
the first generation submodule is used for carrying out underground target substance migration experiments by utilizing the environmental data of the target area and generating a plurality of groups of underground target substance migration experiment data.
And the second generation submodule is used for inverting parameters in the underground target substance migration control equation by utilizing a plurality of groups of underground target substance migration experimental data to generate an underground target substance migration simulation model.
In some alternative embodiments, when the target area is the first area, the constructing module 1101 includes:
and the first determination submodule is used for determining the vertical space position parameter of the first area, the concentration of heavy metal salt in the direction from the ore layer to the ground surface and the flow parameter of the heavy metal salt in the first area based on the environmental data.
And the second determination submodule is used for determining the attenuation period of the radioactive elements in the heavy metal salt in the first area and the adsorption efficiency of the target plant to the heavy metal salt based on the attribute data of the underground target substance.
And a third determination submodule for determining a first attenuation coefficient based on the attenuation period of the radioactive element in the heavy metal salt in the first region and the adsorption efficiency of the target plant on the heavy metal salt.
The first construction submodule is used for acquiring an initial time derivative order, an initial space derivative order and an initial super-diffusion coefficient, and constructing an underground target substance migration control equation corresponding to the first area based on the first attenuation coefficient, a vertical space position parameter of the first area, the concentration of heavy metal salt in the direction from the mineral seam to the earth surface, a flow parameter of the heavy metal salt in the first area, the initial time derivative order, the initial space derivative order and the initial super-diffusion coefficient.
In some alternative embodiments, the first build sub-module includes a subsurface target material migration control equation corresponding to the first region as shown in the following relationship:
（2）
wherein,first attenuation coefficient, < >>Vertical spatial position parameter representing the first area, < ->Representing a time parameter->Indicating the concentration of heavy metal salts in the direction from the seam to the surface,/->Representing the time hausdorff derivative, +.>Representing the spatial hausdorff derivative,/>Initial time derivative order, ++>Initial spatial derivative order, ++>Flow parameters of heavy metal salts in the first region, < >>Representing the initial superdiffusion coefficient.
In some alternative embodiments, when the target area is the second area, the constructing module 1101 includes:
And the fourth determining submodule is used for determining a vertical space position parameter of the second area, a seam position parameter of the second area, a pressure parameter of the second area, the concentration of the heavy metal salt in the direction from the disposal warehouse to the ground surface, the concentration of the radioactive heavy metal salt of the heavy metal salt solute in the direction from the seam to the ground surface, a flow parameter of heavy metal salt gas migration in the second area and a flow parameter of heavy metal salt solute migration in the second area based on the environmental data.
And a fifth determination submodule for determining a second attenuation coefficient based on the attribute data of the underground target substance.
The second construction submodule is used for acquiring the time fractional order of the initial heavy metal salt gas migration, the space fractional order of the initial heavy metal salt gas migration and the diffusion coefficient of the initial heavy metal salt gas migration, and constructing a control equation of the heavy metal salt gas migration corresponding to the second area based on the second attenuation coefficient, the vertical space position parameter of the second area, the concentration of the heavy metal salt gas in the direction from the treatment warehouse to the ground surface, the flow parameter of the heavy metal salt gas migration in the second area, the time fractional order of the initial heavy metal salt gas migration, the space fractional order of the initial heavy metal salt gas migration and the flow parameter of the heavy metal salt gas migration.
And the third construction submodule is used for acquiring the time fractional order of initial heavy metal salt solute migration, the space fractional order of initial heavy metal salt solute migration and the diffusion coefficient of initial heavy metal salt solute migration, and constructing a control equation of heavy metal salt solute migration corresponding to the second region based on the second attenuation coefficient, the vertical space position parameter of the second region, the ore layer position parameter of the second region, the pressure parameter of the second region, the concentration of heavy metal salt solute radioactive heavy metal salt in the direction from the ore layer to the earth surface, the flow parameter of heavy metal salt solute migration in the second region, the time fractional order of initial heavy metal salt solute migration, the space fractional order of initial heavy metal salt solute migration and the flow parameter of heavy metal salt solute migration.
And the fourth construction submodule is used for constructing an underground target substance migration control equation corresponding to the second region based on the control equation of heavy metal salt gas migration and the control equation of heavy metal salt solute migration.
In some alternative embodiments, the fourth construction sub-module includes a subsurface target material migration control equation corresponding to the second region as shown in the following relationship:
（4）
（5）
wherein, Vertical spatial position parameter representing the second area, < >>Representing a time parameter->Concentration of gaseous heavy metal salt in the direction from the disposal reservoir to the surface,/->Concentration of heavy metal salt solute radioactive heavy metal salt in direction from ore layer to ground surface, < + >>Representing the pressure parameter of the second region, +.>，/>，/>Representing the time fractional derivative, whereinTime fractional order of initial heavy metal salt gas migration, +.>Time fractional order of initial heavy metal salt solute migration, +.>Spatial fractional order indicative of initial heavy metal salt gas migrationOrder of->Spatial fractional derivative representing Riemann-LiuVil type, ++>The derivative representing this type is the derivative of the Riemann-LiuVil class score order,/I>A seam location parameter indicative of a second zone, < + >>Spatial fractional order of initial heavy metal salt solute migration, +.>Flow parameter indicative of heavy metal salt gas migration in the second zone->A flow parameter indicative of heavy metal salt solute transport in the second zone,/->Diffusion coefficient indicating initial heavy metal salt gas migration, +.>Diffusion coefficient indicating initial heavy metal salt solute migration, +.>Representing a second attenuation coefficient.
In some alternative embodiments, when the target region is the third region, the constructing module 1101 includes:
And a sixth determining submodule, configured to determine a vertical spatial position parameter of the third area, a concentration of heavy metal salt in a direction from the ore layer to the earth surface, a flow parameter of heavy metal salt under the upper cladding layer, and a flow parameter of heavy metal salt under the pore matrix based on the environmental data.
A seventh determination submodule determines a third attenuation coefficient based on attribute data of the subsurface target material.
And a fifth construction submodule, configured to obtain an initial time cutoff fractional order, an initial cutoff coefficient and a diffusion coefficient of heavy metal salt under the initial upper cladding, and construct a heavy metal salt migration control equation under the upper cladding corresponding to the third region based on the third attenuation coefficient, the vertical spatial position parameter of the third region, the concentration of heavy metal salt in the direction from the ore layer to the ground surface, the flow parameter of heavy metal salt under the upper cladding, the initial time cutoff fractional order, the initial cutoff coefficient and the diffusion coefficient of heavy metal salt under the initial upper cladding.
And a sixth construction submodule, configured to obtain an initial time fractional order and a diffusion coefficient of heavy metal salt under the initial pore matrix, and construct a pore matrix heavy metal salt migration control equation corresponding to the third region based on the third attenuation coefficient, the vertical spatial position parameter of the third region, the concentration of heavy metal salt in the direction from the mineral seam to the earth surface, the flow parameter of heavy metal salt under the pore matrix, the initial time fractional order and the diffusion coefficient of heavy metal salt under the initial pore matrix.
And a seventh construction submodule for constructing an underground target substance migration control equation corresponding to the third region based on the heavy metal salt migration control equation under the upper cladding and the heavy metal salt migration control equation under the pore matrix.
In some alternative embodiments, the seventh building sub-module includes a subsurface target material migration control equation corresponding to the third region as shown in the following relationship:
（6）
（7）
wherein,vertical spatial position parameter representing third area, < ->Representing a time parameter->Concentration of heavy metal salts in the direction from the seam to the surface,/->Indicating seam location,/->) Indicating the thickness of the upper cladding>Representing the boundary of the earth's surface,time cut-off fractional derivative representing the carport type, wherein +.>Representing that the derivative of this type is the truncated fractional derivative of the carport class, 0 represents the initial moment,/->Initial time cut-off fractional order, +.>Representing the initial truncated coefficient of the block,time fractional derivative representing the type of carport, wherein +.>Representing that the derivative of this type is the fractional derivative of the carport class, 0 represents the initial moment, ++>Initial time fraction order, +.>Representing an overlay containingThe flow parameters of the heavy metal salt under the layer,indicating the flow parameters of the heavy metal salt under the pore matrix, < + > >Indicating the diffusion coefficient of the heavy metal salt under the initial coating,representing the diffusion coefficient of the heavy metal salt under the original pore matrix, < >>Representing a third attenuation coefficient.
In some alternative embodiments, when the target area is a landfill, the constructing module 1101 includes:
and an eighth determination submodule for determining a spatial position parameter of the landfill, a concentration of leachate of the landfill and a flow parameter of the leachate of the landfill based on the environmental data.
And a ninth determination submodule, configured to determine an attenuation coefficient of landfill leachate in the anaerobic environment based on the attribute data of the underground target material.
And the eighth construction submodule is used for acquiring an initial fractional order capacity coefficient, an initial cut-off coefficient, an initial fractional order and an initial diffusion coefficient, and constructing an underground target substance migration control equation in a target area based on the attenuation coefficient of landfill leachate in an anaerobic environment, the spatial position parameter of the landfill, the concentration of the landfill leachate, the flow parameter of the landfill leachate, the initial fractional order capacity coefficient, the initial cut-off coefficient, the initial fractional order and the initial diffusion coefficient.
In some alternative embodiments, the eighth build sub-module includes a subsurface target material migration control equation in the target area as shown in the following relationship:
（8）
Wherein,spatial location parameters representing the target area, +.>Representing a time parameter->Concentration of landfill leachate, +.>Representing the total duration of the simulation, +.>Initial fractional capacity coefficient, ++>Representing the initial truncated coefficient of the block,initial fractional order, ++>Flow parameters of landfill leachate, < ->The initial diffusion coefficient is set to be,attenuation coefficient of landfill leachate in anaerobic environment, < ->Representing the sign of the fractional derivative>Short for the truncated fractional derivative +.>Representing the truncated coefficient>Indicating the moment at which the fractional derivative starts memorization.
In some alternative embodiments, the analog module 1103 includes:
the solving submodule is used for obtaining the initial concentration of the underground target substance, solving the underground target substance migration simulation model by utilizing a finite difference method based on the initial concentration of the underground target substance, and generating the time-space distribution of the underground target substance in the target area.
In some alternative embodiments, further comprising:
and the evaluation module is used for evaluating the migration process of the underground target substance by using the simulation data and determining the environmental pollution control strategy of the target area based on the evaluation result.
Further functional descriptions of the above respective modules and units are the same as those of the above corresponding embodiments, and are not repeated here.
The simulation of the migration process of the target substance in this embodiment is presented in the form of functional units, where the units refer to ASIC (Application Specific Integrated Circuit ) circuits, processors and memories executing one or more software or fixed programs, and/or other devices that can provide the above-described functionality.
The embodiment of the invention also provides computer equipment, which is provided with the simulation device of the migration process of the target substance shown in the figure 11.
Referring to fig. 12, fig. 12 is a schematic structural diagram of a computer device according to an alternative embodiment of the present invention, as shown in fig. 12, the computer device includes: one or more processors 10, memory 20, and interfaces for connecting the various components, including high-speed interfaces and low-speed interfaces. The various components are communicatively coupled to each other using different buses and may be mounted on a common motherboard or in other manners as desired. The processor may process instructions executing within the computer device, including instructions stored in or on memory to display graphical information of the GUI on an external input/output device, such as a display device coupled to the interface. In some alternative embodiments, multiple processors and/or multiple buses may be used, if desired, along with multiple memories and multiple memories. Also, multiple computer devices may be connected, each providing a portion of the necessary operations (e.g., as a server array, a set of blade servers, or a multiprocessor system). One processor 10 is illustrated in fig. 12.
The processor 10 may be a central processor, a network processor, or a combination thereof. The processor 10 may further include a hardware chip, among others. The hardware chip may be an application specific integrated circuit, a programmable logic device, or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general-purpose array logic, or any combination thereof.
Wherein the memory 20 stores instructions executable by the at least one processor 10 to cause the at least one processor 10 to perform a method for implementing the embodiments described above.
The memory 20 may include a storage program area that may store an operating system, at least one application program required for functions, and a storage data area; the storage data area may store data created according to the use of the computer device, etc. In addition, the memory 20 may include high-speed random access memory, and may also include non-transitory memory, such as at least one magnetic disk storage device, flash memory device, or other non-transitory solid-state storage device. In some alternative embodiments, memory 20 may optionally include memory located remotely from processor 10, which may be connected to the computer device via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
Memory 20 may include volatile memory, such as random access memory; the memory may also include non-volatile memory, such as flash memory, hard disk, or solid state disk; the memory 20 may also comprise a combination of the above types of memories.
The computer device further comprises input means 30 and output means 40. The processor 10, memory 20, input device 30, and output device 40 may be connected by a bus or other means, for example in fig. 12.
The input device 30 may receive input numeric or character information and generate key signal inputs related to user settings and function control of the computer apparatus, such as a touch screen, a keypad, a mouse, a trackpad, a touchpad, a pointer stick, one or more mouse buttons, a trackball, a joystick, and the like. The output means 40 may include a display device, auxiliary lighting means (e.g., LEDs), tactile feedback means (e.g., vibration motors), and the like. Such display devices include, but are not limited to, liquid crystal displays, light emitting diodes, displays and plasma displays. In some alternative implementations, the display device may be a touch screen.
The embodiments of the present invention also provide a computer readable storage medium, and the method according to the embodiments of the present invention described above may be implemented in hardware, firmware, or as a computer code which may be recorded on a storage medium, or as original stored in a remote storage medium or a non-transitory machine readable storage medium downloaded through a network and to be stored in a local storage medium, so that the method described herein may be stored on such software process on a storage medium using a general purpose computer, a special purpose processor, or programmable or special purpose hardware. The storage medium can be a magnetic disk, an optical disk, a read-only memory, a random access memory, a flash memory, a hard disk, a solid state disk or the like; further, the storage medium may also comprise a combination of memories of the kind described above. It will be appreciated that a computer, processor, microprocessor controller or programmable hardware includes a storage element that can store or receive software or computer code that, when accessed and executed by the computer, processor or hardware, implements the methods illustrated by the above embodiments.
Although embodiments of the present invention have been described in connection with the accompanying drawings, various modifications and variations may be made by those skilled in the art without departing from the spirit and scope of the invention, and such modifications and variations fall within the scope of the invention as defined by the appended claims.
Claims (15)
1. A method of simulating a migration process of a target substance, the method comprising:
acquiring attribute data and environment data of an underground target substance corresponding to a target area, and constructing an underground target substance migration control equation based on the attribute data and the environment data of the underground target substance corresponding to the target area;
performing parameter inversion on the underground target substance migration control equation to generate an underground target substance migration simulation model;
and simulating the migration process of the underground target substance by using the underground target substance migration simulation model to generate target substance migration process simulation data.
2. The method of claim 1, wherein the parametric inversion of the subsurface target material migration control equation generates a subsurface target material migration simulation model, comprising:
carrying out underground target substance migration experiments by utilizing the environmental data of the target area to generate a plurality of groups of underground target substance migration experiment data;
and inverting parameters in the underground target substance migration control equation by utilizing the plurality of groups of underground target substance migration experimental data to generate an underground target substance migration simulation model.
3. The method of claim 1, wherein constructing a subsurface target substance migration control equation based on the attribute data and the environmental data of the subsurface target substance corresponding to the target region comprises: when the target area is the first area,
determining a vertical spatial position parameter of the first area, the concentration of heavy metal salt in the direction from the ore layer to the ground surface and a flow parameter of the heavy metal salt in the first area based on the environmental data;
determining the attenuation period of the radioactive element in the heavy metal salt in the first area and the adsorption efficiency of the target plant on the heavy metal salt based on the attribute data of the underground target substance;
determining a first attenuation coefficient based on the attenuation period of the radioactive element in the heavy metal salt in the first region and the adsorption efficiency of the target plant on the heavy metal salt;
and acquiring an initial time derivative order, an initial space derivative order and an initial super-diffusion coefficient, and constructing an underground target substance migration control equation corresponding to a first area based on the first attenuation coefficient, the vertical space position parameter of the first area, the concentration of heavy metal salt in the direction from the ore layer to the earth surface, the flow parameter of the heavy metal salt in the first area, the initial time derivative order, the initial space derivative order and the initial super-diffusion coefficient.
4. The method of claim 3, wherein the constructing a subsurface target material migration control equation corresponding to a first zone based on the first attenuation coefficient, the vertical spatial location of the first zone, the concentration of heavy metal salts in the mineral seam to surface direction, the flow parameters of heavy metal salts in the first zone, the initial time derivative order, the initial spatial derivative order, and the super diffusion coefficient, wherein the subsurface target material migration control equation corresponding to the first zone is represented by the following relationship:
wherein,representing a first attenuation coefficient, ">Vertical spatial position parameter representing the first area, < ->A time parameter is represented by a time parameter,/>indicating the concentration of heavy metal salts in the direction from the seam to the surface,/->Representing the time hausdorff derivative, +.>Representing the spatial hausdorff derivative,/>Representing the initial time derivative order,/->Representing the initial spatial derivative order,/->Indicating the flow parameters of the heavy metal salt in the first zone, < >>Representing the initial superdiffusion coefficient.
5. The method of claim 1, wherein constructing a subsurface target substance migration control equation based on the attribute data and the environmental data of the subsurface target substance corresponding to the target region comprises: when the target area is the second area,
Determining a vertical spatial position parameter of the second area, a seam position parameter of the second area, a pressure parameter of the second area, a concentration of a gas heavy metal salt in a direction from a disposal warehouse to the ground surface, a concentration of a heavy metal salt solute radioactive heavy metal salt in a direction from the seam to the ground surface, a flow parameter of heavy metal salt gas migration in the second area and a flow parameter of heavy metal salt solute migration in the second area based on the environmental data;
determining a second attenuation coefficient based on attribute data of the subsurface target material;
acquiring the time fractional order of initial heavy metal salt gas migration, the space fractional order of initial heavy metal salt gas migration and the diffusion coefficient of initial heavy metal salt gas migration;
constructing a control equation of heavy metal salt gas migration corresponding to a second region based on the second attenuation coefficient, the vertical spatial position parameter of the second region, the concentration of the heavy metal salt gas in the direction from the treatment warehouse to the ground surface, the flow parameter of heavy metal salt gas migration in the second region, the time fractional order of the initial heavy metal salt gas migration, the spatial fractional order of the initial heavy metal salt gas migration and the flow parameter of the heavy metal salt gas migration;
Acquiring a time fractional order of initial heavy metal salt solute migration, a space fractional order of initial heavy metal salt solute migration and a diffusion coefficient of initial heavy metal salt solute migration;
constructing a control equation of heavy metal salt solute migration corresponding to a second region based on the second attenuation coefficient, the vertical spatial position parameter of the second region, the seam position parameter of the second region, the pressure parameter of the second region, the concentration of heavy metal salt solute radioactive heavy metal salt in the direction from the seam to the earth surface, the flow parameter of heavy metal salt solute migration in the second region, the time fractional order of initial heavy metal salt solute migration, the spatial fractional order of initial heavy metal salt solute migration and the flow parameter of heavy metal salt solute migration;
and constructing an underground target substance migration control equation corresponding to the second area based on the heavy metal salt gas migration control equation and the heavy metal salt solute migration control equation.
6. The method of claim 5, wherein the controlling equation for migration of the heavy metal salt gas and the controlling equation for migration of the heavy metal salt solute are based on constructing an underground target substance migration controlling equation corresponding to a second region, and the underground target substance migration controlling equation corresponding to the second region is represented by the following relational expression:
Wherein,vertical spatial position parameter representing the second area, < >>Representing a time parameter->Indicating the concentration of gaseous heavy metal salts in the direction from the treatment reservoir to the surface,/->Represents the concentration of heavy metal salt solute radioactive heavy metal salt in the direction from the ore deposit to the ground surface,representing the pressure parameter of the second region, +.>，/>，/>Representing the time fractional derivative, wherein +.>Time fractional order representing initial heavy metal salt gas migration, +.>Represents the time fractional order of initial heavy metal salt solute migration,spatial fractional order indicative of initial heavy metal salt gas migration, +.>Spatial fractional derivative representing Riemann-LiuVil type, ++>The derivative representing this type is the derivative of the Riemann-LiuVil class score order,/I>A seam location parameter indicative of a second zone, < + >>Spatial fractional order indicative of initial heavy metal salt solute migration,/->Flow parameter indicative of heavy metal salt gas migration in the second zone->A flow parameter indicative of heavy metal salt solute transport in the second zone,/->Diffusion coefficient indicating initial heavy metal salt gas migration, +.>Diffusion coefficient indicating initial heavy metal salt solute migration, +.>Representing a second attenuation coefficient。
7. The method of claim 1, wherein constructing a subsurface target substance migration control equation based on the attribute data and the environmental data of the subsurface target substance corresponding to the target region comprises: when the target area is the third area,
Determining a vertical spatial position parameter of a third area, the concentration of heavy metal salt in the direction from the ore layer to the ground surface, a flow parameter of heavy metal salt under the upper coating layer and a flow parameter of heavy metal salt under the pore matrix based on the environmental data;
determining a third attenuation coefficient based on attribute data of the subsurface target material;
acquiring an initial time cut-off fractional order, an initial cut-off coefficient and an initial diffusion coefficient containing heavy metal salt under an upper coating layer;
constructing an upper-cladding-containing lower heavy metal salt migration control equation corresponding to a third region based on the third attenuation coefficient, the vertical space position parameter of the third region, the concentration of the heavy metal salt in the direction from the ore layer to the ground surface, the flow parameter of the upper-cladding-containing lower heavy metal salt, the initial time cut-off fractional order, the initial cut-off coefficient and the initial upper-cladding-containing lower heavy metal salt diffusion coefficient;
acquiring an initial time fractional order and a diffusion coefficient of heavy metal salt under an initial pore matrix;
constructing a pore matrix heavy metal salt migration control equation corresponding to a third region based on the third attenuation coefficient, the vertical space position parameter of the third region, the concentration of heavy metal salt in the direction from the ore layer to the earth surface, the flow parameter of heavy metal salt under the pore matrix, the initial time fractional order and the diffusion coefficient of heavy metal salt under the initial pore matrix;
And constructing an underground target substance migration control equation corresponding to the third region based on the heavy metal salt migration control equation under the upper coating and the heavy metal salt migration control equation under the pore matrix.
8. The method of claim 7, wherein the constructing an underground target substance migration control equation corresponding to a third zone based on the under-coating heavy metal salt migration control equation and the under-pore-matrix heavy metal salt migration control equation, wherein the underground target substance migration control equation corresponding to the third zone is represented by the following relationship:
wherein,vertical spatial position parameter representing third area, < ->Representing a time parameter->Indicating the concentration of heavy metal salts in the direction from the seam to the surface,/->Indicating seam location,/->) Indicating the thickness of the upper cladding>Representing the boundary of the earth's surface->Time cut-off fractional derivative representing the carport type, wherein +.>Representing that the derivative of this type is the truncated fractional derivative of the carport class, 0 represents the initial moment,/->Representing the order of the initial time cut-off fraction, < >>Representing the initial truncated coefficient, ++>Time fractional derivative representing the type of carport, wherein +.>Representing that the derivative of this type is the fractional derivative of the carport class, 0 represents the initial moment, ++ >Representing the initial time fraction order, +.>Indicating the flow parameters of the heavy metal salt under the upper coating, < >>Indicating the flow parameters of the heavy metal salt under the pore matrix, < + >>Indicating the diffusion coefficient of the heavy metal salt under the initial upper coating layer,/-, for the first coating layer>Representing the diffusion coefficient of the heavy metal salt under the original pore matrix, < >>Representing a third attenuation coefficient.
9. The method of claim 1, wherein constructing a subsurface target substance migration control equation based on the attribute data and the environmental data of the subsurface target substance corresponding to the target region comprises: when the target area is a landfill site,
determining space position parameters of the landfill, concentration of landfill leachate and flow parameters of the landfill leachate based on the environmental data;
determining the attenuation coefficient of landfill leachate in an anaerobic environment based on the attribute data of the underground target substances;
acquiring an initial fractional order capacity coefficient, an initial cut-off coefficient, an initial fractional order and an initial diffusion coefficient;
and constructing an underground target substance migration control equation in a target area based on the attenuation coefficient of the landfill leachate in an anaerobic environment, the spatial position parameter of the landfill, the concentration of the landfill leachate, the flow parameter of the landfill leachate, the initial fractional capacity coefficient, the initial cutoff coefficient, the initial fractional order and the initial diffusion coefficient.
10. The method of claim 9, wherein the constructing a target zone underground target material migration control equation based on the landfill leachate attenuation coefficient in an anaerobic environment, the landfill spatial location parameter, the landfill leachate concentration, the landfill leachate flow parameter, the initial fractional capacity coefficient, the initial cutoff coefficient, the initial fractional order, and the initial diffusion coefficient, wherein the target zone underground target material migration control equation is represented by the following relationship:
wherein,spatial location parameters representing the target area, +.>Representing a time parameter->Indicating the concentration of landfill leachate, < + >>Representing the total duration of the simulation, +.>Representing the initial fractional capacity coefficient, +.>Representing the initial truncated coefficient, ++>Representing the initial fractional order, +.>Indicating the flow parameters of the landfill leachate, <' > for the landfill leachate>Representing the initial diffusion coefficient, +.>Representing the attenuation coefficient of landfill leachate in anaerobic environment, < + >>Representing the sign of the fractional derivative>Short for the truncated fractional derivative +.>Representing the truncated coefficient>Indicating the moment at which the fractional derivative starts memorization.
11. The method of claim 1, wherein simulating the migration process of the subsurface target material using the subsurface target material migration simulation model to generate simulation data comprises:
and obtaining the initial concentration of the underground target substance, solving the underground target substance migration simulation model by using a finite difference method based on the initial concentration of the underground target substance, and generating the time-space distribution of the underground target substance in the target area.
12. The method as recited in claim 1, further comprising:
and evaluating the migration process of the underground target substance by using the simulation data, and determining the environmental pollution control strategy of the target area based on the evaluation result.
13. A simulation apparatus for a migration process of a target substance, the apparatus comprising:
the construction module is used for acquiring attribute data and environment data of the underground target substances corresponding to the target area and constructing an underground target substance migration control equation based on the attribute data and the environment data of the underground target substances corresponding to the target area;
the inversion module is used for carrying out parameter inversion on the underground target substance migration control equation to generate an underground target substance migration simulation model;
And the simulation module is used for simulating the migration process of the underground target substance by using the underground target substance migration simulation model and generating target substance migration process simulation data.
14. A computer device, comprising:
a memory and a processor in communication with each other, the memory having stored therein computer instructions which, upon execution, perform the method of simulating a migration process of a target substance according to any one of claims 1 to 12.
15. A computer-readable storage medium having stored thereon computer instructions for causing a computer to perform the simulation method of the target substance migration process according to any one of claims 1 to 12.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311309905.2A CN117059185B (en) | 2023-10-11 | 2023-10-11 | Simulation method, device, equipment and medium for migration process of target substance |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311309905.2A CN117059185B (en) | 2023-10-11 | 2023-10-11 | Simulation method, device, equipment and medium for migration process of target substance |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117059185A true CN117059185A (en) | 2023-11-14 |
| CN117059185B CN117059185B (en) | 2024-01-26 |
Family
ID=88669593
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311309905.2A Active CN117059185B (en) | 2023-10-11 | 2023-10-11 | Simulation method, device, equipment and medium for migration process of target substance |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117059185B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190085525A1 (en) * | 2016-12-14 | 2019-03-21 | Hohai University | Simple method for determining breakthrough time of anti-seepage liners in landfill |
| US20210164345A1 (en) * | 2019-05-06 | 2021-06-03 | Xi'an Huaxian Petroleum Technology Co., Ltd | A Flow Simulation and Transient Well Analysis Method Based on Generalized Tube Flow and Percolation Coupling |
| CN113919141A (en) * | 2021-09-22 | 2022-01-11 | 中国矿业大学 | Coal mine area storage yard soil heavy metal risk management and control system and migration inversion method |
| CN115048883A (en) * | 2022-07-27 | 2022-09-13 | 中国长江三峡集团有限公司 | Dynamic simulation method and device for soil-vegetation soluble colloid migration process |
-
2023
- 2023-10-11 CN CN202311309905.2A patent/CN117059185B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190085525A1 (en) * | 2016-12-14 | 2019-03-21 | Hohai University | Simple method for determining breakthrough time of anti-seepage liners in landfill |
| US20210164345A1 (en) * | 2019-05-06 | 2021-06-03 | Xi'an Huaxian Petroleum Technology Co., Ltd | A Flow Simulation and Transient Well Analysis Method Based on Generalized Tube Flow and Percolation Coupling |
| CN113919141A (en) * | 2021-09-22 | 2022-01-11 | 中国矿业大学 | Coal mine area storage yard soil heavy metal risk management and control system and migration inversion method |
| CN115048883A (en) * | 2022-07-27 | 2022-09-13 | 中国长江三峡集团有限公司 | Dynamic simulation method and device for soil-vegetation soluble colloid migration process |
Non-Patent Citations (1)
| Title |
|---|
| 万朔阳 等: "基于Hydrus-1D对西坝镇农田土壤重金属迁移模拟及空间解析", 科学技术与工程, vol. 20, no. 2, pages 854 - 859 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117059185B (en) | 2024-01-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rolle et al. | Nernst‐Planck‐based description of transport, coulombic interactions, and geochemical reactions in porous media: Modeling approach and benchmark experiments | |
| Li et al. | Analytical solutions for advective–dispersive solute transport in double‐layered finite porous media | |
| Pu et al. | Assessment of consolidation-induced contaminant transport for compacted clay liner systems | |
| Li et al. | Gas pressure model for layered municipal solid waste landfills | |
| Finsterle et al. | Sealing of a deep horizontal borehole repository for nuclear waste | |
| Dutova et al. | The HydroGeo software package and its usage | |
| Wendling et al. | Gas transport modelling at different spatial scales of a geological repository in clay host rock | |
| Wei et al. | A COMSOL‐PHREEQC Coupled Python Framework for Reactive Transport Modeling in Soil and Groundwater | |
| Montarnal et al. | Presentation and use of a reactive transport code in porous media | |
| Michel et al. | A novel conceptual framework for long‐term leaching of autochthonous soil particles during transient flow | |
| CN117059185B (en) | Simulation method, device, equipment and medium for migration process of target substance | |
| Yu et al. | Numerical investigation on light non–aqueous phase liquid flow in the vadose zone considering porosity effects on soil hydraulic properties | |
| Healy et al. | VS2DRTI: Simulating heat and reactive solute transport in variably saturated porous media | |
| National Research Council et al. | Science and technology for environmental cleanup at Hanford | |
| Larson | Development of the conceptual models for chemical conditions and hydrology used in the 1996 performance assessment for the Waste Isolation Pilot Plant | |
| Freshley et al. | Advanced Simulation Capability for Environmental Management (ASCEM) Phase II Demonstration | |
| Hudson et al. | Coupled THMC modeling for safety assessment of geological disposal of radioactive wastes: The DECOVALEX project (1992–2015) | |
| Toran et al. | A two‐tiered approach to reactive transport: Application to Sr mobility under variable pH | |
| Keshari | Recent trends in flow through porous media | |
| Jiang et al. | Approximate analytical solution to the axisymmetric model for PVD-enhanced in-situ flushing of double-layered contaminated soils | |
| Çokça et al. | A course delivery method by using a software to predict contaminant migration through a soil layer | |
| Štefaňák et al. | Assessment of the Landfill Barrier System through Numerical Analysis: Rehabilitation and Expansion of Belgrade Landfill Case Study | |
| Alves et al. | Stochastic and deterministic models to evaluate the critical distance of a near surface repository for the disposal of intermediate and low level radioactive wastes | |
| Fryar et al. | Hanford Site Composite Analysis: Comparison of One-and Three-Dimensional Simulation of Vadose Zone Transport and Transfer to Groundwater for Selected Waste Sites | |
| Souza et al. | Review of historic forensic geotechnical engineering |
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 |