CN105678454A - Method for determining soil restoration amount of polluted field on the basis of condition simulation - Google Patents

Abstract

The invention provides a method for determining soil restoration amount of a polluted field on the basis of condition simulation. The method comprises following steps: obtaining soil data of a study area and collecting samples to obtain sample data; performing normal state conversion to the sample data and establishing a semi-variation function model; performing condition simulation to obtain the restoration scope of risk loss; determining the area of polluted area according to the restoration scope and calculating the restoration earth volume in combination with the thickness of soil layers. By use of the method, soil condition of the study area can be comprehended; the scope and earth volume of the study area to be restored are estimated on the basis of condition simulation; the risk brought by soil restoration estimation result uncertainty is quantitatively estimated; the relation between decision results quantized by a transfer function and the risk loss is introduced to determine the restoration scope under the principle of risk loss minimization so that the problem of uncertainty existing in determining polluted scope by spatial interpolation method based on discrete sample points is solved.

Description

A kind of contaminated site soil remediation method for determination of amount based on condition simulation

Technical field

The present invention relates to ambient soil analysis technical field, particularly relate to a kind of contaminated site soil remediation method for determination of amount based on condition simulation.

Background technology

In recent years; iterative method along with " move back two and enter three " policy that Industry Structure adjusts; in many cities, industrial undertaking in succession closes down or transfers; leave over substantial amounts of industrial pollution place; and the soil in these places is also affected by pollution in various degree, its improvement reparation problem is also more and more prominent. China provides foundation in the determination that " contaminated site risk assessment technology directive/guide " (HJ25.3-2014) that 2014 promulgate is contaminated site soil risk controlling value, but to how still to lack science according to determined reparation target delimitation reparation border and estimation earth volume to be repaired in specific construction practice. Repairing border and delimit inaccurate, the fixed too low or too high meeting of earth volume to be repaired causes that repair faces environmental risk or economic risk.

Summary of the invention

The technical problem to be solved in the present invention is to provide a kind of contaminated site soil remediation method for determination of amount based on condition simulation, it is possible to the repair coverage delimited fast and accurately, it is determined that soil remediation amount.

For solving above-mentioned technical problem, embodiments of the invention provide a kind of contaminated site soil remediation method for determination of amount based on condition simulation, and the described contaminated site soil remediation method for determination of amount based on condition simulation includes:

Obtain survey region soil information data and carry out sample collecting acquisition sample data;

Sample data is carried out normal state change, and sets up semivariable function model;

Carry out condition simulation, show that the repair coverage of risk of loss delimited;

Repair coverage according to delimiting determines Polluted area area, in conjunction with the thickness of place layer and then calculate reparation earth volume.

Preferably, described acquisition survey region soil information data include: at least one in regional soil type, infiltration coefficient, natural moisture content, the content of organic matter, polluter and pollution distribution.

Preferably, described acquisition survey region soil information data also carry out sample collecting acquisition sample data, including:

By lattice method, survey region is carried out scene lay sampled point and sample;

Dissimilar ground granule is carried out stratified sampling by binding regional soil data.

Preferably, described sample data is carried out normal state change, including:

Adopting Box-Cox changing method, sample data is carried out normal state change, makes data change into normal distribution, its formula is as follows:

