張 慧，馬鑫鵬，蘇 航，隋虹均，史曉磊，李昕陽，楊 歡，鄭志志

(東北農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院，哈爾濱 150030)

土壤背景值通常是指在未受或者很少受到人類活動影響的條件下土壤本身的物質(zhì)組成及含量，是反映土壤屬性的重要參考依據(jù)[1]。土壤重金屬作為土壤環(huán)境的組成部分，其含量狀況能夠直接影響農(nóng)作物的產(chǎn)量和質(zhì)量，并且能夠通過食物鏈進(jìn)入人體，影響人體健康[2]。因此合理確定土壤背景值并以此評估土壤中重金屬的含量狀況具有重要的現(xiàn)實意義[3]。

當(dāng)前對土壤重金屬含量狀況或污染評價研究采用的評價標(biāo)準(zhǔn)包括國家土壤環(huán)境質(zhì)量一級標(biāo)準(zhǔn)[4]、省市土壤背景值[5]、區(qū)域土壤背景值[6]、深層土壤重金屬含量均值[7]等等，大部分研究選取土壤背景值作為評價標(biāo)準(zhǔn)。不同的學(xué)者從各自的角度對土壤背景值進(jìn)行了研究，主要集中在土壤重金屬背景值空間分布狀況，不同城市、地質(zhì)條件、地貌類型等條件下的背景值含量狀況，土壤重金屬背景值的確定方法等方面[8-12]，但是對于土壤重金屬背景值進(jìn)行系統(tǒng)、全面的研究還相對較少，并且對背景值確定方法的研究通常是以一種重金屬作為研究對象，研究結(jié)果的可靠性以及對不同重金屬的適用性有待進(jìn)一步考究。

本文以松嫩平原位于黑龍江省南部的區(qū)域作為研究區(qū)，以多目標(biāo)區(qū)域地球化學(xué)調(diào)查成果作為數(shù)據(jù)基礎(chǔ)，采用置信區(qū)間法、累積頻率法、迭代剔除法[13]對生物毒性明顯的5種重金屬(As、Hg、Cd、Cr、Pb)以及具有一定毒性的3種重金屬(Cu、Ni、Zn)的背景值進(jìn)行研究，確定8種重金屬的背景值。在此基礎(chǔ)上采用地統(tǒng)計學(xué)、普通克里金插值法等對土壤重金屬的空間分布進(jìn)行探究[14]，并根據(jù)確定的重金屬背景值，采用地累積指數(shù)法對研究區(qū)土壤重金屬累積、污染程度進(jìn)行評價，以期為土壤重金屬背景值研究提供參考，為松嫩平原土壤重金屬累積特征、污染評價等提供理論指導(dǎo)，同時也為松嫩平原農(nóng)業(yè)生產(chǎn)、糧食安全、土壤保護(hù)等提供實踐基礎(chǔ)。

1 研究區(qū)概況

研究區(qū)主要是松嫩平原位于黑龍江省南部的區(qū)域，地理坐標(biāo)為122°20′E～128°00′E，44°40′N～48°00′N，包括齊齊哈爾市南部、大慶市、綏化市南部、哈爾濱市西部等28個市(縣)，覆蓋區(qū)域面積約為8.15萬km2，占黑龍江省總面積的17.23%(見圖1)。地貌景觀以平原為主，地勢平坦開闊，嫩江、松花江、拉林河、呼蘭河等水系及其支流穿過。地貌類型主要包括褶斷剝蝕低山、侵蝕-冰水臺地、沖積湖積低平原、剝蝕沖積湖積平原、沖積平原、沖積-洪積平原、風(fēng)蝕風(fēng)積沙地等。研究區(qū)地處中溫帶，夏季溫?zé)岫嘤辏竞渎L。土壤以黑土、黑鈣土、草甸土為主，土壤肥力較高，主要以農(nóng)業(yè)為主，是我國重要的商品糧生產(chǎn)基地。

2 材料與方法

2.1 數(shù)據(jù)來源及處理

研究數(shù)據(jù)以多目標(biāo)區(qū)域地球化學(xué)調(diào)查成果為基礎(chǔ)，樣本采集按照多目標(biāo)區(qū)域地球化學(xué)調(diào)查技術(shù)規(guī)范要求進(jìn)行。采用雙層網(wǎng)格化采樣方法，使用GPS并結(jié)合地形圖定點采集表層和深層土壤樣品，并記錄采樣點的經(jīng)緯度、海拔高度等信息，在研究區(qū)表層(0～20 cm)按照2 km×2 km的格網(wǎng)采集四個土樣混合成一個土壤樣本，共采集20923個土壤樣本；在研究區(qū)深層(180～200 cm)按照4 km×4 km的格網(wǎng)采集四個土樣混合成一個土壤樣本，共采集5295個土壤樣本。

