黃河上游蘭州水文站2004—2018年水沙特性分析

馬 雁，賈生海，張彥洪

(甘肅農業大學水利水電工程學院，甘肅 蘭州 730070)

黃河作為西北地區重要的淡水資源，在地區社會經濟發展中扮演著重要的角色。黃河泥沙問題由來已久，相關學者為此做了大量調查研究，其中趙陽等研究了黃河流域近70年水沙現狀及成因分析[1]，姚文藝等通過黃河上游進出水庫的水沙資料分析了水庫運行對徑流泥沙的調節及影響[2]，彭少明等利用IHA(Indicators of Hydrologic Alteration)指標體系分析了梯級水庫群不同運行方式對黃河上游水文變化和對河流生態的影響[3- 4]，同時人類活動和氣候變化在一定程度上影響著流域水沙運移規律[5- 7]。

黃河上游通常是指內蒙古托克托縣河口鎮以上河段，全長3472km，占總長的63.5%，流域面積約為38.6萬km2，占全河總流域面積的51.3%左右。黃河上游天然徑流量的99%來自蘭州以上地區，泥沙主要來源區間位于蘭州至青銅峽地段[8]。黃河上游水體泥沙含量小，水能資源蘊藏量大，沿途建造的水利樞紐工程大大推動了區域社會經濟和地方水利電力行業的發展。但因黃河含沙量大，人類活動及水利工程的修建，改變了河道的水沙沖淤平衡，河岸面積縮減、部分寬谷河道主河槽淤積嚴重，河道健康和生態功能進一步減退，調控洪水的能力減弱，可能導致洪水淹沒及人民財產的損失[9- 11]。因此，研究黃河上游水沙特性變化規律，對水資源合理配置、維持地方經濟社會發展，保障人民財產安全具有一定的現實意義。

1 黃河上游泥沙來源分析

黃河主要支流有湟水、白河、黑河、洮河、大夏河、大黑河、祖厲河、清水河等，多年平均天然徑流量達330億m3，占黃河天然徑流總量的57.2%。黃河中的泥沙主要來源于支流泥沙的匯入，上游典型支流為洮河，控制水文站為紅旗站[12]。黃河上游部分河段經過黃土高原地區，該地區汛期降雨集中、降水強度大，而植被稀疏矮小，結構疏松的表土在雨水的強力沖刷下，極易引發水土流失；內蒙古河口鎮以上黃河段流經騰格里沙漠和毛烏素沙漠，地表沙土在水流侵蝕及風蝕等條件下每年將有大量的泥沙進入黃河。因此，降雨侵蝕和水土流失是黃河上游泥沙的主要來源，這一觀點同時也被相關學者證實[13]。近年來，土地利用類型及城市化等人類活動在一定程度上改變了黃河上游流域下墊面狀況，使流域水循環過程發生改變，對流域泥沙的產生也有較大影響[14]。

2 蘭州水文站不同年份水沙特征分析

蘭州水文控制站位于黃河干流之上，其控制流域面積達22.26萬km2。該站2004—2018年年徑流量、年輸沙量及年平均含沙量來源于水利部《中國河流泥沙公報》發布數據。

2.1 蘭州水文站不同年份年徑流量變化特征

2004—2018年蘭州水文控制站15年間年徑流量變化特征如圖1所示。由圖1可知，2004—2018年，年徑流量在235.6億～441.8億m3之間變化，平均徑流量為302.8億m3，年徑流量呈先增大后減小再增大的變化趨勢。2012年是15年來第一個年徑流量達到高峰的年份，為380.4億m3，而第二次年徑流量高峰出現在2018年，達到441.8億m3，為2004年來的最大徑流量，較其他年份平均流量增加50.85%。經統計2012年該水文站所在流域降雨、徑流和2018年黃河流域洪水等資料發現[15—17]，2012年黃河上游產生了一次較大規模的洪水過程，7月30日10:18洪峰流量達到3860m3/s，是1986年以來的最大洪峰流量；這次洪水過程中日均流量大于2500m3/s的天數為29d，日均流量超過3000m3/s的天數為9d，因此導致2012年年徑流量高于除2018年的其他年份。2018年汛期黃河流域暴雨過程頻繁，共出現3次編號洪水，中雨以上降水過程共21次，其中在第2號洪水過程中，蘭州水文控制站7月23日8:18洪峰流量達3610m3/s，此次洪水量達到3.05億m3；在第3號洪水過程中，9月26日23:00洪峰流量達3590m3/s，總水量為191億m3，洪水歷時87d，本次日均流量大于3000m3/s的天數達30d，3次洪水過程導致2018年年徑流量明顯大于其他年份。

2.2 蘭州水文站不同年份年輸沙量變化特征

