精河流域徑流變化特征及其對降水變化的響應

董煜 張立山 陳學剛

摘要：利用精河水文站1964年-2012年徑流、降雨數(shù)據(jù)，采用M-K檢驗、累積曲線模型及交叉小波分析等方法，對艾比湖流域精河的徑流、降水的變化特征及多時間尺度相關進行分析，并定量分析了降水變化對徑流變化的貢獻率。精河年徑流量總體上徑流呈增長趨勢，但趨勢不明顯；降雨總體上呈增加趨勢，且增加趨勢明顯。精河年徑流量與降水量在多時間尺度上存在顯著相關關系，且以正相關關系為主，顯著相關區(qū)域集中在1986年-1991年在2～4 a的周期帶、1985年-1990年的7 a周期帶、2002年-2011年的 4～5 a周期帶。分析認為，近50年來該流域降水變化對徑流增加的貢獻率平均為68.4%。

關鍵詞：徑流；降水；M-K檢驗；交叉小波分析；精河流域

中圖分類號：P332 文獻標志碼：A 文章編號：1672-1683（2016）04-0060-05

Abstract：Based on the annual runoff and precipitation data of Jinghe Hydrological Station from 1964 to 2012 in Jinghe River，the paper studied the changing trend and multi-time scale correlation between the annual runoff and precipitation by using the methods of Mann-Kendall test，cross-wavelet test and double mass curve，and the effects of the precipitation change on runoff was also quantitatively analyzed.The results showed that runoff increased insignificantly and precipitation increased significantly，respectively.The multi time scale correlativity was significant between the annual runoff and precipitation in Jinghe River，which was dominated by positive correlation.The significantly correlated regions were the 2～4-year periods during 1986-1991，the 7-year periods during 1985-1990，and the 4～5 year periods during 2002-2011.The result of double mass analysis showed that the increases of runoff in JingheRiver was mainly attributed to precipitation change，and its contribution was 68.4%.

Key words：runoff；precipitation；Mann-Kendall test；cross-wavelet analysis；Jinghe River

河川徑流的時空變化直接影響區(qū)域社會經(jīng)濟的發(fā)展及生態(tài)環(huán)境安全[1-3]，其受氣象條件、下墊面性質(zhì)、地形地貌以及社會經(jīng)濟活動等多種要素的影響 [4-5]，同時具有確定性變化和隨機性變化特性。水資源短缺已成為我國西北干旱地區(qū)社會可持續(xù)發(fā)展和生態(tài)環(huán)境保護的主要限制因素。已有研究表明，該地區(qū)氣候由暖干向暖濕化變化的趨勢[6-7]，已對區(qū)內(nèi)徑流產(chǎn)生了明顯的影響[ 8-10]。因此在氣候暖濕變化的背景下，有必要研究干旱地區(qū)地表徑流的變化特征及氣候因素對徑流的影響程度，這對于干旱地區(qū)合理利用水資源、有效改善當?shù)氐纳鷳B(tài)環(huán)境、建立先進的現(xiàn)代水資源管理觀念、建設合理符合實際的干旱地區(qū)流域綜合管理體系具有現(xiàn)實意義。

艾比湖作為新疆最大的咸水湖，流入的河流因中上游生產(chǎn)、灌溉而大量引水，除博爾塔拉河、精河存在入湖徑流外，其他入湖河流如奎屯河、四棵樹河、古爾圖河等已斷流 [11]。入湖水量減少，導致湖面積萎縮，湖西北部的湖底干枯，使得鹽塵影響到整個艾比湖流域生態(tài)安全，對流域農(nóng)牧業(yè)、工業(yè)生產(chǎn)及人體健康造成危害[12]。近年針對艾比湖流域入湖徑流開展眾多研究，主要涉及博爾塔拉河、精河徑流與氣溫、降雨關系、年徑流變化等問題[13-19]。但對于精河月徑流的變化特征趨勢及是否存在突變變化、徑流對降水的定量響應研究相對較少。本文采用累積距平、M-K趨勢檢驗、交叉小波分析等方法，對徑流與降雨在年際多時間尺度以及徑流對降水變化的水文響應進行分析研究。

1 研究區(qū)域及資料來源