采集的土壤樣本進(jìn)行自然風(fēng)干，剔除植物根系、砂礫、動物殘體等雜質(zhì)，經(jīng)研磨后過20目尼龍篩，然后對土樣進(jìn)行消解，重金屬As采用硝酸-鹽酸-高氯酸消解，重金屬Hg采用硝酸-硫酸-五氧化二釩消解，重金屬Cr采用硫酸-硝酸-氫氟酸消解，重金屬Pb采用鹽酸-硝酸-氫氟酸-高氯酸消解，重金屬Cd、Cu、Ni、Zn采用鹽酸-硝酸-高氯酸消解。重金屬土樣進(jìn)行消解后，在實驗室測定元素含量。

圖1 研究區(qū)位置示意圖Fig.1 Map of study area

Cr、Cu、Ni、Pb、Zn元素含量采用原子吸收光譜法測定；Cd元素含量采用石墨爐原子吸收光譜法測定；Hg元素含量采用冷原子吸收光譜法測定；As元素含量采用硼氫化鉀-硝酸銀分光光度法測定，同時采用插入國家標(biāo)準(zhǔn)物質(zhì)和加標(biāo)回收兩種方式來保證測試的準(zhǔn)確度。

2.2 研究方法

2.2.1 土壤重金屬背景值確定方法

(1)置信區(qū)間法。置信區(qū)間法在土壤背景值確定的應(yīng)用中通常是按照數(shù)據(jù)正態(tài)分布類型的不同，采用符合數(shù)據(jù)特征的計算方法確定區(qū)域土壤元素背景值的上限。此方法對元素含量數(shù)據(jù)的要求較高，數(shù)據(jù)需要服從正態(tài)分布或者對數(shù)正態(tài)分布，一般情況下采用95%置信上限來代表區(qū)域土壤元素背景值。

(2)累積頻率法。累積頻率法是在相對累積頻率和元素含量分布圖中，曲線拐點所對應(yīng)的元素含量通常是該元素背景含量和異常值的分界點，高于拐點的樣點代表了異常值，低于拐點對應(yīng)含量的樣本才能反映出區(qū)域土壤重金屬背景值狀況，元素含量低于拐點的樣本含量均值可代表該元素的土壤背景值，同時該方法對數(shù)據(jù)也有要求，低于拐點的樣本含量數(shù)據(jù)應(yīng)該屬于正態(tài)分布或者近似正態(tài)分布。

(3)迭代剔除法。迭代剔除法基于數(shù)據(jù)均值和標(biāo)準(zhǔn)差將異常值剔除，對剩余的數(shù)據(jù)進(jìn)行計算以求取土壤重金屬背景值[15]。

2.2.2 土壤重金屬空間分布研究方法 采用普通克里金插值法。普通克里金插值法從變量自身特點出發(fā)，考慮觀測點的整體空間分布情況，考慮樣本點的空間相互位置關(guān)系與待估計點的空間位置關(guān)系，對待估點進(jìn)行一種無偏最優(yōu)估計，并且能給出估計精度[16]。在ArcGIS 10.2軟件的支撐下，運(yùn)用普通克里金插值法對樣點進(jìn)行空間插值，進(jìn)而研究重金屬在空間上的分布特征。

2.2.3 土壤重金屬污染評價方法 采用地累積指數(shù)法。地累積指數(shù)法通常稱為Muller指數(shù)，由德國海德堡大學(xué)沉積物研究所的科學(xué)家Muller在1969年提出，最初廣泛應(yīng)用于定量研究水環(huán)境沉積物中重金屬污染程度，近年來逐漸引入到土壤重金屬污染評價中[17]，其表達(dá)式為：

Igeo=log2[Cn(k×Bn)]

其中：Igeo為地累積指數(shù)，Cn為樣品中元素n含量的實測值(mg·kg-1)；Bn為元素n的化學(xué)背景值(mg·kg-1)；k為修正系數(shù)，一般取1.5。

地累積指數(shù)一般分為7個等級(0～6級)，污染程度從無污染到極強(qiáng)污染，具體等級劃分見表1。

3 結(jié)果與分析

3.1 土壤重金屬背景值

