龍羊峽水庫調度對徑流的影響及蓄補水規律

尚文繡，許明一，尚弈，嚴登明

(1.黃河勘測規劃設計研究院有限公司，鄭州 450003；2.華北水利水電大學，鄭州 450045)

徑流是河流的關鍵組成要素和河流生態服務功能的重要影響因素[1-3]。河流天然徑流過程一般隨時間表現出明顯的豐枯變化并具有周期性特征[4]，這些變化特征發揮著傳輸物質、傳遞信息、塑造棲息地、維持生物生活史模式等作用[5-6]。受到人類活動影響，全球河長1 000 km以上的河流中約2/3處于非自然流動狀態[7]，很多河流的徑流過程已經發生了顯著變化[8-9]。

水庫通過調蓄河川徑流實現供水、防洪、發電等功能，是改變河流天然徑流過程的主要原因之一[10-12]。隨著人類社會對供水安全保障和電力的需求不斷增長，水庫數量快速增長[13]，世界上已建在建壩高超過30 m的大壩已超過1.5萬座[14]。2019年我國已建成大中型水庫4 722 座，總庫容超過8 000 億m3。水庫建設運行直接改變了河流徑流過程，進而改變物質輸移、生物遷徙、水文循環等過程，對河流生態健康造成復雜影響[15-17]。量化水庫對徑流過程的影響是識別水庫生態影響、提升水庫生態保護功能的基礎。當前研究[9,18-20]中，評價水庫對徑流過程影響的主要方法是量化水庫運行前和運行后兩個時段的徑流過程關鍵特征，將關鍵特征的差異作為水庫的影響。但是河流水文情勢受到氣候變化、工程調節、取水、退水等多因素的復合影響[21-23]，對比水庫運行前后的徑流過程變化反映的是多因素的共同作用，并不能得到水庫調度這單一因素發揮的作用。

位于黃河上游的龍羊峽水庫是黃河干流上唯一一座多年調節水庫，控制黃河天然徑流量的40%以上，調節庫容193.6 億m3，對黃河徑流具有很強的調節能力[24]。黃河天然徑流量和取用水已經發生明顯變化，天然年徑流量從580 億m3(1919—1975年系列)衰減至490 億m3(1956—2016年系列)，進入21世紀以來黃河取水量增長80 億m3以上，說明黃河徑流過程變化受到多因素影響。2000—2016年龍羊峽水庫入庫水文站實測年徑流量比20世紀60—90年代減少了16.24%，導致難以通過對比水庫運行前、后兩個時段徑流過程變化來準確量化龍羊峽水庫的作用。本文對比龍羊峽水庫運行后同一時段內入庫徑流過程和出庫徑流過程關鍵特征的差異，結合突變檢驗，量化龍羊峽水庫在黃河徑流過程變化中發揮的作用，并通過顯著性檢驗識別龍羊峽水庫蓄補水規律。

1 研究方法與數據

1.1 研究區域

研究對象為位于青海省共和縣、貴南縣交界處的龍羊峽水庫。水庫壩址控制黃河流域面積13.1萬km2，占全流域面積的17.4%。龍羊峽水庫1978年開工建設，1986年開始蓄水。水庫的開發任務以發電為主，兼有防洪、灌溉、防凌、養殖、旅游等綜合效益。

通過對比同一時段內龍羊峽水庫出庫徑流過程相對于入庫徑流過程的變化，可以揭示水庫調度這一因素對河流徑流過程的影響。選用的水文站監測斷面和評價指標見表1。1986—2016年唐乃亥斷面天然年徑流量186.53億m3，貴德斷面天然年徑流量187.71億m3，兩斷面天然年徑流量僅相差0.64%，說明兩斷面間的區間來水對徑流影響很小。

表1 情景設置Tab.1 Scenario setting

1.2 水文情勢變化評價方法

采用水文變化指標(indicators of hydrologic alteration，IHA)[20]評價龍羊峽水庫出庫斷面水文情勢相對于入庫斷面水文情勢的變化情況，評價時段為1986—2016年。IHA指標體系包含5組33個評價指標(表2)，分別描述了流量、頻率、發生時機、持續時間和變化率等水文情勢關鍵特征。計算1986—2016年龍羊峽水庫出庫斷面各IHA指標值相對于入庫斷面相應IHA指標值的變化幅度(即偏離度)。

表2 IHA評價指標體系Tab.2 Indicators of Hydrologic Alteration

采用變化范圍法(range of variability approach，RVA)評價水文情勢變化程度[25]。將近天然時段各IHA指標值的75%分位數至25%分位數作為RVA目標區間[19]，水文改變度為

