謝阿廷，程 艷，趙晨曦，胡 霞，張 健，王 悅

(1.新疆農(nóng)業(yè)大學(xué)草業(yè)與環(huán)境科學(xué)學(xué)院，烏魯木齊 830052；2.新疆環(huán)境污染控制與風(fēng)險預(yù)警重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830011；3.新疆環(huán)境保護(hù)科學(xué)研究院，烏魯木齊 830011；4.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院，烏魯木齊 830046)

0 引 言

【研究意義】水質(zhì)評價是對水環(huán)境質(zhì)量狀況的定性分析與定量評定、是關(guān)于水環(huán)境質(zhì)量現(xiàn)狀和變化的理論、方法及其應(yīng)用技術(shù)的科學(xué)，對水體的污染程度進(jìn)行等級劃分和確定污染類型的過程，通過對水質(zhì)監(jiān)測數(shù)據(jù)的合理評價，才能制定科學(xué)的整治規(guī)劃采取有效的措施[1]。隨著水資源的匱乏，對水資源、水環(huán)境的優(yōu)劣也日益關(guān)注，水質(zhì)評價的重要性也凸顯出來了，對于流域附近的居民，河流水質(zhì)的狀況影響著生活環(huán)境的舒適，對河流水質(zhì)的評價也具有一定的實(shí)際意義。近年來對伊犁河流域的研究，大部分集中于伊犁河生物群落結(jié)構(gòu)的分析研究，鮮有對伊犁河流域的綜合水質(zhì)進(jìn)行評價，而對該流域水質(zhì)的評價具有重要意義，為水環(huán)境污染防控提供科學(xué)依據(jù)和技術(shù)支撐?！厩叭搜芯窟M(jìn)展】有關(guān)水質(zhì)評價的方法很多，主要分為單因子評價法[2]、綜合污染指數(shù)法、模糊評價法、灰色評價法、人工神經(jīng)網(wǎng)絡(luò)法和主成分分析法[3]。單因子評價法方法簡單，是以最差水質(zhì)指標(biāo)所屬類別作為綜合水質(zhì)類別，但卻無法判斷水體的整體污染情況；污染指數(shù)法無法對斷面水質(zhì)類別進(jìn)行定性判斷，卻能反映水體的相對污染程度。其他綜合評價方法是基于水質(zhì)分析是一個多因素、多指標(biāo)的復(fù)雜過程，從而采用一些合理的數(shù)學(xué)模式，對這些諸多因素進(jìn)行綜合，以便做出合理評價或分析，其結(jié)果往往反映的是對象的綜合性特征，及其個體在整體中的綜合表現(xiàn)，是對單因子分析、綜合指數(shù)分析等方法的合理補(bǔ)充。因此，采用多方法進(jìn)行水質(zhì)評價是全面判別水體污染狀況的必然趨勢。例如，李曉航等[4]采用單因子污染指數(shù)法、內(nèi)梅羅污染指數(shù)法和綜合污染指數(shù)法對艾比湖流域水質(zhì)進(jìn)行了水質(zhì)評價，并探討了艾比湖流域的總體污染狀況。上述綜合性方法中，主成分分析法具有從大而復(fù)雜的數(shù)據(jù)庫獲得更好的關(guān)于水質(zhì)和分析的信息分析，并能有效地識別染源[5-6]。伊元榮等[7]采用多元統(tǒng)計(jì)主成分分析方法對伊犁河進(jìn)行水質(zhì)評價，結(jié)果表明，主成分分析方法能真實(shí)、準(zhǔn)確分析出城市河流各水質(zhì)指標(biāo)的內(nèi)在關(guān)系。周豐等[8]提出了一種基于多元統(tǒng)計(jì)分析和RBFNNs的水質(zhì)評價方法，對各組樣本進(jìn)行水質(zhì)評價，評價結(jié)果與傳統(tǒng)單項(xiàng)指數(shù)法的結(jié)果基本一致。【本研究切入點(diǎn)】不同的方法有各自不同的優(yōu)缺點(diǎn)，其評價結(jié)果所表達(dá)的水環(huán)境特征也不盡相同。目前，研究所涉及到的河段僅有2個常規(guī)水質(zhì)監(jiān)測點(diǎn)，而河段內(nèi)城鎮(zhèn)生活、工業(yè)、農(nóng)業(yè)、養(yǎng)殖等各類污染源分布較為密集，該河段夏季的水污染較為復(fù)雜和顯著，現(xiàn)有監(jiān)測點(diǎn)位和資料遠(yuǎn)不能合理反映區(qū)域水質(zhì)的現(xiàn)狀特征，因此，選取多個有效斷面及水質(zhì)評價方法的聯(lián)用對伊犁河水污染管理與防控具有重要意義。【擬解決的關(guān)鍵問題】通過在夏季7月對研究河段進(jìn)行水質(zhì)加密監(jiān)測并運(yùn)用單因子評價法、內(nèi)梅羅指數(shù)法、重金屬水質(zhì)質(zhì)量指數(shù)法、主成分分析法，對伊犁河河段水環(huán)境特征進(jìn)行評價分析，全面掌握研究河段的水環(huán)境質(zhì)量現(xiàn)狀特征，為水環(huán)境污染防控提供科學(xué)依據(jù)和技術(shù)支撐。

