焦作市城區(qū)土壤重金屬空間變異特征及污染評(píng)價(jià)

李艷利，肖春燕，王守全，劉坤朋

（1.河南理工大學(xué) 資源環(huán)境學(xué)院，河南 焦作454000；2.河海大學(xué) 環(huán)境學(xué)院，江蘇 南京210098；3.蘇州科技學(xué)院 環(huán)境與科學(xué)工程學(xué)院，江蘇 蘇州215009）

城市土壤是城市生態(tài)系統(tǒng)的重要組成部分，其環(huán)境質(zhì)量狀況直接關(guān)系到人類的健康安全［1－2］。伴隨著工業(yè)化和城市化的發(fā)展，來(lái)自工業(yè)排放、大氣沉降、生活垃圾和汽車尾氣排放等方式造成的重金屬在土壤中累積。且城市土壤重金屬污染問(wèn)題也引起了國(guó)內(nèi)外學(xué)者的關(guān)注［3－10］，東南沿海經(jīng)濟(jì)發(fā)達(dá)城市和中西部等大城市紛紛開展了對(duì)土壤重金屬的研究［11－17］。目前國(guó)內(nèi)關(guān)于土壤重金屬的研究主要集中在兩個(gè)方面，一是對(duì)土壤重金屬的污染分布及來(lái)源進(jìn)行分析［18－21］，另一方面是對(duì)土壤重金屬進(jìn)行污染評(píng)價(jià)和修復(fù)［14－15］。土地利用方式、人類干擾方式和距離污染源遠(yuǎn)近等的不同均會(huì)使城市土壤重金屬含量存在顯著的空間差異。

本文采用GIS空間分析技術(shù)和SPSS技術(shù)，與綜合污染指數(shù)評(píng)價(jià)法相結(jié)合的方法，在焦作市開展表層土壤重金屬污染程度、空間變異特征以及分布規(guī)律的研究。為焦作市土壤重金屬防治提供一定理論基礎(chǔ)。

1 數(shù)據(jù)來(lái)源與研究方法

1.1 研究區(qū)概況

焦作市北依太行與山西省接壤，東臨新鄉(xiāng)市，西臨濟(jì)源市，南臨黃河與鄭州、洛陽(yáng)相望，地跨北緯35°10′—35°21′，東 經(jīng) 113°4′—113°26′，東 西 長(zhǎng) 約32.5km，南北寬約19.7km。焦作市城區(qū)地轄解放區(qū)、中站區(qū)、馬村區(qū)、山陽(yáng)區(qū)4個(gè)區(qū)，面積4.24×105km2，市區(qū)人口約88萬(wàn)人。焦作市屬暖溫帶大陸性季風(fēng)型氣候，春旱多風(fēng)、夏熱多雨、秋高氣爽、冬寒少雪。平原地區(qū)年平均氣溫14.9℃，山區(qū)比平原低2～3℃，年均降雨量603.5mm，無(wú)霜期231d，年日照時(shí)數(shù)2 422.7h。

1.2 研究方法

1.2.1 樣品采集與測(cè)定 2012年4月在焦作市3個(gè)城區(qū)（解放區(qū)、山陽(yáng)區(qū)、中站區(qū)）采集0—20cm表層土壤。按照網(wǎng)格法（2km×2km）均勻布點(diǎn)，采樣點(diǎn)的選取兼顧了各種土地類型，具有空間代表性、經(jīng)濟(jì)性和可操作性，使用GPS定位，并準(zhǔn)確記錄經(jīng)緯度。樣品采用五分法采集，混合均勻后留下1kg土樣裝入塑料袋，共采集了44個(gè)樣品，其中山陽(yáng)區(qū)23個(gè)，解放區(qū)12個(gè)，中站區(qū)9個(gè)。將采集的樣品放置于通風(fēng)處在室溫下自然風(fēng)干，避免陽(yáng)光直接曝曬和其他污染。用研缽反復(fù)研磨，過(guò)20目篩子得到粗樣，再次研磨，過(guò)100目篩子得到細(xì)樣，然后低溫保存、備用。所有樣品經(jīng)標(biāo)準(zhǔn)方法進(jìn)行消化處理后，Cr，Cu的測(cè)定采用火焰原子吸收分光光度法，Cd，Pb的測(cè)定采用石墨爐原子吸收分光光度法進(jìn)行測(cè)定。

