寧夏引黃灌區(qū)排水溝重金屬分布特征及風(fēng)險評估

李延林，鄭 燦，邱小琮，尹 亮，楊子超

(1.寧夏大學(xué)土木與水利工程學(xué)院，銀川 750021；2.寧夏大學(xué)生命科學(xué)學(xué)院，銀川 750021)

近年來，隨著工業(yè)化和城市化進程的加快，大量水域被重金屬污染。人類活動排放的重金屬通過多種方式進入水體，被水中懸浮顆粒物和底部沉積物所吸附，并在其中富集，最終主要儲存于水體沉積物中[1]。重金屬具有不能被微生物降解、潛伏期長、毒性大、難去除、易富集的特征，水體中的重金屬被生物體吸收、富集之后可沿食物鏈傳遞，放大富集進入人體，危害人身體健康[2]。然而在人類活動、降雨、生物擾動、水環(huán)境化學(xué)條件變化等因素的影響下，沉積物中的重金屬又被釋放到水體中，造成水體的二次污染[3]。因此，研究水體中重金屬的含量及分布特征是十分必要的，可為研究區(qū)域的重金屬污染狀況提供第一手的數(shù)據(jù)資料，有利于開展區(qū)域性重金屬污染防治工作。

農(nóng)業(yè)生產(chǎn)中農(nóng)藥、劣質(zhì)化肥等的不合理使用是重金屬的重要來源[5]，殘留在土壤中的農(nóng)藥、化肥等污染物及所含的重金屬在地表水、降水沖刷淋溶等作用下被運輸?shù)胶恿髦幸鹞廴綶6]；化工、鋼鐵等企業(yè)生產(chǎn)過程中產(chǎn)生的“三廢”中含有Hg、Pb、Cu、Zn、As、Cd、Cr等重金屬會對環(huán)境造成污染[7]；城市污水中除含有大量的有機物外，還有Pb、Zn、Cd、Cr等重金屬，污水不處理或處理不完整直接排入河流會增加水體和沉積物中重金屬的濃度[6,8]。張晴雯[9]等分析了寧夏黃灌區(qū)濕地沉積物剖面中TOC、TN、TP、Zn、Cu含量的垂向分布特征，并與濕地周圍農(nóng)田、排水溝和豬場排污口沉積物相應(yīng)層次元素的含量進行了比較。朱嵬[10]等調(diào)查了寧夏黃河流域9個濕地底泥中重金屬(As、Cd、Cr、Cu、Hg、Ni、Pb、Zn)的含量與富集情況，并濕地底泥中重金屬的富集系數(shù)和生態(tài)危害系數(shù)以及各采樣點的生態(tài)危害指數(shù)進行探討。夏修杰[11]等為了解黃河中下游河道淤積泥沙重金屬污染狀況，沿河道采集15 個深層淤積泥沙樣本，對樣本中的Pb、Cu、Cd、Hg、Zn、Ni、Cr共7種重金屬和As 的含量進行了檢測，選擇單因子質(zhì)量指數(shù)法和地累積指數(shù)法開展了評價。

