他維媛 周書宇 金盛華 趙東旭 孟昭君 梁麗華,3#

(1.陜西省環(huán)境調(diào)查評估中心,陜西 西安 710054;2.西北大學城市與環(huán)境學院,陜西 西安 710127;3.陜西省地表系統(tǒng)與環(huán)境承載力重點實驗室,陜西 西安 710127)

近年來,我國各類礦山采選行業(yè)發(fā)展迅猛,由礦產(chǎn)資源的開發(fā)和冶煉工業(yè)引起的土壤重金屬污染問題因其嚴重性而受到國內(nèi)外的廣泛關注[1]。尾礦庫屬于高環(huán)境風險設施場所,其中通常含有選礦后遺留的重金屬和其他有害物質(zhì),大量的尾礦堆積不僅會侵占土地資源,尾礦滲濾液的滲漏還可能導致區(qū)域性土壤重金屬含量超標[2],對周圍地區(qū)的生態(tài)環(huán)境和居民的健康造成極大威脅,因此,對礦區(qū)土壤重金屬污染現(xiàn)狀及其來源進行分析具有重大意義。

商洛市礦產(chǎn)種類多、儲量大,目前已有一些學者對該地區(qū)土壤重金屬的污染做了研究。例如楊維鴿等[3]20的研究表明商州區(qū)尾礦庫周邊耕地土壤中Pb和Cd為主要污染物。劉夢云等[4]對商洛市煙草用地重金屬污染狀況進行了評價,結果表明土壤中存在不同程度的Cd污染。高利峰[5]66對商洛市不同功能區(qū)土壤重金屬污染特征進行研究,結果表明當?shù)赝寥乐饕蹾g和Cd的污染?，F(xiàn)有的對商洛市土壤重金屬的研究往往只使用一種評價方法,而各單一評價方法存在的局限性通常導致評價結果不夠全面。因此在對土壤重金屬進行污染評價時,需要選用不同的評價方法,綜合考慮重金屬的累積程度、生物毒性、遷移特性、污染源特點等方面,才能全面客觀地反映土壤重金屬的污染情況。本研究以商洛市13個尾礦庫周邊土壤為研究對象,選取Cd、Cr、As、Pb、Cu、Ni共計6種重金屬,采用多種評價方法對其污染特征和來源進行研究,以期為該地土壤重金屬污染的治理與管控提供數(shù)據(jù)支撐和科學依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

商洛市位于陜西省南部,市內(nèi)礦產(chǎn)資源豐富,已發(fā)現(xiàn)各類礦產(chǎn)資源60多種,全市共有尾礦庫118座,占陜西省的37.7%。

此次研究涉及的尾礦庫主要礦種包括鉬礦、釩礦、銅礦、浮選金礦、涉氰金礦、非硫鐵礦,現(xiàn)存庫容從12萬m3至116萬m3不等,運行狀況多為停用和閉庫。研究區(qū)內(nèi)13個尾礦庫的基本情況如表1所示。

表1 尾礦庫基本信息Table 1 Basic information of tailings ponds

1.2 樣品采集與分析

按照《土壤環(huán)境監(jiān)測技術規(guī)范》(HJ/T 166—2004)的要求,在商洛市13個尾礦庫內(nèi)布設采樣點(見圖1),著重對土壤有污染跡象或可能受污染的區(qū)域進行布點采樣。采樣時剔除粒徑較大的石塊等雜質(zhì),份樣量大于1 kg。

圖1 研究區(qū)采樣點分布示意圖Fig.1 Schematic diagram of sampling sites in the study area

采樣完成后,參照《工業(yè)固體廢物采樣制樣技術規(guī)范》(HJ/T 20—1998)制樣及保存。Cr、Pb、Cu、Ni采用火焰原子吸收分光光度法檢測,Cd采用電感耦合等離子體質(zhì)譜法檢測,As采用微波消解—原子熒光法檢測。

1.3 評價方法

