韓江干流潮安站泥沙變化分析

張郁坡

摘 要：研究了1951—2011年韓江干流潮安站泥沙的基本變化特征。結果表明，潮安站年輸沙量年內分配不均，月際變化明顯，輸沙量主要集中在4—9月，且前汛期大于后汛期。1994年以后，由于上游多個梯級電站、水庫和樞紐的建成，年平均含沙量及年輸沙量呈現明顯減少的趨勢。M-K突變檢驗顯示，2001年以后輸沙量快速減少。

關鍵詞：輸沙量；汛期；河床；水資源

中圖分類號：P333 文獻標識碼：A DOI：10.15913/j.cnki.kjycx.2016.15.078

1 流域概況

韓江流域位于粵東、閩西南地區，地理位置在東經115°13′～117°09′，北緯23°17′～26°05′，是廣東省僅次于珠江流域的第二大流域。該流域包括廣東、福建和江西三省，干流全長470 km（自河源至東溪口），落差920 m，平均坡度為0.4‰，流域面積為30 112 km2。梅江是韓江的主流，發源于廣東省紫金縣的七星崠，在三河壩與汀江匯合后，經潮州市進入韓江三角洲，分北、東、西溪出南海。

潮安水文站是國家重點基本水文站，其控制流域面積為29 077 km2，占總流域面積的96.6%，是韓江干流的控制站。本文通過1955—2011年潮安站的輸沙量資料，分析了韓江輸沙量的年內、年際變化，可以為該流域河床變化、河沙開采及水資源的發展與合理利用提供一定的參考。

2 泥沙的年內變化

2.1 分析方法

2.1.1 不均勻系數

不均勻系數Cu（Unevenness Coefficient）又稱變差系數，是反應時間序列分配不均勻程度的一個綜合性指標，能用于多種性質的時間序列的分析研究，因其具有良好的時間空間可比性，在我國得到了廣泛應用。本文用來分析輸沙量年內分布的不均勻性，其計算公式為：

式（3）中：Wsi 為第i個月輸沙量；WS為全年輸沙量。CI取值介于8.33～100.當CI=8.33時，年輸沙量均勻地分配在12個月內；當8.3320時，輸沙量月際變化明顯；當CI=100時，年輸沙量全部集中在一個月內。

2.2 潮安站泥沙年內變化情況分析

2.2.1 不均勻性分析和集中指數分析

通過式（1）（3）分別計算逐年Cu和CI值，并點繪逐年變化過程線，如果圖1所示。從圖1可以看出，潮安站輸沙量年內分布的Cu值基本上都大于1，平均值為1.38，最大值為2.43（2007）；CI值均大于10，大部分年份大于20，平均值為25.1，最大值為57.5（2007）。由此可見，輸沙量月際變化明顯。

計算不同時間段輸沙量年內分配系數Cu，并點繪逐時段變化過程線，如圖2所示。由圖2可知，1975—1984年的Cu值是最小的，之后逐漸變大。由此可見，潮安站輸沙量年內分配不均勻現象越來越明顯。

2.2.2 泥沙月分配分析

圖3為潮安站各月輸沙量和徑流量的分布情況。從圖3中可以看出，與徑流的月分布相同，潮安站的輸沙量主要集中在4—9月，占全年輸沙量的87%以上，表現出該站泥沙年內分布集中、極不均衡的特點。但相比徑流量的分布，沙量的分布更加不均勻，且出現了前汛期沙量大于后汛期的情況。

給定顯著性水平a，查正態分布表得到臨界值Ua，如果|UFk|>Ua，則表明序列存在明顯的趨勢變化。將所有UFk組成一條曲線C1，通過信度檢驗可知其是否具有趨勢。同理，將此方法引用到反序列，將反序列Xn，Xn-1，…，X1表示為X1，X2，…，Xn，重復上述計算過程，并使計算值乘以-1，得到UBk，其在圖中表示為C2。分析繪出的UF和UB曲線圖，如果UF或UB>0，則表明序列呈上升趨勢；如果UF或UB<0，則表明序列呈下降趨勢。當其超過臨界值線（信度線）時，表明上升或下降的趨勢明顯。超過臨界線的范圍確定為出現突變的時間區域。如果UF和UB出現兩線相交，且交點在臨界區域之內，則交點對應的時刻便是突變開始的時間。