1.2.2 土壤污染評(píng)價(jià)方法 土壤污染評(píng)價(jià)采用單因子污染指數(shù)法和內(nèi)梅羅綜合指數(shù)法。單因子污染指數(shù)法適合單一因子污染的評(píng)價(jià)，內(nèi)梅羅綜合指數(shù)法能全面、綜合地反映土壤的污染程度。

（1）單因子指數(shù)法：

式中：Pi——土壤中污染物i的環(huán)境質(zhì)量指數(shù)；ci——污染物i的實(shí)測(cè)濃度（mg／kg）；Si——污染物i的評(píng)價(jià)標(biāo)準(zhǔn)（mg／kg），采用河南省土壤重金屬環(huán)境背景值［22］。下同。

（2）內(nèi)梅羅綜合指數(shù)法：

式中：P綜——土壤中重金屬元素的綜合污染指數(shù)；n——參與評(píng)價(jià)的重金屬的種類總數(shù)。根據(jù)P值變幅，結(jié)合作物受害程度和污染物積累狀況劃分土壤質(zhì)量分級(jí)：Pi≤1為非污染，1＜Pi≤2為輕度污染，2＜Pi≤3為中度污染，Pi≥3為重度污染；PN≤0.7為安全等級(jí)，0.7＜PN≤1為警戒線，1＜PN≤2為輕度污染，2＜PN≤3為中度污染，PN≥3為重度污染。

1.2.3 數(shù)據(jù)分析方法 數(shù)據(jù)分析采用SPSS 16.0軟件完成，半變異函數(shù)擬合、樣點(diǎn)重金屬濃度分布圖以及重金屬的插值圖用ArcGIS 9.3軟件完成。

2 結(jié)果與分析

2.1 土壤重金屬含量的統(tǒng)計(jì)分析

土壤表層樣品重金屬含量測(cè)定結(jié)果如表1所示。土壤中Cr，Cu，Cd和Pb的濃度范圍分別為：133.13～783.13，20.88～72.13，0.007～1.773和7.45～45.49mg／kg，土壤樣品中 Cr，Cu，Cd和 Pb這4種重金屬平均含量分別為：378.86，36.26，0.33和20.23mg／kg。變異系數(shù)反映了所有采樣點(diǎn)的平均變異程度，研究區(qū)域Cr，Cu，Cd和Pb的變異系數(shù)分別為：63.23%，28.21%，107.21%和32.97%，這4種重金屬元素的平均變異程度由大到小依次為：Cd，Cr，Pb和Cu。說(shuō)明Cr和Cd受人類活動(dòng)的強(qiáng)烈干預(yù)，Cu和Pb受人類干擾相對(duì)較小。相應(yīng)研究采用單樣本的Kolmogorov—Smirnov檢驗(yàn)，Cd和Pb符合正態(tài)分布，Cr和Cu表現(xiàn)出偏態(tài)分布特征。按照地統(tǒng)計(jì)學(xué)空間變異分析的要求，對(duì)Cr和Cu進(jìn)行對(duì)數(shù)轉(zhuǎn)換，使其符合正態(tài)分布，用于空間變異分析。

2.2 土壤重金屬污染狀況

對(duì)所有土壤樣品重金屬分析數(shù)據(jù)取其平均值并通過(guò)單因子污染指數(shù)進(jìn)行評(píng)價(jià)，計(jì)算結(jié)果如表2所示。Cr污染指數(shù)大于1，其他元素均小于1。Cr，Cu，Cd和 Pb的污染指數(shù)分別為：1.52，0.18，0.56和0.06。土壤中Cr和Cd的最大值分別為土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)的3.13和2.96倍。而Cu和Pb的最大值分別為環(huán)境質(zhì)量標(biāo)準(zhǔn)的0.36和0.13倍。這也說(shuō)明Cr和Cd在空間分布上具有較為強(qiáng)烈的空間變異性。

表1 焦作市城區(qū)表層土壤重金屬含量的統(tǒng)計(jì)特征（n＝44）

表2 焦作市城區(qū)土壤重金屬單因子污染指數(shù)

