基于Mann-Kendall 檢驗(yàn)和重標(biāo)極差分析的水質(zhì)變化趨勢(shì)

董　旭,梅　琨,商　栩,黃樹輝,黃　宏①

(1.溫州市水利局珊溪管理局,浙江 溫州　325000;2.溫州醫(yī)科大學(xué)浙南水科學(xué)研究院,浙江 溫州　325035;3.浙江省流域水環(huán)境與健康風(fēng)險(xiǎn)研究重點(diǎn)實(shí)驗(yàn)室,浙江 溫州　325035;4.溫州醫(yī)科大學(xué)環(huán)境與公共衛(wèi)生學(xué)院環(huán)境科學(xué)系,浙江 溫州　325035)

?

基于Mann-Kendall 檢驗(yàn)和重標(biāo)極差分析的水質(zhì)變化趨勢(shì)

董旭1,梅琨2,商栩3,黃樹輝4,黃宏2①

(1.溫州市水利局珊溪管理局,浙江 溫州325000;2.溫州醫(yī)科大學(xué)浙南水科學(xué)研究院,浙江 溫州325035;3.浙江省流域水環(huán)境與健康風(fēng)險(xiǎn)研究重點(diǎn)實(shí)驗(yàn)室,浙江 溫州325035;4.溫州醫(yī)科大學(xué)環(huán)境與公共衛(wèi)生學(xué)院環(huán)境科學(xué)系,浙江 溫州325035)

摘要：將Mann-Kendall檢驗(yàn)和重標(biāo)極差(R/S)分析相結(jié)合，建立了一套綜合分析過(guò)去和未來(lái)水質(zhì)變化趨勢(shì)的方法。首先，用季節(jié)性Mann-Kendall檢驗(yàn)識(shí)別水質(zhì)時(shí)間序列的過(guò)去變化趨勢(shì)；然后，用R/S分析方法計(jì)算水質(zhì)時(shí)間序列的Hurst指數(shù)(IH)，0

關(guān)鍵詞：水質(zhì);變化趨勢(shì);Mann-Kendall檢驗(yàn);重標(biāo)極差分析;長(zhǎng)程相關(guān)性

由于水環(huán)境保護(hù)長(zhǎng)期滯后于經(jīng)濟(jì)發(fā)展,水污染特別是飲用水污染已成為制約許多國(guó)家和地區(qū)可持續(xù)發(fā)展的突出問(wèn)題,水環(huán)境特別是水源地環(huán)境質(zhì)量的改善已成為事關(guān)國(guó)計(jì)民生的重大事項(xiàng)。由于水環(huán)境系統(tǒng)的復(fù)雜性,水質(zhì)數(shù)據(jù)資料往往具有非正態(tài)性、季節(jié)變異性、出現(xiàn)漏測(cè)值和小于檢測(cè)限值等特征,使得常規(guī)參數(shù)統(tǒng)計(jì)分析方法的應(yīng)用受到限制[1]。因此,目前被廣泛應(yīng)用于水質(zhì)變化趨勢(shì)分析的方法主要是Mann-Kendall檢驗(yàn)[2]、Sen′s slope estimates[3]等非參數(shù)統(tǒng)計(jì)分析方法。然而,這些方法側(cè)重于定性或半定量地分析時(shí)間序列資料在過(guò)去的變化趨勢(shì),而不能直接推斷水環(huán)境質(zhì)量的未來(lái)變化趨勢(shì)。20世紀(jì)50年代以來(lái),HURST[4]發(fā)現(xiàn)水文時(shí)間序列具有長(zhǎng)程相關(guān)性(持續(xù)性)或長(zhǎng)程反相關(guān)性(反持續(xù)性),即時(shí)間序列的未來(lái)變化趨勢(shì)跟過(guò)去變化趨勢(shì)一致或相反。近年來(lái)研究發(fā)現(xiàn),河流水質(zhì)也具有長(zhǎng)程相關(guān)性或長(zhǎng)程反相關(guān)性[5]。識(shí)別水質(zhì)時(shí)間序列的長(zhǎng)程相關(guān)性或長(zhǎng)程反相關(guān)性有助于了解水質(zhì)動(dòng)態(tài)變化規(guī)律,也為在趨勢(shì)性分析的基礎(chǔ)上進(jìn)一步推斷未來(lái)水質(zhì)變化提供依據(jù)。

