城市表層土壤重金屬污染分析

1、城市表層土壤重金屬污染分析摘 要隨著城市工業(yè)化水平的加快，除了醫(yī)療器械、冶金、化工生產(chǎn)、電子等諸多行業(yè)有大量重金屬排出，人們的日常生活，農(nóng)業(yè)用藥以及交通尾氣等也使排污量的大量增加，汞、鉛、鉻、銅等重金屬污染已成為一個(gè)世界性難題。土壤中重金屬的富集會(huì)導(dǎo)致環(huán)境污染，對(duì)植物生長(zhǎng)安全及人體健康等構(gòu)成直接危害。因此研究城市土壤中重金屬的分布與污染問(wèn)題，有著越來(lái)越重要的現(xiàn)實(shí)意義。針對(duì)問(wèn)題一，利用Matlab軟件插值繪制出城市功能區(qū)分布圖、海拔圖和重金屬元素等濃度曲線圖，得到了各重金屬元素的空間分布；運(yùn)用內(nèi)梅羅綜合指數(shù)結(jié)合國(guó)家土壤重金屬含量二級(jí)標(biāo)準(zhǔn)，對(duì)比分析得出結(jié)論：工業(yè)區(qū)、交通區(qū)的重金屬污染程度較大。針對(duì)

2、問(wèn)題二，基于統(tǒng)計(jì)學(xué)方法，使用平均值、標(biāo)準(zhǔn)差、變異系數(shù)、偏度和峰度，分析8種元素在各功能區(qū)的污染情況和分布特點(diǎn)，得到重金屬污染的主要原因是工廠廢物排放和汽車尾氣的排放。針對(duì)問(wèn)題三，基于氣體擴(kuò)散理論，建立重金屬污染物的指數(shù)分布模型；運(yùn)用ArcGIS 軟件繪制各重金屬元素濃度的直方圖、QQ-PLOT圖、空間趨勢(shì)圖以及半變異函數(shù)圖像，分析出重金屬污染物的傳播特征，從而選定出指數(shù)范普通克里格插值為重金屬元素的最優(yōu)插值模型。利用Matlab軟件對(duì)插值生成的數(shù)據(jù)進(jìn)行分析，畫出梯度矢量圖，準(zhǔn)確確定出8種重金屬污染源的位置。針對(duì)問(wèn)題四，在問(wèn)題三的模型基礎(chǔ)上引入時(shí)間變量，通過(guò)高斯公式來(lái)研究城市地質(zhì)環(huán)境的演變模型，

3、擬需要不同時(shí)間段的數(shù)據(jù)材料。關(guān)鍵詞：土壤重金屬污染 空間分布 地統(tǒng)計(jì)分析 ArcGIS 內(nèi)梅羅指數(shù) 克里格插值一 問(wèn)題重述隨著城市經(jīng)濟(jì)的快速發(fā)展和城市人口的不斷增加，人類活動(dòng)對(duì)城市環(huán)境質(zhì)量的影響日顯突出。對(duì)城市土壤地質(zhì)環(huán)境異常的查證，以及如何應(yīng)用查證獲得的海量數(shù)據(jù)資料開展城市環(huán)境質(zhì)量評(píng)價(jià)，研究人類活動(dòng)影響下城市地質(zhì)環(huán)境的演變模式，日益成為人們關(guān)注的焦點(diǎn)。按照功能劃分，城區(qū)一般可分為生活區(qū)、工業(yè)區(qū)、山區(qū)、主干道路區(qū)及公園綠地區(qū)等，分別記為1類區(qū)、2類區(qū)、5類區(qū)，不同的區(qū)域環(huán)境受人類活動(dòng)影響的程度不同。將所考察的城區(qū)劃分為間距1公里左右的網(wǎng)格子區(qū)域，按照每平方公里1個(gè)采樣點(diǎn)對(duì)表層土（010 厘米深

4、度）進(jìn)行取樣、編號(hào)，并用GPS記錄采樣點(diǎn)的位置。應(yīng)用專門儀器測(cè)試分析，獲得了每個(gè)樣本所含的多種化學(xué)元素的濃度數(shù)據(jù)。另一方面，按照2公里的間距在那些遠(yuǎn)離人群及工業(yè)活動(dòng)的自然區(qū)取樣，將其作為該城區(qū)表層土壤中元素的背景值。請(qǐng)運(yùn)用數(shù)學(xué)建模的分析方法來(lái)完成以下任務(wù)：(1) 給出8種主要重金屬元素在該城區(qū)的空間分布，并分析該城區(qū)內(nèi)不同區(qū)域重金屬的污染程度。(2) 通過(guò)數(shù)據(jù)分析，說(shuō)明重金屬污染的主要原因。(3) 分析重金屬污染物的傳播特征，由此建立模型，確定污染源的位置。(4) 分析你所建立模型的優(yōu)缺點(diǎn)，為更好地研究城市地質(zhì)環(huán)境的演變模式，還應(yīng)收集什么信息？有了這些信息，如何建立模型解決問(wèn)題？二 模型假設(shè)1

5、、忽略海拔對(duì)重金屬濃度的影響。2、所給數(shù)據(jù)無(wú)壞值。3、取樣點(diǎn)分布均勻，各區(qū)域的個(gè)數(shù)可以代表區(qū)域面積。4、該城市內(nèi)的土壤ph值呈中性且為旱地。5、重金屬元素隨廢氣、粉塵傳播的速率為常數(shù)。三 符號(hào)說(shuō)明符號(hào)意義第種重金屬污染元素污染物的環(huán)境實(shí)測(cè)值污染物的評(píng)價(jià)標(biāo)準(zhǔn)P綜j第j個(gè)功能區(qū)的內(nèi)梅羅指數(shù)土壤污染中污染指數(shù)的最大值土壤污染中污染指數(shù)的平均值C重金屬元素的濃度擴(kuò)散系數(shù)grad梯度q單位時(shí)間通過(guò)單位法向長(zhǎng)度的流量Q總流量p0內(nèi)污染物濃度的單位面密度k重金屬元素污染的傳播速率四 模型的建立與求解4.1任務(wù)一的分析與解決4.1.1任務(wù)一的分析任務(wù)一的關(guān)鍵是如何給出8種主要重金屬元素在該城區(qū)的空間分布，并分

