基于數(shù)值分析的藻類生長監(jiān)測因子關(guān)聯(lián)性研究

周斌，楊利航，張曉

（1.浙江省環(huán)境監(jiān)測中心，浙江杭州310012；2.浙江環(huán)境監(jiān)測工程有限公司，浙江杭州310012）

基于數(shù)值分析的藻類生長監(jiān)測因子關(guān)聯(lián)性研究

周斌1，2，楊利航1，2，張曉1，2

（1.浙江省環(huán)境監(jiān)測中心，浙江杭州310012；2.浙江環(huán)境監(jiān)測工程有限公司，浙江杭州310012）

連續(xù)監(jiān)測湖州市與太湖入湖口交界斷面的10項水質(zhì)指標(biāo)，分別采用Pearson相關(guān)分析法和聚類分析法分析藻類暴發(fā)期和非暴發(fā)期各項水質(zhì)指標(biāo)之間的關(guān)聯(lián)性。兩種方法分析結(jié)果表明：在藻類非暴發(fā)期，藻類密度與葉綠素a、CODMn、總磷呈正相關(guān)關(guān)系，與氮磷比呈負(fù)相關(guān)關(guān)系；在藻類暴發(fā)期，藻類密度與葉綠素a、CODMn、pH、DO呈顯著正相關(guān)關(guān)系；與氨氮、總氮、氮磷比呈顯著負(fù)相關(guān)關(guān)系。通過分析不同時期各藻類生長監(jiān)測因子間的關(guān)聯(lián)性，得出不同時期藻類預(yù)警監(jiān)測應(yīng)重點關(guān)注的監(jiān)測指標(biāo)，為太湖流域藻類預(yù)警監(jiān)測與防控工作提供技術(shù)參考。

藻類；監(jiān)測因子；聚類分析；相關(guān)性

0 前言

近年來，水體環(huán)境污染造成的水體富營養(yǎng)化程度日益嚴(yán)重，水體的富營養(yǎng)化導(dǎo)致了藻類的異常增殖，尤其是由水體富營養(yǎng)化所引起的藍(lán)藻水華現(xiàn)象已成為國內(nèi)外重點關(guān)注的環(huán)境問題。藻類的暴發(fā)不僅破壞水生生態(tài)系統(tǒng)的平衡，還嚴(yán)重影響了飲用水源地的用水安全。太湖是長江三角洲地區(qū)社會經(jīng)濟(jì)發(fā)展的重要水源，由于人口密度大，社會經(jīng)濟(jì)活動頻繁，受生產(chǎn)、生活污水的不斷排入影響，導(dǎo)致太湖自身污染凈化能力有所降低，加之年均氣溫高，降雨豐沛，藻類生物量普遍偏高，存在發(fā)生較大藻類水華的風(fēng)險[1-2]。2007年，太湖藍(lán)藻大面積暴發(fā)，無錫全城自來水?dāng)喙１M管經(jīng)過近十年的治理，太湖流域水環(huán)境有了明顯改善，但太湖的生態(tài)現(xiàn)狀仍然脆弱，每年隨著夏季高溫的到來，仍時有大面積的藍(lán)藻暴發(fā)。因此，太湖水污染防治一直是各級政府水環(huán)境綜合治理規(guī)劃的重點。浙江省湖州市位于太湖南岸，市內(nèi)河網(wǎng)水系與太湖交換密切，其與太湖交界監(jiān)測斷面是太湖藻類預(yù)警應(yīng)急監(jiān)測的重要點位。研究湖州與太湖交界斷面水質(zhì)變化情況，了解藻類暴發(fā)期與非暴發(fā)期間藻類生物量與各項水質(zhì)因子之間的關(guān)聯(lián)性，找出藻類預(yù)警應(yīng)急監(jiān)測期間應(yīng)重點關(guān)注的水質(zhì)監(jiān)測指標(biāo)，可以為浙江省太湖流域藻類暴發(fā)防治監(jiān)測工作提供技術(shù)參考。

1 對象與方法

1.1 采樣點布設(shè)

