曹 丹

(遼寧省葫蘆島水文局，遼寧 葫蘆島 125000)

渾河大伙房以上水體為大伙房控制流域內(nèi)的全部河流和水庫，主要包括大伙房水庫及渾河上游、社河和蘇子河3條主要入庫支流[1]。大伙房水庫為遼寧中東部等城市的飲用水源地，控制流域內(nèi)水體具有連續(xù)的空間性，隨著社會經(jīng)濟快速發(fā)展，水庫控制流域內(nèi)的水體污染程度發(fā)生變化，對區(qū)域的飲水安全產(chǎn)生不同程度的影響[2]。為加大大伙房水庫飲用水源地保護力度，亟需結(jié)合大伙房水庫控制流域水體水質(zhì)監(jiān)測數(shù)據(jù)，對水庫以上水體水質(zhì)變化特征進行分析，并對水體主要污染來源進行識別，從而為大伙房水庫水體保護規(guī)劃和具體措施的制定提供支撐依據(jù)[3]。近些年來，對于大伙房水庫以上控制水體水污染特征取得一定研究成果[4-14]，但這些研究存在樣點少、研究時間早、采樣時段窄等問題，污染源空間和時間分布特征不清晰。為提高水源地飲水安全和利用價值，本文結(jié)合大伙房水庫以上控制水體的監(jiān)測數(shù)據(jù)，并結(jié)合主成分分析方法[15]和水質(zhì)指標之間、水質(zhì)指標與環(huán)境因素之間相關性分析，建立水體污染源識別技術，對大伙房水庫以上控制水體的污染源及時空分布特征進行識別。成果為《大伙房飲用水水源保護條例》的落實和最嚴格水資源管理制度實施提供技術支撐。

1 污染影響因素識別方法

1.1 識別方法

針對水體污染識別的方法研究，對不同水質(zhì)監(jiān)測指標與影響因子之間采用主成分分析方法進行關聯(lián)分析，結(jié)合區(qū)域污染特征進行動態(tài)識別。通過對線性變化的多個變量進行主要因子分析是主成分分析的常用的多元統(tǒng)計方法。其主要優(yōu)點在于對多個變量進行主成分分析，采用“降維、去相關”對原始指標主成分綜合分析。將一定關聯(lián)的因子進行重新組合分析，確定其綜合指標是主成分分析的主要方式，該方法將原始隨機變量進行正交變換成隨機新變量。采用協(xié)方差將原隨機變量進行對角轉(zhuǎn)換，新的正交坐標系對原坐標系進行幾何表現(xiàn)，使得具有p個正交方向的樣點分布，采用降維方式對多變量進行處理，新變量中具有較大方程的可以對多個原始變量進行主要信息的綜合反映，使其包含原始變量信息，主要變量構建方程為：

Zp=up1X1+up2X2+…uppXp

(1)

式中，Zp—p個主成分；Xp—p個變量；u—權重。

1.2 主成分分析的基本步驟

(1)對原有變量進行坐標變化，可得：

Zp=u1px1+u2px2+…uppxp

(2)

其中：

(3)

(2)提出主成分，z1稱為第一主成分，其滿足條件：

u′1u1=1

(4)

var(zi)=maxvar(u′x)

(5)

z2稱為第二主成分，其滿足條件：

cov(z1，z2)=0

(6)

u′2u2=1

(7)

var(z2)=maxvar(Ux)

(8)

其余主成分所滿足的條件以此類推。

2 大伙房水庫以上控制水體污染特征

采用GB 3838—2002《地表水環(huán)境質(zhì)量標準》，對大伙房水庫以上水體水質(zhì)數(shù)據(jù)進行對比分析，對其主要污染因子進行時空分布特征的識別，對控制流域內(nèi)主要入庫河流富營養(yǎng)化等級評價。

2.1 全部樣點水質(zhì)參數(shù)統(tǒng)計分析

選取的18個水質(zhì)采樣點分別位于渾河、社河和蘇子河作為大伙房水庫以上控制水體內(nèi)，各采樣點分布見表1。按照《水和廢水監(jiān)測分析方法》，進行29項水質(zhì)指標統(tǒng)計分析，結(jié)果見表2。

