吳轉(zhuǎn)璋 耿天召 伍震威 王歡

摘要 采用單因子評價法、模糊綜合評價法和主成分分析法對青弋江流域7個監(jiān)測斷面2017—2019年6項水質(zhì)指標的監(jiān)測數(shù)據(jù)進行分析，識別了流域的主要水質(zhì)指標，得出青弋江流域的綜合水質(zhì)狀況和水質(zhì)時空變化趨勢。結(jié)果表明：青弋江流域整體水質(zhì)為優(yōu)，主要水質(zhì)指標為NH 3-N和TP;下游因工業(yè)污染源較多，水質(zhì)最差，上游主要污染源為林地腐殖質(zhì)等帶來的N、P污染，隨著水體自凈中游水質(zhì)優(yōu)于上游;汛期由于雨水沖刷帶來的農(nóng)業(yè)面源污染加重，水質(zhì)劣于其他月份，隨著對流域污染源排放管控的加強，流域水質(zhì)隨時間推移呈優(yōu)化趨勢。

關(guān)鍵詞 青弋江;單因子評價;模糊綜合評價;主成分分析;時空變化趨勢

中圖分類號 X 824? 文獻標識碼 A? 文章編號 0517-6611（2021）18-0079-05

doi：10.3969/j.issn.0517-6611.2021.18.020

開放科學(xué)（資源服務(wù)）標識碼（OSID）：

Water Quality Evaluation and Spatiotemporal Variation Characteristics of Qingyi River Basin

WU Zhuan-zhang1，GENG Tian-zhao2，WU Zhen-wei1 et al （1.Anhui Province Eco-Environmental Monitoring Center，Hefei， Anhui 230071;2.Suzhou Ecological Environment Bureau，Suzhou，Anhui 234000）

Abstract The single factor evaluation method， fuzzy comprehensive evaluation method and principal component analysis method were used to analyze the monitoring data of six water quality indicators of seven sections in Qingyi River Basin from 2017 to 2019. The main water quality indexes of the basin were identified. The comprehensive water quality and the spatiotemporal trends of water quality were obtained. The results showed that the overall water quality of Qingyi River Basin was excellent， and the main water quality indicators were NH 3-N and TP; the downstream water quality was the worst due to the large number of industrial pollution sources， and the main pollution source in the upstream was N and P pollution caused by forest humus. With the self purification of water body， the water quality in the middle reaches was better than that in the upper reaches; the water quality in flood season was worse than that in other months due to the agricultural non-point source pollution caused by rainwater， and the water quality of the basin tends to be optimized over time with strengthening the management of pollution discharge.

Key words Qingyi River;Single factor evaluation;Fuzzy comprehensive evaluation;Principal component analysis;Spatiotemporal trends

基金項目 第二次全國污染源普查工業(yè)源普查報表制度及普查數(shù)據(jù)質(zhì)量控制實施技術(shù)支持項目（22110399005）。

作者簡介 吳轉(zhuǎn)璋（1981—），女，安徽太湖人，工程師，碩士，從事環(huán)境監(jiān)測和評價方面的研究。*通信作者，正高級工程師，從事環(huán)境管理研究。

收稿日期 2020-12-16

“十三五”以來，隨著“水十條”以及相關(guān)污染防治攻堅戰(zhàn)行動計劃的發(fā)布實施，我國水污染防治工作取得了明顯成效，水環(huán)境質(zhì)量得到顯著改善，但水環(huán)境質(zhì)量形勢依然嚴峻，水生態(tài)環(huán)境保護工作依然任重道遠[1]。水環(huán)境質(zhì)量評價是水生態(tài)環(huán)境保護中的一項基礎(chǔ)性工作，通過對監(jiān)測數(shù)據(jù)的合理分析，摸清水環(huán)境質(zhì)量的優(yōu)劣和變化特征，能夠為制定科學(xué)的水生態(tài)環(huán)境保護措施提供有效參考[2]。因此，采用合適的方法對水環(huán)境質(zhì)量現(xiàn)狀進行評價具有重要意義。

