李 哲 卿，張 紅 艷，許 景 璇，白 凱 華，代 俊 峰，4

(1.桂林理工大學(xué) 環(huán)境科學(xué)與工程學(xué)院，廣西 桂林 541004； 2.廣西巖溶地區(qū)水污染控制與用水安全保障協(xié)同創(chuàng)新中心，廣西 桂林 541004； 3.桂林水文中心，廣西 桂林 541000； 4.廣西環(huán)境污染控制理論與技術(shù)重點(diǎn)實(shí)驗(yàn)室，廣西 桂林 541004)

0 引 言

河流生態(tài)環(huán)境日益引起社會(huì)的關(guān)注，河流水量水質(zhì)的變化規(guī)律分析對(duì)于河流生態(tài)環(huán)境保護(hù)十分重要。目前，對(duì)河流水量、水質(zhì)變化的單一研究較多，水量水質(zhì)綜合變化分析也逐漸引起重視。孟瑩瑩等[1]在對(duì)上海市3個(gè)典型分流制雨水系統(tǒng)的旱流污水(合流制排水系統(tǒng)在晴天時(shí)輸送的污水)水量水質(zhì)進(jìn)行綜合分析之后，再結(jié)合流程圖法、CMBM法在雨水系統(tǒng)旱流溯源中的應(yīng)用，得到了理想結(jié)果。Cao等[2]將水量水質(zhì)轉(zhuǎn)化為評(píng)價(jià)指標(biāo)，提出線性可加性指標(biāo)，結(jié)合水量水質(zhì)進(jìn)行了水資源量的計(jì)算。然而，豐枯不同水期的河流水量水質(zhì)綜合變化研究相對(duì)較少。在豐水期，徑流增加對(duì)河流水質(zhì)的影響包括兩個(gè)方面：一是增大河流的污染物稀釋能力和自凈能力，二是增大流域侵蝕能力，可能帶來更多的污染物進(jìn)入河流。這兩方面影響的綜合結(jié)果如何，在不同地區(qū)和不同河流可能呈現(xiàn)不同的結(jié)果。

靈渠作為漓江的重要支流，溝通了長(zhǎng)江支流湘江和珠江支流桂江(漓江)，在漓江流域的發(fā)展中發(fā)揮了重要作用，2018年靈渠成功入選世界灌溉工程遺產(chǎn)名錄。由于經(jīng)濟(jì)社會(huì)的快速發(fā)展，靈渠周圍村鎮(zhèn)的生活污水和農(nóng)業(yè)生產(chǎn)污水的排放，給靈渠水質(zhì)帶來了污染風(fēng)險(xiǎn)。處于漓江流域的靈渠，水資源相對(duì)豐富，但豐枯季節(jié)變異大，豐枯水期對(duì)靈渠水質(zhì)的影響如何，值得深入研究。本文基于收集的長(zhǎng)序列水量水質(zhì)數(shù)據(jù)，分析豐枯水期靈渠水量水質(zhì)的變化及其關(guān)系，以為靈渠水環(huán)境保護(hù)和生態(tài)環(huán)境改善提供依據(jù)。

1 資料與方法

1.1 研究區(qū)概況

靈渠位于漓江流域上游(見圖1)，屬于亞熱帶季風(fēng)氣候，降雨量的年際年內(nèi)分配受季風(fēng)活動(dòng)的影響較大，雨季的降雨量高達(dá)1 100～1 400 mm，而干季的降雨量?jī)H為370～430 mm，年平均蒸發(fā)量約為1 260.8～1 792.2 mm。本次研究選取了靈渠水文站1960～2019年徑流數(shù)據(jù)和2004～2019年的水質(zhì)數(shù)據(jù)，進(jìn)行靈渠豐枯水期水量水質(zhì)的變化分析。

圖1 靈渠水文站位置Fig.1 The location of hydrological station in Lingqu

1.2 研究方法

采用統(tǒng)計(jì)分析法、Mann-Kendall法和Spearman秩相關(guān)系數(shù)法分析靈渠水量和水質(zhì)的變化特點(diǎn)。