44個(gè)采樣點(diǎn)中，Cr的污染指數(shù)超過(guò)1的有21個(gè)，超標(biāo)率為47.73%，其中大多數(shù)為中度污染。Cr污染較重的區(qū)域主要集中在山陽(yáng)區(qū)。土壤中Cu的單項(xiàng)污染指數(shù)均小于1，說(shuō)明焦作市土壤中Cu濃度未超過(guò)國(guó)家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)Ⅱ級(jí)標(biāo)準(zhǔn)值。土壤中Cd的污染指數(shù)超過(guò)1的有10個(gè)，超標(biāo)率為22.73%，焦作市污水處理廠附近、焦作師范?？茖W(xué)校、萬(wàn)方小區(qū)等采樣點(diǎn)的土壤單因子污染指數(shù)均大于1，說(shuō)明這些地區(qū)已受到不同程度的輕度污染。土壤Pb的污染指數(shù)均小于1，濃度值接近河南省土壤背景值22.3mg／kg，表明焦作市土壤中Pb污染程度很低。

焦作市城區(qū)土壤重金屬的綜合污染指數(shù)超過(guò)1的有23個(gè)，占52.7%，表明焦作市城區(qū)土壤重金屬污染整體較為嚴(yán)重。其中44個(gè)樣點(diǎn)中6.82%的樣點(diǎn)綜合污染指數(shù)在1～2，屬于輕度污染。40.18%的樣點(diǎn)綜合污染指數(shù)在2～3，屬于中度污染。2.27%的樣點(diǎn)綜合污染指數(shù)≥3，屬于重度污染。說(shuō)明不同地區(qū)污染積累程度差異明顯。

2.3 焦作市不同城區(qū)土壤重金屬分布特征及累積狀況

2.3.1 焦作市不同城區(qū)土壤重金屬空間分布特征焦作市3個(gè)城區(qū)土壤重金屬空間分布如表3所示。元素Cr和Cd的最大值均出現(xiàn)在山陽(yáng)區(qū)，元素Pb和Cu的最大值均出現(xiàn)在解放區(qū)，Pb的最小值出現(xiàn)在山陽(yáng)區(qū)，Cu的最小值出現(xiàn)在中站區(qū)。

Cr和Pb均在解放區(qū)呈最大變異程度，說(shuō)明解放區(qū)Cr和Pb受人類活動(dòng)的干擾程度最強(qiáng)。Cd在山陽(yáng)區(qū)和解放區(qū)均呈極強(qiáng)變異程度，這說(shuō)明Cd在山陽(yáng)區(qū)和解放區(qū)受人類干擾程度均很強(qiáng)。而Cu在3個(gè)城區(qū)的變異程度均較低，說(shuō)明Cu受人類活動(dòng)干擾程度相對(duì)較低。

表3 焦作市不同城區(qū)土壤重金屬含量描述性統(tǒng)計(jì)分析

解放區(qū)土壤中Pb和Cu的平均值在3個(gè)城區(qū)中最高。其中Pb（20.89mg／kg）分別是山陽(yáng)區(qū)的1.04和1.07倍。Cu（37.88mg／kg）分別是山陽(yáng)區(qū)的1.06和1.07倍。Pb低于環(huán)境背景值，其他元素均超過(guò)環(huán)境背景值，山陽(yáng)區(qū)土壤Cr含量是背景值的7.75倍，解放區(qū)土壤Cr含量是背景值的3.59倍，中站區(qū)土壤Cr含量是背景值的4.71倍。山陽(yáng)區(qū)土壤Cd含量是背景值的6.91倍，解放區(qū)土壤Cd含量是背景值的2.91倍，中站區(qū)土壤Cd含量是背景值的3.25倍。3個(gè)城區(qū)Cr和Cd的分布特征總體上存在山陽(yáng)區(qū)＞中站區(qū)＞解放區(qū)的趨勢(shì)，Pb和Cu的分布特征總體上存在解放區(qū)＞山陽(yáng)區(qū)＞中站區(qū)的趨勢(shì)。