目前水環(huán)境質(zhì)量評價的方法有很多，比較常用的有單因子評價法、綜合污染指數(shù)法、模糊綜合評價法和主成分分析法等。其中單因子評價法計算簡單，可操作性強，可以直接根據(jù)評價因子的監(jiān)測值與評價標準之間的關(guān)系，得出各評價因子的達標情況、超標因子以及超標倍數(shù)等結(jié)果，能清晰明了地判斷出水體主要污染因子及超標斷面位置，因此單因子評價法適用于快速判斷水質(zhì)類別，但是該方法用最差的水質(zhì)指標等級作為該斷面的水質(zhì)類別，評價結(jié)果過于保守，難以對同一等級水質(zhì)類別進行更詳細的優(yōu)劣排序。綜合污染指數(shù)法是將各評價因子的監(jiān)測值與其評價標準之比作為單項污染指數(shù)，通過賦權(quán)綜合計算各單項污染指數(shù)得到綜合污染指數(shù)，該方法可以對整體水質(zhì)作出定量描述，但不同的賦權(quán)方法計算出的評價結(jié)果差異較大。模糊綜合評價法通過建立隸屬函數(shù)和權(quán)重集得到水質(zhì)的綜合評價結(jié)果，評價結(jié)果更加客觀全面，適用于分析水質(zhì)是否達到水質(zhì)目標要求，但目前隸屬度由人為主觀確定得出，評價結(jié)果會出現(xiàn)較大差異。主成分分析法能夠考慮到各指標間的信息關(guān)聯(lián)，在最大限度地保留原始信息的同時，對高維變量進行降維篩選出主要指標，且能根據(jù)各斷面的綜合得分值對流域水質(zhì)進行時空變化特征分析，但主成分分析法難以對水質(zhì)等級進行判定。以上評價方法側(cè)重點不同，各有優(yōu)劣，相互結(jié)合可以優(yōu)勢互補[3-9]。

為此，筆者以長江的重要一級支流青弋江作為研究對象，根據(jù)2017—2019年青弋江流域的水質(zhì)監(jiān)測資料，結(jié)合單因子評價法、模糊綜合評價法和主成分分析法對青弋江流域水質(zhì)進行評價和分析，全面客觀地反映流域水質(zhì)狀況與變化特征，從而為流域水環(huán)境管理提供科學(xué)依據(jù)。

1 材料與方法

1.1 數(shù)據(jù)來源

青弋江是長江的重要一級支流，干流全長309 km，流域面積約7 195 km2，起源于黃山市黟縣，經(jīng)石臺縣、黃山區(qū)，于周家坦注入陳村水庫（太平湖），出陳村水庫流經(jīng)涇縣、宣城、南陵、蕪湖等地，于蕪湖市區(qū)入長江[10-11]。

該研究數(shù)據(jù)來源于青弋江流域7個常規(guī)監(jiān)測斷面2017—2019年數(shù)據(jù)，斷面設(shè)置自上游至下游分別為溪口村、城關(guān)上游、琴溪橋、涇南交界、百園新村、海南渡、寶塔根。根據(jù)青弋江地理位置特征可知，溪口村和城關(guān)上游為上游斷面，琴溪橋、涇南交界和百園新村為中游斷面，海南渡和寶塔根為下游斷面。其區(qū)位分布見圖1。

1.2 數(shù)據(jù)預(yù)處理

根據(jù)對原始數(shù)據(jù)的初步分析，剔除未檢出數(shù)據(jù)以及水溫、pH等空值較多的項目，選取溶解氧（DO）、高錳酸鹽指數(shù)（COD Mn）、五日生化需氧量（BOD 5）、氨氮（NH 3-N）、化學(xué)需氧量（COD Cr）和總磷（TP）6項指標進行評價分析。

1.3 評價方法

1.3.1 單因子評價法。

