許 淑 敏

(天津大學(xué)，天津 300072)

0 引 言

水箱模型是日本的菅原正于上世紀(jì)50年代提出，并不斷發(fā)展成為一種被世界各國廣泛采用的流域水文模型。最早提出的水箱模型是黑箱模型，后經(jīng)眾多學(xué)者研究現(xiàn)已發(fā)展成為一種概念性模型。在我國南方地區(qū)的濕潤季節(jié)，往往采用串聯(lián)水箱模型。

在使用水箱模型進(jìn)行計(jì)算時(shí)，必須確定每個(gè)水箱的初始蓄水深度(底水)。它的大小，對次雨產(chǎn)生的徑流量和過程有密切的關(guān)系，即初始蓄水深高，反映前期濕潤，次雨的產(chǎn)流快，損失小，徑流量較大，反之則反。因此確定底水是水箱模型在應(yīng)用中一個(gè)重要問題。

當(dāng)起漲流量較大時(shí)(此時(shí)各層水箱底水均超過邊孔高度)可采用起漲流量確定各層水箱的初始蓄水深度(底水)；久旱之后的次雨模擬則可認(rèn)為第一層水箱的底水為零，甚至第二層也可認(rèn)為為零。但是對于起漲時(shí)處于退水階段(此時(shí)初始蓄水深度(底水)低于水箱的邊孔高度)，這種情況并未有好的解決辦法。實(shí)際工作中為了避開并解決這一問題，采用的辦法是將底水參數(shù)化[1]。

將底水參數(shù)化固然使得問題得到簡化，但是缺點(diǎn)也是顯而易見：將各次降雨初始的初始蓄水深度(底水)作為不變的參數(shù)處理使得該參數(shù)失去本來具有的物理意義，同時(shí)也使得優(yōu)選后的其他參數(shù)也偏離其物理意義，進(jìn)而使得作為概念性模型的水箱模型再次淪為黑箱模型。

1 利用次雨初損確定底水及模擬成果

如圖1所示：對于第一層水箱，底水=最大初損-次雨初損，因而對于起漲時(shí)處于退水階段這種情況可采用此方法確定第一層水箱的初始蓄水深度(底水)。

圖1 底水示意圖Fig.1 bottom water schematic

1.1 資料整理

為檢驗(yàn)此法效果，本文選用湖北清江恩施流域1996-1999年4年資料，據(jù)該流域產(chǎn)流特點(diǎn)，確定采用三層串聯(lián)水箱模型。為模擬地面徑流的非線性效果，第一層水箱開2個(gè)邊孔，第三層水箱不開底孔(流域閉合)，邊孔位于底部。同時(shí)為了模擬調(diào)蓄作用并聯(lián)一個(gè)線性水庫。

由上述4年資料中選取起漲時(shí)處于退水階段這種情況的洪水，可按如下原則選取：①為了保證起漲流量較低(處于退水段)，選擇在次雨之前10天左右沒有降雨的洪水。②首選峰形

較好的單峰洪水，若滿足條件的洪水較少，則可選擇有多個(gè)洪峰的洪水，但本文將有針對性的研究主要洪峰。

按上述2原則共選出了10場(960418，960503，960514，960603，960722，970606，970706，990522，990626，990716)洪水。

1.2 底水的確定

將各場洪水起漲點(diǎn)之前的降雨累加作為次雨初損，選其最大者作為最大初損；根據(jù)底水=最大初損-次雨初損求得第一層水箱底水。

因挑選的洪水均在長期無雨退水段產(chǎn)生的，上兩層水箱底水均低于出流孔高度，故可認(rèn)為起漲流量是由第三層水箱引起的；根據(jù)Qg=αx[2]求得第三層水箱的底水。

本次研究中將第二層水箱底水參數(shù)化，具體取值由參數(shù)優(yōu)選確定。各場次洪水次雨初損及水箱底水取值如表1。

表1 洪水初損、第一層水箱底水、起漲流量、第三層水箱底水值Tab.1 The value of elemantary losing of the flood,bottom water of the first layer of water bank,the rising discharge and the third layer of water bank

1.3 率定驗(yàn)證