1 材料與方法

1.1 材 料

研究河段位于伊犁河伊寧市段，河段總長約40 km，自東向西流，水深較淺，水流湍急，水體自凈能力較強(qiáng)，河道南北側(cè)均有河道或排水口匯入，南岸以休閑生活為主，多為漁業(yè)養(yǎng)殖、農(nóng)業(yè)灌溉和景觀生活用水；北岸則以城市、工業(yè)園區(qū)為主，工業(yè)用水多由城西污水處理廠集中處理后排入河內(nèi)。為較全面評價研究河段水環(huán)境質(zhì)量狀況，考慮到河道兩側(cè)污染來源的差異性，分別于河段南北岸各選取5個斷面，北岸斷面為A1、A4、A8、A9、A10、南岸斷面為A2、A3、A5、A6、A7。

2015年7月20～25日，按照《水質(zhì)采樣技術(shù)指導(dǎo)》(HJ 494-2009)要求，現(xiàn)場用PH-B-1型酸度計(jì)和YSI-85-25型便攜式多功能水質(zhì)分析儀對pH值、水溫、溶氧、電導(dǎo)率、進(jìn)行現(xiàn)場測定；采集0～50 cm的上層水，水樣封存于不同材質(zhì)的瓶子中，依次編號，加適量保護(hù)劑后，低溫保存(4℃)。圖1

圖1 采樣點(diǎn)位置
Fig.1 Location of sampling points

1.2 方 法

1.2.1 水樣指標(biāo)的選擇

根據(jù)西北地區(qū)水體污染情況和有關(guān)水質(zhì)評價的文獻(xiàn)[9-11]，選取具有該河段特征污染性質(zhì)的氨氮、總氮、總磷、高錳酸鹽指數(shù)、化學(xué)需氧量、陰離子表面活性劑、溶解氧、氯化物、以及重金屬六價鉻、銅、鉛、鋅、鎘、錳、砷、鐵等指標(biāo)為該河流段的水質(zhì)評價指標(biāo)。

1.2.2 評價標(biāo)準(zhǔn)的確定

該河段水質(zhì)目標(biāo)為《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB3838-2002)中的Ⅱ類標(biāo)準(zhǔn)，則以此標(biāo)準(zhǔn)作為研究河段水質(zhì)評價的依據(jù)標(biāo)準(zhǔn)。

1.2.3 評價方法

1.2.3.1 單因子評價法

單因子評價法是最直接、簡單明了的水質(zhì)評價方法，是《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB3838-2002)中推薦的水質(zhì)評價方法之一。其方法為在所有參與綜合水質(zhì)評價的水質(zhì)指標(biāo)中選擇水質(zhì)最差的單項(xiàng)指標(biāo)所屬類別來確定所屬水域綜合水質(zhì)類別。

計(jì)算公式如下：

式中，Sij：單項(xiàng)水質(zhì)評價因子i在第j取樣點(diǎn)的標(biāo)準(zhǔn)指數(shù)；Cij：水質(zhì)評價因子i在第j取樣點(diǎn)的濃度，(mg/L)；Csi：水質(zhì)評價因子i的評價標(biāo)準(zhǔn)，(mg/L)。

1.2.3.2 內(nèi)梅羅指數(shù)

內(nèi)梅羅指數(shù)是兼顧極值和平均值的一種計(jì)權(quán)型多因子環(huán)境質(zhì)量指數(shù)評價法，能夠全面反映水體污染現(xiàn)狀并甄別出主要污染物，是水質(zhì)評價的常用方法之一。計(jì)算表達(dá)式為：