本研究采用4個單一指數(shù)——污染因子、單因子污染指數(shù)、潛在生態(tài)風險因子和地累積指數(shù)和3個綜合指數(shù)——污染負荷指數(shù)、內(nèi)梅羅綜合污染指數(shù)和潛在生態(tài)風險指數(shù)進行評價。單一指數(shù)用于篩選土壤中主要污染元素,綜合指數(shù)用于評估尾礦庫和研究區(qū)的污染情況[6]。在采用多種污染指數(shù)對重金屬的污染特征進評價的基礎上,對不同的指標和評價方法也進行了比較,以期為評價方法的選擇提供參考。

污染因子和污染負荷指數(shù)分別用于評估單個重金屬和研究區(qū)整體污染狀況[7]。污染因子和污染負荷指數(shù)的計算方法見式(1)至式(3)。

(1)

(2)

(3)

式中:FC,i為第i種重金屬的污染因子;ci為第i種重金屬的質(zhì)量濃度,mg/kg;Bi為第i種重金屬的背景值,mg/kg,選用陜西省土壤重金屬元素背景值[8];IPL,j為采樣點j處的污染負荷指數(shù);n為重金屬種類總數(shù);IPL為研究區(qū)的污染負荷指數(shù);m為采樣點總數(shù)。

單因子污染指數(shù)法是以污染物的環(huán)境質(zhì)量標準為基準的一種評價方法,一般作為其他評價方法的基礎[9]。單因子污染指數(shù)的計算方法見式(4)。

(4)

式中:Pi為第i種重金屬的單因子污染指數(shù);Si為第i種重金屬的的標準值,mg/kg,本研究選用《土壤環(huán)境質(zhì)量 農(nóng)用地土壤污染風險管控標準(試用)》(GB 15618—2018)中的風險篩選值。

內(nèi)梅羅綜合污染指數(shù)法基于單因子污染指數(shù)法,可綜合反映各重金屬對環(huán)境的作用,同時可以突出高濃度污染物對環(huán)境的影響[10]3。內(nèi)梅羅綜合污染指數(shù)的計算方法見式(5)。

(5)

式中:PN為內(nèi)梅羅綜合污染指數(shù);Pi,ave為Pi均值;Pi,max為Pi最大值。

潛在生態(tài)風險因子考慮了毒性影響,潛在生態(tài)風險指數(shù)用于綜合評價土壤或沉積物中重金屬污染程度及潛在生態(tài)風險[11],計算公式見式(6)和式(7)。

FER,i=FC,i×Ti

(6)

IR=∑FER,i

(7)

式中:FER,i為第i種重金屬的潛在生態(tài)風險因子;Ti為第i種重金屬的毒性響應系數(shù),Cd、Cr、As、Pb、Ni、Cu的毒性響應系數(shù)分別取30、2、10、5、5、5[12];IR為潛在生態(tài)風險指數(shù)。

地累積指數(shù)可以反映重金屬污染水平和人為活動的影響,計算公式見式(8):

(8)

式中:Igeo為地累積指數(shù);1.5為修正指數(shù),通常用于表征沉積特征和其他人為影響[13]。

污染因子、污染負荷指數(shù)處于(0,1)、[1,3)、[3,6)、[6,+∞)區(qū)間分別代表輕污染、中污染、重污染和極重污染;單因子污染指數(shù)、內(nèi)梅羅綜合污染指數(shù)處于(0,0.7)、[0.7,1.0)、[1.0,2.0)、[2.0,3.0)、[3.0,+∞)區(qū)間分別代表無污染、警戒線、輕污染、中污染和重污染;潛在生態(tài)風險因子處于(0,40)、[40,80)、[80,160)、[160,320)、[320,+∞)區(qū)間分別代表低風險、中風險、較高風險、高風險和極高風險;潛在生態(tài)風險指數(shù)處于(0,150)、[150,300)、[300,60)、[600,+∞)區(qū)間分別代表低風險、中風險、較高風險和高風險;地累積指數(shù)處于(-∞,0)、[0,1)、[1,2)、[2,3)、[3,4)、[4,5)、[5,+∞)區(qū)間分別代表無污染、無至輕污染、輕污染、輕至中污染、中污染、中至重污染和重污染。

