楊怡青，胡 劍，陳龍贊

(寧波市水利水電規(guī)劃設(shè)計(jì)研究院，浙江 寧波 315192)

在流域防洪體系中，大中型水庫一直承擔(dān)攔洪蓄洪、削峰錯(cuò)峰的工程任務(wù)，是整個(gè)防洪體系的重要組成部分[1]。近年來在極端氣候頻發(fā)、城市建設(shè)規(guī)模持續(xù)擴(kuò)張、流域規(guī)劃防洪格局改變等新形勢下，流域整體洪量增加，流域干流泄洪壓力明顯增大[2]。因此，提升大中型水庫防洪能力、挖掘其工程應(yīng)用價(jià)值是緩解干流防洪壓力、減少洪水入侵平原的重要舉措，是提升水利防災(zāi)減災(zāi)能力、支撐新形勢下經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展的內(nèi)在要求。

國內(nèi)外關(guān)于水庫防洪能力的研究起步于20世紀(jì)90年代，馮平等通過概率組合方法估算了水庫的防洪能力，得出的結(jié)果比僅依賴設(shè)計(jì)洪水過程給出的防洪標(biāo)準(zhǔn)更能反映水庫的實(shí)際防洪作用[3]。仲志余等基于動(dòng)庫容調(diào)洪的基本原理并考慮了泥沙淤積對水庫防洪能力的影響，計(jì)算分析了三峽水庫的防洪能力[4]。胡秀英利用多庫聯(lián)合調(diào)度模型計(jì)算分析多水庫聯(lián)合防洪能力[5]。S J Bennett等結(jié)合美國11000座水庫的實(shí)測資料，研究了累積沉積物(如泥沙)等對水庫防洪能力的影響[6]。Niewiadomska等基于多庫聯(lián)合調(diào)度模型研究了實(shí)時(shí)水庫防洪與預(yù)報(bào)方法[7]?？傮w來說，國內(nèi)外關(guān)于水庫防洪能力的研究經(jīng)歷了從靜態(tài)到動(dòng)態(tài)[8]，從單庫防洪能力到多庫聯(lián)合[9]防洪能力的變化[10]。

國內(nèi)外現(xiàn)有的研究均是從水庫上游來水的角度出發(fā)，結(jié)合水庫的運(yùn)行狀況分析優(yōu)化水庫的防洪能力[11]。本文將從水庫現(xiàn)有防洪能力、洪水特性分析等方面出發(fā)，結(jié)合下游興奉橋斷面安全泄量，優(yōu)化水庫運(yùn)行調(diào)度方案，從而達(dá)到提升水庫防洪能力[12]的目的。

1 研究區(qū)域

1.1 甬江流域

甬江流域位于全國海岸帶中段，地勢總體呈西南高、東北低。西南部為天臺山支脈四明山脈，總面積達(dá)2521km2，約占甬江流域總面積的43%；平原位于流域中部及北部沿海，姚江、奉化江、甬江“三江”橫貫其間。

甬江流域?qū)賮啛釒Ъ撅L(fēng)氣候區(qū)，四季分明，溫暖濕潤，光照充足，雨量豐沛。秋冬季，受北方冷高壓氣團(tuán)控制，以晴冷干燥天氣為主。春夏季，北方冷高壓氣團(tuán)逐漸衰退，雨水增多。

甬江流域主要水系包括甬江、同為甬江河源的奉化江和姚江，在寧波市區(qū)三江口匯合后稱為甬江。甬江流域具體包括了奉化江流域和姚江流域。

甬江流域基本位于寧波市境內(nèi)，是寧波市社會(huì)經(jīng)濟(jì)活動(dòng)最重要的區(qū)域。流域面積占寧波市總陸域面積的60%以上，人口和GDP均占全市的70%以上，是寧波市經(jīng)濟(jì)社會(huì)發(fā)展的核心區(qū)域。

1.2 流域水庫概況

近年來，隨著流域水利治理工作的不斷推進(jìn)，甬江流域總體構(gòu)建了“上蓄、中疏、下排”的防洪排澇工程體系，流域防洪排澇能力不斷增強(qiáng)[14]。“上蓄”水庫工程情況如下：

奉化江流域有4座大型水庫，1座中型水庫，分別為周公宅水庫、皎口水庫、亭下水庫、橫山水庫以及葛岙水庫(中型)。

姚江流域有1座大型水庫，4座中型水庫，分別為四明湖水庫(大型)、梁輝水庫、陸埠水庫、雙溪口水庫以及溪下水庫。

1.3 流域防洪排澇能力現(xiàn)狀

