多風(fēng)險(xiǎn)因素影響下的水庫防洪調(diào)度風(fēng)險(xiǎn)綜合評(píng)估研究

管新建，胡 棟，孟 鈺

(鄭州大學(xué) 水利與環(huán)境學(xué)院，鄭州 450001)

0 引 言

水庫是重要的防洪興利工程，合理分配水庫防洪與興利庫容，是提高水資源利用率的重要手段。在水庫防洪調(diào)度過程中存在較多水文、水力等不確定性因素以及調(diào)控汛限水位可能導(dǎo)致不同程度洪災(zāi)風(fēng)險(xiǎn)。因此，如何全面考慮防洪風(fēng)險(xiǎn)因素，綜合衡量調(diào)度風(fēng)險(xiǎn)，合理調(diào)控汛限水位，協(xié)調(diào)水庫防洪與興利關(guān)系，對(duì)水庫防洪調(diào)度風(fēng)險(xiǎn)研究具有重要意義。

關(guān)于水庫防洪調(diào)度風(fēng)險(xiǎn)方面的研究，國外研究主要集中在風(fēng)險(xiǎn)損失評(píng)估方面。早期，為評(píng)估汛期水庫調(diào)整蓄水和棄水帶來的風(fēng)險(xiǎn)，在設(shè)計(jì)水庫汛期限制水位時(shí)開始引入風(fēng)險(xiǎn)理論[1]。國外對(duì)于洪水風(fēng)險(xiǎn)影響的評(píng)估主要采用由洪水造成的經(jīng)濟(jì)損失量化值進(jìn)行表示，設(shè)置多種汛限水位方案，對(duì)比不同汛限水位所帶來的損失或收益，優(yōu)選出合理的汛限水位值[2]。此外，部分研究將潰壩風(fēng)險(xiǎn)損失與水庫大壩風(fēng)險(xiǎn)率相結(jié)合[3-5]。國內(nèi)研究主要集中在水庫防洪風(fēng)險(xiǎn)因素模擬與風(fēng)險(xiǎn)率量化方面。目前，國內(nèi)學(xué)者通過分析水庫防洪調(diào)度中相關(guān)水文、水力等不確定性因素，總結(jié)了水庫風(fēng)險(xiǎn)分析的對(duì)象、目標(biāo)、及常見的幾類風(fēng)險(xiǎn)因素[6-7]，主要采用隨機(jī)算法對(duì)風(fēng)險(xiǎn)因素進(jìn)行模擬，常用方法包括蒙特卡洛模擬法[8-10]與貝葉斯理論等[11-13]，并采用水庫、庫區(qū)及下游多類風(fēng)險(xiǎn)率指標(biāo)對(duì)風(fēng)險(xiǎn)模擬結(jié)果進(jìn)行評(píng)估[10,17]。隨后，基于Copula函數(shù)的洪水過程多類水文要素的聯(lián)合概率分布研究逐漸受到重視與推廣[14,15]。可見，國內(nèi)在水庫防洪調(diào)度風(fēng)險(xiǎn)識(shí)別與模擬方法上的研究較為成熟，但對(duì)風(fēng)險(xiǎn)結(jié)果的評(píng)估主要還是以單風(fēng)險(xiǎn)率指標(biāo)為主，也有部分研究提出了多目標(biāo)風(fēng)險(xiǎn)評(píng)估指標(biāo)體系[8]，但基于綜合模型對(duì)風(fēng)險(xiǎn)綜合評(píng)估的研究仍然相對(duì)較少。