精河起源于天山山脈婆羅科努山北麓，由南向北流入艾比湖[17]，整個流域位于精河縣境內(nèi)[20]，流域面積2 150 km²，河流全長114 km，多年平均年徑流量4.75×10⁸m³，流域范圍介于東經(jīng)81°46′-83°51′E、北緯44°02′-45°10′N。烏圖精河與冬都精河是精河的主要支流，流域地形由南向北部傾斜，主要由山地、洪積地貌、沖積-湖積地貌構(gòu)成，精河縣農(nóng)業(yè)生產(chǎn)區(qū)主要位于沖積-湖積平原（圖1）。該流域位于中緯度地帶，屬于溫帶干旱荒漠類型的大陸性氣候，春季氣溫回升快且不穩(wěn)定，多低溫降雪天氣，秋季降溫迅速。大風天氣多出現(xiàn)在春秋兩季，給工農(nóng)業(yè)生產(chǎn)、社會生活帶來諸多影響。

研究數(shù)據(jù)來源于精河水文站1964年-2012年的逐月徑流實測數(shù)據(jù)及降水量逐月實測數(shù)據(jù)（精河水文站位于44°22′N，82°55′E，海拔為619.2 m）。

2 研究方法

Mann-Kendal是一種非參數(shù)統(tǒng)計檢驗方法，亦稱無分布檢驗，主要特點是不需要樣本遵從一定的分布，不受少數(shù)異常值的干擾，計算簡單、檢測范圍寬、定量化程度高，是突變檢測方法中常用的方法，具體計算過程見文獻[21]。利用在Mann-Kendall進行單調(diào)趨勢檢驗時，統(tǒng)計量Zc為正值時，表示研究序列呈增加趨勢，Zc為負值時表示研究序列呈減少趨勢。

運用Morlet小波變化研究精河徑流與降雨的多時間尺度的變化特征，小波能量譜定義為

交叉小波分析由交叉譜分析與小波變換兩種方法結(jié)合而產(chǎn)生的一種在時頻域中分析兩個信號相關性的分析方法[24]，通過降雨和徑流序列的小波系數(shù)相乘計算交叉小波譜，可以將降雨和徑流的相關性關系在不同的時間尺度上顯示出來[25]。

運用交叉小波變換譜分析，進一步從多時間尺度的角度研究徑流與降雨在時頻域中的相關性，并運用置信水平為95%的紅噪聲總體譜進行顯著性檢驗。

流域降水與徑流變化是否一致可以通過徑流-降水雙累積曲線直觀表示，雙累積曲線呈現(xiàn)為直線，表示降水與徑流變化存在一致性，如果直線發(fā)生偏移，則表示研究流域的降水徑流特征發(fā)生了變化[26]。

3 結(jié)果分析

3.1 精河徑流量與降水量的趨勢分析

精河年均徑流量（圖2）為4.75×10⁸m³，最大徑流量出現(xiàn)在1988年，為6.01×10⁸m³，最小徑流值出現(xiàn)在1992年，為3.68×10⁸m³，絕對變化幅度為2.33×10⁸m³。1964年至2012年間，年徑流量呈波動變化，但總體上呈增加趨勢，年徑流增加率為0.176 m³/（s·10a）。流域內(nèi)降水量呈波動增加趨勢，平均每年增加0.990 8 mm。多年平均降雨量為150.4 mm，最大降水量出現(xiàn)在2003年（249.8 mm），最小值出現(xiàn)在2008年（98.8 mm），絕對變化幅度為151 mm。通過圖2的降雨量及徑流量線性趨勢變化分析，精河徑流量的變化與降水量的變化趨勢大致相符。

根據(jù)M-K趨勢分析，年徑流量的檢驗統(tǒng)計量Zc =0.844 74Z0.05=1.96，增加趨勢明顯。

在圖3中，精河9月份徑流量檢驗統(tǒng)計量大于0.05水平下的顯著性檢驗，這表明在1964年至2012年49年序列中9月份徑流變化呈顯著增加趨勢。同理，11月份徑流變化與9月相同。2月、7月、8月的徑流量的Zc值處于于0與-1.96之間，表示這三個月徑流量雖然呈減少趨勢，但趨勢變化不顯著。而其他月份的檢驗統(tǒng)計量值介于0與1.96之間，表明徑流量呈現(xiàn)增加趨勢但不顯著。對于降水量，2月份降雨在49年序列中表現(xiàn)為顯著增加，其他月份在49年序列沒有顯著變化。