黃河寧夏段位于中衛(wèi)下沿河至石嘴山二水文站之間，沿黃河兩岸地形呈“J”型帶狀分布，該段是黃河上游主要的灌溉用水和水污染嚴(yán)重的河段之一。引黃灌區(qū)涉及青銅峽市、永寧縣、銀川市、賀蘭縣、平羅縣、惠農(nóng)縣、石嘴山市、中衛(wèi)縣、中寧縣、吳忠市、靈武市等11個縣(市)和20多個國營農(nóng)、林、牧場，年均引黃河水70 億m3，其中農(nóng)業(yè)用水占93%～95%，年退入黃河水量25 億m3左右。在黃河寧夏段自流灌區(qū)內(nèi)部，有大小排水溝近200條，主要排水干溝32條，排水總量約25 億m3[12]。近幾年來，引黃灌區(qū)排水溝不僅接納農(nóng)田退水，同時，大量的生活污水及工業(yè)廢水等排入排水溝，使入黃排水溝的水環(huán)境質(zhì)量呈逐年惡化趨勢。胡宜剛[13]等研究了寧蒙引黃灌區(qū)農(nóng)田排水溝渠水質(zhì)特征，發(fā)現(xiàn)寧蒙引黃灌區(qū)農(nóng)田排水污染以氮含量超標(biāo)為主，尤其以農(nóng)溝和干溝為主，而磷濃度較低，基本不造成污染；劉國強[12]等對寧夏引黃灌區(qū)農(nóng)田退水污染現(xiàn)狀進行分析，結(jié)果表明化肥與規(guī)模養(yǎng)殖是灌區(qū)水體TN(總氮)和TP(總磷)的主要污染源，工業(yè)廢水和城市生活污水對水體質(zhì)量仍然具有重要影響；黃仲冬[13]等對寧夏引黃灌區(qū)典型區(qū)域農(nóng)田水質(zhì)特征進行研究，結(jié)果顯示農(nóng)田退水中銨態(tài)氮和總磷的濃度均超過《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》，分別屬于中度污染和重度污染。目前已有的研究主要集中在引黃灌區(qū)排水溝農(nóng)田退水中TN和TP等的研究，對引黃灌區(qū)排水溝重金屬的研究較少?；诖?，本研究在搜集引黃灌區(qū)各排水溝2011-2015年重金屬因子指標(biāo)基礎(chǔ)上，采用主成分分析法分析影響排水溝水質(zhì)的主要重金屬，采用內(nèi)梅羅綜合污染指數(shù)法進行重金屬污染風(fēng)險評估，對寧夏引黃灌區(qū)各排水溝內(nèi)梅羅綜合污染指數(shù)與社會經(jīng)濟指標(biāo)灰關(guān)聯(lián)法分析，確定影響各排水溝的社會經(jīng)濟指標(biāo)，旨在為寧夏引黃灌區(qū)排水溝重金屬污染綜合整治提供一定的理論依據(jù)。

1 材料與方法

1.1 數(shù)據(jù)采集

寧夏引黃灌區(qū)排水溝重金屬檢測數(shù)據(jù)取自《寧夏回族自治區(qū)環(huán)境狀況公報》，包括中衛(wèi)第四排水溝、吳忠南干溝、銀川靈武溝、吳忠清水溝、銀川銀新干溝、銀川四二干溝及石嘴山第三、五排水溝匯合后入黃口7條排水溝。按照引黃灌區(qū)排水溝在寧夏回族自治區(qū)的位置，分別收集中衛(wèi)市、吳忠市、銀川市、石嘴山市的各項社會經(jīng)濟指標(biāo)的原始數(shù)據(jù)，數(shù)據(jù)來源于2011-2015年《寧夏回族自治區(qū)統(tǒng)計年鑒》和《寧夏回族自治區(qū)水資源公報》。

1.2 污染風(fēng)險評價方法

1.2.1 主成分分析法

主成分分析法是一種將多因子納入同一系統(tǒng)進行分析，從而找出關(guān)鍵影響因子的一種統(tǒng)計分析方法。主成分分析法應(yīng)用于水環(huán)境因子分析主要有兩方面：一是構(gòu)建綜合評價指標(biāo)體系，評價水質(zhì)的污染程度，并進行分級；二是評價各單項指標(biāo)的作用，從中選取主要指標(biāo)，確定影響水質(zhì)的主要因子[14]。

1.2.2 灰關(guān)聯(lián)法

灰關(guān)聯(lián)法就是以各排水溝內(nèi)梅羅綜合污染指數(shù)為母序列，影響因子為子序列，計算各排水溝內(nèi)梅羅綜合污染指數(shù)與影響因子的關(guān)聯(lián)度，根據(jù)關(guān)聯(lián)度的大小確定影響因子的作用大小[15,16]。

1.2.3 內(nèi)梅羅綜合污染指數(shù)法