該研究的目的是建立一套能同時(shí)分析過(guò)去和未來(lái)水質(zhì)變化趨勢(shì)的方法。首先,基于非參數(shù)檢驗(yàn)方法分析在過(guò)去一段時(shí)期內(nèi)水質(zhì)的變化趨勢(shì);然后,基于長(zhǎng)程相關(guān)性分析方法判斷水質(zhì)的未來(lái)變化趨勢(shì)與過(guò)去是否一致;最后,結(jié)合流域污染控制規(guī)劃和進(jìn)展綜合推斷水質(zhì)的未來(lái)變化趨勢(shì)。此外,還以浙江省溫州市珊溪水庫(kù)及其入庫(kù)支流總氮(TN)濃度和高錳酸鹽指數(shù)(CODMn)的變化趨勢(shì)為例開展實(shí)例研究。

1研究方法

1.1季節(jié)性Mann-Kendall檢驗(yàn)

季節(jié)性Mann-Kendall檢驗(yàn)是Mann-Kendall檢驗(yàn)[6]的一種推廣,其零假設(shè)H0為隨機(jī)變量與時(shí)間獨(dú)立,且全年 12 個(gè)月的水質(zhì)資料具有相同的概率分布[7-8]。假設(shè)對(duì)某一個(gè)斷面進(jìn)行n年逐月監(jiān)測(cè),對(duì)該水質(zhì)時(shí)間序列數(shù)據(jù)進(jìn)行季節(jié)性Mann-Kendall檢驗(yàn)的計(jì)算過(guò)程可簡(jiǎn)要概括為3個(gè)主要步驟。

(1)計(jì)算第i月的差值統(tǒng)計(jì)量Si和方差VSi。對(duì)第i月的歷年水質(zhì)時(shí)間序列進(jìn)行比較,如果后面的值大于前面的值記為正號(hào),否則記為負(fù)號(hào)。令正負(fù)號(hào)之和為Si:

(1)

式(1)中,n為年數(shù)；xij和xik為水質(zhì)數(shù)據(jù)；sgn為函數(shù)符號(hào)。

(2)

在零假設(shè)下,均值ESi服從正態(tài)分布:

ESi=0，

(3)

VSi=ni(ni-1)(2ni+5)/18。

(4)

(2)計(jì)算全年的差值統(tǒng)計(jì)量S和方差VS。全年的S和VS分別為各月之和。如果季節(jié)數(shù)和年份數(shù)足夠大,如3a以上的逐月水質(zhì)監(jiān)測(cè)數(shù)據(jù)[6-7],S服從正態(tài)分布,標(biāo)準(zhǔn)方差Z為

(5)

(3)顯著性趨勢(shì)檢驗(yàn)。Kendall 檢驗(yàn)統(tǒng)計(jì)量τ定義為:τ=S/m(m為可作比較的差值數(shù)據(jù)組個(gè)數(shù)之和)。τ>0或τ<0分別指示水質(zhì)時(shí)間序列具有上升或下降趨勢(shì),τ=0則指示無(wú)趨勢(shì)。在雙邊的趨勢(shì)檢驗(yàn)中,在給定的α置信水平上,如果 |Z|≥Z1-α/2,則拒絕原假設(shè)H0。通常取顯著性水平α為 0.10 和 0.01。α≤0.01或0.01<α≤0.10分別指示檢驗(yàn)具有極顯著或顯著水平。

1.2重標(biāo)極差法(R/S)分析

