武水流域水文資料“三性”及徑流特性趨勢分析

張英智

（湘西水文水資源勘測中心，湖南 吉首 416000）

1 基本情況

1.1 河流概況

武水上游稱峒河（河溪以上干流稱峒河，于河溪納入沱江后稱武水），是沅江中游一級支流，流域在東經109°20′～109°55′，北緯28°01′～28°29′之間，發源于花垣縣雅酉鎮坡腳村老人山，流經禾排、大興寨、矮寨、寨陽、吉首市區、河溪、潭溪、洗溪，于瀘溪縣武溪鎮城北社區匯入沅江，左岸主要支流有小龍洞河、德夯溪、洽比河、文溪河(分為浪頭河和龍舞河)、司馬河、丹青河等；右岸主要支流有萬溶江、沱江、能溪河等，河流總長150 km，干流平均坡降1.95‰，流域面積為3 691 km2（其中湖南省占3 614 km2）。流域內植被良好，流域內最高點為花垣蓮花山，海拔高1 179 m。武水流域水系情況見圖1。

圖1 武水流域水系及水文站網分布圖

1.2 地形地貌

流域地貌屬中低山峽谷地貌，地勢西南高、東北低，山勢起伏險峻，群峰連綿，若干條峽谷縱橫其間，兩岸山峰聳立，溪河沿山環繞，形成溝壑萬千、山高谷深，地面被切割成大小不等的小盆地、臺地、高峰、溝谷和陡坡、相互交錯、千姿百態的奇特地勢。

1.3 徑流洪水

流域地處亞熱帶季風濕潤區，產流方式以蓄滿產流為主；多年平均降水量1 400 mm左右，多年平均降雨日數在160天左右；多年平均徑流深787 mm。

洪水主要源于暴雨，暴雨和洪水在時空分布上基本一致，大暴雨洪水的主要天氣系統高空是低渦、切變線，地面是冷鋒或靜止鋒。根據吉首、河溪站50多年實測資料統計，洪水一般發生在4—9月，尤以5—8月較多，武水屬典型的山溪性河流，河道坡降大，洪水陡漲陡落，一次洪水過程多為2～5天，單峰多于復峰。

1.4 站網布置

武水流域干流主要布設有吉首、鳳凰、河溪3個國家基本水文站；有雅酉等13個國家基本雨量站。

鳳凰水文站位于鳳凰縣沱江鎮沙灣社區，設立于1962年1月，1966年1月上遷1 100 m改為鳳凰（二）水文站，是武水一級支流沱江中游的重要控制站。現測驗項目有水位、流量、降水量。

吉首水文站位于吉首市峒河街道新橋路，設立于1953年8月，屬省級重要水文站，為武水上游峒河區域代表站，處在沅江水系一級支流武水上游。測驗項目有水位、流量、降水量。

河溪水文站位于吉首市河溪鎮河溪村，設立于1956年11月，1957年1月開始觀測，屬于沅江一級支流武水控制站，在峒河與沱江匯合口下游1.3 km處，現測驗項目有水位、流量、泥沙、降水量和蒸發量。武水流域站網情況見圖1和表1。

表1 武水流域國家基本水文站情況表

2 流域水文資料“三性”分析

水文資料是水文分析計算的依據和基礎，資料的質量好壞關乎水文信息預報、水文分析計算、水利工程設計的合理性，因此，在使用水文資料之前，必須進行水文資料的可靠性、一致性、代表性分析。本文重點對徑流資料的“三性”進行分析。

2.1 可靠性分析

水文資料的可靠性，其實就是指水文資料的正確與否；對于徑流資料來說主要應從徑流資料的來源、測驗方法、整編精度和水量平衡等方面進行檢查。主要檢查洪水沿河長演進是否相應，洪峰流量沿河長變化及其發生時間是否相應合理，洪水總量是否平衡；河槽蓄水量與出水量是否大致相等等。檢查上下游日平均流量變化是否相應；檢查上下游水量是否平衡、徑流模數是否合理，以及降水徑流關系的檢查與對照。