式中，Ci為項(xiàng)目i的實(shí)測濃度；Cs為項(xiàng)目i的標(biāo)準(zhǔn)值；max(Pi)為項(xiàng)目單因子污染指數(shù)的最大值；ave(Pi)為各項(xiàng)目單因子污染指數(shù)的平均值。內(nèi)梅羅指數(shù)的評價標(biāo)準(zhǔn)[11]：當(dāng)PN≤0.7時，表示無污染；當(dāng)0.72時，表示強(qiáng)污染。

1.2.3.3 重金屬污染水質(zhì)質(zhì)量指數(shù)

以1.2.2中的Ⅱ類水質(zhì)標(biāo)準(zhǔn)作為該河流域水中重金屬污染評價參照標(biāo)準(zhǔn)，使用水質(zhì)質(zhì)量指數(shù)法對水中的重金屬污染進(jìn)行評價：

式中，Ci為重金屬i的實(shí)測濃度；Qi為重金屬i對應(yīng)的水質(zhì)標(biāo)準(zhǔn)；Ai為重金屬i的污染指數(shù)；WQI為水中重金屬的水質(zhì)質(zhì)量指數(shù)。

根據(jù)WQI的值，將河水中的重金屬污染情況分為4個等級[12]，當(dāng)WQI≤1時，表示無污染；當(dāng)13時，表示強(qiáng)污染。

1.2.3.4 主成分分析法

是英國生物學(xué)家Kart Pearson在20世紀(jì)初提出來的。能夠綜合評價指標(biāo)，分析多種水污染因子在各斷面的污染程度，并能有效的識別污染源。即能減少參與數(shù)據(jù)分析的變量個數(shù)，同時也不會造成統(tǒng)計(jì)信息的大量浪費(fèi)和丟失。在進(jìn)行因子分析之前先對數(shù)據(jù)進(jìn)行KMO和巴特利特檢驗(yàn)，KMO檢驗(yàn)是為了看數(shù)據(jù)是否適合因子分析，結(jié)果為KMO值為0.624，一般認(rèn)為在0.6左右適合做因子分析，Sig.值為0小于0.001，表明監(jiān)測數(shù)據(jù)適合做因子分析。計(jì)算步驟如下：

(1)將原始數(shù)據(jù)xij標(biāo)椎化，以后的變量為z1，z2，……zp，標(biāo)準(zhǔn)化數(shù)據(jù)矩陣，以消除變量間在數(shù)量級和量綱上的不同。

(2)求標(biāo)準(zhǔn)化數(shù)據(jù)的相關(guān)矩陣rij及相關(guān)矩陣的特征值λ和特征向量。

(3)計(jì)算方差貢獻(xiàn)率與累積方差貢獻(xiàn)率。

(4)確定因子：設(shè)F1，F(xiàn)2，....，F(xiàn)p為p個主成分，其中前m個因子包含的數(shù)據(jù)信息總量(及累計(jì)貢獻(xiàn)率)不低于85%，從而選取主成分個數(shù)。

(5)用原指標(biāo)的線性組合來求各主成分得分F，以各主分的方差貢獻(xiàn)率為權(quán)，由各因子的線性組合得到綜合評價指標(biāo)函數(shù)。

2 結(jié)果與分析

2.1 重金屬污染水質(zhì)質(zhì)量指數(shù)

通過對重金屬指標(biāo)進(jìn)行評價，各采樣點(diǎn)位水體中重金屬的水質(zhì)質(zhì)量指數(shù)值較低，WQI范圍為：0.001～0.63，表明研究河段不存在重金屬污染。表1

表1 重金屬污染水質(zhì)質(zhì)量指數(shù)
Table 1 Ater quality index of heavy metal pollution

指標(biāo)Index六價鉻Cr6+銅Cu鋅Zn鎘Cd錳Mn砷AsWQI0.120.0010.0010.50.630.06

2.2 單因子評價法

研究表明，氨氮除A8斷面處于Ⅱ類水體外，其它斷面均在Ⅰ類水平?？偟狝3、A4屬于Ⅱ類水體，A1、A2、A10、屬于Ⅲ類水體，A7、A9為Ⅳ類水體，A8、A5為Ⅴ類水體，A6為劣Ⅴ類水體，90%的斷面未達(dá)到標(biāo)準(zhǔn)要求，總氮污染較為突出?？偭讓儆冖蝾愃w的斷面有A3、A7、A8、A10，屬于Ⅲ類水體的有A1、A2、A4、A6、A9，A5達(dá)到Ⅳ類水體，存在超標(biāo)現(xiàn)象。高錳酸鹽除A5斷面Ⅱ類水體外，其它斷面均為Ⅰ類水體；化學(xué)需氧量A3、A4、A7屬于Ⅲ類水體，A1、A2、A9、A8、A6、A10屬于Ⅳ類水體，A5屬于Ⅴ類水體，超標(biāo)較為普遍。陰離子表面活性劑均處于Ⅰ類水平。