(1) 統(tǒng)計(jì)分析法，用來分析徑流的年內(nèi)與年際分布特征，主要是計(jì)算出相應(yīng)的參數(shù)，年內(nèi)分布特征是用不均勻系數(shù)Cn和完全調(diào)節(jié)系數(shù)[3]Cr來表示，年際分布特征是用極值比K與變差系數(shù)Cv來反映。Cn與Cr越大則說明不同月份的徑流量之間的差異就越大，則年內(nèi)分布不均勻[4]；K值與Cv越大，則年際分布越明顯[5]。

(2) Mann-Kendall方法是一種非參數(shù)檢驗(yàn)方法[6]，計(jì)算簡(jiǎn)便，不受異常數(shù)值的影響，該方法常被用來分析水文序列的變化趨勢(shì)[7]，其公式如下：

(1)

(2)

j

S方差計(jì)算公式為

(3)

式中：m，n表示年時(shí)間序列中具有相同值的變量數(shù)目；tp是第p組的相同值數(shù)目。

統(tǒng)計(jì)量Z值的公式如下：

(4)

若Z值是正值，表明該序列有上升的傾向；負(fù)值則有下降傾向[8-9]。

(3) Spearman秩相關(guān)系數(shù)是一個(gè)非參數(shù)指標(biāo)，用來衡量?jī)蓚€(gè)變量之間的相互依賴性。利用單調(diào)方程表示變量之間的相關(guān)性，若兩變量完全單調(diào)相關(guān)時(shí)，Spearman秩相關(guān)系數(shù)為1或-1，Spearman相關(guān)系數(shù)又被認(rèn)為是等級(jí)變量之間的Pearson相關(guān)系數(shù)[10-11]。

(4) 綜合污染指數(shù)是評(píng)價(jià)水環(huán)境質(zhì)量的一種重要方法[12]，具體公式如下：

(5)

(6)

式中：P是綜合污染指數(shù)；Pi是i污染物的污染指數(shù)；n是污染物的種類；Ci是i污染物的實(shí)測(cè)濃度平均值；Si是i污染物評(píng)價(jià)標(biāo)準(zhǔn)值。

2 不同水期內(nèi)的徑流、水質(zhì)指標(biāo)變化分析

2.1 徑流變化分析

本文依照靈渠流量變化將年內(nèi)時(shí)段劃分為兩個(gè)水文期：豐水期(3～8月份)與枯水期(1～2月份，9～12月份)。

2.1.1年內(nèi)變化特征

靈渠各月徑流大小差異較大(見圖2)，豐水期流量占總流量的79%，5月份出現(xiàn)最大值，占全年流量的20%；枯水期流量占總流量的21%，最小值則出現(xiàn)在12月份。

圖2 徑流年內(nèi)分配Fig.2 Annual distribution of runoff

為了方便分析，將1960～2019年徑流以10 a為時(shí)段計(jì)算Cn和Cr值(見表1)。繪制Cn，Cr曲線，見圖3～4。結(jié)果表明：枯水期的Cn，Cr均值大于豐水期，且都處于波動(dòng)上升狀態(tài)，意味著枯水期徑流量的變化比豐水期明顯。

表1 1960～2019年不同時(shí)段年內(nèi)分配結(jié)果Tab.1 Annual distribution results in different periods from 1960 to 2019

圖3 豐水期徑流年內(nèi)分配參數(shù)曲線Fig.3 Annual distribution parameter curve in wet season

圖4 枯水期徑流年內(nèi)分配參數(shù)曲線Fig.4 Annual distribution parameter curve in dry season

2.1.2年際變化特征

不同水文期的年際徑流變化如圖5所示，靈渠豐水期的徑流呈現(xiàn)上升趨勢(shì)，但幅度較緩，線性斜率k為0.02，豐水期徑流的變差系數(shù)Cv和年極值比K分別為0.60和4.30?？菟诘哪觌H徑流呈現(xiàn)下降趨勢(shì)，線性斜率k為-0.04，Cv值為0.65，K值為17.76，兩者均大于豐水期。