本文針對(duì)水庫汛限水位調(diào)控及洪水資源化問題，提出水庫防洪調(diào)度風(fēng)險(xiǎn)綜合評(píng)估總體框架，確定水庫防洪調(diào)度中涉及的四種風(fēng)險(xiǎn)因素，采用蒙特卡羅模擬技術(shù)進(jìn)行隨機(jī)模擬；應(yīng)用水庫調(diào)洪演算模擬模型獲取調(diào)洪演算過程與結(jié)果；通過分析水庫防洪調(diào)度過程中可能出現(xiàn)的風(fēng)險(xiǎn)情況，構(gòu)建風(fēng)險(xiǎn)綜合評(píng)價(jià)指標(biāo)體系；采用層次分析法對(duì)水庫前、后汛期不同汛限水位方案進(jìn)行綜合評(píng)價(jià)。以陸渾水庫為例，對(duì)陸渾水庫進(jìn)行了防洪調(diào)度風(fēng)險(xiǎn)綜合分析與評(píng)價(jià)，為水庫防洪調(diào)度風(fēng)險(xiǎn)綜合評(píng)估研究提供了技術(shù)支撐，為陸渾水庫實(shí)際防洪調(diào)度風(fēng)險(xiǎn)研究提供了科學(xué)參考。

1 水庫防洪調(diào)度風(fēng)險(xiǎn)綜合評(píng)估總體框架

水庫防洪調(diào)度風(fēng)險(xiǎn)綜合評(píng)估總體框架包括三個(gè)部分：①風(fēng)險(xiǎn)因素模擬。識(shí)別影響水庫防洪調(diào)度的風(fēng)險(xiǎn)因素，采用蒙特卡羅模擬技術(shù)對(duì)各風(fēng)險(xiǎn)因素進(jìn)行隨機(jī)模擬；②水庫調(diào)洪演算模擬模型。依據(jù)水庫調(diào)洪演算基本理論，建立水庫調(diào)洪演算模擬模型；③風(fēng)險(xiǎn)綜合評(píng)估模型。通過分析水庫防洪調(diào)度過程中可能出現(xiàn)的風(fēng)險(xiǎn)情況，構(gòu)建風(fēng)險(xiǎn)綜合評(píng)價(jià)指標(biāo)體系；采用層次分析法對(duì)多風(fēng)險(xiǎn)因素影響下的水庫風(fēng)險(xiǎn)進(jìn)行綜合評(píng)估。總體框架流程如圖1所示。

圖1 水庫防洪調(diào)度風(fēng)險(xiǎn)綜合評(píng)估總體框架Fig.1 The framework of comprehensive evaluation of reservoir flood control scheduling

2 基于蒙特卡羅技術(shù)的水庫防洪調(diào)度風(fēng)險(xiǎn)因素模擬

通過分析水庫在防洪調(diào)度過程中可能出現(xiàn)的風(fēng)險(xiǎn)情況，識(shí)別出四種主要的風(fēng)險(xiǎn)因素，分別是設(shè)計(jì)洪水過程、洪水預(yù)報(bào)誤差、水庫調(diào)度滯時(shí)、水位-庫容關(guān)系。根據(jù)各風(fēng)險(xiǎn)因素特點(diǎn)確定其概率分布函數(shù)，選用蒙特卡羅技術(shù)對(duì)各風(fēng)險(xiǎn)因素進(jìn)行隨機(jī)模擬。

2.1 設(shè)計(jì)洪水過程的推求

(1)典型洪水選擇的原則：①從實(shí)測(cè)洪水資料中選取歷時(shí)不小于5 d的洪水過程；②選取洪水時(shí)要求一日洪量包含洪峰，洪量長(zhǎng)時(shí)段包含短時(shí)段的原則。

(2)采用同頻率放大法推求設(shè)計(jì)洪水過程線[16]：①由實(shí)測(cè)洪水資料利用P-Ⅲ曲線可得出設(shè)計(jì)洪峰和不同時(shí)段的設(shè)計(jì)洪量；②根據(jù)典型洪水過程計(jì)算其洪峰及最大洪量，將典型洪水過程按照洪峰、最大一日洪量、最大三日洪量的不同倍比分別進(jìn)行放大計(jì)算。

2.2 預(yù)報(bào)誤差隨機(jī)模擬