主要污染指標(biāo)的空間分布情況為：總氮處于中度污染，其最大值(4.10)出現(xiàn)在采樣點(diǎn)A6，為劣Ⅴ類水體；總磷處于輕度污染，其最大值(2.40)出現(xiàn)在采樣點(diǎn)A5，為Ⅳ類水體；化學(xué)需氧量處于輕度污染，其最大值(2.12)出現(xiàn)在采樣點(diǎn)A5，為Ⅴ類水體。圖2

圖2 單因子指數(shù)
Fig.2 Single factor index

2.3 內(nèi)梅羅指數(shù)

研究表明，河段水體綜合水平良好，整體表現(xiàn)為耗氧性有機(jī)污染；水體中主要污染物為總氮、總磷、化學(xué)需氧量，其中總氮為強(qiáng)污染物。研究河段總氮、總磷、化學(xué)需氧量的污染與農(nóng)業(yè)用水中作物沒吸收的化肥和農(nóng)藥的殘留及居民污水和牲畜糞便等流入流域的污染源息息相關(guān)，也與岸邊植被樹木凋落的腐殖物質(zhì)有關(guān)。表2

表2 評價因子內(nèi)梅羅指數(shù)
Table 2 Nemerow index evaluation factor

評價因子Factor氨氮NH4+總氮TN總磷TP化學(xué)需氧量COD陰離子表面活性劑An-ionic surfactan溶氧Do氯化物Cl高錳酸鹽指數(shù)KMno4Pn0.562 53.276 91.881 81.824 20.348 70.731 30.088 10.482 3污染程度無污染強(qiáng)污染中污染中污染無污染低污染無污染無污染

內(nèi)梅羅指數(shù)的評價結(jié)果與單因子評價法結(jié)果具有一致性。研究河段總氮為強(qiáng)污染，總磷、化學(xué)需氧量為低污染。表明研究河段水體以營養(yǎng)鹽污染和有機(jī)污染為主。

2.4 主成分分析

為進(jìn)一步明確影響不同監(jiān)測斷面的主要污染指標(biāo)，及其主要污染因子的空間分布綜合特征，用主成分分析對10個斷面的8個水質(zhì)指標(biāo)進(jìn)行因子分析，提取特征值大于1的若干個主成分，同時得到各組分貢獻(xiàn)率與各指標(biāo)的關(guān)系。依據(jù)Kaiser-Harris準(zhǔn)則提取特征值大于1的主成分，所以只選擇了前兩個主成分，第一個主成分的方差貢獻(xiàn)率是71.269%，前兩個主成分方差占所有主成分方差的87.849%，由此可見，選前兩個主成分足夠替代原來的變量，幾乎涵蓋了原變量的全部信息。其中第一個主成分在氨氮、高錳酸鹽指數(shù)、化學(xué)需氧量變量上有較大的載荷，所以其反映的是這些變量的信息，第二個因子反映的是總氮的變量信息。表3

根據(jù)監(jiān)測指標(biāo)數(shù)據(jù)的標(biāo)椎化與成分得分系數(shù)，由此直接寫出各公因子的表達(dá)式

F1=0.252×Z氨氮-0.166×Z總氮+0.182×Z總磷+0.089×Z化學(xué)需氧量-0.311×Z溶氧+0.118×Z高錳酸鹽指數(shù)-0.174×Z氯化物+0.251×陰離子表面活性劑

F2=-0.095×Z氨氮+0.489×Z總氮+0.002×Z總磷+0.12×Z化學(xué)需氧量+0.215×Z溶氧+0.122×Z高錳酸鹽指數(shù)+0.511×Z氯化物-0.093×Z陰離子表面活性劑

F=[λ1/(λ1+λ2)]F1[λ2/(λ1+λ2)]F2.

表3 F綜合得分
Table 3 F Comprehensive Score

水樣編號Sample numberF1得分F1 scoreF1排名F1 rankingF2得分F2 scoreF2排名F2 ranking綜合得分FComprehensive scoreF排名F rankingA1-0.1215-0.4547-0.1835A2-0.1286-0.5408-0.2056A3-0.2678-0.97210-0.4008A4-0.0074-0.8099-0.1584A52.16011.54212.0431A60.03931.06020.2323A7-0.2547-0.2226-0.2487A8-0.945100.4263-0.68710A90.05020.18640.6932A10-0.5279-0.2185-0.4699

