吳會(huì)民 張韋 李文雯 蔡超 王永辰 李楠 繳建華

摘要：運(yùn)用主成分分析法對(duì)三口池塘水質(zhì)和沉積物共20個(gè)指標(biāo)的檢測(cè)結(jié)果進(jìn)行了綜合評(píng)價(jià)。結(jié)果表明，在養(yǎng)殖前期水質(zhì)最差，養(yǎng)殖中期水質(zhì)最好，養(yǎng)殖后期逐漸下降。

關(guān)鍵詞：濱海型鹽堿地池塘;理化因子;沉積物;主成分分析

中圖分類(lèi)號(hào)：X824? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：0439-8114（2019）17-0043-03

DOI：10.14088/j.cnki.issn0439-8114.2019.17.011? ? ? ? ? ?開(kāi)放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Abstract： In order to evaluate the quality of pond of coastal saline-alkaline land by principal component analysis（PCA）， 20 indicators of water and sediment were detected. The results showed that the cultural water quality was the worst at the early stage， but it turned better at the middle， and it got worse at the last stages of the breeding.

Key words： pond of coastal saline-alkaline land; physical-chemical factors; sediment; principal component analysis （PCA）

目前水環(huán)境質(zhì)量評(píng)價(jià)的方法有很多種，主要有兩大類(lèi)，即單因子評(píng)價(jià)法和綜合評(píng)價(jià)法[1]。單因子評(píng)價(jià)法計(jì)算較為簡(jiǎn)便，但局限性也較明顯，即缺乏對(duì)水質(zhì)總體的評(píng)價(jià)。綜合評(píng)價(jià)法主要包括模糊綜合評(píng)價(jià)法[2]、人工神經(jīng)網(wǎng)絡(luò)評(píng)價(jià)法[3]、層次分析法[4]、灰色聚類(lèi)法[5]和投影尋蹤法[6]等，但是這些方法在判斷水環(huán)境優(yōu)劣時(shí)需要與相應(yīng)的標(biāo)準(zhǔn)值進(jìn)行比較，然而很多水環(huán)境指標(biāo)沒(méi)有標(biāo)準(zhǔn)值，尤其是在漁業(yè)水質(zhì)評(píng)價(jià)時(shí)，可供利用的有標(biāo)準(zhǔn)值的指標(biāo)很少，這給水質(zhì)評(píng)價(jià)帶來(lái)很多不便。主成分分析法（PCA）也稱(chēng)為主變量分析法，旨在利用降維的思想對(duì)高維變量系統(tǒng)進(jìn)行最佳綜合和簡(jiǎn)化，將較多的指標(biāo)轉(zhuǎn)化為少數(shù)幾個(gè)綜合指標(biāo)[7]。由于該方法不需要參考指標(biāo)的標(biāo)準(zhǔn)值，因此，可以根據(jù)任何可以量化的指標(biāo)進(jìn)行評(píng)價(jià)比較，尤其適用于標(biāo)準(zhǔn)值很少的漁業(yè)水體的評(píng)價(jià)比較。本試驗(yàn)以天津市水產(chǎn)研究所淡水試驗(yàn)站的3口養(yǎng)殖池塘為研究對(duì)象，研究了池塘環(huán)境各指標(biāo)的相關(guān)性，并首次利用主成分分析法對(duì)養(yǎng)殖池塘不同階段的環(huán)境質(zhì)量進(jìn)行了比較，以期為該類(lèi)型池塘選擇合適的調(diào)控措施提供參考，并為建立漁業(yè)水體環(huán)境評(píng)價(jià)方法奠定基礎(chǔ)。

1? 材料與方法

1.1? 研究對(duì)象及樣品采集

