污染場(chǎng)地土壤重金屬污染空間特征分析
——以某搬遷電鍍廠為例

耿治鵬，宋頡，王春林，蔣衛(wèi)國(guó)*，陳征

1.遙感科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京師范大學(xué)地理科學(xué)學(xué)部

2.環(huán)境遙感與數(shù)字城市北京市重點(diǎn)實(shí)驗(yàn)室, 北京師范大學(xué)地理科學(xué)學(xué)部

3.生態(tài)環(huán)境部土壤與農(nóng)業(yè)農(nóng)村生態(tài)環(huán)境監(jiān)管技術(shù)中心

隨著我國(guó)環(huán)保事業(yè)的發(fā)展，位于城市內(nèi)部的污染企業(yè)相繼關(guān)停、搬遷，但這類企業(yè)在生產(chǎn)過(guò)程中遺留在場(chǎng)地土壤中的污染并未消失，仍然存在著較大的風(fēng)險(xiǎn)隱患[1-3]。自然環(huán)境本身存在的重金屬，一般為穩(wěn)定的化合態(tài)，對(duì)周邊生物和環(huán)境的影響較小，而工業(yè)生產(chǎn)會(huì)導(dǎo)致周邊土壤游離態(tài)重金屬濃度增加，易對(duì)周邊生物和環(huán)境產(chǎn)生較大影響[4]。重金屬污染造成的危害具有滯后性和隱蔽性[5]。我國(guó)人口基數(shù)大，工業(yè)發(fā)展較快，工業(yè)生產(chǎn)產(chǎn)生的固體廢物的堆放和傾倒、有害廢水的滲透、有害氣體及飄塵隨雨水降落等過(guò)程導(dǎo)致或加劇了土壤重金屬污染[6]。原環(huán)境保護(hù)部與原國(guó)土資源部2014年公布的《全國(guó)土壤污染狀況調(diào)查公報(bào)》顯示，我國(guó)受Ni、Cr、Pb等重金屬污染的耕地面積約占總耕地面積的1/5，耕地土壤的重金屬污染會(huì)導(dǎo)致農(nóng)作物的污染和減產(chǎn)、品質(zhì)下降，同時(shí)誘發(fā)癌癥危害人體的健康，污染鄰近地塊的生態(tài)環(huán)境[6-8]。

目前，國(guó)內(nèi)外對(duì)于重金屬污染的空間分析做了許多研究，如Mapanda等[9]分析了津巴布韋首都某長(zhǎng)期用污水灌溉的耕地，發(fā)現(xiàn)污水灌溉致使土壤中Zn、Cu、Cd有明顯的聚集；Benhaddya等[10]采用克里金插值法分析阿爾及利亞某地區(qū)的土壤重金屬空間分布特征，并進(jìn)行了健康風(fēng)險(xiǎn)評(píng)價(jià)；Lopes等[11]分析巴西某冶煉廠周邊土壤重金屬的污染特征，結(jié)果顯示Cd和Zn空間聚集特征顯著；任加國(guó)等[12]將半變異函數(shù)理論與地統(tǒng)計(jì)學(xué)方法相結(jié)合，大尺度分析土壤重金屬空間相關(guān)性；路一帆等[13]將空間克里金插值法與健康風(fēng)險(xiǎn)評(píng)價(jià)相結(jié)合，分層分析小尺度下污染場(chǎng)地的空間特征和健康風(fēng)險(xiǎn)；劉麗麗等[14]采用主成分分析法解析重金屬污染源主要成分，結(jié)合空間分布分析華南某重金屬污染場(chǎng)地的污染源；王梅霞等[15]利用半變異函數(shù)和克里金插值，研究了復(fù)雜地形條件下土壤砷空間分布特征，揭示了復(fù)雜地形區(qū)土壤砷的空間變異性；王春光等[16]結(jié)合污染源分布特征，利用反距離權(quán)重法(IDW)分析全區(qū)整體的重金屬濃度空間分布規(guī)律及原因，并對(duì)污染指數(shù)進(jìn)行評(píng)價(jià)。以上研究并未從重金屬污染的空間特征的角度分析小區(qū)域場(chǎng)地污染的相關(guān)性。

