李瑞娟,周冰
(1. 西安航空學院能源與建筑學院,陜西 西安 710077;2.青海省地質礦產測試應用中心,青海 西寧 810000)
隨著礦產資源的不斷開采,對當地自然環境造成嚴重的生態破壞,其中以土壤重金屬污染現象尤為嚴重。目前,國內針對礦區土壤金屬污染評價的方法日趨成熟,包括內梅羅指數法、地累積指數法和潛在生態危害指數法等[1],但對于土壤污染的來源及分布研究并不全面,因此在運用以上方法進行金屬污染評價的基礎上,對污染元素的來源及空間分布特征進行研究,對查明重度污染區域及今后的污染防治具有十分重要的現實意義。
安徽省某金屬尾礦庫位于銅陵市東北部,三面環山,北有小峰山,南有烏木山,東北壩下有鐘團村。該尾礦庫主要包括銅陵銅尾礦、楊山沖尾礦、銅官山尾礦和金口嶺尾礦,設計總庫容約為540萬m3,目前尾砂堆存約670 m3,占地面積約44萬m2。其中銅陵銅尾礦廢棄地是一處存在重金屬Cu毒害的原生裸地,含有大量金屬硫化物和有毒物質,不僅對周圍環境造成了嚴重污染與破壞,還通過食物鏈對人體健康造成危害。
目前該尾礦廢棄地下游仍住有居民約260戶,農田約70 km2,且布置有公路和鐵路設施。由于該地區屬于亞熱帶季風氣候,夏季多雨,年平均降水量約1400 mm,尾礦重金屬元素易通過降雨排水和地表徑流等方式沖刷至下游[2],對下游居民安全及交通設施造成危害。故本文針對銅陵銅尾礦廢棄地及周邊55 km2范圍內土壤進行金屬污染評價及綜合利用研究,具有一定的現實意義,也為其他類似礦區提供參考價值。
為查明土地質量現狀,了解土壤污染狀況,采集研究區0~20 cm表層樣品。表層土壤樣點均勻布置,布設密度采用20個點/km2,根據實際情況,共采集表層土壤樣品754個,其采樣具體方法如下:
(1)樣品采集時,刮去地表植物凋落物,以刻槽方式采集地表至20 cm深處的土柱,并去除雜草、礫石、肥料團塊等雜物[3]。
(2)在每處采樣時,要在周圍10 m范圍內至少采集3處以上土壤樣品組合成一件樣品,采樣點坐標記錄中間采樣點GPS顯示的坐標。
(3)所采集的樣品原始重量應不低于2000 g,以保證自然干燥過篩后,截取的粒級樣品(<0.8 mm)重達1000 g[4]。
(4)采樣用的布袋要經過洗滌,比較濕的和污染較重的樣品要用塑料口袋裝好隔離,避免與其他樣品相互污染。
為了充分考察土壤污染情況,查明土壤中污染元素的含量及分布特征,故對土壤樣品中污染元素進行測試分析,包括砷、鎘、鉻、鉛、汞、鎳、銅和鋅共計8項,其測試方法及儀器見表1。

表1 測試方法及儀器統計Table 1 Test methods and instrument statistics
土壤中存在的重金屬污染大多表現為綜合污染,故應在傳統單因子指數法的基礎上將同個樣品中不同污染物的污染指數和同個污染物中不同樣品的污染指數兩者進行綜合,從而達到對單個污染物或采樣點綜合評價[5]。其中,內梅羅指數法是國內對土壤環境質量進行綜合污染指數評價普遍方法之一,它能夠全面考慮不同重金屬元素對土壤環境質量的影響,綜合反映出不同污染物對土壤環境的不同作用[6]。其計算公式如下:

式中:PN-土壤綜合污染指數;Ci?Si -土壤中第i種污染物的單項污染指數;(Ci?Si )max-土壤污染物中的最大單項污染指數;(Ci?Si )ave-土壤中各類污染物的單項污染指數的平均值。
根據內梅羅綜合污染指數可將土壤劃分為5個等級。
隨著GIS 技術的不斷發展,可運用地理信息系統軟件Arc GIS 10.2分析土壤中金屬的空間分布狀況[7]。其中影響該軟件對土壤中金屬元素空間分布預測精度的因素包括采樣點的數量,布點的空間結構和空間插值方法。
本次研究共采集樣點數量754個,滿足實際需求;布點方法可基于GIS 技術準確快速地直觀表達,同時將土壤屬性數據與樣點空間數據相關聯,從而建立樣點數據庫便于數據保存與實時更新;空間插值方法選用普通克里金插值法對研究區土壤中金屬含量進行內插,把已知數據點之間的未知區域部分通過插值來進行估算并賦值,通過這些離散點的特征數據得出連續的空間數據集,然后按不同濃度值重新分類后分級著色得出各金屬元素的空間分布狀況圖,從而對研究區土壤的重金屬含量和污染狀況全面把握。
目前,銅陵銅尾礦土地貧瘠,土壤養分匱乏,富含重金屬,遺留大片廢地,僅在邊緣區域存在少數植物存活,而中央區域幾乎無植物生長。
根據樣品測試結果,運用SPSS 22.0統計分析軟件對上述754個土壤樣品的測試指標進行描述性分析,分別統計其極小值、極大值、中位數、平均值、標準差和變異系數,并參考全國土壤污染狀況調查中當地背景值[8],見表3。
由表3可得,研究區土壤中Hg元素變異系數最大,高達233%,具有較強的變異性,而其他金屬元素的變異系數較低,均在[0,1]范圍,具有中等變異性,這說明人為活動對礦區內土壤化學組成成分影響有限。同時根據各金屬元素的平均含量與當地背景值對比可得,目前研究區內Pb、Cu、Zn和Hg四種元素已經在土壤中開始累積。
運用內梅羅指數法對研究區754 個土壤樣品中的砷、鎘、鉻、鉛、汞、鎳、銅和鋅八種金屬元素含量計算其單項污染指數,以國家土壤環境質量標準(GB 15618-2018)中pH≤6.5 的二級標準為評價標準,按公式(1)計算樣品中污染因子的污染指數,其評價結果見表4。并根據表2對樣品劃分污染等級程度,其結果見表5。