3.2 精河徑流及降雨的多尺度相關分析

3.2.1 徑流量與降水量年內(nèi)分配特征

從精河域徑流的年內(nèi)分配特征來看（圖4），徑流主要集中在夏季的6、7、8三個月，占到年徑流量的63%；在1至5月、9至12月的徑流較少，最大月徑流量為7月。

精河流域降水的年內(nèi)分配特征顯示，降水主要集中在4月-8月，5個月的降水占到年降水量的69%，1月-3月、9月-12月的降水量較少，與精河最大徑流量月相比，兩者的月分配特征變化存在差異，月降水變化最大值比徑流的變化提前。

3.2.2 徑流量與降水量多尺度特征

應用交叉小波對精河徑流與降水的變化特征及二者不同時間尺度相關關系進行分析，進一步加深對精河的水文特征了解。

對精河流域徑流、降水采用Morlet小波函數(shù)進行小波變換，得到年徑流及降水的連續(xù)小波譜，見圖5。圖中深色區(qū)為95%置信區(qū)，表示功率譜與置信水平為95%的紅噪聲標準譜的比值大于1。

從圖5中可以看出，在95%的置信度下：年均徑流量存在3個顯著周期，分別為1～5 a左右（1964年-1973年）、1～6 a左右（2001年-2012年）和1～7 a左右（1982年-1996年）的周期；降雨只存在1～7 a左右的顯著周期（1977年-2012年）。徑流與降雨在時頻域中都存在著顯著的、不同尺度的周期。徑流在1～6 a及1～7 a尺度的周期顯著性最高，周期尺度范圍較大。降水的顯著周期與徑流存在一定的差異，通過置信度95%的顯著性的區(qū)域，降水周期較徑流周期長，但兩者高能量區(qū)表現(xiàn)出比較好的一致性，且降雨周期尺度涵蓋了徑流的大部分尺度的周期，在小波分析中高能量譜的范圍基本包含了徑流和降水的高能量區(qū)，說明流域徑流變化對降水變化具有積極的響應。

圖6中，用實線和虛線分別表示正、負相關關系，點劃線表示正相關與負相關關系界線，反映精河年徑流量與降水量之間在不同時間尺度下相關性關系。深色區(qū)域表示小波功率譜與置信水平為95%的紅噪聲總體譜的比值大于1，表示徑流與降雨呈顯著性相關關系。

從交叉小波譜看，1964年-2012年精河年降水量與徑流量在不同時間尺度上以正相關關系為主。

3.3 精河降雨變化對徑流的影響分析

圖7顯示，自1980年開始，精河徑流-降水雙累積曲線發(fā)生偏移。可以將曲線發(fā)生偏移前的時間設定為基準期，表示精河徑流受自然因素影響，為天然徑流量。確定基準期為1964年至1980年，并根據(jù)累積曲線對精河徑流劃分不同時段。

依據(jù)基準期年降水和年徑流資料，建立基準期內(nèi)的徑流序列R和降水序列P的相關方程，其方程為：R=0.9152P +178.43，r=0.962。

由年降水量與年徑流量的相關方程，計算得出精河不同時段的徑流量，可近似代表天然徑流量。用基準期實測徑流數(shù)據(jù)與不同時段的天然徑流量的差值，為此時段降水變化對徑流變化的影響值；基準期實測值與各時段實測值的差值減去降水變化的影響值即為其他因素對徑流變化的影響值。分析表1可以看出，降水是徑流變化的主要影響因素，影響率平均大于68%，這與相關研究成果一致[17，20]。因為數(shù)據(jù)來源的水文站為山口水文站，其上游受人類影響相對綠洲帶要較小。其他影響因素可能是溫度、冰川積雪、地表覆被、蒸發(fā)量等原因。

4 結(jié)語

（1）精河年徑流量呈增長趨勢，但趨勢不明顯；年降水量總體上呈增加趨勢，且增加趨勢明顯。在1964年至2012年49年中，徑流在9月、11月，降雨在2月呈現(xiàn)增加趨勢，徑流降雨在其余月份未有顯著增加或減少變化。

（2）精河年徑流量與降水量在多時間尺度存在顯著相關關系，并以正相關為主，顯著相關區(qū)域主要集中在1986年-1991年在2～4 a的周期帶、1985年-1990年的7 a周期帶、2002年-2011年的4～5 a周期帶。

