[發(fā)明專利]圍墾地土壤重金屬污染綜合評(píng)價(jià)圖生成方法

無(wú)效

• 申請(qǐng)?zhí)枺?/td>: 201210124951.0
• 申請(qǐng)日: 2012-04-25
• 公開(kāi)（公告）號(hào): CN102636632A
• 公開(kāi)（公告）日: 2012-08-15
• 發(fā)明（設(shè)計(jì)）人: 申廣榮;徐敬敬;錢(qián)振華;黃秀梅;張慶龍
• 申請(qǐng)（專利權(quán)）人: 上海交通大學(xué)
• 主分類號(hào): G01N33/24
• 分類號(hào): G01N33/24
• 代理公司: 上海科盛知識(shí)產(chǎn)權(quán)代理有限公司 31225
• 代理人: 趙繼明

1.一種圍墾地土壤重金屬污染綜合評(píng)價(jià)圖生成方法，其特征在于，該方法包括以下步驟：

1)獲取待測(cè)圍墾地的采樣點(diǎn)土壤樣本；

2)對(duì)上述采樣點(diǎn)土壤樣本的重金屬元素含量數(shù)值進(jìn)行數(shù)據(jù)預(yù)處理，分別獲得各重金屬元素的擬插值子樣本；

3)對(duì)擬插值子樣本進(jìn)行探索性空間數(shù)據(jù)分析處理，獲得采樣點(diǎn)土壤樣本的各向異性參數(shù)、步長(zhǎng)和步數(shù)；

4)對(duì)采樣點(diǎn)土壤樣本重金屬元素含量數(shù)值進(jìn)行函數(shù)擬合，獲得適用于重金屬元素含量插值的理論半方差函數(shù)模型；

5)利用包括局域離群值的采樣點(diǎn)進(jìn)行空間變異預(yù)測(cè)結(jié)果分布圖的繪制；

6)重金屬污染綜合評(píng)價(jià)圖的生成：首先計(jì)算各重金屬元素的污染指數(shù)，然后計(jì)算多因子綜合污染指數(shù)，根據(jù)分類標(biāo)準(zhǔn)生成重金屬污染綜合評(píng)價(jià)圖。

2.根據(jù)權(quán)利要求1所述的一種圍墾地土壤重金屬污染綜合評(píng)價(jià)圖生成方法，其特征在于，所述的步驟2)中的數(shù)據(jù)預(yù)處理具體為：

2.1)通過(guò)重金屬元素含量數(shù)值的偏峰度計(jì)算及檢驗(yàn)，考察其是否符合標(biāo)準(zhǔn)正態(tài)分布條件H₀，若符合標(biāo)準(zhǔn)正態(tài)分布條件H₀，則進(jìn)行步驟3)；否則，執(zhí)行步驟2.2)；

2.2)計(jì)算樣本均值μ、標(biāo)準(zhǔn)差σ，去除區(qū)間(μ-1.96σ，μ+1.96σ)以外的值后再次執(zhí)行步驟2.1)，若符合標(biāo)準(zhǔn)正態(tài)分布條件H₀則進(jìn)行步驟3)，否則執(zhí)行步驟2.3)；

2.3)對(duì)采樣點(diǎn)土壤樣本的數(shù)據(jù)進(jìn)行對(duì)數(shù)變換，再次執(zhí)行步驟2.2)，重新計(jì)算變換后樣本的μ和σ，若符合標(biāo)準(zhǔn)正態(tài)分布條件H₀則進(jìn)行步驟3)，否則執(zhí)行步驟2.4)；

2.4)反復(fù)執(zhí)行步驟2.3)，若符合標(biāo)準(zhǔn)正態(tài)分布條件H₀則進(jìn)行步驟3)，否則采用Box-Cox變換后再執(zhí)行步驟2.2)，若不符合標(biāo)準(zhǔn)正態(tài)分布條件H₀，則停止計(jì)算，重新返回步驟1)，否則執(zhí)行步驟3)。

3.根據(jù)權(quán)利要求2所述的一種圍墾地土壤重金屬污染綜合評(píng)價(jià)圖生成方法，其特征在于，所述的符合標(biāo)準(zhǔn)正態(tài)分布條件H₀是指：重金屬元素含量數(shù)值的三階矩g₁和四階矩g₂滿足如下式的條件時(shí)，表明數(shù)據(jù)為標(biāo)準(zhǔn)正態(tài)分布：

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_α為在顯著性水平為α?xí)r的u值，x_i(i＝1，2，...，n)為第i個(gè)樣本的重金屬元素含量，n為樣本個(gè)數(shù)，為n個(gè)采樣點(diǎn)重金屬元素含量的平均值。