隨著流域防洪排澇工程體系建設(shè)步伐的不斷推進(jìn)，奉化江干支流及姚江干流洪水不再進(jìn)入平原，平原澇水外排能力明顯增強(qiáng)。

奉化江干支流上游山區(qū)相對集中且集雨面積較大?，F(xiàn)狀條件下，奉化江干支流防洪壓力仍然較大，奉化江干支流沿線最高洪水位普遍較堤防設(shè)計(jì)水位高0.03～0.19m，除東江、鄞江及縣江上游堤段以及奉化江鄞州大橋下游堤段防洪能力仍基本滿足20年一遇、50年一遇及100年一遇防洪標(biāo)準(zhǔn)要求外，奉化江“鄞州大橋上游”、東江“東江大閘下游”、鄞江“百梁橋下游”以及剡江干流沿線堤防安全超高不足，干流防洪安全缺乏有效保障。

2 案例研究分析

綜合考慮水文資料、地質(zhì)地貌資料以及水庫工程狀況等，本文選定甬江流域子流域奉化江流域上游橫山水庫作為研究對象并結(jié)合具體實(shí)例進(jìn)行分析研究。

2.1 研究思路

結(jié)合奉化江流域水文地質(zhì)狀況和橫山水庫運(yùn)行調(diào)度現(xiàn)狀，研究思路如下：

(1)了解橫山水庫的概況，明確橫山水庫現(xiàn)有的運(yùn)行調(diào)度方案。

(2)進(jìn)行甬江流域的設(shè)計(jì)暴雨和設(shè)計(jì)洪水計(jì)算，分析不同重現(xiàn)期下橫山水庫與其下游受保護(hù)斷面之間區(qū)域(非受控山區(qū))的洪水狀況。

(3)根據(jù)上述分析計(jì)算結(jié)果分析現(xiàn)有調(diào)度方案的作用和存在的問題，嘗試優(yōu)化現(xiàn)有調(diào)度方案并與之前調(diào)度方案的防洪能力進(jìn)行比較分析。

2.2 水文分析

2.2.1基礎(chǔ)資料

流域及區(qū)域水文分析計(jì)算共采用54處國家雨量站(包括甬江流域周邊的雨量站)。其中，流域內(nèi)雨量站點(diǎn)44處，附近雨量站點(diǎn)10處，各測站雨量資料起始年份和觀測年限不一。綜合各雨量站具體情況，選用1956—2016年共計(jì)61年的資料作為水文分析的基礎(chǔ)數(shù)據(jù)。全部61年雨量資料均通過“一致性”“可靠性”和“代表性”審查。

2.2.2設(shè)計(jì)暴雨

(1)設(shè)計(jì)暴雨計(jì)算

流域暴雨選樣采用年最大值法，具體步驟如下：①根據(jù)各雨量測站實(shí)測降雨資料，繪制雨量分布等值線；②采用網(wǎng)格法計(jì)算各研究流域長系列逐日面降雨量；③通過統(tǒng)計(jì)得出各時(shí)段的年最大面雨量；④進(jìn)行面雨量的頻率計(jì)算并繪制年最大雨量的頻率分布曲線；⑤通過各時(shí)段年最大雨量的頻率分布曲線查得不同頻率下設(shè)計(jì)雨量。

通過以上計(jì)算可以得到不同頻率下甬江流域、姚江流域、奉化江流域年最大1d、3d降雨量，成果見表1。

(2)設(shè)計(jì)暴雨成果的合理性檢查

甬江流域設(shè)計(jì)面雨量頻率參數(shù)與當(dāng)?shù)厮氖謨灾薪o定的值相當(dāng)，未出現(xiàn)明顯偏差。查證甬江流域現(xiàn)有的面雨量時(shí)間序列與設(shè)計(jì)暴雨計(jì)算中得到的設(shè)計(jì)結(jié)果大體相當(dāng)，沒有出現(xiàn)量級上、分布上相悖的情況，總體考慮設(shè)計(jì)暴雨的結(jié)果較為合理[15]。

2.2.3設(shè)計(jì)洪水

(1)產(chǎn)流計(jì)算

根據(jù)甬江流域及其子流域奉化江流域的歷史暴雨空間分布資料，本次選定“1963年12號臺風(fēng)”降雨作為奉化江流域典型降雨。

山區(qū)及平原均采用初損后損法進(jìn)行產(chǎn)流計(jì)算。產(chǎn)流以前的總損失水量稱為初損，以流域平均水深表示；后損主要是流域產(chǎn)流以后的下滲損失，以平均下滲率表示。一次降雨形成的徑流深計(jì)算方法如下式：