表2 土壤質量分級標準Table 2 Standard table of soil quality classification

表4 污染評價結果Table 4 Pollution evaluation results

表5 污染等級結果Table 5 Pollution level results
由表4可得,研究區土壤中Hg和Cu元素的超標現象最嚴重,各采樣點中Hg和Cu的綜合污染指數分別為3.822和3.641,達到重度污染;其次是 Cd,其綜合污染指數為 1.376,構成輕度污染;Ni和Zn兩種元素除部分存在輕微污染外,其綜合污染指數表明研究區總體上并未受到污染;As、Cr和Pb三種元素均處于安全水平,沒有造成污染。結合綜合污染指數可得,研究區土壤采集樣品中19.8%的采樣點已達到重度污染,這將嚴重威脅當地環境安全。根據樣點分布及現場勘查推測重度污染區域主要分布在采區和尾礦廢棄地附近。
利用Arc GIS 10.2中的 Geostatistical Analysis 地統計分析模塊,采用普通克里金(Kriging)插值法繪制研究區污染元素在空間上的分布狀況,見圖1。

圖1 污染元素分布Fig. 1 Distribution of pollution elements
圖1展示了研究區內土壤8種污染元素含量的空間分布。由圖1可知,砷、汞、銅、鉻和鋅元素具有相似的空間分布特征,呈現以采區為中心向四周遞減的趨勢,元素含量較高的區域集中分布在采區和尾渣堆積地,其他區域含量隨距采區距離的增大而遞減。其中,研究區中含量較高的污染元素為銅元素,最大值為198.670 mg/kg,最小值為13.258 mg/kg,遠大于其他污染元素的含量,并通過圖中空間分布可得銅元素集中分布于采區和尾礦廢棄地,而其他區域含量幾乎為0,分析其原因可能和礦區不充分開采有關,造成大量銅元素殘留,另外尾礦廢棄地自身金屬元素含量偏高,并伴隨著降雨、刮風和人為活動等因素不斷向周圍遷移。
根據上述研究區目前存在元素Pb、Cu、Zn和Hg累積,借鑒曹紅梅對香根草的研究[9],在重度污染區域隨機選擇10個樣點處種植香根草以期研究香根草對重金屬的富集吸收及轉運能力。根據表1中測試方法及儀器分別測試土壤樣品、香根草根部、莖部內Pb、Cu、Zn和Hg四種金屬元素的含量,取其平均值統計結果可見表6。

表6 香根草富集系數與轉運系數統計Table 6 Statistics of the enrichment coefficient and transport coefficient of vetiver
富集系數BCF[10]是評價植物對于重金屬富集能力的重要指標,轉運系數TCF[11]是體現植物將重金屬離子由根部運輸至地上部的主要指標。由表6可知,香根草對Pb、Cu、Zn和Hg四種元素中,僅有Cu的富集系數大于1,其他均小于1,說明香根草根系對Cu的富集能力很強,可以作為研究區生態恢復的備用植物,但其轉運系數較低,如何在不破壞根系的前提下有效利用香根草將Cu轉移出土壤進一步利用仍待研究。
(1)通過測試指標分析和內梅羅污染評價可得,目前研究區土壤中累積污染元素主要為Pb、Cu、Zn和Hg,其中尤以Hg和Cu金屬污染最為嚴重,其綜合污染指數分別達3.822和3.641,且19.8%的采樣點已達到重度污染程度。
(2)通過Arc GIS及普通克里金插值法,得出As、Hg、Cr和Zn元素具有相似的空間分布特征:以采區為中心向四周遞減。其中含量較高的污染元素為Cu,集中分布于采區和尾礦廢棄地,并未遷移至下游居民區。這與現場勘查重度污染樣點分布地點一致,說明了引入Arc GIS 對污染元素空間分布狀況研究的可行性,為后續綜合治理做下了鋪墊。
(3)針對研究區土壤銅污染問題,可選擇香根草作為生態恢復植物。香根草對Cu的富集系數為1.31,具有很強的富集能力,可有效降低土壤中Cu元素的含量,但其轉運系數較低,如何在不破壞根系的前提下有效利用香根草將Cu轉移出土壤進一步利用仍待研究。