在表1的基礎(chǔ)上進(jìn)行參數(shù)優(yōu)選，將十場流量資料分成兩組，第一組為：960418，960503，960514，960603，960722場洪水，用來率定模型；第二組為：970606，970706，990522，990626，990716場洪水，用來驗(yàn)證模型。經(jīng)過反復(fù)調(diào)整參數(shù)，最終得到的參數(shù)成果如表2。

1.4 效果檢驗(yàn)

據(jù)水文情報(bào)預(yù)報(bào)規(guī)范，洪峰及峰現(xiàn)時(shí)間精度規(guī)定，降雨徑流預(yù)報(bào)接受的許可誤差為：實(shí)測洪峰流量的20%；峰現(xiàn)時(shí)間接受的許可誤差為預(yù)報(bào)根據(jù)時(shí)間至實(shí)測洪峰出現(xiàn)時(shí)間之間時(shí)距的30%。本次模擬研究的評(píng)定結(jié)果如表3所示。

由表3可知，在十場洪水中，若以洪峰流量作為評(píng)價(jià)標(biāo)準(zhǔn)，則有8場合格預(yù)報(bào)，2場不合格預(yù)報(bào)；若以峰現(xiàn)時(shí)間為評(píng)價(jià)標(biāo)準(zhǔn)，也有8場合格預(yù)報(bào)，2場不合格預(yù)報(bào)。合格率均為80%，屬于預(yù)報(bào)項(xiàng)目精度等級(jí)中的乙級(jí)。對于不同工程，安全要求是不同的，無論是更關(guān)注洪峰流量還是峰現(xiàn)時(shí)間，本研究方法都可以滿足精度要求。但本研究還存在一個(gè)問題，若某工程對洪峰流量和峰現(xiàn)要求都有較高要求時(shí)，本法精度尚低，合格率為60%，預(yù)報(bào)項(xiàng)目精度為丙級(jí)，有待進(jìn)一步研究。

表2 三層串聯(lián)水箱模型優(yōu)選參數(shù)成果表Tab.2 The result table of optimizing parameters ofthree layer cascade tank model

表3 洪峰及峰現(xiàn)時(shí)間精度評(píng)定表Tab.3 Precision evaluation form of peak discharge and the time of peak appearing

最終的模擬效果如圖2。其中點(diǎn)狀圖為實(shí)測流量數(shù)據(jù)，藍(lán)色線狀圖為模擬數(shù)據(jù)?？v坐標(biāo)Q表示洪水不同時(shí)刻流量，橫坐標(biāo)t表示時(shí)間段數(shù)，本文是以3 h作為單位時(shí)間段間隔。綜合來看，模擬結(jié)果流量過程線與實(shí)測數(shù)據(jù)變化趨勢較一致，整體模擬效果比較好。

2 結(jié) 論

(1)將水箱模型的底水與次雨聯(lián)系起來，使參數(shù)更符合其物理意義。

(2)以洪峰流量或峰現(xiàn)時(shí)間作為評(píng)定標(biāo)準(zhǔn)，本研究方法的合格率較高，可達(dá)80%，預(yù)報(bào)項(xiàng)目精度等級(jí)為乙級(jí)。

(3)以峰現(xiàn)時(shí)間和峰現(xiàn)時(shí)間作為評(píng)定標(biāo)準(zhǔn)，本研究方法的合格率為60%，預(yù)報(bào)項(xiàng)目精度等級(jí)為丙級(jí)，有待進(jìn)一步提高。

圖2 恩施流域模擬效果對比圖Fig.2 Simulation results of enshi basin

3 結(jié) 語

文中第一層水箱底水確定采用次雨初損確定使得該參變量比起將其參數(shù)化處理更加符合其物理意義，然模擬壤中流的水箱底水仍然作為參數(shù)處理未免略顯不足，故而筆者認(rèn)為模擬壤中流的水箱底水的確定值得進(jìn)一步探討，以期得到符合其物理意義的確定方法。

□

[1] 徐宗學(xué).水文模型[M].北京：科學(xué)出版社，2008.

[2] 湯成友，項(xiàng)祖?zhèn)?水箱模型在大尺度流域?qū)崟r(shí)洪水預(yù)報(bào)模型研制中的應(yīng)用[J].水文，2007，27(5).