研究區(qū)表、深層8種土壤重金屬元素含量描述性統(tǒng)計分析結(jié)果見表2和表3。采用置信區(qū)間法、累積頻率法、迭代剔除法三種方法，應(yīng)用SPSS 19.0軟件對深層土壤(180～200 cm)重金屬樣本數(shù)據(jù)進(jìn)行統(tǒng)計，結(jié)果見表4。從表中可以看出，采用置信區(qū)間法和累積頻率法計算出的土壤重金屬背景值基本上相一致，而采用迭代剔除法計算出的背景值含量相對于前兩種方法偏差較大，主要是由于在進(jìn)行迭代剔除時去除了個別元素大量的樣本數(shù)據(jù)，致使個別元素樣本數(shù)據(jù)量下降，在一定程度上不能很好地反映區(qū)域背景值整體情況，而置信區(qū)間法采用95%的置信區(qū)間，累積頻率法通過分布曲線將自然背景值分離出來，兩種方法測定的背景值相一致，因此采用置信區(qū)間法和累積頻率法計算求得的土壤重金屬元素背景值求均值的結(jié)果作為松嫩平原黑龍江省南部土壤重金屬背景值(見表5)。

與國家《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB15618-1995)一級標(biāo)準(zhǔn)相比，8種重金屬元素背景值都遠(yuǎn)低于國家土壤環(huán)境質(zhì)量一級標(biāo)準(zhǔn)，與黑龍江省土壤背景值相比，除了重金屬As、Ni的背景值略高于黑龍江省土壤背景值外，其余重金屬元素背景值均低于黑龍江省土壤背景值，與1986年測得的松嫩平原土壤背景值相比，除了重金屬Hg的背景值略低于1986年松嫩平原背景值外，其余重金屬背景值均高于1986年測得的松嫩平原土壤重金屬背景值。

表1 地累積指數(shù)分級標(biāo)準(zhǔn)

表2 表層土壤重金屬描述性統(tǒng)計

表3 深層土壤重金屬描述性統(tǒng)計

表4 不同方法計算的土壤重金屬元素背景值

表5 土壤重金屬背景值比較

注: A,松嫩平原黑龍江省南部土壤背景值； B,國家土壤環(huán)境質(zhì)量一級標(biāo)準(zhǔn)； C,黑龍江省土壤背景值； D,1986年松嫩平原土壤背景值。

Note: A,Soil background level in southern Songnen Plain of Heilongjiang Province; B,National soil environmental quality standards; C,Soil background level in Heilongjiang Province; D,Soil background level of Songnen Plain in 1986.

3.2 土壤重金屬空間分布

3.2.1 變異函數(shù)擬合 地統(tǒng)計學(xué)是以區(qū)域化變量理論作為依據(jù)，能夠反映區(qū)域化變量的許多特征[18]。采用GS+9.0軟件對松嫩平原黑龍江省南部深層土壤(180～200 cm)樣本含量數(shù)據(jù)進(jìn)行擬合，確定擬合的最優(yōu)模型，擬合結(jié)果及相關(guān)參數(shù)見表6。重金屬As、Hg、Cd、Cr、Zn變異函數(shù)擬合的結(jié)果為指數(shù)模型，重金屬Cu、Ni變異函數(shù)擬合的結(jié)果為球狀模型，重金屬Pb變異函數(shù)擬合的結(jié)果為高斯模型。塊金值C0反映了區(qū)域化變量隨機(jī)性大??；基臺值(C0+C)反映了區(qū)域化變量變化的幅度；塊金系數(shù)C0/(C0+C)通常能反映區(qū)域化變量的空間自相關(guān)性；變程反映了區(qū)域化變量的影響范圍；擬合系數(shù)(r2)反映了模型擬合的精度[19]，擬合系數(shù)越接近1表示擬合的精度越高。

重金屬As、Hg、Cr、Cu、Ni、Zn塊金系數(shù)均小于0.25，具有強(qiáng)烈的空間自相關(guān)性，重金屬Cd、Pb的塊金系數(shù)介于0.25～0.75之間，具有中等程度的空間自相關(guān)性，表明8種重金屬元素含量均主要受結(jié)構(gòu)性因素(地質(zhì)、地貌類型條件等)的影響。

3.2.2 普通克里金插值 在變異函數(shù)擬合的基礎(chǔ)上，根據(jù)擬合確定的最優(yōu)模型，應(yīng)用ArcGIS 10.2軟件對深層土壤(180～200 cm)重金屬樣本含量數(shù)據(jù)進(jìn)行空間插值，插值的結(jié)果見圖2。從圖2中可以看出重金屬As、Cd、Cr、Cu、Ni、Zn在研究區(qū)的分布相對一致，呈現(xiàn)出東北高西南低的空間分布格局，特別是重金屬Cr、Cu東北部和西南部差異更為明顯；重金屬Hg、Pb在研究區(qū)分布相對較為均勻，重金屬Hg僅僅在研究區(qū)東部小范圍內(nèi)含量較高，重金屬Pb在研究區(qū)西北部小范圍內(nèi)含量較高。