(1)

式中：ξt為標(biāo)準(zhǔn)正態(tài)分布變量；rt為0～1的均勻隨機(jī)數(shù)。

(2)

(3)

(4)

(5)

采用式(6)計(jì)算可得對(duì)數(shù)正態(tài)分布隨機(jī)序列：

xt=a+exp(σyξt+μy) (t=1,2,…)

(6)

2.3 調(diào)度滯時(shí)隨機(jī)模擬

調(diào)度滯時(shí)是由于洪水調(diào)度決策及調(diào)度方案實(shí)施程序中出現(xiàn)不確定性等情況，從而導(dǎo)致實(shí)際調(diào)度實(shí)施時(shí)間滯后的時(shí)間。調(diào)度滯時(shí)風(fēng)險(xiǎn)因素通常采用三角分布來進(jìn)行描述，三角分布涉及三個(gè)參數(shù)：調(diào)度滯時(shí)可能的最大值c、最可能值b和最小值a，計(jì)算公式如式(7)所示：

(7)

式中，rt為0～1的均勻隨機(jī)數(shù)；xt為三角分布隨機(jī)數(shù)。

2.4 水位-庫容關(guān)系隨機(jī)模擬

由于使用的測(cè)量器材不同，以及在發(fā)生洪水時(shí)產(chǎn)生的沖擊程度不同等原因，導(dǎo)致水位-庫容關(guān)系的發(fā)生變化[17]。根據(jù)以往經(jīng)驗(yàn)統(tǒng)計(jì)，水位-庫容函數(shù)V(Z)通常符合正態(tài)分布，可按公式(8)計(jì)算得到V(Z)：

V(Z)=X(Z)+X(Z)×0.1×U

(8)

式中：V(Z)為實(shí)際水位-庫容函數(shù)；X(Z)為設(shè)計(jì)水位-庫容函數(shù)；Z為水位；U為標(biāo)準(zhǔn)正態(tài)分布隨機(jī)數(shù)。

3 水庫調(diào)洪演算模擬模型

水庫調(diào)洪模擬演算基本原則是依據(jù)水庫調(diào)洪演算基本原理，結(jié)合實(shí)際水庫防洪調(diào)度規(guī)則對(duì)入庫洪水進(jìn)行調(diào)洪模擬計(jì)算。水庫調(diào)洪演算基本原理主要包括水量平衡原理和動(dòng)力平衡原理[16]。水量平衡原理可用圣維南方程組表示，但通常難以得出精確的分析解，所以假定水庫水位與庫容在ΔT時(shí)段內(nèi)呈線性變化，水量平衡方程可簡(jiǎn)化為式(9)：

(9)

式中：Q1、Q2分別代表時(shí)段初、末的入庫流量；q1、q2分別代表時(shí)段初、末的出庫流量；V1、V2分別代表時(shí)段初、末的水庫蓄水量；ΔT為時(shí)段長(zhǎng)。

動(dòng)力平衡原理可由水庫蓄泄平衡方程或蓄泄曲線來表示[16]。其中，水庫泄水量 與壩前水位q有關(guān)，這種關(guān)系隨防洪調(diào)度過程中所采用的泄洪設(shè)施變化而變化。

本文所編制的水庫調(diào)洪演算模擬模型是依據(jù)水庫調(diào)洪演算基本原理，結(jié)合水庫實(shí)際防洪調(diào)度規(guī)則等條件，采用Visual Basic 6.0語言編寫模擬模型，包括輸入模塊、計(jì)算模塊、輸出模塊[18]。其中輸入模塊是實(shí)現(xiàn)程序基本數(shù)據(jù)的輸入，包括設(shè)置水庫初始起調(diào)水位，添加洪水入庫流量、水位-庫容關(guān)系、水位-泄流量關(guān)系等條件；計(jì)算模塊是依據(jù)水量平衡方程結(jié)合水庫防洪調(diào)度規(guī)則對(duì)洪水入庫流量與水庫泄流量關(guān)系進(jìn)行分析計(jì)算等。輸出模塊是輸出水庫模擬調(diào)洪演算的結(jié)果，包括洪水入庫時(shí)間、入庫流量、水庫下泄量、各時(shí)段水庫水位及庫容等。