3.2 潮安站泥沙年際變化情況分析

3.2.1 年際變化趨勢分析

用M-K趨勢分析了1955—2011年泥沙序列。趨勢檢驗結果表明，潮安站年平均含沙量和年輸沙量的Z值分別為-5.14和-3.91，Z的絕對值均大于2.32，表示通過了置信度99%的顯著性檢驗。由此可見，潮安站年平均含沙量和年輸沙量減少的趨勢十分明顯。具體如表1所示。

繪制了潮安站年平均含沙量差積累計曲線，如圖4所示。在圖4中可以看出，1987年以前呈上升趨勢，1994年以后呈下降趨勢，原因為90年代以后，韓江中上游多個梯級電站、水庫和樞紐建成，使上游來沙在區間沉積，大量泥沙被攔截，因此，韓江中下游含沙量及年輸沙量明顯減少。

3.2.2 年際變化突變檢測

利用Mann-Kendall法對潮安站的年輸沙量數據進行了突變檢驗，趨勢檢驗結果如圖5所示。由圖5可見，潮安站的年輸沙量從1994年以來有明顯的減少趨勢，20世紀末，這種減少的趨勢通過了a=0.05的顯著性檢驗，甚至通過了a=0.001的檢驗。這表明潮安站的年輸沙量的減少趨勢是十分明顯的，再通過UFk和UBk兩條曲線的交點位置，可以確定年輸沙量的減少屬于突變現象，具體是從2001年開始的。其原因為2000年以后潮州供水樞紐的建設改變了水的流速，進而影響了含沙量和年輸沙量。

3.2.3 輸沙量徑流量累積分析

水沙累積關系線表示的是徑流量與輸沙量的關系。如果流域的水沙關系發生系統變化，則水沙累積關系線上將出現明顯的轉折，即累積曲線斜率發生變化。潮安站徑流量輸沙量累積關系線點如圖6所示。

從圖6可以發現，潮安站在1993年以前的累積關系線基本呈一條直線，表明輸沙量與徑流量的關系未發生系統變化；1994年以后，累積線偏向徑流量軸，表明1994年以后潮安站年輸沙量相比于徑流量出現了系統減少。

4 結論

河流泥沙的產生是降水徑流、自然地理因素和人類活動等綜合作用的結果，泥沙的變化過程也是極其復雜的。通過本文的分析，韓江干流潮安站的年輸沙量月際變化顯明顯，輸沙量主要集中在4—9月，且前汛期大于后汛期。1994年以前，輸沙量與徑流量的關系是相對穩定的；1994年以后，年輸沙量相比于徑流量出現了系統減少，這是因為90年代以后，上游多個梯級電站、水庫和樞紐的建成攔截了部分泥沙。此外，上游河段的水土保持工作減少了泥沙的沖刷量，使潮安站年平均含沙量及年輸沙量呈現明顯減少的趨勢。M-K突變檢驗顯示，2001年以后年輸沙量快速減少。

參考文獻

[1]張信寶，文安邦.長江上游干流和支流河流泥沙近期變化及其原因[J].水利學報，2004（04）.

[2]聶文晶.洮河流域沙特性分析[J].甘肅水利水電技術，2012（01）.

[3]杜文印.珠江三角洲近30年水沙分配變化分析[J].廣東水利水電，2012（06）.

[4]陸建宇，王秀慶.徑流年內分配不均勻性的度量指標及其應用[J].水利發電，2015（11）.

〔編輯：張思楠〕