6、析該城區(qū)內(nèi)不同區(qū)域重金屬的污染程度。論文中運(yùn)用ArcGIS軟件插值繪制出8種重金屬元素在該城區(qū)的空間分布圖像，并按劃分好的5種功能區(qū)，對(duì)各區(qū)域中的重金屬污染物程度進(jìn)行了定量描述。為了全面反映各污染物對(duì)土壤的不同作用，突出高濃度污染物對(duì)環(huán)境質(zhì)量的影響，克服平均值法中各個(gè)污染物分擔(dān)的缺陷，論文中采用了內(nèi)梅羅綜合污染指數(shù)法，能較全面地反映環(huán)境質(zhì)量水平。4.1.2任務(wù)一的解決論文中利用已知的取樣點(diǎn)坐標(biāo)值及其所在的功能區(qū)域編號(hào)，運(yùn)用Matlab軟件繪制出該城市的功能區(qū)域分布圖像及海拔圖（見附錄一），圖中用不同顏色表示了城市中的5 個(gè)功能區(qū)域（如下圖）。通過(guò)對(duì)功能區(qū)分布圖和8種元素的濃度分布圖進(jìn)行比較分析

7、，我們可以得到各元素在市內(nèi)的分布情況，具體分析如下： 圖1 功能區(qū)分布圖圖3 Cd濃度分布圖圖2 As濃度分布圖由上圖我們可以看出受As元素影響較大的區(qū)域主要集中在山區(qū)、工業(yè)區(qū)以及交通區(qū)，并且從圖像中還可以得知As的分布范圍比較廣泛；Cd元素影響較大的區(qū)域主要集中在山區(qū)、交通區(qū)，部分生活區(qū)由于靠近交通區(qū)也受到了較大影響。此外，從曲線的分布范圍來(lái)看，該元素分布范圍十分廣泛，幾乎遍及了這個(gè)城市大部分地區(qū)?？梢酝茰y(cè)這兩種元素的污染源主要是非點(diǎn)源污染，而且這兩種元素?cái)U(kuò)散性較強(qiáng)。圖5 Cu濃度分布圖圖4 Cr濃度分布圖 上圖顯示，Cr元素主要集中在部分交通區(qū)與生活區(qū)交界地帶，山區(qū)和交通區(qū)的部分地方也受到

8、了Cr元素的影響。圖像中元素濃度曲線較為集中，說(shuō)明此種元素在城市內(nèi)的分布主要集中在特定的某幾塊區(qū)域，屬于非均勻分布；Cu元素主要集中在部分交通區(qū)與生活區(qū)交界地帶以及工業(yè)區(qū)和交通區(qū)的局部。與Cr的分布情況相類似，該元素的分布也主要集中在某幾個(gè)特定的區(qū)域，呈現(xiàn)一種較不均勻的分布狀態(tài)。可以推測(cè)這兩種元素的污染源主要是點(diǎn)源污染，擴(kuò)散性較弱。圖7 Ni濃度分布圖圖6 Hg濃度分布圖 該圖顯示，受高濃度Hg影響的區(qū)域較多，各個(gè)功能區(qū)均出現(xiàn)被污染跡象，但該元素的分布主要集中在功能區(qū)的部分區(qū)域，沒有出現(xiàn)大范圍的傳播或擴(kuò)散現(xiàn)象；Ni元素在山區(qū)、生活區(qū)和交通區(qū)的交界地帶有較高濃度的分布。元素的污染濃度曲線分布范圍

9、十分廣泛，說(shuō)明這種元素?cái)U(kuò)散能力較強(qiáng)。圖9 Zn濃度分布圖圖8 Pb濃度分布圖上圖顯示Pb的分布范圍較為廣泛，受其影響較大的主要是交通區(qū)、工業(yè)區(qū)和生活區(qū)；Zn主要富集在交通區(qū)和工業(yè)區(qū)，在人為活動(dòng)影響較小的山區(qū)分布不是很多。 以上分別定性討論了8種重金屬元素空間分布，該城市的環(huán)境綜合污染程度與8種重金屬污染物質(zhì)相關(guān)。因此，需要一種同時(shí)考慮多種污染物綜合污染水平的多因子評(píng)價(jià)方法。我們引入了內(nèi)梅羅綜合污染歸一化指數(shù)評(píng)價(jià)體系，它兼顧了單因子污染指數(shù)的平均值和最高值，能較全面地反映環(huán)境質(zhì)量水平，而且可以突出污染較重的污染物的作用，其計(jì)算公式如下： ，式中為污染物的環(huán)境實(shí)測(cè)值（ug·g-1）；為污

10、染物的評(píng)價(jià)標(biāo)準(zhǔn)（ug·g-1）；為土壤污染中污染指數(shù)的最大值；為土壤污染中污染指數(shù)的平均值1-3。通過(guò)上述公式計(jì)算得到的值越大，說(shuō)明該區(qū)域受所有重金屬污染的綜合影響程度越大。我們根據(jù)內(nèi)梅羅評(píng)價(jià)體系利用Excel軟件計(jì)算出每一個(gè)采樣點(diǎn)的P綜值，具體見附錄三。依據(jù)得出的結(jié)果，分區(qū)域計(jì)算出每一個(gè)功能區(qū)的P綜平均值，記為P綜j，所得結(jié)果如表1：功能區(qū) 功能區(qū)域P綜j1生活區(qū)4.692工業(yè)區(qū)15.343山 區(qū)1.684交通區(qū)12.425公 園3.85 表1 各功能區(qū)內(nèi)梅羅指數(shù)由表1的結(jié)果可以得出以下兩點(diǎn)結(jié)論：（1）工業(yè)區(qū)和交通區(qū)受重金屬污染程度較大；（2）山區(qū)是五個(gè)功能區(qū)中環(huán)境質(zhì)量最好的區(qū)域。

11、工業(yè)區(qū)由于集中了大量化工企業(yè)，這些企業(yè)向環(huán)境中排放了大量重金屬污染物，對(duì)環(huán)境的影響很大；在交通區(qū)，汽車排放的尾氣是重金屬元素的主要來(lái)源，其附近環(huán)境也受到較大影響；由于山區(qū)的人為干擾因素最小，而且綠化覆蓋率高,環(huán)境自我調(diào)節(jié)能力強(qiáng),所以受到的污染物影響也就最小。為了進(jìn)一步分析該城市受重金屬污染的主要原因，我們需要對(duì)5個(gè)功能區(qū)的土壤受污染程度進(jìn)行分級(jí)，為此論文中引入了國(guó)家土壤重金屬含量二級(jí)標(biāo)準(zhǔn)（GB15618-1995）。如表2所示：元素As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)標(biāo)準(zhǔn)含量303003002005