4 水庫防洪調(diào)度風(fēng)險(xiǎn)綜合評(píng)估模型

4.1 風(fēng)險(xiǎn)綜合評(píng)價(jià)指標(biāo)體系

通過對(duì)水庫庫區(qū)與下游在分期洪水調(diào)度過程中可能的風(fēng)險(xiǎn)情況進(jìn)行分析與匯總，建立風(fēng)險(xiǎn)綜合評(píng)價(jià)指標(biāo)體系，其中分為兩類風(fēng)險(xiǎn)指標(biāo)：庫區(qū)風(fēng)險(xiǎn)控制指標(biāo)與下游風(fēng)險(xiǎn)控制指標(biāo)；前者包括水庫征地水位L1與水庫移民水位L2兩個(gè)風(fēng)險(xiǎn)指標(biāo)；后者包括最大安全泄量Q、平均淹沒人口P、平均淹沒耕地S三個(gè)風(fēng)險(xiǎn)指標(biāo)，具體的風(fēng)險(xiǎn)綜合評(píng)價(jià)指標(biāo)體系如圖2所示。

圖2 風(fēng)險(xiǎn)綜合評(píng)價(jià)指標(biāo)體系Fig.2 Comprehensive evaluation index system of flood risk

根據(jù)各風(fēng)險(xiǎn)指標(biāo)特點(diǎn)，確定其計(jì)算公式如下：

(1)庫區(qū)與下游風(fēng)險(xiǎn)超標(biāo)率。庫區(qū)風(fēng)險(xiǎn)指標(biāo)水庫征地水位L1、移民水位L2及下游風(fēng)險(xiǎn)指標(biāo)最大安全泄量Q均采用超標(biāo)率進(jìn)行估算，計(jì)算公式如(10)所示：

(10)

式中：R1、R2分別為水庫征地水位L1與移民水位L2的超標(biāo)率；R3為水庫最大安全泄量Q的超標(biāo)率；Mi為第i類水位或泄量的超標(biāo)次數(shù)；N為模擬場(chǎng)次數(shù)。

(2)平均淹沒人口。根據(jù)水庫實(shí)際歷史調(diào)查資料，采用同倍比方法估算[19]平均淹沒人口P，計(jì)算公式如(11)所示：

(11)

式中：P為平均淹沒人口；hi為單場(chǎng)次調(diào)洪最高庫水位；h1為水庫移民水位；h2為歷史實(shí)測(cè)情況下發(fā)生人口淹沒現(xiàn)象的實(shí)測(cè)水位；ph為歷史實(shí)測(cè)情況下發(fā)生人口淹沒現(xiàn)象的淹沒人口數(shù)；m為移民水位的超標(biāo)次數(shù)；N為模擬場(chǎng)次數(shù)。

(3)平均淹沒耕地。根據(jù)水庫實(shí)際歷史調(diào)查資料，采用同倍比方法估算[19]平均淹沒耕地S，計(jì)算公式如(12)所示：

(12)

式中：S為平均淹沒耕地；wi為單場(chǎng)次調(diào)洪最高庫水位；w1為水庫征地水位；w2為歷史實(shí)測(cè)情況下發(fā)生淹沒耕地現(xiàn)象的實(shí)測(cè)水位；qh為歷史實(shí)測(cè)情況下發(fā)生淹沒耕地現(xiàn)象的淹沒耕地面積；m為征地水位的超標(biāo)次數(shù)；N為模擬場(chǎng)次數(shù)。

4.2 綜合評(píng)價(jià)方法

