李曉新

（1.山西省防汛抗旱指揮部辦公室，山西 太原 030002；2.太原理工大學(xué)，山西 太原 030024）

1 引言

漳澤水庫承擔(dān)著漳澤電廠、長治鋼鐵廠等十余家大中型企業(yè)用水以及1.32萬hm2耕地的灌溉用水。設(shè)計年供水量1億m3，目前漳澤水庫實際工業(yè)供水量達(dá)4 691萬m3，隨著社會經(jīng)濟的發(fā)展，需水量將逐步增加并接近水庫設(shè)計供水量，批復(fù)的工業(yè)供水指標(biāo)達(dá)7 070萬m3。與此同時，自20世紀(jì)80年代以來，漳澤水庫平均入庫徑流量明顯減少，平均徑流量僅7 950萬m3，水資源供需矛盾比較突出，如果遭遇連續(xù)枯水年，水庫的供水形勢將更加嚴(yán)峻。因此，在水庫防汛與興利矛盾日漸突出的形勢下，除采取工程措施外，有必要研究經(jīng)濟的非工程措施，即改變傳統(tǒng)調(diào)度方式，充分發(fā)揮水庫的調(diào)蓄潛力。以漳澤水庫目前的汛期汛限水位靜態(tài)分期控制方式進(jìn)行計算，26年來，汛期多年平均棄水量為4 200萬m3，占來水量的26.4%。由此可見，現(xiàn)行的汛限水位控制方案棄水量較大，洪水資源仍有利用空間。為更好地利用汛期洪水資源來提高水庫興利效益，在現(xiàn)有基礎(chǔ)上進(jìn)行漳澤水庫汛限水位動態(tài)控制研究是十分必要的。同時漳澤水庫大壩壩體穩(wěn)定、封閉條件及工程條件良好；區(qū)域降水量、洪水具有顯著的季節(jié)性特征，開展汛限水位動態(tài)控制研究是可行的。

2 漳澤水庫汛限水位動態(tài)控制域分析

2.1 降雨預(yù)報信息的可利用性分析

經(jīng)過多年的研究和實踐，汛限水位控制相關(guān)理論已有很大發(fā)展。汛限水位的控制方法主要經(jīng)歷了三個階段[1，2]：單一固定水位控制、分期固定水位控制及分期動態(tài)水位控制。水庫汛限水位的單一控制只利用了洪水的統(tǒng)計信息；分期固定水位控制利用了暴雨、洪水季節(jié)變化規(guī)律，改變了整個汛限水位單一固定不變的控制方式；水庫汛限水位動態(tài)控制利用氣象及水文預(yù)報的信息和當(dāng)前洪水信息，在有效、可靠的預(yù)見期內(nèi)和即時洪水情勢條件下，適時調(diào)整水庫洪水起調(diào)水位（汛限水位），確保在大洪水來臨前將水庫水位安全降至原設(shè)計汛限水位，且保證在發(fā)生設(shè)計或校核標(biāo)準(zhǔn)洪水時水庫大壩安全不受影響，達(dá)到在不降低水庫防洪標(biāo)準(zhǔn)的前提下增加水庫效益的目的。水庫汛限水位動態(tài)控制方法在一定程度上解決了汛期防洪與興利的矛盾。該方法利用洪水預(yù)報信息或降雨預(yù)報信息，根據(jù)具體情況實時將水庫水位控制在動態(tài)控制域范圍內(nèi)。相對于以往設(shè)定一個汛限水位而言，“控制域”的提出更為科學(xué)合理。針對漳澤水庫的實際情況，由于其水庫洪水預(yù)報方案精度低、預(yù)報結(jié)果不準(zhǔn)確，在對水庫相關(guān)流域降雨預(yù)報信息可利用性進(jìn)行分析的基礎(chǔ)上，采用預(yù)泄能力約束法、考慮與設(shè)計條件差別的方法綜合確定汛限水位動態(tài)控制域。

根據(jù)目前的氣象預(yù)報發(fā)展水平與預(yù)報精度，對長治氣象臺未來24 h和未來24～48 h的短期降雨預(yù)報的精度、預(yù)報誤差分布規(guī)律進(jìn)行分析發(fā)現(xiàn)：

第一，長治市氣象臺未來24 h和未來24～48 h各預(yù)報降雨量對應(yīng)的實際降雨量的經(jīng)驗頻率分布屬于偏態(tài)型，空報率大于漏報率，即除了無雨預(yù)報以外，其他量級有“預(yù)報偏大”的傾向。

第二，長治市氣象臺預(yù)報未來24 h和未來24～48 h為無雨、小雨和中雨時，發(fā)生大雨及以上量級的頻率較低。

第三，根據(jù)現(xiàn)有資料，在大雨以上量級降雨預(yù)報的發(fā)布次數(shù)上，統(tǒng)計得到的大雨及以上量級降雨預(yù)報的頻率、空報率和漏報率波動較大。