(1)

式中：Di是第i個IHA指標的水文改變度；No是評價時段內落到RVA目標區間的年數；Ne是評價時段內落到RVA目標區間的預期年數，根據本文選擇的RVA目標區間，預期年數為評價時段總年數的一半；i是IHA指標編號，i=1,2,…,33。水文改變度分為3個等級：Di<33%時代表無改變或低度改變；33%≤Di≤67%時代表中度改變；Di>67%時代表高度改變。

基于IHA-RVA方法，Shiau等[26]提出了整體改變度為

(2)

式中：DA是整體改變度。

1.3 流量過程變異程度評價方法

通過流量過程變異程度表征1980—2016年龍羊峽水庫入庫斷面和出庫斷面實測月徑流相對于天然月徑流的變化程度

(3)

表3 流量過程變異程度賦分標準Tab.3 Flow deviation indicator scoring standard table

1.4 顯著性檢驗方法

采用t檢驗判斷龍羊峽水庫蓄補水量與黃河天然徑流量、黃河流域地表水耗水量占天然徑流量的比例之間是否存在顯著的線性關系。

提出原假設H0：β=0。β是因變量y與自變量x之間的線性回歸系數，β為0代表y與x之間沒有線性相關關系。構造檢驗統計量

(4)

1.5 突變檢驗方法

采用Mann-Kendall(M-K)檢驗法[28]分析龍羊峽出庫斷面流量過程變異程度的突變點。對于一組序列長度為n的獨立序列yf(f=1,2,…,n)，構造統計量df和Ff

(5)

(6)

式中：當yi>yj時，ri=1，否則ri=0，j=1,2,…,i；E(df)和V(df)分別為df的均值和方差；令F1=0。給定顯著性水平α=0.05，當|Ff|>Uα時，表明序列有顯著的趨勢變化。Uα為概率超過1-α/2時標準正態分布的值。

將yf逆序排列，按照式(5)和式(6)得到一組新的Ff，再令統計量Bf=-Ff。如果Ff和Bf隨f的變化曲線有交點，且交點位于±Uα之間，那么交點為突變發生的時間。

1.6 數據來源

龍羊峽水庫1980—2016年入庫和出庫日徑流數據來自相應水文站提供的監測斷面實測日徑流資料。1980—2016年入庫和出庫天然月徑流數據采用水利部黃河水利委員會的徑流還原成果。黃河流域歷年地表水耗水量來自《黃河水資源公報》。

2 龍羊峽水庫調度對徑流的影響

2.1 對年際和年內徑流過程的影響

1986—2016年，龍羊峽水庫入庫最大年徑流量327.74億m3，最小年徑流量105.77億m3，極值比3.10，年徑流量變異系數0.26。經過水庫調蓄后，出庫最大年徑流量252.88億m3，最小年徑流量108.61億m3，極值比2.33，年徑流量變異系數0.19。這說明龍羊峽水庫調度減小了年徑流量的變化范圍和離散程度。

對于年內徑流過程，1986—2016年龍羊峽水庫入庫與出庫月均流量見圖1。水庫調度后：6—10月的月均流量減小6.11億～14.73億m3；11月—次年5月的月均流量增加，其中11月和5月增幅較小，12月—次年4月的月均流量增加5.09億～8.58億m3。

圖1 1986—2016年龍羊峽水庫入庫與出庫月均流量對比Fig.1 Comparison of monthly average inflow and outflow of Longyangxia Reservoir from 1986 to 2016

2.2 對水文情勢的影響

對比1986—2016年龍羊峽水庫出庫徑流和入庫徑流的IHA指標值，結果見圖2。指標1～12顯示：與入庫徑流相比，出庫徑流6—10月的月均流量減小，降幅18.19%～43.15%；11—5月的月均流量增加，其中12—3月增幅超過100%。指標13～21反映了極端高/低流量事件的變化：與入庫徑流相比，出庫徑流的最大1、3、7、30、90 d平均流量減小34.74%～47.30%；最小1、3、7、30、90 d平均流量均增大，且隨著天數增加，增幅從38.99%增加到144.43%。龍羊峽水庫入庫及出庫徑流均沒有發生過零流量(指標23)。出庫徑流的基流指標比入庫徑流增加80.24%(指標24)。指標25和26反映了與入庫徑流相比，出庫徑流的最大1 d平均流量發生日期略推遲，而最小1 d平均流量發生日期大幅提前。指標27～30反映了高/低流量脈沖的變化:與入庫徑流相比，出庫徑流高流量脈沖發生次數增加，但平均持續時間縮短，總持續時間縮短約2/3；低流量脈沖的發生次數與平均持續時間均減小，總持續時間減小90%以上。指標31～33反映了日流量的變化率，與入庫徑流相比，出庫徑流的平均上升速率和平均下降速率均顯著增大，流量翻轉次數大幅增加。