單因子評價法是用水體各監(jiān)測因子的監(jiān)測結(jié)果對比該因子的水質(zhì)分類標準，確定該因子水質(zhì)類別;在所有因子的水質(zhì)類別中選取水質(zhì)類別最差的類別作為

該水體的水質(zhì)類別。水質(zhì)分類標準值根據(jù)《地表水環(huán)境質(zhì)量標準》（GB 3838—2002）確定。

1.3.2 模糊綜合評價法。

模糊評價法通過各因子的濃度值，除以對應(yīng)的質(zhì)量標準中各類別的平均值，得出權(quán)重值。再依據(jù)各監(jiān)測因子在地表水環(huán)境質(zhì)量標準各類別中的歸屬，得出隸屬度。最后通過矩陣復(fù)合運算，確定總體水環(huán)境的類別[12]。

1.3.3 主成分分析法。

主成分分析法是設(shè)法將原來變量重新組合成一組新的相互無關(guān)的綜合變量，同時根據(jù)實際需要從中取出幾個較少的但盡可能多地反映原來變量信息的綜合變量的統(tǒng)計方法。主要計算步驟如下[13-14]：①原始數(shù)據(jù)標準化;②計算各指標之間的相關(guān)系數(shù)矩陣;③計算相關(guān)系數(shù)矩陣的特征值與特征向量;④計算主成分貢獻率及累積貢獻率，確定主成分;⑤將標準化后的指標變量轉(zhuǎn)換為主成分;⑥對主成分進行加權(quán)求和，即得最終評價值并用其評判水質(zhì)，權(quán)數(shù)為每個主成分的方差貢獻率。

2 結(jié)果與分析

2.1 單因子評價結(jié)果

首先采用《地表水環(huán)境質(zhì)量評價辦法（試行）》（環(huán)辦〔2011〕22號文）規(guī)定的劃分水質(zhì)類別的單因子評價法來評價溪口村、城關(guān)上游、琴溪橋、涇南交界、百園新村、海南渡和寶塔根7個監(jiān)測斷面水質(zhì)狀況。評價結(jié)果顯示，2017—2019年青弋江流域各監(jiān)測斷面水質(zhì)類別均穩(wěn)定在Ⅱ類（表1），水質(zhì)狀況為優(yōu)，確定青弋江流域水質(zhì)類別主要指標是TP。

2.2 模糊綜合評價結(jié)果

青弋江流域水質(zhì)模糊綜合評價顯示（表1），2017—2019年青弋江流域各監(jiān)測斷面水質(zhì)類別均穩(wěn)定在Ⅰ類，水質(zhì)狀況為優(yōu)。通過對權(quán)重指標分析可知，2017—2019年7個監(jiān)測斷面的DO權(quán)重值均最大，分別為025、0.27、0.33、0.26、0.26、0.27、0.36、0.33、0.36、0.26、0.28、0.29、0.32、031、0.29、0.24、0.27、0.31、0.24、0.24、027，說明DO是影響青弋江流域水質(zhì)的主要指標。

2.3 主成分分析結(jié)果

分別對各監(jiān)測斷面不同監(jiān)測指標年均值和月均值進行2個獨立的主成分分析。采用各監(jiān)測斷面不同監(jiān)測指標年均值進行流域主要污染指標的篩選以及水環(huán)境質(zhì)量的時空變化趨勢分析，采用各監(jiān)測斷面不同監(jiān)測指標月均值進行各監(jiān)測斷面水環(huán)境質(zhì)量的時間變化趨勢分析。

運用 KMO 和 Bartlett球體檢驗進行統(tǒng)計，KMO檢驗是對原始變量之間的簡相關(guān)系數(shù)和偏相關(guān)系數(shù)的相對大小進行檢驗，Bartlett球性檢驗用于檢驗各變量是否各自獨立[15-16]。該研究計算得到年均值和月均值數(shù)據(jù)組KMO檢驗統(tǒng)計值分別為0.52和0.56，均大于0.50，Bartlett球形檢驗顯著性概率P值均小于 0.05，表明以上數(shù)據(jù)組指標間相互不獨立，相關(guān)性較強，適宜進行主成分分析。