表1 水質(zhì)采樣點位置

由表2可知，除總氮外，其他指標均低于Ⅱ類地表水環(huán)境標準，表明大伙房水庫控制流域內(nèi)水體整體具有良好的水質(zhì)狀況。超標嚴重的2個指標為總氮和糞大腸菌群，2個指標最低值均大于對應的Ⅱ類水質(zhì)環(huán)境標準值，全年總氮在整個控制流域水體內(nèi)污染指數(shù)較高。渾河和蘇子河2個支流是超標樣本糞大腸菌群的主要來源，全年各月份均超過Ⅱ類水質(zhì)標準。大伙房水庫、紅升水庫、腰堡水庫在7月份有10個樣本pH值超標，其他各指標最高值均不超標。變異系數(shù)高于100%的水質(zhì)指標分別為氨氮、汞、揮發(fā)酚、糞大腸菌群、硫酸鹽、鐵、錳、葉綠素、懸浮物、渾濁度、色度，表明各指標具有較為明顯的異質(zhì)性時空分布特征。

表2 所有樣品29項水質(zhì)指標統(tǒng)計性描述

2.2 總氮變化特征

考慮到大伙房水庫控制流域內(nèi)總氮和糞大腸菌群超標嚴重，結(jié)合各采樣點總氮監(jiān)測數(shù)據(jù)，對大伙房水庫控制流域內(nèi)河流2014—2018年總氮時空變化特征進行分析，如圖1所示。

由圖1可知，各采樣點總氮濃度均值為2.76mg/L，濃度區(qū)間主要在0.63～9.25之間變化。有98%的樣本超過了地表水Ⅲ類水質(zhì)標準(1.0mg/L)，各樣本中濃度超過1.5mg/L(Ⅳ類水質(zhì)標準)的比例達到92%，超過2.0mg/L(Ⅴ類標準)的樣本比例達到77%，總氮污染程度較高。從總氮污染空間分布特征可看出，渾河上游總氮濃度均值為3.35mg/L，超標樣本個數(shù)達到71個，總氮污染程度最高，其次是蘇子河，其總氮濃度均值為2.95mg/L，超標樣本的個數(shù)為114個，社河總氮污染程度在3條主要河流中最低，其總氮濃度均值為2.45mg/L，有132個樣本超過標準值。

圖1 全部樣點總氮濃度分布圖

2.3 糞大腸菌群變化特征

大伙房水庫控制流域全部樣點糞大腸菌群濃度如圖2所示，共454個樣本，糞大腸菌群濃度范圍為0～198000個/L，其中33%滿足地表水Ⅰ類水質(zhì)標準(200個g/L)，59%滿足Ⅱ類水質(zhì)標準(2000個g/L)，76%滿足Ⅲ類標準(10000個g/L)，81%滿足Ⅳ類標準(20000個g/L)，89%滿足Ⅴ類標準(40000個g/L)。整體渾河上游糞大腸菌群濃度均值為121226個/L，84個樣本超過標準值，為濃度最高的河流，其次為蘇子河，其糞大腸菌群濃度均值為17371個/L，超標樣本的個數(shù)為152個，社河及庫區(qū)范圍內(nèi)糞大腸菌群濃度最低，其糞大腸菌群濃度均值分別為6165個/L和149個/L，超標樣本的個數(shù)90個和130個。

圖2 全部樣點糞大腸菌群濃度分布圖

2.4 富營養(yǎng)化指數(shù)變化特征

結(jié)合綜合加權指數(shù)法對大伙房水庫控制流域內(nèi)2014—2017年富營養(yǎng)指數(shù)變化特征進行分析，如圖3所示。