層次分析法(AHP)是由Saaty提出的確定指標(biāo)相對(duì)權(quán)重的常用方法[20]。其基本原理是：首先根據(jù)評(píng)價(jià)者對(duì)各風(fēng)險(xiǎn)指標(biāo)的重要程度構(gòu)造判斷矩陣，然后對(duì)各判斷矩陣進(jìn)行一致性檢驗(yàn)；若通過一致性檢驗(yàn)，則求得各判斷矩陣的最大特征值及其對(duì)應(yīng)的特征向量即為各指標(biāo)權(quán)重值；否則，需要重新構(gòu)造判斷矩陣進(jìn)行計(jì)算。具體分為以下3個(gè)步驟：

(1)構(gòu)造判斷矩陣。按圖2所示風(fēng)險(xiǎn)綜合評(píng)價(jià)指標(biāo)體系結(jié)構(gòu)關(guān)系構(gòu)造判斷矩陣如式(13)所示：

(13)

式中：aij表示上一層指標(biāo)ui對(duì)uj的相對(duì)重要程度，i,j=1,2,…,n。

(2)一致性檢驗(yàn)。通過判斷矩陣偏離一致性C.I.與判斷矩陣隨機(jī)一致性R.I.的比值C.R.來檢驗(yàn)判斷矩陣的一致性。根據(jù)判斷矩陣按式(14)可計(jì)算其最大特征值λmax，然后按式(15)可計(jì)算偏離一致性比值C.I.：

|A-λE|=0

(14)

式中：A為待求的各判斷矩陣；λ為所求特征值；λmax為所求特征值λ中的最大值；

C.I.=(λmax-n)/(n-1)

(15)

通過查表可確定隨機(jī)一致性R.I.，按式(16)可計(jì)算 ：

C.R.=C.I./R.I.

(16)

當(dāng)C.R.<0.1時(shí)，認(rèn)為總排序滿足一致性條件，否則需重新構(gòu)造判斷矩陣。

(3)風(fēng)險(xiǎn)綜合值計(jì)算。按式(17)對(duì)各風(fēng)險(xiǎn)指標(biāo)計(jì)算值進(jìn)行規(guī)范化處理：

(17)

根據(jù)各指標(biāo)的權(quán)重wj和規(guī)范化值rij按式(18)計(jì)算其風(fēng)險(xiǎn)綜合值A(chǔ)i：

(18)

5 應(yīng)用實(shí)例

5.1 研究區(qū)基本情況

陸渾水庫位于黃河支流伊河上，控制流域面積3 492 km2，水庫壩頂高程333 m，設(shè)計(jì)總庫容約13.2 億m3，具體設(shè)計(jì)標(biāo)準(zhǔn)如表1所示。水庫前汛期7月1日-8月31日，汛限水位為317.0 m，相應(yīng)庫容5.7 億m3；后汛期9月1日-10月31日，汛限水位為317.5 m，相應(yīng)庫容5.9 億m3。該工程是以防洪為主，兼顧灌溉、發(fā)電的綜合利用工程，汛期時(shí)通過配合黃河下游三門峽、故縣水庫削減洪峰流量以減輕下游防洪壓力。目前，陸渾水庫汛限水位相對(duì)較低，可在汛期適當(dāng)調(diào)整以充分利用洪水資源，增加水庫興利效益。

表1 陸渾水庫設(shè)計(jì)洪水標(biāo)準(zhǔn)Tab.1 Design flood standards for Luhun Reservoir

5.2 主要風(fēng)險(xiǎn)因素模擬與調(diào)洪演算方案