分別對各水質(zhì)指標年均值和月均值數(shù)據(jù)組進行降維處理，結(jié)果見表2和圖2。從年均值數(shù)據(jù)組降維結(jié)果來看，方差貢獻率大于10%的主成分有4個，對應(yīng)的數(shù)值分別為3865%、27.12%、13.13%和10.55%，前4個主成分的累積方差貢獻率為89.45%，高于80%。主成分1與TP和NH 3-N呈顯著正相關(guān)關(guān)系，主成分2與DO呈顯著正相關(guān)關(guān)系，主成分3與COD Mn和BOD 5呈顯著正相關(guān)關(guān)系，主成分4與DO和COD Cr呈顯著正相關(guān)關(guān)系。以上主成分涉及所有水質(zhì)指標。

從月均值數(shù)據(jù)組降維結(jié)果來看，方差貢獻率大于10%的主成分也有4個，對應(yīng)的數(shù)值分別為31.30%、20.33%、1723%和13.58%，前4個主成分的累積方差貢獻率為8244%，高于80%。主成分1與TP、NH 3-N和COD Mn呈顯著正相關(guān)關(guān)系，主成分2與DO呈顯著正相關(guān)關(guān)系，主成分3與BOD 5呈顯著正相關(guān)關(guān)系，主成分4與COD Cr呈顯著正相關(guān)關(guān)系。以上主成分涉及所有水質(zhì)指標。

對比分析青弋江流域各監(jiān)測斷面年均值數(shù)據(jù)組和月均值數(shù)據(jù)組降維結(jié)果，可知2組數(shù)據(jù)的主成分分析結(jié)果基本一致，篩選出青弋江流域主要水質(zhì)指標為TP和NH 3-N。

青弋江流域年均值數(shù)據(jù)組水質(zhì)主成分綜合評價結(jié)果見表3，月均值數(shù)據(jù)組水質(zhì)主成分綜合評價結(jié)果見表4。

2.4 水質(zhì)時空變化特征分析

對青弋江流域年均值和月均值數(shù)據(jù)組水質(zhì)主成分綜合評價結(jié)果進行分析，得出青弋江流域水質(zhì)時空變化趨勢及各斷面水質(zhì)時間變化趨勢（圖3）。

2017—2019年青弋江流域年均值數(shù)據(jù)組水質(zhì)主成分綜合評價結(jié)果表明（表3）：青弋江流域下游水質(zhì)最差，中游優(yōu)于上游;2019年各斷面水質(zhì)均優(yōu)于2017年水質(zhì)，青弋江流域整體水質(zhì)隨時間推移呈優(yōu)化趨勢;2018年整體水質(zhì)變化趨勢不明顯，溪口村和琴溪橋2018年水質(zhì)在2017—2019年中最優(yōu)，百園新村和寶塔根水質(zhì)在3年中最差。