2 結果與討論

2.1 重金屬含量統(tǒng)計分析

商洛市土壤重金屬的描述性統(tǒng)計見表2。由表2可知,Cd和Pb均值為相應的標準值1.7、1.1倍,表明Cd和Pb含量超出了全國尺度上的污染物容許量;Cd、Cr、As、Pb、Ni和Cu均值分別為相應背景值的11.6、1.1、2.1、8.7、1.1、1.9倍,表明這些重金屬在研究區(qū)土壤中存在累積效應,這與湯波等[14]的研究結果一致,說明該地區(qū)采礦活動已經(jīng)對當?shù)赝寥佬纬闪宋廴?。與高利峰[5]652015年對商洛市農(nóng)業(yè)區(qū)土壤重金屬的研究結果對比,本次研究中Cd和Pb均值均高于2015年,這可能與該區(qū)域尾礦庫曾發(fā)生尾礦泄漏事故有關。利用背景值和標準值得出的評價結果的差異性,研究區(qū)土壤對外源重金屬具有一定的緩沖能力,即在特定的土壤標準下,土壤可以容納一定量重金屬而不超標[15]。國家標準一般只能起篩選作用[16],由于各地土壤類型和理化性質(zhì)差異巨大,土壤污染具有明顯的區(qū)域特征,因此在區(qū)域尺度上利用背景值判斷土壤中重金屬的累積程度更為合理。

表2 研究區(qū)土壤重金屬統(tǒng)計分析結果Table 2 Results of statistical analysis of soil heavy metal content in the study area

通常用變異系數(shù)來度量概率分布的離散程度[10]4,從表2可以看出,13個尾礦庫附近土壤中Cd、As、Pb變異系數(shù)較高,說明整個研究區(qū)這3種重金屬濃度空間分散程度較大,表明受人為干擾的可能性較大。值得注意的是,Cr和Ni的變異系數(shù)分別為19.80%和17.70%,屬于中等變異,然而土壤Cr和Ni均值卻與當?shù)乇尘爸到咏?這表明土壤中Cr和Ni可能主要受大氣沉降、成土因素等自然源的影響。

尾礦庫重金屬質(zhì)量濃度見圖2。總體上看,各尾礦庫中Pb含量最高(1～6號和12號尾礦庫尤為明顯)。楊維鴿等[3]18對商洛市尾礦庫周邊耕地土壤污染的研究也表明Pb為首要污染元素,與本研究結果一致。

圖2 尾礦庫重金屬質(zhì)量濃度Fig.2 Heavy metal content of tailings ponds

2.2 重金屬污染特征分析

2.2.1 污染因子和污染負荷指數(shù)

由圖3可以看出,重金屬污染因子均值表現(xiàn)為Cd(11.57)>Pb(8.74)>As(2.08)>Cu(1.93)>Ni(1.11)>Cr(1.06)。土壤中Cd和Pb的污染水平最高,部分采樣點達到極重污染,為主要污染元素;研究區(qū)的污染負荷指數(shù)為1.99,為中污染水平,與鄧文博等[17]的研究結果一致。其中5號尾礦庫污染負荷指數(shù)最大(3.39),為重污染水平。

圖3 研究區(qū)土壤重金屬的污染因子和尾礦庫的污染負荷指數(shù)Fig.3 Contamination factors of soil heavy metals and pollution load index of tailings ponds in the study area

2.2.2 單因子污染指數(shù)和內(nèi)梅羅綜合污染指數(shù)