12、0050300250 表2 國(guó)家土壤重金屬含量二級(jí)標(biāo)準(zhǔn)土壤污染綜合指數(shù)分級(jí)標(biāo)準(zhǔn)為：綜合污染指數(shù)>3為重污染，23為中污染，12為輕污染，0.71為警戒級(jí)，0.7為安全級(jí)。利用單項(xiàng)污染指數(shù)評(píng)價(jià)方法，對(duì)5大功能區(qū)中的各重金屬元素的污染狀況進(jìn)行量化。最后，由內(nèi)梅羅評(píng)價(jià)方法給出各功能區(qū)的綜合污染指數(shù)。根據(jù)計(jì)算得出的單因子指數(shù)以及土壤污染綜合指數(shù)分級(jí)標(biāo)準(zhǔn)，分別得到各功能區(qū)中輕度、中度和重度的污染比例，由此列出下表所示的數(shù)據(jù)： As (g/g)Cd (ng/g)Cr (g/g)輕度中度重度輕度中度重度輕度中度重度10.00%0.00%0.00%31.82%4.55%2.27%0.00%2.27%0.

13、00%20.00%0.00%0.00%50.00%5.56%8.33%0.00%0.00%0.00%30.00%0.00%0.00%4.55%0.00%0.00%0.00%0.00%0.00%40.72%0.00%0.00%37.68%10.87%2.17%0.72%0.00%0.72%50.00%0.00%0.00%20.00%2.86%5.71%0.00%0.00%0.00%平均0.31%0.00%0.00%29.47%6.27%2.82%0.31%0.31%0.31%Cu (g/g)Hg (ng/g)Ni (g/g)輕度中度重度輕度中度重度輕度中度重度12.27%0.00%0.00%2.

14、27%0.00%0.00%0.00%0.00%0.00%22.78%0.00%2.78%8.33%2.78%8.33%0.00%0.00%0.00%30.00%0.00%0.00%0.00%0.00%0.00%1.52%0.00%0.00%42.17%0.00%0.72%2.17%0.00%3.62%0.00%0.72%0.00%50.00%0.00%0.00%0.00%2.86%0.00%0.00%0.00%0.00%平均1.57%0.00%0.63%2.19%0.63%2.51%0.31%0.31%0.00%Pb (g/g)Zn (g/g)8種重金屬綜合污染程度輕度中度重度輕度中度重度輕度

15、中度重度12.27%0.00%0.00%13.64%2.27%4.55%20.45%4.55%4.55%25.56%0.00%0.00%8.33%2.78%8.33%27.78%11.11%8.33%30.00%0.00%0.00%0.00%0.00%0.00%4.55%0.00%0.00%40.00%0.00%0.00%17.39%3.62%3.62%28.99%7.97%5.80%50.00%0.00%0.00%5.71%0.00%2.86%11.43%8.57%2.86%平均0.94%0.00%0.00%10.97%2.19%3.45%20.69%6.27%4.39%表3 各功能區(qū)單因子

16、污染程度通過(guò)對(duì)表中數(shù)據(jù)的比較可以得到如下結(jié)果：（1）8種重金屬元素對(duì)城市的污染程度依次為：Zn>Cd>Hg>Cu， As, Ni，Cr, Pb的污染程度很小；（2）5個(gè)功能區(qū)受污染程度的大小依次為：工業(yè)區(qū)>交通區(qū)>生活區(qū)>公園綠地區(qū)>山區(qū)。4.2任務(wù)二的分析與解決4.2.1 城市工農(nóng)業(yè)生產(chǎn)中土壤內(nèi)重金屬的主要來(lái)源土壤中重金屬的來(lái)源是有多種途徑的。人類工農(nóng)業(yè)生產(chǎn)活動(dòng)是造成土壤重金屬污染的重要原因。土壤中重金屬的主要來(lái)源有：（1）大氣沉降。大氣中的污染物主要來(lái)源于工業(yè)生產(chǎn)、汽車尾氣排放及汽車輪胎磨損產(chǎn)生的大量含重金屬的有害氣體和粉塵等。大氣中的大多數(shù)重金屬

17、經(jīng)自然沉降和雨淋沉降進(jìn)入土壤，他們主要分布在工礦的周圍和公路、鐵路的兩側(cè)，并向四周擴(kuò)散，隨著距城市的距離加大而降低4。(2)工業(yè)活動(dòng)。人類的工業(yè)活動(dòng)是導(dǎo)致環(huán)境中重金屬增加的重要因素之一，有典型的以污染源為中心的同心圈分布現(xiàn)象5。工業(yè)區(qū)由于污染源分布較為集中，因此，在該區(qū)域內(nèi)土壤重金屬含量一般會(huì)高于正常含量，進(jìn)而導(dǎo)致區(qū)域性土壤重金屬污染。工業(yè)區(qū)建立的時(shí)間越長(zhǎng)，其周圍土壤中富集的重金屬含量也就越大，對(duì)環(huán)境的危害程度也就越大6。工業(yè)活動(dòng)導(dǎo)致的重金屬污染通常與煤、石油等化石燃料的使用有關(guān)。金屬冶煉造成的土壤重金屬污染通常表現(xiàn)為復(fù)合污染7，如Cu冶煉廠附近土壤中不僅Cu的含量增加，同時(shí)Ni的含量也有所增

18、加8。(3)交通運(yùn)輸。機(jī)動(dòng)車尾氣既是城市大氣的主要污染源，也是公路兩側(cè)土壤重金屬污染的主要來(lái)源。汽車尾氣排放、輪胎添加劑中的重金屬元素的釋放均可影響到土壤中鉛、鋅、銅的含量，且這些元素的積累量都與交通流量呈正相關(guān)性。（4）污泥施肥。污泥中含有大量的有機(jī)質(zhì)和氮、磷、鉀等營(yíng)養(yǎng)元素，但同時(shí)也含有大量的金屬元素。隨著大量的市政污泥進(jìn)入農(nóng)田，農(nóng)田中的重金屬含量不斷增高。污泥施肥不僅可導(dǎo)致土壤中Cd、Hg、Cr、Cu、Zn、Ni、Pb含量的增加，而且還會(huì)導(dǎo)致農(nóng)作物中重金屬的積累，對(duì)人體健康產(chǎn)生不利影響。（5）污水灌溉。污水灌溉一般指使用經(jīng)過(guò)一定處理的城市污水灌溉農(nóng)田和草地。城市污水包括工業(yè)廢水、生活污水和