水體中往往是多種重金屬同時存在的，內(nèi)梅羅綜合污染指數(shù)法不僅能夠反映水體重金屬污染現(xiàn)狀，還可以根據(jù)各種重金屬對復(fù)合污染的不同貢獻甄別出主要污染物，是水體重金屬污染評價的常用方法[17,18]。其計算方法為：

單因子污染指數(shù)：

Pi=Ci/Si

(1)

多因子綜合污染指數(shù)：

(2)

式中：Ci為重金屬i的實測濃度；Si為相應(yīng)的水質(zhì)指標(biāo)，寧夏引黃灌區(qū)排水溝的水質(zhì)目標(biāo)為Ⅳ類水，因灌區(qū)的水要排放到黃河，故采用GB 3838-2002 中Ⅲ類水質(zhì)標(biāo)準(zhǔn)作為參比標(biāo)準(zhǔn)[19]；max(Pi)為重金屬單因子污染指數(shù)的最大值；ave(Pi)為各金屬單因子污染指數(shù)的平均值。

單因子污染指數(shù)和多因子污染指數(shù)的評價標(biāo)準(zhǔn)如表1所示[20]。

2 結(jié)果與分析

2.1 寧夏引黃灌區(qū)排水溝重金屬時空分布特征

將寧夏引黃灌區(qū)7條排水溝Hg等8種重金屬濃度的平均值與我國地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB 3838-2002)中的Ⅴ類水質(zhì)標(biāo)準(zhǔn)進行對比，Hg等8種重金屬濃度的平均值均未超過地表水Ⅴ類水質(zhì)量標(biāo)準(zhǔn)。圖1～圖8為寧夏引黃灌區(qū)7條排水溝中Hg等8種重金屬濃度平均值在2011-2015年內(nèi)變化曲線。

圖1 重金屬Hg時空分布Fig.1 Temporal and spatial distribution of heavy metal Hg

圖2 重金屬Pb時空分布Fig.2 Temporal and spatial distribution of heavy metal Pb

圖3 重金屬Se時空分布Fig.3 Temporal and spatial distribution of heavy metal Se

圖4 重金屬Zn時空分布Fig.4 Temporal and spatial distribution of heavy metal Zn

圖5 重金屬As時空分布Fig.5 Temporal and spatial distribution of heavy metal As

圖6 重金屬Cd時空分布Fig.6 Temporal and spatial distribution of heavy metal Cd

圖7 重金屬Cu時空分布Fig.7 Temporal and spatial distribution of heavy metal Cu

中衛(wèi)第四排水溝、銀川銀新干溝、銀川四二干溝和石嘴山第三、五排水溝匯合后入黃口、Hg濃度平均值在2011-2015年總體變化趨勢較平穩(wěn)，無明顯的變化；吳忠南干溝和銀川靈武溝Hg濃度平均值在2011-2015年內(nèi)變化幅度較大，均在2012年達到最高；吳忠清水溝Hg濃度平均值在2011-2013年逐漸升高，并在2013年達到最高，2013-2015年逐漸降低。

中衛(wèi)第四排水溝Pb濃度平均值在2011-2015年保持不變；石嘴山第三、五排水溝匯合后入黃口Pb濃度平均值在2011-2015年內(nèi)變化幅度較大，2011-2012年降低，2012-2013升高，2013-2015年呈逐漸降低趨勢；其余各排水溝Pb濃度平均值在2011-2105年呈逐漸降低趨勢。

中衛(wèi)第四排水溝和吳忠清水溝Se濃度平均值在2011-2015年保持不變；銀川銀新干溝和石嘴山第三、五排水溝匯合后入黃口Se濃度平均值在2001-2015年逐漸降低；其余各排水溝Se濃度平均值在2011-2015年變化趨勢較平穩(wěn)，無明顯的變化。

吳忠清水溝和吳忠南干溝Zn濃度平均值2011-2015年保持不變；銀川靈武溝Zn濃度平均值在2011-2014年變化趨勢平穩(wěn)，2014-2015年突然升高；中衛(wèi)第四排水溝和銀川銀新干溝Zn濃度平均值2011-2015年變化幅度較大，中衛(wèi)第四排水溝Zn濃度平均值在2013年達到最高，銀川銀新干溝Zn濃度平均值2015年最高；銀川四二干溝在2011-2015年總體上呈逐漸升高趨勢；石嘴山第三、五排水溝匯合后入黃口Zn濃度平均值2011-2015年總體上呈下降趨勢。

