劉　琦，孫冰心

(1.雞西水文局，黑龍江 雞西 158100；2.牡丹江水文局，黑龍江 牡丹江 157000)

?

M-K和滑動(dòng)T確定穆棱河流域徑流變化突變時(shí)間

劉琦1，孫冰心2

(1.雞西水文局，黑龍江 雞西 158100；2.牡丹江水文局，黑龍江 牡丹江 157000)

摘要：Mann-Kendall和滑動(dòng)T檢驗(yàn)法是近年來被廣泛應(yīng)用的特征分析方法，為了提高研究精度和準(zhǔn)確度，將兩種方法結(jié)合起來，以M法確定突變節(jié)點(diǎn)，用滑動(dòng)T檢驗(yàn)M法結(jié)論的準(zhǔn)確性。穆棱河作為黑龍江省的十大河流之一，近年來流域附近對(duì)水資源的開采程度日益加大，因此，研究該流域徑流變化具有十分重要的意義。

關(guān)鍵詞：徑流；Mann-Kendall；滑動(dòng)T；穆棱河流域

穆棱河總長834 km，流域面積18 427 km2，位于黑龍江省東部。河寬25～140 m，水深0.8～2.2 m。多年平均結(jié)冰期為11月中旬至次年4月上旬，平均為144 d，平均最大冰厚0.9 m[1]，河道總落差699 m，多年平均徑流量為23.5億m3，折合徑流深為133.5 mm。該流域?qū)嶋H供水量23.821億m3，其中地表水供水12.290億m3，占實(shí)際供水量的51.6%，地下水供水11.531億m3，占48.4%[1]。

由于穆棱河在黑龍江省水域系統(tǒng)占據(jù)的重要地位，研究地表徑流量變化規(guī)律，預(yù)測未來流量趨勢變化，對(duì)該區(qū)農(nóng)業(yè)節(jié)水灌溉及水資源可持續(xù)利用與管理具有十分重要意義。論文以穆棱河流域上游至下游水文站徑流實(shí)測資料為基礎(chǔ)，依據(jù)Mann-Kendall秩次檢驗(yàn)法和滑動(dòng)T檢驗(yàn)法兩種方法結(jié)合應(yīng)用，分析徑流變化發(fā)生突變的時(shí)間點(diǎn)，確定突變的發(fā)生時(shí)間，以便更好的掌握水情變化，為人民的生命、財(cái)產(chǎn)安全和防汛抗旱工作提供技術(shù)支持。

1理論方法概述

1.1Mann-Kendall秩次檢驗(yàn)法