浙江省湖州市的新塘港、新港口、小梅口、城西大橋、港湖大橋5個監(jiān)測點位為浙江省太湖流域藻類預(yù)警應(yīng)急監(jiān)測的關(guān)鍵點位，在夏秋季太湖藻類暴發(fā)期，受太湖水倒灌等影響，藍(lán)藻積聚現(xiàn)象發(fā)生頻繁，將采樣點位布設(shè)在這5個監(jiān)測點位，具有代表性。

1.2 采樣頻次

于某年1～12月進(jìn)行一年的連續(xù)定位采樣，按照歷年的藻類暴發(fā)時間把6月上旬至10月下旬劃分為藻類暴發(fā)期，每隔1～2 d在設(shè)定采樣點進(jìn)行采樣1次；其余時間劃分為非藻類暴發(fā)期，每周進(jìn)行采樣1次。

1.3 分析項目與方法

測定分析10個水質(zhì)因子，每個樣品每個因子重復(fù)3次測量。其中10個因子包括理化因子與生物因子，理化因子包括：水溫、pH、溶解氧（DO）、總氮、氨氮、總磷、氮磷比、高錳酸鉀指數(shù)（CODMn），評價及分析方法依據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》（GB 3838-2002）；生物因子包括：葉綠素a（Chla）、藻類密度，藻類密度采用YSI-6600便攜式水質(zhì)監(jiān)測儀現(xiàn)場監(jiān)測，葉綠素a測定采用丙酮萃取比色法[4]。

1.4 數(shù)據(jù)結(jié)構(gòu)

將藻類密度做ln（X+1）變量轉(zhuǎn)換使其符合正態(tài)分布，再按時間順序建立該采樣點理化因子和生物因子的數(shù)據(jù)庫。

1.5 數(shù)據(jù)分析方法與原理

1.5.1 聚類分析

聚類分析（Cluster Analysis）是根據(jù)事物本身的特征研究個體分類的方法。聚類分析的基本思想認(rèn)為所研究的樣品或因子之間存在著程度不同的相似性，將相似性較大的因子聚合為一類，把另外一些相似程度較大的因子又聚合為另一類。這樣將所有因子分成許多類別，同一類中的個體有較大的相似性，不同類中的個體差異性很大。

聚類分析根據(jù)實際的需要可分為兩類：按樣品進(jìn)行聚類分析，稱為Q型聚類法；按因子進(jìn)行聚類分析，稱為R型聚類法。本研究是對水質(zhì)因子進(jìn)行分類，采用的是R型聚類法。

1.5.2 數(shù)據(jù)的標(biāo)準(zhǔn)化

試驗測定的水質(zhì)因子的量綱不同且數(shù)量級相差很大，為了去除因子的量綱對計算的影響，首先對因子的原始數(shù)據(jù)進(jìn)行了標(biāo)準(zhǔn)化，采用的標(biāo)準(zhǔn)化公式如下：

式中y1為因子X的標(biāo)準(zhǔn)化數(shù)值；Xi為因子X的第i個樣品值；n為樣品數(shù)目；X為因子X的平均值，以下在聚類分析提供的數(shù)據(jù)都是標(biāo)準(zhǔn)化后的數(shù)據(jù)。

1.5.3 聚類分析中的距離、相似系數(shù)

假設(shè)有n個樣品，每個樣品測得p項因子，原始資料陣為

其中Xij（i=1，…n；j=1，…；p）為第i個樣品的第j個因子的觀測數(shù)據(jù)。第i個樣品Xi為矩陣X的第i行所描述，所以任何兩個樣品XK與XL之間的相似性，可以通過矩陣X中的第K行與第L行的相似程度來判斷；任何兩個變量XK與XL之間的相似性，可以通過第K列與第L列的相似程度來判斷。

p個因子之間相似性的定義與樣品相似性定義類似，但此時是在n維空間中來研究的，變量之間的相似性是通過原始資料矩陣X中p列間相似關(guān)系來研究的。

這時，每個因子都可以看成是m維空間的一個點，n個因子就組成m維空間的n個點，用各點之間的距離與相似系數(shù)來判斷各因子之間的靠近程度。

設(shè)dpq為任意兩個因子p、q之間的距離，把因子xp和xq之間的距離簡記為dpq。歐式距離（Euclidean distance，即二階Minkowski距離）dpq的計算公式如下：