由圖4可以看出,單因子污染指數(shù)均值表現(xiàn)為Cd(1.73)>Pb(1.10)>As(0.93)>Cu(0.41)>Cr(0.26)>Ni(0.17)。研究區(qū)土壤中As、Cd和Pb均有采樣點達重污染,其中Cd和Pb污染情況較重,為主要污染元素。從尾礦庫污染情況來看,5號和13號尾礦庫為重污染,內(nèi)梅羅綜合污染指數(shù)分別為10.22和3.78。內(nèi)梅羅綜合污染指數(shù)法具有格外強調(diào)最不利因素的特點,5號尾礦庫Cd的單因子污染指數(shù)最大,是導致5號尾礦庫土壤內(nèi)梅羅綜合污染指數(shù)偏大的主要原因。針對毒性強、危害大的土壤重金屬進行研究時,使用內(nèi)梅羅綜合污染指數(shù)法適當強調(diào)最不利因素的影響可以更好地引起人們的警覺和重視。

圖4 研究區(qū)土壤重金屬的單因子污染指數(shù)和尾礦庫的內(nèi)梅羅綜合污染指數(shù)Fig.4 Single factor pollution index of soil heavy metals and Nemerow pollution index of tailings ponds in the study area

2.2.3 潛在生態(tài)風險因子和潛在生態(tài)風險指數(shù)

由圖5可以看出,土壤中重金屬之間潛在生態(tài)風險因子差異很大,其中Cd的潛在生態(tài)風險因子均值達到了極高風險水平,為主要生態(tài)風險元素;這與崔雅紅等[18]對安康市石煤礦區(qū)土壤以及李堆淑等[19]對陜南某冶煉區(qū)周邊蔬菜地土壤的研究結論一致。從尾礦庫污染情況來看,尾礦庫之間潛在生態(tài)風險指數(shù)介于80.38～2 896.54,差異也較大。其中5號尾礦庫為高風險。研究區(qū)土壤整體呈現(xiàn)出低至中風險,土壤重金屬防控的重點在于Cd污染的治理。

圖5 研究區(qū)土壤重金屬的潛在生態(tài)風險因子和尾礦庫的潛在生態(tài)風險指數(shù)Fig.5 Potential ecological risk factor of soil heavy metals and ecological risk index of tailings ponds in the study area

2.2.4 地累積指數(shù)

由圖6可以看出,Cd和Pb在個別采樣點為重污染和中至重污染,兩者的Igeo均值呈現(xiàn)輕污染,表明研究區(qū)土壤中Cd和Pb出現(xiàn)富集,為主要污染元素。從尾礦庫的Igeo來看,5號尾礦庫的Cd為主要污染貢獻元素;對于1號和6號尾礦庫,Pb為主要污染貢獻元素。其余尾礦庫也存在不同程度的Cd和Pb污染。圖6表明,研究區(qū)土壤整體呈現(xiàn)出無至輕污染。CHEN等[20]采用地累積指數(shù)法對商洛市某礦區(qū)土壤重金屬進行污染評價時,也得出Cd在研究區(qū)土壤中污染現(xiàn)象普遍的結論。

圖6 研究區(qū)土壤重金屬和尾礦庫的地累積指數(shù)Fig.6 Geoaccumulation index of soil heavy metals and tailings ponds in the study area

2.3 重金屬來源分析

單一指數(shù)評價結果表明,該區(qū)域土壤中Cd和Pb為主要污染元素,研究區(qū)內(nèi)尾礦庫礦種以鉬礦為主,鉬礦的主要伴生元素為Cd和Pb[21-22],這造成了土壤中Cd和Pb的富集;交通運輸和燃煤燃燒也是Pb的主要來源之一[23-25]。此外,研究區(qū)內(nèi)大量堆積的尾礦也與土壤中Cd和Pb的富集密切相關。