19、商業(yè)污水。由于城市化的迅速發(fā)展，大量的被污染廢水被排入河道，致使城市污水中含有的許多重金屬離子隨著污水灌溉而進(jìn)入土壤。在分布上，往往是緊靠污染源頭和城市工業(yè)區(qū)的土壤污染嚴(yán)重，遠(yuǎn)離污染源頭和城市工業(yè)區(qū)土壤鮮有污染9。近年來(lái)污水灌溉已成為農(nóng)業(yè)灌溉用水的重要組成部分，我國(guó)自60年代至今，污灌面積迅速擴(kuò)大，以北方旱作地區(qū)污灌最為普遍，約占全國(guó)污灌面積的90以上。南方污灌面積僅占6，其余在西北和青藏10。污灌會(huì)導(dǎo)致土壤重金屬Hg、Cd、Cr、As、Cu、Zn、Pb等含量增加。4.2.2 數(shù)據(jù)處理方法論文分別對(duì)8種土壤重金屬元素在5個(gè)功能區(qū)濃度的平均值、標(biāo)準(zhǔn)差、變異系數(shù)、偏度以及峰度進(jìn)行了計(jì)算，以各元素濃

20、度在每個(gè)功能區(qū)的平均值作為衡量污染程度的基本標(biāo)準(zhǔn)；利用標(biāo)準(zhǔn)差來(lái)描述同種元素在五個(gè)功能區(qū)中的分布情況，標(biāo)準(zhǔn)差越高表示該元素在相應(yīng)功能區(qū)中的濃度分布越不均勻；變異系數(shù)定義為標(biāo)準(zhǔn)值與平均值之比，其值反映了采樣總體中各采樣點(diǎn)之間的平均變異程度。統(tǒng)計(jì)值之中的偏度反映了正態(tài)分布雙尾特征，標(biāo)準(zhǔn)正態(tài)分布偏度為0。峰度是用來(lái)衡量分布的集中程度或分布曲線的尖峭程度的指標(biāo)，反映了樣本的集中程度，標(biāo)準(zhǔn)正態(tài)分布峰度為3。當(dāng)峰度大于3時(shí)，表示分布的形態(tài)比正態(tài)分布更尖峭，為尖峰分布，當(dāng)峰度小于3時(shí)，比正態(tài)分布更平緩，為平峰分布。在分析8個(gè)重金屬元素對(duì)各功能區(qū)的影響程度時(shí)，我們利用了單項(xiàng)污染指數(shù)評(píng)價(jià)方法，其計(jì)算公式為： ，式

21、中為土壤中污染物的環(huán)境質(zhì)量指數(shù)；為污染物的實(shí)測(cè)值（ug·g-1）；為污染物的評(píng)價(jià)標(biāo)準(zhǔn)（ug·g-1）。>1表示人類活動(dòng)能引起該元素增加；值越大，人為影響程度越大。4.2.3 數(shù)據(jù)處理及分析1. As元素As平均ug·g-1標(biāo)準(zhǔn)差ug·g-1變異系數(shù)偏度峰度評(píng)價(jià)16.272.15020.34290.40172.6971.74179292927.25144.24430.58531.60396.04322.01427533.609071.5911190.4408671.8212586.204251.1233594168.788.124180.522372

22、1.0054623.593181.58556855.7080433.2383290.5673274.60593931.95451.739921 表4 各功能區(qū)As統(tǒng)計(jì)量數(shù)據(jù)分析：a.從平均值來(lái)看，交通區(qū)的As元素含量（168.7）明顯大于其他區(qū)，說(shuō)明As主要來(lái)自交通區(qū)；b.從變異系數(shù)來(lái)看，各個(gè)區(qū)的數(shù)值（均<1.0）都比較小，說(shuō)明該元素含量在各區(qū)波動(dòng)較小，樣本值都比較接近；c.交通區(qū)偏度在1左右，不是很大，表明樣本分布接近于正態(tài)分布，含量在該區(qū)內(nèi)分布比較均勻；d.交通區(qū)的峰度較大，因此樣本比較集中，進(jìn)一步說(shuō)明分布比較均勻；e交通區(qū)的評(píng)價(jià)值最高，表明對(duì)于As元素，在交通區(qū)人為影響程度比較大；

23、結(jié)論： 重金屬元素As污染主要集中于交通區(qū)，并在該區(qū)呈均勻分布，可能來(lái)源于汽車的尾氣。2. Cd元素Cd平均ng·g-1標(biāo)準(zhǔn)差ng·g-1變異系數(shù)偏度峰度評(píng)價(jià)1289.9614183.68080.63352.01188.35732.2304720282393.1111237.57640.60431.64075.41263.0239316243168.788.124180.5223721.0054623.593181.1716899774360.0145243.3920.6760621.8761618.3271482.7693422525280.5429235.84290.8

24、406661.9328285.8838392.158021978表5 各功能區(qū)Cd統(tǒng)計(jì)量數(shù)據(jù)分析：a.從平均值來(lái)看，Cd在除山區(qū)之外的各功能區(qū)中都有較多地分布，表明各區(qū)內(nèi)均會(huì)產(chǎn)生Cd污染；b.其在各功能區(qū)中的變異系數(shù)較小，說(shuō)明Cd元素含量在各區(qū)波動(dòng)比較小，樣本值比較接近；c.各功能區(qū)偏度均小于2，相對(duì)較小，表明該元素分布大致接近正態(tài)分布，該元素在各區(qū)含量分布均勻；d.各功能區(qū)峰度（均大于3）較大，表明樣本值比較集中，該元素分布比較均勻；e.評(píng)價(jià)指數(shù)中工業(yè)區(qū)最大，表明Cd元素污染最嚴(yán)重的區(qū)域在工業(yè)區(qū)；結(jié)論：重金屬元素Cd在除山區(qū)之外的各功能區(qū)分布均勻，推測(cè)污染主要原因是工業(yè)區(qū)中電鍍、采礦、冶煉、

25、染料、電池和化學(xué)工業(yè)等排放的廢水，交通運(yùn)輸中尾氣排放。生活區(qū)與公園綠地區(qū)Cd則主要來(lái)源于工業(yè)區(qū)排放廢氣顆粒的大氣沉降或廢水的傳播。3. Cr元素Cr平均ug·g-1標(biāo)準(zhǔn)差ug·g-1變異系數(shù)偏度峰度評(píng)價(jià)169.0184107.89131.56325.782636.58262.226400293253.409244.00220.82394.229822.86791.722876344339.3902327.899310.708282.89170613.57761.256764418458.0539181.605531.4056858.97875692.78811.872706