現(xiàn)有的長(zhǎng)程相關(guān)性分析方法有重標(biāo)極差法(R/S)[9]、去趨勢(shì)波動(dòng)分析法(detrended fluctuation analysis,DFA)[5]等。R/S分析方法的特點(diǎn)是能揭示時(shí)間序的分形特征,能明確時(shí)間序列未來(lái)的變化特征與過(guò)去是相同的還是相反的[9]。假設(shè)對(duì)某一個(gè)斷面進(jìn)行多年逐月水質(zhì)監(jiān)測(cè),對(duì)該數(shù)據(jù)序列進(jìn)行R/S分析的主要原理和步驟為:將水質(zhì)數(shù)據(jù)按監(jiān)測(cè)日期先后順序排列,構(gòu)成一個(gè)水質(zhì)時(shí)間序列數(shù)據(jù){ξt=x1,x2，…,xn}。對(duì)于任意正整數(shù)k≥1,其均值系列為

(6)

累積離差(Ftk)為

(7)

極差(Rk)為

Rk=maxFtk-minFtk,1≤t≤k。

(8)

標(biāo)準(zhǔn)差(Sk)為

(9)

如果存在以下關(guān)系:

Rk/Sk∝kIH。

(10)

則說(shuō)明時(shí)間序列存在Hurst現(xiàn)象,式(10)中IH為Hurst指數(shù)。IH=0.5指示時(shí)間序列完全獨(dú)立,0

1.3水質(zhì)變化趨勢(shì)綜合判斷

Mann-Kendall檢驗(yàn)和R/S分析分別都被廣泛應(yīng)用于水質(zhì)的趨勢(shì)分析和長(zhǎng)程相關(guān)性分析,然而它們各有側(cè)重點(diǎn),單獨(dú)依靠其中任何一種方法都不足以判斷未來(lái)的水質(zhì)變化趨勢(shì)。整體而言,Mann-Kendall檢驗(yàn)側(cè)重于分析過(guò)去一段時(shí)期內(nèi)時(shí)間序列的趨勢(shì)性,R/S分析方法側(cè)重于分析時(shí)間序列的長(zhǎng)程相關(guān)(持續(xù)性)和長(zhǎng)程反相關(guān)性(反持續(xù)性)。將Mann-Kendall檢驗(yàn)和R/S分析相結(jié)合,不僅能鑒別在過(guò)去一段時(shí)間內(nèi)水質(zhì)的變化趨勢(shì)上升抑或下降,還能明確未來(lái)水質(zhì)變化趨勢(shì)跟過(guò)去變化趨勢(shì)相同抑或相反,進(jìn)而可以結(jié)合流域污染控制規(guī)劃和進(jìn)展推斷未來(lái)水質(zhì)的變化趨勢(shì)。

2研究區(qū)域與數(shù)據(jù)收集

珊溪水利樞紐(27.46°～27.99° N,119.62°～120.27° E)位于飛云江上游,由珊溪水庫(kù)、趙山渡水庫(kù)和配套引水渠系3部分組成,流域集水面積約2 300 km2,被稱為溫州人民的“大水缸”(圖1)。

圖1　水質(zhì)監(jiān)測(cè)點(diǎn)位置

流域集水區(qū)內(nèi)人口總數(shù)約 54 萬(wàn)人,土地利用類型以林地為主,約占總面積的 71%;其次為耕地,約占 20%。珊溪水庫(kù)壩址以上控制流域總面積約1 529 km2。多年平均入庫(kù)徑流量18.6億m3,正常蓄水位 142 m,相應(yīng)庫(kù)容和水域面積分別為12.91億m3和36 km2。流域集水區(qū)內(nèi)工業(yè)已完成搬遷和轉(zhuǎn)型,但受農(nóng)村生活污水、畜禽養(yǎng)殖廢棄物和農(nóng)業(yè)化肥等的影響,局部庫(kù)灣曾發(fā)生藻類異常增殖現(xiàn)象。隨著水源地環(huán)境綜合整治的不斷推進(jìn),入庫(kù)支流水質(zhì)惡化的趨勢(shì)得到有效遏制,庫(kù)區(qū)水體整體優(yōu)良,但個(gè)別月份TN和CODMn等參數(shù)屬于GB 3838—2002《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》中Ⅲ類水質(zhì)。