對(duì)豐枯水期的年徑流序列進(jìn)行非參數(shù)統(tǒng)計(jì)，1960～2019年豐枯水期的統(tǒng)計(jì)量Z值分別為0.23和-2.96，且枯水期的統(tǒng)計(jì)值通過了99%的顯著檢驗(yàn)，表明枯水期的徑流下降趨勢(shì)顯著。

圖5 不同水文期內(nèi)的徑流變化Fig.5 Runoff variation in different hydrological periods

2.2 水質(zhì)變化分析

2.2.1水質(zhì)變化趨勢(shì)

2004～2019年靈渠高錳酸鹽指數(shù)的濃度變化范圍是0.5～4.0 mg/L，2004年以外的其他年份均處于Ⅰ類標(biāo)準(zhǔn)(參照GB 3838-2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》)；氨氮濃度指標(biāo)的變化范圍是0.01～0.62 mg/L，含量基本達(dá)標(biāo)，處于Ⅰ～Ⅱ類水之間；總磷指標(biāo)變化范圍是0.01～0.5 mg/L，僅有2007年處于Ⅲ類水，其余年份均處于Ⅰ～Ⅱ類水之間，總體情況較好。

靈渠屬于Ⅲ類水質(zhì)保護(hù)標(biāo)準(zhǔn)，據(jù)此計(jì)算靈渠水質(zhì)污染綜合指數(shù)變化(見圖6)。結(jié)果表明：2004～2019年水質(zhì)狀況總體好轉(zhuǎn)，污染綜合指數(shù)呈下降趨勢(shì)，趨勢(shì)線的斜率為-0.01。

圖6 水質(zhì)污染綜合指數(shù)Fig.6 Comprehensive index of water pollution

2.2.2污染物濃度變化

應(yīng)用Mann-Kendall法計(jì)算豐枯水期污染物濃度的統(tǒng)計(jì)值Z(見表2)，不同污染物濃度豐枯水期的變化關(guān)系見圖7～8。結(jié)果顯示：豐枯水期的高錳酸鹽指數(shù)濃度統(tǒng)計(jì)值Z均為正值，但并未通過顯著檢驗(yàn)，處于不明顯的上升趨勢(shì)；氨氮、總磷濃度的統(tǒng)計(jì)值在豐枯水期都為負(fù)值，表明氨氮、總磷濃度在豐枯水期內(nèi)整體呈現(xiàn)下降趨勢(shì)，且氨氮通過了99%的顯著檢驗(yàn)。

表2 不同水期污染物濃度統(tǒng)計(jì)值ZTab.2 Statistical values Z of pollutant concentrations in different water periods

圖7 豐水期不同污染物濃度均值變化Fig.7 Variation of mean concentration of different pollutants in wet season

圖8 枯水期不同污染物濃度均值變化Fig.8 Variation of mean concentration of different pollutants in dry season

2.2.3污染物通量變化

根據(jù)資料的詳細(xì)程度，本文分析污染物通量的影響因素時(shí)，忽略了地下水交換的影響及底質(zhì)背景對(duì)污染物通量的影響，主要分析水量和污染物濃度對(duì)污染物通量的影響。污染物通量是用豐(枯)水期的月均流量乘以豐(枯)水期各月實(shí)測(cè)污染物濃度的平均值，再乘以時(shí)間轉(zhuǎn)化系數(shù)來計(jì)算的，具體公式如下：

L=KQC

(7)

式中：L是污染物通量，t；K是時(shí)間轉(zhuǎn)化系數(shù),按一個(gè)月30 d計(jì)，K=2.59×106s；Q是豐(枯)水期的月均流量,m3/s；C是豐(枯)水期各月實(shí)測(cè)污染物濃度的平均值,mg/L。