26、872543.63614.840130.3400891.5196176.6721211.407612903表6 各功能區(qū)Cr統(tǒng)計(jì)量數(shù)據(jù)分析：a從平均值來(lái)看，Cr在各功能區(qū)的含量比較接近，其中生活區(qū)、交通區(qū)和工業(yè)區(qū)較多，公園綠地和山區(qū)含量在背景值范圍內(nèi)；b從變異系數(shù)來(lái)看，生活區(qū)與交通區(qū)的變異系數(shù)值較大（均大于1），表明樣本值的波動(dòng)較大，存在偏差較大的樣本點(diǎn)，可能是由于某些特殊的樣本值；c各功能區(qū)的偏度都比較大，尤其是生活區(qū)與交通區(qū)，表明Cr含量在該功能區(qū)內(nèi)分布很不均勻；d各功能區(qū)的峰度也都比較大，尤其是交通區(qū)，表明Cr含量在交通區(qū)樣本值比較集中，分布相對(duì)均勻；e從評(píng)價(jià)指數(shù)看出，生活區(qū)的Cr含量受

27、人為影響較大；結(jié)論:重金屬元素Cr分布不均勻，Cr含量較高的相對(duì)集中在某區(qū)域，推測(cè)該區(qū)域附近存在有關(guān)金屬加工、電鍍的工業(yè)部門，產(chǎn)生大量Cr元素，由于大氣沉降等原因，對(duì)生活區(qū)污染十分嚴(yán)重。4. Cu元素Cu平均ug·g-1標(biāo)準(zhǔn)差ug·g-1變異系數(shù)偏度峰度評(píng)價(jià)149.403247.16290.95472.33549.08273.7426652892127.5358414.94133.25355.607332.94499.661805556315.191866.9753260.4591490.6275592.4890971.311914601462.21493120.22221

28、.9323699.406373101.49794.713252086530.1917122.681960.7512643.78356319.129692.287251082表7 各功能區(qū)Cu統(tǒng)計(jì)量數(shù)據(jù)分析：a從平均值來(lái)看，工業(yè)區(qū)的Cu含量明顯較高，表明Cu主要來(lái)源于工業(yè)區(qū)；b從變異系數(shù)來(lái)看，工業(yè)區(qū)的變異系數(shù)（大于1）較大，表明樣本值的波動(dòng)較大，存在偏差較大的樣本點(diǎn)，Cu在工業(yè)區(qū)含量高可能是由于某些特殊的樣本值；c工業(yè)區(qū)Cu含量偏度較大，表明Cu在工業(yè)區(qū)內(nèi)分布不均勻；d交通區(qū)Cu含量峰度比較大，表明交通區(qū)內(nèi)樣本值比較集中，Cu含量在交通區(qū)內(nèi)分布相對(duì)；e從評(píng)價(jià)指數(shù)看出，工業(yè)區(qū)的Cu含量受人為影響較

29、為嚴(yán)重；結(jié)論：綜合說(shuō)明Cu元素分布是不均勻的，推測(cè)造成污染主要原因是：在工業(yè)區(qū)某處存在銅冶煉廠等部門，產(chǎn)生大量Cu元素；交通區(qū)內(nèi)Cu含量也較多，分布較均勻，主要是由于汽車尾氣的排放。5. Hg元素Hg平均ng·g-1標(biāo)準(zhǔn)差ng·g-1變異系數(shù)偏度峰度評(píng)價(jià)193.0407102.90241.1062.550510.59162.6583051952642.35532244.0753.4935115.44389531.6770818.35300794337.9955818.628720.4902860.3847292.3856711.170173164446.82252180.2

30、74.8794986.47576243.5913212.766358185114.9917224.27981.9503994.84977126.918043.285477551表8 各功能區(qū)Hg統(tǒng)計(jì)量數(shù)據(jù)分析：a從平均值來(lái)看，Hg主要分布在工業(yè)區(qū)和交通區(qū)，生活區(qū)與公園綠地區(qū)含量也較多，山區(qū)Hg含量在背景值范圍內(nèi)；b工業(yè)區(qū)、交通區(qū)Hg含量變異系數(shù)較大，表明樣本值波動(dòng)較大，存在偏差較大的樣本點(diǎn)；生活區(qū)與公園綠地區(qū)變異系數(shù)較小，樣本值比較接近；c工業(yè)區(qū)與交通區(qū)Hg含量的偏度較大，表明Hg在工業(yè)區(qū)與交通區(qū)內(nèi)分布不均勻；d交通區(qū)的峰度較大，表明樣本相對(duì)集中，看出交通區(qū)內(nèi)Hg含量分布相對(duì)均勻；e從評(píng)價(jià)指數(shù)

31、看出，工業(yè)區(qū)與交通區(qū)Hg含量受人為影響很大；結(jié)論：綜合看出，城區(qū)內(nèi)Hg呈現(xiàn)不均勻分布，Hg污染現(xiàn)象較為明顯。特別是在工業(yè)區(qū)和交通區(qū)中污染的尤為嚴(yán)重，推測(cè)工業(yè)區(qū)內(nèi)某幾個(gè)區(qū)域建有關(guān)于電子、電池的加工廠，排放的廢物中含Hg含量很高；交通區(qū)、生活區(qū)由于海拔較低會(huì)較大程度的受到大氣程度的影響，并且生活區(qū)內(nèi)大量使用節(jié)能燈時(shí)也會(huì)產(chǎn)生汞污染。我們認(rèn)為公園內(nèi)會(huì)有較高的Hg污染可能是應(yīng)為采用了污水灌溉的方法來(lái)灌溉公園內(nèi)的植物，使得公園中土壤的Hg含量超標(biāo)。6. Ni元素Ni平均ug·g-1標(biāo)準(zhǔn)差ug·g-1變異系數(shù)偏度峰度評(píng)價(jià)118.34235.66230.30870.46513.11381.

32、491241685219.81178.37020.42250.61832.93811.610704607316.071412.059810.750393.05743714.104941.256405519417.617111.786980.6690648.72783992.80871.432284671515.289714.9742960.3253360.9989843.6511521.243066202表9 各功能區(qū)Ni統(tǒng)計(jì)量數(shù)據(jù)分析：a從平均值來(lái)看，Ni在各功能區(qū)分布比較均勻，且都在所給背景值的范圍之內(nèi)，表明該城區(qū)Ni含量基本上沒有受到人為因素的影響；b各區(qū)Ni含量變異系數(shù)較小，表明樣本值