水質(zhì)時(shí)間序列的長(zhǎng)度對(duì)趨勢(shì)檢驗(yàn)具有很大影響,過(guò)短的時(shí)間(如2～3 a)序列不能準(zhǔn)確判斷是否存在趨勢(shì),過(guò)長(zhǎng)的時(shí)間序列則可能導(dǎo)致一種趨勢(shì)掩蓋或抵消另一種趨勢(shì)[7]。一般認(rèn)為,用季節(jié)性Mann-Kendall檢驗(yàn)判斷水質(zhì)趨勢(shì)時(shí),序列長(zhǎng)度一般以5～8 a為宜。研究數(shù)據(jù)從當(dāng)?shù)厮块T獲取,為2010—2014年珊溪水庫(kù)及主要支流入庫(kù)斷面的逐月監(jiān)測(cè)數(shù)據(jù)(圖2)。根據(jù)庫(kù)區(qū)水質(zhì)現(xiàn)狀,選取TN和CODMn參數(shù)作為研究對(duì)象。TN濃度的測(cè)定采用堿性過(guò)硫酸鉀消解-紫外分光光度法(GB 11894—1989),CODMn的測(cè)定采用滴定法(GB/T 11892—1989)。

3結(jié)果與分析

3.1季節(jié)性Mann-Kendall檢驗(yàn)結(jié)果

2010—2014年,珊溪水庫(kù)庫(kù)區(qū)和主要入庫(kù)支流逐月TN濃度和CODMn時(shí)間序列見圖2～3,季節(jié)性Mann-Kendall檢驗(yàn)結(jié)果見表1。司前溪TN濃度明顯上升,黃坦坑和莒江溪明顯下降,峃作口溪、三插溪和洪口溪均下降但不明顯(表1)。司前溪CODMn明顯上升,三插溪和洪口溪均上升但不明顯,峃作口溪、黃坦坑和莒江溪均明顯下降 (表1)。跟其他支流相比,黃坦坑TN濃度和CODMn下降趨勢(shì)最明顯(圖2～3)。黃坦坑流域內(nèi)的黃坦鎮(zhèn)曾經(jīng)是文成縣最大的生豬養(yǎng)殖基地,2011年生豬50頭以上的養(yǎng)殖場(chǎng)達(dá)302個(gè)。大量畜禽糞便未經(jīng)治理直排入河,導(dǎo)致水質(zhì)嚴(yán)重惡化,長(zhǎng)期屬于劣V類,在枯水期水體甚至發(fā)黑發(fā)臭。為保護(hù)珊溪水庫(kù)水質(zhì),管理部門在黃坦坑小流域開展水土保持生態(tài)清潔型小流域建設(shè)工程,申報(bào)并實(shí)施農(nóng)村沼氣國(guó)債項(xiàng)目和浙江省“811”環(huán)境整治項(xiàng)目,至2012年基本完成[10],并取得了明顯成效。珊溪水庫(kù)庫(kù)區(qū)水體TN濃度和CODMn分別為上升和下降(表1)。2012年以來(lái),溫州市大力推進(jìn)庫(kù)區(qū)農(nóng)房改造和跨區(qū)域統(tǒng)籌集聚,最大限度減少一、二級(jí)水源保護(hù)區(qū)人口數(shù)量,同時(shí)全面實(shí)施水源保護(hù)“五大”工程建設(shè),包括生活污水治理、生活垃圾治理、畜禽養(yǎng)殖污染治理、主要支流生態(tài)保護(hù)與修復(fù)、水質(zhì)自動(dòng)在線監(jiān)測(cè)和預(yù)警應(yīng)急體系等水源保護(hù)“五大”工程建設(shè),最大限度減少入庫(kù)污染物總量。為進(jìn)一步評(píng)估水源地綜合整治效果,對(duì)各監(jiān)測(cè)斷面2012年前后TN濃度和CODMn進(jìn)行Mann-Whitney檢驗(yàn) (表2)。