從圖3中可以看出，從2014年起，大伙房水庫控制流域內(nèi)水體富營養(yǎng)化指數(shù)逐年遞減，29%的樣本在2014年富營養(yǎng)化程度處于輕度-中度之間，7%的樣本在2017年處于輕度等級。除2017年外，每年7、8、9月富營養(yǎng)化指數(shù)高，5、10、11月低，夏季富營養(yǎng)化指數(shù)總體較高，秋冬季營養(yǎng)化指數(shù)較低，富營養(yǎng)化指數(shù)在2017年總體處于中等程度，月際沒有顯著差異。其中2017年5月富營養(yǎng)化指數(shù)整體偏高，接近富營養(yǎng)化警戒線，主要是總氮濃度超標導致，濃度高樣點渾73是大伙房水庫內(nèi)渾河入庫后的樣點，北雜木位于渾河入庫口。渾73與北雜木2點的富營養(yǎng)化指數(shù)高于其他樣點，處于輕度富營養(yǎng)化水平，并且在大部分時間下，北雜木富營養(yǎng)程度低于渾73水質(zhì)監(jiān)測采樣點。位于蘇子河和社河入庫口的古樓、臺溝2個采樣點處于中等富營養(yǎng)化程度，因此大伙房水庫庫區(qū)主要受渾河上游來水導致的富營養(yǎng)化可能性較高。樣點渾37左和渾37右位于水庫中部的同一監(jiān)測斷面，富營養(yǎng)化指數(shù)大小及變化趨勢一致。2014—2015年渾37左富營養(yǎng)化指數(shù)高于渾37右，是由于靠近渾37左的岸邊農(nóng)業(yè)種植更加密集。撫順取水口和壩前2點位于大伙房水庫出庫口附近，除2014年7月和2017年5月外，其他時間均處于中營養(yǎng)水平。

圖3 庫區(qū)及入庫口富營養(yǎng)化指數(shù)分布圖

3 大伙房水庫以上控制水體面源污染識別

3.1 污染識別結(jié)果

分別選取渾河上游、蘇子河、社河的5個采樣點，進行總氮、總磷、糞大腸菌群等18個指標污染源識別，將354個樣本的18個指標標準化數(shù)據(jù)進行主成分分析，分析矩陣見表3。

表3 主成分分析矩陣

第一個主成分解釋了總體方差的20.8%，包括總磷、氟化物、鐵、錳等指標，推測為地下水污染與底泥釋放源，來自于地下水補給及河流底泥沉積的污染物釋放。總磷和鐵的Spearman相關系數(shù)為0.513，具有顯著相關性，撫順地區(qū)地下水質(zhì)中根據(jù)《撫順市地質(zhì)環(huán)境監(jiān)測成果報告》超標的指標主要有氟化物、鐵、錳，地表水和地下水的互相補給過程中會使得超標污染物相互遷移。同時通過調(diào)查顆粒形態(tài)磷沉積在底泥之中主要來源于農(nóng)業(yè)面源，底泥在厭氧條件下其鐵、磷等污染物能夠與氫氧化三鐵進行吸附和釋放，因此地下水污染及底泥釋放是其第一個主成分的主要污染源。

氣溫、水溫和溶解氧在第二主成分中總體方差為15.5%，其中負貢獻率指標為溶解氧。溶解氧和水溫按照Spearman相關性分析具有顯著相關，水中溶解的氧因為升高的水溫使得其濃度較低，水溫和氣溫屬于氣候自然變化，而不是污染來源，因此氣候因素是第二主成分。

總氮和硝酸鹽氮在第三主成分中總體方差為9.9%，且具有顯著相關性。大氣沉降、化肥流失、動物排泄、生活污水等是水體氮素的主要來源，化肥流失和動物糞便是硝酸鹽氮的主要來源，因此第三主成分中農(nóng)田流失是總氮和硝酸鹽氮的主要來源，控制流域內(nèi)尿素用量較高，其耕地主要分布在流域兩岸丘陵地區(qū)，夏季化肥受降水和水土流失影響流失風險較高，因此農(nóng)業(yè)種植源是第三主成分污染主要來源。

