張善紅+李堆淑

摘要 土壤重金屬污染問(wèn)題是當(dāng)今學(xué)術(shù)界、政府和公眾廣泛關(guān)注的問(wèn)題之一。本文總結(jié)了GIS在土壤重金屬污染評(píng)價(jià)中制作環(huán)境質(zhì)量專題圖、空間分析重金屬污染、指導(dǎo)土壤采樣方案以及管理土壤數(shù)據(jù)庫(kù)等方面的應(yīng)用，并預(yù)測(cè)了GIS在土壤重金屬污染評(píng)價(jià)應(yīng)用方面的發(fā)展趨勢(shì)，以供相關(guān)人員參考。

關(guān)鍵詞 重金屬；GIS；污染評(píng)價(jià)；研究進(jìn)展

中圖分類號(hào) X825 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1007-5739（2017）22-0151-01

重金屬是一類典型的環(huán)境污染物，土壤重金屬污染具有危害周期長(zhǎng)、范圍廣、毒性大、生物難降解、污染隱蔽等特點(diǎn)[1]。在自然條件下，除了部分土壤的重金屬含量較高之外，一般土壤中有害重金屬含量都小于構(gòu)成污染的臨界值。但由于受人類生產(chǎn)活動(dòng)的影響，土壤污染加重。加上近幾十年城市化、工業(yè)化的發(fā)展，土壤中重金屬污染物含量不斷提高，土壤污染問(wèn)題越來(lái)越突出，重金屬污染問(wèn)題已引起學(xué)術(shù)界、政府和公眾的廣泛關(guān)注[2]。

地理信息系統(tǒng)（GIS）技術(shù)具有強(qiáng)大的空間分析功能，集空間數(shù)據(jù)的采集、分析、管理、可視化輸出為一體[3]，可在可視化的環(huán)境下對(duì)土壤重金屬含量空間分布進(jìn)行研究[4]。基于GIS的土壤重金屬研究，一般利用GIS軟件的計(jì)算和圖形顯示功能。利用克里格等插值方法，將采樣點(diǎn)得到的點(diǎn)源數(shù)據(jù)轉(zhuǎn)化成面源數(shù)據(jù)[5]，同時(shí)將土壤圖數(shù)字化后建立起相應(yīng)的空間與屬性的數(shù)據(jù)庫(kù)，從而繪制出土壤重金屬空間分布圖、重金屬污染評(píng)價(jià)圖等。該方法可以重構(gòu)土壤重金屬的空間分布，實(shí)現(xiàn)土壤重金屬污染的可視化；利用GIS技術(shù)的空間分析功能，可以判斷土壤重金屬的空間分布特征及其與周圍環(huán)境之間的相關(guān)關(guān)系，有助于查明土壤環(huán)境質(zhì)量在區(qū)域空間內(nèi)的變化結(jié)果、原因及趨勢(shì)，為研究者或政府提供一個(gè)客觀、準(zhǔn)確的分析依據(jù)。本文從以下4個(gè)方面介紹GIS在土壤重金屬污染評(píng)價(jià)中的應(yīng)用。

1 制作環(huán)境質(zhì)量專題圖

GIS的基本功能是制圖，即根據(jù)已知點(diǎn)狀數(shù)據(jù)的屬性，利用不同的空間插值方法繪制等值線圖和空間分布圖，并按照不同的數(shù)值范圍進(jìn)行設(shè)色，更加直觀、形象地反映土壤重金屬污染的空間分布特點(diǎn)。目前，使用GIS中的克里格插值法制圖展示重金屬在土壤中的空間分布已經(jīng)被廣泛應(yīng)用[6]。郝麗虹等[7]運(yùn)用克里格插值方法中的普通插值法，對(duì)海南島農(nóng)用地中的Hg、Cd、Cr、Pb、As 5種重金屬含量空間分布進(jìn)行可視化處理，研究結(jié)果表明海南島工業(yè)區(qū)外源性重金屬污染最為嚴(yán)重；王曉輝等[8]采用普通克里格和析取克里格插值方法制作了淮南一礦區(qū)的重金屬空間分布圖，結(jié)果表明第1類元素的超標(biāo)部分主要位于2個(gè)重礦區(qū)，第2類元素超標(biāo)部分位于老城區(qū)；李春芳等[9]基于GIS空間插值相結(jié)合的方法研究了龍口市9種重金屬的空間結(jié)構(gòu)和分布特征；劉曉雙[10]運(yùn)用GIS中普通克里格方法和綜合指數(shù)法研究了云浮硫鐵礦區(qū)6種重金屬的空間分布狀況以及污染狀況評(píng)價(jià)圖；王學(xué)軍等[11]分析了北京東郊污灌區(qū)農(nóng)用土壤表層重金屬含量，并制作了Cu、Zn、Cr、Pb 4種重金屬含量的三維空間分布圖及污染評(píng)價(jià)指數(shù)分級(jí)圖。

