大通河流域近60年徑流變化特征和趨勢分析

王大超，杜麗芳，路 賀，史偉明，張 鈺

（蘭州大學資源環境學院，甘肅蘭州 730000）

1 研究背景

20世紀80年代以來，由于人類活動和氣候變化的影響，水資源問題日益突出。中國水資源時空分布不均衡，西北地區水資源較為缺乏。河川徑流決定了地表水資源的數量，是最易開發利用的水資源［1］。河流取水、水庫蓄水、跨流域調水等人類活動和全球氣溫升高、降水量變化等氣候變化，會影響流域徑流的時空分布規律，也使得流域水文要素序列產生一定程度的變異［2-3］；導致流域水文特征頻率不一致和時空不對應，水文規律失真，使水文水利計算和水資源分析評價的成果出現偏差［4］。

大通河是西北地區少有水資源豐富的河流，水資源總量達28.95億m3，目前流域開發利用方式主要是建設梯級水電站和跨流域調水。研究表明，1990年以來大通河流域開發利用主要對生態環境和水文過程造成了負面影響［5］，流域自上至下影響程度逐漸遞增［6］；近20 a來該流域徑流時空分布規律發生了變化［7］；梯級電站建設對流域洪水過程影響較大［6］，外調水工程直接導致上游尕大灘水文站徑流減少［8］，下游享堂水文站徑流量也呈現遞減趨勢［9］，流域徑流深總體表現為減少趨勢［8］。

分析大通河流域徑流變化規律和趨勢，對其流域水資源開發利用和區域水利規劃具有重要的現實意義。目前關于大通河干流徑流特征的已有研究中，未對整個干流進行全面分析，因此本文選用該河干流天堂和享堂水文站1956～2015年60 a長系列天然徑流資料，采用滑動平均法、累計距平、M-K檢驗和Morlet小波分析等方法，分析了該河干流徑流的年內年際變化、周期、突變特征和趨勢。

2 大通河概況

大通河是湟水一級支流，位于甘肅省和青海省境內，在上游祁連和剛察縣稱默勒河，在中游門源縣名為浩門河。大通河流域海拔介于1 700～5 000 m，呈西北—東南走向，集水面積15 130 km2，其中青海省內12 943 km2。大通河發源于天峻縣托勒山的木里山泉，流經青海省祁連、剛察、海晏、門源、互助縣和甘肅省天祝、永登縣，入河口位于青海省民和縣享堂鎮，河流全長560 km，河道比降5.67%。以尕大灘站、天堂站和享堂站為控制站點（見表1），可將大通河流域劃分為上、中、下游（見圖1）。

大通河流域主要為大陸性高寒半濕潤山地氣候，徑流以大氣降水和冰川消融為主要補給來源，年平均氣溫約-4℃～-1℃，年平均降水量約400～600 mm，降水主要集中在5～10月，平均年徑流量27.7億m3（1956～2015年）。大通河流域水利工程開發較晚，從1958年就開始“引大入秦”工程的規劃勘察，到1985年開始實施水電站規劃建設［10］，至2015年已建成梯級電站42座，規劃了6處跨流域調水工程，已建成或正在建設“引大入秦”“引大濟湟”“引硫濟金”和“引大濟西”［5］，后期擬建的有“引大濟黑”和“引大濟湖”調水工程［11］。

表1 大通河流域上、中、下游水文站信息

圖1 大通河流域控制水文站分布

3 研究資料與方法

3.1 研究資料

本次主要研究大通河流域的水文特征，根據水文站控制情況，從甘肅省水文水資源局和黃河上游局收集到大通河干流中游天堂水文站、下游享堂水文站1956～2015年的長系列實測徑流資料，經檢驗，資料代表性和可靠性良好。

3.2 研究方法

3.2.1 年際變化和趨勢性分析方法

本文采用線性傾向分析法、滑動平均法和距平累計曲線法來分析大通河徑流量年際趨勢性變化。線性傾向分析法在水文氣象趨勢分析中比較常用，數理統計中通過線性回歸方程來確定自變量和因變量之間關系［12］。滑動平均法是序列經過滑動平均后削弱序列中的短周期，有效修正差異較大值引起的誤差，滑動均值隨時間的變化可以反映序列的變化趨勢［13］。距平累積曲線又稱為差積曲線，根據繪制出的曲線變化情形判斷顯著的持續性變化及趨向，可以大致判斷出序列突變的時間［14-15］。

3.2.2 小波分析方法