2017—2019年青弋江流域月均值數(shù)據(jù)組水質(zhì)主成分綜合評價結(jié)果表明（表4）：溪口村斷面2017年5、6月水質(zhì)最差，10、12月份水質(zhì)最優(yōu);2018年4、9月水質(zhì)最差，10—12月水質(zhì)最優(yōu);2019年2、6、9月水質(zhì)最差，5、10—12月水質(zhì)最優(yōu)。城關(guān)上游斷面2017年7、8月水質(zhì)最差，1月次之，其他月份曲線基本平穩(wěn);2018年9月水質(zhì)最差，6月次之，2月水質(zhì)最優(yōu);2019年5—8、11月水質(zhì)最差，1月水質(zhì)最優(yōu)。琴溪橋斷面2017年9月水質(zhì)最差，11、12月水質(zhì)最優(yōu);2018年水質(zhì)曲線基本平穩(wěn);2019年6月水質(zhì)最差，12月水質(zhì)最優(yōu)。涇南交界斷面2017年9月水質(zhì)最差，其他月份曲線基本平穩(wěn);2018年10月水質(zhì)最差，7、12月次之，其他月份曲線基本平穩(wěn);2019年8、9月水質(zhì)最差，1月最優(yōu)。百園新村斷面連續(xù)3年水質(zhì)變化曲線基本一致，7月份最差，其他月份曲線基本平穩(wěn)。海南渡斷面2017年7月水質(zhì)最差，12月次之，其他月份基本平穩(wěn);2018年9月水質(zhì)最差，12月次之，7月最優(yōu);2019年9月水質(zhì)最差，6月次之，1月最優(yōu)。寶塔根斷面連續(xù)3年水質(zhì)變化基本一致，6—9月最差，12月次之，其他月份基本平穩(wěn)。綜合分析各斷面水質(zhì)變化，得出青弋江流域汛期水質(zhì)明顯劣于其他月份。

3 討論

3.1 評價結(jié)果對比分析

該研究的3種水質(zhì)評價方法側(cè)重點不同，評價結(jié)果也略有差異。單因子評價法用最差的水質(zhì)等級代表評價結(jié)果，根據(jù)各采樣點位的TP濃度評價青弋江水質(zhì)為Ⅱ類，評價結(jié)果不夠全面，但計算結(jié)果安全性高。模糊綜合評價結(jié)果與單因子評價結(jié)果相差1個等級，評價結(jié)果顯示青弋江水質(zhì)為Ⅰ類，模糊綜合評價法充分考慮了各評價因子的綜合影響，評價結(jié)果更加客觀全面，但模糊綜合評價對指標權(quán)重矢量的確定主觀性較強，會掩蓋一些重要水質(zhì)指標的影響。主成分分析法難以對水質(zhì)等級進行判定。

影響單因子評價結(jié)果的主要水質(zhì)指標是TP，模糊綜合評價中DO權(quán)重值最大，DO對模糊綜合評價結(jié)果的影響最大，主成分分析法結(jié)果顯示第1個主成分對應(yīng)水質(zhì)指標為TP和NH 3-N，第2個主成分對應(yīng)水質(zhì)指標為DO，主成分分析法與單因子評價法和模糊綜合評價法的評價結(jié)果基本一致，且包含單因子和模糊綜合評價的主要水質(zhì)指標，因此主成分分析法更為客觀、全面。

單因子和模糊綜合評價法只能對水體進行定性評價，主成分分析法能對水體質(zhì)量進行定量描述，且能根據(jù)描述結(jié)果進行時空變化特征分析。

3.2 評價結(jié)果與實際調(diào)查情況對比分析

青弋江流域上游主要污染源是林地和農(nóng)村生活源，中游主要為農(nóng)村生活源和農(nóng)業(yè)面源，下游進入蕪湖市段主要為工業(yè)污染源，因此整個流域以富營養(yǎng)污染為主，流域下游由于工業(yè)污染源較多，水質(zhì)最差，上游主要污染源為林地腐殖質(zhì)等帶來的N、P污染，隨著水體自凈，中游水質(zhì)優(yōu)于上游，汛期由于雨水沖刷帶來的農(nóng)業(yè)面源污染加重水質(zhì)劣于其他月份;隨著對流域污染源排放管控的加強，流域水質(zhì)隨時間推移呈優(yōu)化趨勢。以上調(diào)查結(jié)果均與評價結(jié)果一致。

4 結(jié)論

（1）單因子、模糊綜合和主成分分析相結(jié)合的評價方法可用于青弋江流域水質(zhì)評價和時空變化特征的綜合分析。結(jié)果顯示，青弋江流域整體水質(zhì)狀況穩(wěn)定在Ⅱ類以上，主要水質(zhì)指標為NH 3-N和TP;青弋江流域下游水質(zhì)最差，中游水質(zhì)優(yōu)于上游;2019年各斷面水質(zhì)均優(yōu)于2017年，青弋江流域整體水質(zhì)隨時間的推移呈優(yōu)化趨勢，青弋江流域汛期水質(zhì)明顯劣于其他月份。

