馮菊霞,王桂榮,周 瑤
(1.宜川縣水利工作隊,陜西 延安 716200;2.寶雞市地下水管理監(jiān)測中心,陜西 寶雞 721001;3.陜西省渭河生態(tài)區(qū)保護中心,陜西 西安 710004)
劉莊水庫位于仕望河上游西川河支流丁盤溝上的劉莊村。仕望河發(fā)源于黃龍縣圪臺街孟家山,從河源至宜川縣城段為自南向北流,在縣城至黨灣鄉(xiāng)之間納入西川河折向東流,在秋林以西納入交里河,向東流入黃河,屬黃河一級支流,全長112.8 km,面積2 356 km2,全河平均比降9.7‰,仕望河上游森林密布,植被較好,清水長流,含沙量較小,河道多為土基河床,縣城以下河床下切較深,基巖出露,形成裸露的石峽谷,基巖出露區(qū)植被差,水土流失嚴(yán)重。西川河于縣城附近與仕望河交匯,全長52.1 km,面積519.6 km2,平均比降7.7‰。丁盤溝河長24.1 km,面積125.3 km2,平均比降11.7‰。莊水庫以上屬黃土高原稍林區(qū),森林茂密,植被覆蓋率高,劉莊水庫控制流域面積100 km2,壩址以上河流長度15.0 km,河道比降為10‰。
劉莊水庫地處暖溫帶半濕潤大陸性季風(fēng)氣候區(qū),冬春寒冷少雨,夏季多暴雨天氣,秋季常有連陰雨,其多年平均值為604.5 mm。多年平均氣溫9.9℃,1月平均氣溫-5.7℃,7月平均氣溫23.6℃,極端最高氣溫39.9℃,極端最低氣溫-22.4℃,無霜期178 d,早霜多出現(xiàn)在10月中旬,晚霜終于4月中旬。最大凍土深0.5 m,最大風(fēng)速18 m/s。
劉莊水庫位于黃龍山以北的黃土林區(qū),水庫以上森林植被好,河川徑流年內(nèi)分配較均勻,洪水徑流系數(shù)小,下滲損失大。年降水量的70%集中于6-9月,由于地處大陸腹地,遠離水汽源地,故降水年際變化大,且年內(nèi)分配不均勻,因而徑流洪水的年際變化大,洪水多發(fā)生在7、8兩月,連陰雨多發(fā)生在9月。但劉莊水庫壩址處無實測水文資料,水庫自建成以來未進行系統(tǒng)的庫區(qū)水文觀測。只有劉莊水庫所在河流—仕望河建有大村水文站,該站位于水庫下游,集水面積2 141 km2,為對水庫徑流計算方法進行研究探討,選取水文序列比較成熟的大村水文站1959-2000年共42 a的實測水文資料系列,綜合水庫附近最能表征流域降水特征的英旺雨量站的1961-2000年共40年的實測降水資料和流域內(nèi)宜川氣象站1952-2000年49 a的系列資料進行計算。選取的水文資料可靠性、代表性和一致性滿足規(guī)范要求,能作為計算的依據(jù)。
對收集到的仕望河大村水文站1959-2000年42 a實測年徑流資料系列進行審查,未發(fā)現(xiàn)不合理現(xiàn)象。
大村水文站1959-1980年徑流資料的還原直接采用《延安地區(qū)水文手冊》的計算成果,1981-2000年徑流資料還原按水資源評價細(xì)則進行。經(jīng)計算,大村水文站42 a還原水量平均為0.03×104萬 m3,占實測年徑流量的3.8%,占天然年徑流量的3.7%。
將還原后的徑流系列,用算術(shù)平均法求均值,用矩法求變差系數(shù)Cv,經(jīng)計算機程序計算,目估適線,大村站年徑流多年平均值為0.807億 m3,Cv值為0.57,Cs/Cv=3.0。計算結(jié)果見圖1。計算可知 Cv=0.57,Cs/Cv=3.0,平均值=0.807。
圖1 大村站年徑流量頻率曲線
基于流域資料的完整程度和資料審查合格的基礎(chǔ)上,本次徑流計算選用水文比擬法與等值線圖法分別進行,最后對兩種方法計算的徑流成果進行比選,得到最合理的成果。
2.4.1 水文比擬法
劉莊水庫所在河流下游設(shè)有大村水文站,流域面積2 141 km2,多年平均徑流量8 070×104萬 m3,面平均雨量610 mm。用下式為計算劉莊水庫面平均雨量580 mm時的徑流量:
W劉莊=W參×F劉/F參×H劉/H參
(1)
式中:W劉、W參為劉莊水庫、參證站多年平均年徑流量(104m3);F劉、F參為劉莊水庫、參證站控制集水面積(km2); H劉、H參為劉莊水庫、參證站面平均雨量(mm)。
經(jīng)計算,劉莊水庫多年平均年徑流量為358×104m3。
2.4.2 等值線圖法
通過查閱“九五”國家重點科技攻關(guān)項目:“陜西省地表水資源調(diào)查評價研究”的多年平均徑流深等值線圖可知,劉莊水庫重心處的徑流深為35 mm,則劉莊水庫多年平均年徑流量為350×104m3。
從以上2種方案的計算結(jié)果可以看出,水文比擬法和等值線圖法算得劉莊水庫設(shè)計年徑流量結(jié)果,基本一致,相差很小,說明結(jié)果可信。由于大村水文站屬于劉莊水庫本流域站,故選用該值。即劉莊水庫壩址處多年平均年徑流量為358×104m3。
由陜西省“九五”國家重點科技攻關(guān)項目查得,劉莊水庫年徑流變差系數(shù)Cv=0.45,Cs/Cv=3.0,求得劉莊水庫設(shè)計頻率年徑流量見表1。
表1 劉莊水庫設(shè)計頻率年徑流量成果表 104m3
根據(jù)大村站的頻率計算成果,在理論頻率線上初選頻率為95%、75%、50%年份相應(yīng)的徑流量,再根據(jù)還原后的實際徑流量,比選較為安全合理的年份,最終確定為1995、1973、1962年為相應(yīng)年份。按水文比擬法的原理,據(jù)其年內(nèi)分配百分比計算劉莊水庫設(shè)計頻率年徑流的月年分配結(jié)果,見表2。
表2 劉莊水庫設(shè)計頻率年徑流量年內(nèi)分配 104m3
劉莊河水庫所在的流域主要由黃土林區(qū)和石山林區(qū)相間,面積基本接近,缺洪水資料,對設(shè)計洪水的計算分析可采用以下所述的計算方法,參數(shù)用黃土林區(qū)和石山林區(qū)參數(shù)的平均值。
洪水總量計算可采用《延安地區(qū)實用水文手冊》提供的經(jīng)驗公式(2):
WP=KnQPn
(2)
式中:WP為洪水總量(萬m3);QP為洪峰流量(m3/s);Kn、n為經(jīng)驗參數(shù)。
查《延安地區(qū)實用水文手冊》,石山林區(qū)Kn=0.25,n=1.36;黃土林區(qū)Kn=0.18,n=1.27,參數(shù)按照黃土林區(qū)和石山林區(qū)參數(shù)的平均值,則Kn=0.215,n=1.315。
WP=0.215QP1.315
(3)
洪峰流量計算,同樣采用《延安地區(qū)實用水文手冊》提供的經(jīng)驗公式,并用綜合指數(shù)法與匯水面積相關(guān)法兩種方法進行校核。
3.2.1 綜合指數(shù)法
用式(4)計算重現(xiàn)期為N年的年最大流量QMN
QMN=C·Nα·Fβ·ψγ·H3Nη
(4)
式中:H3N為重現(xiàn)期為N年的3h面雨量;F=100 km2;Ψ = F / L2;L為河道長度,L=15 km。
查《延安地區(qū)實用水文手冊》各分區(qū)參數(shù)指數(shù)見表3。對各參數(shù)實際選取按照黃土林區(qū)和石山林區(qū)參數(shù)的平均值選取,然后進行設(shè)計洪峰流量計算,結(jié)果見表4。
表3 參數(shù)指數(shù)結(jié)果
表4 綜合指數(shù)法設(shè)計洪峰流量
3.2.2 匯水面積相關(guān)法
利用公式(5)計算重現(xiàn)期為N年的年最大流量QMN
QMN=KN·Fn
(5)
式中:KN、n為重現(xiàn)期為N年的經(jīng)驗參數(shù),查《延安地區(qū)實用水文手冊》洪峰面積相關(guān)參數(shù)表,500 a一遇石山林區(qū)KN=22.0,n=0.74,黃土林區(qū)KN=6.32,n=0.67;50 a一遇石山林區(qū)KN=12.0 ,n=0.74,黃土林區(qū)KN=3.58 ,n=0.67。參數(shù)按照黃土林區(qū)和石山林區(qū)參數(shù)的平均值,則500 a一遇Kn=14.16,n=0.705;50 a一遇KN=7.79,n=0.705。
按照公式(5)計算本流域年最大流量可得:500 a一遇洪峰流量364 m3/s;50 a一遇洪峰流量200 m3/s。
從設(shè)計洪水計算成果綜合表5可看出,500 a一遇年最大流量299~364 m3/s;50 a一遇年最大流量146~200 m3/s,考慮安全,采用匯水面積相關(guān)法,即500 a一遇洪峰流量為364 m3/s,洪水總量為504萬 m3;50 a一遇洪峰流量為200 m3/s,洪水總量為242萬 m3。
表5 洪水計算成果綜合表
洪水過程采用5點概化過程線,設(shè)計洪水流量過程見表6和表7。
表6 頻率P=0.2%設(shè)計洪水流量過程
表7 頻率P=2.0%設(shè)計洪水流量過程
對于一些缺乏水文實測資料的地區(qū),其徑流及設(shè)計洪水計算成果的精度和合理性非常重要,尤其對水資源開發(fā)利用和水利工程設(shè)計關(guān)系重大。通過以劉莊水庫為例,針對壩址處無實測水文資料和水庫自建成以來未進行系統(tǒng)的庫區(qū)水文觀測的實際,借鑒類似工程的計算方法,對徑流及設(shè)計洪水計算過程進行了闡述,從計算成果看,合理可行,符合劉莊水庫實際。