任建營張新建任 婕

（1.山西省漳澤水庫管理局 山西長治 046021；2.山西省防汛抗旱指揮部辦公室 山西太原 030002）

1 緣起

漳澤水庫位于山西省長治市北郊，是海河流域漳衛(wèi)南運河水系濁漳南源干流上的一座以工業(yè)、城市供水、防洪、灌溉為主兼顧水產(chǎn)養(yǎng)殖、旅游等綜合利用的大（2）型水利工程。流域面積3 176 km2，總庫容4.273億m3。

近十多年來，漳澤水庫管理局大力進(jìn)行水庫管理信息化建設(shè)，在水庫防洪調(diào)度和日常興利調(diào)度方面于省內(nèi)最早開發(fā)使用調(diào)洪計算軟件[1]，該軟件中設(shè)置了泄水洞和溢洪道閘門開度計算功能，用于水庫日常泄水的計算。隨著水庫管理水平的日益提高，對水庫防洪和興利調(diào)度計算精度方面也提出了更高的要求。漳澤水庫位于太行山上黨盆地，為平原水庫，大壩長2 514 m，正常高水位時水面面積31 km2，每厘米的水位差就有幾十萬立方米庫容，原有50 cm的水位數(shù)據(jù)級差難以滿足日益精細(xì)管理的需求。

2010年，水庫管理局以1 cm級差的水位—庫容—面積數(shù)據(jù)為基礎(chǔ)，升級改造水庫調(diào)洪計算軟件。在此過程中，發(fā)現(xiàn)并解決了泄水洞閘門流量計算精度偏低和溢洪道孔堰流邊界系數(shù)精度偏低帶來的問題。

2 泄水洞流量計算精度的提高

漳澤水庫泄水洞位于大壩壩軸樁號0+800處，為壩內(nèi)埋管式半有壓洞，洞首為高25.5 m的進(jìn)水塔，內(nèi)裝3.4m×3.7m（寬×高）平板鋼閘門，洞身斷面為3.2m×3.575 m（寬×高）的馬蹄型，設(shè)計限制最大下泄流量為90 m3/s。

在水庫泄水流量不大時，一般多通過泄水洞泄水。在庫水位不同時，同樣的泄水流量閘門開度是不同的。根據(jù)泄水洞流量計算公式①，已知上游水頭H和閘門開度e手工計算求泄水流量Q比較容易，但由于閘門開度e不能變換為以上游水頭H和泄水流量Q為自變量的顯函數(shù)形式，故已知上游水頭H和泄水流量Q不能直接求出閘門開度e。以往一般采用試算法或者查曲線圖的方法來求解閘門的開度，不僅計算速度慢，精度也不高。因此，開發(fā)計算機(jī)軟件用數(shù)值解法解決閘門開度計算的問題非常必要，它既提高了計算精度，又提高了工作效率。

該功能的使用方法是：以人機(jī)對話方式輸入庫水位和泄水流量，選擇泄水洞為泄水建筑物后，便可計算出泄水洞閘門開度。由于閘門開度是按泄水開始時的庫水位計算的，故不能考慮由于泄水引起庫水位下降后，對泄水流量的影響。對于較大流量、較長時間的泄水，應(yīng)在庫水位發(fā)生變化后，重新計算閘門開度。

2.1 泄水洞流量計算問題的發(fā)現(xiàn)

在閘門開度計算應(yīng)用實踐中，發(fā)現(xiàn)計算機(jī)軟件泄水洞流量計算存在較大誤差，例如：上游水頭8.7 m時，要求計算下泄流量為15 m3/s，通過軟件計算可得閘門開度為0.512 m，如果把閘門開度為0.512 m的計算結(jié)果帶入到公式①，得出泄水洞流量14.33 m3/s，計算精度為95.53%，誤差較大。

式中：Q——泄水洞流量，m3/s；

m——流量系數(shù)，無量綱；

b——閘寬，b=3.2 m；

e——閘門開度，m；

H——上游水頭，m，H為庫水位高程減泄水洞進(jìn)口底高程后之值；

ε——垂直收縮系數(shù)。

我們先是改進(jìn)計算方法來提高計算精度，但效果不明顯。

2.2 問題的解決

在泄水洞流量計算中，m和ε是重要的計算參數(shù)，但它們卻不是常數(shù)，而是隨著上游水頭H和閘門開度e的比值e/H而變化的量，這種變化規(guī)律無法用精確的數(shù)學(xué)公式表達(dá)，只是兩個數(shù)值表，計算者只能通過插值計算的方法來確定m和ε的取值。表1就是漳澤水庫泄水洞原設(shè)計提供的ε取值表。

表1 漳澤水庫泄水洞設(shè)計提供的ε取值表

通過仔細(xì)觀察發(fā)現(xiàn)，原來提供的公式計算參數(shù)中閘門垂直收縮系數(shù)ε表數(shù)值分布不均、精度偏低，小數(shù)點后的有效數(shù)字只有2位，這有可能是計算精度不高的原因，但一時我們也找不到更高精度的參數(shù)表來替代。

后來通過查閱《水力學(xué)計算手冊》[2]發(fā)現(xiàn)了泄水洞流量計算閘門垂直收縮系數(shù)參數(shù)表，于是采用了以下表2新的ε表格數(shù)據(jù)。

表2 泄水洞流量計算閘門垂直收縮系數(shù)ε參數(shù)表

表2比表1的參數(shù)分布均勻，閘門垂直收縮系數(shù)ε值提高了一個量級。通過閘門開度計算軟件測試，采用新的參數(shù)后，通過一系列數(shù)據(jù)測算，泄水洞流量計算精度有所提高，見表3。

表3 泄水洞閘門垂直收縮系數(shù)ε取值變化前后流量計算對比表