3.3 表層土壤重金屬污染程度研究

采用地累積指數(shù)法對松嫩平原黑龍江省南部表層土壤(0～20 cm)重金屬污染程度進(jìn)行研究，以確定的松嫩平原黑龍江省南部土壤重金屬背景值作為評價標(biāo)準(zhǔn)，8種重金屬處于不同污染程度的樣本數(shù)量見表7。重金屬As、Hg、Cd、Cr、Cu、Ni、Pb、Zn處于無污染程度的樣本數(shù)量分別占總樣本數(shù)量的97.72%、59.92%、83.45%、97.63%、99.24%、98.86%、99.02%、99.26%，只有重金屬Hg、Cd無污染的樣本量占全部樣本量的比例低于95%；重金屬As、Hg、Cd、Cr、Cu、Ni、Pb、Zn處于無污染到中度污染的樣本數(shù)量占全部樣本量的比例分別為2.27%、36.00%、15.61%、2.37%、0.76%、1.14%、0.95%、0.73%，除重金屬Hg、Cd污染樣本占全部樣本的比例較高之外，其余重金屬所占比例都很低；重金屬Hg、Cd處于中度污染的樣本占全部樣本量的比例分別為3.18%、0.77%，除此之外，重金屬處于中度污染到強(qiáng)污染、強(qiáng)污染、強(qiáng)污染到極強(qiáng)污染、極強(qiáng)污染的樣本量都很少，可以忽略不計。

表6 變異函數(shù)擬合結(jié)果及相關(guān)參數(shù)

圖2 深層土壤重金屬元素含量插值圖Fig.2 Interpolation of heavy metal content in deep soil

污染程度 Pollution levelAsHgCdCrCuNiPbZn無污染Uncontaminated2044612537174612042720763206842071820769無污染到中度污染Uncontaminated to moderately contaminated47575323267496158239199153中度污染Moderately contaminated266516202061中度污染到強(qiáng)污染Moderately toseverly contaminated01272200000強(qiáng)污染Severelycontaminated040900000強(qiáng)污染到極強(qiáng)污染Severelyto extremely contaminated014000000極強(qiáng)污染Extremely contaminated08200000

總體而言，松嫩平原黑龍江省南部表層土壤重金屬As、Cr、Cu、Ni、Pb、Zn處于無污染程度，是相對清潔的，但重金屬Hg和Cd存在著一定程度的輕微污染，應(yīng)當(dāng)引起相關(guān)部門的注意，同時也需要加強(qiáng)對重金屬Hg和Cd的污染控制與管理。

4 結(jié)論與討論

4.1 結(jié)論

(1)采用置信區(qū)間法和累積頻率法得到的松嫩平原黑龍江省南部8種土壤重金屬元素As、Hg、Cd、Cr、Cu、Ni、Pb、Zn的背景值分別為9.33、0.016、0.079、50.82、18.74、24.16、22.70、57.34(mg·kg-1)。

(2)普通克里金插值的結(jié)果表明，深層土壤重金屬As、Cd、Cr、Cu、Ni、Zn在研究區(qū)的分布相對一致，呈現(xiàn)出東北高西南低的空間分布格局，重金屬Hg、Pb在研究區(qū)分布相對較為均勻。

(3)地累積指數(shù)法評估的結(jié)果表明，松嫩平原黑龍江省南部表層土壤總體而言是相對清潔的，但重金屬Hg、Cd出現(xiàn)了一定程度的輕微污染。

4.2 討論

背景值作為土壤重金屬污染評價的重要參考標(biāo)準(zhǔn)，其值的準(zhǔn)確性決定了土壤污染狀況判定的精度，尤其是對于較小范圍區(qū)域內(nèi)的土壤重金屬污染狀況進(jìn)行研究時，更應(yīng)該選取準(zhǔn)確的背景值以刻畫區(qū)域土壤重金屬含量狀況，在實際條件允許的前提下，應(yīng)當(dāng)采集深層土壤樣本以確定區(qū)域背景值。同時應(yīng)該在本研究的基礎(chǔ)上考慮土壤重金屬背景值的影響因素，探索不同地質(zhì)、地貌等條件下的土壤重金屬背景值，以實現(xiàn)從多角度、全面分析土壤重金屬背景值特征，為土壤重金屬背景值、污染評價等研究提供參考。