Mann-Kendall秩次檢驗(yàn)法是由Mann在1945年提出的，Mann-Kendall檢驗(yàn)方法屬于非參數(shù)統(tǒng)計(jì)檢驗(yàn)的方法，最初由Mann和Kendall[2]提出，由于該檢驗(yàn)方法不要求樣本遵從一定的分布，且檢驗(yàn)結(jié)果基本不受少數(shù)異常值的干擾，所以被廣泛應(yīng)用于分析降水、氣溫、徑流等要素隨時(shí)間序列的趨勢變化。其基本原理是在假設(shè)原假設(shè)Ho在時(shí)間序列數(shù)據(jù)為同一系列數(shù)據(jù)，沒有變化的情況下，設(shè)此連續(xù)的時(shí)間序列數(shù)據(jù)為X1，X2，X3，……，Xi表示第i個(gè)樣本的累計(jì)數(shù)，其中樣本中的Xi>Xj(1

它的檢驗(yàn)原理是假設(shè)有一時(shí)間序列如下：X1，X2，X3，……，Xn，構(gòu)造一秩序列mi，mi表示Xi>Xj(1≤j≤i)的樣本累積數(shù)。構(gòu)造統(tǒng)計(jì)量[5]：

(1)

式中：mi為某一元素xi大于前面xi-1，xi-2，…，x2，x1的累積數(shù)。

在原序列隨機(jī)獨(dú)立的情況下，dk的均值、方差分別為：

(2)

(3)

然后將dk標(biāo)準(zhǔn)化，構(gòu)造統(tǒng)計(jì)量：

(4)

其中，UF1=0，UFi服從標(biāo)準(zhǔn)正態(tài)分布。

取顯著水平α，如果|UFi|>Uα，說明序列趨勢性顯著。再計(jì)算反序列，重復(fù)式(1)至(4)，得到UBk，k=n，n-1，…，1。同樣，UB1=0。

計(jì)算步驟：①計(jì)算順序時(shí)間序列的秩序列dk，并按式(4)計(jì)算UFk。②計(jì)算逆序時(shí)間序列的秩序列dk，也按式(4)計(jì)算出UBk。③給定顯著性水平=0.05，那么臨界值U0.05=±1.96。將UFk和UBk兩條統(tǒng)計(jì)量序列曲線和臨界值±1.96繪于同一圖上。

分析統(tǒng)計(jì)量序列曲線。當(dāng)UFk或UBk大于0時(shí)，時(shí)間序列呈上升趨勢，當(dāng)UFk或UBk小于0時(shí)，時(shí)間序列呈下降趨勢。如果曲線超過臨界值線，則序列趨勢性顯著。當(dāng)統(tǒng)計(jì)量序列曲線出現(xiàn)交點(diǎn)，且交點(diǎn)在臨界線之間時(shí)，此時(shí)交點(diǎn)所在的年份即為可能發(fā)生突變的時(shí)間。

1.2滑動(dòng)T檢驗(yàn)法

滑動(dòng)T檢驗(yàn)法是利用T檢驗(yàn)法的原理，對(duì)序列逐點(diǎn)進(jìn)行T檢驗(yàn)。T檢驗(yàn)法的原理[6]是：設(shè)滑動(dòng)點(diǎn)前后，兩個(gè)序列總體的分布函數(shù)各為F1(x)和F2(x)，從總體F1(x)和F2(x)中分別抽取容量為n1和n2的兩個(gè)樣本，要求檢驗(yàn)原假設(shè)：F1(x) =F2(x)，則

(5)

其中，

(6)

(7)

(8)

T服從t(n1+n2-2)分布，選擇顯著性水平α，查t分布表得到臨界值tα/2，當(dāng)T>tα/2時(shí)，拒絕原假設(shè)，說明其存在顯著性差異；當(dāng)T

2徑流特征趨勢變化實(shí)例分析

2.1資料選取

水文資料是水文分析和計(jì)算的依據(jù)，它直接影響著整個(gè)工程的設(shè)計(jì)精度和安全。所以對(duì)水文資料的選擇必須嚴(yán)格依照可靠性、一致性、代表性的審查標(biāo)準(zhǔn)。

水文數(shù)據(jù)主要包括穆棱河上游—下游水文站的逐日徑流量數(shù)據(jù)。水文資料主要源于水文年鑒(松嫩流域卷)。水文站的位置及水文資料年限如表1所示。

特別要說明的是：密山橋水文站1979年建站，因此，水量資料是從1979年開始的。

表1　四個(gè)水文站的基本情況

2.2實(shí)例分析

采用Mann-Kendall檢驗(yàn)法對(duì)穆棱河流域四個(gè)水文站的天然多年徑流序列進(jìn)行趨勢性分析和檢驗(yàn)。選擇顯著性水平，對(duì)穆棱河流域徑流序列進(jìn)行滑動(dòng)T檢驗(yàn)法分析計(jì)算。計(jì)算結(jié)果如圖1、圖2所示。

圖1　穆棱河流域年徑流序列突變點(diǎn)M-K分析圖

圖2　滑動(dòng)T檢驗(yàn)法分析圖

由圖1可知，1975—1982年徑流量的統(tǒng)計(jì)量UF值和UB值超過了置信水平0.05相應(yīng)的臨界值(Y=1.96和Y=-1.96)，表明徑流量變化趨勢較顯著。UF和UB曲線相交于1963年附近，且交點(diǎn)位于兩臨界線之間，表明徑流量在1963年附近發(fā)生了一次減少的突變。UF和UB曲線在1998年再一次相交，徑流又經(jīng)歷了一次突變。

圖1中可見，在臨近點(diǎn)±1.96范圍內(nèi)出現(xiàn)的UF和UB交匯點(diǎn)分別在1954年、1959年、1964年和1998年附近。根據(jù)M-K檢驗(yàn)法規(guī)定，統(tǒng)計(jì)量序列曲線出現(xiàn)交點(diǎn)，且交點(diǎn)在臨界線之間時(shí)，此時(shí)交點(diǎn)所在的年份即為可能發(fā)生突變的時(shí)間。然后，根據(jù)圖1得到的突變時(shí)間點(diǎn)經(jīng)滑動(dòng)T驗(yàn)證，見圖2，圖2中可以看出，發(fā)生超出臨界水平的時(shí)間點(diǎn)出現(xiàn)在1998年附近，結(jié)合圖1，突變時(shí)間點(diǎn)是1998年，突變區(qū)間在1998年附近。

滑動(dòng)T是用來檢驗(yàn)M-K方法的檢驗(yàn)法，運(yùn)用滑動(dòng)T驗(yàn)證M-K方法適用于水量變化比較大的年份，1998年是穆棱河水量突增年份，百年一遇洪水年份，在兩種方法結(jié)合下得到確定性分析。通過M-K分析得到的其它突變年份，盡管在滑動(dòng)T檢驗(yàn)中沒有得到印證，但是查閱相關(guān)年鑒資料，1954年、1959年是小水年，1964年是中水年，均有不同明顯程度的水量增加趨勢。

徑流趨勢檢驗(yàn)方法Mann-Kendall和滑動(dòng)T兩種方法的搭配使用，既說明了出現(xiàn)突變時(shí)間，又印證了突變時(shí)間的準(zhǔn)確性。M-K檢驗(yàn)法觀察到1998年附近為突變時(shí)間區(qū)域，滑動(dòng)T檢驗(yàn)法得到1998年附近確實(shí)存在突變時(shí)間點(diǎn)。由M-K法獲得的其它小水或中水年份，通過年鑒資料得到了充分的驗(yàn)證，M-K劃分出的突變時(shí)間點(diǎn)均屬有效突變時(shí)間點(diǎn)，在未來的水量預(yù)測、防洪抗汛、水資源優(yōu)化利用等方面，M-K可以做為定性分析方法應(yīng)用于穆棱河中長期徑流特征分析中。

3結(jié)語

徑流趨勢檢驗(yàn)方法Mann-Kendall和滑動(dòng)T兩種方法的搭配使用，既說明了出現(xiàn)突變時(shí)間，又印證了突變時(shí)間的準(zhǔn)確性。M-K檢驗(yàn)法和滑動(dòng)T共同分析得到1998年附近為突變時(shí)間區(qū)域，其它突變時(shí)間點(diǎn)均在年鑒中得到論證。通過穆棱河流域上測站的水量資料和檢驗(yàn)方法的搭配使用，M-K應(yīng)用于實(shí)際生產(chǎn)生活中具有科學(xué)性、可操作性和實(shí)際意義。

參考文獻(xiàn)：

[1]王鳳梅.穆棱河文化的起源和發(fā)展[J].戲劇之家(上半月)，2011(12):81.

[2]郭巧玲，楊云松，暢祥生，等.1957—2008年黑河流域徑流年內(nèi)分配變化[J].地理科學(xué)進(jìn)展，2011，30(5):550-557.

[3]尹云鶴，吳紹洪，陳剛.1961—2006年我國氣候變化趨勢與突變的區(qū)域差異[J].自然資源學(xué)報(bào)，2009，24(12):2147-2157.

[4]陳學(xué)君，姜夢晶.甘肅定西地區(qū)干旱氣候突變分析[J].干旱區(qū)資源與環(huán)境，2012，26(11):68-75.

[5]王玉新．嫩江流域徑流演變規(guī)律及其歸因分析[D]．長春：吉林大學(xué)，2012.

[6]呂琳莉，李朝霞.雅魯藏布江中下游徑流變異性識(shí)別[J].水力發(fā)電，2013，39 (5)：13-15.

Ensure the abrupt runoff change time in Muling river byM-K and sliding T Method

LIU Qi1，SUN Bingxin2

(1.HydrologicalBureauofJixi，Jixi158100，China；2.HydrologicalBureauofMudanjiang，Mudanjiang158100，China)

Abstract：The characteristic analysis method of M-K and sliding T are both widely used in recent years. These two methods are combined in order to improve the research accuracy. The method of M-K is used to make sure the abrupt runoff change time point, and then the results accuracy of which is validated by method of slide T. Muling River is one of the most important river in Heilongjiang province, the water resource exploitation degree of which increases daily in recent years. So, it is extremely significant to research this basin water resource.

Key words：runoff；M-K method；sliding T；Muling river basin

作者簡介：劉琦(1989-)，女，助理工程師，長期從事水文情報(bào)預(yù)報(bào)工作。E-mail：627291058@qq.com

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

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

文章編號(hào)：2096-0506(2015)02-0035-03