除了制作土壤重金屬的空間分布圖、污染評(píng)價(jià)圖外，也有很多學(xué)者基于GIS將土地利用圖、地形圖及相關(guān)地理底圖與專題地圖進(jìn)行比較分析，用于揭示重金屬的空間分布與周圍地理環(huán)境的相關(guān)關(guān)系。如Simasuwannarong B等[12]將土壤重金屬空間分布圖與工廠分布、交通線路分布、農(nóng)業(yè)用地等地理底圖進(jìn)行分析比較，揭示農(nóng)業(yè)用地重金屬污染與工廠廢物排放、汽車尾氣排放的相關(guān)關(guān)系；蔣詩(shī)泉等[13]將地貌特征圖與重金屬空間分布圖相結(jié)合制作了三維地貌上重金屬的空間分布，更精確地探究重金屬污染的原因。

克里格插值方法通常是將有限的點(diǎn)源數(shù)據(jù)轉(zhuǎn)化成面源數(shù)據(jù)，為土壤重金屬的空間分布預(yù)測(cè)提供了一種線性、無(wú)偏、最優(yōu)估值，反映土壤污染物的二維分布情況。但克里格插值方法受到變異函數(shù)、待估樣點(diǎn)的幾何性質(zhì)、已知樣點(diǎn)數(shù)據(jù)結(jié)構(gòu)等的影響，此外還要求使用該方法的數(shù)據(jù)呈正態(tài)分布，且數(shù)據(jù)之間具有一定的相關(guān)性。由于重金屬空間差異性受人為影響較大，如汽車尾氣排放、工廠“三廢”排放、垃圾填埋場(chǎng)位置等，因而利用該方法做出的重金屬空間分布圖容易模糊重金屬污染污染區(qū)域的邊界，將污染變異區(qū)域擴(kuò)大化。近年來(lái)，GIS與數(shù)學(xué)模型的集成在土壤重金屬污染評(píng)價(jià)中已相當(dāng)普遍，并成為重金屬污染評(píng)價(jià)的一個(gè)新趨勢(shì)，如模糊數(shù)學(xué)、序貫高斯、回歸模型、高精度曲面建模等，這些方法考慮因素越來(lái)越多，使評(píng)價(jià)更加靈活、客觀。

2 進(jìn)行土壤重金屬污染空間分析

利用GIS的疊加分析、緩沖區(qū)分析、數(shù)字地形分析等空間分析功能，將土壤重金屬污染的數(shù)據(jù)與研究區(qū)的其他背景數(shù)據(jù)進(jìn)行疊加，從而產(chǎn)生新的不同屬性數(shù)據(jù)，且這些數(shù)據(jù)綜合了原來(lái)多個(gè)要素的性質(zhì)，能更好地反映土壤污染與周圍環(huán)境以及各種因素之間的關(guān)系。例如趙彥鋒等[14]將工廠分布圖、鎮(zhèn)區(qū)緩沖圖與重金屬空間分布圖進(jìn)行疊加，得到土壤Pb含量隨與城鎮(zhèn)距離的增加呈持續(xù)下降趨勢(shì)，土壤Cu、Zn、Cd含量隨與工廠距離的增加呈下降趨勢(shì)，但土壤重金屬含量與距離公路的遠(yuǎn)近、土地利用類型的關(guān)系不明顯；史文嬌等[15]將雙城市空間分布圖與工廠、道路、城鎮(zhèn)、土壤類型、土壤理化性質(zhì)、高程等數(shù)據(jù)進(jìn)行疊加，分析重金屬空間分異的原因；謝小進(jìn)[16]將土壤空間分布圖與不同類型工業(yè)區(qū)、不同級(jí)別公路緩沖區(qū)進(jìn)行疊加，分析重金屬污染的原因。