近年來(lái)，伴隨著產(chǎn)業(yè)結(jié)構(gòu)的優(yōu)化調(diào)整，許多污染產(chǎn)業(yè)搬遷或停產(chǎn)，但受其污染的地塊被遺留下來(lái)，制約了城市化的發(fā)展進(jìn)程。分析土壤重金屬污染空間分布特征有助于從污染途徑角度探究污染原因，分析其聚集特征并與場(chǎng)地本身的要素相關(guān)聯(lián)，有助于從污染源角度探究污染原因，對(duì)污染地塊的治理具有重要意義。筆者以D市某電鍍廠為研究區(qū)，分析場(chǎng)地區(qū)域的重金屬空間分布特征，以期為城市內(nèi)部污染場(chǎng)地的治理修復(fù)提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)電鍍廠位于D市東北部，廠區(qū)地塊總面積為29 076 m2，由電鍍區(qū)、污水處理區(qū)、供電房、食堂、水泥地面組成。廠區(qū)位于北回歸線以南，屬于亞熱帶季風(fēng)氣候，年降水量豐富，年平均氣溫為22.6℃。D市制造業(yè)實(shí)力雄厚，電鍍廠所屬的化工業(yè)是支柱型產(chǎn)業(yè)。

1.2 數(shù)據(jù)的采集

研究所用數(shù)據(jù)來(lái)自生態(tài)環(huán)境部土壤與農(nóng)業(yè)農(nóng)村生態(tài)環(huán)境監(jiān)管技術(shù)中心提供的污染場(chǎng)地調(diào)查報(bào)告[17]。場(chǎng)地土壤數(shù)據(jù)于2016年3月采集。由于廠區(qū)內(nèi)部存在混凝土等硬質(zhì)地面，在數(shù)據(jù)采集的過(guò)程中使用專業(yè)鉆探設(shè)備，對(duì)水泥地面和廠房區(qū)內(nèi)的采樣點(diǎn)進(jìn)行破碎和土壤采樣。采集襯管內(nèi)的土壤樣品，將取樣管割開(kāi)，劃去表面土壤，均勻采集樣品用于測(cè)定土壤理化性質(zhì)和重金屬濃度。所有樣品采集后放入低溫保溫箱中，并及時(shí)送檢。土壤理化性質(zhì)由上海市環(huán)境測(cè)試服務(wù)部實(shí)驗(yàn)室測(cè)定。

研究區(qū)內(nèi)共設(shè)置33個(gè)土壤樣本點(diǎn)，按照5點(diǎn)采樣法，每個(gè)點(diǎn)位采集5個(gè)土壤樣本，共165個(gè)土壤樣本。采樣深度為0～0.5 m，采樣點(diǎn)在空間上的分布如圖1所示。

圖1 研究區(qū)采樣點(diǎn)位Fig.1 Sampling points in the study area

共測(cè)定Cd、Hg、Pb、Ag、Cr、Ni、Cu、Zn 8種重金屬，將各點(diǎn)位重金屬濃度參照GB 15618——2018《土壤環(huán)境質(zhì)量 農(nóng)用地土壤污染風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)（試行）》中的標(biāo)準(zhǔn)值，分析研究區(qū)域內(nèi)重金屬污染特征。為進(jìn)行空間分析，在ArcGIS 10.6軟件中將研究區(qū)場(chǎng)地信息按照WGS 1984坐標(biāo)系矢量化（圖1），場(chǎng)地共分為電鍍區(qū)、污水處理區(qū)、供電房、食堂及廠區(qū)道路5個(gè)部分。

1.3 研究方法