1.5.4 離差平方和系統(tǒng)聚類法

選用常用的系統(tǒng)聚類法，令n個樣品自成一類，計算出相似性測度，此時類間距離與樣品間距離是等價的，把測度最小的兩個類合并；然后按照某種聚類方法計算類間的距離，再按最小距離準(zhǔn)則并類；這樣每次減少一類，持續(xù)下去直到所有樣品都?xì)w為一類為止。離差平方和法（Ward’s method）是系統(tǒng)聚類法的一種，是沃爾德（Ward，1963）首先提出的，其基本思想是同一類內(nèi)各案例的離差平方和應(yīng)該較小，不同類的離差平方和應(yīng)該較大。求解過程是首先使每個案例自成一類，每一步使離差平方和增加，最小的兩類合并為一類，直到所有的案例合并為一類為止。

用G表示類，XG為X中元素的均值，則G的直徑為：

用Dp，Dq分別表示Gp和Gq的直徑，Dp+q用DpYDq表示大類Dw（p，q）的直徑，Dw（p，q）表示Gp和Gq的距離，則：

2 分析與結(jié)論

2.1 藻類非暴發(fā)期水質(zhì)因子的線性相關(guān)分析

采用Pearson相關(guān)分析法對藻類非暴發(fā)期監(jiān)測的10個水質(zhì)因子進(jìn)行分析，得到相關(guān)系數(shù)與伴隨系數(shù)（P值）矩陣表。經(jīng)計算可得Pearson相關(guān)系數(shù)矩陣，為10×10階對稱矩陣。由于矩陣為對稱矩陣，所以只列出它的上三角部分，詳見表1。

2.1.1 藻類密度與相關(guān)因子的關(guān)系

目前我國尚無慢阻肺首次確診前錯失早期診斷時間的研究報道。因此本研究重點驗證慢阻肺確診前錯失早期診斷時間是否可以反映慢阻肺的嚴(yán)重程度，探討與疾病嚴(yán)重程度相關(guān)的臨床危險因素，從而進(jìn)一步探索建立有利于慢阻肺早期確診的方法。

藻類密度是反映水體中浮游植物量最重要的綜合因子，可以較直觀地描述太湖水體富營養(yǎng)的情況，故選它作為指示因子。由表1數(shù)據(jù)分析可得，在藻類非暴發(fā)期，藻類密度與葉綠素a、CODMn、總磷之間均存在顯著正相關(guān)關(guān)系，與氮磷比呈顯著負(fù)相關(guān)關(guān)系，但相關(guān)性程度不高（相關(guān)系數(shù)均小于0.5）。

表1 藻類非暴發(fā)期各水質(zhì)指標(biāo)Pearson相關(guān)系數(shù)矩陣Table 1 The Pearson correlation coefficient matrix of non-outbreak period ofalgae

2.2 藻類非暴發(fā)期水質(zhì)因子的聚類分析

表2為對非暴發(fā)期10個水質(zhì)因子進(jìn)行聚類分析的詳細(xì)過程，合并10個水質(zhì)因子一共用了9步。圖1為反映聚類全過程的樹形圖。聚類分析結(jié)果反映出太湖水質(zhì)環(huán)境的重要特征：在藻類非暴發(fā)期，10項因子總體可以分為四類：總氮、氮磷比、氨氮；水溫、pH；總磷、CODMn、DO；葉綠素a、藻類密度。每類之間因子的相關(guān)關(guān)系較大，且有相同的變化趨勢。

從圖1橫坐標(biāo)的垂直距離可以看出，藻類密度與葉綠素a的相關(guān)性最大，這與用Pearson相關(guān)分析的結(jié)果一致。葉綠素a通常在水質(zhì)監(jiān)測中作為藻類生物量的代表性參數(shù)，并且有些研究直接以葉綠素a值直接表征藻類暴發(fā)情況。藻類在非暴發(fā)時期基本處于休眠的狀態(tài)，藻類的生長對湖體的水文參數(shù)與理化性質(zhì)的影響較小。

表2 藻類非暴發(fā)期水質(zhì)因子聚類過程Table 2 The clustering process of water quality factors in non-outbreak period of algae