選取天津市水產(chǎn)研究所淡水試驗(yàn)站3口養(yǎng)殖池塘，每口面積為0.267 hm2，池底土質(zhì)，池深3.0 m，有效水深2.5 m。每口池塘配備一臺(tái)1.5 kW的葉輪式增氧機(jī)和一臺(tái)120 W的投餌機(jī)。養(yǎng)殖模式為魚(yú)蝦混養(yǎng)，主養(yǎng)凡納濱對(duì)蝦，養(yǎng)殖密度為120萬(wàn)尾/hm2，放養(yǎng)規(guī)格為1.0～1.5 cm/尾，混養(yǎng)種類(lèi)主要為烏克蘭鱗鯉（125 g/尾，7 500尾/hm2）、異育銀鯽（50 g/尾，3 000尾/hm2）、草魚(yú)（250 g/尾，450尾/hm2）、鰱（125 g/尾，2 250尾/hm2）、鳙（250 g/尾，750尾/hm2）。養(yǎng)殖過(guò)程中對(duì)凡納濱對(duì)蝦和魚(yú)類(lèi)進(jìn)行分別投喂。魚(yú)類(lèi)投喂人工配合飼料，日投餌4次;凡納濱對(duì)蝦日投喂2次，上下午各一次，在魚(yú)類(lèi)集中攝食時(shí)沿池塘四周投喂。每天中午及晚上開(kāi)啟增氧機(jī)，中午開(kāi)啟時(shí)間為12：00～14：00，晚上開(kāi)啟時(shí)間根據(jù)實(shí)際情況而定，另外，高溫季節(jié)及陰雨天氣加長(zhǎng)開(kāi)啟時(shí)間。養(yǎng)殖過(guò)程中主要是通過(guò)池塘的清淤殺菌、苗種檢疫殺菌、定期投喂藥餌、水質(zhì)殺菌、底質(zhì)改良等措施預(yù)防病害的發(fā)生。整個(gè)養(yǎng)殖周期沒(méi)有從外部引水，補(bǔ)充水源主要為自然降雨。

在養(yǎng)殖季節(jié)（5—10月）每月采樣1次。每次采樣都選擇在晴天上午9：00～10：00進(jìn)行。用于主成分分析的指標(biāo)有20個(gè)。其中，水質(zhì)的測(cè)定指標(biāo)有pH、懸浮物、化學(xué)需氧量、氨氮、亞硝酸鹽氮、硝酸鹽氮、總氮、活性磷酸鹽、總磷、銅、鋅、總硬度、總堿度、鹽度;沉積物的測(cè)定指標(biāo)有總磷、總氮、銅、鋅、硫化物、有機(jī)碳。各測(cè)定指標(biāo)依據(jù)相應(yīng)國(guó)家標(biāo)準(zhǔn)進(jìn)行測(cè)定。

1.2? 研究方法

主成分分析法（PCA）也稱(chēng)為主變量分析法，旨在利用降維的思想對(duì)高維變量系統(tǒng)進(jìn)行最佳綜合和簡(jiǎn)化，將較多的指標(biāo)轉(zhuǎn)化為少數(shù)幾個(gè)綜合指標(biāo)。主成分個(gè)數(shù)確定的原則：①只選取特征根>1時(shí)對(duì)應(yīng)的主成分;②累積貢獻(xiàn)率達(dá)到80%以上對(duì)應(yīng)的主成分;③根據(jù)特征根變化的突變點(diǎn)決定主成分變量。

用SPSS13.0對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行主成分分析。

2? 結(jié)果與分析

2.1? 水質(zhì)、沉積物指標(biāo)的相關(guān)性

用相關(guān)分析研究了水質(zhì)和沉積物各因子之間的相關(guān)關(guān)系，結(jié)果顯示，沉積物銅含量、水體活性磷酸鹽含量和總堿度、總硬度與5個(gè)指標(biāo)存在顯著相關(guān)性，水體pH、氨氮含量、鋅含量及總磷與4個(gè)指標(biāo)存在顯著相關(guān)性，沉積物鋅含量、有機(jī)碳含量和水體化學(xué)耗氧量與3個(gè)指標(biāo)存在顯著相關(guān)性，底泥總氮含量和水體總硬度與2個(gè)指標(biāo)存在顯著相關(guān)性，水體亞硝酸鹽氮含量與1個(gè)指標(biāo)存在相關(guān)性。由此可見(jiàn)，各指標(biāo)之間存在著錯(cuò)綜復(fù)雜的關(guān)系。

以沉積物總氮為自變量，以其他全部指標(biāo)為因變量通過(guò)逐步回歸分析方法建立了回歸方程。由表1可知，沉積物中總氮（TN沉積物）與水質(zhì)鋅（Zn+）、亞硝酸鹽氮（NO2-）、懸浮物（SS）的回歸方程為T(mén)N沉積物 =-7.39 Zn++0.435 NO2--2.523 SS。

2.2? 主成分分析計(jì)算