1.3.1 空間插值分析

空間插值是通過(guò)采樣點(diǎn)所獲得的離散數(shù)據(jù)，按照插值公式原理估算未知區(qū)域的連續(xù)數(shù)據(jù)的過(guò)程，即根據(jù)周邊相鄰采樣點(diǎn)的值來(lái)求得未知點(diǎn)的值[18]。IDW以插值點(diǎn)與采樣點(diǎn)之間的距離為權(quán)重，離插值點(diǎn)距離越近的采樣點(diǎn)的權(quán)重越大[19]。各采樣點(diǎn)值對(duì)插值點(diǎn)值作用的權(quán)重大小是成比例的，這些權(quán)重的總和為1。其公式為：

式中：Z(x0)為樣點(diǎn)x0處的預(yù)測(cè)值；n為用于插值樣點(diǎn)的數(shù)量；Z(xi)為樣點(diǎn)xi處的觀測(cè)值； λi為第i個(gè)樣點(diǎn)對(duì)預(yù)測(cè)點(diǎn)作用的權(quán)重，其公式如下：

式中：di為預(yù)測(cè)點(diǎn)與各樣點(diǎn)間的距離；p為距離的冪，在預(yù)測(cè)過(guò)程中，權(quán)重隨著樣點(diǎn)與預(yù)測(cè)點(diǎn)之間距離的增加而減小。

1.3.2 空間自相關(guān)分析

空間自相關(guān)性指位于同一分布區(qū)內(nèi)的觀測(cè)數(shù)據(jù)的不同變量之間存在的潛在相關(guān)性，莫蘭指數(shù)（Moran'sI）用于描述這些變量在空間上的聚集程度?？臻g自相關(guān)分為全局空間自相關(guān)（Global Moran'sI）和局部空間自相關(guān)（Local Moran'sI）[20]。全局空間自相關(guān)是在整個(gè)研究區(qū)范圍內(nèi)對(duì)待分析要素的地理屬性的聚集或離散程度進(jìn)行分析；局部空間自相關(guān)適用于當(dāng)存在全局空間自相關(guān)時(shí)，分析是否存在空間異質(zhì)性并確定空間異常值點(diǎn)位置，在不存在全局空間自相關(guān)時(shí)，尋找可能被掩蓋的局部空間自相關(guān)的位置[21]。

式中：Wij為空間權(quán)重矩陣；Xi和Xj分別為采樣點(diǎn)i和j上的測(cè)量值；為Xi的平均值；S2為Xi的方差；I為全局莫蘭指數(shù)。當(dāng)I＞0時(shí)，表示所分析的變量在空間上呈正相關(guān)，其值越大，空間相關(guān)性越明顯；當(dāng)I=0時(shí)，表示空間上不具有相關(guān)性；當(dāng)I＜0時(shí)，表示變量在空間上呈負(fù)相關(guān)，其值越小，空間差異越明顯。局部空間自相關(guān)分析可得出局部關(guān)聯(lián)分析的Lisa圖，可區(qū)分空間各點(diǎn)位之間高值（HH）聚類、低值（LL）聚類、由高值圍繞的低值異常值（LH）聚類以及由低值圍繞的高值異常值（HL）聚類等特征。

空間自相關(guān)分析中常用標(biāo)準(zhǔn)化統(tǒng)計(jì)量的臨界值（Z）檢查所分析要素是否在空間上具有顯著性相關(guān)性，其計(jì)算公式如下：

式中：E(I)為待分析要素也就是全局自相關(guān)的期望值；VAR(I)表示所分析變量全局莫蘭指數(shù)的方差。Z越大說(shuō)明所分析的要素在空間上的相關(guān)性越高，當(dāng)|Z|＞2.58時(shí)，則說(shuō)明所分析的變量在空間上差異性顯著，不存在隨機(jī)分布的可能性。

1.3.3 半變異函數(shù)分析