為更好地將降雨預(yù)報用于水庫調(diào)度中，并使汛限水位動態(tài)控制建立在安全可靠的基礎(chǔ)上，在接到不同量級降雨預(yù)報信息后能夠判斷是否應(yīng)該加大泄量，采用頻率分析法進(jìn)一步分析某一量級降雨預(yù)報條件下實際可能發(fā)生降雨量的頻率分布規(guī)律，給決策者提供“外延”信息。

表1 非參數(shù)檢驗成果表

經(jīng)非參數(shù)假設(shè)檢驗（表1），長治市氣象臺無雨、小雨和中雨預(yù)報對應(yīng)的實際降雨量分布屬于偏態(tài)型，符合PIII型分布。以各級降雨預(yù)報對應(yīng)的實際降雨量為樣本，利用矩法估計其統(tǒng)計參數(shù)（式1—2）初估參數(shù)x和CV，初步選定CS（式3）與CV的比值，然后以PIII型曲線作為理論線型，目估適線法確定統(tǒng)計參數(shù)，由于大雨以上量級樣本容量太小，未對其進(jìn)行頻率計算，只對未來24 h和未來24～48 h無雨、小雨和中雨預(yù)報相應(yīng)的實際降雨量進(jìn)行頻率計算。計算結(jié)果表明：未來24 h預(yù)報無雨、小雨和中雨時，發(fā)生大暴雨的頻率均低于0.01%；未來24～48 h預(yù)報中雨時，發(fā)生大暴雨的頻率約為0.2%～0.1%。當(dāng)預(yù)測未來24 h、未來24～48 h預(yù)報無雨、小雨和中雨時，發(fā)生大暴雨的頻率均非常小。因此，長治市氣象臺未來24 h和未來24～48 h無雨、小雨和中雨預(yù)報在漳澤水庫汛限水位動態(tài)控制中可以應(yīng)用。

2.2 汛限水位動態(tài)控制域上限的確定

2.2.1 確定步驟

汛限水位動態(tài)控制域上限的確定分4個步驟：

第一，計算有效預(yù)見期t1。計算公式為：

式中：tf——短期降雨預(yù)報可利用的預(yù)見期；

tc——預(yù)報信息傳遞時間；

tj——預(yù)泄決策調(diào)令傳播時間；

tz——閘門操作時間。

第二，求有效預(yù)見期t1期間平均出流量q出。假定經(jīng)t1時期后浮動水位增值ΔZ1，則絕對水位為Zt1=Z0+ΔZ1，根據(jù)Zt1與Z0，查水位—泄量關(guān)系得出qt1和q0，近似求出t1期間平均出流q出＝（qt1+q0）/2。

第三，求Q入。如有可利用的洪水預(yù)報信息，則根據(jù)面臨時刻t0預(yù)報的入庫洪水過程推求t1期間平均入庫流量，用梯形法或矩形法可求出Q入；若無洪水預(yù)報信息，利用降雨預(yù)報信息，則t1期間平均入庫流量則為發(fā)布大雨以上量級降雨預(yù)報信息后，有效預(yù)見期內(nèi)的平均入庫流量。

第四，汛限水位動態(tài)控制域上限計算。

預(yù)泄水量（q出－Q入）×t1＝ΔV1

上浮后的庫容Vt1＝V0+ΔV1。

其中，V0為面臨時刻t0水位Z0所相應(yīng)庫容。由Vt1查求庫容水位關(guān)系得到Z′t1，若Z′t1≈Zt1則試算停止，否則重設(shè)Zt1，重復(fù)第二步至第四步步驟，直到滿足要求為止。Zt1即為允許上浮的汛限水位。

2.2.2 確定方法

第一，采用預(yù)泄能力約束法來確定漳澤水庫汛限水位動態(tài)控制域上限。

漳澤水庫預(yù)泄能力主要受四方面約束：一是受水庫下游安全泄量300 m3/s的約束。二是受預(yù)泄時間，即有效預(yù)見期的約束。根據(jù)短期降雨預(yù)報可利用性分析成果，未來48 h短期預(yù)報信息在實際調(diào)度中是可利用的，假定實際降雨發(fā)生在最不利時刻，即24～48 h降雨預(yù)報的24 h之初，則降雨預(yù)報可利用的預(yù)見期為24 h?？紤]到信息傳遞時間、預(yù)泄決策調(diào)令傳達(dá)及實施時間、閘門操作時間等在實際調(diào)度中可能有所差異，有效預(yù)見期取20 h或22 h。三是受預(yù)泄時的水庫入庫流量約束。此次漳澤水庫汛限水位動態(tài)控制分析采用的是降雨預(yù)報信息，即水庫相關(guān)流域氣象臺發(fā)布未來48 h大雨、暴雨、大暴雨降雨預(yù)報信息后，水庫開始預(yù)泄。因此，預(yù)泄時入庫流量即為發(fā)布24～48 h大雨及以上量級降雨預(yù)報信息后的入庫流量。鑒于安全考慮，在確定漳澤水庫汛限水位控制域時，將有效預(yù)見期內(nèi)的最大入庫流量作為水庫預(yù)泄時的入庫流量。統(tǒng)計漳澤水庫流域2004—2009年發(fā)布24～48 h大雨以上量級降雨預(yù)報信息，有效預(yù)見期內(nèi)相應(yīng)實測入庫流量為1.1～28.1 m3/s，平均流量5.4 m3/s，最大入庫流量28.1 m3/s，為安全起見，取30 m3/s。四是受城市供水與灌溉平均出流約束。根據(jù)漳澤水庫設(shè)計年供水量，有效預(yù)見期內(nèi)水庫出流取5 m3/s。