2.3.2 焦作市不同城區(qū)土壤重金屬累積狀況 表4是不同城區(qū)土壤的重金屬污染指數(shù)。其中山陽(yáng)區(qū)、解放區(qū)和中站區(qū)土壤Cd，Cu和Pb的污染指數(shù)均小于1，說(shuō)明焦作市城區(qū)土壤Cu，Cd和Pb處于非污染狀態(tài)。解放區(qū)Cr污染指數(shù)小于1，說(shuō)明解放區(qū)Cr處于非污染狀態(tài)，山陽(yáng)區(qū)和中站區(qū)Cr污染指數(shù)均大于1，達(dá)到污染等級(jí)，屬于輕污染狀態(tài)。從4種重金屬元素綜合污染指數(shù)來(lái)看，山陽(yáng)區(qū)土壤重金屬達(dá)到中度污染，中站區(qū)土壤重金屬屬輕度污染，解放區(qū)的處于尚清潔狀態(tài)。

表4 焦作市不同城區(qū)土壤重金屬污染指數(shù)

山陽(yáng)區(qū)工業(yè)相對(duì)較集中，商業(yè)發(fā)達(dá)以及公路網(wǎng)密集、車流量大等，導(dǎo)致大量重金屬的排放，累積于土壤中，這可能是造成其中度污染的原因。解放區(qū)建區(qū)較早，但其屬于尚清潔狀態(tài)，主要?dú)w功于其工業(yè)相對(duì)較少，僅電廠和火車站周圍的采樣點(diǎn)污染指數(shù)大于2。中站區(qū)農(nóng)田以土壤中重金屬污染指數(shù)均大于2，可能是造成中站區(qū)達(dá)到輕度污染的原因，且該區(qū)分布著礦務(wù)局電廠、多氟多化工廠等，也可能導(dǎo)致附近土壤中重金屬的累積。

2.4 焦作市城區(qū)土壤重金屬來(lái)源分析

統(tǒng)計(jì)分析中的相關(guān)分析可有效揭示土壤重金屬之間的關(guān)系和來(lái)源。相關(guān)性分析表明（表5），土壤中Cd與Cr，Pb具有顯著的空間相關(guān)性，Pb和Cu具有極顯著的相關(guān)性，表明土壤重金屬之間表現(xiàn)為同源性或復(fù)合污染。Pb，Cu主要來(lái)源于汽車輪胎磨損及油的泄漏［23－24］，因此Pb和Cu可能主要來(lái)源于交通排放污染物和鋼鐵廠等。一般認(rèn)為Cr主要來(lái)源于成土母質(zhì)［25］，但Cr和Cd也有相同來(lái)源，因此可推斷Cr和Cd具有復(fù)合污染的特征，其來(lái)源包括交通排放污染物、工業(yè)和成土母質(zhì)等。

表5 焦作市城市土壤重金屬的相關(guān)分析

2.5 土壤重金屬污染空間變異結(jié)構(gòu)分析

應(yīng)用球形模型對(duì)重金屬元素的半變異函數(shù)進(jìn)行了擬合，表6是各參數(shù)半變異函數(shù)的相應(yīng)參數(shù)。塊金值表示變量在小尺度上的非連續(xù)變異，主要來(lái)源于抽樣尺度的隨機(jī)變異和測(cè)量誤差，反映了隨機(jī)效應(yīng)的程度。基底效應(yīng)揭示了區(qū)域化變量的空間相關(guān)程度，如果其值小于25%，說(shuō)明系統(tǒng)具有強(qiáng)烈的空間相關(guān)性；在25%～75%表明系統(tǒng)具有中等程度的空間相關(guān)性；大于75%表明系統(tǒng)空間相關(guān)性很弱。

表6 焦作市城區(qū)土壤重金屬半方差函數(shù)理論模型及相關(guān)參數(shù)

Cu和Cd的基底效應(yīng)值分別為100%和85.41%，空間相關(guān)性很弱，以隨機(jī)變異為主，主要受外源污染（交通運(yùn)輸、電廠以及管理措施等人為活動(dòng)）的影響，反映了人類活動(dòng)強(qiáng)烈影響著土壤Cu和Cd的空間變異性。Cu和Cd的空間相關(guān)性的最大距離分別為15.05和17.21km。Pb的基底效應(yīng)值為56.98%，呈中等程度的空間相關(guān)性，Cr的基底效應(yīng)值為1.4%，具有強(qiáng)烈的空間相關(guān)性，表明Cr主要受內(nèi)源污染（如成土母質(zhì)、氣候變化、土壤特性等）的控制。而Pb受到內(nèi)源和外源因素的綜合影響。Pb和Cr空間自相關(guān)的最大平均距離分別為6.33和2.62km。

