陳利粉,陳 龍,宋 盈,馬建茹,楊衛(wèi)芳

(河南省安陽生態(tài)環(huán)境監(jiān)測中心,河南 安陽 455000)

0 引言

彰武水庫是位于海河流域洹河干流上的一座中型水庫，其控制流域面積為120 km2，總庫容7 063萬m3，興利庫容1 703萬m3[1]。彰武水庫水源主要來自于南海泉，并與上游的大型水庫小南海水庫聯(lián)合調度運用，流入洹河，橫穿安陽市境內，最終匯入衛(wèi)河，其主要作用是保障安鋼、大唐電廠、安化等工業(yè)用水及下游18萬hm2的農田引蓄灌溉、水產養(yǎng)殖，同時也是安陽市重要的防洪屏障[2-3]。作為安陽市主要的工、農業(yè)供水水源地和地下水補給源之一，通過水質監(jiān)測數據，分析彰武水庫水質狀況對于庫區(qū)周邊及下游工、農業(yè)生產以及開展庫區(qū)水污染防治工作具有十分重要的意義，為此，該水庫的水質評價也受到了廣泛關注。郭學士[3]采用Daniel檢驗法對彰武水庫和小南海水庫的關鍵因子進行了趨勢分析；唐敏等[2]采用綜合營養(yǎng)狀態(tài)指數法(TLI)對彰武水庫的營養(yǎng)化狀態(tài)進行評判；李中原等[4]通過單因子水質標識指數法、綜合污染指數法、有機污染指數法對彰武水庫水環(huán)境狀況進行了評估；張培等[5]采用輸出系數法和等標污染負荷法，對彰武水庫流域18個鄉(xiāng)鎮(zhèn)農村生活面源污染物進行了評估研究。但這些研究方法多采用單因子評價法、綜合指數法和水質綜合營養(yǎng)狀態(tài)指數法，而應用主成分分析法、加拿大水質指數法(CCME WQI)對彰武水庫水質進行評價的研究卻少見報道。

主成分分析法是一種多元統(tǒng)計技術，利用降維的方式對數據進行處理，將較多的指標轉化為少數幾個綜合指標，同時又設法保留原始數據中存在的關系，而且具有優(yōu)于原變量的性質，使原問題更易于分析，被廣泛的應用于水環(huán)境質量評價中[6-14]。CCME WQI是目前國際上應用較為廣泛的水質評價方法之一，該方法同時考慮了超標范圍F1、超標頻率F2、超標幅度F3等3個因素對水質狀況的影響，與其他形式的水質指數相比，其具有基準選擇的靈活性和對缺失數據的較好的容忍度等優(yōu)勢，因此，近年來，國內外較多學者從不同角度對CCME WQI進行了研究并應用于水質評價中[15-18]。

本研究以彰武水庫出庫口作為研究對象，利用Excel 2020軟件、Origin 2021軟件、SPSS 24.0軟件，綜合評價斷面水質狀態(tài)，分析2011—2020年年際及每年豐、平、枯水期不同時期水質指標的變化特征，應用主成分分析法、綜合營養(yǎng)狀態(tài)指數法(TLI)和加拿大水質指數法(CCME WQI)3種不同評價方法對彰武水庫水質進行評價的差異性，并采用主成分分析方法，提取水質污染因子并分析污染源，擬為安陽市彰武水庫污染防治工作提供參考。

1 材料與方法

1.1 水質指標選擇及測定方法

水樣的采集及分析工作均由河南省安陽生態(tài)環(huán)境監(jiān)測中心完成，監(jiān)測點位為彰武水庫出庫口，研究收集了2011—2020年的水質監(jiān)測數據。為保證監(jiān)測數據的穩(wěn)定性和可比性，于全年不同水期選取代表月份(9月代表豐水期，1月代表平水期，5月代表枯水期)作為研究時段，每月月初采集1次樣品，在采樣點下0.5 m深位置處采集水樣進行分析。監(jiān)測指標包括水溫(T)、pH值、溶解氧(DO)、透明度(SD)、高錳酸鹽指數(CODMn)、生化需氧量(BOD)、氨氮(NH3-N)、總磷(TP)、總氮(TN)和葉綠素a(chla)，其中pH值、溶解氧、透明度在現場進行測定，其余指標均選用《地表水環(huán)境質量標準》(GB 3838-2002)[19]中推薦的方法進行監(jiān)測，根據Ⅲ類標準進行水質評價。