五日生化需氧量、氨氮、總磷和糞大腸菌群等指標在第四主成分中方差總體較高，畜禽養(yǎng)殖和農(nóng)村生活是其污染主要來源，糞大腸菌群濃度較高的采樣點主要分布在大伙房水庫控制流域內(nèi)農(nóng)村居民人口較多的區(qū)域。糞大腸菌群和人口數(shù)按照Spearman相關性分析具有較高的關聯(lián)度，區(qū)域內(nèi)畜禽養(yǎng)殖廠主要分布在大伙房水庫控制流域河流沿岸，排放的污染物易流失進入水體。因此定義PC4為農(nóng)村生活及畜禽養(yǎng)殖源。

第五個主成分解釋了總體方差的7.4%，包括高錳酸鹽指數(shù)和化學需氧量，代表有機物污染。高錳酸鹽指數(shù)與水體中各指標均具有較高的關聯(lián)度，其主要污染為工業(yè)、生活廢水。

砷、硫酸鹽和氯化物在第六主成分中總體解釋方差為6.3%，其主要污染指標主要受地質(zhì)影響。

pH值為第七主成分，其主要和鉀鈉鈣鎂等離子具有較高的關聯(lián)度，其主要受歷史積累所影響。

3.2 主要識別污染源時空分布

對識別出2大超標水質(zhì)的污染源時空變化特征進行分析，確定污染源的高風險時間與樣點位置。

(1)農(nóng)業(yè)種植源

對大伙房水庫控制流域內(nèi)的農(nóng)業(yè)種植源得分進行時空分布統(tǒng)計，如圖4所示。

從圖4中可看出從2014—2018年大伙房水庫控制流域內(nèi)農(nóng)業(yè)種植源呈現(xiàn)先遞增后遞減的變化特征，各監(jiān)測點從2017年開始春季主成分值高，到夏季有所減小，到秋季再逐步增加的變化分布，這主要受地區(qū)農(nóng)業(yè)施肥方式有關，春季施肥量較高，污染流失率較大，夏季由于作物較為茂盛減低了降雨影響的化肥流失濾。農(nóng)作物在秋季收割后由于裸露地表增加了其流失率。渾河上游北雜木、蘇子河紅升水庫、永陵和古樓等樣點由于有面積較大的耕地，其農(nóng)業(yè)種植源影響下的污染物流失率均較高。

圖4 農(nóng)業(yè)種植源得分分布圖

(2)農(nóng)村生活及畜禽養(yǎng)殖源

對大伙房水庫控制流域內(nèi)的農(nóng)村生活及畜禽養(yǎng)殖源得分進行時空分布統(tǒng)計，如圖5所示。

圖5 農(nóng)村生活及畜禽養(yǎng)殖源得分分布圖

北雜木、木奇和永陵采樣點的農(nóng)村生活及畜禽養(yǎng)殖源得分均高于其他水質(zhì)采樣點，渾河上游和蘇子河是采樣點的主要分布河流，農(nóng)村生活及畜禽養(yǎng)殖是渾河上游和蘇子河較為嚴重的污染來源。通過對渾河、蘇子河、社河入庫口家禽調(diào)查，其數(shù)量超過7萬只，2016—2017各采樣點農(nóng)村生活及畜禽養(yǎng)殖源得分具有明顯遞減變化，這主要是因為嚴格劃定了畜禽禁養(yǎng)區(qū)。

4 結(jié)論及建議

(1)大伙房水庫控制流域2個污染因子具有相似的空間變化特征，即渾河上游>蘇子河>社河和庫區(qū)?？偟獫舛瓤傮w呈現(xiàn)逐年降低趨勢，年內(nèi)變化為春季高、夏季下降、秋冬季升高。糞大腸菌群濃度年際間無差異，同年內(nèi)先增加再降低，峰值常出現(xiàn)在7—9月。

(2)庫區(qū)水體整體處于中營養(yǎng)等級，營養(yǎng)水平為磷限制型，且富營養(yǎng)化指數(shù)呈逐年下降趨勢。應在庫區(qū)上游設置污染修復區(qū)，以降低上游污染對入庫水質(zhì)影響。

(3)本文未對大伙房水庫控制流域內(nèi)面源污染調(diào)控措施進行深入分析，存在不足，在后續(xù)研究中還應對流域內(nèi)面源污染調(diào)控措施進行深入分析。