半變異函數(shù)屬于地統(tǒng)計(jì)分析中的函數(shù)的一種，能夠較為準(zhǔn)確地描述區(qū)域化變量的空間變化特性[22]。半變異函數(shù)分析可以直接對(duì)變量的空間相關(guān)性進(jìn)行分析，常用于研究變量的空間性質(zhì)隨距離變化的趨勢(shì)，得出研究變量的空間分布規(guī)律[23]。計(jì)算公式如下：

式中：h為采樣點(diǎn)之間步長(zhǎng)的變化量；r(h)為采樣點(diǎn)之間的半變異值；N(h)為在空間中步長(zhǎng)向量均為h時(shí)的離散點(diǎn)對(duì)數(shù)量；Y(xi)為位置在x上的區(qū)域化變量。r(h)半變異函數(shù)值隨著h的增大而增加，在半變異函數(shù)趨于平穩(wěn)時(shí)，離散點(diǎn)對(duì)間隔為0，步長(zhǎng)范圍越小。

2 結(jié)果與討論

2.1 土壤重金屬特征統(tǒng)計(jì)及相關(guān)性

2.1.1 土壤重金屬濃度

各類重金屬濃度統(tǒng)計(jì)見(jiàn)表1。為保證場(chǎng)地污染的篩選值的科學(xué)性，研究采用的標(biāo)準(zhǔn)限值結(jié)合了多個(gè)場(chǎng)地標(biāo)準(zhǔn)文件。由表1可知，重金屬Hg、Pb、Cd、Ag濃度的平均值小于當(dāng)?shù)乇尘爸?，最大值不超過(guò)標(biāo)準(zhǔn)限值，在研究區(qū)場(chǎng)地不存在污染情況。Ni、Cu、Zn濃度的平均值均超過(guò)的標(biāo)準(zhǔn)限值，Cr濃度的平均值雖未超過(guò)標(biāo)準(zhǔn)限值，但Ni、Cr、Cu、Zn濃度的最大值均遠(yuǎn)高于當(dāng)?shù)乇尘爸?，最大超?biāo)點(diǎn)位的超標(biāo)倍數(shù)分別為9.55、1.35、5.94、10.67，在研究區(qū)場(chǎng)地造成污染。這4種重金屬濃度標(biāo)準(zhǔn)差較大，說(shuō)明其在土壤中的分布具有顯著的差異，空間特征明顯。

表1 土壤重金屬濃度統(tǒng)計(jì)Table 1 Statistics of heavy metal contents in soil

2.1.2 土壤超標(biāo)重金屬之間的相關(guān)性

對(duì)各采樣點(diǎn)位的重金屬濃度進(jìn)行相關(guān)性分析，分析各重金屬濃度平均值的關(guān)聯(lián)性。Ni、Cr、Cu、Zn濃度之間的相關(guān)系數(shù)如表2所示。從表2可以得出，Ni與Cr濃度之間具有較強(qiáng)的空間相關(guān)性，相關(guān)系數(shù)為0.68，說(shuō)明Ni濃度高的采樣點(diǎn)位，Cr的濃度也較高，二者在研究區(qū)的分布特征較為相似；同理，Ni與Cu、Cu與Cr、Cr與Zn濃度之間同樣具有相關(guān)性，但相關(guān)性較弱，相關(guān)系數(shù)分別為0.56，0.45，0.32。而Ni與Zn、Cu與Zn濃度之間沒(méi)有相關(guān)性，各采樣點(diǎn)重金屬濃度的差異較大。

表2 超標(biāo)重金屬之間的相關(guān)系數(shù)Table 2 Correlation coefficient of over-standard heavy metals

2.2 土壤超標(biāo)重金屬空間分布特征分析

2.2.1 全局趨勢(shì)分析