1.2 數據分析方法

湖庫綜合營養(yǎng)狀態(tài)指數(TLI)評價是指通過與湖泊營養(yǎng)狀態(tài)相關的一系列指標，以及這些指標之間的相互關系，對湖泊的營養(yǎng)狀態(tài)進行評價的一種方法，通常選取葉綠素a(chla)、總磷(TP)、總氮(TN)、透明度(SD)、高錳酸鹽指數(CODMn)等指標進行評價[20]。

(1)

(2)

式中：TLIj、Wj、rij分別代表第j種參數的營養(yǎng)狀態(tài)指數、營養(yǎng)狀態(tài)指數的相關權重及其與基準參數chla的相關系數。通常采用0～100的一系列數字把水庫營養(yǎng)狀態(tài)分為5級，包括貧營養(yǎng)(TLI<30)、中營養(yǎng)(30≤TLL≤50)、輕度富營養(yǎng)(5070)。

主成分分析(PCA)是一種利用降維的方式，對批量數據進行處理，將較多的指標在保留原始數據之間關系的同時轉化為少數幾個綜合指標，可避免評價因子選擇的隨意性，使評價結果更為客觀。應用SPSS軟件對各水質指標進行降維處理，提取特征值>1的主成分因子，將主成分載荷矩陣中的系數與標準化后的監(jiān)測數據對應相乘后求和，得到各主成分的表達式Fi，再將這些主成分的值乘以各自對應特征值的方差貢獻率后，再相加，得到主成分綜合得分F，進一步計算各因子得分和綜合得分[21]。

CCME WQI是由加拿大環(huán)境部長理事會提出的一種水質評價方法，計算過程如下：

(3)

(4)

(5)

(6)

(7)

式中F1、F2、F3分別代表超標范圍、超標頻率和超標幅度，CCME WQI取值介于0～100之間，根據其數值大小，可將水體的水質狀況分為5個等級，見表1[18]。

表1 CCME WQI水質等級劃分及說明Tab.1 Classification and explaination of Canadian water quality index

2 結果與分析

2.1 水質理化指標描述性統(tǒng)計

由表2可見：彰武水庫水溫呈季節(jié)性變化，pH值在7.25～8.16之間，水體全年呈弱堿性；與平水期相比，除葉綠素a外，各監(jiān)測因子均值豐水期均小于平水期，其中溶解氧、生化需氧量與總氮最大值與均值明顯下降，氨氮最大值明顯下降；與平水期相比，除氨氮外，各監(jiān)測因子均值枯水期均小于平水期，其中溶解氧最大值與均值明顯下降，氨氮均值與標準差明顯上升；溶解氧的最小值出現在豐水期，而葉綠素a的均值和最大值均出現在豐水期，這可能是由于降雨的增多，藻類生長在水體表面形成遮光阻氣層，使水體光合作用受阻，影響水體氧平衡；氨氮均值在枯水期最大，平水期次之，豐水期最小，水體氨氮質量濃度介于Ⅰ類水和Ⅱ類水之間；各水期總氮均值均大于地表水Ⅴ類標準，而總磷濃度均較低，且變化趨勢不明顯，可見，彰武水庫屬于磷限制水體，而控氮則是其主要任務之一[22]；豐水期生化需氧量與化學需氧量的比值為0.209(小于0.3)，平水期與枯水期均大于0.3(平水期0.311，枯水期0.339)，表明彰武水庫在豐水期可生化性不高，有機污染與水體富營養(yǎng)化較為嚴重[23]。

