水箱模型應(yīng)用研究

水箱模型將復(fù)雜的降水徑流過程簡單慨化為流域的蓄水與出流關(guān)系進行模擬.針對皂河小流域的降水、徑流和水循環(huán)等水文特性,建立了4層串聯(lián)水箱模型,模擬流域的降水徑流過程,計算流域徑流量,并結(jié)合河流水質(zhì)模型,綜合模擬河流流量及主要水質(zhì)指標(biāo)變化.結(jié)果表明:各項指標(biāo)的模擬值與實測值吻合較好,說明水箱模型用于小流域的徑流模擬較為準(zhǔn)確,可與水質(zhì)模型耦合使用,綜合模擬小流域的水文水質(zhì)變化。 ?

在水箱模型的擬線性狀態(tài)基礎(chǔ)上建立擴展的卡爾曼濾波方程，構(gòu)造了水箱模型的實時校正模型。實例表明，水箱模型的實時校正模型較不帶校正的水箱模型，預(yù)報精度有很大的提高。研究表明，水箱模型的實時校正模型對不良的母體模型參數(shù)有較強的適應(yīng)性或容錯性。 ?

應(yīng)用水箱模型建立三峽入庫洪水實時連續(xù)預(yù)報模型的方法。按照河段流量傳播時間，將入庫站(寸灘)以上干、支流劃分為若干子河段，各子河段按照計算時段長分成若干單元河段，各單元河段區(qū)間降雨徑流預(yù)報采用水箱模型，河道流量演算采用連續(xù)馬斯京根法。對河系洪水預(yù)報模型的檢驗表明，預(yù)報精度在85%以上，能夠滿足實時洪水預(yù)報的要求。 ?

根據(jù)現(xiàn)有資料，采用二層直列式水箱模型。在參數(shù)率定時，由于流域中僅有福才站以上流域的降雨徑流觀測記錄，而相應(yīng)時段松濤水庫的庫水位記錄缺少，無法還原庫周區(qū)（即福才站至大壩區(qū)間流域）洪水過程，因此無法模擬庫周區(qū)洪水的入庫過程，相應(yīng)的模型參數(shù)不能正確率定。由于庫周區(qū)與入庫站（福才站）控制流域的地形、地貌等極相似，并且流域面積較小，故在預(yù)報時庫周區(qū)洪水預(yù)報模型中的參數(shù)直接引用福才站的參數(shù)。這對預(yù)報精度會有一定的影響，但以后隨資料的積累這一問題將會逐漸得到解決。

水箱模型的基本思想是假定流域中的出流及下滲量，是流域相應(yīng)蓄水深的函數(shù)。從這點出發(fā)，將流域的雨洪轉(zhuǎn)化過程的各個環(huán)節(jié)（產(chǎn)流、坡面匯流、河道匯流等）用若干個彼此相聯(lián)的水箱進行模擬。以水箱中的蓄水深度為控制，計算流域的產(chǎn)流、匯流及下滲過程。較小的流域，可用若干個串聯(lián)的直列式水箱模型模擬出流和下滲過程?？紤]降雨和產(chǎn)、匯流的不均勻，需要分區(qū)計算的較大流域，可用若干個并聯(lián)組合的水箱模型結(jié)構(gòu)，模擬整個流域的雨洪轉(zhuǎn)化過程。直列式水箱模型，每層水箱的側(cè)邊有出流孔，底部有下滲孔。上層水箱的入流為流域面上的降雨，下層水箱的入流為上層水箱的下滲量。各層水箱的出流量可理解為流域各蓄水層形成的不同水源的徑流量。對二層水箱模型來說，上層水箱的出流量相當(dāng)于地表徑流，下層水箱的出流量相當(dāng)于地下徑流。

流域上的降水，輸入第一層水箱，形成蓄水深Z1，當(dāng)Z1大于出流孔高H11時，開始出流，同時它的下滲水量注入第 二層水箱，使第二層水箱的蓄水深Z2開始上升，當(dāng)Z2大于出流孔高H21時，第二層水箱開始出流。若降雨繼續(xù)，則各出流孔繼續(xù)出流，相應(yīng)時刻各出流量疊加，即得流域的總出流過程。

第一層水箱設(shè)有二個出流孔，孔高為H11、H12，下滲系數(shù)為R10。假如在T時 該層的蓄水深為Z(T)，則T時的出流量Q(T)和下滲量F(T) 的計算如下:

Q(T)=0 當(dāng)Z(T)

Q(T)=[Z(T)-H11]×R11 當(dāng)H11

Q(T)=[Z(T)-H11]×R11+[E(T)-H12]×R12 當(dāng)Z(T)>H12

F(T)=Z(T)×R10 ，或當(dāng)?shù)卓赘叱龅酌鍴10時,F(T)=[Z(T)-H10]×R10

T時段末的剩余蓄水深Z`(T)為:

Z`(T)=Z(T)-Q(T)-F(T)

T+1時段的蓄水深為:

Z(T+1)=Z(T)+P(T)-E(T)

P(T)為T時段降雨,E(T)為T時段蒸發(fā)等損失。一次降雨過程的出流和下滲過程可連續(xù)計算。

第二層以下等水箱的出流和下滲量計算,除了以上一層水箱的下滲量為入流外,其余計算同上。

該結(jié)構(gòu)上層水箱有二個出流孔,一個下滲孔，下層水箱一個出流孔,一個下滲孔?？紤]該地區(qū)處于南北氣候過渡地帶，氣候變化大，對上層水箱的下滲孔高出底面一個H10高度（即死藏水），這樣該結(jié)構(gòu)包含著9個待定系數(shù)。這9個系數(shù)確定后，便可由降雨預(yù)報出流過程。

模型參數(shù)雖有一定的物理意義，如上層水箱的H11有初損值的意義，但與初損值又不完全相同，出流系數(shù)R12、R11與退水曲線的指數(shù)方程的指數(shù)有一定的聯(lián)系，但又不能直接確定。因此，這些參數(shù)的確定，都需經(jīng)過分析和反復(fù)試算，以計算與實測洪水過程擬合最好為確定模型參數(shù)的準(zhǔn)則。 ?

該模型最早由日本菅原正巳博士在40年代提出 ，主要發(fā)展時期是 1 95 1～ 1 95 8年 ,6 0年代已應(yīng)用于日本各河流。 1 971～ 1 974年 ,主要是70年代初，水箱模型則普遍應(yīng)用于世界。

水箱模型又稱坦克（Tank）模型 ，也稱黑箱模型。水箱模型是通過降雨過程計算徑流過程的一種降雨徑流模型。水箱模型是應(yīng)用較為廣泛的流域水文模型。 ?