2.3 藻類暴發(fā)期水質(zhì)因子的線性相關(guān)分析

采用Pearson相關(guān)分析法對藻類暴發(fā)期監(jiān)測的10個水質(zhì)因子進(jìn)行分析，得到相關(guān)系數(shù)與伴隨系數(shù)（P值）矩陣表。經(jīng)計算可得Pearson相關(guān)系數(shù)矩陣，為10×10階對稱矩陣。由于矩陣為對稱矩陣，所以只列出它的上三角部分，見表3。

2.3.1 藻類密度與相關(guān)因子的關(guān)系

2.4 藻類暴發(fā)期水質(zhì)因子的聚類分析

表4為藻類暴發(fā)期10個水質(zhì)因子進(jìn)行聚類分析的詳細(xì)過程，合并10個水質(zhì)因子一共用了9步。圖2為反映藻類暴發(fā)期聚類全過程的樹形圖。聚類結(jié)果反映出太湖水質(zhì)環(huán)境的重要特征：在藻類暴發(fā)期，10項水質(zhì)因子總體可以分為四類：總氮、氮磷比、水溫；氨氮、總磷；CODMn、藻類密度、葉綠素a、pH、DO。每類因子之間的相關(guān)關(guān)系較大，且有相同的變化趨勢。從圖2橫坐標(biāo)的垂直距離與Pearson相關(guān)性分析可以看出，藻類密度與CODMn、葉綠素a的相關(guān)性最大，其次是pH、DO，且這些因子均與藻類密度呈顯著正相關(guān)關(guān)系；另外，藻類密度與總氮、氮磷比、氨氮相關(guān)性大且呈顯著負(fù)相關(guān)。

表3 藻類暴發(fā)期Pearson相關(guān)系數(shù)矩陣Table 3 The Pearson correlation coefficient matrix of outbreak period of algae

分析結(jié)果表明藻類密度、CODMn與葉綠素a之間呈現(xiàn)顯著的正相關(guān)關(guān)系。高錳酸鹽指數(shù)（CODMn）主要是指水體中能被高錳酸鉀氧化的有機(jī)物含量。在藻類暴發(fā)期，太湖藻類的生物量主要呈現(xiàn)不斷增長的趨勢，藻類聚積量會為水體提供有機(jī)質(zhì)源，因此藻類密度與有機(jī)污染物的含量變化有同步特點，在藻類生物量高的水體中，可以被氧化的有機(jī)物量也高[5]。葉綠素是藻類的重要組成成分之一，葉綠素a的含量與水體中藻類的種類、數(shù)量等密切相關(guān)，通過測定水體中葉綠素a含量能在一定程度上反映水體的藻類污染狀況[6]。

Pearson相關(guān)性和聚類分析表明藻類密度與pH、溶解氧相關(guān)性較大，且呈顯著正相關(guān)關(guān)系，這與張澎浪，鄭曉紅等人[7－8]的研究結(jié)果一致。在地表水pH范圍內(nèi)，碳酸鹽體系中的碳主要以HCO3-形式存在，結(jié)合水體碳酸鹽的化學(xué)平衡與藻類光合作用原理，在藻類密度較高的水體中，藻類的光合作用向水中釋放了O2和OH-，從而使水中的溶解氧濃度和pH值升高。這說明該港口水體的pH、溶解氧受藻類的大量生長影響較大，在藻類防治預(yù)警監(jiān)測工作中，建議對水體的pH值和溶解氧指標(biāo)的變化趨勢進(jìn)行關(guān)注。