經檢查分析，本次采用的武水流域干支流各代表站資料，均來源于水文機構經過整理、校核、復核、整編、復審、驗收并統一進行刊印后的《水文年鑒》成果；其資料均進行了單站和綜合合理性檢查；本次分析檢查，鳳凰、吉首、河溪站的徑流模數分別為0.022 m3/s·km2、0.031 m3/s·km2、0.026 m3/s·km2，由于各站控制斷面以上多年平均面降雨、下墊面、產匯流條件等差異，徑流模數相差較大（如吉首站以上多年平均面降雨及流域坡降較大，致使徑流模數較鳳凰站偏大40.9%）；河溪站多年平均年徑流量大于吉首站與鳳凰站的多年平均年徑流量之和，即上下游水量平衡；檢查年最大洪水（取典型年2014年7月14—17日洪水）過程，洪峰流量分別是鳳凰站984 m3/s（7月15日）、吉首站1 500 m3/s（7月16日）、河溪站5 060 m3/s（7月16日），洪水沿河長演進、洪峰流量沿河長變化及其發生時間是相應且合理的，且洪水總量平衡；各水文站多年平均流量與集水面積關系合理（見圖2）；流域各站歷年斷面變化較小，水位流量關系較穩定；因此，武水流域徑流資料可靠性較高。

圖2 武水各站年平均流量與集雨面積關系圖

2.2 一致性分析

水文資料的一致性分析，是水文水資源的基礎性工作內容。對于某些資料，從年際來看，一是由于受觀測、測驗方法的改變導致所觀測量的相應改變；二是因為在系列時間范圍內有人類活動影響（如歷史上測站斷面位置發生遷移、存在有跨流域引水和上游建有影響較大的水利工程等），也就是說表征產流條件下墊面發生變化（或下墊面條件改變不在容許范圍內），且變化影響很大時，系列的一致性受到影響。這些變化通常在水文刊印的資料中難以發現，因此，在使用水文資料時，必須檢查歷年資料的一致性。對存在明顯的某些不合理的變化成分，應進行物理成因分析和統計分析，查明原因，加以排除；對有人類活動影響的情況，要根據情況進行還原或修正。對于徑流量資料一致性分析，一般采用有序聚類法、降雨徑流相關法進行檢查分析。

有序聚類法，就是以有序分類來推估最可能的突變點，實質是尋求最優分割點，使同類間的離差平方和較小而類與類之間的離差平方和較大。一般地，若序列有兩個明顯的階段性過程，則總離差平方和的時序變化呈現單谷底現象；若有幾個明顯的階段性過程，則總離差平方和的時序變化會有幾個谷底，這樣可以根據谷底發生的時間劃分序列變化的階段。本次分析，將武水流域鳳凰、吉首、河溪3個站1957—2020年水文序列采用有序聚類法分析流域年徑流序列發生變化的突變點，詳見圖3、圖4、圖5。由圖3、圖4、圖5可知，3個站年徑流離差平方和總體趨于平穩狀態，雖然在2012年吉首、河溪站年徑流離差平方和達到最小值，但兩站徑流序列沒有出現顯著變化。

圖3 鳳凰站年徑流量離差平方和變化過程

圖4 吉首站年徑流量離差平方和變化過程

圖5 河溪站年徑流量離差平方和變化過程

由此可知，武水流域雖建有長潭崗水庫等水利工程，但經分析，水利工程建設對武水年徑流系列一致性的影響不大。

降雨徑流相關法檢查。本文將1957—2020年鳳凰、吉首、河溪站年降雨量與年徑流量進行相關分析，如圖6～圖8。鳳凰站降雨徑流關系點據稍顯散亂，相關系數0.687 3；吉首站、河溪站降雨徑流關系點據稍好，吉首站降雨徑流相關系數0.861 2；河溪站降雨徑流相關系數0.7993；3站降雨徑流關系無明顯系統偏離；關系點據的偏離大小呈隨機變化，無明顯趨勢性。