33、比較接近；c交通區(qū)偏度相對(duì)較大，表明Ni在交通區(qū)含量和其他區(qū)相比，分布不均勻；d交通區(qū)峰度很大，表明樣本很集中，分布比較均勻；e從評(píng)價(jià)指數(shù)看出，該城區(qū)Ni元素含量受人為因素影響較??；結(jié)論：Ni元素在該城區(qū)分布基本是均勻的，在背景值范圍之內(nèi)，基本不受人為因素的影響；但交通運(yùn)輸中汽車尾氣的排放是Ni的主要來(lái)源，存在潛在的威脅。7. Pb元素Pb平均ug·g-1標(biāo)準(zhǔn)差ug·g-1變異系數(shù)偏度峰度評(píng)價(jià)169.106472.32521.04664.305823.63732.229237537293.040885.36730.91752.940511.48643.00131720433

34、8.8376720.663680.5320522.4276248.8629251.179222874463.534232.527260.5119651.3276284.9785092.049490416560.7085745.842980.7551322.1694867.2381081.958341014表 10 各功能區(qū)Pb統(tǒng)計(jì)量數(shù)據(jù)分析：a從平均值來(lái)看，Pb主要分布在工業(yè)區(qū)，在生活區(qū)、交通區(qū)和公園綠地均有分布，山區(qū)該元素含量在背景值范圍內(nèi)；b從變異系數(shù)來(lái)看，各功能區(qū)的值都比較小，表明樣本波動(dòng)小，樣本點(diǎn)比較接近；c從偏度考慮，生活區(qū)偏度較大，表明Pb在生活區(qū)和其他區(qū)比較，分布相對(duì)不均勻；d從

35、峰度來(lái)看，工業(yè)區(qū)、生活區(qū)峰度較大，表明樣本比較集中，工業(yè)區(qū)、生活區(qū)中Pb含量分布均勻；e從評(píng)價(jià)指數(shù)看出，Pb污染最嚴(yán)重的區(qū)域在工業(yè)區(qū)；結(jié)論：綜合看出，鉛呈現(xiàn)均勻分布，主要是工業(yè)污染，但生活中鉛污染無(wú)處不在。交通區(qū)中，汽車尾氣是主要來(lái)源；生活區(qū)中食物、學(xué)習(xí)用品玩具、電池、油漆等等都是鉛污染來(lái)源。8. Zn元素Zn平均ug·g-1標(biāo)準(zhǔn)差ug·g-1變異系數(shù)偏度峰度評(píng)價(jià)1237.0086443.63841.87185.166930.83663.4349077732277.9275350.83091.26232.822810.17934.027934783376.7553536.2

36、77850.4726432.2158469.2925161.0622353974242.8549384.78481.5844236.52307154.870653.5196355815154.2423230.91991.4971244.62133824.893812.235395445表11 各功能區(qū)Zn統(tǒng)計(jì)量數(shù)據(jù)分析：a從平均值來(lái)看，Zn元素主要分布在工業(yè)區(qū)、生活區(qū)和交通區(qū)，公園綠地含量相對(duì)較少，山區(qū)含量在背景值范圍之內(nèi)；b各功能區(qū)（除山區(qū)）變異系數(shù)都比較大，表明樣本波動(dòng)大，存在偏差較大的樣本點(diǎn)；c生活區(qū)和交通區(qū)偏度較大，相對(duì)于其他區(qū)，Zn的分布不均勻；d生活區(qū)和交通區(qū)峰度比較大，表明樣本值

37、比較集中，分布比較均勻；工業(yè)區(qū)中峰度相對(duì)較小，表明分布不均勻；e從評(píng)價(jià)指數(shù)看出，Zn污染最嚴(yán)重的區(qū)域在工業(yè)區(qū)；結(jié)論：重金屬元素Zn呈現(xiàn)不均勻分布，Zn在幾個(gè)地區(qū)會(huì)集中分布，推測(cè)主要來(lái)源于鋅礦開采、冶煉加工等工業(yè)部門的排放。交通區(qū)中汽車輪胎磨損也是Zn污染的重要來(lái)源。綜合上述分析，論文得到了8種重金屬元素在城市內(nèi)的分布情況及各污染元素的主要來(lái)源，如表12所示：重金屬元素污染主要原因As汽車尾氣（交通區(qū)）Cd工廠廢物（工業(yè)區(qū)）汽車尾氣（交通區(qū)）大氣沉降、污水灌溉（生活區(qū)與公園綠地區(qū)）Cr工業(yè)廢物大氣沉降Cu銅冶煉廠廢物（工業(yè)區(qū)）Hg工廠廢物（工業(yè)區(qū)）、大氣沉降Ni汽車尾氣（交通區(qū)）Pb鉛工業(yè)（工業(yè)

38、區(qū)）、汽車尾氣（交通區(qū)）日常生活（生活區(qū)）Zn工廠廢物（工業(yè)區(qū)）、輪胎磨損（交通區(qū)）表12 各污染元素的主要來(lái)源4.3模型的分析、建立與求解經(jīng)過(guò)任務(wù)二的分析與求解，我們知道了8種重金屬污染元素在各功能區(qū)的分布情況，其中As、Cd、Ni、Pb屬于均勻分布，Cr、Cu、Hg、Zn屬于不均勻分布。 但是，從整個(gè)城市來(lái)看，每一種重金屬污染物的分布都是非均勻的，通過(guò)問(wèn)題二可以看出，污染源主要集中在工業(yè)區(qū)、交通區(qū)和生活區(qū)。為了確定污染源的位置，必須清楚重金屬污染物的傳播特征，對(duì)319個(gè)采樣點(diǎn)的地理參數(shù)和各點(diǎn)的污染物濃度含量進(jìn)行插值，因此，建立合適的插值方法是解決本問(wèn)題的關(guān)鍵。圖10模型建立的流程圖4.3.