R=P-I0-ftR-P0

(1)

式中，P—一次降雨量，mm；I0—初損，mm；f—平均下滲率，mm/h；tR—產(chǎn)流歷時(shí)，h；P0—降雨后期不產(chǎn)流的雨量，mm。

具體成果見表2。

(2)匯流計(jì)算

匯流計(jì)算主要包括山區(qū)坡面匯流、山區(qū)河道匯流和平原河網(wǎng)匯流。當(dāng)集水面積大于50km2時(shí)，無資料地區(qū)坡面匯流計(jì)算采用“浙江省瞬時(shí)單位線法”，瞬時(shí)單位線是指流域上分布均勻、歷時(shí)趨于無窮小、強(qiáng)度趨于無窮大、總量為一個(gè)單位的地面凈雨在流域出口斷面形成的地面徑流過程線，其數(shù)學(xué)表達(dá)式為

表1 甬江流域設(shè)計(jì)暴雨成果

表2 奉化江流域不同重現(xiàn)期產(chǎn)流量計(jì)算成果 單位：萬m3

(2)

式中，n—線性水庫的個(gè)數(shù)；Γ(n)—n的伽瑪函數(shù)；K—線性水庫的調(diào)蓄系數(shù)，反映流域的匯流時(shí)間。按照這兩參數(shù)的匯流模式，只要求出n和K兩參數(shù)，即可求出瞬時(shí)單位線。

在實(shí)際應(yīng)用中，一般用S曲線轉(zhuǎn)化成時(shí)段為Δt的無因次單位線U(Δt，t)，即

(3)

U(Δt，t)=S(t)-S(t-Δt)

(4)

用瞬時(shí)單位線計(jì)算設(shè)計(jì)洪水時(shí)，應(yīng)考慮單位線對雨強(qiáng)的非線性反應(yīng)，采用不同時(shí)段凈雨強(qiáng)度條件下使用不同單位線的方法，即變雨強(qiáng)單位線法。

當(dāng)集水面積小于50km2時(shí)，采用“浙江省推理公式法”，即

Qmp=0.278×ip×F

(5)

(6)

式中，F(xiàn)—水庫壩址以上集雨面積，km2；iP—凈雨強(qiáng)度，mm/h；hR—τ時(shí)段內(nèi)凈雨，mm；τ—集流時(shí)間，h；n—暴雨指數(shù)。

計(jì)算得出橫山水庫非受控山區(qū)(橫山—興奉橋)各重現(xiàn)期的洪峰流量，以及按現(xiàn)有調(diào)度方案下的橫山水庫相應(yīng)出流量，詳見表3。

2.2.4結(jié)果分析

綜合2.2.3節(jié)的計(jì)算成果和2.3.2節(jié)調(diào)度方案可以發(fā)現(xiàn)，在流域20～50年一遇的暴雨情況下，因非受控山區(qū)洪峰未超過505.0、625.0m3/s，

表3 橫山水庫非受控山區(qū)不同重現(xiàn)期洪峰流量及水庫相應(yīng)出流 單位：m3/s

造成下游斷面發(fā)生洪峰時(shí)，橫山水庫均存在相應(yīng)出流，水庫出流會(huì)增加干流洪峰、抬升上下游沿線水位，對下游河道防洪造成一定不利影響。

2.3 水庫概況與運(yùn)行調(diào)度方案

2.3.1水庫概況

橫山水庫位于奉化江支流縣江上游。水庫壩址以上集水面積為150.8km2，總庫容為1.108億m3，是一座以防洪、供水、灌溉為主，并具備發(fā)電、養(yǎng)殖等綜合性功能的大型水利工程。水庫大壩按百年一遇防洪標(biāo)準(zhǔn)設(shè)計(jì)，萬年一遇防洪標(biāo)準(zhǔn)校核。

橫山水庫調(diào)度控制斷面興奉橋，位于奉化主城區(qū)，興奉橋與橫山壩址的距離約為14.5km。興奉橋斷面與橫山水庫之間存在74km2的產(chǎn)流區(qū)域，區(qū)間來水不受水庫控制，又稱為非受控區(qū)間，如圖2所示。

圖2 橫山水庫山區(qū)產(chǎn)流分區(qū)

2.3.2水庫運(yùn)行調(diào)度規(guī)則

興奉橋堤防防洪標(biāo)準(zhǔn)為50年一遇，20年一遇安全泄量為505.0m3/s，50年一遇安全泄量為625.0m3/s。橫山水庫正常按興奉橋斷面安全泄量進(jìn)行補(bǔ)償調(diào)節(jié)，具體調(diào)度方式如下：

