CN110826236B - Determination method for soil remediation pollution range - Google Patents
Determination method for soil remediation pollution range Download PDFInfo
- Publication number
- CN110826236B CN110826236B CN201911095336.XA CN201911095336A CN110826236B CN 110826236 B CN110826236 B CN 110826236B CN 201911095336 A CN201911095336 A CN 201911095336A CN 110826236 B CN110826236 B CN 110826236B
- Authority
- CN
- China
- Prior art keywords
- pollution
- partition
- range
- thiessen
- layer
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 239000002689 soil Substances 0.000 title claims abstract description 58
- 238000005067 remediation Methods 0.000 title claims abstract description 28
- 238000005192 partition Methods 0.000 claims abstract description 81
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 78
- 231100000719 pollutant Toxicity 0.000 claims abstract description 78
- 238000005070 sampling Methods 0.000 claims abstract description 49
- 238000004088 simulation Methods 0.000 claims abstract description 21
- 238000000638 solvent extraction Methods 0.000 claims abstract description 8
- 230000001502 supplementing effect Effects 0.000 claims abstract description 7
- 150000001875 compounds Chemical class 0.000 claims abstract description 4
- 238000012545 processing Methods 0.000 claims abstract description 4
- 239000010410 layer Substances 0.000 claims description 77
- 238000009826 distribution Methods 0.000 claims description 48
- 238000010276 construction Methods 0.000 claims description 9
- 238000009412 basement excavation Methods 0.000 claims description 6
- 239000002356 single layer Substances 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 4
- 230000001788 irregular Effects 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000013461 design Methods 0.000 description 10
- 238000011109 contamination Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 239000003802 soil pollutant Substances 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 238000007254 oxidation reaction Methods 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 125000001997 phenyl group Chemical class [H]C1=C([H])C([H])=C(*)C([H])=C1[H] 0.000 description 1
- 125000005575 polycyclic aromatic hydrocarbon group Chemical group 0.000 description 1
- 238000003900 soil pollution Methods 0.000 description 1
- 238000005527 soil sampling Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Algebra (AREA)
- Evolutionary Biology (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Processing Of Solid Wastes (AREA)
Abstract
The invention relates to the technical field of soil remediation and discloses a judging method of a soil remediation pollution range. The invention comprises the following steps: data acquisition and processing; partitioning is performed using Thiessen polygons: correspondingly projecting the pollution point location data on a field ground plane, and carrying out plane partition; performing site layering, and defining a pollution level known by each layer of partition; after each layer of the field is divided into Thiessen polygons, projecting sampling sample points belonging to the corresponding depth layer to a plane, and obtaining the pollution degree in the corresponding Thiessen polygons according to the pollutant data of the sampling points in the corresponding layer; performing Kriging interpolation simulation on pollutants of each layer, and comparing and supplementing unknown partition pollution levels according to interpolation results; counting the amount of the repairing compound; and (5) deriving inflection point coordinates. The invention has high judging accuracy and can effectively avoid the phenomenon of staggered layers at the overlapping position of the upper layer pollution range and the lower layer pollution range.
Description
Technical Field
The invention relates to the technical field of soil remediation, in particular to a judging method of a soil remediation pollution range.
Background
In the design and implementation of soil remediation projects, an important problem to be faced first is the definition of the remediation range. The current main flow range judging mode is to simulate the distribution condition of soil pollutants on each depth layer by adopting an inverse distance interpolation mode according to the sampling point distribution condition in the field, and to respectively define the pollution range of each layer.
However, compared with the actual situation, the pollution range defined by the method is often larger, so that the repair engineering quantity is increased, and excessive repair can be generated on the field; in addition, as each depth layering is independently carried out when dividing the repair boundary, the longitudinal relation among the depth layering is not considered, the overlapping part of the upper and lower layer pollution range boundaries is always staggered and mixed, the phenomenon of staggered layers after the different depth layer boundaries are overlapped on the plane is serious, further technical scheme design and actual construction design are not facilitated, the workload of subsequent work is increased, and certain adverse effects are caused to the links of demonstration, construction, evaluation, acceptance and the like of the repair project.
Disclosure of Invention
The invention provides a method for judging the soil restoration pollution range, which has high judgment accuracy and can effectively avoid the phenomenon of staggered layer at the overlapping position of the upper and lower pollution ranges.
The technical problems to be solved are as follows: the pollution range judged by the existing method is larger than that of the actual situation, excessive repair is easy to cause, funds and resources are wasted, serious staggered layers are generated at the boundary of the upper layer and the lower layer, and site construction is difficult.
In order to solve the technical problems, the invention adopts the following technical scheme:
the invention discloses a judging method of a soil remediation pollution range, which comprises the following steps:
step one, data acquisition and processing: extracting information from the collected original site data and carrying out statistical analysis;
step two, partitioning by utilizing Thiessen polygons: correspondingly projecting the pollution point location data on a field ground plane, and carrying out plane partition;
step three, field layering is carried out, and known pollution levels of all layers of subareas are defined; after each layer of the field is divided into Thiessen polygons, projecting sampling sample points belonging to the corresponding depth layer to a plane, and obtaining the pollution degree in the corresponding Thiessen polygons according to the pollutant data of the sampling points in the corresponding layer;
step four, adopting Kriging interpolation simulation to pollutants of each layer, and comparing and supplementing unknown partition pollution levels according to interpolation results;
step five, counting the amount of the repairing compound: counting the repairing formula quantity, the total repairing formula quantity and the pollution partition formula quantity of each layer;
and step six, deriving coordinates of the inflection points, and generating excavation inflection points of the to-be-repaired area for subsequent soil repair construction.
In the second step, the planar partition mode is that Thiessen polygons are divided according to all existing sampling points in the red line range of the field to obtain Thiessen partitions in the field, then the soil is deeply layered according to the longitudinal distribution rule of pollutants, and each deep layer is divided according to the same Thiessen partition.
The invention relates to a judging method of a soil remediation pollution range, and further, the specific division of Thiessen subareas in the step two comprises the following steps:
2.1, importing pollution point data into ArcGIS software in sequence from left to right and from top to bottom in the red line range of the field;
2.2, dividing an irregular triangle network meeting the Delong criterion in all points;
2.3, longitudinally layering the soil to be repaired in the field according to the factors such as the distribution depth of the pollutants, the distribution characteristics of the concentration of the pollutants, the geological structure of the soil, the distribution characteristics of underground water and the like, and designating the depth range of each layer;
and 2.4, carrying out single-layer partition after the longitudinal layered partition is completed, wherein the partition principle is to partition each layer of land according to the Thiessen polygonal network divided by all points.
The invention relates to a judging method of a soil remediation pollution range, which further comprises the following steps of:
3.1, determining the concentration of the pollutants in each Thiessen zone;
according to the Thiessen polygon algorithm, only one sampling point exists in each Thiessen partition, all sample data in the sampling points are taken as the pollutant concentration of the area, and the maximum concentration value of each pollutant in the sampling points is taken as the pollutant concentration of the area;
3.2, determining the pollution level of each Thiessen partition;
judging pollution levels according to the concentration of all pollutants at the sampling point, wherein the pollution levels are specifically classified into pollution-free, light-pollution and heavy-pollution;
and 3.3, comparing the respective pollution levels of all pollutants in the sampling points in each Thiessen partition, and determining the highest value in the pollution levels of all pollutants as the pollution level of the Thiessen partition.
The invention relates to a judging method of a soil remediation pollution range, and further relates to a simulation supplementing method in the step four, which comprises the following steps:
4.1, performing a Kriging interpolation simulation on each pollutant in each plane layering by adopting a Kriging interpolation method, and simulating the concentration distribution gradient of the pollutant in a full field area;
4.2, summarizing the pollutant levels determined in the third step aiming at the simulated concentration gradient distribution of each pollutant in each layer, and dividing each pollutant into light and heavy pollution ranges according to the respective pollution levels; then, overlapping the light gradient range and the heavy gradient range of all pollutants in each plane layering, and fusing the range boundaries with the heaviest pollution levels in the pollutants to obtain the total pollution gradient distribution;
and 4.3, superposing the total pollution gradient distribution map of each layer with the pollution gradient distribution of the Thiessen partition in the step 3.3.
In the method for judging the soil remediation pollution range, in step 4.3, when the total pollution gradient distribution map is compared with the Thiessen subareas obtained in step 3.3 and judging results in two partial subareas are inconsistent, the sampling results of the area and the peripheral area are called for manual comparison judgment.
The invention relates to a method for judging a soil remediation pollution range, which comprises the following steps of further, in the fourth step, sampling points of partial Thiessen subareas at the depth without sampling samples, and lack of sampling data for judging the pollution level in the area, estimating the pollution level in the area by adopting a Kriging interpolation simulation mode, respectively interpolating and superposing target pollutants in each layered area to obtain a total pollution gradient distribution map, and comprehensively judging the pollution level of the Thiessen subareas lacking point position information.
The invention relates to a judging method of a soil remediation pollution range, and further relates to a determining method of each parameter in the step five, which comprises the following steps:
5.1, determining the pollution range of pollutants according to the Thiessen subareas obtained in the second step and the pollutant pollution level distribution of each layer obtained in the fourth step; the contaminated areas above the repair target are considered to be within repair scope;
5.2, obtaining the earthwork quantity in each polluted partition range according to the depth range and the restoration range of each plane layering;
5.3, determining the total polluted earthwork in the polluted level according to the earthwork in each polluted partition range obtained in the step 5.2;
and 5.4, superposing the pollution ranges in each pollution level to obtain a total plane restoration range and total polluted soil earthwork.
The invention relates to a judging method of a soil restoration pollution range, and further, the specific determining method of the earthwork in the step 5.2 is that the area of each pollution zone is multiplied by the thickness of the plane layer.
Compared with the prior art, the judging method of the soil remediation pollution range has the following beneficial effects:
the invention provides a definite soil restoration range judging mode and provides a set of corresponding soil restoration range judging implementation method, and the total soil restoration range and the soil restoration range on each depth layer can be defined by analyzing the detection data, the restoration target value and the depth layer information of the soil sampling points, and the soil restoration range can comprise the information of pollutant types, pollution concentrations and the like, so that the technical scheme and the construction scheme are convenient to carry out process design in the compiling stage.
The method is basically used for partitioning Thiessen polygons, optimizing coordinates of inflection points and avoiding too complex excavation boundaries. And the boundary optimization work is finished when the repair boundary is determined, so that disputes are avoided when the boundary optimization is performed in the later period.
The pollution range is judged by combining Thiessen polygons and Kerling interpolation: based on the Thiessen polygon partition, the coordinates of inflection points are optimized, and the excavation boundary is prevented from being too complex; in the area where the pollution degree cannot be determined through the Thiessen polygon partition, performing pollutant distribution simulation in a Kriging interpolation mode, judging according to the Thiessen polygon inner point position data by adopting a main land block, performing auxiliary judgment on other unknown land blocks in a Kriging interpolation and inverse distance interpolation combined mode, and avoiding excessive determination of the repair quantity; and the simulation precision is enhanced compared with other interpolation methods by combining judgment of the distribution degree of pollutants in the corresponding region.
According to the invention, each layer of soil is partitioned in a Thiessen polygonal mode: different from the traditional pollution range judging method, after judging the range of each layer respectively, the method is divided into separate areas, and finally the serious interleaving situation among the range boundaries of each layer is caused; when judging the soil restoration range, the method adopts a unified Thiessen polygon partition mode to uniformly partition each layer, then judges the pollution degree of each partition in each layer, has high restoration boundary overlap ratio between different layers, and finally obviously improves the staggered layer phenomenon between the obtained layer ranges, thereby being convenient for further process design and construction design.
The following describes a method for judging the soil remediation pollution range according to the present invention with reference to the accompanying drawings.
Drawings
FIG. 1 is a Thiessen partition within a field of the present embodiment;
FIG. 2 is a Thiessen polygonal pollution level of a first soil partition in this example;
FIG. 3 is a graph showing the simulation of pollution zones by Kriging interpolation for the first soil zone in this example;
fig. 4 is a block diagram of the first layer comprehensive simulation result in this embodiment after the Thiessen polygon pollution block is supplemented.
Detailed Description
The method combines the Thiessen polygon and the Kriging interpolation in the pollution range judgment, improves the simulation precision, and can be more suitable for judging the soil pollution range and the soil quantity in practice. Meanwhile, as a uniform partition mode is adopted for each layer in the field partition process, the phenomenon of staggered boundary layer at the overlapping position of the upper layer pollution range and the lower layer pollution range can be effectively avoided, and the subsequent process design and engineering implementation are greatly facilitated. In addition, after the pollution scope is judged in the mode of the invention, the information such as the pollutant types, the pollutant concentrations and the like in the corresponding areas can be reserved in each partition, so that the design of a targeted repair scheme can be conveniently carried out according to the characteristics of each partition in the subsequent process design.
According to the method for judging the soil remediation pollution range, the reasonable field remediation range and pollution partition can be judged through the Thiessen polygon partition and inverse distance difference value and the Kriging interpolation combination simulation mode according to the field investigation result and the field remediation target value, regular partition division is carried out on each depth layer, the occurrence of partition staggered layer phenomenon is avoided, and the operability of partition design is improved.
In this embodiment, a coking site in Shanxi province is taken as a test area, and the site is mainly contaminated by benzene series, polycyclic aromatic hydrocarbon and petroleum hydrocarbon, and is a domestic typical organic pollution site, and the method specifically comprises the following steps:
step one, data acquisition and processing:
the original site data are obtained by sampling soil at different pollution points and performing experimental analysis; extracting detailed information of each pollution point data according to the original site data, and carrying out statistical analysis on the initial data;
the information extracted from each pollution point data includes: sampling point location, sample name, sampling coordinates, sampling depth and concentration of each target pollutant.
The sampling point distribution is based on the following principle: firstly, performing site history analysis, judging suspicious pollution sources, performing point distribution sampling at suspicious positions such as a production factory, a storage workshop and a tank body, then performing encryption point distribution in a whole factory according to a sample detection analysis result, and properly encrypting the point distribution in a pollution area by combining site characteristics according to the point distribution density requirement in site environment investigation technical Specification.
In the example, more than 100 sampling points are distributed, and the distribution principle is that the pollution-free areas in the whole factory are distributed at intervals of 40-80 m, and the pollution areas are distributed at intervals of 20m from the periphery.
Step two, partitioning by utilizing Thiessen polygons: correspondingly projecting the pollution point location data on a field ground plane, and carrying out plane partition;
specifically, the arranged pollution point location data is imported into ArcGIS software or other software with the same function, and then sampling points are projected onto a plane according to sampling coordinates for partitioning.
The planar partition method comprises the steps of carrying out Thiessen polygon division according to all existing sampling points in a red line range of a field to obtain Thiessen partitions in the field range, then carrying out depth layering on soil according to a longitudinal distribution rule of pollutants, and carrying out region division on each depth layering according to the same Thiessen partition.
In this embodiment, the specific division includes the following steps:
2.1, importing pollution point data into ArcGIS software in sequence from left to right and from top to bottom in the red line range of the field;
2.2, dividing an irregular triangular net meeting the Delong criterion (namely, each triangle circumscribing circle does not contain other points) in all points;
2.3, longitudinally layering the soil to be repaired in the field according to the factors such as the distribution depth of the pollutants, the distribution characteristics of the concentration of the pollutants, the geological structure of the soil, the distribution characteristics of underground water and the like, and designating the depth range of each layer;
under the normal condition, the surface soil is relatively heavy in pollution and large in variation range, and the depth range of each layer is 1-2 m; the concentration of the deep soil pollutants is relatively low, the distribution is uniform, and the depth range of each layer can be properly increased, for example, 2-3 m;
and 2.4, carrying out single-layer partition after the longitudinal layered partition is completed, wherein the partition principle is to partition each layer of land parcels according to Thiessen polygonal networks divided by all points, and the planar partition modes of the single layers after the partition by the method are unified, so that the pollution level judgment in the third step is convenient.
Thiessen zones within the field in this embodiment are shown in FIG. 1.
Step three, field layering is carried out, and known pollution levels of all layers of subareas are defined; after each layer of the field is divided into Thiessen polygons, projecting sampling sample points belonging to the corresponding depth layer to a plane, and obtaining the pollution degree in the corresponding Thiessen polygons according to the pollutant data of the sampling points in the corresponding layer;
the specific layering and grading mode specifically comprises the following steps:
3.1, determining the concentration of the pollutants in each Thiessen zone;
according to the Thiessen polygon algorithm, only one sampling point exists in each Thiessen partition, one or more sample data can exist in the sampling point, and when a plurality of sample data exist, the maximum concentration value of each pollutant is taken as the pollutant concentration of the area;
3.2, determining the pollution level of each Thiessen partition;
and judging the pollution level according to the concentration of all pollutants at the sampling point, such as no pollution, light pollution and heavy pollution, wherein the specific level division standard is judged according to the site restoration target value, the pollutant concentration and superscalar multiple, the restoration process treatment difficulty and the restoration process selection.
In this example, the pollution-free criterion is that the pollutant is undetected or the concentration is lower than the repair target; the light pollution judging standard is that the exceeding standard multiple of the pollutant is within 3 times of the restoration target, and the pollutant can be treated to reach the standard by adopting a chemical oxidation process and adding a restoration medicament; the standard exceeding multiple of the severe pollution judgment is more than 3 times, and the severe pollution judgment is repaired by adopting an ectopic thermal desorption process.
And 3.3, comparing the respective pollution levels of all pollutants in the sampling points in each Thiessen partition, and determining the highest value in the pollution levels of all pollutants as the pollution level of the Thiessen partition.
The contamination level of the upper first level partition in this embodiment is shown in fig. 2.
Step four, adopting Kriging interpolation simulation to pollutants of each layer, and comparing and supplementing unknown partition pollution levels according to interpolation results;
because the method is used for avoiding the staggered layer behavior of each layer boundary when carrying out Thiessen polygon partitioning, each layer partitioning principle is based on the spatial relationship of all sampling points; when pollution level determination is carried out in each layer according to the data of the corresponding depth sampling point positions, for partial Thiessen partition sampling point positions, no sample is sampled at the depth, and the situation that sampling data is used for determining the pollution level in the region is lacking; in addition, the pollution degree in the region is estimated by adopting a Kriging interpolation simulation mode, each target pollutant in each layered region is respectively interpolated and overlapped to obtain a total pollution gradient distribution map, and then the pollution level of the Thiessen region lacking point location information is comprehensively judged.
In this embodiment, a specific analog supplementing method includes the following steps:
4.1, performing a Kriging interpolation simulation on each pollutant in each plane layering by adopting a Kriging interpolation method, and simulating the concentration distribution gradient of the pollutant in a full field area;
in this example, there are 12 kinds of site target pollutants, and 5 plane layers are divided, so 12 times of kriging interpolation simulation are performed in each layer, and 60 times of interpolation simulation are performed in total;
4.2, summarizing the pollutant levels determined in the third step aiming at the simulated concentration gradient distribution of each pollutant in each layer, and dividing each pollutant into light and heavy pollution ranges according to the respective pollution levels; then, overlapping the light gradient range and the heavy gradient range of all pollutants in each plane layering, and fusing the range boundaries with the heaviest pollution levels in the pollutants to obtain the total pollution gradient distribution;
4.3, superposing the total pollution gradient distribution map of each layer with the pollution gradient distribution of the Thiessen partition in the step 3.3;
comparing the total pollution gradient distribution map with the Thiessen subareas obtained in the step 3.3, and if two different judging results exist in part of subareas, calling the sampling results of the area and the peripheral area to carry out manual comparison judgment;
4.4, areas in which the pollution level cannot be judged through the step 3.3 due to lack of sampling points possibly exist in the Thiessen subareas, and the areas adopt a total pollution gradient distribution map to carry out auxiliary judgment so as to complement the pollution levels of all Thiessen subareas in the areas;
in this embodiment, the first layer of partition is shown in fig. 3, and the partition after the Thiessen polygon pollution partition is supplemented by the comprehensive simulation result is shown in fig. 4.
Step five, counting the amount of the repairing compound: counting the repairing formula quantity, the total repairing formula quantity and the pollution partition formula quantity of each layer;
the specific statistical method is as follows:
5.1, determining the pollution range of pollutants according to the Thiessen subareas obtained in the second step and the pollutant pollution level distribution of each layer obtained in the fourth step; contaminated areas above the repair target (mild, and severe contamination) are considered to be within repair scope;
5.2, obtaining the earthwork quantity in each polluted partition range according to the depth range and the restoration range of each plane layering;
the specific determination method of the earthwork quantity is that the area of each polluted area is multiplied by the thickness of the plane layer;
taking the first layer as an example in the example, the depth of the first layer ranges from 0 to 2m, the thickness is 2m, and the light pollution area is about 3.5 ten thousand m after the first four steps are judged 2 The area of heavy pollution is about 3.8 ten thousand m 2 The method comprises the steps of carrying out a first treatment on the surface of the Thus the light pollution dosage is about 7 ten thousand meters 3 The area of heavy pollution is about 7.6 ten thousand m 3 ；
5.3, determining the total polluted earthwork in the polluted level according to the earthwork in each polluted partition range obtained in the step 5.2;
in this embodiment, the total contaminated soil volume in the upper first layer of the contaminated soil layer is 14.6 ten thousand meters 3 ；
5.4, superposing the pollution ranges in each pollution level to obtain a total plane restoration range and total polluted soil earthwork quantity;
in this example, the total contaminated soil had an earth volume of about 29 ten thousand meters 3 。
Step six, deriving inflection point coordinates, and generating an excavation inflection point of the region to be repaired for subsequent soil repair construction;
and (3) deriving inflection point coordinates of each process partition through ArcGIS software, and generating an excavation inflection point of the to-be-repaired area for subsequent soil repair construction.
The above examples are only illustrative of the preferred embodiments of the present invention and are not intended to limit the scope of the present invention, and various modifications and improvements made by those skilled in the art to the technical solution of the present invention should fall within the scope of protection defined by the claims of the present invention without departing from the spirit of the present invention.
Claims (8)
1. A judging method for a soil remediation pollution range is characterized by comprising the following steps: the method comprises the following steps:
step one, data acquisition and processing: extracting information from the collected original site data and carrying out statistical analysis;
step two, partitioning by utilizing Thiessen polygons: correspondingly projecting the pollution point location data on a field ground plane, and carrying out plane partition;
step three, field layering is carried out, and known pollution levels of all layers of subareas are defined; after each layer of the field is divided into Thiessen polygons, projecting sampling sample points belonging to the corresponding depth layer to a plane, and obtaining the pollution degree in the corresponding Thiessen polygons according to the pollutant data of the sampling points in the corresponding layer;
step four, adopting Kriging interpolation simulation to pollutants of each layer, and comparing and supplementing unknown partition pollution levels according to interpolation results; for a region with a part of Thiessen partition sampling points and no sampling sample at the depth, lack of sampling data is used for judging the pollution level in the region, the pollution level in the region is estimated by adopting a Kriging interpolation simulation mode, each target pollutant in each layered region is respectively interpolated and overlapped to obtain a total pollution gradient distribution map, and then the pollution level of the Thiessen partition lacking point location information is comprehensively judged;
step five, counting the amount of the repairing compound: counting the repairing formula quantity, the total repairing formula quantity and the pollution partition formula quantity of each layer;
and step six, deriving coordinates of the inflection points, and generating excavation inflection points of the to-be-repaired area for subsequent soil repair construction.
2. The method for determining a soil remediation pollution range according to claim 1, wherein: in the second step, the planar partition mode is that Thiessen polygons are divided according to all existing sampling points in the red line range of the field to obtain Thiessen partitions in the field range, then the soil is deeply layered according to the longitudinal distribution rule of pollutants, and each deep layer is divided according to the same Thiessen partition.
3. The method for determining a soil remediation pollution range according to claim 2, wherein: the specific division of the Thiessen partition in the second step comprises the following steps:
2.1, importing pollution point data into ArcGIS software in sequence from left to right and from top to bottom in the red line range of the field;
2.2, dividing an irregular triangle network meeting the Delong criterion in all points;
2.3, longitudinally layering the soil to be repaired in the field according to the distribution depth of the pollutants, the distribution characteristics of the concentration of the pollutants, the geological structure of the soil and the distribution characteristics of underground water, and designating the depth range of each layer;
and 2.4, carrying out single-layer partition after the longitudinal layered partition is completed, wherein the partition principle is to partition each layer of land according to the Thiessen polygonal network divided by all points.
4. The method for determining a soil remediation pollution range according to claim 1, wherein: the specific layering and grading mode in the third step comprises the following steps:
3.1, determining the concentration of the pollutants in each Thiessen zone;
according to the Thiessen polygon algorithm, only one sampling point exists in each Thiessen partition, all sample data in the sampling points are taken as the pollutant concentration of the area, and the maximum concentration value of each pollutant in the sampling points is taken as the pollutant concentration of the area;
3.2, determining the pollution level of each Thiessen partition;
judging pollution levels according to the concentration of all pollutants at the sampling point, wherein the pollution levels are specifically classified into pollution-free, light-pollution and heavy-pollution;
and 3.3, comparing the respective pollution levels of all pollutants in the sampling points in each Thiessen partition, and determining the highest value in the pollution levels of all pollutants as the pollution level of the Thiessen partition.
5. The method for determining a soil remediation pollution range according to claim 1, wherein: the simulation supplementing method in the fourth step comprises the following steps:
4.1, performing a Kriging interpolation simulation on each pollutant in each plane layering by adopting a Kriging interpolation method, and simulating the concentration distribution gradient of the pollutant in a full field area;
4.2, summarizing the pollutant levels determined in the third step aiming at the simulated concentration gradient distribution of each pollutant in each layer, and dividing each pollutant into light and heavy pollution ranges according to the respective pollution levels; then, overlapping the light gradient range and the heavy gradient range of all pollutants in each plane layering, and fusing the range boundaries with the heaviest pollution levels in the pollutants to obtain the total pollution gradient distribution;
and 4.3, superposing the total pollution gradient distribution map of each layer with the pollution gradient distribution of the Thiessen partition in the step 3.3.
6. The method for determining a soil remediation pollution range of claim 5 wherein: in step 4.3, when the total pollution gradient distribution map is compared with the Thiessen subareas obtained in step 3.3 and the judging results of the two subareas are inconsistent, the sampling results of the area and the peripheral area are called for manual comparison judgment.
7. The method for determining a soil remediation pollution range according to claim 1, wherein: the method for determining each parameter in the fifth step comprises the following steps:
5.1, determining the pollution range of pollutants according to the Thiessen subareas obtained in the second step and the pollutant pollution level distribution of each layer obtained in the fourth step; the contaminated areas above the repair target are considered to be within repair scope;
5.2, obtaining the earthwork quantity in each polluted partition range according to the depth range and the restoration range of each plane layering;
5.3, determining the total polluted earthwork in the polluted level according to the earthwork in each polluted partition range obtained in the step 5.2;
and 5.4, superposing the pollution ranges in each pollution level to obtain a total plane restoration range and total polluted soil earthwork.
8. The method for determining a soil remediation pollution range of claim 7 wherein: the specific method of determining the amount of earth in step 5.2 is the area of each contaminated zone multiplied by the thickness of the planar layer.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911095336.XA CN110826236B (en) | 2019-11-11 | 2019-11-11 | Determination method for soil remediation pollution range |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911095336.XA CN110826236B (en) | 2019-11-11 | 2019-11-11 | Determination method for soil remediation pollution range |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110826236A CN110826236A (en) | 2020-02-21 |
| CN110826236B true CN110826236B (en) | 2023-11-21 |
Family
ID=69553888
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911095336.XA Active CN110826236B (en) | 2019-11-11 | 2019-11-11 | Determination method for soil remediation pollution range |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110826236B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111476504A (en) * | 2020-04-24 | 2020-07-31 | 中科鼎实环境工程有限公司 | Refined site investigation method based on restoration efficiency |
| CN111707490B (en) * | 2020-06-24 | 2023-12-26 | 湘潭大学 | Agricultural land soil pollution investigation staged partition sampling method |
| CN111760897B (en) * | 2020-07-07 | 2022-02-01 | 北京建工环境修复股份有限公司 | Operating well arrangement method and system for conduction type in-situ thermal desorption repair area |
| CN112308290B (en) * | 2020-09-30 | 2023-07-28 | 湘潭大学 | Method for investigating and distributing points in polluted site |
| CN112686497A (en) * | 2020-12-05 | 2021-04-20 | 辽宁大学 | Method for quantizing field soil environmental damage identification object based on kriging interpolation |
| CN113176395A (en) * | 2021-04-26 | 2021-07-27 | 北京市环境保护科学研究院 | Method for determining site soil pollutant exposure concentration in layering manner |
| CN113159454A (en) * | 2021-05-18 | 2021-07-23 | 辽宁大学 | Encryption optimization point distribution method for identifying and evaluating soil environment damage physical quantity |
| CN113486553A (en) * | 2021-07-02 | 2021-10-08 | 浙江大学 | Complex equipment reliability analysis method based on Thiessen polygon area division |
| CN113692797A (en) * | 2021-07-23 | 2021-11-26 | 江苏澳墨环境科技有限公司 | Soil pollution remediation method |
| CN115318811B (en) * | 2022-08-11 | 2023-08-15 | 广东省重工建筑设计院有限公司 | Polluted land block utilizing method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105678454A (en) * | 2016-01-04 | 2016-06-15 | 北京师范大学 | Method for determining soil restoration amount of polluted field on the basis of condition simulation |
| CN106446435A (en) * | 2016-10-09 | 2017-02-22 | 北京建工环境修复股份有限公司 | Method for optimizing contaminated soil restoration scope and boundary |
| CN108268977A (en) * | 2018-01-12 | 2018-07-10 | 武汉智博创享科技股份有限公司 | Reduce the method and system that the extremely extra-high value of pollutant in soil remediation influences |
| CN109754182A (en) * | 2018-12-29 | 2019-05-14 | 上海立昌环境工程股份有限公司 | A kind of calculation method and system of contaminated site soil remediation amount |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2645091A1 (en) * | 2007-11-26 | 2009-05-26 | Tanya L. Michailuck | Assessment and remediation process for contaminated sites |
-
2019
- 2019-11-11 CN CN201911095336.XA patent/CN110826236B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105678454A (en) * | 2016-01-04 | 2016-06-15 | 北京师范大学 | Method for determining soil restoration amount of polluted field on the basis of condition simulation |
| CN106446435A (en) * | 2016-10-09 | 2017-02-22 | 北京建工环境修复股份有限公司 | Method for optimizing contaminated soil restoration scope and boundary |
| CN108268977A (en) * | 2018-01-12 | 2018-07-10 | 武汉智博创享科技股份有限公司 | Reduce the method and system that the extremely extra-high value of pollutant in soil remediation influences |
| CN109754182A (en) * | 2018-12-29 | 2019-05-14 | 上海立昌环境工程股份有限公司 | A kind of calculation method and system of contaminated site soil remediation amount |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110826236A (en) | 2020-02-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110826236B (en) | Determination method for soil remediation pollution range | |
| CN107525907B (en) | Multi-objective optimization method for underground water pollution monitoring network | |
| CN102254349B (en) | Method for constructing three-dimensional physical model of sedimentary stratum system by using drilling data | |
| Hugelius et al. | Estimated stocks of circumpolar permafrost carbon with quantified uncertainty ranges and identified data gaps | |
| Zhu et al. | Building 3D solid models of sedimentary stratigraphic systems from borehole data: an automatic method and case studies | |
| CN102609782B (en) | Bian Caibian based on coal-mining subsidence simulation and second ploughing rate answers the method for optimizing on opportunity | |
| CN109754182B (en) | Method and system for calculating soil remediation quantity of contaminated site | |
| Abdeslam et al. | Application of drastic method for determining the vulnerability of an alluvial aquifer: Morsott-El Aouinet north east of Algeria: using arcgis environment | |
| Lelliott et al. | A structured approach to the measurement of uncertainty in 3D geological models | |
| He et al. | Reservoir simulation using smart proxy in SACROC unit-Case study | |
| CN112308290A (en) | Method for investigating and stationing polluted site | |
| CN115828508B (en) | Automatic ground water critique prediction method based on GIS platform | |
| CN113720639B (en) | Method for investigating and supplementing sampling of polluted site | |
| CN112697849B (en) | Dynamic pollution source positioning method | |
| CN111459955A (en) | Three-dimensional geological structure model automatic updating method and system based on GIS platform | |
| Tao et al. | Three-dimensional delineation of soil pollutants at contaminated sites: Progress and prospects | |
| Shaorui et al. | Predictions of overbreak blocks in tunnels based on the wavelet neural network method and the geological statistics theory | |
| Wilschut et al. | Joint history matching of well data and surface subsidence observations using the ensemble Kalman filter: a field study | |
| CN113298419A (en) | Optimal selection method of typical condition underground water pollution plume spatial distribution interpolation technology | |
| Li et al. | Discrimination of effective source rocks and evaluation of the hydrocarbon resource potential in Marsel, Kazakhstan | |
| CN110687603B (en) | Geological modeling method for internal seepage barrier of offshore oilfield reservoir | |
| Wang et al. | A real-time online structure-safety analysis approach consistent with dynamic construction schedule of underground caverns | |
| CN104035985A (en) | Mining method for abnormal data of basic geographic information | |
| CN114974442A (en) | Multiphase extraction and remediation method for organic contaminated site | |
| Sambo et al. | An integrated reservoir modelling and evolutionary algorithm for optimizing field development in a mature fractured reservoir |
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 |