39、1 理論分析傳播特征由于重金屬元素主要來(lái)源于大氣中廢氣、粉塵顆粒；因此我們采用氣體擴(kuò)散模型11來(lái)得出重金屬元素污染的傳播特征方程，以污染源為原點(diǎn)建立新的平面直角坐標(biāo)系OXY。根據(jù)氣體擴(kuò)散模型我們引入時(shí)間變量t，結(jié)合元素濃度關(guān)于地理位置(x，y)的變化，建立重金屬污染濃度的平面矢量場(chǎng)。污染源的位置坐標(biāo)為（0，0）。記t時(shí)刻時(shí)，平面任一點(diǎn)重金屬元素的濃度為,根據(jù)假設(shè)設(shè)單位時(shí)間通過(guò)單位法向長(zhǎng)度的流量與濃度C成正比，有：是擴(kuò)散系數(shù)，grad表示梯度，負(fù)號(hào)表示由濃度高向濃度低的地方擴(kuò)散。先考察平面區(qū)域，其面積為S，包圍的曲線為l ，曲線l構(gòu)成圓形區(qū)域， l的外法線向量為，則在(t, t + Dt)內(nèi)通過(guò)

40、的流量為：內(nèi)污染物濃度的增量為：污染源在Dt 時(shí)間內(nèi)產(chǎn)生的濃度總量為:其中p0為常數(shù)，表示內(nèi)污染物濃度的單位面密度。根據(jù)質(zhì)量守恒定律和連續(xù)性原理，單位時(shí)間內(nèi)通過(guò)所選曲線l 的向外傳播的元素濃度與l 曲線內(nèi)濃度增量之和，等于污染源在單位時(shí)間產(chǎn)生的污染物濃度。有：。即：又根據(jù)曲線積分的Stokes公式：（div表示散度）所以：其中k為常數(shù)，表示重金屬元素污染的傳播速率。由以上兩式得：即為：根據(jù)參考文獻(xiàn)11解得：由于此題中重金屬元素濃度為固定值，傳播的濃度矢量場(chǎng)為穩(wěn)定的，不隨時(shí)間變化，因此對(duì)擴(kuò)散模型進(jìn)行改進(jìn)，令kt=m，m表示傳播的距離。取m為常數(shù)，則將擴(kuò)散模型改進(jìn)重金屬元素濃度傳播途徑的函數(shù)。即：

41、由此看出重金屬污染物的傳播隨地域呈指數(shù)變化。為了尋找到合理的插值方法，我們采用地統(tǒng)計(jì)學(xué)方法，結(jié)合ArcGIS軟件中的地統(tǒng)計(jì)分析擴(kuò)展模塊，分析土壤中8種重金屬的空間變異性和相關(guān)性。4.3.2 基于ArcGIS軟件的重金屬元素分布特征分析（1）以Cu為例，利用ArcGIS軟件對(duì)該元素進(jìn)行空間數(shù)據(jù)探索分析，利用直方圖顯示數(shù)據(jù)的概率分布特征。由圖中右上角數(shù)據(jù)可以看出，其偏度遠(yuǎn)大于0，峰度遠(yuǎn)大于3，可以說(shuō)明Cu在城市中的空間分布不是正態(tài)分布。圖11 Cu濃度直方圖 對(duì)Cu的采樣點(diǎn)濃度進(jìn)行對(duì)數(shù)轉(zhuǎn)換后得到下圖，可以看出中值接近均值，峰度有很大程度的下降且接近于3，近似于正態(tài)分布。說(shuō)明Cu的分布是指數(shù)形式的分

42、布。 圖12 對(duì)數(shù)變換后的Cu濃度直方圖(2)利用QQ-PLOT圖來(lái)檢驗(yàn)數(shù)據(jù)的正態(tài)分布情況。QQ-PLOT是根據(jù)變量分布的分位數(shù)對(duì)所指定的理論分布分位數(shù)繪制的圖形，是一種用來(lái)檢驗(yàn)樣點(diǎn)數(shù)據(jù)分布的統(tǒng)計(jì)圖。如果被檢驗(yàn)的樣點(diǎn)數(shù)符合所指定的分布，則代表樣點(diǎn)的點(diǎn)簇在一條直線上。直線表示正態(tài)分布，從下圖中可以看出散點(diǎn)與直線擬合地并不是很好，不符合正態(tài)分布。 圖13 Cu濃度QQ-PLOT圖對(duì)Cu的采樣點(diǎn)濃度進(jìn)行對(duì)數(shù)轉(zhuǎn)換后得到下圖，看到散點(diǎn)與直線擬合地非常好，說(shuō)明轉(zhuǎn)換后的Cu濃度值符合正態(tài)分布，也說(shuō)明Cu的分布是一種類似于指數(shù)形式的分布。 圖14 對(duì)數(shù)變換后的Cu濃度QQ-PLOT圖（3）空間趨勢(shì)反映了空間物

43、體在空間區(qū)域上變化的主題特征，它主要揭示了空間物體的總體規(guī)律，而忽略局部的變異。通過(guò)趨勢(shì)分析，可以將研究區(qū)平面圖上的點(diǎn)轉(zhuǎn)化為感興趣的屬性值為高度的三維視圖。將這些點(diǎn)按照兩個(gè)方向投影到與地圖平面正交的平面上，每個(gè)方向可用一個(gè)多項(xiàng)式來(lái)擬合。如果擬合曲線為平直的，說(shuō)明沒有全局趨勢(shì)；如果擬合的曲線為確定的曲線（二次函數(shù)曲線或指數(shù)函數(shù)曲線等），則存在某種全局趨勢(shì)。一旦存在全局趨勢(shì)，則可以通過(guò)確定性內(nèi)進(jìn)行插值。 圖15 Cu濃度的趨勢(shì)分析圖上圖中藍(lán)線與綠線基本平直，說(shuō)明Cu的濃度不具有全局趨勢(shì)，在插值時(shí)不考慮空間相關(guān)性。 （4）半變異函數(shù) 該圖可以反應(yīng)數(shù)據(jù)的空間相關(guān)程度，只有數(shù)據(jù)空間相關(guān)，才有必要進(jìn)行空間

44、插值法。圖表的橫坐標(biāo)表示任兩點(diǎn)的空間距離，縱標(biāo)表示該兩點(diǎn)的半變異函數(shù)值。根據(jù)距離越近越相似的原理，因而x值越小，y值應(yīng)該越小。圖16 Cu濃度的半變異函數(shù)圖由圖像看出，分布成明顯的分層現(xiàn)象.圖中分為3層，底層表示背景值，上面兩層表示污染源。由圖可得，Cu有兩個(gè)以上的污染源。其他7種元素的分析方法與Cu相同。4.3.3 插值方法的確定8種污染物都是用指數(shù)范普通克里格（Universal Cokriging）插值。使用此插值，可以確定函數(shù)以指數(shù)形式擬合。進(jìn)行分析時(shí)首先分析數(shù)據(jù)中的變化趨勢(shì)，獲得具體的擬合模型；其次對(duì)殘差數(shù)據(jù)（即原始數(shù)據(jù)減去趨勢(shì)數(shù)據(jù)）進(jìn)行克里格分析；最后，將趨勢(shì)面分析和殘差分析的克里