小波分析是20世紀80年代出現的一門新興數學分支領域，適用性較好，對具有多個時間尺度變化和非平穩變化特征的序列，能夠比較清楚地體現序列時域和頻域的局部變化特點［16］，是研究水文特征的較為可靠的方法之一。本文選用常用的非正交小波——Morlet小波對大通河流域年徑流進行周期性分析，對序列連續進行小波變換可以得出小波變換系數，進而繪制出小波分析等值線圖，從中識別水文序列在時間上的周期性特征［17］；通過繪制小波方差圖進行檢驗并識別徑流序列主周期，小波方差圖上波峰對應尺度為主周期。

3.2.3 Mann-Kendall方法

Mann-Kendall（M-K）秩次檢驗是一種非參數檢驗方法，對類型變量和順序變量具有較強的適用性和檢驗能力，可以檢測徑流時間序列的變化趨勢，以及確定趨勢起始位置，指出突變區域［18-19］。M-K檢驗一般取顯著性水平α=0.05，臨界值U0.05=±1.96，比較計算徑流時間序列UF和UB值的正負情況，得知序列呈現上升或下降趨勢［20］。當UF和UB曲線超過臨界線時，表明趨勢顯著；通過分析統計序列構成UF和UB的兩條曲線的交點是否位于臨界線之間，判斷變異是否存在，突變開始的時間就是交點對應的時刻。

圖2 1956～2015年天堂站、享堂站徑流量月平均分配曲線

4 結果與分析

4.1 年內徑流變化特征

大通河流域徑流主要補給來源是大氣降水，以降雨補給為主，融雪補給為輔。從大通河干流天堂站、享堂站徑流年內分配曲線（見圖2）可以看出，兩站的年內徑流量分配規律基本一致，各季節徑流量大小關系都呈現為：夏季>秋季>春季>冬季。大通河徑流年內分配極不均勻，受夏季和秋季降雨影響，汛期徑流主要集中在5～10月，這6個月徑流量總量占全年流量的比例約82%，11月至次年4月的徑流量總量僅占全年流量的約18%。

4.2 年際徑流變化特征和趨勢分析

根據天堂、享堂站1956～2015年徑流資料，計算得到兩站多年平均徑流量分別為24.75億m3、27.41億m3；在1989年天堂和享堂站都出現最大年徑流量，分別為40.04億、50.19億m3；天堂、享堂站最小年徑流量分別為12.68億（1973年），13.22億m3（2004年）；兩站年徑流量極值比分別為3.19，3.80，變差系數Cv值較為接近，分別為0.30，0.31，流域年徑流量序列變化較大。

20世紀90年代，大通河流域開始水利工程開發，主要在天堂站和享堂站之間。1956～1990年，大通河天堂站和享堂站多年平均徑流量分別為23.66億，29.21億m3，90年代以后兩站多年平均徑流量為24.90億，25.05億m3，天堂站多年平均徑流量相對增加了5.24%，而享堂站多年平均徑流量相對減少了14.24%。繪制大通河干流天堂站、享堂站年徑流過程及趨勢線，由圖3可知，天堂站徑流趨勢線斜率k1為0.053 9，享堂站徑流趨勢線斜率k2為-0.094 5，表明大通河上中游徑流量60 a來呈現緩慢增長的趨勢，而下游徑流量出現衰減的趨勢。

圖3 1956～2015年大通河流域天堂站、享堂站年徑流趨勢和5 a滑動曲線

由圖3天堂和享堂水文站年徑流5 a滑動曲線可以看出，1995年以前大通河流域中游和下游年徑流量豐枯水年對應，變化趨勢基本一致，1956～1980年和1990～1995年徑流量總體呈現減少趨勢，1981～1989年呈現增加趨勢；1995年以后，天堂站1996～2012年徑流量總體呈現增加趨勢，而享堂站年徑流量呈現明顯的減少趨勢。近20 a來下游影響徑流變化的因素較多，享堂站年徑流過程線比天堂站年徑流過程線變化更加劇烈。

采用累積距平法繪制近60 a來大通河干流天堂站和享堂站逐年徑流距平累計曲線，見圖4。由圖4可知，天堂站和享堂站變化趨勢差異較大，根據天堂站逐年徑流距平累計曲線，1956～1961年大通河中上游流域累計距平值持續增加，表明該時段距平值持續為正，處于豐水期；1962～1982年呈減小趨勢，處于枯水期；1982～1990年徑流量不斷增加，處于豐水期；1991～2004年，徑流量波動頻繁；2005～2013年徑流量不斷增加，處于豐水期；2014～2015年徑流量減少，處于枯水期。