①水庫臺汛汛限水位105.17m。

②當(dāng)庫水位105.17m

③當(dāng)114.77m

④當(dāng)Z>118.17m時(shí)，泄洪洞及泄洪閘全開敞泄，控制下泄流量不超過入庫流量。

2.4 水庫防洪能力現(xiàn)狀及提升分析

2.4.1水庫防洪能力現(xiàn)狀

根據(jù)現(xiàn)有的橫山水庫運(yùn)行調(diào)度方案，當(dāng)奉化江流域遭遇20～50年一遇的設(shè)計(jì)暴雨時(shí)對橫山水庫進(jìn)行調(diào)洪計(jì)算并分析得出水庫削峰與蓄洪效益。

當(dāng)橫山水庫遭遇流域20～50年一遇設(shè)計(jì)暴雨時(shí)，水庫削峰、蓄洪效果較好，削峰率達(dá)53.9%～55.3%，攔蓄洪率達(dá)48.5%～57.8%；且因水庫前期按505.0、625.0m3/s補(bǔ)償調(diào)度，預(yù)留一定庫容應(yīng)對中后期大洪水，50年一遇水庫調(diào)洪效果優(yōu)于20年一遇的調(diào)洪效果。

2.4.2水庫運(yùn)行規(guī)則優(yōu)化

根據(jù)2.3.4節(jié)和2.4.1節(jié)中的結(jié)果分析可以發(fā)現(xiàn)，橫山水庫現(xiàn)有的調(diào)度方案可以較好地滿足削峰蓄洪的任務(wù)，但是依舊存在值得優(yōu)化的地方。根據(jù)橫山水庫現(xiàn)有的運(yùn)行調(diào)度方案以及產(chǎn)流域的水文分析結(jié)果，考慮提高下游防洪能力，現(xiàn)擬定新的方案如下：

①水庫汛限水位105.17m。

②當(dāng)庫水位Z<117.77m(20年一遇)時(shí)，按興奉橋安全泄量210.0m3/s進(jìn)行補(bǔ)償調(diào)節(jié)。

③當(dāng)117.77m

④當(dāng)118.11m

⑤當(dāng)Z>118.97m時(shí)，泄洪閘全開，控制下泄流量不超過入庫流量。

2.4.3結(jié)果分析

計(jì)算奉化江流域遭遇20～50年一遇設(shè)計(jì)暴雨下，按新的調(diào)度方案(即按210.0m3/s補(bǔ)償調(diào)度)橫山水庫的調(diào)洪演算結(jié)果，詳見表4。

由表4的計(jì)算結(jié)果可以發(fā)現(xiàn)：

①當(dāng)遭遇流域20年一遇設(shè)計(jì)暴雨時(shí)，若橫山水庫前期按興奉橋210.0m3/s補(bǔ)償調(diào)度，水庫削峰率達(dá)80.5%，攔蓄洪率達(dá)72.3%，較按505.0、625.0m3/s補(bǔ)償調(diào)度(削峰率53.9%，攔洪率48.5%)的削峰、蓄洪效果提升明顯。

②當(dāng)遭遇流域50年一遇設(shè)計(jì)暴雨時(shí)，若橫山水庫前期按興奉橋210.0m3/s補(bǔ)償調(diào)度，中后期調(diào)蓄能力有所減弱，削峰率僅45.7%，較按505.0、625.0m3/s補(bǔ)償調(diào)度(削峰率55.3%)的削峰效果有所降低，但因前期調(diào)蓄增強(qiáng)，攔蓄洪率達(dá)65.5%，較按505.0、625.0m3/s補(bǔ)償調(diào)度(攔蓄洪率57.8%)蓄洪效果有所提升。

表4 橫山水庫削峰與蓄洪效益統(tǒng)計(jì)結(jié)果

3 結(jié)語

本文從調(diào)整水庫補(bǔ)償調(diào)度流量角度出發(fā)，優(yōu)化奉化江流域上游橫山水庫的防洪調(diào)度方案，有效增強(qiáng)了水庫削峰、攔洪蓄洪能力，緩解下游河道的防洪壓力，水庫防洪能力得到明顯提升。

限于篇幅和基礎(chǔ)資料情況，本文研究在內(nèi)容的深度和廣度上存在不足。事實(shí)上，除去水庫補(bǔ)償調(diào)度流量這一影響因素，從擴(kuò)充庫容、調(diào)整預(yù)泄量、更改平均起調(diào)水位等角度出發(fā)，分析相應(yīng)的水庫優(yōu)化調(diào)度方案，提升水庫防洪能力也具有很好的研究和實(shí)用價(jià)值。