鄢 波，夏自強(qiáng)，3，周艷先，王景才，陳起川

(1.河海大學(xué)水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室，2.河海大學(xué)水文水資源學(xué)院，3.河海大學(xué)國際河流研究所，江蘇南京 210098)

位于我國東北的黑龍江(俄羅斯稱阿穆爾河)流經(jīng)蒙古東部、中國東北地區(qū)、俄羅斯遠(yuǎn)東地區(qū)南部，是東亞和東北亞最大的河流之一，是一條國際河流，黑龍江全長4 370 km，它的中段為中俄界河，流經(jīng)黑龍江省長達(dá)2 900 km［1］。近年來由于邊境地區(qū)經(jīng)濟(jì)的迅速發(fā)展，中俄界河正扮演越來越重要的角色，中俄界河的水環(huán)境問題、水資源安全問題、跨境生態(tài)問題受到越來越多的重視［2］。近年來，各國學(xué)者對黑龍江流域的降水以及徑流等開展了大量的研究，取得了一些重要成果［1，3-5］。

哈巴羅夫斯克站(以下簡稱為哈巴站)地理位置重要，研究哈巴站的徑流變化特征、規(guī)律，對于分析黑龍江下游的水文情勢具有重要的意義，可為界河水資源的開發(fā)和保護(hù)提供科學(xué)的依據(jù)。該站的長期徑流系列是全流域各站展延的重要依據(jù)，尤其是對于資料相對短缺的國內(nèi)各水文站點(diǎn)。由于近代黑龍江流域的人為干擾相對較少，依據(jù)該站的長期徑流系列可以進(jìn)一步分析研究全球氣候變化對黑龍江流域水資源的影響。

1 研究區(qū)概況

哈巴站位于東經(jīng)135°10'、北緯48°31'，地處黑龍江干流中下游的分界處，集水面積163.0 萬km2，占黑龍江流域面積的87.9%，是黑龍江下游重要的控制站，也是黑龍江干流接納結(jié)雅河、布列亞河、松花江和烏蘇里江等主要支流后的主要控制站，在流域中所處的位置見圖1 所示。徑流資料選用哈巴站1897—2005 年共109 a 年徑流資料和1897—1985年共89 a 月徑流資料。

該水文站所在的黑龍江流域?qū)俅箨懶约撅L(fēng)氣候，冬季氣候嚴(yán)寒，氣溫低于0℃的時間長達(dá)6 個月之久，河流封凍時間170 ～200 d，北部有永久性島狀凍土分布區(qū)［6］。流域內(nèi)夏季短暫，降水集中，流域產(chǎn)匯流條件良好，徑流量集中且豐富，徑流系數(shù)較大。徑流的補(bǔ)給以降水為主，約占65% ～80%;融雪徑流為輔，約占15% ～20%;地下徑流補(bǔ)給不足10%［6］。傳統(tǒng)的四季劃分不能體現(xiàn)黑龍江流域的季節(jié)特性，也不利于各個季節(jié)的水文分析。采用吳瓊等提出的四季劃分來分析黑龍江干流各個季節(jié)的徑流情況，該劃分為:冬季漫長，11 月—次年3 月;夏季適中，6—8月;春秋過渡季節(jié)短促，分別為4—5 月，9—10 月［7］。

圖1 哈巴站在黑龍江流域的位置

2 徑流的變化特征分析

2.1 徑流年際變化分析

2.1.1 年際和年代際變化特征

選用Mann-Kendall 法對哈巴站年徑流進(jìn)行趨勢檢驗(yàn)。Mann-Kendall 法最初由Mana［8］于1945 年提出，是檢驗(yàn)序列變化趨勢的有效工具，被廣泛應(yīng)用于氣象參數(shù)和水文序列的分析［9］。該檢驗(yàn)法不需要樣本遵循一定的分布，也不受個別異常值的干擾，具體公式見文獻(xiàn)［10］。

年際變化上，綜合圖2 及表1，可以看出哈巴站多年平均徑流量為2 604.06 億m3，年徑流量最豐年為1897 年的4 221.7 億m3，最枯年為1979 年的1352.9 億m3，最豐、最枯年之比為3.12，年徑流量的變差系數(shù)Cv為0.225。通過5 a 滑動平均過程線可以較為直觀地看出哈巴站的年徑流趨勢變化。從趨勢線來看，哈巴站的年徑流在這個統(tǒng)計時段內(nèi)呈現(xiàn)微弱下降趨勢，減少速率為3.068 億m3/a。年徑流量經(jīng)過M-K 趨勢檢驗(yàn)Z = -1.482，未通過95%顯著性檢驗(yàn)，減少趨勢不顯著。

