楊海燕,陶興玲,蔣金銘,萬 芳
(1.四川省達(dá)州水文水資源勘測中心,四川 達(dá)州 635000;2.華北水利水電大學(xué)水利學(xué)院,鄭州 450046)
四川省江口水庫位于達(dá)州市前河、中河、后河匯口處宣漢縣城,以發(fā)電為主,兼有防洪、灌溉、發(fā)電等功能,為州河流域僅有一定調(diào)蓄能力的水庫,對州河下游城鎮(zhèn)(特別是達(dá)州市城區(qū))的防洪安全具有十分重要的作用。同時(shí)江口水庫以下人煙稠密,主要有羅江鎮(zhèn)、達(dá)州市城區(qū)、達(dá)州機(jī)場等,而州河為典型的山區(qū)河流,流域暴雨洪水頻繁,城市防洪能力較低,災(zāi)害嚴(yán)重,威脅人民生命財(cái)產(chǎn)安全,制約經(jīng)濟(jì)發(fā)展。因此,對江口水庫入庫流量預(yù)報(bào)及水庫調(diào)度的模擬研究,為州河流域水庫的實(shí)時(shí)洪水調(diào)度提供有力有效的技術(shù)支撐,促進(jìn)水庫的安全經(jīng)濟(jì)運(yùn)行,避免或盡量減少因洪水對人民群眾生命和財(cái)產(chǎn)造成的損失,同時(shí),也有利于實(shí)現(xiàn)水資源的優(yōu)化配置。本文基于近期投入使用的“四川省水旱災(zāi)害防御決策支持系統(tǒng)”,選取州河江口水庫為例,對江口水庫的入庫流量預(yù)報(bào)及調(diào)度進(jìn)行動(dòng)態(tài)模擬研究。
針對流域水旱災(zāi)害特點(diǎn)和流域水旱災(zāi)害防御決策指揮需求,近年來四川省水利廳投資建設(shè)的新平臺——“四川省水旱災(zāi)害防御決策支持系統(tǒng)”已于2022年逐步投入使用。該平臺融合了水務(wù)及水文等多重模塊,融合了實(shí)時(shí)雨水情、氣象網(wǎng)格數(shù)值預(yù)報(bào)、流域洪水預(yù)報(bào)、動(dòng)態(tài)洪水風(fēng)險(xiǎn)編制、流域防洪調(diào)度、洪水調(diào)查等內(nèi)容,為進(jìn)一步分析研判汛情風(fēng)險(xiǎn)、預(yù)測災(zāi)害發(fā)展趨勢及水庫調(diào)度決策指揮提供多層次的信息服務(wù)和多種手段支持。
“四川省水旱災(zāi)害防御決策支持系統(tǒng)”的水文模塊中除了涵蓋原老平臺的雨水情查詢功能外,還添加了預(yù)報(bào)、調(diào)度、簡報(bào)發(fā)布等多重功能模塊。其中“預(yù)報(bào)調(diào)度”模塊通過在平臺中構(gòu)建相應(yīng)的洪水預(yù)報(bào)模型,結(jié)合未來氣象數(shù)值預(yù)報(bào),即可預(yù)報(bào)水庫的入庫流量。在平臺中錄入水利工程基本信息后,便可在入庫流量預(yù)報(bào)的基礎(chǔ)上實(shí)現(xiàn)水庫調(diào)度的動(dòng)態(tài)模擬。水庫的入庫流量預(yù)報(bào)及調(diào)度模擬將為流域水旱災(zāi)害防御工作提供科學(xué)的技術(shù)支撐,充分發(fā)揮水庫防洪蓄水治澇的能力,為防洪減災(zāi)和水資源合理利用發(fā)揮重要作用。
州河屬渠江左岸一級支流,上源分前、中、后河三支,均發(fā)源于大巴山南麓,河長310 km,控制流域面積11 165 km2,經(jīng)宣漢、達(dá)縣至三匯鎮(zhèn)匯入渠江。前河為主流,發(fā)源于城口縣燕麥鄉(xiāng)光頭山,自東北向西南流,至宣漢縣城附近與后河匯合,主河道長178 km,集水面積2 249 km2;中河發(fā)源于重慶城口縣白芷鄉(xiāng)的白芷山西南,中河流域面積1 402 km2,從源頭至匯合口河長101.5 km;后河起源于萬源市皮窩鄉(xiāng)趙家山下鐵鋪,于宣漢縣城東的江口處與前河相匯,注入州河,集水面積2 987 km2,河長154 km。
州河流域?qū)賮啛釒駶櫦撅L(fēng)氣候,流域多年平均降水量1 308.4 mm。年內(nèi)降水主要集中在5-10月,降水量占全年總降水量的80%以上;而暴雨發(fā)生在7月、9月的概率較大,因此,年內(nèi)變化呈雙峰過程。
徑流在年內(nèi)的變化與降水的年內(nèi)變化相應(yīng)。5-11月為豐水期,該時(shí)段內(nèi)的多年平均徑流量占年徑流量的86.8%;其中7-9月最豐,占年水量的49.5%。暴雨洪水匯流時(shí)間短,州河洪水具有陡漲陡落、峰高量大、過程線尖瘦等特點(diǎn);過程線形狀單峰洪水歷時(shí)2~3 d,復(fù)峰歷時(shí)約5~6 d。
江口水庫樞紐工程系州河干流梯級開發(fā)中的第一級,地處宣漢縣城區(qū)上游,前、后河匯合口下游約300 m州河源頭處,距離達(dá)州城區(qū)約51 km,是一座以發(fā)電為主,兼有防洪、灌溉的骨干水利工程。集雨面積6 364 km2,總庫容2.77億m3,有效庫容1.48億m3,調(diào)洪庫容1.33億m3,具有不完全年調(diào)節(jié)性能[1]。
江口水庫擋水和泄水建筑物設(shè)計(jì)洪水標(biāo)準(zhǔn)為100年一遇,洪峰流量12 500 m3/s;校核洪水標(biāo)準(zhǔn)為1000年一遇,洪峰流量為15 800 m3/s。水庫正常蓄水位330.0 m,設(shè)計(jì)洪水位330.62 m,校核洪水位332.70 m,主汛期汛限水位327.00 m。
水文站網(wǎng)是短期洪水預(yù)報(bào)的基礎(chǔ)。州河上游分為前河、中河、后河主要三支流,江口水庫位于前河、中河、后河匯口處,前河上有南壩水文站,中河有黃金口水文站,后河有毛壩水文站,以及眾多的雨量監(jiān)測站。截至目前,達(dá)州水文中心在江口水庫以上流域基本建成了水情遙測站126處(其中基本水文站4處,中小河流水文站2處,中小河流水位站3處,雨量站117處)。流域水系及站點(diǎn)分布見圖1,站點(diǎn)分布全流域,站網(wǎng)密度大,滿足水庫入庫流量預(yù)報(bào)需求。
圖1 江口水庫所在流域水系及水文站網(wǎng)分布
河系預(yù)報(bào)主要為江口水庫以上流域支流各主要控制站點(diǎn)逐級演算至江口水庫處,并疊加區(qū)間流量過程得到江口電站的入庫流量。其中,河系方案中各一級預(yù)報(bào)站點(diǎn)及區(qū)間來水采用新安江三水源蓄滿產(chǎn)流模型(SMS_3)和三水源滯后演算模型(LAG_3),各站點(diǎn)逐級演算采用馬斯京根河道流量演算法[2]。
本次方案主要借用“中國洪水預(yù)報(bào)系統(tǒng)”分別率定前、中、后河各支流控制站點(diǎn)南壩、黃金口、毛壩三個(gè)基本水文站(即本河系預(yù)報(bào)方案中一級預(yù)報(bào)站點(diǎn))的新安江模型和馬斯京根河道演算各參數(shù),區(qū)間方案參數(shù)借用本流域具有相似流域下墊面條件站點(diǎn)已率定的新安江模型參數(shù),再通過“四川省水旱災(zāi)害防御決策支持系統(tǒng)”平臺中“預(yù)報(bào)調(diào)度”模塊構(gòu)建江口水庫的河系方案。河系預(yù)報(bào)方案概化圖見圖2。
圖2 江口水庫河系預(yù)報(bào)方案
由于江口水庫沒有長系列實(shí)測水文資料,本次無法對該方案進(jìn)行檢驗(yàn)評定,故只對本河系預(yù)報(bào)方案中一級預(yù)報(bào)站點(diǎn)的新安江模型方案進(jìn)行檢驗(yàn)評定,檢驗(yàn)評定見表1。檢驗(yàn)評定表明,各一級預(yù)報(bào)站的方案精度均可達(dá)到乙級以上。
表1 基本水文站新安江模型方案檢驗(yàn)評定
本次河系預(yù)報(bào)方案依據(jù)《水文情報(bào)預(yù)報(bào)規(guī)范》(GB/T 22482-2008),利用江口水庫的多場次實(shí)測洪水入流過程對預(yù)報(bào)方案進(jìn)行檢驗(yàn)后再做精度評定[3]:方案精度達(dá)乙級以上,可做預(yù)報(bào)及滾動(dòng)預(yù)報(bào)使用。
4.1.1 2021年“8·29”洪水
2021年8月28~29日,達(dá)州轄區(qū)普降大到暴雨、局部大暴雨,本次暴雨時(shí)空分布較均勻,持續(xù)時(shí)間較長,暴雨中心主要集中在萬源市、宣漢縣、達(dá)川區(qū)、開江縣一帶;最大1h雨量為22.0 mm,最大日雨量為215.0 mm(28日);最大過程雨量為293.0 mm。受暴雨影響,28~29日,達(dá)州轄區(qū)境內(nèi)各干支流發(fā)生了不同程度的洪水過程。其中:江口水庫以上各支流發(fā)生較大洪水過程,其最大入庫流量8 600 m3/s,最高庫水位超汛限水位2.55 m。
4.1.2 2022年“10·5”洪水
2022年10月4~6日,達(dá)州市普降大到暴雨、局部大暴雨。暴雨中心主要集中在達(dá)川區(qū)、萬源市、宣漢縣、通川區(qū)一帶,暴雨時(shí)間主要集中在4日0時(shí)至6日8時(shí)。達(dá)州轄區(qū)最大過程雨量為298.0 mm,最大日雨量為170.5 mm(4日),最大1 h雨量為4日16時(shí)35.0 mm。受強(qiáng)降雨影響,江口水庫以上河流發(fā)生了一次一般洪水過程,江口水庫最大入庫流量4 800 m3/s,最高庫水位超汛限水位2.39 m。
根據(jù)渠江流域洪水預(yù)報(bào)多年經(jīng)驗(yàn),本文模擬實(shí)際預(yù)報(bào)過程,以檢驗(yàn)方案預(yù)報(bào)調(diào)度是否合理可用,故預(yù)報(bào)時(shí)間為降雨開始后、各一級預(yù)報(bào)站點(diǎn)未出洪峰之前?!八拇ㄊ∷禐?zāi)害防御決策支持系統(tǒng)”平臺的具體操作流程為:首先,用“預(yù)報(bào)調(diào)度”中“預(yù)見期降雨設(shè)置”模塊獲取未來氣象數(shù)據(jù),同時(shí),可對氣象數(shù)據(jù)進(jìn)行人工修正,并保存雨量方案;其次,在“預(yù)報(bào)調(diào)度”模塊中選擇已構(gòu)建好的“江口水庫”的河系方案,在獲取雨量方案及其他數(shù)據(jù)后進(jìn)行作業(yè)預(yù)報(bào),同時(shí),可根據(jù)實(shí)際情況對預(yù)報(bào)結(jié)果進(jìn)行檢查和人工修正;最后,在江口水庫入庫流量預(yù)報(bào)結(jié)果的基礎(chǔ)上,通過模擬不同的出庫方式,來演算水庫入庫流量、水庫水位等變化。
江口水庫動(dòng)態(tài)模擬中,“8·29”洪水以2021年8月28日22時(shí)為預(yù)報(bào)依據(jù)時(shí)間,“10·5”洪水以2022年10月4日20時(shí)為預(yù)報(bào)依據(jù)時(shí)間,前期預(yù)熱期為30 d,預(yù)見期時(shí)段長為24個(gè)時(shí)段。利用已構(gòu)建的江口水庫河系預(yù)報(bào)方案,在獲取未來氣象數(shù)值天氣數(shù)據(jù)的前提下,進(jìn)行江口水庫的入庫流量預(yù)報(bào)。為了檢驗(yàn)水庫調(diào)度模擬是否與實(shí)際過程相符,本次的出庫流量模擬按照實(shí)際出庫流量設(shè)置,“8·29”“10·5”入庫流量及水庫調(diào)度動(dòng)態(tài)模擬見表2、表3。
表2 “8·29”洪水江口入庫流量及水庫調(diào)度動(dòng)態(tài)模擬
表3 “10·5”洪水江口入庫流量及水庫調(diào)度動(dòng)態(tài)模擬
由江口水庫入庫流量及水庫調(diào)度動(dòng)態(tài)模擬結(jié)果分析:入庫流量動(dòng)態(tài)模擬中洪峰預(yù)見期可達(dá)7~10 h,入庫洪峰流量數(shù)值誤差小于5%,時(shí)間誤差為1h,誤差滿足預(yù)報(bào)要求;同時(shí),在入庫流量過程預(yù)報(bào)結(jié)果的基礎(chǔ)上,根據(jù)實(shí)際出庫流量模擬得到的江口水庫水位動(dòng)態(tài)變化過程與實(shí)際水位過程較吻合,效果良好。
準(zhǔn)確可靠的水庫入庫流量及調(diào)度動(dòng)態(tài)模擬,可為及時(shí)、科學(xué)的水庫群聯(lián)合調(diào)度提供有力的技術(shù)支撐,以達(dá)到削減洪峰流量、降低洪峰水位、實(shí)現(xiàn)水庫精細(xì)化調(diào)度、水資源科學(xué)配置等最終目的?;凇八禐?zāi)害防御決策支持系統(tǒng)”平臺建立的江口水庫河系預(yù)報(bào)方案,實(shí)現(xiàn)了江口水庫的入庫流量預(yù)報(bào)及水庫調(diào)度進(jìn)行動(dòng)態(tài)模擬,為水庫的實(shí)時(shí)洪水調(diào)度提供數(shù)據(jù)支撐及直觀的圖形顯示,為類似水庫的入庫流量預(yù)報(bào)調(diào)度提供參考?;凇八禐?zāi)害防御決策支持系統(tǒng)”平臺模擬未來水情信息時(shí),對未來數(shù)字天氣預(yù)報(bào)依賴程度較高;同時(shí),該方法須基于已率定好模型參數(shù)的基礎(chǔ)上構(gòu)建方案,因此,須在構(gòu)建河系方案前做大量的模擬檢驗(yàn)率定工作,以得到較高預(yù)報(bào)精度方案參數(shù)。