陸渾水庫設(shè)計(jì)防洪運(yùn)用方式為：當(dāng)入庫流量未達(dá)到1 000 m3/s時(shí)，原則上按進(jìn)出庫平衡方式運(yùn)用；否則，按下泄流量1 000 m3/s控制運(yùn)用；當(dāng)水庫水位達(dá)20年一遇洪水位321.5 m時(shí)，若入庫流量未達(dá)到20年一遇洪水位相應(yīng)的泄流能力2 560 m3/s，原則上按進(jìn)出庫平衡方式運(yùn)用；否則，按敞泄運(yùn)用。目前，陸渾水庫前汛期汛限水位為317 m，后汛期汛限水位317.5 m，以現(xiàn)狀運(yùn)行水位為基礎(chǔ)，按照其前、后汛期時(shí)段的分期，設(shè)定如表2所示的不同汛限水位方案進(jìn)行水庫調(diào)洪模擬演算。

表2 陸渾水庫調(diào)洪演算方案Tab.2 Flood control scheme for LuHun reservoir

6 風(fēng)險(xiǎn)綜合評(píng)價(jià)結(jié)果與分析

6.1 指標(biāo)權(quán)重計(jì)算結(jié)果

基于4.1小節(jié)構(gòu)建的風(fēng)險(xiǎn)綜合評(píng)價(jià)指標(biāo)體系，采用層次分析法進(jìn)行指標(biāo)權(quán)重及綜合風(fēng)險(xiǎn)估計(jì)。根據(jù)各風(fēng)險(xiǎn)目標(biāo)遞階結(jié)構(gòu)，結(jié)合專家經(jīng)驗(yàn)，分層構(gòu)造與修正判斷矩陣，包括風(fēng)險(xiǎn)綜合評(píng)估總目標(biāo)Z、庫區(qū)風(fēng)險(xiǎn)X、下游風(fēng)險(xiǎn)Y判斷矩陣，對(duì)各判斷矩陣進(jìn)行最大特征值計(jì)算與一致性檢驗(yàn)，結(jié)果如表3所示。

根據(jù)上述過程計(jì)算并繪制各風(fēng)險(xiǎn)指標(biāo)權(quán)重值圖如圖3所示，由各風(fēng)險(xiǎn)指標(biāo)權(quán)重值可以看出，決策者在協(xié)調(diào)庫區(qū)風(fēng)險(xiǎn)和下游風(fēng)險(xiǎn)的情況下，更側(cè)重于庫區(qū)風(fēng)險(xiǎn)，權(quán)重值結(jié)果是合理的。

6.2 風(fēng)險(xiǎn)綜合評(píng)價(jià)結(jié)果

根據(jù)水庫防洪調(diào)度規(guī)則、水位-庫容關(guān)系、水位-泄流量關(guān)系等初始條件，選取計(jì)算時(shí)段ΔT=2 h，將初始數(shù)據(jù)輸入到水庫調(diào)洪模擬演算模型中，獲取前、后汛期水庫調(diào)洪演算模擬結(jié)果；根據(jù)各風(fēng)險(xiǎn)指標(biāo)權(quán)重值，計(jì)算得到風(fēng)險(xiǎn)綜合評(píng)價(jià)結(jié)果如表4、表5所示。

圖3 各風(fēng)險(xiǎn)指標(biāo)權(quán)重值圖Fig.3 Weights of risk indices

表4 陸渾水庫前汛期調(diào)洪演算和風(fēng)險(xiǎn)綜合評(píng)價(jià)結(jié)果Tab.4 Flood regulation routing and risk comprehensive evaluation of the Luhun Reservoir in pre flood season

表5 陸渾水庫后汛期調(diào)洪模擬演算和綜合評(píng)價(jià)結(jié)果Tab.5 Flood regulation routing and risk comprehensive evaluation of the Luhun Reservoir in post flood season