3 指導(dǎo)土壤采樣方案

采集方案（采集的數(shù)目和密度）對(duì)GIS制作重金屬空間分布圖質(zhì)量有很大影響。一個(gè)理想的采集方案包括合理的采集布點(diǎn)方法與合適的采集密度[17]。較低的采集樣點(diǎn)密度會(huì)降低空間結(jié)構(gòu)的描述能力，增加預(yù)測(cè)誤差。因此，采集樣點(diǎn)越多，重金屬空間分布圖就越準(zhǔn)確?？死锔癫逯捣椒ㄖ辽傩枰?0～100個(gè)采樣點(diǎn)才能獲得穩(wěn)定的半方差函數(shù)[18]，才能更加準(zhǔn)確地描繪空間結(jié)構(gòu)；有的學(xué)者建議至少100～150個(gè)采樣點(diǎn)才能建立一個(gè)平穩(wěn)的變異函數(shù)[19]，至少300個(gè)樣點(diǎn)才能分析各向異性[20]。此外，采集樣點(diǎn)的位置也會(huì)影響插值的準(zhǔn)確性。采集樣點(diǎn)均勻分布可以更好地反映插值要素在空間上的分布，能更精確地建立半方差函數(shù)，提高空間分布圖的準(zhǔn)確性。因此，為了最大限度地節(jié)約野外采集時(shí)間成本和室內(nèi)實(shí)驗(yàn)分析成本，在保證精度的前提下，GIS可以指導(dǎo)土壤采集方案，減少采集樣點(diǎn)的個(gè)數(shù)。

4 土壤數(shù)據(jù)庫(kù)管理

長(zhǎng)時(shí)間監(jiān)測(cè)以及大區(qū)域土壤重金屬污染研究需要建立專門的土壤數(shù)據(jù)庫(kù)，土壤數(shù)據(jù)庫(kù)是進(jìn)行土壤環(huán)境質(zhì)量動(dòng)態(tài)分析、監(jiān)測(cè)和評(píng)價(jià)的基礎(chǔ)。近些年有學(xué)者開(kāi)始建立數(shù)據(jù)庫(kù)，并進(jìn)行相應(yīng)的研究。如李衛(wèi)江等[21]構(gòu)建了一個(gè)集采樣點(diǎn)監(jiān)測(cè)數(shù)據(jù)與元數(shù)據(jù)、空間數(shù)據(jù)為一體的上海市基本農(nóng)田土壤重金屬數(shù)據(jù)庫(kù)，對(duì)土壤環(huán)境的長(zhǎng)期管理、農(nóng)產(chǎn)品安全生產(chǎn)與管理提供重要數(shù)據(jù)依據(jù)。

5 結(jié)語(yǔ)

GIS在土壤重金屬污染評(píng)價(jià)中的應(yīng)用主要集中在制作專題地圖、進(jìn)行土壤重金屬污染空間分析、指導(dǎo)土壤采樣方案及土壤數(shù)據(jù)庫(kù)管理等。隨著GIS技術(shù)的普及以及相關(guān)技術(shù)的發(fā)展，在可視化環(huán)境下，從時(shí)空方面對(duì)土壤重金屬污染研究及土壤重金屬污染擴(kuò)散研究是今后研究的重點(diǎn)。

6 參考文獻(xiàn)

[1] 顧濟(jì)滄，趙娟.云南省土壤重金屬污染現(xiàn)狀及治理技術(shù)研究[J].環(huán)境科學(xué)導(dǎo)刊，2010，29（5）：68-71.

[2] 周燕，李德志，宋云，等.崇明島土壤重金屬含量的空間分布及污染評(píng)價(jià)[J].華東師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2009，5（5）：75-84.

[3] TAM N F Y，WONG Y S.Spatia lvariation of heavy metals insurface sediments of Hong Kong Mangrove swamps[J].Environrnental Pollution，2000，110（2）：195- 200.

