CN110826236A - Method for judging soil remediation pollution range - Google Patents
Method for judging soil remediation pollution range Download PDFInfo
- Publication number
- CN110826236A CN110826236A CN201911095336.XA CN201911095336A CN110826236A CN 110826236 A CN110826236 A CN 110826236A CN 201911095336 A CN201911095336 A CN 201911095336A CN 110826236 A CN110826236 A CN 110826236A
- Authority
- CN
- China
- Prior art keywords
- pollution
- layer
- partition
- thiessen
- range
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 57
- 238000005067 remediation Methods 0.000 title claims abstract description 38
- 238000005192 partition Methods 0.000 claims abstract description 81
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 77
- 231100000719 pollutant Toxicity 0.000 claims abstract description 77
- 238000005070 sampling Methods 0.000 claims abstract description 46
- 238000004088 simulation Methods 0.000 claims abstract description 21
- 238000000638 solvent extraction Methods 0.000 claims abstract description 16
- 230000001502 supplementing effect Effects 0.000 claims abstract description 5
- 238000012545 processing Methods 0.000 claims abstract description 4
- 239000010410 layer Substances 0.000 claims description 85
- 238000009826 distribution Methods 0.000 claims description 46
- 238000010276 construction Methods 0.000 claims description 9
- 238000010586 diagram Methods 0.000 claims description 7
- 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
- 239000013589 supplement Substances 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 6
- 239000003802 soil pollutant Substances 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000004939 coking Methods 0.000 description 1
- 239000000356 contaminant Substances 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000011066 ex-situ storage Methods 0.000 description 1
- 238000011835 investigation 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
- 239000005416 organic matter Substances 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
- 238000012795 verification Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- 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 method for judging a soil remediation pollution range. The invention comprises the following steps: data acquisition and processing; partitioning by utilizing Thiessen polygons: correspondingly projecting the pollution point location data on a field ground plane, and carrying out plane partitioning; carrying out site layering and defining known pollution levels of each layer partition; after the Thiessen polygon partition is completed in each layer of the site, projecting sampling sample point positions belonging to the corresponding depth layer to a plane, and obtaining the pollution degree in the corresponding Thiessen partition according to the sampling point pollutant data in the corresponding layer; performing Krigin interpolation simulation on each layer of pollutants, comparing according to interpolation results, and supplementing unknown partition pollution levels; counting the repairing amount; and deriving inflection point coordinates. The invention has high judgment accuracy and can effectively avoid the layer error phenomenon generated at the overlapping part of the upper and lower pollution ranges.
Description
Technical Field
The invention relates to the technical field of soil remediation, in particular to a method for judging a soil remediation pollution range.
Background
In the design and implementation process of the soil remediation project, an important problem to be faced first is the definition of the remediation range. At present, the mainstream range judgment mode is that according to the sampling point distribution condition in a field, the distribution condition of soil pollutants on each depth layer is simulated by adopting an inverse distance interpolation mode, and the pollution ranges of all layers are respectively defined.
However, the pollution range defined by the method is often larger than the actual situation, so that the repair engineering amount is increased, and the site is possibly over-repaired; in addition, because each depth layer is independently carried out when the repairing boundary is divided, the longitudinal relation between the depth layers is not considered, finally, the overlapping parts of the upper and lower pollution range boundaries are often staggered and mixed, the staggered layer phenomenon is serious after the boundaries of different depth layers are overlapped to a plane, the design of a further technical scheme and the design of actual construction are not facilitated, the workload of subsequent work is increased, and certain adverse effects are caused on links such as the demonstration, construction, evaluation and acceptance of repairing projects.
Disclosure of Invention
The soil remediation pollution range judgment method provided by the invention has high judgment accuracy and can effectively avoid the staggered layer phenomenon generated at the overlapping part of the upper and lower pollution ranges.
The technical problem to be solved is that: the pollution range judged by the existing method is larger than the actual situation, excessive repair is easily caused, capital and resource waste is caused, and the serious layer-staggered phenomenon is generated at the boundary of an upper layer and a lower layer, so that the field construction is difficult.
In order to solve the technical problems, the invention adopts the following technical scheme:
the invention relates to a method for judging the soil remediation pollution range, which comprises the following steps:
step one, data acquisition and processing: extracting information from the collected original field data and performing statistical analysis;
step two, partitioning by utilizing a Thiessen polygon: correspondingly projecting the pollution point location data on a field ground plane, and carrying out plane partitioning;
thirdly, carrying out site layering and defining known pollution levels of partitions of each layer; after the Thiessen polygon partition is completed in each layer of the site, projecting sampling sample point positions belonging to the corresponding depth layer to a plane, and obtaining the pollution degree in the corresponding Thiessen partition according to the sampling point pollutant data in the corresponding layer;
fourthly, performing Krigin interpolation simulation on each layer of pollutants, and comparing and supplementing unknown partition pollution levels according to interpolation results;
step five, counting the repairing formula: counting the repairing amount of each layer, the total repairing amount and the pollution partition amount of each layer;
and sixthly, deriving inflection point coordinates, and generating an excavation inflection point of the area to be repaired for subsequent soil repair construction.
The method for judging the soil remediation pollution range further comprises the step two of carrying out Thiessen polygonal division according to all existing sampling point positions within the range of the red line of the field in a plane division mode to obtain Thiessen sections within the range of the field, then carrying out deep layering on the soil according to the longitudinal distribution rule of pollutants, and carrying out regional division on each deep layering according to the same Thiessen section.
The invention relates to a method for judging the soil remediation pollution range, and further the step two, the specific division of the Thiessen zones comprises the following steps:
2.1, importing the pollution point data into ArcGIS software from left to right and from top to bottom in a field red line range;
2.2, dividing an irregular triangular net meeting the delaunay criterion in all points;
2.3, longitudinally layering the soil to be repaired in the field according to factors such as pollutant distribution depth, pollutant concentration distribution characteristics, soil geological structure and underground water distribution characteristics, and specifying the depth range of each layer;
and 2.4, carrying out single-layer partition after the longitudinal layered partition is finished, wherein the partition principle is that each layer of land is partitioned according to a Thiessen polygonal network partitioned by all points.
The invention relates to a method for judging the soil remediation pollution range, and further comprises the following steps:
3.1, determining the concentration of the pollutants in each Thiessen subarea;
according to the Thiessen polygon algorithm, only one sampling point is arranged in each Thiessen partition, and the maximum concentration value of each pollutant in all sample data in the sampling point is taken as the pollutant concentration of the area;
3.2, determining the pollution level of each Thiessen partition;
judging the pollution level according to the concentration of all pollutants at the sampling point, and specifically classifying the pollution level into no pollution, light pollution and severe 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 each pollution level of the pollutants as the pollution level of the Thiessen partition.
The invention relates to a method for judging the soil remediation pollution range, and further the simulated supplement method in the fourth step comprises the following steps:
4.1, performing once Krigin interpolation simulation on each pollutant in each plane layer by adopting a Krigin interpolation method, and simulating the concentration distribution gradient of the pollutant in the whole field area;
4.2, simulating the concentration gradient distribution simulated by each pollutant in each layer by adopting the pollutant levels determined in the third step, and dividing each pollutant into light and heavy pollution ranges according to the respective pollution levels; then overlapping the light and heavy gradient ranges of all pollutants in each plane layer, and fusing the range boundaries with the heaviest pollution levels in the pollutants to obtain total pollution gradient distribution;
and 4.3, superposing the total pollution gradient distribution diagram of each layer with the pollution gradient distribution of the Thiessen partition in the step 3.3.
The method for judging the soil remediation pollution range further comprises the step 4.3 of calling the sampling results of the area and the peripheral area to perform manual comparison judgment when the total pollution gradient distribution map is compared with the Thiessen subarea obtained in the step 3.3 and the judgment results of the two modes in partial subareas are inconsistent.
The invention relates to a method for judging a soil remediation pollution range, which comprises the following steps of, in the fourth step, estimating the pollution degree in a region by adopting a kriging interpolation simulation mode for the situation that part of Thiessen partition sampling points are not sampled at the depth and lack of sampling data for judging the pollution level in the region, respectively interpolating and overlapping target pollutants in each layered region to obtain a total pollution gradient distribution diagram, and then comprehensively judging the pollution level of the Thiessen partition lacking point information.
The invention relates to a method for judging the soil remediation pollution range, and further comprises the following steps:
5.1, determining the pollution range of pollutants according to the Thiessen partition obtained in the step two and the pollutant pollution level distribution of each layer obtained in the step four; contaminated areas above the remediation target are all considered to be within the remediation range;
5.2, obtaining the earth volume in each pollution subarea range according to the depth range and the restoration range of each plane layering;
5.3, determining the total polluted earthwork amount in the pollution level according to the earthwork amount in each pollution subarea range obtained in the step 5.2;
and 5.4, overlapping the pollution ranges in the pollution levels to obtain the total plane restoration range and the total polluted soil volume.
The invention relates to a method for judging the soil remediation pollution range, and further, in the step 5.2, the concrete determination method of the earth volume is to multiply the area of each pollution subarea by the thickness of the plane layer.
Compared with the prior art, the method for judging the soil remediation pollution range has the following beneficial effects:
the invention provides a definite soil restoration range judgment mode and provides a set of corresponding soil restoration range judgment implementation method, 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 point position, the information such as pollutant types and pollution concentrations can be included, and the technical design is conveniently carried out at the technical scheme and construction scheme establishment stage.
The method is basically divided into Thiessen polygon areas, the coordinates of inflection points are optimized, and the condition that the excavation boundary is too complex is avoided. And determining to repair the boundary and completing the boundary optimization work, thereby avoiding disputes generated when the boundary optimization is carried out at the later stage.
The invention utilizes the combination of Thiessen polygon and Critical interpolation to judge the pollution range: optimizing inflection point coordinates based on Thiessen polygon subareas to avoid excessively complex excavation boundaries; in an area where the pollution degree cannot be determined in the Thiessen polygon partition, performing pollutant distribution simulation in a Critical interpolation mode, determining the main land parcel according to point data in the Thiessen polygon partition, and performing auxiliary determination on the other unknown land parcels in a manner of combining Critical interpolation and reverse distance interpolation to avoid excessive determination of the repair volume; compared with other interpolation methods, the simulation precision is enhanced by combining and judging the distribution degree of the pollutants in the corresponding area.
The soil of each layer is uniformly partitioned in a Thiessen polygonal mode: different from the situation that the boundary staggering condition between the range boundaries of each layer is serious finally because each layer is separately partitioned after the range is respectively determined in the traditional pollution range determination method; when the soil restoration range is judged, the uniform partition is firstly carried out on each layer in a uniform Thiessen polygon partition mode, then the pollution degree of each partition in each layer is judged, the restoration boundary overlap ratio between different layers is high, the finally obtained layer staggering phenomenon between each layer range is obviously improved, and further process design and construction design are facilitated.
The method for determining the soil remediation pollution range according to the present invention will be further described with reference to the accompanying drawings.
Drawings
FIG. 1 is a Thiessen section within the field in this example;
FIG. 2 is a Thiessen polygonal contamination level for the first layer soil partition of the present example;
FIG. 3 is a Kriging interpolation simulation pollution zone of the first layer soil zone in this embodiment;
fig. 4 shows the partition after the first layer comprehensive simulation result supplements the thiessen polygon contamination partition in this embodiment.
Detailed Description
The method combines the Thiessen polygon and the Critical interpolation in the pollution range judgment, improves the simulation precision, and can judge the soil pollution range and the amount in a more practical way. Meanwhile, because a uniform partition mode is adopted for each layer in the field partition process, the staggered layer phenomenon of boundary at the overlapping part of the upper and lower pollution ranges can be effectively avoided, and the follow-up process design and engineering implementation are greatly facilitated. In addition, after the pollution range is judged by the method, the pollutant types, the pollutant concentrations and other information in the corresponding areas can be reserved in each subarea, and the targeted repair scheme design can be conveniently carried out according to the characteristics of each subarea in the subsequent process design.
According to the method for judging the soil remediation pollution range, the reasonable field remediation range and the pollution partition can be judged in a combined simulation mode of Thiessen polygonal partition, reverse distance difference and Crimen interpolation according to the field investigation result and the field remediation target value, regular partition division is performed on each depth layer, the partition layer staggering phenomenon is avoided, and the operability of partition design is improved.
In the embodiment, a coking field in Shanxi province is taken as a test area, main pollutants in the field are benzene series, polycyclic aromatic hydrocarbon and petroleum hydrocarbon, the field is a domestic typical organic matter pollution field, and the field verification test specifically comprises the following steps:
step one, data acquisition and processing:
the original site data is 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 performing statistical analysis on the initial data;
the information extracted for each pollution point data includes: sampling point location, sample name, sampling coordinates, sampling depth and concentration of each target pollutant.
The distribution of sampling points is based on the following principle: the method comprises the steps of firstly, carrying out site history analysis, judging suspicious pollution sources, carrying out point distribution sampling at suspicious positions such as production plants, storage plants, tank bodies and the like, then carrying out encrypted point distribution in the whole plant area according to a sample detection and analysis result, referring to the point distribution density requirement in site environment survey technical specification, and properly encrypting in a polluted area by combining site characteristics.
In the embodiment, more than 100 sampling points are distributed in the whole plant area, and the distribution principle is that points are distributed in the pollution-free area at intervals of 40-80 m and in the pollution area at intervals of 20m from the periphery.
Step two, partitioning by utilizing a Thiessen polygon: correspondingly projecting the pollution point location data on a field ground plane, and carrying out plane partitioning;
specifically, the sorted pollution point location data is imported into ArcGIS software or other software with the same function, and then the sampling point location is projected onto a plane according to sampling coordinates for partitioning.
The idea of plane partitioning is that within a field red line range, Thiessen polygon partitioning is carried out according to all existing sampling point positions to obtain Thiessen partitions within the field range, then, according to a longitudinal distribution rule of pollutants, soil is subjected to deep partitioning, and each deep partition is subjected to region partitioning according to the same Thiessen partition.
Specifically, in this embodiment, the specific division includes the following steps:
2.1, importing the pollution point data into ArcGIS software from left to right and from top to bottom in a field red line range;
2.2, dividing an irregular triangular net meeting the delaunay criterion (namely, each triangle circumcircles does not contain other points) in all the points;
2.3, longitudinally layering the soil to be repaired in the field according to factors such as pollutant distribution depth, pollutant concentration distribution characteristics, soil geological structure and underground water distribution characteristics, and specifying the depth range of each layer;
in general, because the surface soil is relatively heavily polluted and the variation range is large, the depth range of each layer is 1-2 m; the concentration of the deep soil pollutants is relatively low, the deep soil pollutants are distributed uniformly, and the depth range of each layer can be properly increased, such as 2-3 m;
and 2.4, carrying out single-layer partition after the longitudinal layered partition is finished, wherein the partition principle is that each layer of land is partitioned according to the Thiessen polygonal network partitioned by all points, and the plane partition modes of all single layers are unified after the partition is carried out by the method, so that the pollution level judgment in the third step is facilitated.
The Thiessen zones within the field in this embodiment are shown in FIG. 1.
Thirdly, carrying out site layering and defining known pollution levels of partitions of each layer; after the Thiessen polygon partition is completed in each layer of the site, projecting sampling sample point positions belonging to the corresponding depth layer to a plane, and obtaining the pollution degree in the corresponding Thiessen partition according to the sampling point pollutant data 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 subarea;
according to the Thiessen polygon algorithm, only one sampling point exists in each Thiessen partition, one or more sample data possibly 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, wherein the specific level division standard is judged according to the site restoration target value, the pollutant concentration and standard exceeding multiple, the restoration process treatment difficulty and the restoration process selection if no pollution, light pollution and severe pollution exist.
In this example, the non-contamination determination criterion is that no contaminant is detected or the concentration is lower than the repair target; the judgment standard of the light pollution is that the exceeding multiple of the pollutant is within 3 times of the restoration target, and the pollutant can reach the standard by adopting a chemical oxidation process and adding a restoration agent; the standard exceeding multiple of the severe pollution is more than 3 times, and an ex-situ thermal desorption process is required for repairing.
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 each pollution level of the 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.
Fourthly, performing Krigin interpolation simulation on each layer of pollutants, and comparing and supplementing unknown partition pollution levels according to interpolation results;
in the method, when the Thiessen polygon partition is carried out, the layer staggering behavior of each layer boundary is avoided, so that each layer partition principle is based on the spatial relationship of all sampling point positions; when pollution levels are judged according to data of sampling points of corresponding depths in each layer, sampling points of partial Thiessen sub-areas are not sampled at the depth, and sampling data are lacked to judge the condition of the pollution levels in the areas; and as a supplement, the pollution degree in the region is estimated by adopting a kriging interpolation simulation mode, each target pollutant in each layered region is interpolated and superposed to obtain a total pollution gradient distribution diagram, and then the pollution level of the Thiessen region lacking point location information is comprehensively judged.
In this embodiment, a specific simulation supplementing method includes the following steps:
4.1, performing once Krigin interpolation simulation on each pollutant in each plane layer by adopting a Krigin interpolation method, and simulating the concentration distribution gradient of the pollutant in the whole field area;
in the example, 12 field target pollutants are divided into 5 plane layers, so that the kriging interpolation simulation is carried out for 12 times in each layer, and the interpolation simulation is carried out for 60 times in total;
4.2, simulating the concentration gradient distribution simulated by each pollutant in each layer by adopting the pollutant levels determined in the third step, and dividing each pollutant into light and heavy pollution ranges according to the respective pollution levels; then overlapping the light and heavy gradient ranges of all pollutants in each plane layer, and fusing the range boundaries with the heaviest pollution levels in the pollutants to obtain total pollution gradient distribution;
4.3, superposing the total pollution gradient distribution diagram 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 the judgment results of the two modes are different in partial subareas, calling the sampling results of the area and the peripheral area for manual comparison judgment;
4.4, areas which cannot be used for judging the pollution levels through the 3.3 steps due to lack of sampling points possibly exist in the Thiessen zones, the areas are subjected to auxiliary judgment by adopting a total pollution gradient distribution diagram, and the pollution levels of all Thiessen zones in the layer are supplemented;
in the first layer of partitions in this embodiment, the polluted partitions subjected to auxiliary judgment by pollutant difference value simulation are shown in fig. 3, and the partitions supplemented by the comprehensive simulation result for the thiessen polygonal polluted partitions are shown in fig. 4.
Step five, counting the repairing formula: counting the repairing amount of each layer, the total repairing amount and the pollution partition amount of each layer;
the specific statistical method is as follows:
5.1, determining the pollution range of pollutants according to the Thiessen partition obtained in the step two and the pollutant pollution level distribution of each layer obtained in the step four; contaminated areas above the target of remediation (mild, and severe contamination) were considered to be within the scope of remediation;
5.2, obtaining the earth volume in each pollution subarea range according to the depth range and the restoration range of each plane layering;
the specific determination method of the earth volume is that the area of each pollution subarea is multiplied by the thickness of the plane layer;
in this example, the first layer is taken as an example, the depth of the first layer is 0-2 m, the thickness of the first layer is 2m, and the area of the first layer with light pollution is 3.5 ten thousand m after the first four steps2The area of severe pollution is about 3.8 ten thousand meters2(ii) a Therefore, the light pollution amount is about 7 ten thousand meters3The area of severe pollution is about 7.6 ten thousand meters3;
5.3, determining the total polluted earthwork amount in the pollution level according to the earthwork amount in each pollution subarea range obtained in the step 5.2;
in this example, the total contaminated soil volume in the upper first-tier contamination level was 14.6 km3;
5.4, overlapping the pollution ranges in the pollution levels to obtain a total plane restoration range and a total polluted soil volume;
in this example, the total soil volume of the contaminated soil was about 29 km3。
Deriving inflection point coordinates, and generating an excavation inflection point of the area to be repaired for subsequent soil repair construction;
and deriving inflection point coordinates of each process partition through ArcGIS software, and generating an excavation inflection point of the area to be repaired for subsequent soil repair construction.
The above-mentioned embodiments are merely illustrative of the preferred embodiments of the present invention, and do not limit the scope of the present invention, and various modifications and improvements of the technical solution of the present invention by those skilled in the art should fall within the protection scope defined by the claims of the present invention without departing from the spirit of the present invention.
Claims (9)
1. A method for judging the 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 field data and performing statistical analysis;
step two, partitioning by utilizing a Thiessen polygon: correspondingly projecting the pollution point location data on a field ground plane, and carrying out plane partitioning;
thirdly, carrying out site layering and defining known pollution levels of partitions of each layer; after the Thiessen polygon partition is completed in each layer of the site, projecting sampling sample point positions belonging to the corresponding depth layer to a plane, and obtaining the pollution degree in the corresponding Thiessen partition according to the sampling point pollutant data in the corresponding layer;
fourthly, performing Krigin interpolation simulation on each layer of pollutants, and comparing and supplementing unknown partition pollution levels according to interpolation results;
step five, counting the repairing formula: counting the repairing amount of each layer, the total repairing amount and the pollution partition amount of each layer;
and sixthly, deriving inflection point coordinates, and generating an excavation inflection point of the area to be repaired for subsequent soil repair construction.
2. The method for determining the soil remediation pollution zone of claim 1, wherein: and in the second step, the plane partitioning mode is that within the range of the red line of the field, Thiessen polygon partitioning is carried out according to all existing sampling point positions to obtain Thiessen partitions within the range of the field, then, according to the distribution rule of pollutants in the longitudinal direction, deep layering is carried out on the soil, and region partitioning is carried out on each deep layering according to the same Thiessen partition.
3. The method for determining the soil remediation pollution zone of claim 2, wherein: the specific division of the Thiessen subareas in the second step comprises the following steps:
2.1, importing the pollution point data into ArcGIS software from left to right and from top to bottom in a field red line range;
2.2, dividing an irregular triangular net meeting the delaunay criterion in all points;
2.3, longitudinally layering the soil to be repaired in the field according to factors such as pollutant distribution depth, pollutant concentration distribution characteristics, soil geological structure and underground water distribution characteristics, and specifying the depth range of each layer;
and 2.4, carrying out single-layer partition after the longitudinal layered partition is finished, wherein the partition principle is that each layer of land is partitioned according to a Thiessen polygonal network partitioned by all points.
4. The method for determining the soil remediation pollution zone of 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 subarea;
according to the Thiessen polygon algorithm, only one sampling point is arranged in each Thiessen partition, and the maximum concentration value of each pollutant in all sample data in the sampling point is taken as the pollutant concentration of the area;
3.2, determining the pollution level of each Thiessen partition;
judging the pollution level according to the concentration of all pollutants at the sampling point, and specifically classifying the pollution level into no pollution, light pollution and severe 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 each pollution level of the pollutants as the pollution level of the Thiessen partition.
5. The method for determining the soil remediation pollution zone of claim 1, wherein: the analog supplement method in the fourth step comprises the following steps:
4.1, performing once Krigin interpolation simulation on each pollutant in each plane layer by adopting a Krigin interpolation method, and simulating the concentration distribution gradient of the pollutant in the whole field area;
4.2, simulating the concentration gradient distribution simulated by each pollutant in each layer by adopting the pollutant levels determined in the third step, and dividing each pollutant into light and heavy pollution ranges according to the respective pollution levels; then overlapping the light and heavy gradient ranges of all pollutants in each plane layer, and fusing the range boundaries with the heaviest pollution levels in the pollutants to obtain total pollution gradient distribution;
and 4.3, superposing the total pollution gradient distribution diagram of each layer with the pollution gradient distribution of the Thiessen partition in the step 3.3.
6. The method for determining the soil remediation pollution zone of claim 5, wherein: and 4.3, when the total pollution gradient distribution map is compared with the Thiessen subareas obtained in the step 3.3, and judging results of two modes in partial subareas are inconsistent, calling sampling results of the area and the peripheral area for manual comparison and judgment.
7. The method for determining the soil remediation pollution zone of claim 5, wherein: and in the fourth step, for the situation that part of Thiessen partition sampling points have no sampling sample at the depth and lack of sampling data to judge the pollution level in the area, inferring the pollution degree in the area by adopting a Krigin interpolation simulation mode, respectively interpolating and superposing target pollutants in each layered area to obtain a total pollution gradient distribution diagram, and then comprehensively judging the pollution level of the Thiessen partition lacking point position information.
8. The method for determining the soil remediation pollution zone of 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 partition obtained in the step two and the pollutant pollution level distribution of each layer obtained in the step four; contaminated areas above the remediation target are all considered to be within the remediation range;
5.2, obtaining the earth volume in each pollution subarea range according to the depth range and the restoration range of each plane layering;
5.3, determining the total polluted earthwork amount in the pollution level according to the earthwork amount in each pollution subarea range obtained in the step 5.2;
and 5.4, overlapping the pollution ranges in the pollution levels to obtain the total plane restoration range and the total polluted soil volume.
9. The method for determining the soil remediation pollution zone of claim 8, wherein: the specific determination of the amount of earth in step 5.2 is the area of each contaminated partition 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 true CN110826236A (en) | 2020-02-21 |
CN110826236B 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) |
Cited By (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 |
CN111707490A (en) * | 2020-06-24 | 2020-09-25 | 湘潭大学 | Method for staged and zoned sampling of agricultural land soil pollution survey |
CN111760897A (en) * | 2020-07-07 | 2020-10-13 | 北京建工环境修复股份有限公司 | Operating well arrangement method and system for conduction type in-situ thermal desorption repair area |
CN112308290A (en) * | 2020-09-30 | 2021-02-02 | 湘潭大学 | Method for investigating and stationing polluted site |
CN112686497A (en) * | 2020-12-05 | 2021-04-20 | 辽宁大学 | Method for quantizing field soil environmental damage identification object based on kriging interpolation |
CN113159454A (en) * | 2021-05-18 | 2021-07-23 | 辽宁大学 | Encryption optimization point distribution method for identifying and evaluating soil environment damage physical quantity |
CN113176395A (en) * | 2021-04-26 | 2021-07-27 | 北京市环境保护科学研究院 | Method for determining site soil pollutant exposure concentration in layering manner |
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 |
CN115318811A (en) * | 2022-08-11 | 2022-11-11 | 广东省重工建筑设计院有限公司 | Method for utilizing polluted land blocks |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090142137A1 (en) * | 2007-11-26 | 2009-06-04 | Pioneer Professional Services Group Ltd. | Assessment and remediation process for contaminated sites |
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 |
-
2019
- 2019-11-11 CN CN201911095336.XA patent/CN110826236B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090142137A1 (en) * | 2007-11-26 | 2009-06-04 | Pioneer Professional Services Group Ltd. | Assessment and remediation process for contaminated sites |
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 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111476504A (en) * | 2020-04-24 | 2020-07-31 | 中科鼎实环境工程有限公司 | Refined site investigation method based on restoration efficiency |
CN111707490A (en) * | 2020-06-24 | 2020-09-25 | 湘潭大学 | Method for staged and zoned sampling of agricultural land soil pollution survey |
CN111707490B (en) * | 2020-06-24 | 2023-12-26 | 湘潭大学 | Agricultural land soil pollution investigation staged partition sampling method |
CN111760897A (en) * | 2020-07-07 | 2020-10-13 | 北京建工环境修复股份有限公司 | 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 |
CN112308290A (en) * | 2020-09-30 | 2021-02-02 | 湘潭大学 | Method for investigating and stationing 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 |
CN115318811A (en) * | 2022-08-11 | 2022-11-11 | 广东省重工建筑设计院有限公司 | Method for utilizing polluted land blocks |
CN115318811B (en) * | 2022-08-11 | 2023-08-15 | 广东省重工建筑设计院有限公司 | Polluted land block utilizing method |
Also Published As
Publication number | Publication date |
---|---|
CN110826236B (en) | 2023-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110826236A (en) | Method for judging 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 | |
Bianchi et al. | Integrating deterministic lithostratigraphic models in stochastic realizations of subsurface heterogeneity. Impact on predictions of lithology, hydraulic heads and groundwater fluxes | |
CN109754182B (en) | Method and system for calculating soil remediation quantity of contaminated site | |
Alberti et al. | Null-space Monte Carlo particle tracking to assess groundwater PCE (Tetrachloroethene) diffuse pollution in north-eastern Milan functional urban area | |
Lelliott et al. | A structured approach to the measurement of uncertainty in 3D geological models | |
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 | |
CN103345566B (en) | Based on the geochemical anomaly discrimination and evaluation method of Geological Connotation | |
CN106446435A (en) | Method for optimizing contaminated soil restoration scope and boundary | |
Kearsey et al. | Testing the application and limitation of stochastic simulations to predict the lithology of glacial and fluvial deposits in Central Glasgow, UK | |
CN112308290A (en) | Method for investigating and stationing polluted site | |
Ouellon et al. | Hydraulic conductivity heterogeneity of a local deltaic aquifer system from the kriged 3D distribution of hydrofacies from borehole logs, Valcartier, Canada | |
Turner et al. | Improved understanding of groundwater flow in complex superficial deposits using three-dimensional geological-framework and groundwater models: an example from Glasgow, Scotland (UK) | |
CN115828508B (en) | Automatic ground water critique prediction method based on GIS platform | |
CN113720639A (en) | Method for investigation, supplement and sampling of polluted site | |
Collins et al. | Horizontal and radial collector wells: simple tools for a complex problem | |
Colombo et al. | Null-space Monte Carlo particle backtracking to identify groundwater tetrachloroethylene sources | |
CN113298419A (en) | Optimal selection method of typical condition underground water pollution plume spatial distribution interpolation technology | |
Lukjan et al. | Assessment of alluvial aquifer heterogeneity and development of stochastic hydrofacies models for the Hat Yai Basin in Southern Thailand | |
CN114974442A (en) | Multiphase extraction and remediation method for organic contaminated site | |
CN109307690B (en) | Method and device for investigating petroleum hydrocarbon soil pollution and electronic equipment | |
WO2021142895A1 (en) | Design method for sve process parameters in petroleum contamination site | |
DeAngelo et al. | Geothermal play fairway analysis, part 2: GIS methodology | |
CN107729629B (en) | Environmental hydrogeological zoning application method for polluted site |
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 |