中衛(wèi)第四排水溝As濃度平均值在2011-2013年變化趨勢平穩(wěn)，2013-2015年呈降低趨勢；吳忠南干溝 As濃度平均值在2011-2013年逐年升高趨勢，2013-2015年逐漸降低；銀川靈武溝和銀川四二干溝As濃度平均值2011-2015年變化幅度較大，均在2013年達到最高；吳忠清水溝As濃度平均值2011-2015年總體上呈上升趨勢；銀川銀新干溝和石嘴山第三、五排水溝匯合后入黃口As濃度平均值2011-2015年呈逐漸降低趨勢。中衛(wèi)第四排水溝、吳忠南干溝和吳忠清水溝Cd濃度平均值在2011-2013年逐年降低，2013-2015年保持不變；石嘴山第三、五排水溝匯合后入黃口Cd濃度平均值在2011-2015年保持不變；其余各排水溝Cd濃度平均值在2011-2015年變化趨勢較平穩(wěn)。

中衛(wèi)第四排水溝、吳忠南干溝和吳忠清水溝Cu濃度平均值在2011-2015年保持不變；銀川靈武溝和銀川四二干溝Cu濃度平均值2011-2014年變化幅度較小，2014-2015年突然上升；銀川銀新干溝Cu濃度平均值在2011-2015年變化幅度大，在2015年達到最高；石嘴山第三、五排水溝匯合后入黃口Cu濃度平均值在2011-2015年變化不大。

中衛(wèi)第四排水溝和吳忠清水溝Cr濃度平均值在2011-2015年保持不變；吳忠南干溝、銀川靈武溝、銀川銀新干溝和銀川四二干溝在2011-2015年變化不大；石嘴山第三、五排水溝匯合后入黃口Cr濃度平均值在2011-2015年變化幅度不大，2014-2015年突然降低。

2.2 基于主成分分析法的寧夏引黃灌區(qū)排水溝重金屬分析

應(yīng)用DPS數(shù)據(jù)處理系統(tǒng)對寧夏引黃灌區(qū)排水溝重金屬進行主成分分析，運用方差最大正交旋轉(zhuǎn)法對因子載荷矩陣進行旋轉(zhuǎn),按照85%的累積方差貢獻率提取主成分,然后選擇旋轉(zhuǎn)后載荷值大于0.6的指標(biāo)作為主要因子進行分析[14,21,22]。各因子的主成分得分值與對應(yīng)的方差貢獻率乘積的和即為各因子的綜合得分[23],計算各重金屬因子的綜合得分，按照分值大小排序，確定影響寧夏引黃灌區(qū)各排水溝的主要重金屬因子。

對各排水溝重金屬進行主成分分析，結(jié)果見表2，計算各排水溝第一主成分重金屬因子的綜合得分，得到影響中衛(wèi)第四排水溝水質(zhì)主要重金屬因子為Cd、Pb；影響吳忠南干溝水質(zhì)主要重金屬因子為Pb、Cu；影響銀川靈武溝水質(zhì)主要重金屬因子為Hg、Zn；影響吳忠清水溝水質(zhì)主要重金屬因子為As、Zn；影響銀川銀新干溝水質(zhì)主要重金屬因子為Cu、As；影響銀川四二干溝水質(zhì)主要重金屬因子為Zn、As；影響石嘴山第三、五排水溝匯合后入黃口水質(zhì)主要重金屬因子為Pb、Cr。

表2 各排水溝主成分分析結(jié)果Tab.2 Principal component analysis results of each drainage ditch

2.3 基于內(nèi)梅羅綜合污染指數(shù)寧夏引黃灌區(qū)排水溝污染評價

2.3.1 內(nèi)梅羅綜合污染指數(shù)