對(duì)高錳酸鹽指數(shù)、氨氮、總磷通量變化使用Mann-Kendall法進(jìn)行分析，計(jì)算豐枯水期污染物通量的統(tǒng)計(jì)值Z(見表3)，豐枯水期高錳酸鹽指數(shù)、氨氮、總磷通量變化趨勢(shì)見圖9～11。結(jié)果顯示：CODMn通量在豐水期間變化比枯水期劇烈，CODMn通量在豐水期呈上升趨勢(shì)，在枯水期呈下降趨勢(shì)。氨氮通量統(tǒng)計(jì)值Z在豐水期和枯水期都為負(fù)值，通過了99%的顯著檢驗(yàn)，說明氨氮通量在豐水期和枯水期都呈現(xiàn)明顯的下降趨勢(shì)?？偭淄吭谪S水期和枯水期都呈現(xiàn)下降趨勢(shì)，枯水期總磷通量的變化小于豐水期。

表3 不同水期污染物通量統(tǒng)計(jì)值ZTab.3 Statistical value Z of pollution flux in different water periods

圖9 CODMn通量均值變化Fig.9 Mean change of CODMn flux

圖10 氨氮通量均值變化Fig.10 Mean change of ammonia nitrogen flux

圖11 總磷通量均值變化Fig.11 Mean change of total phosphorus flux

3 不同水期內(nèi)徑流與水質(zhì)指標(biāo)的相關(guān)分析

3.1 徑流與水質(zhì)指標(biāo)濃度、污染物通量的相關(guān)關(guān)系

應(yīng)用Spearman秩相關(guān)系數(shù)法[13-14]對(duì)2004～2019年的徑流量年平均值與水質(zhì)指標(biāo)濃度、污染物通量進(jìn)行相關(guān)性分析，豐水期和枯水期徑流與水質(zhì)指標(biāo)間的相關(guān)系數(shù)見表4。結(jié)果顯示：在豐水期和枯水期間，高錳酸鹽指數(shù)濃度與徑流相關(guān)系數(shù)都為負(fù)值；氨氮、總磷與豐枯水期徑流的相關(guān)系數(shù)均為正值。豐水期、枯水期徑流與高錳酸鹽指數(shù)、氨氮、總磷通量間的相關(guān)系數(shù)均為正值，說明當(dāng)徑流量增加時(shí)，污染物通量也會(huì)隨之上升，徑流量大小對(duì)污染物通量起了決定作用。

當(dāng)處在枯水期時(shí)，徑流與高錳酸鹽指數(shù)、氨氮和總磷污染物濃度、污染物通量間的相關(guān)系數(shù)均大于豐水期，表明枯水期徑流對(duì)污染物的影響更大。

表4 徑流與污染物的相關(guān)性結(jié)果分析Tab.4 Analysis of correlation results between runoff and pollutants

3.2 河流污染物的來源分析

一般情況下，河流污染來自點(diǎn)源污染和面源污染兩個(gè)主要類型，可以通過分析水質(zhì)濃度變化與徑流間的關(guān)系來判斷河流污染的來源。當(dāng)徑流量增加，水質(zhì)濃度呈現(xiàn)上升趨勢(shì)時(shí)，此時(shí)污染物主要由面源污染引起；反之，隨著徑流量的增加，水質(zhì)濃度出現(xiàn)下降的狀況，污染物則主要由點(diǎn)源污染引起，此時(shí)的徑流表現(xiàn)出對(duì)污染物的稀釋作用[15]。

表4顯示，隨著徑流的增加，高錳酸鹽指數(shù)濃度下降，表明可氧化污染物主要來自點(diǎn)源污染；當(dāng)徑流增加時(shí)，氨氮與總磷濃度隨之上升，說明氮磷污染以面源污染為主。

3.3 靈渠污染防控對(duì)策

實(shí)地調(diào)查結(jié)果顯示，靈渠周邊沒有工業(yè)企業(yè)點(diǎn)源污染，污染來源主要是未經(jīng)深度處理的農(nóng)村生活污水，規(guī)模以下畜禽養(yǎng)殖場(chǎng)廢水排放，和水稻、柑橘等種植業(yè)肥料流失。結(jié)合本文的分析計(jì)算結(jié)果，提出靈渠的污染防控對(duì)策：

