[发明专利]围垦地土壤重金属污染综合评价图生成方法

权利要求书

1.一种围垦地土壤重金属污染综合评价图生成方法，其特征在于，该方法包括以下步骤：

1)获取待测围垦地的采样点土壤样本；

2)对上述采样点土壤样本的重金属元素含量数值进行数据预处理，分别获得各重金属元素的拟插值子样本；

3)对拟插值子样本进行探索性空间数据分析处理，获得采样点土壤样本的各向异性参数、步长和步数；

4)对采样点土壤样本重金属元素含量数值进行函数拟合，获得适用于重金属元素含量插值的理论半方差函数模型；

5)利用包括局域离群值的采样点进行空间变异预测结果分布图的绘制；

6)重金属污染综合评价图的生成：首先计算各重金属元素的污染指数，然后计算多因子综合污染指数，根据分类标准生成重金属污染综合评价图。

2.根据权利要求1所述的一种围垦地土壤重金属污染综合评价图生成方法，其特征在于，所述的步骤2)中的数据预处理具体为：

2.1)通过重金属元素含量数值的偏峰度计算及检验，考察其是否符合标准正态分布条件H₀，若符合标准正态分布条件H₀，则进行步骤3)；否则，执行步骤2.2)；

2.2)计算样本均值μ、标准差σ，去除区间(μ-1.96σ，μ+1.96σ)以外的值后再次执行步骤2.1)，若符合标准正态分布条件H₀则进行步骤3)，否则执行步骤2.3)；

2.3)对采样点土壤样本的数据进行对数变换，再次执行步骤2.2)，重新计算变换后样本的μ和σ，若符合标准正态分布条件H₀则进行步骤3)，否则执行步骤2.4)；

2.4)反复执行步骤2.3)，若符合标准正态分布条件H₀则进行步骤3)，否则采用Box-Cox变换后再执行步骤2.2)，若不符合标准正态分布条件H₀，则停止计算，重新返回步骤1)，否则执行步骤3)。

3.根据权利要求2所述的一种围垦地土壤重金属污染综合评价图生成方法，其特征在于，所述的符合标准正态分布条件H₀是指：重金属元素含量数值的三阶矩g₁和四阶矩g₂满足如下式的条件时，表明数据为标准正态分布：

g1=nΣi=1n(xi-x‾)3(n-1)(n-2)(Σi=1n(xi-x‾)2n-1)32≤uα6n(n-1)(n+1)(n-2)(n+3)]]>

g2=n(n+1)Σi=1n(xi-x‾)4(n-1)(n-2)(n-3)(Σi=1n(xi-x‾)2n-1)2-3(n-1)2(n-2)(n-3)≤uα24n(n-1)(n-3)(n-2)(n+3)(n+5)]]>

式中u_α为在显著性水平为α时的u值，x_i(i＝1，2，...，n)为第i个样本的重金属元素含量，n为样本个数，为n个采样点重金属元素含量的平均值。