3 溢洪道相對開度精度的提高

溢洪道位于大壩右端，為正槽開敞式，閘室設(shè)胸墻，設(shè)4孔9.2 m×6.6 m（寬×高）的弧形鋼閘門，采用2×400 kN的油壓啟閉機(jī)。溢流堰為低駝峰堰，堰前鋪蓋高程比堰頂高程低2.0 m，胸墻底高程比堰頂高程高6.6 m。閘門高6.6 m，單孔寬9.2 m，4孔全寬36.8 m。出口為挑流消能，溢洪道全長304.0 m，最大流量2 100 m3/s。

3.1 問題的發(fā)現(xiàn)

溢洪道流量計算公式如下：

（1）當(dāng)開度e＜0.7H，為孔口出流，采用公式②

式中：Q——溢洪道流量，m3/s；

μ0——流量系數(shù)，μ0=0.685-0.19e/H；

e——閘門開度，m；

B——堰長，B＝9.2×4＝36.8 m；

H0——溢洪道前包括行進(jìn)流速在內(nèi)的堰頂水頭值，H0=H+（α·υ02）/（2g）,由于溢洪道閘前開闊，流速υ0很小，可忽略不計，即H0=H。

（2）當(dāng)開度e≥0.7H，為堰流，采用公式③

式中：Q——溢洪道流量，m3/s；

ε——側(cè)收縮系數(shù)，ε=0.9；

m——流量系數(shù)，當(dāng)P/H＞0.34時，m=0.452（P/H)-0.032；當(dāng)P/H≤0.34時，m=0.385+0.224×（P/H)0.934；

B——堰長，B=9.2×4＝36.8 m；

H——堰頂水頭，m；

P——堰高，P=2 m。

庫水位級差為1 cm計算溢洪道泄流能力時，發(fā)現(xiàn)溢洪道在堰頂水頭9.43 m處變?yōu)榭琢?，計算流量值比堰頂水頭9.42 m時減小32.7 m3/s，出現(xiàn)了在堰流變孔流時兩個公式計算結(jié)果銜接不好，得到的計算結(jié)果見表4。

結(jié)合該區(qū)間附近的數(shù)據(jù)，得到以下局部誤差放大圖。

3.2 問題的解決

在堰頂水頭9.4 m附近，正好是兩個計算公式的銜接之處，即堰流轉(zhuǎn)變?yōu)榭琢鞯霓D(zhuǎn)變點。在堰頂水頭[9.39，9.47]之間按堰頂水頭級差1 cm計算溢洪道泄流能力時，發(fā)現(xiàn)在堰頂水頭9.43 m處數(shù)據(jù)發(fā)生異常變化，這在原來按堰頂水頭水位級差50 cm時是難以發(fā)現(xiàn)的。但公式都是依據(jù)經(jīng)典的水力學(xué)理論，本身并無問題，只能考慮公式的參數(shù)選擇是否合適。

表4 水位級差1 cm溢洪道部分泄流能力計算表

圖1 修正前溢洪道水位—流量曲線銜接放大圖

通過仔細(xì)分析研究并反復(fù)計算驗證，發(fā)現(xiàn)原孔堰流變流系數(shù)精度低是導(dǎo)致發(fā)生錯誤的原因，在將原來公式②、③中：當(dāng)開度e≥0.7H改變?yōu)椋篹≥0.721 3H，即相對開度≥0.721 3時，過渡為堰流，此時，將堰流變孔流公式的分界點堰頂水頭前移到9.15 m時，得到的計算結(jié)果銜接正常，見表5。

表5 修正后的溢洪道部分泄流能力計算表

采取以上的修正計算公式參數(shù)的方法，得到了銜接合理的計算結(jié)果，但這樣的方法必須在理論上找到合理的根據(jù)。根據(jù)《水力學(xué)計算手冊》[2]，對于寬頂堰底坎，判別孔堰流相對開度值為e≥0.65H，實用堰底坎，判別孔堰流相對開度值為e≥0.75H。通過對漳澤水庫溢洪道的駝峰堰分析，該駝峰堰既不是寬頂堰，又不同于實用堰，而是介于兩者之間，即界于[0.65，0.75]之間，取判別孔堰流相對開度值e=0.721 3H是符合實際情況的。而原來判別孔堰流相對開度e=0.7H可以認(rèn)為是低精度參數(shù)，高精度的水位計算必須有高精度的參數(shù)與之配套，這樣，可依此對原來的溢洪道泄流能力數(shù)據(jù)做一個精確的修正。

修改參數(shù)后的計算結(jié)果的圖形如圖2：

圖2 修正后溢洪道堰頂水頭—泄流量曲線銜接放大圖

4 結(jié)論

綜上所述，由于計算工具的局限，20世紀(jì)中期設(shè)計建造的許多水庫，其泄水建筑物流量計算公式參數(shù)普遍存在精度偏低的情況。在水資源日益緊缺的今天，這種情況已經(jīng)不能適應(yīng)水資源精細(xì)管理的現(xiàn)實。漳澤水庫通過引入新的參數(shù)表，使水庫泄水洞流量計算的精度有所提高，通過提高溢洪道相對開度精度，使水庫溢洪道流量公式參數(shù)得以修正的做法，為解決此類問題提供了有益的借鑒。

［1］張新建.山西省漳澤水庫防洪調(diào)度模擬系統(tǒng)［M］//水利部信息化工作領(lǐng)導(dǎo)小組辦公室，水利部水利信息中心編著.水利信息化技術(shù)和建設(shè)成果研究.北京：中國水利水電出版社，2005:413-420.

［2］武漢大學(xué)水利水電學(xué)院水力學(xué)流體力學(xué)教研室.李 煒主編.水力學(xué)計算手冊［M］.北京：中國水利水電出版社，2006:101、75.