計算寧夏引黃灌區(qū)2011-2015年排水溝中8種重金屬的綜合污染指數(shù)，結(jié)果如圖9所示。從各排水溝的Pn值來看，中衛(wèi)第四排水溝污染等級在2011年處于警戒線，2012-2015年污染整體處于尚清潔水平；吳忠南干溝污染等級在2012年、2013年和2014年達到警戒線，其余各年污染處于安全等級；吳忠清水溝污染等級在2013年達到輕度污染，其余各年污染水平處于尚清潔，但2012年和2014年P(guān)n值分別為0.57和0.5污染等級接近警戒線；銀川靈武溝2011-2015年污染整體處于清潔水平，但2011年、2012年和2013年的Pn值分別為0.5、0.5和0.64接近警戒線等級；銀川銀新干溝和銀川四二干溝在2011-2015年污染整體處于清潔水平；石嘴山第三、五排水溝匯合后入黃口在2011年的Pn=0.58接近警戒線等級，其余各年均污染整體處于清潔水平。

圖9 寧夏引黃灌區(qū)2011-2015年排水溝重金屬的綜合污染指數(shù)Fig.9 Comprehensive pollution index of heavy metals in drainage ditch in ningxia irrigation district of Yellow River diversion from 2011-2015

2.3.2 各排水溝內(nèi)梅羅綜合污染指數(shù)與社會經(jīng)濟指標(biāo)的灰關(guān)聯(lián)分析

選取各市農(nóng)業(yè)耗水總量、農(nóng)用化肥使用量、農(nóng)業(yè)總產(chǎn)值、工業(yè)總產(chǎn)值、工業(yè)廢水排放總量和城鎮(zhèn)生活污水排放總量社會經(jīng)濟指標(biāo)，用DPS進行各排水溝綜合污染指數(shù)與社會經(jīng)濟指標(biāo)灰關(guān)聯(lián)分析，分析結(jié)果見表3。

影響中衛(wèi)第四排水溝水質(zhì)的主要重金屬是Cd、Cr，與內(nèi)梅羅綜合污染指數(shù)關(guān)聯(lián)度高的是工業(yè)廢水排放總量、農(nóng)業(yè)總產(chǎn)值、城鎮(zhèn)生活污水排放量，污染等級在2011年處于警戒線，2012-2015年污染整體處于尚清潔水平，主要由于對工業(yè)廢水排放總量進行了控制，強化對重點行業(yè)、重點企業(yè)主要污染物控制；影響吳忠南干溝水質(zhì)的主要重金屬是Pb、Cu，與內(nèi)梅羅綜合污染指數(shù)關(guān)聯(lián)度高的是城鎮(zhèn)生活污水排放量、工業(yè)總產(chǎn)值、農(nóng)業(yè)耗水量，在2011-2013年水質(zhì)惡化，主要由于排入?yún)侵?/p>

表3 各排水溝綜合污染指數(shù)主要影響因子Tab.3 Main influence factors of comprehensive pollution index of each drainage ditch