基于上述基本條件分析，可得出如下結(jié)論：

一是若有效預(yù)見期取22h，則可預(yù)泄水量2178萬m3。對主汛期而言，在原汛限水位（901.0 m）對應(yīng)的庫容基礎(chǔ)上增加2178萬m3，此時對應(yīng)的水位為901.8m，即預(yù)泄能力約束方法確定的主汛期水庫汛限水位動態(tài)控制域上限；對前后汛期而言，在水庫現(xiàn)行分期汛限水位901.5 m對應(yīng)的庫容基礎(chǔ)上增加2 178萬m3，此時對應(yīng)的水位為902.3 m，即預(yù)泄能力約束方法確定的前后汛期水庫汛限水位動態(tài)控制域上限。

二是若有效預(yù)見期取20 h，則可預(yù)泄水量1 980萬m3。對主汛期而言，在原汛限水位（901.0 m）對應(yīng)的庫容基礎(chǔ)上增加1 980萬m3，此時對應(yīng)的水位為901.75 m，即預(yù)泄能力約束方法確定的主汛期水庫汛限水位動態(tài)控制域上限；對前后汛期而言，在水庫現(xiàn)行分期汛限水位901.5 m對應(yīng)的庫容基礎(chǔ)上增加1 980萬m3，此時對應(yīng)的水位為902.2 m，即預(yù)泄能力約束方法確定的前后汛期水庫汛限水位動態(tài)控制域上限。綜上所述，漳澤水庫以預(yù)泄能力為基礎(chǔ)的汛限水位動態(tài)控制域上限計算成果見表2。

表2 漳澤水庫預(yù)泄能力方法下汛限水位動態(tài)控制域

第二，考慮與設(shè)計條件差別的方法確定漳澤水庫汛限水位動態(tài)控制域上限。

漳澤水庫上游流域內(nèi)各支流有6座中型水庫，根據(jù)近年來歷史資料及現(xiàn)場實地調(diào)研考查，部分水庫現(xiàn)狀與《1984年漳澤水庫改建工程初步設(shè)計報告》中得出汛限水位901.0 m的水庫狀況有較大差別。其中莊頭水庫流域內(nèi)大部分為石灰?guī)r地區(qū)，植被較差，水土流失嚴(yán)重，完全處于空庫狀態(tài)；陶清河水庫為空庫度汛；西堡水庫常年無水。根據(jù)1984年前后上述3座水庫入庫徑流對比，假定發(fā)生較大量級洪水時3座水庫處于死水位狀態(tài)，經(jīng)計算，3座水庫與設(shè)計條件相比，可多截留的水量為3 353萬m3，相當(dāng)于漳澤水庫入庫洪水量減少了3 353萬m3，以此作為水庫汛限水位上浮原設(shè)計汛限水位多出的庫容值。據(jù)此計算，對主汛期而言，在原汛限水位901.0 m對應(yīng)的庫容基礎(chǔ)上增加了3 353萬m3，此時對應(yīng)的水位為902.2 m。對前后汛期而言，對應(yīng)的水位上限為902.4 m。即考慮與設(shè)計條件差別方法確定的前后汛期水庫汛限水位動態(tài)控制域上限。

3 結(jié)語

對漳澤水庫流域未來24 h和未來24～48 h短期降雨預(yù)報精度進(jìn)行分析，表明未來24 h和未來24～48 h無雨、小雨、中雨預(yù)報信息均可用于漳澤水庫汛限水位動態(tài)控制。采用預(yù)泄能力約束法、考慮與設(shè)計條件差別的方法，得出了三種汛限水位動態(tài)控制域方案。在保障安全的前提下、在有效預(yù)見期內(nèi)，須將水庫水位降至原設(shè)計汛限水位及以下，確保水庫大壩安全且汛限水位逐步調(diào)整為原則，并考慮漳澤水庫所處流域的特點及調(diào)度經(jīng)驗，結(jié)合水庫汛期分期研究成果，確定漳澤水庫汛限水位動態(tài)控制域為主汛期901.0～901.75 m，前后汛期901.5～902.2 m。