2.6 土壤重金屬的空間分布分布格局

由圖1可看出，焦作市城市土壤Cr濃度普遍較高，由南到北，濃度逐漸遞減，東西方向上由中間向兩邊逐漸遞增，呈帶狀分布。Cu元素含量由東南往西北逐漸下降，高濃度區(qū)主要集中在高新區(qū)西南部、解放區(qū)中部以及山陽(yáng)區(qū)偏東北部，低含量區(qū)主要位于焦作市西部。Cu污染高值區(qū)集中的高新區(qū)西南部是農(nóng)業(yè)區(qū)，土壤中重金屬含量易受農(nóng)田施肥的影響，且由實(shí)際調(diào)查知道該地區(qū)鄉(xiāng)鎮(zhèn)企業(yè)較多，農(nóng)田多采用企業(yè)污水灌溉，也對(duì)土壤中Cu含量做出一定的貢獻(xiàn)。位于解放區(qū)中部的焦作電廠也是高值區(qū)域，電廠會(huì)產(chǎn)生大量的粉煤灰、飄塵，大氣沉降亦造成主導(dǎo)風(fēng)向（東北方向）上重金屬含量的高值，對(duì)土壤重金屬污染也起到一定的積累作用。Cd元素含量從東向西逐漸降低，且中北部區(qū)域濃度最低。高值區(qū)域主要集中在人口比較集中、交通流量高的區(qū)域，例如焦作師范?？茖W(xué)校、人民廣場(chǎng)、萬(wàn)方小區(qū)、火車站、中州鋁廠家屬院等區(qū)域附近，因此可能交通污染、商業(yè)活動(dòng)對(duì)土壤中Cd的積累起到一定的作用。農(nóng)業(yè)活動(dòng)如農(nóng)田施肥及污水灌溉也可能加劇重金屬污染。同時(shí)，舊城區(qū)改造前該區(qū)域建有許多化工廠、冶煉廠等，雖然均已拆除，但重金屬帶來(lái)的威脅可能仍然存在。Pb元素含量普遍較低，解放區(qū)原化工三廠附近濃度最高，呈環(huán)狀向南逐漸降低，不同區(qū)域差異性較低，Pb的變異程度較低上也從側(cè)面證明了這點(diǎn)。同時(shí)考慮到焦作市季風(fēng)風(fēng)向?yàn)闁|北風(fēng)，解放區(qū)擁有容量為1.2×106kW的焦作電廠、焦西礦、原化工三廠這些重污染企業(yè)，化工三廠的位置恰好位于Pb濃度最高值點(diǎn)，可以初步推斷，Pb的污染源可能是位于解放區(qū)北部的工礦企業(yè)，有較大可能性為大氣污染擴(kuò)散點(diǎn)源。因此Pb濃度分布受土壤屬性和工業(yè)污染共同影響。

圖1 焦作市城區(qū)土壤重金屬含量空間分布

3 討論

研究結(jié)果顯示，相較于國(guó)內(nèi)外其他城市Cr，Pb，Cd和Cu的污染狀況，焦作城區(qū)土壤Cr平均含量比其他城市普遍都高，而Pb含量普遍較其他城市低。Cr平均含量較高可能是因?yàn)榻棺鞒菂^(qū)土壤中Cr的基底值較高的原故；Pb含量較低的原因可能是因?yàn)槠囕喬ツp和油的泄露是Pb的一個(gè)重要來(lái)源。焦作市作為中小城市，交通運(yùn)輸量相對(duì)較小因而Pb含量普通較低。分析中發(fā)現(xiàn)，焦作城區(qū)總體樣本中各采樣點(diǎn)Cr具有強(qiáng)烈的空間相關(guān)性，且變異程度較高，也充分說(shuō)明了Cr的分布和污染程度差異較大，其主要受內(nèi)源污染（成土母質(zhì)、成土過(guò)程、土壤特性等）的影響。而焦作城區(qū)Cu和Cd的空間相關(guān)性很弱，主要以隨機(jī)變異為主，說(shuō)明Cu和Cd受外源污染的影響。焦作城區(qū)總體樣本中Cd的分布和污染程度差異極大，說(shuō)明其主要受人類活動(dòng)的影響，其高值區(qū)域主要集中在人口比較集中、交通流量高的區(qū)域，可能商業(yè)活動(dòng)和交通流量對(duì)其污染分布貢獻(xiàn)較大。Cu的高值區(qū)主要集中在工農(nóng)業(yè)企業(yè)集中的區(qū)域，因?yàn)槠滢r(nóng)田多采用污水灌溉，電廠產(chǎn)生大量的粉煤灰、飄塵，也對(duì)Cu的積累做出了一定貢獻(xiàn)。城市Pb呈中等程度的空間相關(guān)性，且總體樣本中所有采樣點(diǎn)Pb的分布和污染程度差異也較小。土壤Cr，Cd和Pb之間以及Cu和Pb之間呈顯著相關(guān)，說(shuō)明重金屬污染具有復(fù)合污染和同源污染的特征。其中外界重金屬來(lái)源，包括大氣沉降、城市垃圾，污水灌溉，工業(yè)廢渣和農(nóng)用化學(xué)品等對(duì)焦作城區(qū)土壤重金屬的含量有一定貢獻(xiàn)。