聚類分析結(jié)果顯示，藻類密度的增大伴隨著水體總氮與氨氮含量的減少、氮磷比數(shù)值的降低。湖庫浮游藻類的生長是由營養(yǎng)鹽、水文氣象條件等因素控制的，氮、磷營養(yǎng)鹽類是藻類生長的限制因子[9]。藻類的生長、繁殖需要吸收水體中大量的可溶性磷鹽與氮鹽，從而導(dǎo)致水體這些營養(yǎng)鹽含量的減少。太湖四周企業(yè)與農(nóng)田較多，地表特別是農(nóng)田地區(qū)受暴雨徑流、灌溉、河流沖刷的作用會產(chǎn)生土壤流失，農(nóng)田流失的土壤攜帶有大量磷、氮化合物，使得水體富含磷、氮營養(yǎng)鹽類，即通常所稱非點源污染。湖泊水體中的浮游植物主要利用氨氮和硝態(tài)氮進(jìn)行生產(chǎn)，而由于生物化學(xué)的原因氨氮更容易被藻類吸收。另外分析結(jié)果表明：在該港口藻類暴發(fā)期，磷營養(yǎng)鹽對藻類的影響小于氮對藻類生長的影響。

表4 藻類暴發(fā)期水質(zhì)因子聚類過程Table 4 The clustering process of water quality factors in the outbreak period of algae

圖2 藻類爆發(fā)期聚類樹狀圖Fig.2 The clustering dendrogram of outbreak period ofalgae

3 結(jié)論

（1）在藻類非暴發(fā)期，藻類密度與葉綠素a、CODMn、總磷之間均存在顯著正相關(guān)關(guān)系，與氮磷比呈顯著負(fù)相關(guān)關(guān)系，但相關(guān)性程度不高（相關(guān)系數(shù)均小于0.5）；其他因子與藻類密度的關(guān)聯(lián)性不顯著。非暴發(fā)期藻類的生長對湖體的水文參數(shù)與理化性質(zhì)的影響較小，建議該時期對藻類密度、葉綠素a、高錳酸鹽指數(shù)、總磷等監(jiān)測指標(biāo)進(jìn)行關(guān)注。

（2）在藻類暴發(fā)期，藻類密度除了與葉綠素a、CODMn關(guān)系密切，與水體的pH、溶解氧濃度也呈現(xiàn)顯著的正相關(guān)關(guān)系；與總氮、氮磷比、氨氮呈顯著負(fù)相關(guān)關(guān)系；但是與藻類非爆發(fā)期相比，與總磷的相關(guān)性不大。因此，我們可以在太湖流域藻類暴發(fā)預(yù)警監(jiān)測工作中，通過監(jiān)測水體中的藻類密度、pH值、溶解氧濃度、葉綠素a、CODMn的變化，對藻類暴發(fā)的趨勢進(jìn)行預(yù)測和風(fēng)險評估，以便于提前采取相應(yīng)的有效措施，破壞藻類生長能力，緩解藻類的暴發(fā)，防止水質(zhì)的進(jìn)一步惡化。

（3）在太湖流域藻類防治工作中，營養(yǎng)鹽的削減仍是最可行的控制途徑。分析結(jié)果表明該區(qū)域湖體流域氮營養(yǎng)鹽對藻類生長的影響大于磷的影響，可能與該處的氮污染嚴(yán)重程度有關(guān)，政府部門與其它相關(guān)職能部門需注意該港口湖體流域污水處理廠改造，提高流域脫氮能力；減少源頭氮等污染物排放，控制農(nóng)田人工復(fù)合肥料量使用。

參考文獻(xiàn)：

[1]吳攀，鄧建明，秦伯強(qiáng)，等.水溫和營養(yǎng)鹽增加對太湖冬、春季節(jié)藻類生長的影響[J].環(huán)境科學(xué)研究，2013，26（10）：1064-1071.

[2]朱廣偉，金穎薇，任杰，等.太湖流域水庫型水源地硅藻水華發(fā)生特征及對策分析[J].湖泊科學(xué)，2016，28（1）：9-21.

[3]于海燕，周斌，胡尊英.浙江省太湖交界斷面藻類預(yù)警應(yīng)急監(jiān)測指標(biāo)特征分析[J].環(huán)境污染與防治，2010，32（2）：59-62.

[4]國家環(huán)境保護(hù)總局，水和廢水監(jiān)測分析方法編委會.水和廢水監(jiān)測分析方法[M].第4版（增補(bǔ)版）.北京：中國環(huán)境科學(xué)出版社，2008：701-705，722.

[5]張巍，王學(xué)軍，江耀慈，等.太湖水質(zhì)指標(biāo)相關(guān)性與富營養(yǎng)化特征分析[J].環(huán)境污染與防治，2002，24（1）：50-53.

