李 萍

（昆明理工大學環(huán)境科學與工程學院，云南 昆明 650000）

0 引言

洱海是云南省第二大淡水湖泊，大理市飲用水源地，蒼山洱海國家級自然保護區(qū)和風景名勝區(qū)，具有提供工農(nóng)業(yè)生產(chǎn)用水和保持生物多樣性等多種功能，是大理市乃至大理州經(jīng)濟社會可持續(xù)發(fā)展的基礎(chǔ)[1]。

近年來由于洱海周圍人口急劇增加，農(nóng)業(yè)、旅游業(yè)迅速發(fā)展，大量營養(yǎng)物質(zhì)排入洱海，致使洱海富營養(yǎng)化嚴重[2]。面對日益嚴峻的洱海保護形勢，開啟了洱海保護搶救模式，通過“七大行動”、流域轉(zhuǎn)型發(fā)展等各項措施，洱海水質(zhì)下滑趨勢得到遏制，但洱海水質(zhì)尚未得到根本性改善，水生態(tài)環(huán)境仍然脆弱，仍處在敏感的富營養(yǎng)化轉(zhuǎn)型期[1]。隨著洱海流域人口增長、空間開發(fā)強度加大、水資源利用量增加，洱海的保護治理仍然面臨著嚴峻挑戰(zhàn)。

本文選取洱海2011—2020年的監(jiān)測數(shù)據(jù)，采用主成分分析法進行水質(zhì)評價，分析研究洱海水質(zhì)時空變化特征，對洱海保護治理工作提供一定的參考。

1 材料和方法

1.1 采樣點位置

根據(jù)洱海的地理情況、入出湖特點，在洱海北部、中部、南部設(shè)置了11個測點，其中北部設(shè)置桃源（631）、湖心0（632）、雙廊（633）、喜洲（280）、湖心1（281）共5個測點，中部設(shè)置龍龕（283）、湖心（284）、塔村（285）3個測點，南部設(shè)置小關(guān)邑（286）、湖心3（287）、石房子（288）3個測點，詳見圖1。

圖1 洱海采樣布點圖

1.2 監(jiān)測時間及監(jiān)測指標

2011—2020年每月1—10日對洱海11個測點進行采樣監(jiān)測，監(jiān)測指標為pH、DO、CODMn、CODCr、BOD5、NH3-N、TP、TN、chla、SD等10項。

1.3 主成分分析法

主成分分析法本質(zhì)上是一種簡化數(shù)據(jù)和降低數(shù)據(jù)維數(shù)的方法，廣泛應(yīng)用于包括水質(zhì)評價在內(nèi)的多個領(lǐng)域。對于十幾種水質(zhì)因子變量描述的數(shù)據(jù)結(jié)果，變量多可以更加細致地描述水質(zhì)狀況，但變量多，處理難度增大，增加分析的復雜程度。主成分分析通過先找到含信息量最大的水質(zhì)變量，然后考察該水質(zhì)變量包含的信息量，當其不足以表達水質(zhì)數(shù)據(jù)包含的信息時，繼續(xù)找下一個變量，直到滿足相關(guān)表達要求為止[3-6]。

2 分析結(jié)果與討論

2.1 數(shù)據(jù)分析

（1）標準化處理及相關(guān)系數(shù)計算

通過SPSS 23.0對洱海11個測點監(jiān)測數(shù)據(jù)進行標準化處理[7]，對標準化后的數(shù)據(jù)進行主成分分析，進行變量相關(guān)性檢驗（見表1）、KMO和Bartlett檢驗（見表2）。由表2可知KMO 取樣適切性量數(shù)為0.669，＞0.5， Bartlett形度檢驗顯著性均＜0.05，表明各個變量具有相關(guān)性，適合進行主成分分析。

表1 監(jiān)測指標相關(guān)性矩陣

表2 KMO和BartIett檢驗結(jié)果

（2）主成分的確定

計算各成分特征值、方差百分比和累計方差百分比（見表3）。確定了三個主成分，累計方差百分比達到74.377%，可以對洱海水質(zhì)空間特征進行評價。

表3 主成分分析總方差分布表

分析發(fā)現(xiàn)在初始主成分載荷矩陣中，各因子差異不明顯，載荷反映信息過多。為了更好地解釋主成分，采用方差最大旋轉(zhuǎn)（ Varimax） 的方法對初始載荷矩陣進行旋轉(zhuǎn)，得到旋轉(zhuǎn)后的主成分載荷矩陣（見表4）。

從表4來看，與第一主成分密切相關(guān)的是CODCr、NH3-N、TN，與第二主成分密切相關(guān)的是pH、TP，與第三主成分密切相關(guān)的是DO、SD?？蓪⒍Ｎ廴疚镏饕秩?，一類為有機物污染，另一類為富營養(yǎng)化物質(zhì)，然后以pH、DO、SD為代表的因子，將其視為藻類生境指標。

（3）主成分得分

利用旋轉(zhuǎn)后的成分矩陣除以各主成分對應(yīng)特征值的平方根，求得相應(yīng)特征向量，將標準化后的數(shù)據(jù)和特征向量相乘[8]，得到如下的主成分表達式：