圖2　2010—2014年珊溪水庫(kù)和支流逐月TN濃度

圖3　2010—2014年珊溪水庫(kù)和支流逐月CODMn

結(jié)果表明,司前溪2012年以后TN濃度顯著大于2012年以前,黃坦坑、洪口溪和莒江溪?jiǎng)t相反,峃作口溪和三插溪TN濃度在2012年前后差異不顯著。三插溪、司前溪和洪口溪CODMn在2012年前后差異不顯著,峃作口溪、黃坦坑和莒江溪2012年以后則顯著小于2012年以前。庫(kù)區(qū)TN濃度在2012年前后幾乎沒(méi)有變化,而CODMn在2012年以后有所下降,但差異不顯著?？梢?庫(kù)區(qū)和大部分支流水質(zhì)惡化的勢(shì)頭得到了有效遏制,說(shuō)明流域環(huán)境污染綜合整治工作富有成效。然而,司前溪TN濃度和CODMn都具有顯著上升趨勢(shì),原因可能是司前溪監(jiān)測(cè)斷面位于司前鎮(zhèn),隨著近年來(lái)城鎮(zhèn)規(guī)模的擴(kuò)大和人口的增加,大量生活污水未經(jīng)處理直接排放入河。季節(jié)性Mann-Kendall檢驗(yàn)明確了2010—2014年期間水質(zhì)變化趨勢(shì),為結(jié)合R/S分析和流域污染控制規(guī)劃判斷水質(zhì)未來(lái)變化趨勢(shì)奠定了基礎(chǔ)[9]。

表12010—2014年珊溪水庫(kù)和支流TN濃度和CODMn季節(jié)性Mann-Kendall檢驗(yàn)和R/S分析結(jié)果

Table 1Seasonal Mann-Kendall test and R/S analysis of monthly TN and CODMnconcentrations of the Shanxi Reservoir and its tributaries during 2010-2014

水質(zhì)參數(shù)支流和水庫(kù)Mann-Kendall檢驗(yàn)結(jié)果R/S分析結(jié)果τα過(guò)去變化趨勢(shì)IH長(zhǎng)程(反)相關(guān)性未來(lái)水質(zhì)變化趨勢(shì)TN濃度峃作口溪-0.160.18下降但不顯著0.68較強(qiáng)持續(xù)性不會(huì)惡化黃坦坑-0.72<0.01很顯著下降0.89很強(qiáng)持續(xù)性改善三插溪-0.030.83下降但不顯著0.58較弱持續(xù)性不會(huì)惡化司前溪0.42<0.01很顯著上升0.61較弱持續(xù)性惡化洪口溪-0.170.16下降但不顯著0.63較弱持續(xù)性不會(huì)惡化莒江溪-0.310.01很顯著下降0.68較強(qiáng)持續(xù)性改善珊溪水庫(kù)0.030.83上升但不顯著0.72較強(qiáng)持續(xù)性不會(huì)惡化CODMn峃作口溪-0.39<0.01很顯著下降0.59較弱持續(xù)性改善黃坦坑-0.87<0.01很顯著下降0.96很強(qiáng)持續(xù)性改善三插溪0.110.35上升但不顯著0.71較強(qiáng)持續(xù)性不會(huì)惡化司前溪0.190.10很顯著上升0.77強(qiáng)持續(xù)性惡化洪口溪0.060.62上升但不顯著0.77強(qiáng)持續(xù)性不會(huì)惡化莒江溪-0.220.06顯著下降0.56較弱持續(xù)性改善珊溪水庫(kù)-0.250.04顯著下降0.72較強(qiáng)持續(xù)性改善