全局趨勢(shì)指的是所要分析的變量在空間整體上的變化趨勢(shì)[27]。將研究區(qū)二維平面上的采樣點(diǎn)結(jié)合4種重金屬濃度轉(zhuǎn)化成三維視圖，運(yùn)用多項(xiàng)式分析投影擬合，確保采樣點(diǎn)數(shù)據(jù)最大限度地分布在趨勢(shì)面上。趨勢(shì)面分析可以在三維立體的角度更直觀地反映重金屬在采樣點(diǎn)范圍內(nèi)的分布特征。4種重金屬的空間趨勢(shì)三維視圖如圖2所示。由圖2可知，在場(chǎng)地大部分區(qū)域Ni濃度為500～1 000 mg/kg，超過(guò)標(biāo)準(zhǔn)限值；在場(chǎng)地大部分區(qū)域Cr濃度均為0～400 mg/kg，低于標(biāo)準(zhǔn)限值，只在場(chǎng)地的東北位置，也即是北側(cè)的電鍍中心區(qū)域濃度超過(guò)標(biāo)準(zhǔn)限值；Cu在場(chǎng)地大部分區(qū)域濃度均處于標(biāo)準(zhǔn)限值（500 mg/kg）以上，而且場(chǎng)地整體的變化幅度較大；Zn在場(chǎng)地只有2處區(qū)域濃度遠(yuǎn)超標(biāo)準(zhǔn)限值（700 mg/kg），其余區(qū)域均低于標(biāo)準(zhǔn)限值。

圖2 重金屬空間趨勢(shì)面圖像Fig.2 Spatial trend analysis image of heavy metals

全局趨勢(shì)分析的結(jié)果可知，Ni的濃度在空間上呈北高南低，東西方向變化不明顯；Cr的濃度在場(chǎng)地東北高、西南低；Cu的濃度在場(chǎng)地整體變化幅度大；Zn的濃度屬于中間高四周低。

2.2.2 空間預(yù)測(cè)分析

將采樣點(diǎn)的4種超標(biāo)重金屬濃度平均值作為插值進(jìn)行空間插值分析，得出超標(biāo)重金屬的插值圖像（圖3）。從圖3可知，Ni濃度高值區(qū)主要集中在北側(cè)電鍍區(qū)域、污水處理廠區(qū)及廠區(qū)道路，根據(jù)電鍍廠的生產(chǎn)方式，初步推斷電鍍生產(chǎn)、污水處理不達(dá)標(biāo)和運(yùn)輸過(guò)程泄露是Ni的主要污染來(lái)源；Cr濃度高值區(qū)主要集中在北側(cè)電鍍區(qū)域并向周邊呈逐漸降低的趨勢(shì)，初步推斷電鍍生產(chǎn)是Cr的主要污染來(lái)源；Cu的富集區(qū)在圖像中較多，在南北側(cè)的電鍍區(qū)、污水處理廠區(qū)和廠區(qū)道路均有分布，而且呈高值的點(diǎn)位較多且集中，說(shuō)明電鍍生產(chǎn)與運(yùn)輸過(guò)程的泄露是導(dǎo)致廠區(qū)內(nèi)Cu濃度超標(biāo)的主要原因，并且超標(biāo)覆蓋的面積較大；Zn的高值區(qū)在南北側(cè)電鍍中心區(qū)域，說(shuō)明是該區(qū)域的電鍍生產(chǎn)造成Zn的超標(biāo)。綜上，電鍍生產(chǎn)是導(dǎo)致Ni、Cr、Cu、Zn超標(biāo)的共同原因，同時(shí)污水處理不達(dá)標(biāo)、運(yùn)輸過(guò)程泄露導(dǎo)致Ni、Cu在廠區(qū)的污水池和水泥道路區(qū)域超標(biāo)。

從圖3看出，4種重金屬濃度的分布各不相同，Ni和Cr濃度高值分布區(qū)域比較近似，核心區(qū)域均在北側(cè)電鍍區(qū)廠區(qū)邊界中間處。結(jié)合4種超標(biāo)重金屬濃度的相關(guān)性分析，Ni與Cr的分布具有較高的相關(guān)性，基本與4種重金屬的空間圖像分布特性相符。

