CN112686497A - Method for quantizing field soil environmental damage identification object based on kriging interpolation - Google Patents
Method for quantizing field soil environmental damage identification object based on kriging interpolation Download PDFInfo
- Publication number
- CN112686497A CN112686497A CN202011408044.XA CN202011408044A CN112686497A CN 112686497 A CN112686497 A CN 112686497A CN 202011408044 A CN202011408044 A CN 202011408044A CN 112686497 A CN112686497 A CN 112686497A
- Authority
- CN
- China
- Prior art keywords
- pollutant
- concentration
- soil
- characteristic
- standard
- 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.)
- Pending
Links
- 239000002689 soil Substances 0.000 title claims abstract description 58
- 230000007613 environmental effect Effects 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 113
- 231100000719 pollutant Toxicity 0.000 claims abstract description 112
- 238000005070 sampling Methods 0.000 claims abstract description 65
- 238000011002 quantification Methods 0.000 claims abstract description 14
- 238000011156 evaluation Methods 0.000 claims abstract description 12
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 238000012216 screening Methods 0.000 claims abstract description 4
- 238000012544 monitoring process Methods 0.000 claims abstract description 3
- 238000003900 soil pollution Methods 0.000 claims abstract description 3
- 239000002356 single layer Substances 0.000 claims description 21
- 238000011109 contamination Methods 0.000 claims description 18
- 239000010410 layer Substances 0.000 claims description 15
- 230000009466 transformation Effects 0.000 claims description 10
- 239000000356 contaminant Substances 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 6
- 238000007726 management method Methods 0.000 claims description 6
- 238000013139 quantization Methods 0.000 claims description 6
- 238000011835 investigation Methods 0.000 claims description 5
- 238000011160 research Methods 0.000 claims description 4
- 239000002344 surface layer Substances 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims 1
- 238000013461 design Methods 0.000 abstract description 4
- 238000012502 risk assessment Methods 0.000 abstract description 3
- 230000006399 behavior Effects 0.000 description 2
- 238000003912 environmental pollution Methods 0.000 description 2
- FBOUIAKEJMZPQG-AWNIVKPZSA-N (1E)-1-(2,4-dichlorophenyl)-4,4-dimethyl-2-(1,2,4-triazol-1-yl)pent-1-en-3-ol Chemical compound C1=NC=NN1/C(C(O)C(C)(C)C)=C/C1=CC=C(Cl)C=C1Cl FBOUIAKEJMZPQG-AWNIVKPZSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001364 causal effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000013501 data transformation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Landscapes
- Processing Of Solid Wastes (AREA)
Abstract
The invention relates to a field soil environmental damage identification object quantification method based on kriging interpolation. The method collects and analyzes relevant basic information of a field to be quantified, designs a sampling, monitoring and identifying evaluation scheme, and obtains geographic information coordinates of sampling points and soil samples; performing laboratory analysis and test on the soil sample to obtain the concentration of the characteristic pollutant, and establishing a soil pollution database of the characteristic pollutant of the sampling point; determining the standard exceeding times of the concentration of the characteristic pollutants of the sampling point by comparing with the characteristic pollutant screening value; and quantifying the pollution area and the earthwork of the site characteristic pollutants according to a kirschner spatial interpolation method, and determining the physical quantity of the soil environment damage of the polluted site. Providing scientific basis for value quantification, health risk assessment, restoration scheme design, pollution site management and control and safety utilization of the pollution site; and a technical support is provided for identifying the environmental damage of the soil in the polluted site.
Description
Technical Field
The invention relates to the field of identification and evaluation of soil environmental damage of a polluted site, in particular to a method for quantifying an identification object of the soil environmental damage of the site based on Krigin interpolation.
Background
In recent years, site pollution has become one of the important environmental issues of concern in countries around the world. The polluted sites of industrial enterprises are distributed all over the world and have the characteristics of wide distribution, high pollutant concentration, long residual time, great environmental threat and the like. With the rapid development of urbanization, in order to improve the reasonable planning of cities and adjust and move polluted enterprises, a large amount of idle polluted land caused by enterprise movement appears in main large and medium-sized cities in China. The region is a non-negligible important pollution source in cities, and brings serious threats to the surrounding environment and the physical health of residents, and simultaneously causes waste of urban land resources. The problem of the polluted site has attracted extensive attention, and environmental supervision and treatment of the polluted site are urgent.
The reasons and consequences of environmental damage caused by polluted sites are complex, administrative supervision is in short supply, the basis for compensation and punishment of ecological environmental damage is lacked, environmental disputes influence social stability, and the safety bottom line of the ecological environment is threatened. The environmental damage appraisal and evaluation is the important content of ecological environment damage compensation, and the environmental damage appraisal and evaluation specifically refers to the process of appraising and evaluating organizations according to the specified procedures and methods, applying scientific technology and professional knowledge, distinguishing pollutant properties, evaluating the range and degree of environmental damage caused by environmental pollution or ecological damage behaviors, judging the causal relationship between the environmental pollution or ecological damage behaviors and the environmental damage, determining the baseline level of ecological environment recovery, and quantifying the amount of the environmental damage and the damage during compensation. The quantification of the soil environmental damage in the polluted site is an important link for identification and evaluation. The physical quantification specifically refers to the calculation of the degree, range and area of the concentration of the different types of site characteristic pollutants exceeding the screening value.
Therefore, a real object quantification method in the field of identification and evaluation of soil environmental damage of a polluted site is urgently provided. The invention provides a method for quantizing site soil environmental damage identification objects based on kriging interpolation, which provides scientific basis for value quantization of a polluted site, health risk assessment, restoration scheme design, management and control of the polluted site and safe utilization; providing technical support for identification and evaluation of soil environmental damage of a polluted site; and a scientific basis and a comprehensive management strategy are provided for reutilization, redevelopment and value guarantee of the polluted site.
Disclosure of Invention
The invention aims to provide a method for quantizing a site soil environmental damage identification object based on kriging interpolation, and solves the specific quantization problem of overproof characteristic pollutants in the field of identification and evaluation of contaminated site soil environmental damage.
The invention provides a method for quantizing a site soil environmental damage identification object based on kriging interpolation, which mainly comprises the following implementation steps:
a field soil environmental damage identification object quantification method based on Krigin interpolation comprises the following steps:
1) collecting and analyzing relevant basic information of a field to be quantified, designing a sampling, monitoring and identifying evaluation scheme, and acquiring geographic information coordinates of sampling points and soil samples;
2) performing laboratory analysis and test on the soil sample to obtain the concentration of the characteristic pollutant, and establishing a soil pollution database of the characteristic pollutant of the sampling point;
3) each sampling point location xiCharacteristic pollutant concentration value and the pollutant sieveSelecting values for comparison, wherein if the concentration is less than the threshold effect concentration, the concentration of the point location characteristic pollutant does not exceed the standard; if the concentration is greater than the threshold effect concentration, the concentration of the point location characteristic pollutant exceeds the standard, and the standard exceeding multiple is calculated;
4) respectively sampling each sampling point position x of all fields to be quantizediPerforming Krigin spatial interpolation on the pollutant concentration, estimating the value of an unknown point by using a sample of the known point, predicting the unknown point by using the known point to form an interpolation area graph of the pollutant concentration from low to high, and obtaining the overproof area of the characteristic pollutant;
5) determining the exceeding depth H of the pollutant according to whether the concentration of the characteristic pollutant detected by each layer of the sampled soil exceeds the standard or not1,H2……HnAnd until the concentration of the characteristic pollutant detected in the last layer does not exceed the standard, the final pollutant exceeding depth is specifically calculated as H ═ H1+H2+……+Hn(ii) a And quantifying to obtain the polluted earthwork amount according to an area map and the sampling depth generated by the Kerrin spatial interpolation, and determining the physical amount of the soil environment damage of the polluted site.
In the step 1), the geographic information coordinates of each sampling point location of the field to be quantized are obtained, and each sampling point location is respectively recorded as each corresponding point xi(i is 1, 2, 3, … …, n), and simultaneously, carrying out layered sampling on the soil of each sampling point, and recording the sampling depth from the surface layer to the deep layer as the pollutant exceeding depth H1,H2,……Hn。
In the step 2), each sampling point x of the field to be quantized is obtainediContaminant concentration value Z (x) (i ═ 1, 2, … …, n)i) For the contaminant concentration value Z (x)i) Performing normal distribution test, and performing logarithmic transformation or power exponential transformation and the like on the data which are not in accordance with the normal distribution to make the data in accordance with the normal distribution; then, the characteristic pollutant value of the known sampling point on the research area A is assumed to be Z (x)i) (i ═ 1, 2, … …, n), then the strip estimates point x0The attribute Z (x) at E A0) Kriging interpolation result Z (x)0) Is the known sample point attribute value Z (x)i) (i ═ 1, 2, … …, n) weighted sumI.e. byIn the formula ofiIs the undetermined weight coefficient, representing each spatial sample point xiConcentration value Z (x) ofi) To be estimated value Z x (x)0) The degree of contribution of.
The Z (xi) has a relation of distance and relative direction change, the object to be researched is called as 'regional variable' by the Krigin interpolation method, and the undetermined weight coefficient lambda can be obtained by obtaining the unbiased condition aiming at the unbiased and minimum variance condition of the Krigin methodi(i ═ 1, 2, … …, n) satisfies the relationship:on the premise of unbiased condition, the solved to-be-determined weight coefficient lambda can be obtained by taking the kriging variance as minimumiThe system of equations of (1):
in the formula, C (x)i,xj) Is two points x within the investigation regioni,xjμ is the lagrange multiplier.
Two points x within the investigation regioni,xjIs specifically based onCalculating to obtain a variation function of xi and xj; where C (xi, xj) is a function of variation of two points xi and xj within the study area, Var represents variance, z (xi), and z (xj) is a characteristic contaminant value (i ═ 1, 2, … …, n) (j ═ 1, 2, … …, n).
In the step 3), each sampling point x is calculatediThe formula of the pollutant concentration standard-exceeding multiple is K ═ xiDetermining standard exceeding multiple K of each sampling point position of the characteristic pollutant in each single-layer soil under different sampling depths of a field to be quantified; wherein K is the exceeding multiple of the soil characteristic pollutant, B is the baseline threshold of the concentration of the characteristic pollutant in the soilThe value is obtained.
In the step 4), each sampling point x of all fields to be quantifiediThe specific method for carrying out spatial interpolation on the pollutant concentration comprises the following steps:
4.1) carrying out normal distribution test on the field pollutant concentration data, and carrying out logarithmic transformation or power exponential transformation on the data which are not in accordance with the normal distribution so as to enable the sampling point data to be in accordance with the normal distribution;
4.2) predicting the spatial distribution form of the overproof pollutants by a kriging interpolation model;
4.3) performing Critical interpolation to obtain the standard exceeding pollutant distribution in the form of regular grid, and setting the size of the grid (G) to be 1 x 1m2The pollutant area calculation formula is S ═ C × G, wherein S is the pollutant overproof area, and C is the total number of grids;
4.4) calculating the standard exceeding area S of each single-layer characteristic pollutant under different sampling depths1……Sn;
4.5) dividing each single layer into S with characteristic pollutant concentration not exceeding the standard according to different pollution degrees10Slight contamination with S11Moderate contamination S12Severe contamination S13To refine the overproof area;
4.6) the overproof area of each monolayer of pollutants is S1=S11+S12+S13The total pollutant overproof area is S ═ S1+……+Sn。
In the step 5), the specific method is as follows:
5.1) whether the concentration of the characteristic pollutant detected by each layer of sampled soil exceeds the standard or not, determining the exceeding concentration of the pollutant through experimental analysis, and further determining the exceeding depth H of the pollutant1,H2……HnUntil the concentration of the pollutants detected in the last layer does not exceed the standard; the final contamination depth is specifically calculated as H ═ H1+H2+……+Hn(ii) a The calculation formula of the polluted earthwork amount of the site to be quantified is V-S-H;
5.2) calculating the polluted soil volume V of each single-layer characteristic pollutant under different sampling depths1,V2……Vn;
5.3) according to different pollution degrees, each single layer can be divided into V with characteristic pollutant concentration not exceeding the standard10Slight contamination V11Moderate pollution V12Severe contamination V13The excessive earth volume is refined, and a scientific basis is provided for the management and control and the safety utilization of the polluted site;
5.4) the over-standard earth volume of each monolayer pollutant is V1=V11+V12+V13. The total volume of the polluted earthwork is V ═ V1+V2+……+Vn(ii) a In the formula, V is the earth volume exceeding the pollution standard, and S is the area exceeding the pollutant standard.
The invention has the beneficial effects that: the invention provides a method for quantizing a site soil environmental damage identification object based on a kriging interpolation, which quantizes the pollution depth, the pollution area and the earth volume of site characteristic pollutants, determines the object amount of the site soil environmental damage, and provides scientific basis for value quantization, health risk assessment, restoration scheme design, and control, management and safety utilization of a polluted site; providing technical support for identification and evaluation of soil environmental damage of a polluted site; and a scientific basis and a comprehensive management strategy are provided for reutilization, redevelopment and value guarantee of the polluted site.
Drawings
FIG. 1 is a real object quantification step of identification and evaluation of soil environmental damage of a contaminated site based on Krigin spatial interpolation;
FIG. 2 is a plot of distribution and partitioning of a coal processing plant;
fig. 3 and 4 are graphs of the exceeding area of the characteristic pollutants of different layers of soil.
Detailed Description
The invention is further described with reference to the accompanying drawings and the specific implementation process.
As shown in fig. 1, a real object quantification step for identifying site soil environmental damage is specifically implemented as follows:
1. basic data collection
1.1) collecting needed relevant data through local departments, organizations and other ways, mainly comprising: topography, soil type map, land type map, etc., as shown in FIG. 2;
1.2) acquiring geographic information coordinates of each sampling point of a field to be quantized; recording each sampling point as each corresponding point xi(i=1,2,……,n);
1.3) carrying out layered sampling on the soil at each sampling point, and recording the sampling depth from the surface layer to the deep layer as H1,H2,……Hn。
2. Establishment of data and model base
2.1) acquiring a pollutant concentration value Z (xi) of each sampling point xi (i ═ 1, 2, … …, n) of the site to be quantified; suppose that the characteristic contaminant value of the known sampling point on the study area A is Z (x)i) (i ═ 1, 2, … …, n), then point x is to be estimated0The attribute Z (x) at E A0) Kriging interpolation result Z (x)0) Is a weighted sum of known sample point attribute values z (xi) (i ═ 1, 2, … …, n), i.e.In the formula ofiIs an undetermined weight coefficient which represents the concentration value Z (xi) at each spatial sampling point xi and the value Z x (x) to be estimated0) The degree of contribution of (c);
2.2) the relation between the distance and the change of the relative direction exists between Z (xi), and the object to be researched is called as 'regional variable' by a Krigin interpolation method;
2.3) obtaining the undetermined weight coefficient lambda aiming at the unbiased condition and the unbiased condition obtained by the minimum variance condition of the Kriging methodi(i ═ 1, 2, … …, n) satisfies the relationship:on the premise of unbiased condition, the solved to-be-determined weight coefficient lambda can be obtained by taking the kriging variance as minimumiThe system of equations of (1):in the formula, C (x)i,xj) Is two points x within the investigation regioni,xjμ is a lagrange multiplier;
2.4) theThe variation function of two points xi, xj in the region of investigation is based in particular onCalculating a variation function of xi and xj, wherein C (xi, xj) is the variation function of xi and xj at two points in the research area, Var represents variance, z (xi), and z (xj) is a characteristic pollutant value (i ═ 1, 2, … …, n) (j ═ 1, 2, … …, n).
3. Exceeding multiple of each sampling point of field to be quantized
3.1) sampling points xiComparing the characteristic contaminant concentration value of (a) with the contaminant screening value; if the concentration is less than the threshold effect concentration, the concentration of the point is not out of standard; if the concentration is greater than the threshold effect concentration, the concentration of the pollutant at the point exceeds the standard;
3.2) calculating the respective sample points xiThe formula of the standard exceeding multiple of the characteristic pollutant concentration is K ═ xiDetermining standard exceeding multiple K of each sampling point position of the characteristic pollutant in each single-layer soil under different sampling depths of a field to be quantified; k is the exceeding multiple of the soil characteristic pollutant, and B is the baseline threshold value of the concentration of the characteristic pollutant in the soil.
4. Kriging spatial interpolation
The application steps of the kriging interpolation method are as follows:
4.1) inputting original data, namely sampling points;
4.2) gridding, namely selecting the range of the area and the size of a grid, and carrying out gridding treatment on the area;
4.3) data inspection and analysis, observing whether the data accord with normal distribution or not, and if not, carrying out data transformation; then removing obvious difference points according to whether the sampling value is in accordance with the actual situation;
4.4) carrying out variation function calculation by using the variation function to know the space structure of the variable; 4.5) Critical interpolation estimation.
5. Area exceeding standard of site pollutant to be quantified
Respectively carrying out the quantization on all sampling points x of all fields to be quantizediThe specific method for carrying out spatial interpolation on the pollutant concentration comprises the following steps:
5.1) carrying out normal distribution test on the field pollutant concentration data, and carrying out logarithmic transformation or power exponential transformation on the data which are not in accordance with the normal distribution so as to enable the sampling point data to be in accordance with the normal distribution;
5.2) predicting the spatial distribution form of the overproof pollutants by a kriging interpolation model;
5.3) as shown in FIGS. 3 and 4: performing kriging interpolation to obtain standard-exceeding pollutant distribution in regular grid form, and setting the size of grid (G) to 1 × 1m2Then the area calculation formula of the pollutants is S ═ C × G; s is the pollutant standard area, C is the total grid number
5.4) calculating the standard exceeding area S of each single-layer characteristic pollutant under different sampling depths1……Sn;
5.5) according to different pollution levels, each single layer can be divided into S with characteristic pollutant concentration not exceeding standard10Slight contamination with S11Moderate contamination S12Severe contamination S13To refine the overproof area;
5.6) the overproof area of each monolayer of pollutants is S1=S11+S12+S13(ii) a The total pollutant overproof area is S ═ S1+……+Sn。
6. To-be-quantified site pollutant pollution earthwork amount
6.1) determining the overproof depth H of the pollutant according to the concentration of the pollutant detected in each layer1,H2……Hn(ii) a And determining the pollutant overproof concentration through experimental detection and analysis according to each layer of sampled soil, and further determining the pollutant overproof depth H1,H2……HnUntil the concentration of the pollutants detected in the last layer does not exceed the standard; the final pollutant overproof depth is specifically calculated as H ═ H1+H2+……+Hn(ii) a The calculation formula of the polluted earthwork amount of the site to be quantified is V-S-H;
calculating the polluted soil volume V of each single-layer characteristic pollutant under different sampling depths1,V2……Vn;
6.3) depending on the degree of contamination, each monolayer can be classified as a characteristic stainV for dye concentration not exceeding standard10Slight contamination V11Moderate pollution V12Severe contamination V13The excessive earth volume is refined, and the cost is saved for subsequent field restoration;
6.4) the over-standard earth volume of each monolayer of pollutants is V1=V11+V12+V13(ii) a The total volume of the polluted earthwork is V ═ V1+V2+……+Vn(ii) a In the formula, V is the earth volume with the pollutant exceeding the standard, S is the area with the pollutant exceeding the standard, and H is the sampling depth.
Claims (8)
1. A field soil environmental damage identification object quantification method based on Krigin interpolation is characterized by comprising the following steps:
1) collecting and analyzing relevant basic information of a field to be quantified, designing a sampling, monitoring and identifying evaluation scheme, and acquiring geographic information coordinates of sampling points and soil samples;
2) performing laboratory analysis and test on the soil sample to obtain the concentration of the characteristic pollutant, and establishing a soil pollution database of the characteristic pollutant of the sampling point;
3) each sampling point location xiComparing the concentration value of the characteristic pollutant with the screening value of the pollutant, and if the concentration is less than the threshold effect concentration, the concentration of the characteristic pollutant of the point position does not exceed the standard; if the concentration is greater than the threshold effect concentration, the concentration of the point location characteristic pollutant exceeds the standard, and the standard exceeding multiple is calculated;
4) respectively sampling each sampling point position x of all fields to be quantizediPerforming Krigin spatial interpolation on the pollutant concentration, estimating the value of an unknown point by using a sample of the known point, predicting the unknown point by using the known point to form an interpolation area graph of the pollutant concentration from low to high, and obtaining the overproof area of the characteristic pollutant;
5) determining the exceeding depth H of the pollutant according to whether the concentration of the characteristic pollutant detected by each layer of the sampled soil exceeds the standard or not1,H2……HnAnd until the concentration of the characteristic pollutant detected in the last layer does not exceed the standard, the final pollutant exceeding depth is specifically calculated as H ═ H1+H2+……+Hn(ii) a And quantifying to obtain the polluted earthwork amount according to an area map and the sampling depth generated by the Kerrin spatial interpolation, and determining the physical amount of the soil environment damage of the polluted site.
2. The kriging interpolation-based field soil environmental damage identification real object quantization method according to claim 1, wherein in step 1), geographic information coordinates of each sampling point location of a field to be quantized are obtained, and each sampling point location is respectively recorded as each corresponding point xi(i is 1, 2, 3, … …, n), and simultaneously, carrying out layered sampling on the soil of each sampling point, and recording the sampling depth from the surface layer to the deep layer as the pollutant exceeding depth H1,H2,……Hn。
3. The kriging interpolation-based field soil environmental damage identification object quantization method according to claim 1, wherein in the step 2), each sampling point x of a field to be quantized is obtainediContaminant concentration value Z (x) (i ═ 1, 2, … …, n)i) For the contaminant concentration value Z (x)i) Performing normal distribution test, and performing logarithmic transformation or power exponential transformation and the like on the data which are not in accordance with the normal distribution to make the data in accordance with the normal distribution; then, the characteristic pollutant value of the known sampling point on the research area A is assumed to be Z (x)i) (i ═ 1, 2, … …, n), then the strip estimates point x0The attribute Z (x) at E A0) Kriging interpolation result Z (x)0) Is the known sample point attribute value Z (x)i) (i ═ 1, 2, … …, n), i.e. a weighted sumIn the formula ofiIs the undetermined weight coefficient, representing each spatial sample point xiConcentration value Z (x) ofi) To be estimated value Z x (x)0) The degree of contribution of.
4. The kriging interpolation-based field soil environmental damage identification object quantification method according to claim 3, wherein the field soil environmental damage identification object quantification method is characterized in thatZ (xi) has a relation of distance and relative direction change, the object to be researched is called as 'regional variable' by the Krigin interpolation method, and the undetermined weight coefficient lambda can be obtained by obtaining the unbiased condition aiming at the unbiased and minimum variance condition of the Krigin methodi(i ═ 1, 2, … …, n) satisfies the relationship:on the premise of unbiased condition, the solved to-be-determined weight coefficient lambda can be obtained by taking the kriging variance as minimumiThe system of equations of (1):
in the formula, C (x)i,xj) Is two points x within the investigation regioni,xjμ is the lagrange multiplier.
5. The method for identifying and evaluating the real object of the soil environmental damage of the polluted site based on the Krigin spatial interpolation as claimed in claim 4, wherein x is two points in the research areai,xjIs specifically based onCalculating to obtain a variation function of xi and xj; where C (xi, xj) is a function of variation of two points xi and xj within the study area, Var represents variance, z (xi), and z (xj) is a characteristic contaminant value (i ═ 1, 2, … …, n) (j ═ 1, 2, … …, n).
6. The field soil environmental damage identification object quantification method based on kriging interpolation according to claim 1, wherein in the step 3), each sampling point x is calculatediThe formula of the pollutant concentration standard-exceeding multiple is K ═ xiDetermining standard exceeding multiple K of each sampling point position of the characteristic pollutant in each single-layer soil under different sampling depths of a field to be quantified; wherein K is the soil characteristic pollutionThe standard exceeding multiple of the substance, B is a baseline threshold value of the concentration of the characteristic pollutant in the soil.
7. The field soil environmental damage identification object quantification method based on kriging interpolation according to claim 1,
in the step 4), each sampling point x of all fields to be quantifiediThe specific method for carrying out spatial interpolation on the pollutant concentration comprises the following steps:
4.1) carrying out normal distribution test on the field pollutant concentration data, and carrying out logarithmic transformation or power exponential transformation on the data which are not in accordance with the normal distribution so as to enable the sampling point data to be in accordance with the normal distribution;
4.2) predicting the spatial distribution form of the overproof pollutants by a kriging interpolation model;
4.3) performing Critical interpolation to obtain the standard exceeding pollutant distribution in the form of regular grid, and setting the size of the grid (G) to be 1 x 1m2The pollutant area calculation formula is S ═ C × G, wherein S is the pollutant overproof area, and C is the total number of grids;
4.4) calculating the standard exceeding area S of each single-layer characteristic pollutant under different sampling depths1……Sn;
4.5) dividing each single layer into S with characteristic pollutant concentration not exceeding the standard according to different pollution degrees10Slight contamination with S11Moderate contamination S12Severe contamination S13To refine the overproof area;
4.6) the overproof area of each monolayer of pollutants is S1=S11+S12+S13The total pollutant overproof area is S ═ S1+……+Sn。
8. The field soil environmental damage identification object quantification method based on kriging interpolation according to claim 1,
in the step 5), the specific method is as follows:
5.1) whether the concentration of the characteristic pollutant detected by each layer of sampled soil exceeds the standard or not, determining the exceeding concentration of the pollutant through experimental analysis, and further determining that the pollutant exceeds the standardStandard depth H1,H2……HnUntil the concentration of the pollutants detected in the last layer does not exceed the standard; the final contamination depth is specifically calculated as H ═ H1+H2+……+Hn(ii) a The calculation formula of the polluted earthwork amount of the site to be quantified is V-S-H;
5.2) calculating the polluted soil volume V of each single-layer characteristic pollutant under different sampling depths1,V2……Vn;
5.3) according to different pollution degrees, each single layer can be divided into V with characteristic pollutant concentration not exceeding the standard10Slight contamination V11Moderate pollution V12Severe contamination V13The excessive earth volume is refined, and a scientific basis is provided for the management and control and the safety utilization of the polluted site;
5.4) the over-standard earth volume of each monolayer pollutant is V1=V11+V12+V13(ii) a The total volume of the polluted earthwork is V ═ V1+V2+……+Vn(ii) a In the formula, V is the earth volume exceeding the pollution standard, and S is the area exceeding the pollutant standard.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011408044.XA CN112686497A (en) | 2020-12-05 | 2020-12-05 | Method for quantizing field soil environmental damage identification object based on kriging interpolation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011408044.XA CN112686497A (en) | 2020-12-05 | 2020-12-05 | Method for quantizing field soil environmental damage identification object based on kriging interpolation |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112686497A true CN112686497A (en) | 2021-04-20 |
Family
ID=75447363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011408044.XA Pending CN112686497A (en) | 2020-12-05 | 2020-12-05 | Method for quantizing field soil environmental damage identification object based on kriging interpolation |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112686497A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113313300A (en) * | 2021-05-25 | 2021-08-27 | 辽宁大学 | Method for identifying soil environment damage physical quantity based on reverse distance weighted interpolation |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006349468A (en) * | 2005-06-15 | 2006-12-28 | Ritsumeikan | Estimation system of pollution distribution of soil, estimation method of pollution distribution of soil and estimation program of pollution distribution of soil |
CN107423561A (en) * | 2017-07-11 | 2017-12-01 | 云南瀚哲科技有限公司 | A kind of evaluation method of soil attribute interpolation |
CN110826236A (en) * | 2019-11-11 | 2020-02-21 | 北京建工环境修复股份有限公司 | Method for judging soil remediation pollution range |
CN111552924A (en) * | 2020-04-22 | 2020-08-18 | 中国科学院地理科学与资源研究所 | Method for evaluating heavy metal pollution characteristics and potential ecological risks of soil on scale of villages and towns |
-
2020
- 2020-12-05 CN CN202011408044.XA patent/CN112686497A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006349468A (en) * | 2005-06-15 | 2006-12-28 | Ritsumeikan | Estimation system of pollution distribution of soil, estimation method of pollution distribution of soil and estimation program of pollution distribution of soil |
CN107423561A (en) * | 2017-07-11 | 2017-12-01 | 云南瀚哲科技有限公司 | A kind of evaluation method of soil attribute interpolation |
CN110826236A (en) * | 2019-11-11 | 2020-02-21 | 北京建工环境修复股份有限公司 | Method for judging soil remediation pollution range |
CN111552924A (en) * | 2020-04-22 | 2020-08-18 | 中国科学院地理科学与资源研究所 | Method for evaluating heavy metal pollution characteristics and potential ecological risks of soil on scale of villages and towns |
Non-Patent Citations (1)
Title |
---|
罗娜;陆安祥;王纪华;: "基于空间插值的土壤重金属污染评估分析系统设计与实现", 食品安全质量检测学报, no. 02 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113313300A (en) * | 2021-05-25 | 2021-08-27 | 辽宁大学 | Method for identifying soil environment damage physical quantity based on reverse distance weighted interpolation |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Dickson et al. | Mapping the probability of large fire occurrence in northern Arizona, USA | |
Sorbino et al. | Susceptibility analysis of shallow landslides source areas using physically based models | |
Shad et al. | Predicting air pollution using fuzzy genetic linear membership kriging in GIS | |
CN112288247B (en) | Soil heavy metal risk identification method based on space interaction relationship | |
CN112308290B (en) | Method for investigating and distributing points in polluted site | |
Liu et al. | Uncertainty in positive matrix factorization solutions for PAHs in surface sediments of the Yangtze River Estuary in different seasons | |
CN113159454A (en) | Encryption optimization point distribution method for identifying and evaluating soil environment damage physical quantity | |
CN113720639B (en) | Method for investigating and supplementing sampling of polluted site | |
Yang et al. | An integrated simulation-assessment approach for evaluating health risks of groundwater contamination under multiple uncertainties | |
Nejadkoorki et al. | The design of long-term air quality monitoring networks in urban areas using a spatiotemporal approach | |
CN114330812A (en) | Landslide disaster risk assessment method based on machine learning | |
CN116187543A (en) | Machine learning-based soil heavy metal content prediction method and application thereof | |
CN114544913A (en) | Soil investigation encrypted sampling stationing optimization method | |
Mirzaei et al. | Model calibration and uncertainty analysis of runoff in the Zayanderood River basin using generalized likelihood uncertainty estimation (GLUE) method | |
CN112950047A (en) | Progressive identification method for suspected contaminated site | |
Xuchao et al. | Simulating soil erodibility in southeastern China using a sequential Gaussian algorithm | |
KR102012007B1 (en) | Method for Estimating Soil Pollution Area and Soil Pollution Quantity Considering Regional Standards by Land category | |
Miháliková et al. | Capability of different interpolation models and pedotransfer functions to estimate soil hydraulic properties in Büyükçay Watershed | |
Ren et al. | Analysis of the spatial characteristics of inhalable particulate matter concentrations under the influence of a three-dimensional landscape pattern in Xi'an, China | |
CN112686497A (en) | Method for quantizing field soil environmental damage identification object based on kriging interpolation | |
Linhua | Calculating Environmental Background Value: A Comparative Study of Statistical Versus Spatial Analyses. | |
CN108364090B (en) | Method for predicting benzene series accumulated load on surface of urban road | |
Parysow et al. | Assessing uncertainty of erodibility factor in national cooperative soil surveys: a case study at Fort Hood, Texas | |
KR101469168B1 (en) | Method for verifying representativeness of sample using spatial correlation analysis | |
CN115879594A (en) | Urban settlement population distribution trend prediction method based on geographic detector |
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 |