τ為Kendall 檢驗(yàn)統(tǒng)計(jì)量,τ>0、τ<0和τ=0分別指示水質(zhì)時(shí)間序列具有上升趨勢(shì)、下降趨和無(wú)趨勢(shì)。α≤0.01或0.01<α≤0.10分別指示檢驗(yàn)具有極顯著或顯著水平。IH=0.5指示時(shí)間序列完全獨(dú)立,0

表2珊溪水庫(kù)和支流2012年前后TN濃度和CODMn描述性統(tǒng)計(jì)和Mann-Whitney檢驗(yàn)結(jié)果

Table 2Descriptive statistics and Mann-Whitney test of monthly TN concentrations and CODMnof the Shanxi Reservoir and its tributaries before and after 2012

支流和水庫(kù)TN濃度CODMn2010—2012年2013—2014年2010—2012年2013—2014年峃作口溪0.63±0.40a0.52±0.23a2.27±0.86a1.96±1.33b黃坦坑8.75±5.41a2.69±2.06b6.39±2.06a3.03±0.70b三插溪0.54±0.52a0.44±0.13a1.34±0.57a1.35±0.38a司前溪0.52±0.24b1.48±0.37a1.39±0.65a1.40±0.35a洪口溪0.54±0.50a0.30±0.11b1.51±0.52a1.46±0.47a莒江溪0.57±0.29a0.30±0.09b1.47±0.67a1.08±0.39b珊溪水庫(kù)0.42±0.12a0.42±0.15a1.82±0.66a1.57±0.44a

同一行數(shù)據(jù)后英文小寫字母不同表示相同取樣點(diǎn)不同年份某指標(biāo)差異顯著(P<0.05)。

3.2R/S分析結(jié)果

基于2010—2014年逐月水質(zhì)監(jiān)測(cè)數(shù)據(jù),珊溪水庫(kù)庫(kù)區(qū)和主要入庫(kù)支流逐月TN濃度和CODMn時(shí)間序列的R/S分析結(jié)果見表1。庫(kù)區(qū)和入庫(kù)支流水質(zhì)時(shí)間序列的IH都大于0.50,指示庫(kù)區(qū)和入庫(kù)支流水質(zhì)都具有長(zhǎng)程相關(guān)性(持續(xù)性),而沒(méi)有長(zhǎng)程反相關(guān)性(反持續(xù)性),即未來(lái)水質(zhì)變化趨勢(shì)跟2010—2014年期間的變化趨勢(shì)一致,區(qū)別僅在于持續(xù)性的強(qiáng)弱。參照以往研究,0.50

3.3水質(zhì)未來(lái)變化趨勢(shì)綜合推斷

水庫(kù)和入庫(kù)支流是連續(xù)的水體,水庫(kù)水質(zhì)跟入庫(kù)支流水質(zhì)密切相關(guān),遏制水庫(kù)水質(zhì)惡化、改善水庫(kù)水質(zhì),務(wù)必要從流域尺度大力推進(jìn)環(huán)境污染綜合整治工作。根據(jù)珊溪水源保護(hù)文件,溫州市將通過(guò)水源地人口統(tǒng)籌集聚及水源保護(hù)“五大工程”建設(shè),力爭(zhēng)在3～5 a內(nèi)恢復(fù)庫(kù)區(qū)生態(tài)平衡。從統(tǒng)計(jì)學(xué)角度,Mann-Kendall檢驗(yàn)和R/S分析指示了珊溪水庫(kù)大部分支流水質(zhì)惡化的勢(shì)頭將得到有效遏制,并且朝著改善的方向發(fā)展。隨著“五水共治”工作不斷推進(jìn),流域集水區(qū)內(nèi)生活污水治理、生活垃圾治理、畜禽養(yǎng)殖污染治理以及主要支流生態(tài)保護(hù)和修復(fù)等工程的不斷推進(jìn),大部分支流的水質(zhì)將呈現(xiàn)穩(wěn)定改善的趨勢(shì)。司前溪TN濃度和CODMn具有顯著上升趨勢(shì),而且具有長(zhǎng)程相關(guān)性,未來(lái)水質(zhì)很可能會(huì)繼續(xù)惡化。要盡快查明司前溪污染來(lái)源,采取針對(duì)性措施,從源頭上治理污染源,早日遏制住水質(zhì)惡化勢(shì)頭。

