馬云峰
遼寧省河庫管理服務(wù)中心(遼寧省水文局),遼寧 沈陽 110003
遼寧省現(xiàn)有大中型水庫106座,總庫容173.10億m3,總調(diào)洪庫容84.60億m3,興利庫容93.70億m3。這些水庫與下游河道堤防相結(jié)合形成的防洪體系,在各個流域的防洪中發(fā)揮著重要作用,可基本保證城市和重要河段的防洪安全。例如:渾河上的大伙房水庫,觀音閣水庫、葠窩水庫,大凌河上的閻王鼻子水庫等,均對下游城市汛期防洪減災(zāi)起著重要作用。
1.2.1 防汛形勢新要求
隨著經(jīng)濟社會的發(fā)展,許多河流兩岸的保護對象越來越集中、越來越重要,流域防洪體系在保護人民生命財產(chǎn)安全中發(fā)揮的作用也越來越大。近些年遼寧省對防洪薄弱環(huán)節(jié)加大投資,開展了一大批水庫除險加固、河道整治工程建設(shè),流域防洪能力得到了明顯提升。但如何更有效地發(fā)揮防洪工程特別是大中型水庫的防洪效益,對水庫防汛工作提出了新的要求。另外,近幾年遼寧省汛期臺風(fēng)頻發(fā),2020年甚至出現(xiàn)了10日內(nèi)3次臺風(fēng)來襲的極端天氣,防汛形勢越發(fā)嚴峻。開展大中型水庫抗暴雨能力計算,提高精準防汛能力,是利用現(xiàn)代化手段進行科學(xué)研判和正確決策的重要手段。
1.2.2 水庫調(diào)度的重要參考
目前,中國氣象部門3日內(nèi)降雨預(yù)報精度已較高,對水庫而言,如果能在強降雨來臨前就計算出當(dāng)前剩余防洪庫容能夠滯納的降雨量,再對比氣象部門發(fā)布的降雨預(yù)報,盡早發(fā)布預(yù)警并適當(dāng)調(diào)整水庫調(diào)洪方式,從而最大限度地減少水庫由于調(diào)度不及時而出現(xiàn)的潰壩或上下游淹沒等問題。通過系統(tǒng)分析全省大中型水庫的抗暴雨能力,為各水庫防汛調(diào)度提供參考,可提高防汛決策的科學(xué)性和時效性。特別是對于一些中型水庫,由于水庫經(jīng)濟效益不高,管理經(jīng)費不足,而且沒有國家專項資金支持,維修養(yǎng)護不及時,更需要參考抗暴雨能力控制運行。
水庫所在流域下墊面情況,即前期降雨量和土壤飽和度,土壤能夠承納的降雨越多,則匯入水庫的降雨量越小,水庫的抗暴雨能力就越大。遼寧省東部山區(qū)屬長山白余脈,植被覆蓋率較高,土壤保水性好,而遼西山區(qū)植被覆蓋差,土壤水分蒸發(fā)快,同等降雨量遼東山區(qū)水庫抗暴雨能力小于遼西地區(qū)。
水庫水位越高,剩余庫容越小,所能承納的降雨量越小,抗暴雨能力就越小。因此,汛期要嚴密監(jiān)視水庫水位,在強降雨前采取措施降低水位,騰出足夠庫容,提高水庫抗暴雨能力。
調(diào)度方式?jīng)Q定了水庫的泄洪能力,泄洪越快則抗暴雨能力越大。但與小型水庫不同,大中型水庫最大特點是大多具有溢洪道閘門,汛限水位一般高于堰頂高程,汛期需要重點考慮調(diào)度方式和起調(diào)水位。另外,水庫的下泄流量需要考慮下游河道的行洪能力,溢洪道泄量直接決定水庫降雨承納能力。
由于不同流域下墊面的植被保水率、土壤特性及含水率不同,不同的降雨過程和降雨量會導(dǎo)致不同的產(chǎn)流過程和產(chǎn)流量。例如:在遼西地區(qū),一些地方雖然土壤含水低較低,但若遇高強度短歷時降雨,土壤可能形成“土蓋”,土壤未飽和便開始產(chǎn)流,形成較大的產(chǎn)流量;而遇長歷時較均勻的降雨過程,則大部分的降雨可能會下滲或蒸發(fā),反而產(chǎn)生較少的徑流量。
根據(jù)大中型水庫抗暴雨能力影響因素,結(jié)合水庫工情,抗暴雨能力計算過程分為以下5步,如圖1所示。
圖1 水庫抗暴雨能力分析技術(shù)路線圖
3.1.1 直接法計算
當(dāng)水庫以上流域內(nèi)長期雨量站分布較密,可根據(jù)水庫以上流域內(nèi)各年的年最大24 h面雨量系列直接進行頻率分析計算,得出各重現(xiàn)期的24 h設(shè)計面雨量。
3.1.2 間接法計算
對于無雨量站或雨量資料短缺的大中型水庫,由《遼寧省中小河流(無資料地區(qū))設(shè)計暴雨洪水計算方法》中查得流域重心24 h暴雨均值、Cv值、點面折減系數(shù),進而計算各重現(xiàn)期設(shè)計24 h面暴雨。
3.2.1 劃分暴雨量級
水庫24 h面暴雨量級從50 mm開始,按50 mm逐級遞增,大部分水庫最大降雨量級做到500 mm,降雨高值地區(qū)的最大降雨量級可做到1 000 mm。
3.2.2 收集洪水過程線
收集水庫初步設(shè)計或除險加固設(shè)計中采用的設(shè)計洪水及設(shè)計洪水過程線,補全各頻率設(shè)計洪水及設(shè)計洪水過程線。
3.2.3 推求量級暴雨對應(yīng)洪水過程線
假定設(shè)計暴雨與設(shè)計洪水同頻,依據(jù)水庫24 h設(shè)計面暴雨、各重現(xiàn)期設(shè)計洪水及洪水過程線成果,用內(nèi)插法推求不同量級24 h面降雨對應(yīng)的洪水過程線。
統(tǒng)計水庫站或水文站歷年場次24 h面降雨、降雨開始時刻土壤含水量及三日洪量,一年可選多場,三日洪量指由本場24 h降雨形成的三日洪量,即由統(tǒng)計的三日洪量減去前次洪水的退水量。其中土壤含水量Pa計算按下式計算:
Pa,t=KPa,t -1+ K2Pa,t -2+ …+ KnPa,t -n
(1)
式(1)中,Pa,t為前期影響雨量;Pa,t -1為上一天的Pa值;n統(tǒng)一取15 d;K統(tǒng)一取0.85。
當(dāng)所計算出來的Pa值大于Imax時,則Pa值取Imax,各水庫及水文站Imax參考《遼寧省中小河流(無資料地區(qū))設(shè)計暴雨洪水計算方法》中各水文分區(qū)數(shù)值。不同飽和情況分為三級,即降雨開始時刻的土壤含水量Pa=0 mm、Pa=0.5Imax及Pa=Imax。飽和狀態(tài)下按24 h降雨全部形成三日洪量考慮,相關(guān)線為45度線。
水庫資料條件好的大型水庫可直接繪制本水庫的降雨徑流相關(guān)圖,無實測流量資料或資料不足的水庫,先繪制附近參證水文站的降雨徑流關(guān)系圖,考慮降雨及徑流的地區(qū)差異性,進行相關(guān)的分析處理,轉(zhuǎn)換成水庫的降雨徑流相關(guān)圖。降雨徑流相關(guān)圖中的Pa等值線根據(jù)實測點據(jù)用目估內(nèi)插法求出。
不同飽和情況不同量級的24 h降雨,由水庫降雨徑流相關(guān)圖查得相應(yīng)的三日洪量,與水庫各重現(xiàn)期3日洪量比較,用內(nèi)插法推算不同飽和狀態(tài)不同量級降雨對應(yīng)的洪水過程線。
水庫洪水調(diào)度方式采用已批復(fù)的水庫設(shè)計報告、除險加固設(shè)計報告或度汛方案中確定的調(diào)度方式。對于單個水庫而言,采用二分迭代法進行水庫調(diào)洪演算,分別按照考慮壩體安全、考慮下游控制泄量情況下,特定工況某一起調(diào)水位下的降雨量級和最高水位關(guān)系圖。為便于成果使用,抗暴雨能力分析,在可能暴雨分析基礎(chǔ)上,采用逐級試算法,給出不同工況下的水庫洪水調(diào)節(jié)成果,并給出最大抗暴雨能力級別。
防汛決策是關(guān)系國計民生的大事,而水庫防洪調(diào)度作為防汛決策的重要部分,具有很強的時效性要求,需要實時預(yù)報,滾動修正調(diào)度。大中型水庫抗暴雨能力計算結(jié)果可作為水庫汛期調(diào)度決策的重要參考,一定程度上彌補防汛預(yù)報調(diào)度短板。在實際應(yīng)用中,一場洪水調(diào)度過程中的水庫最高水位、下游最大泄量、調(diào)度期末水庫水位等都需要綜合考慮。