圖2 黑龍江哈巴站年徑流量變化過程

表1 黑龍江哈巴站各年代際徑流量統(tǒng)計

年代際變化上，只有20 世紀(jì)20、70、80 年代和21 世紀(jì)初的多年均值小于整個序列的多年均值。其中，19 世紀(jì)末徑流量最大，為3 010.07 億m3，其次是20 世紀(jì)50 年代，為2 961.71 億m3，21 世紀(jì)初最小，為1 946.82 億m3，各年代徑流量經(jīng)歷了一個減少—增加—減少的變化過程，而各年代年徑流量變差系數(shù)Cv也處于0.15 ～0.35 這個區(qū)間里。

2.1.2 徑流豐枯分析

所謂連豐年是指該年段內(nèi)各年都屬于特豐水年或偏豐水年，連枯年是指該年段內(nèi)各年都屬于特枯水年或偏枯水年［2］。豐水年相應(yīng)頻率P ＜37.5%，枯水年相應(yīng)頻率P ＞62.5%［2］。根據(jù)標(biāo)準(zhǔn)，計算出連豐年和連枯年的時段均值及其與整個序列多年均值的比值(模比系數(shù))K豐和K枯［2］，統(tǒng)計結(jié)果見表2。由表2 可以看出，哈巴站依次出現(xiàn)3 次連豐年，持續(xù)時間分別為3 a、9 a、3 a，K豐值分別為1.20、1.31、1.14;出現(xiàn)4 次連枯年，持續(xù)時間分別為3 a、4 a、7 a、7 a，K枯值分別為0.84、0.74、0.72、0.77。

表2 黑龍江哈巴站徑流量連豐年和連枯年統(tǒng)計值

2.1.3 徑流變化的階段性分析

為了更好地表示徑流年際變化的階段特征，采用距平累積法［11］進(jìn)行分析。即

式中:Pt為第t 年的距平累積值;Ri為第i 年的徑流量;ˉR 為徑流序列平均值;n 為序列長度。

根據(jù)哈巴站年徑流序列距平累積年際變化過程(圖3)，可將哈巴站年徑流變化過程分為如下階段(表3):3 個豐水段即1928—1938 年、1955—1963年、1983—1985 年;2 個平水段即1940—1953 年、1986—1997 年;3 個枯水段即1916—1927 年、1974—1980 年、1999—2005 年。前面分析的連豐年和連枯年均包含在相應(yīng)的階段當(dāng)中。距平累積曲線不同程度的小波動說明在任何階段豐水年和枯水年的交錯出現(xiàn)的普遍性。

圖3 黑龍江哈巴站年徑流量距平累積變化過程

2.2 徑流年內(nèi)分配分析

2.2.1 徑流年內(nèi)分配過程分析

降水是黑龍江徑流的主要補(bǔ)給來源，流量過程與雨量過程也基本相應(yīng)。黑龍江流域降水受氣候和地形的影響，年內(nèi)降水主要集中在6—9 月份，約占全年降水的70%［7］。冬季僅占10%左右［7］。

黑龍江徑流的年內(nèi)分配四季分明，徑流年內(nèi)分配極不均勻，從哈巴站徑流年內(nèi)分配曲線(圖4)和哈巴站四季徑流變化過程(圖5)可以看出，全年徑流主要集中在夏(6—8 月)、秋(9—10 月)兩季，占全年的76.01%。夏秋兩季徑流呈下降趨勢是導(dǎo)致20 世紀(jì)末至21 世紀(jì)初枯水年份增多的一個重要因素。夏秋兩季是洪水發(fā)生的時期，其中汛期(6—9月)的徑流量占全年的63.85%。冬季(11—次年3月)進(jìn)入冰期，徑流主要依靠地下水補(bǔ)給，各月徑流量都很小，冬季徑流量只占全年的8.26%。春季(4—5 月)進(jìn)入春汛，氣溫回升，進(jìn)入春季后上游各支流先后開河，流域的積雪融化和河網(wǎng)儲冰解凍形成春汛，但春汛水量不大，春季徑流量占全年的15.73%［2］。

圖4 黑龍江哈巴站徑流年內(nèi)分配過程

2.2.2 徑流年內(nèi)分配特征分布