珊溪水庫(kù)是一座大型的多年調(diào)節(jié)深水水庫(kù),2012年以來(lái),為控制水庫(kù)庫(kù)區(qū)藍(lán)藻增殖,采取了濾食性魚類投放和3 a封庫(kù)管理措施。庫(kù)區(qū)CODMn時(shí)間序列呈顯著下降趨勢(shì),而且具有強(qiáng)持續(xù)性,未來(lái)將會(huì)持續(xù)改善。然而,由于水庫(kù)對(duì)養(yǎng)分的滯留效應(yīng)較強(qiáng)[12],而且除了入庫(kù)支流輸送的污染[13],珊溪水庫(kù)TN污染來(lái)源還包括大氣沉降、內(nèi)源污染和消落帶污染等。在這些因素的共同作用下,盡管庫(kù)區(qū)TN濃度已經(jīng)遏制住惡化勢(shì)頭,卻仍未呈現(xiàn)下降趨勢(shì)。可見,水環(huán)境綜合整治和生態(tài)系統(tǒng)恢復(fù)是一項(xiàng)長(zhǎng)期系統(tǒng)的工程,除了繼續(xù)大力推進(jìn)水源保護(hù)各項(xiàng)工程建設(shè)外,還要采取有效措施加強(qiáng)對(duì)消落帶污染和內(nèi)源污染的防控。

4結(jié)論

Mann-Kendall檢驗(yàn)側(cè)重于分析時(shí)間序列的過(guò)去變化趨勢(shì),R/S分析側(cè)重于分析時(shí)間序列的長(zhǎng)程相關(guān)性,將兩者相結(jié)合,不僅能鑒別水質(zhì)的過(guò)去變化趨勢(shì),還能鑒定水質(zhì)的未來(lái)變化趨勢(shì)與過(guò)去相同還是相反,進(jìn)而可以結(jié)合流域污染控制規(guī)劃和進(jìn)展推斷水質(zhì)的未來(lái)變化趨勢(shì)。

對(duì)珊溪水庫(kù)及其入庫(kù)支流水質(zhì)變化趨勢(shì)的實(shí)例研究表明,2010—2014年大部分支流TN濃度和CODMn具有下降趨勢(shì),未來(lái)水質(zhì)將會(huì)停止惡化或持續(xù)改善。然而水庫(kù)TN的滯留效應(yīng)較強(qiáng)且污染來(lái)源較多,TN濃度仍無(wú)明顯下降趨勢(shì),除了繼續(xù)控制入庫(kù)支流污染外,還要控制消落帶污染和內(nèi)源污染。該方法理論明確，操作簡(jiǎn)便，結(jié)果可靠,在水質(zhì)、水文、氣象等的趨勢(shì)分析中具有廣闊的應(yīng)用前景。

參考文獻(xiàn)：

[1]HUANG H,LU J.Identification of River Water Pollution Characteristics Based on Projection Pursuit and Factor Analysis [J].Environmental Earth Science,2014,72(9):3409-3417.

[2]KISI O,AY M.Comparison of Mann-Kendall and Innovative Trend Method for Water Quality Parameters of the Kizilirmak River,Turkey [J].Journal of Hydrology,2014,513(26):362-375.

[3]BOUZA-DEANO R,TERNERO-RODRIGUEZ M,FERNANDEZ-ESPINOSA A J.Trend Study and Assessment of Surface Water Quality in the Ebro River [J].Journal of Hydrology,2008,361 (3/4):227-239.

[4]HURST H.Long Term Storage Capacity of Reservoirs [J].Transactions of the American Society of Civil Engineers,1951,116:770-808.