$Y=\left\{\begin{array}{cc}\frac{X^{\mathrm{\λ}}-1}{\mathrm{\λ}},& \mathrm{\λ}\≠0\\ \ln X,& \mathrm{\λ}=0\end{array}\right.$

Wherein, X is initial data; Y is the data after conversion; λ is transformation parameter.

Preferably, described set up semivariable function model, including:

Following formula is utilized to set up semivariable function model:

$\mathrm{\γ}\left(h\right)=\frac{1}{2N\left(h\right)}\underset{i=1}{\overset{N\left(h\right)}{\Σ}}{\[Z\left(x_i\right)-Z(x_i+h)\]}^2$

Wherein, h is the space length of two sample rooms, and N (h) is the sample logarithm that space length is during h, and Z (xi) and Z (xi+h) respectively regionalized variable are at spatial point xi and the xi+h sample value located.

Preferably, described in carry out condition simulation, draw risk of loss repair coverage delimit, including:

Distribution characteristics according to regionalized variable and semivariable function, order Gauss analogy method carries out condition simulation, obtains contaminated site pollutant and exceedes the probability distribution of certain limit.

Preferably, based on transfer function method, x place, position is labelled unjustifiably for not needing underestimating risk of loss R1 and labelling unjustifiably as needing the risk of loss R2 that over-evaluates caused by restoring area to be characterized as respectively caused by restoring area:

$R_1\left(C\right(x\left)\right)=\left\{\begin{array}{cc}0,& C\left(x\right)\≤C_L\\ {\mathrm{\ω}}_1\left(x\right)\left(C\right(x)-C_L),& C\left(x\right)>C_L\end{array}\right.$

$R_2\left(C\right(x\left)\right)=\left\{\begin{array}{cc}{\mathrm{\&omega;}}_2,& C\left(x\right)<C_L\\ 0,& C\left(x\right)\&GreaterEqual;C_L\end{array}\right.$

Wherein, the concentration value (mg/kg) that C (x) obtains at x place, position for pollutant after condition simulation, C_LFor the reparation target (mg/kg) of pollutant, ω₁For underestimating the health risk loss in the unit grids caused by pollutant levels, ω₂For over-evaluating the unit grids caused by pollutant levels excessively repairs risk of loss.

Preferably, after being repeated condition simulation, what obtain x place, position underestimates risk of loss expected valueWith over-evaluate risk of loss expected valueIt is respectively as follows:

Wherein CⁿX the n-th concentration of analog value (mg/kg) that () obtains at x place, position for pollutant after condition simulation.

Preferably, described in draw risk of loss repair coverage delimit, including:

Should follow, when dividing restoring area, the principle that risk of loss minimizes: if the risk of loss of underestimating at x place, position is expected desirably less than over-evaluating risk of loss, namelyThen illustrate this place be divided into need not restoring area more reasonable; Whereas ifThen this place is divided into needs restoring area more reasonable.

Preferably, described in carry out sample collecting obtain sample data, including:

With reference to the method for testing of regulation in USEPA3060A, USEPA6010C and HJ/T350, sample is analyzed.

Having the beneficial effect that of the technique scheme of the present invention:

In such scheme, by obtaining survey region soil information data and carrying out sample collecting acquisition sample data, survey region soil regime can be fully understanded, based on condition simulation method, study area scope to be repaired and earth volume are estimated, the risk size that the uncertainty of quantitative assessment soil remediation amount estimation result is brought, and introduce the relation between transmission function quantity decision result and risk of loss, the repair coverage being minimised as principle with risk of loss delimited, thus effectively solving the uncertain problem existed based on the space interpolation of discrete sample point when determining pollution range, and then provide technical support for soil environment risk Environmental capacity with control.

Accompanying drawing explanation

Fig. 1 is the contaminated site soil remediation method for determination of amount flow chart based on condition simulation of the present invention.

Detailed description of the invention

For making the technical problem to be solved in the present invention, technical scheme and advantage clearly, it is described in detail below in conjunction with the accompanying drawings and the specific embodiments.

As it is shown in figure 1, a kind of contaminated site soil remediation method for determination of amount based on condition simulation of the embodiment of the present invention, the described contaminated site soil remediation method for determination of amount based on condition simulation includes:

Step 101: obtain survey region soil information data and carry out sample collecting acquisition sample data.

Wherein, described acquisition survey region soil information data may include that at least one in regional soil type, infiltration coefficient, natural moisture content, the content of organic matter, polluter and pollution distribution. Preferably, described acquisition survey region soil information data also carry out sample collecting acquisition sample data, including:

By lattice method, survey region is carried out scene lay sampled point and sample;

Dissimilar ground granule is carried out stratified sampling by binding regional soil data.

Carry out on-the-spot laying sampled point and sampling wherein it is possible to survey region is layouted by professional judgement in conjunction with lattice method. Dissimilar ground granule is carried out stratified sampling by the data collected in conjunction with early stage in the procedures of establishment of point. The sample analysis detection method that layering obtains is with reference to USEPA3060A (USEPA; 1996), USEPA6010C (USEPA; 2007) and HJ/T350 (Environmental Protection Department; 2007) in, the method for testing of regulation is analyzed. when sample and analysis, strictly each link is carried out quality control and guarantee. Determine that the pollutant of concentration over-standard in sample are as object of study (wherein the conventional statistical analysis of soil sample data and conversion use SPSS20.0 software) by analyzing.

In the present embodiment, survey region is carried out early stage environmental survey, it is possible to understand regional environment management and Land_use change risk control requirement.

Step 102: sample data is carried out normal state change, and sets up semivariable function model.

Wherein, described sample data is carried out normal state change, it is possible to including:

Adopting Box-Cox changing method, sample data is carried out normal state change, makes data change into normal distribution, its formula is as follows:

$Y=\left\{\begin{array}{cc}\frac{X^{\mathrm{\&lambda;}}-1}{\mathrm{\&lambda;}},& \mathrm{\&lambda;}\&NotEqual;0\\ \ln X,& \mathrm{\&lambda;}=0\end{array}\right.$

Wherein, X is initial data; Y is the data after conversion; λ is transformation parameter.

Preferably, described semivariable function model is set up, it is possible to including:

Following formula is utilized to set up semivariable function model:

$\mathrm{\&gamma;}\left(h\right)=\frac{1}{2N\left(h\right)}\underset{i=1}{\overset{N\left(h\right)}{\&Sigma;}}{\&lsqb;Z\left(x_i\right)-Z(x_i+h)\&rsqb;}^2$

Wherein, h is the space length of two sample rooms, and N (h) is the sample logarithm that space length is during h, and Z (xi) and Z (xi+h) respectively regionalized variable are at spatial point xi and the xi+h sample value located.

Step 103: carry out condition simulation, show that the repair coverage of risk of loss delimited.

Wherein, described in carry out condition simulation, draw risk of loss repair coverage delimit, it is possible to including:

Distribution characteristics according to regionalized variable and semivariable function, order Gauss analogy method carries out condition simulation, obtains contaminated site pollutant and exceedes the probability distribution of certain limit.

Step 104: the repair coverage according to delimiting determines Polluted area area, in conjunction with the thickness of place layer and then calculate reparation earth volume.

The contaminated site soil remediation method for determination of amount based on condition simulation of the embodiment of the present invention, by obtaining survey region soil information data and carrying out sample collecting acquisition sample data, survey region soil regime can be fully understanded, based on condition simulation method, study area scope to be repaired and earth volume are estimated, the risk size that the uncertainty of quantitative assessment soil remediation amount estimation result is brought, and introduce the relation between transmission function quantity decision result and risk of loss, the repair coverage being minimised as principle with risk of loss delimited, thus effectively solving the uncertain problem existed based on the space interpolation of discrete sample point when determining pollution range, and then provide technical support for soil environment risk Environmental capacity with control.

Preferably, condition simulation is the distribution characteristics according to regionalized variable and semivariable function, a kind of Method of Stochastic of a certain normal distribution of multiple obediences " realization " is produced, it is possible to value is carried out the uncertainty assessment based on theory of probability according to Monte-Carlo method. Algorithm employing order Gauss simulation (SGS) method of condition simulation.

Utilize condition simulation method can obtain a certain pollutant of contaminated site and exceed the probability distribution of certain limit. For the relation between quantity decision result and risk of loss, by transfer function method, x place, position is labelled unjustifiably for need not underestimating risk of loss R1 and labelling unjustifiably as needing the risk of loss R2 that over-evaluates caused by restoring area to be characterized as respectively caused by restoring area:

$R_1\left(C\right(x\left)\right)=\left\{\begin{array}{cc}0,& C\left(x\right)\&le;C_L\\ {\mathrm{\&omega;}}_1\left(x\right)\left(C\right(x)-C_L),& C\left(x\right)>C_L\end{array}\right.$

$R_2\left(C\right(x\left)\right)=\left\{\begin{array}{cc}{\mathrm{\&omega;}}_2,& C\left(x\right)<C_L\\ 0,& C\left(x\right)\&GreaterEqual;C_L\end{array}\right.$

Wherein, the concentration value (mg/kg) that C (x) obtains at x place, position for pollutant after condition simulation, C_LFor the reparation target (mg/kg) of pollutant, ω₁For underestimating the health risk loss in the unit grids caused by pollutant levels, ω₂For over-evaluating the unit grids caused by pollutant levels excessively repairs risk of loss.

Preferably, ω 1 is for underestimating the health risk loss in the unit grids caused by pollutant levels, and its loss is relevant with the concentration of pollutant, and unit is unit/(mg/kg); ω 2 excessively repairs risk of loss for over-evaluating in the unit grids caused by pollutant levels, and its loss is unrelated with the concentration of pollutant, is definite value on certain grid area, and unit is unit.

Exposure people caused due to different pollutant is different, and therefore ω 1 should be tried to achieve by the health risk infringement model that comprises pollutant toxicity parameter; ω 2 also should carry out estimation according to the intended recovery scenario in place and obtain. The method application in contaminated site repair coverage and earth volume to be repaired are determined, the risk of loss of unit grids is decided to be respectively with reference to existing document: ω 1=10 unit/(g/kg); ω 2=100 unit.

Preferably, after being repeated condition simulation, what obtain x place, position underestimates risk of loss expected valueWith over-evaluate risk of loss expected valueIt is respectively as follows:

Wherein CⁿX the n-th concentration of analog value (mg/kg) that () obtains at x place, position for pollutant after condition simulation.

Preferably, described in draw risk of loss repair coverage delimit, including:

Should follow, when dividing restoring area, the principle that risk of loss minimizes: if the risk of loss of underestimating at x place, position is expected desirably less than over-evaluating risk of loss, namelyThen illustrate this place be divided into need not restoring area more reasonable; Whereas ifThen this place is divided into needs restoring area more reasonable.

Preferably, described in carry out sample collecting obtain sample data, including:

With reference to the method for testing of regulation in USEPA3060A, USEPA6010C and HJ/T350, sample is analyzed.

Conventional statistical analysis and the conversion of the soil sample data of this research use SPSS20.0 software; Semivariable function matching, Ordinary Kriging Interpolation interpolation and condition simulation adopt GS+9.0 software; Other definitiveness interpolation models and repair coverage drawing, raster symbol-base use ArcGIS9.3 software.

The soil remediation scope of different soil and earth volume to be repaired all may utilize same method and carry out delimiting and estimating. Assuming that pollutant have identical spatial distribution in same soil layer, the thickness being multiplied by place layer by the area to be repaired of certain aspect calculates this layer of earth volume to be repaired. Analogy method all carries out on the grid of 8m × 8m, considers simulation precision and the two aspect factors that assess the cost, the number realization of condition simulation is decided to be 200 times.

Ignore pollutant vertical characteristics variability in same soil layer, according to place soil thickness it is estimated that the earth volume to be repaired that obtains of different interpolation method, and utilize the cumulative probability distribution situation of soil earth volume to be repaired that condition simulation method obtains. Cumulative probability can be considered as when determining different earth volume to be repaired faced by Relative risk value. earth volume to be repaired is chosen to more big, the risk that the uncertainty being distributed by soil pollutant is brought is more little. and in condition simulation technology actual application, should treat, according to economic rationality and acceptable risk size, the determination repairing earth volume and consider.

The result obtained by n times condition simulation, by transmitting, obtains N number of because drawing underestimating risk of loss and over-evaluating risk of loss distribution caused by repair coverage by mistake; Utilize the raster symbol-base function of ArcGIS, and then obtain study area underestimating risk of loss expectation and over-evaluating risk of loss desired distribution figure based on n times condition simulation. as seen from the figure, having the region underestimating more greatly risk of loss, generally to over-evaluate risk of loss less, and vice versa. to the Spring layer of key pollutants content, be there is similar Spatial Distribution Pattern by underestimating the excessive risk loss region that soil pollution substrate concentration causes.

According to risk of loss minimization principle, to underestimate risk of loss less than the topsoil scope to be repaired over-evaluated risk of loss and delimit for standard. need the region repaired to be predominantly located at and exceed the sample point repairing target and compare the region of concentration; Have in the middle part of the place of the sample point that sporadicly exceeds standard for distribution, risk of loss is underestimated despite certain, but over-evaluate risk of loss less than this region, therefore it being divided into and do not need restoring area. the repair coverage demarcation method that this research adopts is based on actual risk of loss, and result has certain practical significance.

The above is the preferred embodiment of the present invention; it should be pointed out that, for those skilled in the art, under the premise without departing from principle of the present invention; can also making some improvements and modifications, these improvements and modifications also should be regarded as protection scope of the present invention.

Claims (10)

1. the contaminated site soil remediation method for determination of amount based on condition simulation, it is characterised in that the described contaminated site soil remediation method for determination of amount based on condition simulation includes:

Obtain survey region soil information data and carry out sample collecting acquisition sample data;

Sample data is carried out normal state change, and sets up semivariable function model;

Carry out condition simulation, show that the repair coverage of risk of loss delimited;

Repair coverage according to delimiting determines Polluted area area, in conjunction with the thickness of place layer and then calculate reparation earth volume.

2. the contaminated site soil remediation method for determination of amount based on condition simulation according to claim 1, it is characterized in that, described acquisition survey region soil information data include: at least one in regional soil type, infiltration coefficient, natural moisture content, the content of organic matter, polluter and pollution distribution.

3. the contaminated site soil remediation method for determination of amount based on condition simulation according to claim 1, it is characterised in that described acquisition survey region soil information data also carry out sample collecting acquisition sample data, including:

By lattice method, survey region is carried out scene lay sampled point and sample;

Dissimilar ground granule is carried out stratified sampling by binding regional soil data.

4. the contaminated site soil remediation method for determination of amount based on condition simulation according to claim 1, it is characterised in that described sample data is carried out normal state change, including:

Adopting Box-Cox changing method, sample data is carried out normal state change, makes data change into normal distribution, its formula is as follows:

$Y=\left\{\begin{array}{cc}\frac{X^{\mathrm{\&lambda;}}-1}{\mathrm{\&lambda;}},& \mathrm{\&lambda;}\&NotEqual;0\\ \ln X,& \mathrm{\&lambda;}=0\end{array}\right.$

Wherein, X is initial data;Y is the data after conversion; λ is transformation parameter.

5. the contaminated site soil remediation method for determination of amount based on condition simulation according to claim 4, it is characterised in that described set up semivariable function model, including:

Following formula is utilized to set up semivariable function model:

$\mathrm{\&gamma;}\left(h\right)=\frac{1}{2N\left(h\right)}\underset{i=1}{\overset{N\left(h\right)}{\&Sigma;}}{\&lsqb;Z\left(x_i\right)-Z(x_i+h)\&rsqb;}^2$

Wherein, h is the space length of two sample rooms, and N (h) is the sample logarithm that space length is during h, and Z (xi) and Z (xi+h) respectively regionalized variable are at spatial point xi and the xi+h sample value located.

6. the contaminated site soil remediation method for determination of amount based on condition simulation according to claim 5, it is characterised in that described in carry out condition simulation, draw risk of loss repair coverage delimit, including:

Distribution characteristics according to regionalized variable and semivariable function, order Gauss analogy method carries out condition simulation, obtains contaminated site pollutant and exceedes the probability distribution of certain limit.

7. the contaminated site soil remediation method for determination of amount based on condition simulation according to claim 6, it is characterized in that, based on transfer function method, x place, position is labelled unjustifiably for not needing underestimating risk of loss R1 and labelling unjustifiably as needing the risk of loss R2 that over-evaluates caused by restoring area to be characterized as respectively caused by restoring area:

$R_1\left(C\right(x\left)\right)=\left\{\begin{array}{cc}0,& C\left(x\right)\&le;C_L\\ {\mathrm{\&omega;}}_1\left(x\right)\left(C\right(x)-C_L),& C\left(x\right)>C_L\end{array}\right.$

$R_2\left(C\right(x\left)\right)=\left\{\begin{array}{cc}{\mathrm{\&omega;}}_2,& C\left(x\right)<C_L\\ 0,& C\left(x\right)\&GreaterEqual;C_L\end{array}\right.$

Wherein, the concentration value (mg/kg) that C (x) obtains at x place, position for pollutant after condition simulation, C_LFor the reparation target (mg/kg) of pollutant, ω₁For underestimating the health risk loss in the unit grids caused by pollutant levels, ω₂For over-evaluating the unit grids caused by pollutant levels excessively repairs risk of loss.

8. the contaminated site soil remediation method for determination of amount based on condition simulation according to claim 7, it is characterised in that after being repeated condition simulation, what obtain x place, position underestimates risk of loss expected valueWith over-evaluate risk of loss expected valueIt is respectively as follows:

Wherein CⁿX the n-th concentration of analog value (mg/kg) that () obtains at x place, position for pollutant after condition simulation.

9. the contaminated site soil remediation method for determination of amount based on condition simulation according to claim 8, it is characterised in that described in draw risk of loss repair coverage delimit, including:

Should follow, when dividing restoring area, the principle that risk of loss minimizes: if the risk of loss of underestimating at x place, position is expected desirably less than over-evaluating risk of loss, namelyThen illustrate this place be divided into need not restoring area more reasonable; Whereas ifThen this place is divided into needs restoring area more reasonable.

10. the contaminated site soil remediation method for determination of amount based on condition simulation according to claim 3, it is characterised in that described in carry out sample collecting obtain sample data, including:

With reference to the method for testing of regulation in USEPA3060A, USEPA6010C and HJ/T350, sample is analyzed.