由各指標(biāo)兩兩相關(guān)分析可知，各指標(biāo)之間存在著復(fù)雜的相關(guān)關(guān)系，說(shuō)明這些指標(biāo)所反映的池塘的環(huán)境信息存在重疊關(guān)系，適合進(jìn)行主成分分析。主成分分析的特征值和主成分方差貢獻(xiàn)率及累計(jì)方差貢獻(xiàn)率和特征值曲線如表2和圖1所示，由于前5個(gè)主成分累計(jì)方差貢獻(xiàn)率已達(dá)到100.000%，故其他主成分未列出。根據(jù)主成分提取原則，主成分1、2、3、4的特征值都大于1，并且4個(gè)主成分累計(jì)方差貢獻(xiàn)率為95.893%，從圖1也可以看出，特征值曲線在主成分4出現(xiàn)了突變，由此可知，將此4個(gè)主成分作為綜合指標(biāo)，可較為全面地反映濱海型鹽堿地池塘的環(huán)境狀況。因此取主成分的個(gè)數(shù)為4。

主成分載荷矩陣代表的是各指標(biāo)與主成分之間的相關(guān)關(guān)系，某一指標(biāo)與某一主成分的載荷絕對(duì)值越大，說(shuō)明該指標(biāo)與該主成分的相關(guān)性越強(qiáng)，聯(lián)系越緊密。從表3可以看出，氨氮、活性磷酸鹽、總堿度、總磷等指標(biāo)在主成分1上的載荷較大，有機(jī)碳dn、pH、亞硝酸鹽氮等指標(biāo)在主成分2上的載荷較大，硫化物dn、總氮dn、懸浮物、鹽度、總硬度等指標(biāo)在主成分3上的載荷較大，鋅dn、鋅等指標(biāo)在主成分4上的載荷較大。從主成分載荷的分布可看出，主成分1可解釋濱海型鹽堿地池塘水體的氮、磷水平，主成分2可解釋濱海型鹽堿地池塘沉積物的有機(jī)碳水平和水體的亞硝酸鹽氮水平，主成分3可解釋濱海型鹽堿地池塘水體的懸浮物水平，主成分4可解釋濱海型鹽堿地池塘的重金屬鋅水平。

由表4水質(zhì)指標(biāo)的載荷矩陣可得出主成分的得分公式[8]。

式中，Z為主成分載荷值;L為標(biāo)準(zhǔn)化的指標(biāo)值;F為主成分得分;i為采樣次數(shù);g為主成分個(gè)數(shù);

j為指標(biāo)個(gè)數(shù)。

主成分綜合評(píng)價(jià)根據(jù)如下公式計(jì)算[8]：

式中，λ是特征值。

2.3? 評(píng)價(jià)結(jié)果

將各監(jiān)測(cè)批次的標(biāo)準(zhǔn)化數(shù)據(jù)代入主成分公式中，便可得到各批次的得分，說(shuō)明了不同污染物對(duì)各監(jiān)測(cè)批次水體的影響程度。再利用主成分綜合評(píng)價(jià)公式得出各監(jiān)測(cè)批次的綜合得分，主成分綜合得分能夠反映水體的水質(zhì)狀況，得分越高，水質(zhì)污染狀況越嚴(yán)重。從主成分綜合得分（表4）可以看出，在養(yǎng)殖前期水質(zhì)最差，其主成分綜合評(píng)分最高值達(dá)到4.48，養(yǎng)殖中期7月、8月主成分綜合得分較低，在7月29日達(dá)到最低值-5.56，水質(zhì)處于最好狀態(tài)，養(yǎng)殖后期水質(zhì)又有所下降。

3? 討論