由表4可知，在水庫前汛期時(shí)，考慮4種風(fēng)險(xiǎn)因素影響下，水庫在遇5年一遇和20年一遇設(shè)計(jì)洪水時(shí)，在318.0、318.5、318.7 m三種汛限水位方案下的風(fēng)險(xiǎn)綜合值均為0；當(dāng)汛限水位提高至318.8m時(shí)，5年一遇設(shè)計(jì)洪水達(dá)到的最高庫水位超過征地水位319.5 m，其超標(biāo)率為1.94%，平均淹沒耕地0.29 hm2，綜合風(fēng)險(xiǎn)值為0.016；在汛限水位方案為319.0、319.3 m時(shí)，其超標(biāo)率與淹沒范圍更大，其超標(biāo)率分別為14.84%與37.74%，平均淹沒耕地分別為2.17與4.04 hm2，綜合風(fēng)險(xiǎn)值分別為0.191與1，并且在汛限水位方案為319.3 m時(shí)，出現(xiàn)淹沒人口，平均淹沒人口達(dá)到22.3人。對(duì)比6種汛限水位方案，汛限水位方案在318.7 m時(shí)為風(fēng)險(xiǎn)出現(xiàn)的臨界水位；目前，陸渾水庫前汛期汛限水位為317.0 m，相應(yīng)庫容5.7 億m3。根據(jù)洪水的設(shè)計(jì)標(biāo)準(zhǔn)，前汛期汛限水位可適當(dāng)抬高至318.7 m，興利庫容可增加0.64 m3，從而增加水庫興利效益。

由表5可知，在水庫后汛期時(shí)，考慮四種風(fēng)險(xiǎn)因素影響下，水庫在遇5年一遇和20年一遇設(shè)計(jì)洪水時(shí)，318.5、319.0、319.1、319.2 m四種汛限水位方案下的風(fēng)險(xiǎn)綜合值均為0；汛限水位方案為319.4 m時(shí)，水庫征地水位出現(xiàn)超標(biāo)，超標(biāo)率為4.84%，平均淹沒耕地0.75 hm2，綜合風(fēng)險(xiǎn)值為0.177；在汛限水位方案為319.5 m時(shí)，征地與移民水位、下游安全泄量均超標(biāo)，超標(biāo)率分別為8.22%、1.52%、2.02%，平均淹沒耕地?cái)U(kuò)大至1.03 hm2，平均淹沒人口增加至3.76人，綜合風(fēng)險(xiǎn)值為1。對(duì)比6種汛限水位方案，汛限水位方案在319.2 m時(shí)為風(fēng)險(xiǎn)出現(xiàn)的臨界水位。目前，陸渾水庫后汛期汛限水位317.5 m，根據(jù)洪水的設(shè)計(jì)標(biāo)準(zhǔn)，后汛期汛限水位可適當(dāng)抬高至319.2 m，興利庫容可增加0.66 m3，從而增加水庫興利效益。

7 結(jié) 語

水庫防洪調(diào)度過程中存在多種風(fēng)險(xiǎn)情況，綜合考慮水庫可能的風(fēng)險(xiǎn)情況，優(yōu)選汛限水位方案，在控制水庫庫區(qū)與下游整體調(diào)度風(fēng)險(xiǎn)的同時(shí)，能有效提高水庫的興利效益。本文以陸渾水庫為例，考慮設(shè)計(jì)洪水過程、洪水預(yù)報(bào)誤差、水庫調(diào)度滯時(shí)、水位-庫容關(guān)系四種風(fēng)險(xiǎn)因素，評(píng)估水庫不同汛限水位方案的風(fēng)險(xiǎn)綜合情況，結(jié)果表明：

(1)根據(jù)陸渾水庫設(shè)計(jì)標(biāo)準(zhǔn)，前汛期汛限水位可抬高至318.7 m，后汛期汛限水位可抬高至319.2 m，從而將洪水資源化，增加水庫興利效益。

(2)在水庫實(shí)際防洪調(diào)度過程中，預(yù)報(bào)預(yù)泄是更為有效的防洪方式。來水可能超過預(yù)報(bào)洪水，可通過對(duì)水庫下泄流量進(jìn)行實(shí)時(shí)調(diào)整，進(jìn)一步降低洪災(zāi)風(fēng)險(xiǎn)。

□