表2 彰武水庫水質理化指標描述性統(tǒng)計Tab.2 Descriptive statistics of water quality index in Zhangwu Reservoir

2.2 各指標之間相關性分析

彰武水庫主要污染指標的Spearman相關性見圖1：生化需氧量與透明度、高錳酸鹽指數、總氮、葉綠素a呈顯著正相關，高錳酸鹽指數與生化需氧量、總氮、總磷、葉綠素a呈顯著正相關，與氨氮呈顯著負相關，透明度與水溫呈顯著負相關；氨氮和總磷作為浮游植物生長的必需營養(yǎng)物質，并未表現出顯著相關性，可能是水體氮主要由硝酸鹽氮組成，而氨氮濃度過小，水體總磷絕對濃度較低，但N/P比過高，進一步證明了彰武水庫同我國較多水庫相類似，屬于磷限值水體，而控氮則是其主要任務之一[22]。

圖1 各監(jiān)測指標相關性分析Fig.1 Correlation analysis of water quality index in Zhangwu Reservoir注:*P≤0.05。

2.3 主要影響因子分析

2.3.1 2011—2020年水庫水質主要影響因子分析

應用SPSS軟件，根據特征值的大小提取主成分，特征值分析見表3。根據主成分選取原則，特征值>1，提取4個主成分，其方差百分比分別為26.281%、16.348%、13.955%、13.318%，累計方差貢獻率為69.901%，表明這4個主成分反映了原始變量提供的69.901%的信息。各主成分的因子載荷如表5所示，初始因子載荷矩陣由SPSS計算得到，而主成分載荷矩陣則由初始因子載荷矩陣除以主成分對應特征值的平方根計算而來。從表4可以看出：總氮、生化需氧量、高錳酸鹽指數、透明度、水溫在成分1上有較大的載荷，反應水體耗氧污染水平，及水體受工業(yè)污染、農業(yè)面源的污染、有機污染和水體營養(yǎng)狀態(tài)4種不同的環(huán)境問題；氨氮、葉綠素a在成分2上有較大載荷，反應水體的營養(yǎng)狀態(tài)及水體受生活污水污染狀況；pH值、溶解氧在成分3上有較大載荷，可作為解釋水體自凈能力特征的指標[24]；總磷在成分4上有較大載荷，反應水體受農業(yè)面源污染狀況。由此表明，彰武水庫營養(yǎng)狀態(tài)變化的主要驅動因子為總氮、生化需氧量、高錳酸鹽指數、葉綠素a、pH值和總磷，水環(huán)境中有機污染源與可還原污染物并存，這與彰武水庫流域工業(yè)污染、農業(yè)面源污染、附近農地的養(yǎng)分流失、季節(jié)性的地表徑流帶來的水土流失徑流污染等有關[5]。

表3 彰武水庫主成分特征值Tab.3 Principal component eigenvalues of Zhangwu reservoir

表4 主成分矩陣Tab.4 Principal component system matrix

經調查發(fā)現：彰武水庫上游洹河幾乎沒有源頭水，主要接納林州市多家污水處理廠出水及沿線農村生活污水等，污水處理廠出水的不穩(wěn)定，且污水處理廠出水水質標準值均高于地表水Ⅴ類標準限制，從而增加了水體污染負荷；水庫上游無天然徑流的匯入，水體自凈能力較差；水庫上游匯水支流常年斷流，成為周邊居民的納污溝，垃圾污水隨雨水沖入干流匯入彰武水庫，影響其水質；彰武水庫周邊分布10個村莊，雖然建設了分散型污水處理設施，但由于管網不配套、收水量小等問題的存在，使得部分村莊的生活污水仍以面源形式直接排入水庫，對其水質造成一定影響；水庫周邊農業(yè)種植現象，造成一定程度的水土流失；雖然彰武水庫于2016年、2019年開展了庫區(qū)網箱水產養(yǎng)殖、沿岸畜禽養(yǎng)殖場的整治工作，但底泥中仍殘留大量餌料，造成總氮的偏高。綜上，調查得到的結果與主成分分析得到的污染物來源相一致。