南干溝的工業(yè)廢水和城鎮(zhèn)生活污水增多，工業(yè)廢水排放總量和城鎮(zhèn)生活污水排放量增加，導(dǎo)致吳忠南干溝水質(zhì)惡化；影響吳忠清水溝主要重金屬是As、Zn，內(nèi)梅羅綜合污染指數(shù)關(guān)聯(lián)度高的是工業(yè)廢水排放總量、農(nóng)用化肥施用量、城鎮(zhèn)生活污水排放量，在2011-2014年水質(zhì)惡化，主要由于排入?yún)侵夷细蓽虾蛥侵仪逅疁系墓I(yè)廢水和城鎮(zhèn)生活污水增多，工業(yè)廢水排放總量和城鎮(zhèn)生活污水排放量增加，導(dǎo)致吳忠清水溝水質(zhì)惡化；影響銀川靈武溝水質(zhì)的主要重金屬是Hg、Zn，內(nèi)梅羅綜合污染指數(shù)關(guān)聯(lián)度高的是農(nóng)業(yè)耗水量、工業(yè)廢水排放總量、農(nóng)用化肥施用；銀川銀新干溝水質(zhì)的主要重金屬是Cu、As，內(nèi)梅羅綜合污染指數(shù)關(guān)聯(lián)度高的是農(nóng)業(yè)總產(chǎn)值、城鎮(zhèn)生活污水排放量、農(nóng)業(yè)耗水量；銀川四二干溝水質(zhì)的主要重金屬是Zn、As，內(nèi)梅羅綜合污染指數(shù)關(guān)聯(lián)度高的是工業(yè)廢水排放總量、農(nóng)用化肥施用量、農(nóng)業(yè)耗水量；銀川靈武溝在2011-2013年水質(zhì)相對較差、銀川四二干溝和銀川銀新干溝在2011年水質(zhì)較差，主要是由于排入銀川靈武溝、銀川四二干溝和銀川銀新干溝的工業(yè)廢水總量和農(nóng)田退水量增多，隨后水質(zhì)逐漸改善是由于減少了排入排水溝的工業(yè)廢水量、城鎮(zhèn)生活污水量和化肥使用量，使農(nóng)業(yè)面源污染和工業(yè)點源污染減輕；影響石嘴山第三、五排水溝匯合后入黃口水質(zhì)的主要重金屬是Pb、Cr，與內(nèi)梅羅綜合污染指數(shù)關(guān)聯(lián)度高的主要影響因子依次為農(nóng)業(yè)耗水量、工業(yè)總產(chǎn)值、城鎮(zhèn)生活污水排放量，2011年污染等級接近警戒線，各年均污染整體處于清潔水平，2011-2015年石嘴山第三、五排水溝匯合后入黃口水質(zhì)改善是由于對工業(yè)廢水排放總量進行了控制，強化對重點行業(yè)、重點企業(yè)主要污染物控制，同時減少了農(nóng)田退水量和化肥使用量。

3 結(jié) 語

(1)2011-2015年寧夏引黃灌區(qū)7條排水溝中8種重金屬含量的平均值均低于地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB3838-2002)中Ⅴ類水質(zhì)標(biāo)準(zhǔn)，7條排水溝重金屬濃度平均值整體上在2011-2015年內(nèi)變化幅度不大。

(2)運用主成分分析法研究確定影響寧夏引黃灌區(qū)各排水溝水質(zhì)的主要重金屬，分析表明，影響中衛(wèi)第四排水溝水質(zhì)的主要重金屬是Cd、Cr；影響吳忠南干溝水質(zhì)的主要重金屬是Pb、Cu；影響銀川靈武溝水質(zhì)的主要重金屬是Hg、Zn；影響吳忠清水溝水質(zhì)的主要重金屬是As、Zn；影響銀川銀新干溝水質(zhì)的主要重金屬是Cu、As；影響銀川四二干溝水質(zhì)的主要重金屬是Zn、As；影響石嘴山第三、五排水溝匯合后入黃口水質(zhì)的主要重金屬是Pb、Cr。

(3)用內(nèi)梅羅綜合污染指數(shù)法評價寧夏引黃灌區(qū)排水溝污染評，結(jié)果表明中衛(wèi)第四排水溝污染等級在2011年處于警戒線，2012-2015年污染整體處于尚清潔水平；吳忠南干溝污染等級在2012年、2013年和2014年達到警戒線，其余各年污染處于安全等級；吳忠清水溝污染等級在2013年達到輕度污染，其余各年污染水平處于尚清潔，但2012年和2014年P(guān)n值分別為0.57和0.5污染等級接近警戒線；銀川靈武溝2011-2015年污染整體處于清潔水平，但2011年、2012年和2013年的Pn值分別為0.5、0.5和0.64接近警戒線等級；銀川銀新干溝和銀川四二干溝在2011-2015年污染整體處于清潔水平；石嘴山第三、五排水溝匯合后入黃口在2011年的Pn=0.58接近警戒線等級，其余各年均污染整體處于清潔水平。