李慶國(guó),段高云,郭其峰,成鵬飛
(黃河勘測(cè)規(guī)劃設(shè)計(jì)研究院有限公司,河南鄭州450003)
黑河干流全長(zhǎng)928 km,在源流區(qū)分為東、西兩岔[1],在黃藏寺村匯合后折向北流,流域面積為11.6萬(wàn)km2,多年平均水資源總量為28億m3。
黃藏寺水利樞紐是黑河干流骨干調(diào)蓄工程,壩址上距黃藏寺村約11 km,是黑河流域重要的水資源配置工程、生態(tài)保護(hù)工程和扶貧開發(fā)工程。工程地理位置見圖 1。樞紐為Ⅱ等大(2)型工程,死水位為2 580.0 m,正常蓄水位為2 628.0 m,水庫(kù)總庫(kù)容為4.03億m3,調(diào)節(jié)庫(kù)容為2.95億m3。電站裝機(jī)規(guī)模為49 MW,保證出力為 6.2 MW,多年平均發(fā)電量為2.03億kW·h。大壩為碾壓混凝土重力壩,最大壩高為123 m。
圖1 黃藏寺水利樞紐工程位置示意
黃藏寺水利樞紐壩址以上東西兩岔八寶河、黑河上分別建有祁連水文站和扎馬什克水文站[2-3],控制流域面積分別為2 452 km2和4 589 km2;壩址下游80 km處有鶯落峽水文站,控制流域面積為10 009 km2;黃藏寺壩址控制流域面積為7 648 km2。黑河流域主要水文站概況見表1。
表1 黑河流域主要水文測(cè)站概況
祁連站、扎馬什克站和鶯落峽站(以下簡(jiǎn)稱“三站”)多年平均來(lái)水量分別為4.58億、7.41億、16.09億m3,多年平均實(shí)測(cè)輸沙量分別為38.69萬(wàn)、113.55萬(wàn)、205.17萬(wàn) t,多年平均含沙量分別為 0.84、1.53、1.28 kg/m3,最大含沙量分別為321 kg/m3(1979年8月11日)、273 kg/m3(1959 年 7 月 13 日)、140 kg/m3(1992年7月15日),見表2。
流域水沙主要有以下特性:
(1)水量年際變化較小,沙量年際變化較大。從年際變化情況看,來(lái)水量變化較小。三站最大年來(lái)水量分別為6.94億、10.47億、23.07億 m3,分別為最小年來(lái)水量的2.2倍、2.0倍、2.3倍;來(lái)沙量年際變化較大,三站最大年來(lái)沙量分別為201.22萬(wàn)、641.65萬(wàn)、995.08萬(wàn)t,分別為最小年來(lái)沙量的19.2倍、43.2倍、271.1倍。
(2)水沙量年內(nèi)分配不均。各站水量年內(nèi)分布不均,主要集中在6—9月,三站多年平均6—9月來(lái)水量分別占年來(lái)水量的64%、69%、68%;來(lái)沙量的集中程度甚于來(lái)水量,三站多年平均6—9月來(lái)沙量分別為年來(lái)沙量的87%、97%、95%,見圖2、圖3。
表2 黑河流域主要水文站多年平均徑流泥沙特征值
圖2 黑河流域各主要水文站水量年內(nèi)分配
圖3 黑河流域各水文站沙量年內(nèi)分配
(3)來(lái)沙量主要集中在洪水期,非洪水期來(lái)沙量很小。以大水大沙的1996年為例,該年扎馬什克站實(shí)測(cè)最大洪峰流量為1 010 m3/s,為本站實(shí)測(cè)流量最大值;祁連站和扎馬什克站共來(lái)水13.90億m3,占兩站1968—2012年總來(lái)水量的 2.9%;來(lái)沙量為 702.27萬(wàn)t,占1968—2012年總來(lái)沙量的12%。從1996年祁連站和扎馬什克站5—9月日均流量及含沙量過(guò)程線可以看出,洪水期洪峰與沙峰對(duì)應(yīng)較好,其中兩站4場(chǎng)洪水來(lái)水量為5.04億m3,占全年來(lái)水量的36.3%;來(lái)沙量為654.41萬(wàn)t,占全年來(lái)沙量的93.2%,見圖4。
圖4 祁連站和扎馬什克站日均流量、含沙量過(guò)程線
(1)入庫(kù)懸移質(zhì)沙量。工程上游受人類活動(dòng)影響較小,因此黃藏寺水利樞紐入庫(kù)泥沙由祁連、扎馬什克站加區(qū)間來(lái)沙組成,區(qū)間沙量根據(jù)兩站與鶯落峽站區(qū)間輸沙模數(shù)推求,黃藏寺水利樞紐年均入庫(kù)懸移質(zhì)沙量為169.7萬(wàn)t。同理,算得黃藏寺1958—2012年多年平均入庫(kù)水量為12.81億m3。
(2)入庫(kù)推移質(zhì)沙量。黃藏寺水利樞紐懸移質(zhì)來(lái)沙量較小,但河道兩岸有泥石流情況發(fā)生,且?guī)靺^(qū)山體裸露、植被較差。確定黃藏寺水利樞紐推移質(zhì)沙量為懸移質(zhì)沙量的20%[5],即年入庫(kù)推移質(zhì)為33.9萬(wàn)t。利用梅葉-彼得和R·摩勒推移質(zhì)輸沙率計(jì)算公式(見式(1))計(jì)算年均入庫(kù)推移質(zhì)泥沙為35.6萬(wàn)t,與采用推懸比[5]的計(jì)算值比較接近。輸沙率計(jì)算公式為
式中:gs為推移質(zhì)單寬輸沙率,t/(s·m);Ks為河床糙率系數(shù);Kr為河床平整情況下的沙粒阻力系數(shù),Kr=為床沙顆粒級(jí)配曲線中90%較之為小的粒徑,取 0.25 mm;r、rs分別為水、沙容重,r=1 t/m3,rs=2.65 t/m3;h為水深,m;i為坡降;d為粒徑,m。
(3)水沙系列。黃藏寺水利樞紐為大型水庫(kù),泥沙淤積年限確定為50 a。黃藏寺壩址處有1968—2012年45 a實(shí)測(cè)入庫(kù)日均水沙資料,再利用45 a實(shí)測(cè)資料進(jìn)行5 a滑動(dòng)計(jì)算,即可從中選取5 a水沙系列組成50 a水沙系列。
綜合選取1968—2012年+1994—1998年50 a系列作為水庫(kù)泥沙沖淤計(jì)算采用系列,系列多年平均入庫(kù)水量為12.78億m3,多年平均入庫(kù)沙量為172.3萬(wàn)t,水沙量分別為設(shè)計(jì)水沙量的0.998倍、1.015倍,多年平均含沙量為1.35 kg/m3。黃藏寺水利樞紐推移質(zhì)沙量為懸移質(zhì)沙量的20%,則入庫(kù)推移質(zhì)為34.5萬(wàn)t。水沙系列特征值統(tǒng)計(jì)見表3。
表3 入庫(kù)懸移質(zhì)系列多年平均水沙特征值統(tǒng)計(jì)
(1)調(diào)度運(yùn)行方式。黃藏寺水利樞紐在4月上中旬、7月、8月中旬、9月中旬(相機(jī))通過(guò)加大流量向下游輸水,完成國(guó)務(wù)院批復(fù)的分水方案,滿足黑河下游生態(tài)用水要求,在調(diào)水的同時(shí)加強(qiáng)對(duì)中游引水口門的管理;在中游灌區(qū)用水時(shí)段按灌溉用水要求進(jìn)行灌溉調(diào)度,使中游地區(qū)農(nóng)業(yè)用水保證率達(dá)到設(shè)計(jì)保證率50%;冬四月(12月至翌年3月)泄放生態(tài)基流,兼顧發(fā)電運(yùn)用[4]。
當(dāng)水庫(kù)來(lái)水量大于用水需求時(shí),水庫(kù)蓄水運(yùn)用,汛期控制水位不超過(guò)汛限水位2 628 m;當(dāng)水庫(kù)來(lái)水量小于用水需求時(shí),水庫(kù)泄水運(yùn)用。非汛期水庫(kù)最高運(yùn)用水位不高于正常蓄水位2 628 m,最低水位不低于死水位2 580 m,特枯水年水庫(kù)水位不低于極限死水位2 560 m。
(2)水庫(kù)沖淤計(jì)算結(jié)果。根據(jù)選取的泥沙系列和水庫(kù)運(yùn)用方式,采用黃河勘測(cè)規(guī)劃設(shè)計(jì)研究院有限公司研發(fā)的一維水文水動(dòng)力學(xué)泥沙沖淤計(jì)算模型進(jìn)行計(jì)算,模型計(jì)算時(shí)段為日。經(jīng)計(jì)算,黃藏寺水庫(kù)50 a總?cè)霂?kù)懸移質(zhì)沙量為0.663億m3(懸移質(zhì)淤積物干容重采用1.3 t/m3),水庫(kù)50 a累計(jì)懸移質(zhì)淤積量為0.513億m3,懸移質(zhì)多年平均排沙比為22.6%,見圖5。推移質(zhì)泥沙按全部淤積考慮,水庫(kù)運(yùn)用50 a后推移質(zhì)淤積量為0.115億m3(推移質(zhì)淤積物干容重采用1.5 t/m3),則黃藏寺水利樞紐運(yùn)用50 a共淤積泥沙0.628億m3。
圖5 水庫(kù)累計(jì)沖淤過(guò)程
(3)水庫(kù)淤積形態(tài)。根據(jù)已建水庫(kù)的實(shí)測(cè)資料分析、經(jīng)驗(yàn)公式計(jì)算、懸移質(zhì)泥沙沖淤計(jì)算結(jié)果,綜合判斷水庫(kù)淤積形態(tài)為三角洲淤積。按照水庫(kù)淤積沙量,試算確定黃藏寺水利樞紐淤積形態(tài)[5-7]。采用經(jīng)驗(yàn)法進(jìn)行泥沙沿程分布設(shè)計(jì)計(jì)算,水庫(kù)運(yùn)用50 a后壩前淤積高程為2 557 m,尾部段和三角洲段占總淤積量的67.4%,沿程段淤積量占總淤積量的12.1%,壩前淤積量占總淤積量的20.5%。水庫(kù)淤積縱剖面見圖6、圖7。
(4)出庫(kù)含沙量。黃藏寺水利樞紐多年平均出庫(kù)含沙量為0.23 kg/m3,最大日出庫(kù)含沙量為6.84 kg/m3,見表4。
圖6 水庫(kù)運(yùn)用50 a后淤積縱剖面(黑河)
圖7 水庫(kù)運(yùn)用50 a后淤積縱剖面(八寶河)
表4 黃藏寺水庫(kù)出庫(kù)含沙量特征值統(tǒng)計(jì)
黑河流域具有水量年際變化小、沙量年際變化大、水沙年內(nèi)分配不均、沙量主要來(lái)自洪水期等特性。黃藏寺水利樞紐年入庫(kù)懸移質(zhì)沙量為169.7萬(wàn)t、推移質(zhì)沙量為33.9萬(wàn)t,泥沙問(wèn)題相對(duì)不嚴(yán)重。采用50 a水沙系列,應(yīng)用一維水文水動(dòng)力學(xué)泥沙沖淤計(jì)算模型進(jìn)行計(jì)算,水庫(kù)計(jì)算期排沙比為22.6%,50 a共淤積泥沙0.628億m3,最大日出庫(kù)含沙量為6.84 kg/m3。水庫(kù)淤積形態(tài)為三角洲淤積,壩前淤積高程為2 557 m。
黃藏寺水利樞紐是黑河干流骨干調(diào)蓄工程,是黑河流域重要的水資源配置工程、生態(tài)保護(hù)工程和扶貧開發(fā)工程,工程已于2016年開工建設(shè)。通過(guò)分析入庫(kù)水沙特性,選取泥沙系列,采用經(jīng)驗(yàn)分析和數(shù)學(xué)模型計(jì)算,合理確定水庫(kù)工程泥沙設(shè)計(jì)成果,為合理確定工程規(guī)模提供技術(shù)支撐。工程泥沙設(shè)計(jì)是基于選定的水沙條件和設(shè)計(jì)運(yùn)用方式,但水沙條件變化及工程調(diào)度運(yùn)行存在一定的不確定性,因此,在工程實(shí)際運(yùn)用過(guò)程中應(yīng)加強(qiáng)調(diào)度運(yùn)用管控和庫(kù)區(qū)淤積監(jiān)測(cè),保證工程效益正常發(fā)揮。