為進(jìn)一步定量分析黑龍江徑流年內(nèi)分配特征，采用年內(nèi)分配不均勻系數(shù)Cu、完全調(diào)節(jié)系數(shù)Cr［12-14］、集中度Cn［15］等指標(biāo)，從不同角度分析徑流年內(nèi)分配的變化規(guī)律。年內(nèi)分配不均勻系數(shù)越大，表明年內(nèi)各月徑流相差越懸殊，年內(nèi)分配越不均勻，完全調(diào)節(jié)系數(shù)越大，徑流年內(nèi)分配越集中，不均勻程度越高。集中度反映徑流年內(nèi)分配集中的程度。這3 個指標(biāo)的具體計算公式可參考文獻(xiàn)［16-19］。

首先對黑龍江哈巴站徑流逐年代際的月徑流按多年平均值計算，然后按以上3 個指標(biāo)進(jìn)行計算，逐年代際的徑流年內(nèi)分配特征指標(biāo)值見表4。由表4 可見，1897—1985 年哈巴站徑流年內(nèi)分配Cu為0.844，Cr為0.387，Cn為0.553。隨著時間推移，Cu和Cr均呈現(xiàn)增加-減少-增加的過程，兩值變化過程一致，兩值在20 世紀(jì)60 年代均最小，說明60 年代徑流年內(nèi)分配最為緩和，在20世紀(jì)10年代兩值最大，年內(nèi)分配最不均勻;Cn在20 世紀(jì)10 年代最大，20世紀(jì)60 年代最小，也正好對應(yīng)Cu和Cr的最大值和最小值。

表3 黑龍江哈巴站年徑流量變化階段

圖5 黑龍江哈巴站1897—1985 年四季徑年際變化過程

表4 黑龍江哈巴站徑流逐年代際的年內(nèi)分配特征指標(biāo)值

綜上，徑流變化的原因:一方面可能是受氣候因素如流域降水影響的結(jié)果，黑龍江流域降水的年際變化較大且階段性明顯，連續(xù)多雨、連續(xù)少雨交替出現(xiàn);另一方面還可能受到上游人類活動如農(nóng)田灌溉、水庫調(diào)節(jié)、工業(yè)用水等的影響。

3 結(jié) 語

a. 研究分析了黑龍江哈巴站徑流年際變化規(guī)律、徑流年內(nèi)分配特征，可為深入了解黑龍江流域水資源的變化趨勢，為流域的水利工程建設(shè)、生態(tài)環(huán)境保護(hù)和水資源合理利用，以及界河的開發(fā)和保護(hù)提供科學(xué)的依據(jù)。

b. 黑龍江哈巴站徑流連豐年、連枯年交替出現(xiàn)，徑流年際變化的階段性較為明顯:3 個豐水段，即1928—1938、1955—1963、1983—1985 年;2 個平水段，即1940—1953、1986—1997 年;3 個枯水段，即1916—1927、1974—1980、1999—2005 年。

c. 哈巴站徑流年內(nèi)分配極不均勻，全年徑流主要集中在夏秋兩季(6—10 月)，占全年的76.01%，其中汛期(6—9 月)占全年的63.85%。哈巴站在20 世紀(jì)60 年代徑流年內(nèi)分配最不集中、最為均勻，10 年代年內(nèi)分配最為集中、最不均勻。

d. 筆者僅應(yīng)用統(tǒng)計檢驗(yàn)方法對徑流系列的變化做了初步分析研究，有關(guān)后繼成果將陸續(xù)發(fā)表。

［1］葉柏生，李翀，楊大慶，等.我國過去50a 來降水變化趨勢及其對水資源的影響(Ⅰ):年系列［J］. 冰川凍土，2004，26(5):587-594. (YE Baisheng，LI Chong，YANG Daqing，et al. Variation trend of precipitation and its impact on water resources in China during last 50 years(Ⅰ):annual variation［J］. Journal of Glaciology and Geocryology，2004，26(5):587-594.(in Chinese))

［2］陸志華，夏自強(qiáng)，于嵐嵐.松花江佳木斯站徑流變化規(guī)律及演變趨勢分析［J］. 水電能源科學(xué)，2011，29(4):14-17.(LU Zhihua，XIA Ziqiang，YU Lanlan. Analysis of change rules and variation trend of runoff of Jiamusi Station in Songhua River［J］. Water Resources and Power，2011，29(4):14-17.(in Chinese))

［3］葉柏生，李翀，楊大慶，等.我國過去50a 來降水變化趨勢及其對水資源的影響(Ⅱ):月系列［J］. 冰川凍土，2005，27(1):100-105. (YE Baisheng，LI Chong，YANG Daqing，et al. Variation Trend of precipitation and its impact on water resources in China during last 50 years(Ⅱ):monthly variation［J］. Journal of Glaciology and Geocryology，2005，27(1):100-105.(in Chinese))