45、格結(jié)果加和，得到最終結(jié)果。在得出圖形后，可以根據(jù)指標(biāo)，驗(yàn)證差值好壞，從而修改模型，再進(jìn)行差值。所得圖形如任務(wù)一中的8種重金屬污染物分布圖像。最終進(jìn)行分析時(shí)，插值的具體參數(shù)如表13：污染物變換擬合多項(xiàng)式次數(shù)最小相關(guān)點(diǎn)數(shù)最大相關(guān)點(diǎn)數(shù)As對(duì)數(shù)355Cd對(duì)數(shù)345Cr對(duì)數(shù)355Cu對(duì)數(shù)3710Hg對(duì)數(shù)2310Ni對(duì)數(shù)355Pb對(duì)數(shù)2610Zn對(duì)數(shù)355表13 插值參數(shù)擬合后的誤差分析：平均值接近0；標(biāo)準(zhǔn)均方根接近1；，均方根??；且平均標(biāo)準(zhǔn)誤差接近均方根，則插值效果好。我們?cè)诓僮鬟^(guò)程中，主要考慮了品均值接近0和標(biāo)準(zhǔn)均方根接近1兩個(gè)指標(biāo)進(jìn)行差值。這兩個(gè)指標(biāo)也是衡量插值效果的主要指標(biāo)，具體得到誤差如表14：

46、污染物平均值標(biāo)準(zhǔn)均方根均方根平均標(biāo)準(zhǔn)誤差A(yù)s-0.181.373.323.43Cd-10.101.17238.036.40Cr-2.330.9867.623.84Cu-1.451.27169.672.962Hg3.85521.703325.69239.34Ni-0.281.0712.163.94Pb-0.650.8445.59188.54Zn1.611.64350.58339.71表14 擬合后的誤差分析從中發(fā)現(xiàn)，只有Cd的指標(biāo)相比其他污染物的效果較差，其余金屬元素分布的插值效果很好。4.3.4 污染源位置的確定從ArcGIS中導(dǎo)出插值后的數(shù)據(jù)到Matlab。通過(guò)Matlab畫出各個(gè)金屬元素分

47、布圖。編寫程序（見附表）畫出其對(duì)應(yīng)的梯度矢量圖,如圖17：圖17 梯度矢量圖通過(guò)矢量指向，找到金屬含量最高、梯度為零的點(diǎn)即為污染源。一種金屬污染物，可能會(huì)有多個(gè)主要污染源，其坐標(biāo)分布如表15：污染物As Cd Cr Cu x(m)y(m)x(m)y(m)x(m)y(m)x(m)y(m)坐標(biāo)4270720817652339337758433637599820323149173542314972853492239352612625282656401023587298444污染物Hg Ni Pb Zn x(m)y(m)x(m)y(m)x(m)y(m)x(m)y(m)坐標(biāo)243931023537622

48、2160529961072311913447232858841084044345293393754782240810366283071147250349952129182715153949517206424161522698156832729544617545表15 各重金屬污染源坐標(biāo)各重金屬污染源示意圖（圖中深色小點(diǎn)表示污染源位置）：圖18 各重金屬污染源示意圖4.4 任務(wù)4的分析與求解模型的優(yōu)點(diǎn)：1. 模型中運(yùn)用ArcGIS軟件中地理處理與分析模塊能很好的做分析重金屬元素的空間分布2. 模型中運(yùn)用了源的概念，依據(jù)高斯公式，理論分析出金屬元素傳播擴(kuò)散滿足指數(shù)分布；3. 使用了符合實(shí)際的指數(shù)范

49、普通克里格插值，比較符合實(shí)際。模型的缺點(diǎn)：1. 模型中沒有考慮海拔對(duì)各重金屬元素分布的影響；2. 模型中沒有考慮各重金屬元素之間的相關(guān)性；3. 模型中，相比其他金屬元素，Cd的插值效果不好，缺少數(shù)據(jù)。沒有進(jìn)一步分析。模型的推廣：為了更好地研究城市地質(zhì)環(huán)境的演變模式，還應(yīng)收集：1. 不同時(shí)間內(nèi)的重金屬元素含量分布，2. 更多關(guān)于該地區(qū)的地質(zhì)、天氣等方面的數(shù)據(jù)。模型的改進(jìn)：通過(guò)1的數(shù)據(jù)，建立元素分布圖關(guān)于時(shí)間的演變模型：建立元素濃度關(guān)于時(shí)間變化的矢量場(chǎng)，得到；利用高斯公式：，得到污染源的位置和其單位時(shí)間內(nèi)的排污量。使用已經(jīng)得到的污染源數(shù)據(jù)，結(jié)合已經(jīng)建立的演變模型，預(yù)測(cè)地質(zhì)環(huán)境的演變模式。通過(guò)2的數(shù)

50、據(jù)，更全面的分析污染傳播特性，建立更準(zhǔn)確的傳播模型。參考文獻(xiàn)1 劉哲民，寶雞土壤重金屬污染及其防治J.干旱區(qū)資源與環(huán)境， 19（2）：101-104，2005。2 李靜，謝正苗，徐建明，等，杭州市郊蔬菜地土壤重金屬環(huán)境質(zhì)量評(píng)價(jià)J.生態(tài)環(huán)境， 12（3）：277-280，2003 。3 潘劍軍主編，土壤資源調(diào)查與評(píng)價(jià)M，北京：中國(guó)農(nóng)業(yè)出版社，2004。4張輝，南京地區(qū)土壤沉積物中重金屬形態(tài)研究，環(huán)境科學(xué)報(bào)，17（3）：346-351,1997。5Inacio M M,Pereira V,Pinto M S.Mercury contamination in sandy soils surround

51、ings an industrial emission source (Estarreja,Portugal)J.Geoderma ,85(4):325-339,1998.6Pichtel J,Sawyerr H T,Czamowska K.Spatial and temporal distribution of metals in soils in Warsaw,Poland J.Environmental Pollution,98(2):164-174,1997.7曹曉東，孫波，宗良綱，等，丘陵紅壤重金屬?gòu)?fù)合污染的空間變異分析J,土壤，37（2）：140-146,2005。8Brit L S,Tom A, Eirik F.et al.Heavy Metal Survey