谷桂華,李能,吳獻(xiàn)花
(1.云南省水文水資源局 玉溪分局,云南 玉溪 653100;2.玉溪師范學(xué)院 高原湖泊生態(tài)環(huán)境保護(hù)重點(diǎn)實(shí)驗(yàn)室,云南 玉溪 653100)
綠汁江流域地處滇中巖溶高原湖盆區(qū),位于云南省楚雄市祿豐縣、雙柏縣和玉溪市易門縣、峨山縣、新平縣境內(nèi).河源段稱為西河,出祿豐壩子后稱星宿江,流經(jīng)易門縣綠汁鎮(zhèn)后稱祿汁江.流域西部為哀牢山,干流全長(zhǎng)325 km,總面積8 573 km2,其中玉溪市境內(nèi)干流長(zhǎng)168.3 km,流域面積2 573.8 km2.流域內(nèi)復(fù)雜的地貌形態(tài)造成氣候類型多樣,其中以亞熱帶高原季風(fēng)氣候?yàn)橹?作為云南省熱區(qū)資源最為豐富的地區(qū)之一,紅河谷——綠汁江流域具有光熱資源豐富、生態(tài)污染小、民族文化多樣、土地可開發(fā)潛力大等優(yōu)勢(shì).通過對(duì)綠汁江流域1960~2016年降水和徑流變化特征進(jìn)行分析,成果可為綠汁江流域水資源開發(fā)利用、水資源科學(xué)管理調(diào)配和綠汁江熱區(qū)經(jīng)濟(jì)帶建設(shè)提供基礎(chǔ)參考.
選用綠汁江干流鴉勒水文站1960~2016年天然徑流數(shù)據(jù),分析流域內(nèi)徑流變化趨勢(shì);選用該流域玉溪境內(nèi)易門站和鴉勒雨量站1960~2016年逐月降水?dāng)?shù)據(jù),分析流域降水變化特征.其中,鴉勒水文站1976~1989年徑流數(shù)據(jù)由上游江邊站同期數(shù)據(jù)經(jīng)面積比擬法插補(bǔ)得到、1976~1979降水由周邊雨量站同期數(shù)據(jù)建立相關(guān)關(guān)系插補(bǔ)得到;流域平均降水量由代表站雨量用算術(shù)平均法計(jì)算得到.所選用和插補(bǔ)計(jì)算數(shù)據(jù)系列成果經(jīng)分析考證,均具有可靠性、代表性和一致性,可供分析使用.
回歸分析法和滑動(dòng)平均法回歸分析是水文統(tǒng)計(jì)學(xué)中常用的一種方法,主要包含線性回歸和非線性回歸.線性回歸又可分為一元線性回歸和多元線性回歸[1].一元線性回歸法研究某一變量隨時(shí)間的變化趨勢(shì)[2],建立y=bx+a,y為依賴x的變化而變化的某一變量;a、b為回歸方程的截距和斜率,根據(jù)斜率的正負(fù)可以得知某變量的變化趨勢(shì)(若b>0,則當(dāng)前系列為增加趨勢(shì);反之為減小趨勢(shì)).
滑動(dòng)平均法又稱移動(dòng)平均法[3],即沿著全長(zhǎng)N個(gè)數(shù)據(jù),逐個(gè)地取相鄰的2n+1個(gè)數(shù)據(jù)進(jìn)行算術(shù)平均,其值作為滑動(dòng)平均值.其具體的計(jì)算公式如下:
(1)
運(yùn)用此方法可以很好地消除系列數(shù)據(jù)的不穩(wěn)定波動(dòng),顯示出序列變化的平穩(wěn)性.
本文采用一元線性回歸分析法和5年滑動(dòng)平均(即n=5)對(duì)綠汁江流域年徑流和年、季降水序列進(jìn)行趨勢(shì)分析.
Mann-Kendall法Mann-Kendall法(簡(jiǎn)稱M-K法),是用來評(píng)估氣候要素時(shí)間序列趨勢(shì)的檢驗(yàn)方法,以適用范圍廣、人為性少、定量化程度高而著稱[4].
對(duì)時(shí)間序列x1,x2,…,xn(n為樣本數(shù)),所有對(duì)偶觀測(cè)值(xi,xj,j>i)中xi (2) 式中: τ為Kendall秩統(tǒng)計(jì)量; U為Kendall秩次相關(guān)系數(shù); P為系列中所有對(duì)偶觀測(cè)值(xi,xj,i ①對(duì)于無趨勢(shì)的序列E(P)=n(n-1)/4; ②當(dāng)P ③當(dāng)P>E(P)時(shí)表示序列可能有上升趨勢(shì); ④當(dāng)n增加時(shí),統(tǒng)計(jì)量Kendall秩次相關(guān)系數(shù)U很快收斂于標(biāo)準(zhǔn)正態(tài)分布. 給定顯著水平ɑ,其檢驗(yàn)臨界值為Uɑ/2.當(dāng)|U|>Ua/2,序列趨勢(shì)顯著;當(dāng)|U| 累積距平法利用下式來表示序列的累積距平[6]: (3) (4) 其主要思想是判斷離散數(shù)據(jù)對(duì)其均值的離散程度,若曲線呈現(xiàn)上升趨勢(shì),表明累積距平值增大,離散數(shù)據(jù)大于其均值;若曲線呈現(xiàn)下降趨勢(shì),表明累積距平值減小,離散數(shù)據(jù)小于其均值. 累積距平法是一種基于均值的檢驗(yàn)方法,通過觀察差積曲線可判斷數(shù)據(jù)點(diǎn)離散程度,以及時(shí)間序列在長(zhǎng)期內(nèi)的變化趨勢(shì)和突變時(shí)間[7].判斷系列是否存在突變:先計(jì)算回歸檢驗(yàn)值T[8],在給定α的條件下,若|T|>Ta/2,說明資料系列存在突變(跳躍),跳躍點(diǎn)為累積距平曲線的極值點(diǎn)(最大、最小)前后,具體根據(jù)系列累積距平曲線過程確定;若|T| 1960~2016年期間,鴉勒水文站以上綠汁江流域多年平均徑流量11.8億m3,多年平均徑流深166.5 mm,流域玉溪段多年平均降水量741.8 mm.降水和徑流年內(nèi)分配都極為不均,汛期(5~10月)降水占全年降水的84.7%,汛期(5~10月)徑流占全年徑流的81.4%. 根據(jù)方法(1),繪制1960~2016年年徑流、年降水、汛期(5~10月)降水和枯季(1~4月、11~12月)降水隨時(shí)間的變化過程線,建立一元線性回歸方程,統(tǒng)計(jì)各方程斜率b,以此判斷系列趨勢(shì)情況,同時(shí)根據(jù)(1)式對(duì)各要素作5年滑動(dòng)平均過程,統(tǒng)計(jì)其趨勢(shì)結(jié)論,見表1. 表1 綠汁江流域降水、徑流變化趨勢(shì)分析結(jié)果 通過一元線性回歸法和5年滑動(dòng)平均法分析得到,綠汁江流域年徑流深、年降水和汛期降水系列一元線性回歸方程的斜率b均小于0,數(shù)據(jù)系列均為減少趨勢(shì),而枯季降水系列一元線性回歸方程的斜率b大于0,說明枯季降水系列為增加趨勢(shì);5年滑動(dòng)平均分析結(jié)論也與一元線性回歸分析的結(jié)論完全一致. 根據(jù)Mann-Kendall法中(2)式進(jìn)行各分析要素的檢驗(yàn)統(tǒng)計(jì)量U計(jì)算,判斷各要素變化趨勢(shì)的顯著性(取a=0.05),見表2. 表2 綠汁江流域降水、徑流變化趨勢(shì)顯著性分析結(jié)果( Mann-Kendall法) 從統(tǒng)計(jì)分析看出,1960~2016年綠汁江流域年徑流深的減少趨勢(shì)是顯著性的,而年降水和汛期降水的減少趨勢(shì)、枯季降水的增加趨勢(shì)則不顯著. 根據(jù)(3)(4)式用累積距平法計(jì)算分析并繪制各要素累積距平過程線圖,見圖1~圖4. 圖1 1960~2016年徑流系列累積距平曲線 圖2 1960~2016年降水系列累積距平曲線 圖3 1960~2016年汛期降水系列累積距平曲線 圖4 1960~2016年枯季降水系列累積距平曲線 從1960~2016年徑流系列累積距平曲線看出,由于|T|=3.51>T(0.05/2)=1.64,說明徑流系列存在突變,突變點(diǎn)在累積距平曲線T檢驗(yàn)值最高點(diǎn)對(duì)應(yīng)的2003年,即在2003年前后徑流系列均值發(fā)生顯著性跳躍. 從1960~2016年降水系列累積距平曲線看出,由于|T|=1.62 從汛期降水系列累積距平曲線看出,由于|T|=1.76>T(0.05/2)=1.64,說明汛期降水系列存在突變,突變點(diǎn)在累積距平曲線T檢驗(yàn)值最高點(diǎn)對(duì)應(yīng)的2009年,即在2009年前后汛期降水系列均值發(fā)生顯著性跳躍. 從枯季降水系列累積距平曲線看出,由于|T|=0.95 (1)綠汁江流域徑流主要由降水形成,1960~2016年多年平均徑流量11.8億m3、多年平均徑流深166.5 mm,流域玉溪段多年平均降水量741.8 mm. (2)綠汁江流域1960~2016年年徑流深系列為顯著減少趨勢(shì),年降水系列和汛期降水系列為不顯著減少趨勢(shì),枯季降水系列為不顯著增加趨勢(shì). (3)綠汁江流域1960~2016年年徑流深系列在2003年前后發(fā)生突變,汛期降水系列在2009年前后發(fā)生突變,年降水系列和枯季降水系列無突變.2 分析成果
2.1 趨勢(shì)變化
2.2 顯著性分析
2.3 突變分析
3 結(jié) 論