（2）青弋江流域整體水質(zhì)為優(yōu)，但中上游農(nóng)業(yè)面源及蕪湖市段工業(yè)點源污染仍是該流域兩大污染風(fēng)險。因此，該流域仍需加強農(nóng)業(yè)面源污染控制和下游蕪湖市段工業(yè)污染末端治理，注重污染預(yù)防，從源頭和全過程減少污染物的產(chǎn)生。

參考文獻

[1] 馬樂寬，謝陽村，文宇立，等.重點流域水生態(tài)環(huán)境保護“十四五”規(guī)劃編制思路與重點[J].中國環(huán)境管理，2020，12（4）：40-44.

[2] 程衛(wèi)國，李亞斌，蘇燕，等.不同賦權(quán)方法的綜合水質(zhì)標識指數(shù)法對比分析[J].灌溉排水學(xué)報，2019，38（11）：93-99.

[3] 初征.水環(huán)境質(zhì)量評價中的幾種方法[J].有色金屬，2010，62（3）：160-162.

[4] 陸衛(wèi)軍，張濤.幾種河流水質(zhì)評價方法的比較分析[J].環(huán)境科學(xué)與管理，2009，34（6）：174-176.

[5] 孫婷，莫躍爽，李芳，等.花溪區(qū)思雅河水質(zhì)分析與評價[J].環(huán)境監(jiān)測管理與技術(shù)，2020，32（4）：43-47.

[6] 郭晶，王丑明，黃代中，等.洞庭湖水污染特征及水質(zhì)評價[J].環(huán)境化學(xué)，2019，38（1）：152-160.

[7] 王竹，朱士江，劉揚，等.不同水質(zhì)評價方法在灤河下游段的比較應(yīng)用[J].節(jié)水灌溉，2019（10）：68-72，77.

[8] 彭福利，張永勇，李茜，等.官廳水庫入庫斷面水質(zhì)多指標評價與演變特征分析[J].中國環(huán)境監(jiān)測，2020，36（1）：65-74.

[9] 吳岳玲.水質(zhì)綜合評價及預(yù)測研究進展[J].安徽農(nóng)業(yè)科學(xué)，2020，48（2）：23-26.

[10]

李峻，孫世群.基于BP網(wǎng)絡(luò)模型的青弋江水質(zhì)預(yù)測研究[J].安徽工程科技學(xué)院學(xué)報（自然科學(xué)版），2008，23（2）：23-26.

[11] 胡晨琦，胡春生，劉永婷，等.青弋江涇縣段最高階地的ESR年代及其構(gòu)造氣候意義[J].山地學(xué)報，2017，35（4）：469-476.

[12] 楊海江，王建宇，石磊，等.模糊數(shù)學(xué)評價法在蘆河水質(zhì)評價中的應(yīng)用[J].現(xiàn)代農(nóng)業(yè)科技，2016（10）：178-180.

[13] 杜軍凱，傅堯，李曉星.模糊-主成分分析綜合評價法在地下水水質(zhì)評價中的應(yīng)用[J].中國環(huán)境監(jiān)測，2015，31（4）：75-81.

[14] 王紅，姚君蘭，李艷薔，等.基于主成分分析法的梁子湖水質(zhì)評價分析[J].湖北大學(xué)學(xué)報（自然科學(xué)版），2017，39（6）：601-608.

[15] 任芳玲，張亞楠.基于主成分分析法的用電量預(yù)測模型[J].徐州工程學(xué)院學(xué)報（自然科學(xué)版），2016，31（2）：50-53.

[16] 李成，李海波，高丹丹，等.主成分分析在城市大氣環(huán)境質(zhì)量評價中的應(yīng)用[J].湖北大學(xué)學(xué)報（自然科學(xué)版），2016，38（6）：567-571.