2.3.2 2011—2020年水庫水質各水期主要影響因子分析

分別對豐水期、平水期、枯水期各水質指標進行降維處理，根據特征值大小提取主成分，各時期特征值見表5。由表5可知，不同時期提取的特征值大于1的主成分因子均為4個，三者的累積貢獻率分別為86.307%、83.641%、82.511%。

表5 彰武水庫不同時期主成分特征值及主要參數Tab.5 Principal component eigenvalues and main parameters in differentperiods of Zhangwu Reservoir

根據成分得分系數可知，在豐水期，生化需氧量、高錳酸鹽指數、氨氮、pH值、總氮在成分1上有較大載荷，表明污染物主要來源于農業(yè)面源污染和生活污水；平水期，pH值、生化需氧量、溶解氧在成分1上有較大載荷，表明污染物主要來源于生活污水，總氮、透明度、高錳酸鹽指數、總磷在成分2上有較大載荷，表明污染物主要來源于工業(yè)污染和農業(yè)面源污染；枯水期，生化需氧量、總氮、透明度在成分1上有較大載荷，表明污染物主要來源于生活污水和農業(yè)面源污染?？梢钥闯觯诓煌?，水庫受到的主要污染物是不同的，但生化需氧量、總氮是水庫水質污染的重要污染指標，生活污水、農業(yè)面源污染、工業(yè)污染是影響彰武水庫水質的主要因素，這與唐敏等的研究相一致[2]。

2.4 綜合營養(yǎng)狀態(tài)指數(TLI)分析

根據綜合營養(yǎng)狀態(tài)指數(TLI)計算結果可知，2011—2020年彰武水庫各水期水質綜合營養(yǎng)狀態(tài)指數為37.5～57.7，均在中營養(yǎng)與輕度富營養(yǎng)之間。從2017年開始，彰武水庫富營養(yǎng)化程度得到較大改善，由輕度富營養(yǎng)轉為中營養(yǎng)狀態(tài)。這可能是由于彰武水庫在2016年和2019年開展了庫區(qū)網箱水產養(yǎng)殖、沿岸畜禽養(yǎng)殖場整治工作，使得污染物攝入減少，水質得到較大改善，但網箱水產養(yǎng)殖仍有大量餌料殘留在庫區(qū)底泥中，使得總氮成為2017年之后影響水體富營養(yǎng)化的主要因子。

2.5 加拿大水質指數(CCME WQI)評價

根據加拿大水質指數(CCME WQI)計算結果可知，2011—2020年彰武水庫各水期水質指數為70.7～100，2012年5月水質極好外，其余時段均為中等～良好之間。據原始記錄分析，總氮、總磷是主要超標因子，這與主成分分析及綜合營養(yǎng)狀態(tài)指數分析結果是一致的。

2.6 主成分分析與綜合營養(yǎng)狀態(tài)指數分析、加拿大水質指數(CCME WQI)對比

將主成分載荷矩陣中的系數與標準化后的監(jiān)測數據對應相乘后求和，得到4個主成分的表達式F1、F2、F3、F4，再將這4個主成分的值乘以各自對應特征值的方差貢獻率后，再相加，得到主成分綜合得分F，進一步計算各因子得分和綜合得分。

F1=-0.341×水溫+0.157×pH值+0.413×透明度+0.248×溶解氧+0.387×高錳酸鹽指數+0.419×生化需氧量+0.028×氨氮+0.448×總氮+0.182×總磷+0.258×葉綠素a

F2=0.296×水溫-0.001×pH值-0.231×透明度-0.348×溶解氧+0.478×高錳酸鹽指數+0.038×生化需氧量-0.501×氨氮-0.114×總氮+0.133× 總磷+0.476×葉綠素a