其中：X1（pH）、X2（DO）、X3（CODMn）、X4（CODCr）、X5（BOD5）、X6（NH3-N）、X7（TP）、X8（TN）、X9（chla）、X10（SD）。

根據(jù)公式（1）~（4），計算得出2011—2020年洱海水質(zhì)主成分分析得分值（見表5），綜合得分越高，說明樣品污染物含量越高。

表5 2011—2020年洱海水質(zhì)空間變化主成分分析得分值統(tǒng)計表

2.2 水質(zhì)時空特征分析

從2011—2020年洱海11個測點的水質(zhì)單因子評價結(jié)果 （見表6） 來看，水質(zhì)類別在Ⅱ~Ⅲ類波動。由于單因子評價法是通過最差水質(zhì)指標來確定水質(zhì)類別，盡管洱海11個測點的水質(zhì)存在一定的變化，但通過單因子評價法不能充分反映水質(zhì)的動態(tài)變化過程和規(guī)律；根據(jù)主成分分析的結(jié)果（表5），所有水質(zhì)指標的資料已納入主成分分析內(nèi)，可定量反映洱海水質(zhì)的空間特征。從主成分綜合得分值的數(shù)值特征（表5和圖2）來看，洱海11個測點的水質(zhì)綜合得分值排序是: 桃源＞龍龕＞雙廊＞湖心0＞喜洲＞小關(guān)邑＞石房子＞湖心3＞塔村＞湖心1＞湖心。由此可見北部桃源水質(zhì)較差，綜合得分值是0.55，是11個測點中分值最高的；而中部湖心的水質(zhì)較好，綜合得分值是-0.52，反映出洱海中部的水質(zhì)最好，南部次之，北部最差。

表6 2011—2020年洱海水質(zhì)空間綜合評價結(jié)果

圖2 洱海主成分綜合得分均值空間變化圖

洱海流域傳統(tǒng)農(nóng)業(yè)發(fā)達，種植面積較大，隨著經(jīng)濟發(fā)展及人口增加，肥料使用的加大，而灌溉及耕作還在沿用傳統(tǒng)模式，導致農(nóng)業(yè)面源污染成為洱海流域重要污染源。洱海入湖河流集中在北部和南部區(qū)域，北部彌苴河、羅時江、永安江、南部波羅江流經(jīng)大面積人口集中地區(qū)，受農(nóng)業(yè)面源污染、散養(yǎng)畜禽養(yǎng)殖、生活污水等干擾較重，同時北部彌苴河、羅時江、永安江是洱海入湖水量最大的三條河流，大量入湖河水帶入污染物，導致北部水質(zhì)最差，南部次之，中部最好[9-11]。

從圖3中可以看出2011—2014年洱海水體污染變輕，水質(zhì)變好，在2014年值最低，綜合得分-1.57；2015—2018年水質(zhì)污染逐步加重，在2018年達到峰值，綜合得分1.56，為近10年最差年份；隨后2019年、2020年水質(zhì)有持續(xù)改善趨勢。

圖3 2011—2020年洱海主成分綜合得分均值變化趨勢

綜合來看，呈現(xiàn)這一變化趨勢的主要原因是近10年來，井噴式的旅游帶動洱海流域餐飲客棧服務(wù)業(yè)的發(fā)展，也帶來人口聚集和生產(chǎn)生活方式的變化，而洱海流域截污治污設(shè)施建設(shè)投入不足，管網(wǎng)覆蓋率和污水收集處理率低。同時，受農(nóng)業(yè)面源污染、生產(chǎn)生活污水影響，主要入湖河流水質(zhì)較差，清水入湖量大幅減少，湖水不能良性循環(huán)。雖然多次采取整治措施，卻一直不能根除污染源，導致洱海陷入邊治理邊污染、水質(zhì)污染不斷加重的惡性循環(huán)，保護治理已經(jīng)迫在眉睫。通過截污治污、生態(tài)搬遷、流域轉(zhuǎn)型發(fā)展、面源污染控制等措施，洱海流域生活污水得有效處置，餐飲客棧等得到有效管控，流域大蒜全面禁種，化肥使用量減少，面源污染得到一定控制，三庫連通工程實施及蒼山十八溪無序取水整治使得清水補充進入洱海，增加清水入湖量，使洱海2019—2020年水質(zhì)有了一定的改善，但是要達到水功能目標還任重道遠，需要持續(xù)推進洱海保護治理工作[12]。

3 結(jié)論

（1）洱海水質(zhì)污染情況主要取決于CODCr、NH3-N、TN、TP、pH、DO、SD等指標，水體污染主要以有機物污染及富營養(yǎng)化為主。

（2）2011—2020年洱海北部桃源水質(zhì)較差，而中部湖心的水質(zhì)較好，總體呈現(xiàn)洱海中部的水質(zhì)最好，南部次之，北部最差。

（3）2011—2014年洱海水質(zhì)污染逐漸變輕，在2014年值最低；2015—2018年水質(zhì)污染有加重趨勢，在2018年值最高；2019年、2020年水質(zhì)持續(xù)改善。