(1) 對(duì)于高錳酸鹽指數(shù)污染消減，應(yīng)加強(qiáng)對(duì)農(nóng)村生活污水等點(diǎn)源污染的治理；對(duì)于氨氮和總磷污染消減，應(yīng)加強(qiáng)對(duì)稻田、柑橘等種植業(yè)和規(guī)模以下畜禽養(yǎng)殖等面源污染的控制。

(2) 豐水期徑流增加了靈渠水量和河流自凈能力，但其也增大了流域氮磷侵蝕能力，增加了氮磷污染入河量。應(yīng)加強(qiáng)豐水期氮磷污染控制，在田間采用農(nóng)藝措施減少肥料流失，在岸邊設(shè)置植被過濾帶，在主要入河溝渠處建立人工濕地，消減豐水期氮磷污染入河量。

(3) 加強(qiáng)對(duì)靈渠生態(tài)環(huán)境保護(hù)的宣傳力度，采取生態(tài)措施加強(qiáng)農(nóng)業(yè)農(nóng)村污染治理，提升農(nóng)村生活污水治理程度，建立和完善農(nóng)業(yè)污染監(jiān)測(cè)網(wǎng)絡(luò)和監(jiān)管制度，嚴(yán)格把控有機(jī)化肥的使用[16]，提高肥料利用率。

4 結(jié) 論

(1) 漓江上游靈渠的高錳酸鹽指數(shù)濃度在豐、枯水期都呈現(xiàn)不顯著的上升趨勢(shì)，氨氮濃度在豐、枯水期都呈現(xiàn)顯著的下降趨勢(shì)，總磷濃度在豐、枯水期的變化幅度最小，呈現(xiàn)不顯著的下降趨勢(shì)。年際間高錳酸鹽指數(shù)濃度的變化比氨氮、總磷劇烈。

(2) 高錳酸鹽指數(shù)、氨氮、總磷通量在豐水期的變化均大于枯水期。氨氮、總磷通量在豐枯水期都呈現(xiàn)下降趨勢(shì)，并且氨氮通量的下降趨勢(shì)顯著。

(3) 枯水期徑流與高錳酸鹽指數(shù)、氨氮、總磷濃度和污染物通量的相關(guān)性，大于豐水期徑流與污染物濃度和污染物通量的相關(guān)性。豐、枯水期徑流與污染物通量之間的相關(guān)系數(shù)絕對(duì)值大于相應(yīng)水期徑流與濃度之間的相關(guān)系數(shù)絕對(duì)值，表明徑流量與污染物通量之間的相關(guān)性更高。豐、枯水期高錳酸鹽指數(shù)濃度與徑流量都呈現(xiàn)負(fù)相關(guān)，說明徑流量增大對(duì)高錳酸鹽指數(shù)的稀釋作用起主導(dǎo)作用；豐、枯水期氨氮、總磷濃度與徑流量都呈正相關(guān)，說明徑流量增大引起的氨氮和總磷污染物入河增加量，大于徑流量對(duì)河流氮磷污染物的稀釋量(能力)。

(4) 河流水質(zhì)變化會(huì)受點(diǎn)源、面源污染兩方面的影響，豐、枯水期的高錳酸鹽指數(shù)濃度隨著徑流的增加而減小，說明高錳酸鹽指數(shù)污染以點(diǎn)源為主。豐、枯水期氨氮、總磷濃度隨著徑流的增加而增加，說明氮磷污染以面源為主。

(5) 注重豐水期污染防控，采用農(nóng)藝措施減少雨季肥料流失，采用河岸植被過濾帶、人工濕地等生態(tài)措施消減氮磷污染的入河量。

致 謝

本研究的野外采樣和室內(nèi)分析得到了“廣西環(huán)境污染控制理論與技術(shù)重點(diǎn)實(shí)驗(yàn)室科教結(jié)合科技創(chuàng)新基地”和“南方石山地區(qū)礦山地質(zhì)環(huán)境修復(fù)工程技術(shù)創(chuàng)新中心”的支持，在此一并致謝！