F3=-0.057×水溫+0.675×pH值-0.346×透明度+0.500×溶解氧-0.032×高錳酸鹽指數-0.144×生化需氧量-0.255×氨氮+0.003×總氮-0.276×總磷+0.088×葉綠素a

F4=0.400×水溫-0.041×pH值-0.119×透明度-0.086×溶解氧-0.131×高錳酸鹽指數+0.473×生化需氧量+0.348×氨氮+0.202×總氮-0.638×總磷+0.316×葉綠素a

F=0.263×F1+0.163×F2+0.140×F3+0.133×F4

主成分分析(PCA)與綜合營養(yǎng)狀態(tài)指數(TLI)、加拿大水質指數法(CCME WQI)如圖2所示，3種評價方法的得分變化規(guī)律大致相同。TLI最高出現在2015年1月，TLI為57.7，輕度富營養(yǎng)，CCME WQI為75.0，為中等水質，當月總磷濃度為0.1 mg/L，達到湖庫Ⅳ類標準，而平水期總磷是影響其綜合營養(yǎng)狀態(tài)的主要因子，3種方法結果相一致。

同時，又對主成分分析得分與TLI、CCME WQI進行Spearman相關分析，主成分分析得分與TLI呈顯著正相關(r2=0.731,P<0.01,雙尾)，與CCME WQI呈顯著性負相關(r2=-0.620,P<0.01,雙尾)，TLI與CCME WQI呈顯著性負相關(r2=-0.646,P<0.01,雙尾)(表6)。結合主成分分析與TLI得分越高水質狀況越差，而CCME WQI得分越高水質狀況越好可知，主成分分析與TLI、CCME WQI評價結果相一致，均可較好地應用于彰武水庫水質狀態(tài)評價中。

表6 綜合營養(yǎng)狀態(tài)指數、加拿大水質指數與主成分分析得分相關性分Tab.6 Correlation analysis between trophic level index,Canadian water quality index and principal component analysis score

圖2 綜合營養(yǎng)狀態(tài)指數法、加拿大水質指數和主成分分析法得分情況Fig.2 Score of trophic level index method,Canadian water quality index and Principal Component analysis

3 結語

本文以彰武水庫出庫口作為研究對象，分析2011—2020年年際及每年豐、平、枯水期不同時期水質指標的變化特征，并采用主成分分析方法，提取水質污染因子并分析污染源。結果顯示：1)通過主成分分析法，彰武水庫營養(yǎng)狀態(tài)變化的主要驅動因子為總氮、生化需氧量、高錳酸鹽指數、葉綠素a、pH值和總磷，且不同水期，其主要污染物不同，但生化需氧量、總氮始終是水庫水質污染的重要污染指標，這主要與生活污水、農業(yè)面源污染、工業(yè)污染、附近農地的養(yǎng)分流失、地表徑流帶來的水土流失污染等有關;2)2011—2020年彰武水庫各水期水質綜合營養(yǎng)狀態(tài)指數為37.5～57.7，均在中營養(yǎng)與輕度富營養(yǎng)之間，而CCME WQI指數為70.7～100，2012年5月水質極好外，其余時段均為中等～良好之間。通過主成分分析法與綜合營養(yǎng)狀態(tài)指數法、加拿大水質指數法比較，發(fā)現三者有較好的一致性，均可以較好的應用于彰武水庫水質狀態(tài)評價中。

基于彰武水庫主要污染特征，筆者建議：加快污水管網建設、強化工農業(yè)污染減排，深化河道內源治理，加強對廢水直排企業(yè)的監(jiān)管和監(jiān)督監(jiān)測；加快彰武水庫周邊農村環(huán)境的綜合整治，合理規(guī)劃建設農村污水處理設施及配套管網，確保設施的良好運行效果，污水不外排，控制農業(yè)面源污染，改善水庫周邊生態(tài)環(huán)境；加大庫區(qū)及匯入支流垃圾清理、監(jiān)管力度，確保改善庫區(qū)水環(huán)境質量。