圖3 場(chǎng)地區(qū)域超標(biāo)重金屬濃度的空間插值分布Fig.3 Spatial interpolation distribution image of over-standard heavy metals in the site area

2.3 土壤超標(biāo)重金屬空間相關(guān)性分析

2.3.1 半變異函數(shù)分析

根據(jù)空間相關(guān)定律，空間上距離越近的事物具有更大可能的相似性，半變異函數(shù)云分析可以定量化這種相似性。地物的空間異質(zhì)性隨著方向發(fā)生變化，為探索場(chǎng)地超標(biāo)重金屬在各方向的變化特征，搜索 0°（正南北方向）、45°（東北——西南方向）、90°（正東西方向）和135°（西南——東北）方向，角度容差為±22.5°。利用 ArcGIS 10.6軟件中地統(tǒng)計(jì)分析中的半變異函數(shù)分析，對(duì)Ni、Cu、Zn、Cr濃度在不同方向上的空間異質(zhì)性進(jìn)行分析。以Ni為例，結(jié)果如圖4所示。橫坐標(biāo)為步長(zhǎng)即采樣點(diǎn)之間的距離，縱坐標(biāo)為半方差的大小，半方差越小空間相似性越高。由圖4可見(jiàn)，大部分采樣點(diǎn)之間的半方差較小，僅小部分樣點(diǎn)的半方差較大，說(shuō)明Ni濃度在一定的步長(zhǎng)范圍具有較高的空間相關(guān)性。通過(guò)對(duì)比不同方向發(fā)現(xiàn)，Ni濃度的半變異函數(shù)云半方差高值點(diǎn)最少，在該方向Ni濃度的空間相關(guān)性最高。

圖4 超標(biāo)重金屬半變異函數(shù)云分析圖Fig.4 Semivariograms cloud analysis diagram of four over-standard heavy metals

2.3.2 全局空間相關(guān)性分析

利用ArcGIS 10.6軟件中的全局空間關(guān)聯(lián)分析，對(duì)4種超標(biāo)重金屬采樣點(diǎn)的濃度進(jìn)行空間自相關(guān)分析，得出4種重金屬濃度的全局莫蘭指數(shù)（I）及相關(guān)參數(shù)（表3）。由表3可知，研究區(qū)內(nèi)重金屬Ni和Cr濃度在空間上呈正相關(guān)，即空間上距離越相近，其濃度越接近，Cu和Zn在空間上呈負(fù)相關(guān)，空間上距離越相近，其濃度差距越大。按照相關(guān)性大小依次排序?yàn)?Cr＞Cu＞Zn＞Ni；Z表示所分析變量空間隨機(jī)分布的顯著性，只有重金屬Cr的Z＞2.58，說(shuō)明其空間分布存在明顯的聚類特征。

表3 超標(biāo)重金屬的空間自相關(guān)結(jié)果Table 3 Spatial autocorrelation results of four over-standard heavy metals

2.3.3 局部冷熱點(diǎn)相關(guān)性分析