[6]金相燦.中國湖泊環(huán)境[M].北京：海洋出版社，1995.

[7]張澎浪，孫承軍.地表水體中藻類的生長對pH值及溶解氧含量的影響[J].中國環(huán)境監(jiān)測，2004，20（4）：49-50.

[8]鄭曉紅.地表水中總磷和總氮對藻類生長的影響以及藻類生長對pH值和溶解氧含量的影響[J].儀器儀表與分析監(jiān)測，2012，（3）：43-45.

[9]申哲民，張濤，馬晶，等.富營養(yǎng)化與溫度因素對太湖藻類生長的影響研究[J].環(huán)境監(jiān)控與預(yù)警，2011，3（2）：1-4.

Study on Correlations of Algae Growth Monitoring Factors based on Numerical Analysis

ZHOU Bin1，2，YANG Li-hang1，2ZHANG Xiao1，2
（1.Zhejiang Province Environmental Monitoring Center，Hangzhou，Zhejiang 310012，China; 2.Zhejiang Environmental Monitoring Engineering Co.，Ltd.，Hangzhou，Zhejiang 310012，China）

10 water quality indexes were continuously monitored at the cross-section of Tai Lake and ports in Huzhou，and the correlations among the 10 indexes were analyzed by Pearson correlation analysis and cluster analysis respectively.The results showed that the algal density was positively correlated with chlorophyll-a，CODMnand totalphosphorus（TP）in the non-outbreak period of algae，and negatively correlated with the ratio of nitrogen to phosphorus.In the outbreak period of algae，the algal density has a significant positive correlation with chlorophyll-a，CODMn，pH，DO，and a negative correlation with ammonia nitrogen（NH3-N），total nitrogen（TN）and the ratio of nitrogen to phosphorus（N/P）.The focus indexes in different periods could be obtained from the analysis results，which could provide technical reference for the algae early-warning and emergency monitoring in Tai Lake.

algae;monitoring factors;cluster analysis;correlation

1006-4184（2017）8-0037-07

浙江力普石墨粉碎球形化生產(chǎn)線點“墨”成金

（丁文供稿）

2017-05-02

浙江省環(huán)?？蒲许椖浚?016A011）。

周斌（1985-），男，浙江杭州人，碩士，工程師，主要從事環(huán)境監(jiān)測方面的工作研究。E-mail:34299660@qq.com。

石墨銷售價格是按照顆粒大小和純度來決定的，顆粒越大、純度越高的石墨價格越貴，所以選用專業(yè)加工設(shè)備至關(guān)重要。由于我國并不掌握石墨深加工核心技術(shù)，而先進(jìn)石墨加工技術(shù)被美、日、歐盟等壟斷，導(dǎo)致我國石墨資源“低出高進(jìn)”，外國“以購代采”的狀態(tài)長期存在。如何改變這一現(xiàn)狀，將石墨資源優(yōu)勢轉(zhuǎn)化為經(jīng)濟(jì)優(yōu)勢？浙江力普粉碎設(shè)備有限公司一馬當(dāng)先，發(fā)揮自身長期從事粉碎設(shè)備行業(yè)累積的經(jīng)驗優(yōu)勢，根據(jù)用戶實際需求研制的“GCL系列石墨粉碎球化生產(chǎn)線”（省級新產(chǎn)品試制計劃項目，編號：2016D60SA642837）為國內(nèi)的石墨粉碎深加工行業(yè)解決了一大難題，成為行業(yè)的佼佼者。該生產(chǎn)線采用系列微粉機(jī)和大小渦輪分級機(jī)的優(yōu)化配置，解決了生產(chǎn)線上下游產(chǎn)量的銜接問題。研制的負(fù)壓除塵回收系統(tǒng)，采用自動加料系統(tǒng)和風(fēng)網(wǎng)系統(tǒng)，將生產(chǎn)線中的全部尾粉在出口統(tǒng)一收集，提高生產(chǎn)線的自動化程度，防止粉塵污染。日前，這一成果在杭州通過省級新產(chǎn)品鑒定，專家認(rèn)為該技術(shù)處國內(nèi)領(lǐng)先水平。