李 雪
(遼寧省鞍山水文局,遼寧 鞍山)
我國(guó)江河眾多,共有7大流域,包括長(zhǎng)江流域、黃河流域、珠江流域、海河流域、淮河流域、松遼流域以及太湖流域。但是目前我國(guó)水資源短缺,受地形、氣候的影響,河流在地區(qū)上分布不均勻,水質(zhì)受到不同程度的污染,河流水質(zhì)的狀況對(duì)大自然環(huán)境以及生物包括人類有著直接影響,所以對(duì)河流水質(zhì)現(xiàn)狀的分析與評(píng)價(jià)是當(dāng)前全世界研究的重要課題,通過(guò)掌握水質(zhì)狀況、變化特征和發(fā)展趨勢(shì),為水資源的利用、水環(huán)境的治理與保護(hù)、污染源的控制提供依據(jù)[1]。
目前,水質(zhì)評(píng)價(jià)常用的方法有單因子評(píng)價(jià)法、人工神經(jīng)網(wǎng)絡(luò)法、模糊評(píng)價(jià)法、灰色系統(tǒng)法、主成分分析法和層次分析法等[2]。主成分分析法有能夠從眾多因子中篩選出主要成分,對(duì)數(shù)據(jù)結(jié)構(gòu)進(jìn)行簡(jiǎn)化,為后續(xù)工作篩選主要因子提供依據(jù)等優(yōu)點(diǎn),是一種將多維因子納入同一系統(tǒng)進(jìn)行定性、定量化研究的理論比較完善的多元統(tǒng)計(jì)分析方法[3-4],評(píng)價(jià)過(guò)程簡(jiǎn)單,減少主觀誤差,客觀地確定各個(gè)原始指標(biāo)的權(quán)重,對(duì)原始指標(biāo)進(jìn)行線性組合,形成幾個(gè)新的綜合指標(biāo),且這幾個(gè)綜合指標(biāo)互不相關(guān),這就達(dá)到了簡(jiǎn)化數(shù)據(jù)的目的;同時(shí),這幾個(gè)綜合指標(biāo)又能保留原始數(shù)據(jù)中的大部分信息,從而更容易抓住復(fù)雜環(huán)境問(wèn)題中的主要矛盾,提高分析效率[5-6],評(píng)價(jià)結(jié)果準(zhǔn)確,更好地綜合反映水質(zhì)的情況。
主成分分析的基本思想是設(shè)法將原來(lái)眾多具有一定相關(guān)性(如p個(gè)指標(biāo)),重新組合成一組新的互相無(wú)關(guān)的綜合指標(biāo)來(lái)代替原來(lái)的指標(biāo)。進(jìn)行主成分分析主要步驟如下[7]:
(1)運(yùn)用SPSS軟件對(duì)指標(biāo)數(shù)據(jù)標(biāo)準(zhǔn)化處理。假定有n個(gè)樣本,每個(gè)樣本共有p個(gè)變量,構(gòu)成一個(gè)n乘以p階的數(shù)據(jù)矩陣:
經(jīng)過(guò)標(biāo)準(zhǔn)化變換為
(2)計(jì)算相關(guān)系數(shù)矩陣R=[rij]。
(3)計(jì)算特征值與特征向量。解特征方程|λIp-R|=0,求出特征值,并使其按大小順序排列:λ1≥λ2≥…≥λp≥0;然后分別求出對(duì)應(yīng)于特征值λi的特征向量。
(5)計(jì)算主成分載荷lij。
(6)各主成分得分Fm=lm1x1+lm2x2+…+lmpxp。
太子河是遼寧省較大河流之一,位于遼寧省東部,是渾河支流。對(duì)太子河流域水質(zhì)進(jìn)行分析與評(píng)價(jià)可為太子河流域合理開發(fā)利用水資源、優(yōu)化配置產(chǎn)業(yè)結(jié)構(gòu)提供依據(jù),從而實(shí)現(xiàn)太子河流域水資源、人口、社會(huì)經(jīng)濟(jì)和生態(tài)環(huán)境的協(xié)調(diào)可持續(xù)發(fā)展。文章對(duì)太子河流域鞍山境內(nèi)水質(zhì)進(jìn)行評(píng)價(jià),鞍山境內(nèi)太子河流域的支流有南沙河、楊柳河以及海城河等,在南沙河上選取郎家廠監(jiān)測(cè)斷面,楊柳河上選取新堡監(jiān)測(cè)斷面,海城河上選取海城監(jiān)測(cè)斷面,通過(guò)對(duì)入太子河干流前的太子河支流的水質(zhì)進(jìn)行評(píng)價(jià),掌握太子河水質(zhì)現(xiàn)狀以及對(duì)成因進(jìn)行分析。
文章根據(jù)2021年的水質(zhì)監(jiān)測(cè)數(shù)據(jù),排除大部分月份小于檢出限的項(xiàng)目,選取pH值(X1)、溶解氧(X2)、高錳酸鹽指數(shù)(X3)、化學(xué)需氧量(X4)、五日生化需氧量(X5)、氨氮(X6)、總磷(以P計(jì))(X7)、氟化物(X8)、六價(jià)鉻(X9)、鉛(X10)、揮發(fā)酚(X11)11個(gè)項(xiàng)目作為評(píng)價(jià)指標(biāo),利用Excel 2010和SPSS17.0對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)計(jì)算和主成分分析[8]。
對(duì)南沙河上郎家廠監(jiān)測(cè)斷面、楊柳河上新堡監(jiān)測(cè)斷面、海城河上海城監(jiān)測(cè)斷面的2020年12月到2021年11月全年的監(jiān)測(cè)數(shù)據(jù)進(jìn)行分析與評(píng)價(jià),監(jiān)測(cè)數(shù)據(jù)見表1,并對(duì)監(jiān)測(cè)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理。
表1 各監(jiān)測(cè)斷面水質(zhì)監(jiān)測(cè)數(shù)據(jù) 單位:mg/L
由表2可知,郎家廠監(jiān)測(cè)數(shù)據(jù)的前4個(gè)主成分的特征根均大于1,分別為2.950、2.265、1.939和1.204;新堡監(jiān)測(cè)數(shù)據(jù)的前5個(gè)主成分的特征根均大于1,分別為3.118、2.242、1.626、1.491和1.084;海城監(jiān)測(cè)數(shù)據(jù)的前5個(gè)主成分的特征根均大于1,分別為2.953、2.432、1.670、1.460和1.171。表明這些主成分可反映11個(gè)原始指標(biāo)提供的絕大部分信息,其中,第一主成分的方差貢獻(xiàn)率大于其余主成分,說(shuō)明各評(píng)價(jià)斷面的污染物主要體現(xiàn)在第一主成分中[1]。
表2 特征值及主成分貢獻(xiàn)率
由表3可知,郎家廠監(jiān)測(cè)斷面高錳酸鹽指數(shù)、化學(xué)需氧量、五日生化需氧量、氟化物在第一主成分上有較高載荷,說(shuō)明第一主成分主要反映了這些原始指標(biāo)的信息。同理,pH值、溶解氧、總磷,鉛、揮發(fā)酚和氨氮、六價(jià)鉻分別在第二、第三和第四主成分上有較高載荷。
表3 初始因子載荷矩陣
新堡監(jiān)測(cè)斷面pH值、五日生化需氧量、總磷、六價(jià)鉻在第一主成分上有較高載荷,說(shuō)明第一主成分主要反映了這些原始指標(biāo)的信息。同理,高錳酸鹽指數(shù)、化學(xué)需氧量,氨氮、鉛,溶解氧、揮發(fā)酚和氟化物分別在第二、第三、第四和第五主成分上有較高載荷。
海城監(jiān)測(cè)斷面溶解氧、高錳酸鹽指數(shù)、化學(xué)需氧量在第一主成分上有較高載荷,說(shuō)明第一主成分主要反映了這些原始指標(biāo)的信息。同理,五日生化需氧量、總磷、六價(jià)鉻、揮發(fā)酚,pH值、鉛,氟化物和氨氮分別在第二、第三、第四和第五主成分上有較高載荷。
這些指標(biāo)主要為生活和工業(yè)污染因子,其中,溶解氧、化學(xué)需氧量、五日生化需氧量、高錳酸鹽指數(shù)和揮發(fā)酚為水體有機(jī)污染綜合指標(biāo),pH值、氨氮、鉛、總磷、氟化物和六價(jià)鉻為無(wú)機(jī)有毒有害物質(zhì),這表明水體水質(zhì)狀況主要受生活污水和工業(yè)廢水污染影響。
計(jì)算各評(píng)價(jià)斷面的主成分得分和綜合得分,對(duì)水體污染程度進(jìn)行定量化描述,得分越大,表明水質(zhì)狀況越差,結(jié)果如圖1、表4所示。郎家廠監(jiān)測(cè)斷面豐水期8—9月和平水期10—11月水質(zhì)優(yōu)于其他時(shí)期,新堡監(jiān)測(cè)斷面豐水期7—9月和平水期10—11月水質(zhì)優(yōu)于其他時(shí)期,海城監(jiān)測(cè)斷面除了6、7月份水質(zhì)差一些,其余月份水質(zhì)均較好??偟膩?lái)說(shuō),太子河鞍山境內(nèi)的三條支流南沙河、楊柳河以及海城河全年水質(zhì)總體較好。
表4 各評(píng)價(jià)斷面水質(zhì)綜合評(píng)價(jià)結(jié)果
圖1 各監(jiān)測(cè)斷面全年水質(zhì)綜合變化趨勢(shì)
郎家廠是太子河支流南沙河上入太子河干流之前的監(jiān)測(cè)斷面,評(píng)價(jià)結(jié)果表明,豐水期8—9月和平水期10—11月水質(zhì)優(yōu)于其他月份,南沙河主要接納選礦廠、化學(xué)工業(yè)公司以及水泥廠等廢污水排放,化學(xué)需氧量、氨氮、總氮、總磷和氟化物排放量比較大,但從全年評(píng)價(jià)結(jié)果來(lái)看,2月份的污染相對(duì)嚴(yán)重一些,可能是因?yàn)?月份處于枯水期,河流流量減少、流速降低,使得水體中污染物濃度稀釋擴(kuò)散能力減弱,這對(duì)水體的自凈能力不利;其余月份水質(zhì)較好,說(shuō)明河流自凈能力可以承受廢水的排放。
新堡是太子河支流楊柳河上入太子河干流之前的監(jiān)測(cè)斷面,評(píng)價(jià)結(jié)果表明,豐水期7—9月和平水期10—11月水質(zhì)優(yōu)于其他月份,楊柳河主要接納燒結(jié)廠、紡織廠以及畜牧業(yè)等廢污水排放,化學(xué)需氧量、氨氮、總氮、總磷和氟化物排放量比較大,但從全年評(píng)價(jià)結(jié)果來(lái)看,6月份的污染相對(duì)嚴(yán)重一些,6月份是豐水期,可能是因?yàn)樯a(chǎn)生活向河流排放污染物量增多超過(guò)水體自凈能力,其他月份水質(zhì)現(xiàn)狀較好。
海城是太子河支流海城河上入太子河干流之前的監(jiān)測(cè)斷面,評(píng)價(jià)結(jié)果表明,除了6、7月份水質(zhì)差一些,其余月份水質(zhì)均較好。海城河主要接納鎂礦廠、印染廠、紡織廠、食品加工廠、造紙廠及居民生活污水等廢污水排放,化學(xué)需氧量、氨氮、總氮、總磷和氟化物排放量比較大,但從全年評(píng)價(jià)結(jié)果來(lái)看,6、7月份的污染相對(duì)嚴(yán)重一些,可能是因?yàn)樯a(chǎn)生活向河流排放污染物量增多超過(guò)水體自凈能力。
文章利用主成分分析法對(duì)太子河鞍山境內(nèi)的三條支流南沙河、楊柳河以及海城河水質(zhì)現(xiàn)狀進(jìn)行了分析與評(píng)價(jià),評(píng)價(jià)結(jié)果表明,三條支流水質(zhì)較好,每個(gè)監(jiān)測(cè)斷面只有1~2個(gè)月份水質(zhì)稍差,主要是因?yàn)榕欧盼廴疚锪砍^(guò)水體自凈能力,根源在于枯水期水體自凈能力下降以及生產(chǎn)生活廢水的過(guò)量排放。評(píng)價(jià)結(jié)果可為水資源保護(hù)和發(fā)展提供科學(xué)依據(jù),政府各部門要完善相關(guān)的法律法規(guī),根據(jù)豐水期、平水期、枯水期以及季節(jié)變化調(diào)整污廢水排放量,同時(shí)加強(qiáng)對(duì)工業(yè)廢水排放的監(jiān)管力度,多措并舉,減輕河流的污染負(fù)荷,堅(jiān)持綠水青山就是金山銀山的環(huán)保理念。