圖2 1986—2016年龍羊峽水庫出庫徑流相對于入庫徑流的IHA指標值偏離度Fig.2 Deviation degree of IHA value of outflow runoff relative to inflow runoff of Longyangxia Reservoir from 1986 to 2016

與入庫徑流相比，出庫徑流的IHA指標中有6個發生低度改變，分別為5、6、8、11月的月均流量，零流量天數和最大1 d平均流量發生日期；其他27個指標均為高度改變。出庫徑流相對于入庫徑流的整體改變度為95.04%，屬于高度改變。

2.3 對流量過程變異程度的影響

1980—2016年龍羊峽水庫入庫及出庫流量過程變異程度見圖3。龍羊峽水庫入庫徑流的流量過程變異程度始終較低，變化范圍為0.02～0.06，除2002年外流量過程變異程度的賦分均為100。龍羊峽水庫以上河段用水很少，20世紀90年代以來龍羊峽水庫以上河段的地表水耗水量僅1.05億～2.82億m3，不足龍羊峽壩址處天然徑流量的2%，不足黃河流域地表水耗水量的1%，且龍羊峽以上黃河干流水庫少、調蓄能力弱，因此龍羊峽入庫徑流的流量變異程度很小。2002年黃河天然徑流量僅246.16億m3，比多年平均值偏低49.76%，是1956年以來黃河來水最枯的年份，因此，龍羊峽以上河段的人類活動對流量過程的干擾相對明顯，但這一年流量過程變異程度賦分仍高達97.10。

圖3 1980—2016年龍羊峽水庫入庫及出庫流量過程變異程度及賦分Fig.3 Flow deviation indicator value and its score of inflow and outflow of Longyangxia Reservoir from 1986 to 2016

1986年以前龍羊峽水庫出庫流量過程變異程度及賦分與入庫接近，但此后流量過程變異程度迅速大幅增加、賦分迅速大幅降低。1980—1985年龍羊峽水庫出庫流量過程變異程度均值0.05，賦分均值98.29。1987—2016年龍羊峽水庫出庫流量過程變異程度均值2.15，是1980—1985年均值的43倍；賦分均值20.14，比1980—1985年均值減小了79.51%。

M-K檢驗結果顯示，在顯著性水平0.05的情況下，龍羊峽水庫出庫流量過程變異程度的突變點為1986年。突變時間與龍羊峽水庫開始蓄水的年份一致。鑒于1986年前后龍羊峽以上河段取用水沒有發生明顯變化，上游也沒有大型水利工程投入運行，龍羊峽水庫投入運行是導致龍羊峽水庫出庫斷面流量過程變異程度突變的主要原因。

3 龍羊峽水庫蓄補水規律

3.1 水庫蓄補水與來水、耗水的關系

1986—2016年黃河天然徑流量及龍羊峽水庫補水量見圖4。龍羊峽水庫最大年蓄水量95.69 億m3，最大年補水量53.74 億m3，年均蓄水4.61 億m3。龍羊峽水庫補水量與黃河天然徑流量間呈現出反比的變化趨勢，表現為天然徑流量小的年份水庫補水、天然徑流量大的年份水庫蓄水。例如：2005年黃河天然徑流量為579.88億m3，比多年平均值(490.00 億m3，1956—2016年系列，下同)偏高18.34%，龍羊峽水庫年蓄水95.69億m3；2006年黃河天然徑流量407.16 億m3，比多年平均值偏低16.91%，龍羊峽水庫年補水51.83 億m3。對1986—2016年黃河天然徑流量和龍羊峽水庫年補水量進行線性擬合，擬合公式的斜率為-0.25，相關系數為0.55。對相關關系進行顯著性檢驗，t值為-3.52，絕對值高于顯著性水平0.05時的臨界值，說明黃河天然徑流量和龍羊峽水庫年補水量之間線性關系顯著，且為負相關關系。

圖4 1986—2016年黃河天然徑流量及龍羊峽水庫補水量Fig.4 Natural runoff of the Yellow River and replenishment of Longyangxia Reservoir from 1986 to 2016