[4] 劉慶，杜志勇，史衍璽，等.基于GIS的山東壽光蔬菜產(chǎn)地土壤重金屬空間分布特征[J].農(nóng)業(yè)工程學(xué)報(bào)，2009，25（10）：258-263.

[5] 馮經(jīng)昆.西南山區(qū)某垃圾焚燒廠周邊土壤重金屬空間分布及其污染評(píng)價(jià)[D].桂林：廣西師范大學(xué)，2014.

[6] 陳峰，尹春芹，蔣新，等.基于GIS的南京市典型蔬菜基地土壤重金屬污染現(xiàn)狀與評(píng)價(jià)[J].中國(guó)環(huán)境監(jiān)測(cè)，2008，24（2）：40-44.

[7] 郝麗虹，張東明，吳鵬飛，等.海南島農(nóng)用地土壤重金屬含量的空間分布[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2009，17（2）：203-234.

[8] 王曉輝，楊晨.基于GIS和地統(tǒng)計(jì)學(xué)的淮南礦區(qū)土壤重金屬含量與空間分布研究[J].長(zhǎng)江流域資源與環(huán)境，2014，23（11）：60-65.

[9] 李春芳，王菲，曹文濤，等.龍口市污水灌溉區(qū)農(nóng)田重金屬來(lái)源、空間分布及污染評(píng)價(jià)[J].環(huán)境科學(xué)，2017，38（3）：1018-1027.

[10] 劉曉雙.基于GIS的云浮硫鐵礦周邊土壤中重金屬污染空間分布于污染評(píng)價(jià)研究[D].廣州：暨南大學(xué)，2010.

[11] 王學(xué)軍，席爽.北京東郊污灌土壤重金屬含量的克立格插值及重金屬污染評(píng)價(jià)[J].中國(guó)環(huán)境科學(xué)，1997，17（3）：225-228.

[12] SIMASUWANNARONG B，SATAPANAJARU T，KHUNTONG S，et al.Spatial distribution and risk assessment of As，Cd，Pb and Zn in Topsoil at Rayong Province，Thailand[J].Water Air Soil Pollut，2012（223）：1931-1943.

[13] 蔣詩(shī)泉，劉中俠，蔣詩(shī)平.基于綜合主成分模型的城市土壤重金屬空間分布與傳播特征研究[J].華東師范大學(xué)學(xué)報(bào)（自然科學(xué)版），2015（3）：136-145.

[14] 趙彥鋒，史學(xué)正，于東升，等.工業(yè)型城鄉(xiāng)交錯(cuò)區(qū)農(nóng)業(yè)土壤Cu、Zn、Pb和Cd的空間分布及影響因素研究[J].土壤學(xué)報(bào)，2007，44（2）：227-234.

[15] 史文嬌，魏開(kāi)，汪景寬，等.雙城市土壤重金屬空間分異及影響因子分析[J].水土保持學(xué)報(bào)，2007，21（1）：59-64.

[16] 謝小進(jìn).上海地區(qū)土壤重金屬空間分布特征及其成因分析[D].上海：上海師范大學(xué)，2010.

[17] 張軍，余江，陳文清.基于GIS的土壤重金屬污染評(píng)價(jià)研究[J].安徽農(nóng)業(yè)科學(xué)，2011，39（29）：18106-18107.

[18] WEBSTER R，OLIVER M A.Sample adequately to estimate variograms of soil properties[J].Journal of Soil Science，1992，43（1）：177-192.

[19] VOLTZ M，WEBSTER R.A comparison of Kriging，cubic splines and classification for predicting soil properties from sample information[J].European Journal of Soil Science，1990，41（3）：473-490.

[20] ROBINSON T P，METTERNICHT G.Testing the performance of spatial interpolation techniques for mapping soil properties[J].Computers and Electronics in Agriculture，2006，50（2）：97-108.

[21] 李衛(wèi)江，康建成，吳永興.基于GIS的土壤重金屬數(shù)據(jù)庫(kù)構(gòu)建及應(yīng)用[J].測(cè)繪科學(xué)，2010，35（4）：184-186.