（3）通過運用降水-徑流雙累積曲線分析了降水對徑流變化的影響率，說明49年來，精河降水是導致精河徑流發(fā)生變化的主要因素，對徑流影響的貢獻率平均為68.4%。

參考文獻（References）：

[1] ZHANG Q，XU C-Y，TAO H，et al.Climate changes and their impacts on water resources in the arid regions：a case study of the Tarim River Basin，China[J].Stochastic Environmental Research and Risk Assessment，2010，24（3）：349-358.

[2] 山成菊，董增川，付曉花，等.近 50 年來灤河流域徑流演變規(guī)律研究[J].南水北調(diào)與水利科技，2013，11（4）：5-8.（SHAN Cheng-ju，DONG Zeng-chuan，F(xiàn)U Xiao-hua，et al.Rules of runoff variation in Luanhe River Basin in recent 50 years[J].South-to-North Water Transfers and Water Science & Technology，2013，11（4）：5-8.（in Chinese））

[3] YANING C，CHANG C X，XING M H，et al.Fifty-year climate change and its effect on annual runoff in the Tarim River Basin，China[J].Quaternary International，2009，208（1）：53-61.

[4] MA X，YASUNARI T，OHATA T，et al.Hydrological regime analysis of the Selenge River Basin，Mongolia[J].Hydrological Processes，2003，17（14）：2929-2945.

[5] XU D，ZHANG G，YIN X.Runoff varation and its impacting factor in Nenjiang River during 1956-2006[J].Advances in Water Science，2009，20（3）：416-421.

[6] 施雅風，沈永平，李棟梁，等.中國西北氣候由暖干向暖濕轉(zhuǎn)型的特征和趨勢探討[J].第四紀研究，2003，23（2）：152-164.（SHI Ya-feng，SHEN Yong-ping，LI Dong-liang，et al.Discussion on the present climate change from warm-dry to warm-wet in northwest[J].Quaternary Sciences，2003，23（2）：152-164.（in Chinese））

[7] LI Q，CHEN Y，SHEN Y，et al.Spatial and temporal trends of climate change in Xinjiang，China[J].Journal of Geographical Sciences，2011，21（6）：1007-1018.

[8] CHEN Y，TAKEUCHI K，XU C，et al.Regional climate change and its effects on river runoff in the Tarim Basin，China[J].Hydrological Processes，2006，20（10）：2207-16.

[9] XU J，CHEN Y，JI M，et al.Climate change and its effects on runoff of Kaidu River，Xinjiang，China：a multiple time-scale analysis[J].Chinese Geographical Science，2008，18（4）：331-9.

[10] LING H，XU H，SHI W，et al.Regional climate change and its effects on the runoff of Manas River，Xinjiang，China[J].Environmental Earth Sciences，2011，64（8）：2203-13.

[11] 賈春光，王曉峰，金海龍，等.新疆艾比湖湖面動態(tài)變化及其影響研究[J].干旱區(qū)資源與環(huán)境，2006，20（4）：152-156.（JIA Chun-guang，WANG Xiao-feng，JIN Hai-long，et al.Research on the dynamic changes of Ebinur Lake Area and the effects.[J].Journal of Arid Land Resources and Environment，2006，20（4）：152-156.（in Chinese））

[12] 李艷紅，楚新正，金海龍.新疆艾比湖流域水文特征分析[J].水文，2006，26（5）：69-71.（LI Yan-hong，CHU Xin-zheng，JIN Hai-long.Study on changes of hydrological characteristics of Ebinur Lake Basin in Xinjiang.[J].Journal of China Hydrology，2006，26（5）：69-71.（in Chinese））

[13] 吉力力·阿不都萬里，穆桂金.艾比湖干涸湖底塵暴及其災害分析.[J].干旱區(qū)地理，2002，25（2）-149-154.（JILIL Abduvali，MU Gui-jin.Analysis on the dust storms and their disasters in the Lake bed region of Ebinur Lake，XinJiang.[J].Arid Land Geography，2002，25（2）-149-154.（in Chinese））

