孫 卓,魏艷鳳
(黑河水文局,黑龍江黑河164300)
河槽匯流曲線預(yù)報(bào)方法應(yīng)用分析
孫 卓,魏艷鳳
(黑河水文局,黑龍江黑河164300)
河槽匯流曲線法是河段上、下游斷面之間流量預(yù)報(bào)的常用方法。文章以馬斯京根河道匯流曲線和瞬時(shí)河槽匯流曲線基本原理為基礎(chǔ),對(duì)兩種方法進(jìn)行了應(yīng)用分析和探討。兩種途徑得到的河槽匯流單位線,經(jīng)多次洪水驗(yàn)證分析,兩種匯流曲線的預(yù)報(bào)精度無(wú)明顯差異。針對(duì)實(shí)例流域具體情況,在參數(shù)計(jì)算方法選取,預(yù)報(bào)斷面區(qū)間來(lái)水量確定等方面提出了具體建議。
河槽匯流曲線;馬斯京根匯流曲線;瞬時(shí)匯流曲線;對(duì)比分析;預(yù)報(bào)精度
目前使用較廣的有馬斯京根匯流曲線與瞬時(shí)河槽匯流曲線。二者都是從聯(lián)解河段水量平衡方程及槽蓄方程得到的。不同點(diǎn)是二者的槽蓄方程不一。
馬斯京根法是用w=k[xQ上+(1-x)Q下],后者的W=KQ下,實(shí)際上是馬斯京根法在x=0條件下的一個(gè)特例。瞬時(shí)單位線在使用時(shí)要轉(zhuǎn)化為指定時(shí)段的單位線。
瞬時(shí)單位線有參數(shù)n、k。n與k的乘積稱為單位線滯時(shí),反映了一次洪水過(guò)程自上游站至下游站的平均傳播時(shí)間。n愈大,l(特征河長(zhǎng))、k愈小,曲線愈尖瘦,說(shuō)明河槽調(diào)節(jié)能力小,k值愈大,即傳播時(shí)間愈長(zhǎng),河槽調(diào)節(jié)能力也愈大,洪水波就顯著展開。馬斯京根河槽單位線共有x、n(河段分段數(shù))、k 3個(gè)參數(shù),要求計(jì)算時(shí)段△t≈k,x愈小,洪水過(guò)程展開程度則大,河槽調(diào)節(jié)能力較強(qiáng),下面分述兩種匯流曲線的實(shí)例應(yīng)用分析[2]。
洪水波在傳進(jìn)中變形或河槽的調(diào)節(jié)作用,就其內(nèi)因來(lái)說(shuō)是不穩(wěn)定流的水深和比降不斷變化引起的:①水流受河槽斷面形狀、糙率的約制和水流內(nèi)摩阻力造成斷面流速分布不均勻;②從流域各處匯入徑流或閘壩泄水變化使河槽水深、流速等水力要素隨時(shí)間而變,于是流量也隨時(shí)間而變。組成下斷面t時(shí)刻流量Qt的i個(gè)流束的水質(zhì)點(diǎn)是分處河段各處的[1]。它們各自的流速都不相等,但都能在t時(shí)刻傳至下斷面組成了Q下t,即:
2.1 馬斯京根匯流曲線應(yīng)用分析
匯流曲線推求:
馬斯京根匯流曲線推求首先要求對(duì)預(yù)報(bào)河段進(jìn)行分段處理,各河段均采用同一短河段的參數(shù)K、X和△t來(lái)求出河段流量演算公式系數(shù)C0、C1、C2,然后用單位入流求出匯流曲線U(t),計(jì)算實(shí)例見表1及圖1。
表1中,取計(jì)算時(shí)段△t=6h,實(shí)例站湯旺河五營(yíng)站至伊新站河段長(zhǎng)為98 km,洪水傳播時(shí)間約20 h。
表1 五營(yíng)~伊新站河槽匯流曲線計(jì)算表
表2 五營(yíng)~伊新站出流量演算表(馬斯京根匯流曲線)
2.2 瞬時(shí)河槽匯流曲線應(yīng)用分析
匯流曲線推求:
由于選取的時(shí)段長(zhǎng)△t不同,同一次洪水分析出來(lái)的河槽單位線也不相同。為了使分析的單位線不受選取時(shí)段△t長(zhǎng)短的影響,假設(shè)單位入流發(fā)生在t趨于零的極小時(shí)段內(nèi),即單位入流是瞬時(shí)的,相應(yīng)于這個(gè)入流的河段出流過(guò)程就是瞬時(shí)河槽單位線[3]。應(yīng)用時(shí),根據(jù)瞬時(shí)河槽單位線轉(zhuǎn)換成需要的時(shí)段△t單位線。
常用的瞬時(shí)河槽單位線為納須瞬時(shí)流域單位線,其表達(dá)式為:
應(yīng)用上式的關(guān)鍵是求參數(shù)n、k,一般采用求流量過(guò)程的矩得到(公式略),計(jì)算結(jié)果n=2.38,k=5.43。匯流曲線計(jì)算時(shí)根據(jù)n、k值查納須瞬時(shí)單位線S(t)曲線表即可,詳見表3~表6和圖2。
圖1 湯旺河五營(yíng)~伊新站河槽單位線(馬斯京根法)
表3 湯旺河五營(yíng)站流量過(guò)程求矩計(jì)算表(△t=6h)
表4 湯旺河伊新站流量過(guò)程求矩計(jì)算表(△t=6h)
表5 湯旺河五營(yíng)~伊新站河槽瞬時(shí)單位線計(jì)算表(△t=6h)
圖2 湯旺河五營(yíng)~伊新站瞬時(shí)河槽單位線
表6 五營(yíng)~伊新站出流量演算表(瞬時(shí)單位線)
1)以上兩種途徑得到的河槽匯流單位線在應(yīng)用于同一場(chǎng)洪水?dāng)M合預(yù)報(bào)時(shí),洪峰流量誤差分別為15.3%和9.5%,均滿足預(yù)報(bào)精度要求,經(jīng)多次洪水驗(yàn)證分析,兩種匯流曲線的預(yù)報(bào)精度無(wú)明顯差異。
2)馬斯京根河槽匯流曲線推求時(shí)雖比較繁瑣,所需參數(shù)較多,計(jì)算量較大,但比較容易得出合適的單位線;瞬時(shí)河槽匯流曲線推求時(shí)計(jì)算量雖較少,且可以查S(t)曲線表得到匯流曲線,但參數(shù)n、k值由于洪水場(chǎng)次選取等因素,經(jīng)常出現(xiàn)不合理的現(xiàn)象。
3)天然河道預(yù)報(bào)斷面之間一般距離較長(zhǎng),區(qū)間入流往往成為河道流量演算預(yù)報(bào)的重要影響因素,選擇實(shí)例洪水進(jìn)行匯流曲線參數(shù)分析時(shí)如何確定區(qū)間入流非常重要。建議根據(jù)支流測(cè)站來(lái)水按計(jì)算時(shí)段和洪水傳播時(shí)間與上游站流量進(jìn)行疊加來(lái)確定具體上游來(lái)水量;對(duì)于無(wú)測(cè)站控制的區(qū)間入流應(yīng)根據(jù)區(qū)間降雨情況進(jìn)行經(jīng)驗(yàn)預(yù)估來(lái)確定上游來(lái)水。由于實(shí)例流域資料條件限制,上述分析結(jié)論還很膚淺,僅供交流參考。
[1]長(zhǎng)江水利委員會(huì).水文預(yù)報(bào)方法:第二版[M].北京:水利電力出版社,1993:29-70.
[2]葛守西.現(xiàn)代洪水預(yù)報(bào)技術(shù)[M].北京:中國(guó)水利水電出版社,1999:125-145.
[3]水利部.水文情報(bào)預(yù)報(bào)規(guī)范[M].北京:中國(guó)水利水電出版社,2000:9-11.
P333
B
1007-7596(2015)05-0027-04
2014-10-31
孫卓(1979-),女,吉林長(zhǎng)嶺人,工程師;魏艷鳳(1976-),女,黑龍江嫩江人,工程師。