［4］葉柏生，陳鵬，丁永建，等.100 多年來東亞地區(qū)主要河流徑流變化［J］.冰川凍土，2008，30(4):556-561.(YE Baisheng，CHEN Peng，DING Yongjian，et al. Discharge changes of the eight large rivers in East Asia during the last more than 100 years［J］. Journal of Glaciology and Geocryology，2008，30(4):556-561.(in Chinese))

［5］牟金玲，狄娟.黑龍江干流中上游段水文特性［J］.黑龍江水利科技，2007，35(2):101-102. (MOU Jinling，DI Juan. Hydrological characteristics in middle and upper reaches of the Heilongjiang River［J］. Heilongjiang Science and Technology of Water Conservancy，2007，35(2):101-102.(in Chinese))

［6］肖迪芳，張雪峰.黑龍江流域水文水資源特性初析［J］.水文，1992(1):51-55.(XIAO Difang，ZHANG Xuefeng.An analysis of hydrology and water resources in Heilongjiang Basin［J］. Journal of China Hydrology，1992(1):51-55.(in Chinese))

［7］吳瓊，梁桂彥，吳玉影，等.黑龍江省四季劃分及氣候特點(diǎn)分析［J］. 林業(yè)勘查設(shè)計，2009(4):95-96. (WU Qiong，LIANG Guiyan，WU Yuying，et al. Four seasons division and climate characteristic analysis of Heilongjiang Province［J］. Forest Investigation Design，2009(4):95-96.(in Chinese))

［8］Mana H B. Non-parametric test against trend［J］.Econometrika，1945，13:245-259.

［9］蔡明科，粟曉玲.黃土高原水文干旱特征分析［J］.西北水資源與水工程，2002，13(4):18-23.(CAI Mingke，SU Xiaoling.Analysis on hydrologic-drought characteristics for loess plateau［J］. Northwest Water Resources ＆ Water Engineering，2002，13(4):18-23.(in Chinese))

［10］WANG Jingcai，XIA Ziqiang，GUO Lidan，et al.Recognizing and forecasting the hydrologic drought in the upper of Weihe basin［C］//IEEE Computer Society.Proceedings of 2011 International Conference on Electrical and Control Engineering，Yichang:IEEE Computer Society，2011:1677-1680.

［11］穆興民，李靖，王飛，等.黃河天然徑流量年際變化過程分析［J］. 干旱區(qū)資源與環(huán)境，2003，17(2):1-5. (MU Xingmin，LI Jing，WANG Fei，et al.Analysis on the annual natural runoff variety process of the Yellow River［J］.Journal of Arid Land Resources ＆ Environment，2003，17(2):1-5.(in Chinese))

［12］鄭紅星，劉昌明.黃河源區(qū)徑流年內(nèi)分配變化規(guī)律分析［J］.地理科學(xué)進(jìn)展，2003，22(6):585-590. (ZHENG Hongxin，LIU Changming. Changes of annual runoff distribution in the headwater of the Yellow River Basin［J］. Progress in Geography，2003，22(6):585-590. (in Chinese))

［13］馮國章，李瑛，李沛成.河川徑流年內(nèi)分配不均勻性的量化研究［J］. 西北農(nóng)林科技大學(xué)學(xué)報:自然科學(xué)版，2000，28 (2):50-53. (FENG Guozhang，LI Ying，LI Peicheng. Quantification of nonuniformity in annual distribution of stream flows［J］. The Journal of Northwest Agricultural University:Natural Sciences，2000，28(2):50-53.(in Chinese))

［14］燕華云，楊貴林，汪青春.長江源區(qū)徑流年內(nèi)分配時程變化規(guī)律分析［J］. 冰川凍土，2006，28(4):526-529.(YAN Huayun，YANG Guilin，WANG Qingchun. Change of annual runoff distribution in the headwaters of the Yangtze River［J］.Journal of Glaciology and Geocryology，2006，28(4):526-529.(in Chinese))

［15］劉賢趙，李嘉竹，宿慶，等.基于集中度與集中期的徑流年內(nèi)分配研究［J］. 地理科學(xué)，2007，27(6):791-795.(LIU Xianzhao，LI Jiazhu，SU Qing，et al. Interannual runoff distribution based on degree and time of concentration for rivers［J］. Scientia Geographica Sinica，2007，27(6):791-795.(in Chinese))