4 結(jié)論

（1）研究區(qū)域內(nèi)4種土壤重金屬元素含量的統(tǒng)計(jì)結(jié)果表明，Cd的變異系數(shù)最大，超過(guò)了100%，變異程度最強(qiáng)，Cr，Pb和Cu均呈中等變異性。

（2）三城區(qū)Cr和Cd的分布特征總體上存在山陽(yáng)區(qū)＞中站區(qū)＞解放區(qū)的趨勢(shì)，Pb和Cu的分布特征總體上存在解放區(qū)＞山陽(yáng)區(qū)＞中站區(qū)的趨勢(shì)。

（3）土壤污染狀況分析表明，土壤Cr和Cd污染相對(duì)嚴(yán)重，樣點(diǎn)超標(biāo)率在47.73%和22.73%，土壤中Cd與Cr，Pb具有顯著的空間相關(guān)性，Pb和Cu具有極顯著的相關(guān)性，土壤重金屬之間表現(xiàn)為同源性或者復(fù)合污染。土壤重金屬污染程度反映了商業(yè)區(qū)分布、廢棄物排放、工業(yè)區(qū)分布、城市交通等對(duì)城市土壤重金屬分布影響的復(fù)雜性。

（4）半變異函數(shù)模型擬合結(jié)果表明，Cu和Cd存在較弱的空間自相關(guān)性，說(shuō)明這兩種元素主要以隨機(jī)變異為主，受人類活動(dòng)影響和干擾大。Cr存在極強(qiáng)的空間自相關(guān)性，說(shuō)明其主要受內(nèi)在因素影響；Pb存在中等程度的空間自相關(guān)性，說(shuō)明其濃度受內(nèi)在因素和外在因素綜合影響。

（5）采用普通克里格最優(yōu)內(nèi)插法得到了焦作市表層土壤重金屬的含量空間分布圖。Cr的空間分布呈帶狀，由東南向西北遞減；Cu的空間分布由中心向四周輻射狀遞減；Cd的空間分布呈帶狀分布，由東部向西部遞減；Pb的空間分布呈環(huán)狀分布，由北部向南部遞減。

［1］ Miguel D E，Jimenez D G M，Llamas J F.The overlooked contribution of compost application to the trace element load in the urban soil of Madrid（Spain）［J］.The Science of the Total Environment，1998，215（1／2）：113－122.

［2］ 陳燕芳.我國(guó)城市土壤重金屬污染及其治理研究進(jìn)展綜述［J］.中國(guó)人口資源與環(huán)境，2011，21（3）：536－539.

［3］ Kimple C R，Morel J L.Urban soil management：a growing concern［J］.Soil Science，2001，65（1）：31－40.

［4］ Turer D，Maynard J B，Sansalone J J.Heavy Metal Contamination in the Soils of Urban Highways：Comparison Between Runoff and Soil Concentrations at Cincinnati，Ohio［J］.Water，Air，and Soil Pollution，2001，132（3／4）：293－314.

［5］ Madrid L，Diaz Barrientos E.Madrid F.Distribution of heavy metal contents of urban soils in parks of Seville［J］.Chemosphere，2002，49（10）：1301－1308.

［6］ Jennings A A，Petersen E J.Variability of North American regulatory guidance for heavy metal contamination of residential soil［J］.Environmental Engineering and Science，2006，5（6）：485－508.