根據享堂站逐年徑流距平累計曲線可知，1956～1968年處于豐水期；1969～1980年徑流量波動頻繁，呈減小趨勢；1981～1990年徑流量又不斷增加，處于豐水期；1991年開始累計距平值呈減小趨勢，表明該時段距平值持續為負，處于枯水期。

圖4 1956～2015年大通河流域天堂站、享堂站年徑流序列標準距平累計曲線

4.3 大通河干流年徑流周期變化

對天堂、享堂水文站年徑流序列采用Morlet小波分析法進行周期性分析并通過小波方差檢驗，小波分析結果如圖5所示。大通河干流天堂站年徑流量存在3，7 a及25 a左右的振蕩周期，徑流序列最大峰值均出現在25 a，表明25 a的周期振蕩最強，是徑流量序列的第一主周期；7 a和3 a時間尺度分別為徑流量變化的第二和第三主周期。在25 a時間尺度的情況下，正負相位交替明顯，經歷了“豐—枯—豐—枯—豐”5個變化期，其中1956～1966年、1978～1990年、2004～2014年為正相位，為豐水年；1967～1977年、1991～2003年為負相位，為枯水年。在7 a時間尺度下，經歷了8個“豐—枯”轉換期。根據上述分析可知：當時間尺度不同，徑流量豐枯交替變化也會不同，徑流量周期變化以第一周期為主；徑流量等值線圖在2015年以后表現為負相位，且等值線顯示未封閉，表明研究區徑流量將繼續減少；按照25 a的第一主周期特征推算，年徑流量到2028年左右將持續呈現減少趨勢。

在大通河下游享堂站年徑流周期性分析圖中，大通河下游享堂站存在3，6 a及30 a左右的振蕩周期，在30 a時間尺度上經歷了“豐—枯—豐—枯—豐”5個變化期，1956～1965年、1981～1996年和2012～2015 年為偏多期，為豐水年；1966～1980年和1997～2011年為偏少期，為枯水年。下游徑流變化主周期大于上游，小周期轉換要比上游頻繁，大通河上下游徑流量變化周期年份對應關系不好，可能與流域徑流補給形式、水資源開發程度和調蓄能力有關，需要進一步研究這種情況的影響因素。

圖5 大通河流域天堂站、享堂站年徑流量小波分析和小波方差

4.4 大通河干流年徑流突變性分析

運用M-K檢驗法對大通河干流天堂站、享堂站年徑流序列進行突變檢驗，假定年徑流變化穩定、時間序列相互獨立而且保持相同連續分布，取顯著性水平α=0.05，臨界值U0.05=±1.96，天堂、享堂水文站年徑流突變分析見圖6。根據圖6可知天堂、享堂兩站年徑流量變化趨勢為：1990年以前基本一致；1956～1960年UF>0，年徑流量上升；1962～1987年UF＜0，年徑流量下降。天堂水文站年徑流統計量曲線UF、UB交點在置信區間里面，突變發生在1989年左右；1988～2015年UF>0，年徑流量上升趨勢，其中2011～2015年顯著上升。享堂水文站UF、UB曲線在1999年置信區間內相交，表明享堂站徑流量在1999年發生突變；1996～2015年UF＜0，徑流量下降，其中2010～2015年呈顯著下降趨勢。

圖6 大通河流域天堂站、享堂站年徑流突變分析

5 結論

本次分析研究采用滑動平均法、累計距平、小波分析法和M-K秩次檢驗方法，系統分析大通河中游天堂水文站和下游享堂水文站年徑流變化特征和趨勢。結果表明：

（1）大通河流域中游和下游年徑流量變化在1956～1995年變化趨勢相似，豐枯水年基本相對應；1995～2015年兩站徑流變化趨勢相反。近60 a，中上游年徑流量總體呈上升趨勢，但增加趨勢不明顯，以0.053 9億m3/a的幅度增加；下游享堂站年徑流量存在明顯下降趨勢，以0.094 5億m3/a的幅度減少。

（2）采用Morlet小波分析對天堂站和享堂站1956～2015年年徑流量序列進行周期性分析，結果表明天堂站年徑流存在3，7 a及25 a左右的振蕩周期，在25 a時間尺度上經歷了“枯-豐-枯-豐-枯”5個循環交替；下游小周期的轉換要比上游頻繁，享堂站年徑流量存在3，6 a及30 a左右的振蕩周期。

（3）大通河天堂站和享堂站徑流變化1990年以前趨勢基本一致，中游天堂站年徑流量在1989年發生了突變，1989～2015年年徑流量呈現上升趨勢；而下游享堂站徑流在1999年發生突變，1996～2015年年徑流量呈現下降趨勢，2010年后顯著下降。