2004—2018年蘭州水文控制站15年間年輸沙量的變化特征如圖2所示。由圖2可知，2004—2018年，年輸沙量在0.079～0.960億t之間變化，其中2004—2007年年輸沙量平均值為0.257億t，呈逐年遞增趨勢；2008—2011年年輸沙量平均值為0.109億t，處于較低水平，2004—2007年年輸沙量是2008—2011年年輸沙量的2.45倍；2012年由于黃河上游一次較大的洪水過程，輸沙量迅速增大到0.372億t；2013—2017年年輸沙量平均值為0.118億t，與2008—2011年均值接近；2018年輸沙量達到最大值0.960億t，是其他年份平均輸沙量的5.54倍。通過查閱和分析水文站相關資料可知[18]，2018年年徑流量較其他年份偏大，為441.8億m3，同時2018年汛期受降雨影響，黃河上游洪水多、水量大、降雨持續時間長，伴隨洪水過程，受湟水、大通河等支流的影響，黃河上游來水含沙量增大，輸沙量明顯增多。

圖2 2004—2018年年輸沙量變化特征

2.3 蘭州水文站不同年份年平均含沙量變化特征

2004—2018年年平均含沙量的變化特征如圖3所示。由圖3可知，年平均含沙量最小值為0.260kg/m3，最大值為2.170kg/m3，對應的年份分別為2009、2018年。2004—2007年年平均含沙量逐年遞增，平均值為0.898kg/m3，與2012年0.978kg/m3平均含沙量接近；2018年平均含沙量為2.170kg/m3，是2004—2007年平均值的2.42倍；其余年份平均值為0.403kg/m3，是2018年的18.57%，2004—2007年平均值的44.88%?？傮w來看，在統計數據中除2018年外，其他年份年平均含沙量均值為0.585kg/m3，變化幅度較小，2018年較其他年份平均增加了270.67%，這種趨勢與該年度徑流量和輸沙量較其他年份大相一致。同時，由于2018年汛期暴雨洪水增多且歷時長，促使其平均含沙量較其他年份明顯偏多。

圖3 2004—2018年年平均含沙量變化特征

2.4 蘭州水文站不同年份年徑流量與年平均含沙量關系

河道中年徑流量大小與降水量及支流來水等因素有關。2004—2018年以來年徑流量與年平均含沙量之間的變化關系如圖4所示。由圖4可以得出，在不同年份隨徑流量的增加，年平均含沙量呈逐漸增大的趨勢，可用線性函數表示。

f(x)=-1.2091+0.0063x

(1)

式中，f(x)—平均含沙量，kg/m3；x—年徑流量，億m3；R—相關系數，這里取R=0.6831。由此可見，在平均含沙量變化不大的情況下，年徑流量越大，年平均含沙量也就越大。

圖4 2004—2018年年徑流量與年平均含沙量變化關系

2.5 蘭州水文站不同年份年輸沙量與年平均含沙量關系

2004—2018年以來年輸沙量與年平均含沙量之間的變化關系如圖5所示。由圖5可以得出，年輸沙量與年平均含沙量之間呈線性關系，可用線性函數表示。

f(x)=0.2078+2.1409x

(2)

式中，f(x)—平均含沙量，kg/m3；x—年輸沙量，億t，R—相關系數，這里取R=0.9770。由此可見，年輸沙量與年平均含沙量二者相關性強，隨著年輸沙量的增大，年平均含沙量也隨之增大。

圖5 2004—2018年年輸沙量與年平均含沙量變化關系

2.6 蘭州水文站不同年份年徑流量與年輸沙量關系

黃河上游水沙異源，不同時期水沙變化主導因子不同，20世紀60年代以前以降雨為主導因子，之后隨著上游大型水庫的修建及水保措施的實施，水沙變化受自然因素和人類活動共同影響[17]。2004—2018年年徑流量與年輸沙量變化關系如圖6所示，由圖6可以得出，在不同年份隨徑流量的增加，年輸沙量呈逐漸增大的趨勢，可用線性函數表示。

f(x)=-0.7849+0.0033x

(3)

式中，f(x)—年輸沙量，億t；x—年徑流量，億m3；R—相關系數，這里取R=0.7962。

由相關系數可知年徑流量與年輸沙量相關性較差，說明水沙變化受自然因素和人類活動的共同影響較大，同時水沙變化并不完全遵循同步規律，存在時間上的差異性。

圖6 2004—2018年水沙變化關系

3 結語

降雨和洪水對黃河上游徑流量影響較大，同時支流洪水不僅導致洪峰的出現，還攜帶大量泥沙。通過分析黃河上游蘭州水文控制站水沙來源及徑流、泥沙特征，研究認為2004—2018年，年徑流量呈先增大后減小再增大的趨勢，輸沙量和平均含沙量隨年份變化不明顯，2018年受暴雨洪水影響，三項指標均高于其他年份；在不同年份隨徑流量和輸沙量的增加，年平均含沙量均呈線性增長的趨勢，水沙變化并不完全遵循同步規律，存在時間上的差異性。因此，加強洪水預報，減少支流泥沙來源，降低洪水過程中攜帶的泥沙量，有助于減少黃河上游泥沙含量，減少下游泥沙淤積量。