從主成分綜合得分可以看出，在剛開(kāi)始養(yǎng)殖時(shí)池塘環(huán)境質(zhì)量最差，這是因?yàn)樘旖虻臑I海型鹽堿地池塘每年在放苗前和剛放苗時(shí)池塘的枝角類(lèi)生物量都很大，而浮游植物的密度都很低，池塘水體氨氮和亞硝酸鹽氮含量都很高。隨著所放魚(yú)苗和蝦苗的攝食，枝角類(lèi)生物量急劇減少，浮游植物種類(lèi)和生物量逐漸增多，池塘水體質(zhì)量明顯改善，整個(gè)池塘生態(tài)系統(tǒng)達(dá)到一個(gè)較好的平衡狀態(tài)。在本研究中，養(yǎng)殖中后期池塘環(huán)境質(zhì)量最好，這可能是因?yàn)楸狙芯恐兄黟B(yǎng)對(duì)象是凡納濱對(duì)蝦，并且套養(yǎng)了一定比例的鰱鱅等濾食性魚(yú)類(lèi)，池塘養(yǎng)殖結(jié)構(gòu)比較合理，養(yǎng)殖過(guò)程中管理比較規(guī)范。利用主成分分析法對(duì)濱海型鹽堿地養(yǎng)殖池塘環(huán)境的評(píng)價(jià)結(jié)果與池塘水質(zhì)、沉積物和浮游生物的分析結(jié)果基本一致，說(shuō)明主成分分析法是一種比較理想的養(yǎng)殖池塘環(huán)境評(píng)價(jià)方法。主成分分析法已廣泛應(yīng)用于河流湖泊等水環(huán)境質(zhì)量評(píng)價(jià)[9，10]，本方法計(jì)算簡(jiǎn)便，不需要復(fù)雜的計(jì)算公式，用成熟的統(tǒng)計(jì)軟件就可以實(shí)現(xiàn)，因此有很好的實(shí)用性。筆者之前發(fā)表的文章也可以實(shí)現(xiàn)任意指標(biāo)的比較，但是計(jì)算方法較復(fù)雜[11]，因此主成分分析法在漁業(yè)環(huán)境評(píng)價(jià)中有更廣闊的應(yīng)用前景。

參考文獻(xiàn)：

[1] 楊? 旭，曾祥亮.基于單因子評(píng)價(jià)法和污染指數(shù)法的鄭州大學(xué)眉湖水質(zhì)評(píng)價(jià)[J].江蘇科技信息，2014（5）：51-53.

[2] 宋? 科，趙? 晟，張? 力，等.基于模糊綜合評(píng)價(jià)法的東極大黃魚(yú)（Larimichthys crocea）養(yǎng)殖區(qū)海水水質(zhì)評(píng)價(jià)[J].海洋與湖沼，2013， 44（2）：383-388.

[3] 王彤彤，張? 劍，涂? 川，等.IPSO-BP神經(jīng)網(wǎng)絡(luò)在渭河天水段水質(zhì)評(píng)價(jià)中的應(yīng)用[J].環(huán)境科學(xué)與技術(shù)，2013，36（8）：175-181.

[4] 盧文喜，李? 迪，張? 蕾.基于層次分析法的模糊綜合評(píng)價(jià)在水質(zhì)評(píng)價(jià)中的應(yīng)用[J].節(jié)水灌溉，2011（3）：46-49.

[5] 趙天燕，賀方升，侯棋棕，等.一種改進(jìn)的灰關(guān)聯(lián)分析法在水質(zhì)評(píng)價(jià)中的應(yīng)用[J].安全與環(huán)境學(xué)報(bào)，2005，5（1）：26-29.

[6] 康? 明，王麗萍，趙璧奎，等.基于投影尋蹤動(dòng)態(tài)聚類(lèi)法的水庫(kù)水質(zhì)評(píng)價(jià)模型[J].水力發(fā)電，2013，39（1）：16-19.

[7] 邢? 靜，張? 越，陳彥磊，等.基于主成分分析法的黃河流域國(guó)控?cái)嗝嫠|(zhì)評(píng)價(jià)[J].節(jié)水灌溉，2013（10）：31-34.

[8] 劉? 瀟，薛? 瑩，紀(jì)毓鵬，等.基于主成分分析法的黃河口及其鄰近水域水質(zhì)評(píng)價(jià)[J].中國(guó)環(huán)境科學(xué)，2015，35（10）：3187-3192.

[9] 龔? 然，徐力剛，徐? 進(jìn)，等.主成分分析法在城市小型湖泊水質(zhì)評(píng)價(jià)中的應(yīng)用[J].環(huán)境科學(xué)與技術(shù)，2014，37（S2）：486-488，545.

[10] 陳? 榮，李雄華.主因子分析法在東江湖水質(zhì)評(píng)價(jià)中的應(yīng)用[J].湖南文理學(xué)院學(xué)報(bào)（自然科學(xué)版），2015，27（4）：64-68.

[11] 吳會(huì)民，姜巨峰，李春艷，等.三種不同養(yǎng)殖品種池塘水質(zhì)變化特征分析及比較[J].淡水漁業(yè)，2016，46（2）：82-86.