［7］ Imperato M，Adamo P，Naimo D，et al.Spatial distribution of heavy metals in urban soils of Naples City（Italy）［J］.Environmental Pollution，2003，124（2）：247－256.

［8］ Jordao C，Nascentes C，Cecon P，et al.Heavy Metal Availability in Soil Amended with Composted Urban Solid Wastes［J］.Environmental Monitoring and As－sessment，2006，112（1／3）：309－326.

［9］ Linde M，Bengtsson H，Oborn I.Concentrations and Pools of Heavy Metals in Urban Soils in Stockholm，Sweden［J］.Water，Air and Soil Pollution，2001，1（3／4）：83－101.

［10］ 柳云龍，章立佳，韓曉菲，等.上海城市樣帶土壤重金屬空間變異特征及污染評(píng)價(jià)［J］.環(huán)境科學(xué)，2012，33（2）：599－605.

［11］ 王濟(jì)，張浩，曾希柏，等.貴陽(yáng)市城區(qū)土壤重金屬分布特征及污染評(píng)價(jià)［J］.土壤，2010，42（6）：928－934.

［12］ 馬建華，李燦，陳云增.土地利用與經(jīng)濟(jì)增長(zhǎng)對(duì)城市土壤重金屬污染的影響：以開封市為例［J］.土壤學(xué)報(bào)，2011，48（4）：743－750.

［13］ 錢翌，張瑋，冉德超.青島城市土壤重金屬的形態(tài)分布及影響因素分析［J］.環(huán)境化學(xué)，2011，30（3）：652－657.

［14］ 黃敏，楊海舟，余萃，等.武漢市土壤重金屬積累特征及其污染評(píng)價(jià)［J］.水土保持學(xué)報(bào)，2010，24（4）：135－139.

［15］ 劉勇，岳玲玲，李晉昌.太原市土壤重金屬污染及其潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)［J］.環(huán)境科學(xué)學(xué)報(bào)，2011，31（6）：1285－1293.

［16］ 鐘曉蘭，周生路，李江濤，等.長(zhǎng)江三角洲地區(qū)土壤重金屬污染的空間變異特征：以江蘇省太倉(cāng)市為例［J］.土壤學(xué)報(bào)，2007，44（1）：33－39.

［17］ 趙淑蘋，陳立新.大慶地區(qū)不同土地利用類型土壤重金屬分析及生態(tài)危害評(píng)價(jià)［J］.水土保持學(xué)報(bào)，2011，25（5）：195－199.

［18］ 謝小進(jìn)，康建成，李衛(wèi)江，等.上海寶山區(qū)農(nóng)用土壤重金屬分布與來(lái)源分析［J］.環(huán)境科學(xué)，2010，31（3）：768－774.

［19］ 李瑞平，郝英華，李光德，等.泰安市農(nóng)田土壤重金屬污染特征及來(lái)源解析［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2011，30（10）：2012－2017.

［20］ 段雪梅，蔡煥興，巢文軍.南京市表層土壤重金屬污染特征及污染來(lái)源［J］.環(huán)境科學(xué)與管理，2010，35（4）：31－35.

［21］ 陳景輝，盧新衛(wèi)，翟萌.西安城市路邊土壤重金屬來(lái)源與潛在風(fēng)險(xiǎn)［J］.應(yīng)用生態(tài)學(xué)報(bào)，2011，22（7）：1810－1816.

［22］ 邵豐收，周皓韻.河南省主要元素的土壤環(huán)境背景值［J］.能源保護(hù)，2005，20（2）：28－29.

［23］ 張一修，王濟(jì)，秦樊鑫，等.貴陽(yáng)市道路灰塵和土壤重金屬來(lái)源識(shí)別比較［J］.環(huán)境科學(xué)學(xué)報(bào)，2012，32（1）：204－212.

［24］ 曹益金，李曉燕，王干珍，等.貴陽(yáng)市城區(qū)土壤重金屬累積現(xiàn)狀研究［J］.土壤通報(bào)，2012，43（2）：484－489.

［25］ Boruvka L，Vacek O，Jehlicka J.Principal component analysis as a tool to indicate the origin of potentially toxic elements in soils［J］.Geoderma，2005，128（3／4）：289－300.