[14] 李加強，陳亞寧，李衛(wèi)紅，等.天山北坡中小河川降水與徑流變化特征--以精河為例[J].干旱區(qū)地理，2010，33（4）：615-622.（LI Jia-qiang，CHEN Ya-ning，LI Wei-hong，et al.Variation features of precipitation and runoff of the middle-small rivers of Northern Piedmont of Tianshan Mountains：a case of Jinghe River[J].Arid Land Geography，2010，33（4）：615-622.（in Chinese））

[15] 劉世薇，周華榮，梁雪瓊，等.艾比湖流域降水與徑流變化特征分析[J].水土保持學報，2011，25（5）：21-25.（LIU Shi-wei，ZhouHua-rong，Liang Xue-qiong，et al.Trend analysis of the precipitation and runoff in Ebinur Lake Basin[J].Journal of Soiland Water Conservation，2011，25（5）：21-25.（in Chinese））

[16] 喬木，周生斌，盧磊.艾比湖流域年徑流變化特征分析[J].水土保持學報，2010，24（6）：236-239.（QIAO Mu，ZHOU Sheng-bin，LU Lei.Trends in runoff variations of the Ebinur Lake Basin during the last 48 years[J].Journal of Soiland Water Conservation，2010，24（6）：236-239.（in Chinese））

[17] 柯琳，葉茂，凌紅波，等.近 48 年區(qū)域氣候變化對精河地表徑流的影響分析[J].水資源與水工程學報，2013，24（1）：54-59.（KE Lin，YE Mao，LING Hong-bo，et al.Analysis of regional climate change impact on surface runoff in Jinghe River in recent forty-eight years[J].Journal of Water Resources&Water Engineering，2013，24（1）：54-59.（in Chinese））

[18] 司桂琴，劉志輝.新疆精河氣溫，降水與徑流量近 51a 變化分析[J].新疆農(nóng)業(yè)科學，2010，47（009）：1792-1796.（SI Gui-qin，LIU Zhi-hui.Analysis on the change of temperature，precipitationand runoff in Jinghe，Xinjiang during the recent 51years[J].Xinjiang Agricultural Sciences，2010，47（009）：1792-1796.（in Chinese））

[19] 董煜，海米提·依米提.艾比湖流域徑流水文特征及其對降水變化響應-以博爾塔拉河為例[J].水土保持研究，2014，21（2）：94-99.（DONGYu，HAIMITI·Yimiti.Runoff characteristics and responses to precipitation in the Ebinur Lake Basin-A case study of Boertala River.[J].Research of Soil and Water Conservation，2014，21（2）：94-99.（in Chinese））

[20] 高明.精河徑流與降水變化特征及響應[J].新疆水利，2011，（2）：9-12.（GAO Min.Characteristics and responses analysis of runoff and precipitation in JingHe River.[J].Xinjiang Water Resources，2011，（2）：9-12.（in Chinese））

[21] 高宇，袁勇，方宏陽.近50年黑河流域降水變化特征分析[J].南水北調(diào)與水利科技，2013，11（5）：16-21.（GAO Yu，YUAN Yong，F(xiàn)ANG Hong-yang.Analysis of precipitation variation in the Heihe River in recent 50 years[J].South-to- North Water Transfers and Water Science & Technology，2013，11（5）：16-21.（in Chinese））

[22] LABAT D.Cross wavelet analyses of annual continental freshwater discharge and selected climate indices[J].Journal of Hydrology，2010，385（1-4）：269-278.

[23] HAO Y，LIU G，LI H，et al.Investigation of karstic hydrological processes of Niangziguan Springs （North China）using wavelet analysis[J].Hydrological Processes ，2012，26（20）：3062-3069.

[24] TORRENCE C，COMPO G P.A practical guide to wavelet analysis[J].Bulletin of the American Meteorological Society，1998，79（1）：61-78.

[25] Nason G P，Sapatinas T.Wavelet packet transfer function modelling of nonstationary time serves[J].Statistics and Computing，2002，12（1）：45-56.

[26] 王隨繼，聞云霞，顏明，等.皇甫川流域降水和人類活動對徑流量變化的貢獻率分析[J].地理科學，2012，67（3）：388-397.（WANG Sui-ji，YAN Yun-xia，YAN Ming，et al.Contributions of precipitation and human activities to the runoff change of the Huangfuchuan Drainage Basin [J].Acta Geographica Sinica，2012，67（3）：388-397.（in Chinese））