研究區(qū)尾礦渣、尾礦水的監(jiān)測指標見表3,這些指標與土壤中Pb、Cd的相關性見表4。研究區(qū)尾礦水pH為5.84～8.62,尾礦渣pH為3.29～9.16,土壤中Pb含量與尾礦水、尾礦渣的pH呈顯著負相關。偏酸的環(huán)境更有利于重金屬元素從尾礦渣和尾礦水中遷移到土壤中[26-27],然后通過發(fā)生絡合或吸附而長期滯留于土壤[28],造成土壤中Pb的富集。土壤中Cd含量與尾礦水中Pb、尾礦渣中Pb和Cd含量顯著相關,是因為尾礦渣和尾礦水中的重金屬易受大氣、地表徑流、土壤下滲等作用的影響[29]遷移到礦山周圍環(huán)境中,引起土壤重金屬含量升高。因此,尾礦庫土壤中Pb、Cd呈現(xiàn)的含量特征與尾礦渣中Pb、Cd的含量特征基本一致。土壤中Cd和尾礦渣中Zn含量顯著相關,是由于Cd和Zn有相似的化學性質(zhì)和行為,在自然界常常伴生[30],兩者之間發(fā)生的協(xié)同作用會促進重金屬離子的遷移轉(zhuǎn)化,加重土壤中Cd污染。

表3 研究區(qū)尾礦渣、尾礦水的監(jiān)測指標Table 3 Monitoring indicators in tailing slag and tailing water in the study area

表4 研究區(qū)尾礦渣、尾礦水中監(jiān)測指標與土壤中Pb、Cd的相關性1)Table 4 Correlation of monitoring indicators in tailing slag and tailing water with Pb and Cd in soil in the study area

從綜合指數(shù)分析結果來看,尾礦庫的污染分布格局為東北部污染水平較高,西南部、東南部污染水平較低。這主要因為研究區(qū)西北部尾礦庫集中且礦山開發(fā)年限久遠,采礦規(guī)模較大,因此該區(qū)域土壤中重金屬具有時間累積效應。此外,1號尾礦庫在再利用過程中會產(chǎn)生大量富集Mo及其主要伴生元素Pb、Cd的微粒,這些微粒自然沉降后在該區(qū)域土壤中發(fā)生二次富集。進一步調(diào)查發(fā)現(xiàn),5號尾礦庫曾被洪水沖毀管道,造成滲濾液直接外排,而尾礦水中的選礦藥劑可以使尾礦區(qū)域的pH更低[31],更有利于該區(qū)域土壤中重金屬遷移與富集。

不同指標計算得出的評價結果不同是因為重金屬污染評價方法的特點及適用性各不相同。污染因子僅考慮了重金屬的含量與背景值之間的關系,污染負荷指數(shù)僅考慮了重金屬含量與標準值之間的關系,兩者只能簡單說明重金屬在土壤中的富集程度[32];地累積指數(shù)在使用中修正系數(shù)的選擇會對計算結果造成較大影響,也存在一定的局限性[33];內(nèi)梅羅綜合污染指數(shù)突出了污染指數(shù)最大的污染物對環(huán)境的影響和貢獻作用;潛在生態(tài)風險指數(shù)側(cè)重于評價重金屬對環(huán)境的潛在風險,考慮了重金屬的毒性作用,但未考慮重金屬在土壤中的遷移特性[34];污染負荷指數(shù)能直觀地反映各種重金屬對污染的貢獻程度,但未考慮不同污染源所引起的背景差別[35]。綜上,在進行土壤重金屬污染風險評價時,應從土壤重金屬的累積程度、生物毒性、遷移特性、污染源特點等方面考慮,選用合適的評價方法,才能準確地評估重金屬在土壤中的污染情況。

3 結 論

(1) 研究區(qū)土壤中的Cd、Cr、As、Pb、Cu和Ni濃度均高于陜西省土壤背景值,Cd和Pb濃度超過標準值,Cd和Pb為主要污染元素。

(2) 研究區(qū)土壤中的Cd、Pb來源于采礦活動和尾礦廢棄物堆積等人為活動影響,Ni和Cu主要受大氣沉降、成土因素等自然因素影響。

(3) 在進行土壤重金屬污染評價時,應綜合考慮土壤重金屬的累積程度、生物毒性、遷移特性等參數(shù)指標,才能更準確地評估重金屬在土壤中的污染情況。