圖6 鳳凰站降雨徑流相關圖

圖7 吉首站降雨徑流相關圖

圖8 河溪站降雨徑流相關圖

綜上所述，武水流域徑流資料具有一致性，無需進行一致性修正。

2.3 代表性分析

水文資料代表性，就是對資料的穩定性、周期性和變化規律進行分析，揭示系列的代表程度。本文選取武水控制站河溪水文站1957—2020年徑流系列進行代表性分析；年徑流多年變化周期性和穩定性，可以通過年徑流差積、累積流量過程進行分析，如圖9。差積曲線隨時間出現多個峰谷交替，也就是說自1957年來河溪站存在多個豐、枯水年的交替循環現象。如圖10，隨著時間的增長，河溪站年徑流模比系數累積曲線趨近于1，因此，武水河溪站徑流資料代表性較好。

圖9 河溪站年徑流模比系數差積曲線圖

圖10 河溪站年徑流模比系數累積曲線圖

由以上分析表明，徑流系列中，豐、平、枯水分布較勻，基本上能反映徑流系列總體的特性，滿足規范要求，具有代表性。

3 徑流特征趨勢分析

3.1 徑流年際變化

從圖9可以看出，1958—2020年，河溪水文站年徑流量經歷了幾個大的豐枯循環：1958—1966年，徑流量呈遞減趨勢；1967—1977年，徑流量總體呈遞增趨勢；1978—1988年，徑流量總體呈遞減趨勢；1994—1999年，徑流量呈遞增趨勢；2000—2011年，徑流量呈遞減趨勢；2012—2020年，徑流量呈遞增趨勢。點繪5年年徑流滑動平均統計分析圖（圖11），由圖11可知，河溪水文站年徑流量的變化趨勢與圖9的變化趨勢是一致的，雖有豐枯交替，但總體上呈微弱上升的趨勢。

圖11 河溪站年徑流量5年滑動平均曲線圖

3.2 Kendall秩次相關檢驗

采用Kendall法，對河溪水文站年徑流序列秩次相關檢驗，當n=64、置信度水平取α=0.05時，Uα/2=1.96。本文對河溪水文站1957—2020年年徑流量序列進行Kendall法檢驗統計計算，其結果U=0.85，|U|＜Uα/2表示序列變化趨勢不顯著；且U＞0表示序列呈上升趨勢。也就是說河溪水文站年徑流量序列變化趨勢呈現出并不顯著的上升趨勢。

3.3 R/S法分析

根據該流域河溪水文站多年徑流資料的分析，將年徑流量時間序列:{X i}={1，2，3，……64}，應用R/S分析法對河溪水文站年徑流量時序進行分析。對所給序列進行分析，并進行直線擬合，得到年徑流量時間序列（圖12）的回歸方程為ln（R/S）=0.695 9ln i-0.325 6，計算結果H=0.695 9＞0.5，這表明河溪水文站年徑流量的變化具有長程相關性，存在正的持續性，意味著未來的總體趨勢將與過去相同，即年徑流量在未來一段時期有持續增加的趨勢。

圖12 河溪水文站年徑流量時間序列R/S分析結果

4 結 語

本文采用有序聚類法等方法對流域水文站徑流資料“三性”進行分析，并利用Kendall趨勢檢驗法、R/S分析法的理論對代表站河溪站年徑流量時序變化的趨勢特性及豐枯演化規律進行研究，得出如下結論:

1）經采用有序聚類法、降雨徑流相關法等，對徑流資料“三性”分析得出，武水徑流資料具有較高的可靠性、明顯的一致性和較好的代表性。

2）由Kendall趨勢檢驗法分析可知，武水控制站河溪水文站年徑流量呈增加趨勢，但并不顯著。

3）由R/S分析法分析，結果表明武水控制站河溪水文站年徑流量變化存在正的持續性，從平均的觀點而言，可以預測河溪水文站的年徑流量將繼續呈微弱增加的趨勢。