研究表明(1)研究河段整體表現(xiàn)為河段北岸污染高于南岸。(2)由于水體流速的作用，將上游污染物帶到下游水域，上游水質(zhì)好于下游。(3)A5斷面F1得分較高，說明該斷面受氨氮、高錳酸鹽指數(shù)、化學(xué)需氧量影響較大，該斷面主要是受農(nóng)田退水及生活污水的影響，該斷面上游為污水處理廠及工業(yè)集中區(qū)。A5、A6斷面F2得分較高，表示該斷面受總氮影響較大?？偟奈廴九c農(nóng)業(yè)污水和居民污水牲畜糞便有關(guān)，這兩個斷面均在污水處理廠附近。A5斷面的綜合得分第一，且遠(yuǎn)甚于其它斷面，污染較其它斷面來說比較嚴(yán)重，主要的影響因子分別是總氮、氨氮、高錳酸鹽指數(shù)、化學(xué)需氧量，流域污染表現(xiàn)為有機(jī)污染。

主成分分析法的結(jié)果表明，研究河段的A5、A9、A6斷面水質(zhì)較差，主要的影響因子分別是總氮、氨氮、高錳酸鹽指數(shù)、化學(xué)需氧量，表明研究河段水體以有機(jī)污染為主，河段的上游區(qū)域水質(zhì)好于下游區(qū)域。

3 討 論

主成分分析法目前是運(yùn)用較多的評價方法之一，操作簡單，減少主觀誤差這與廉鐵輝等[13](2017)對主成分分析法的評價結(jié)論一致。而除了采用主成分分析法之外，更結(jié)合單因子評價法和綜合污染指數(shù)法聯(lián)合評價河流流域，使結(jié)果更具有準(zhǔn)確性。

孫博等(2016)[14]和王洪濤等(2016)[15]采用主成分分析法對河流沉積物中的重金屬污染進(jìn)行評價，得到污染強(qiáng)度的結(jié)論，而研究采用主成分分析法對水質(zhì)指標(biāo)進(jìn)行評價，因?yàn)樵摿饔虻闹亟饘贌o污染。得到伊犁河流域的氮元素污染最為嚴(yán)重，總氮的污染程度是強(qiáng)污染，但該流域卻并不只是受氮元素的單一污染，除了氮元素污染外，還伴隨著總磷，化學(xué)需氧量等其他污染物質(zhì)，流域變現(xiàn)為典型的有機(jī)污染的結(jié)論。因此，采用多方法對河流進(jìn)行水質(zhì)評價是全面判別水體污染狀況的必然趨勢。

4 結(jié) 論

4.1 通過單因子評價法和內(nèi)梅羅指數(shù)法可以明確該河流段的主要污染因子為總氮、總磷、化學(xué)需氧量，其中總氮最大超標(biāo)倍數(shù)為4.36，超標(biāo)率達(dá)93.8%；總磷最大超標(biāo)倍數(shù)為6.37，超標(biāo)率達(dá)75%；化學(xué)需氧量最大超標(biāo)倍數(shù)為3.2，超標(biāo)率達(dá)100%。表明研究河段水體以營養(yǎng)鹽污染和有機(jī)污染為主，河段的中下游區(qū)域?yàn)槲廴据^為嚴(yán)重的區(qū)域。

4.2 通過重金屬污染水質(zhì)質(zhì)量指數(shù)法對該流域重金屬污染程度進(jìn)行評價，研究河段水體中6種重金屬的水質(zhì)質(zhì)量指數(shù)數(shù)值較低，各指標(biāo)均能達(dá)到Ⅱ類水質(zhì)標(biāo)準(zhǔn)；WQI范圍為：0.001～0.63，整體污染水平為無污染。

4.3 綜合F得分中A5為最大值(2.043)，由于A5斷面水質(zhì)受工業(yè)污水與農(nóng)業(yè)污水的影響，該斷面水質(zhì)在該流域中水質(zhì)最差，A6斷面F綜合得分為0.232，主要與農(nóng)業(yè)污水和居民污水有關(guān)；整體水域表現(xiàn)為河段北面整體污染高于南面，水體污染原因與生活污水和農(nóng)業(yè)活動及工業(yè)污水有關(guān)。