1990—2016年黃河流域地表水耗水量占天然徑流量的比例見圖5。《黃河水資源公報》缺少1996年和1997年地表水耗水量數據，因此不計算這2 a地表水耗水量占天然徑流量的比例，僅對25 a的數據進行分析。對黃河流域地表水耗水量占天然徑流量的比例和龍羊峽水庫補水量之間的線性關系進行顯著性檢驗，t值為2.91，高于顯著性水平0.05時的臨界值，說明兩者線性關系顯著，且為正相關關系。統計結果顯示，當黃河流域地表水耗水量占天然徑流量的比例超過64%時，龍羊峽水庫以補水為主，否則以蓄水為主。25 a中有13 a龍羊峽水庫補水，其中11 a黃河流域地表水耗水量占天然徑流量的比例超過64%；有12 a龍羊峽水庫蓄水，其中9 a黃河流域地表水耗水量占天然徑流量的比例低于64%。

圖5 1990—2016年黃河流域地表水耗水量占天然徑流量的比例及龍羊峽水庫補水量Fig.5 Proportion of surface water consumption in natural runoff in the Yellow River basin and water replenishment of Longyangxia Reservoir from 1990 to 2016

3.2 水庫的蓄豐補枯作用

在多年間蓄豐補枯是多年調節水庫的重要任務，分析龍羊峽水庫在黃河流域年際蓄豐補枯中發揮的作用。分析時段為1990—2016年(不含1996年和1997年)，共計25 a。1990—2016年黃河流域年均地表水耗水量占年均天然徑流量的65.44%，但是部分年份占比極大。例如，2000—2002年連續枯水年占比高達94.88%，最枯的2002年天然徑流量僅246.16億m3，但地表水耗水量達到286.05億m3，是天然徑流量的116.20%。黃河除了供水外，還需要在河道內分配大量水量用于輸沙和生態保護，年均入海水量154.80億m3，即使在最枯的2002年，入海水量也有41.90億m3。1990—2016年，有12 a黃河流域地表水耗水量與入海水量之和超過了天然徑流量，在這些年份必須通過水庫補水才能保障河道內、外供水。在12 a中有11 a龍羊峽水庫補水，累計補水403.19億m3，年均補水36.65億m3。

2000—2002年的年均天然徑流量僅為289.31億m3，比多年平均值偏低40.96%；年均地表水耗水量274.51億m3，年均入海水量45.67億m3，天然徑流量與河道內外耗水量間存在30.87億m3的缺口。龍羊峽水庫在3年間累計補水85.35億m3，年均補水28.45億m3。如果沒有龍羊峽水庫調蓄，且維持河道外地表水耗水量不變，那么2000年、2001年入海水量將分別從48.58億、46.53億m3減少到25.25億、28.98億m3，2002年入海水量將為0。黃河天然徑流量年際變化大，易發生連續枯水年，而受到河流含沙量大的影響，黃河干流中下游修建多年調節水庫的難度大，作為黃河干流唯一一座多年調節水庫，龍羊峽水庫的蓄豐補枯對保障流域供水安全和黃河健康生命意義重大。

4 結 論

本文研究了1986—2016年龍羊峽水庫出庫徑流的IHA、流量過程變異程度等指標值相對于入庫徑流的變化，并識別了水庫蓄補水與來水、耗水等因素的相關關系。得到以下結論：

龍羊峽水庫顯著改變了出庫徑流過程，出庫徑流相對于入庫徑流的整體改變度達到高度改變，突出表現為流量波動范圍減小，包括年徑流量極值比減小24.84%、年徑流量變異系數減小26.92%、極端高流量事件的流量量級減小34.74%～47.30%、極端低流量事件的流量量級增大38.99%～144.43%等，但日流量的變化加快、波動加劇。

龍羊峽水庫蓄水導致出庫斷面徑流過程嚴重偏離天然情況。水庫出庫斷面的流量過程變異程度在水庫初次蓄水的1986年發生突變，1980—1985年和1987—2016年兩個時段出庫斷面流量過程變異程度分別為0.05和2.15，賦分分別為98.29和20.14。

龍羊峽水庫年補水量與黃河天然徑流量、黃河流域地表水耗水量占天然徑流量的比例之間均存在顯著的線性關系，其中水庫補水量與黃河天然徑流量間呈負相關關系。在龍羊峽水庫補水的年份中，有84.62%的年份黃河流域地表水耗水量占天然徑流量的比例超過64%。在龍羊峽水庫蓄水的年份中，有75.00%的年份黃河流域地表水耗水量占天然徑流量的比例低于64%。