通過(guò)全局自相關(guān)分析，對(duì)整個(gè)場(chǎng)地所有采樣點(diǎn)進(jìn)行了空間相關(guān)分析，但對(duì)部分特殊點(diǎn)位的相關(guān)分析和空間特征分析仍不夠清晰。為挖掘可能存在的局部點(diǎn)位之間的關(guān)系，對(duì)4種超標(biāo)重金屬分別進(jìn)行局部自相關(guān)分析，進(jìn)一步探索局部空間分布的特征。利用ArcGIS 10.6軟件Anselin Local Moran'sI分析模塊，并設(shè)置相應(yīng)參數(shù)，得到研究區(qū)土壤重金屬濃度聚類分布（圖5）。從高值聚集點(diǎn)（熱點(diǎn)）來(lái)看，Ni、Cr熱點(diǎn)分布在北側(cè)電鍍區(qū)中間位置；從低值聚類點(diǎn)（冷點(diǎn)）來(lái)看，Ni、Cr、Cu、Zn 均存在，整體上冷點(diǎn)分布在電鍍區(qū)的東南部分；從被高值圍繞的低值異常點(diǎn)來(lái)看，4種重金屬均不存在這類特征點(diǎn)；從被低值圍繞的高值異常點(diǎn)來(lái)看，Ni、Cr 2種重金屬存在，均分布在北側(cè)電鍍區(qū)與熱點(diǎn)相鄰處。4種金屬的采樣點(diǎn)多為無(wú)顯著聚集特征點(diǎn)，只有部分點(diǎn)位存在空間聚集特征。

Ni、Cr、Cu、Zn 均存在聚類特征，這些特征點(diǎn)可能受到場(chǎng)地自身的屬性的影響。圖5中的高值點(diǎn)結(jié)合場(chǎng)地空間屬性，可能是電鍍生產(chǎn)或殘留污染物泄露導(dǎo)致的；低值點(diǎn)位靠近墻體和電鍍區(qū)的拐角，可能是邊角區(qū)域不易產(chǎn)生污染。這些特殊點(diǎn)的存在除與場(chǎng)地在空間屬性相關(guān)外，也可能與地面的損毀程度、土質(zhì)有關(guān)。

圖5 超標(biāo)重金屬局部自相關(guān)圖Fig.5 Local autocorrelation diagram of over-standard heavy metals

3 結(jié)論

（1）研究區(qū)電鍍廠部分點(diǎn)位重金屬 Ni、Cr、Cu、Zn濃度超過(guò)GB 15618——2018中標(biāo)準(zhǔn)限值，對(duì)研究區(qū)造成污染，其中Ni、Cu超標(biāo)點(diǎn)位較多。

（2）電鍍生產(chǎn)區(qū)是 Ni、Cr、Cu、Zn 濃度超標(biāo)的共同區(qū)域，有可能是電鍍生產(chǎn)或遺留污染物的泄露導(dǎo)致的。4種重金屬濃度在場(chǎng)地的分布特征不同，Ni在廠區(qū)北高南低；Cr只在東北角區(qū)域濃度超標(biāo)；Zn屬于中間高四周低；Ni和Cr高值分布區(qū)域均在北側(cè)電鍍區(qū)廠區(qū)邊界中間處。

（3）相關(guān)性分析表明，Ni與Cr濃度在廠區(qū)的空間分布上具有較強(qiáng)的相關(guān)性，相關(guān)系數(shù)為0.68，其余重金屬間的相關(guān)性較弱。研究區(qū)內(nèi)土壤 Ni、Cr濃度呈現(xiàn)空間正相關(guān)，Cu、Zn濃度呈現(xiàn)空間負(fù)相關(guān)，Cr濃度空間分布上存在明顯的聚類特征。

（4）從局部冷熱點(diǎn)分析的結(jié)果可知，4種重金屬均存在空間聚集特征。高值點(diǎn)可能是電鍍生產(chǎn)或殘留污染物泄露導(dǎo)致的，低值點(diǎn)可能由于位于邊角區(qū)域，故未產(chǎn)生污染。

（5）本研究的不足之處在于受限于場(chǎng)地的調(diào)查數(shù)據(jù)，未獲得場(chǎng)外的土壤重金屬數(shù)據(jù)，進(jìn)行的分布預(yù)測(cè)存在誤差；重金屬的分布預(yù)測(cè)只按照濃度進(jìn)行預(yù)測(cè)，并未和生態(tài)、健康風(fēng)險(xiǎn)關(guān)聯(lián)，預(yù)計(jì)在后續(xù)的研究中補(bǔ)充。