[5]ONDERKA M,MRAFKOA L,KREINA,etal.Long-Term Persistence of Stream Nitrate Concentrations (Memory Effect) Inferred From Spectral Analysis and Detrended Fluctuation Analysis [J].Water Air & Soil Pollution,2012,223(1):241-252.

[6]丁文榮.橫斷山區(qū)干旱河谷氣候變化趨勢(shì)研究[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2013,29(6):681-687.

[7]孫松.季節(jié)性Kendall檢驗(yàn)法在黃壇水庫(kù)水質(zhì)趨勢(shì)分析中的應(yīng)用[J].廣州環(huán)境科學(xué),2011,26(2):580-585.

[8]代堂剛.運(yùn)用季節(jié)性Kendall檢驗(yàn)數(shù)學(xué)模型分析漁洞水庫(kù)水質(zhì)變化趨勢(shì)[J].水資源研究,2009,30(2):14-15.

[9]于延勝,陳興偉.R/S和Mann-Kendall法綜合分析水文時(shí)間序列未來(lái)的趨勢(shì)特征[J].水資源與水工程學(xué)報(bào),2008,19(3):41-44.

[10]王冠,張錦娟,莊景.水庫(kù)水源地面源污染防治對(duì)策及實(shí)踐[J].安徽農(nóng)業(yè)科學(xué),2013,41(10):4578-4580.

[11]張國(guó)存,查良松.南京近50年來(lái)氣候變化及未來(lái)趨勢(shì)分析[J].南京師范大學(xué)學(xué)報(bào):自然科學(xué)版,2008,31(6):580-586.

[12]胡春華,樓倩,丁文軍,等.鄱陽(yáng)湖氮、磷營(yíng)養(yǎng)鹽的滯留效應(yīng)研究[J].環(huán)境污染與防治,2012,34(9):1-4.

[13]陳玉東,陳梅,張龍江,等.基于GIS的龍墩水庫(kù)典型小流域面源污染氮磷負(fù)荷研究[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2015,31(4):492-499.

(責(zé)任編輯： 陳昕)

Analysis of Variation Trend of Water Quality Based on Mann-Kendall Test and Rescaled Range Analysis.

DONGXu1,MEIKun2,SHANGXu3,HUANGShu-hui4,HUANGHong2

(1.Wenzhou Water Resources Bureau， Shanxi Authority, Wenzhou 325000, China;2.Southern Zhejiang Water Research Institute, Wenzhou Medical University, Wenzhou 325035, China;3.Key Laboratory of Watershed Environmental Science and Health of Zhejiang Province, Wenzhou 325035, China;4.Department of Environmental Science, School of Environmental Science and Public Health, Wenzhou Medical University, Wenzhou 325035, China)

Abstract:A comprehensive set of methods was established based on Mann-Kendall test and rescaled range (R/S) analysis for analysis of variation trends of water quality in the past and future. Firstly, seasonal Mann-Kendall tests were performed to determine variation trend of water quality time series in the past. And then, R/S analysis was done used to calculate Hurst indexes (IH)of the water quality time series(0

Key words:water quality;variation trend;Mann-Kendall test;rescaled range analysis;long-term correlation

作者簡(jiǎn)介：董旭(1987—),男,浙江溫州人,助理工程師,學(xué)士,主要研究方向?yàn)樗Y源利用與保護(hù)。E-mail: 276205244@qq.com

DOI：10.11934/j.issn.1673-4831.2016.02.016

中圖分類號(hào)：X522

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1673-4831(2016)02-0277-06

通信作者①E-mail: huanghongpanda@163.com

基金項(xiàng)目：溫州市公益性科技計(jì)劃(S20140014);溫州市水體污染控制與治理科技創(chuàng)新項(xiàng)目(S20140039,S20140037);溫州醫(yī)科大學(xué)人才科研啟動(dòng)基金(QTJ14045);浙江省自然科學(xué)基金(LQ16C030004)

收稿日期：2015-11-11