張慶華,刁艷芳,程傳民,劉 巍

(1.山東農(nóng)業(yè)大學(xué)水利土木工程學(xué)院,山東泰安 271018;2.泰安市水利勘察設(shè)計研究院,山東泰安 271000)

梯級橡膠壩塌壩泄流計算仿真系統(tǒng)

張慶華1,刁艷芳1,程傳民2,劉 巍2

(1.山東農(nóng)業(yè)大學(xué)水利土木工程學(xué)院,山東泰安 271018;2.泰安市水利勘察設(shè)計研究院,山東泰安 271000)

針對梯級橡膠壩塌壩泄流計算過程比較復(fù)雜、煩瑣的問題,利用VFP語言編制了梯級橡膠壩塌壩泄流計算仿真系統(tǒng)。計算結(jié)果表明:在輸入橡膠壩的壩高、凈寬、側(cè)收縮系數(shù)、塌壩時間等參數(shù)的基礎(chǔ)上,該仿真系統(tǒng)能精確、快速地推算梯級橡膠壩不同塌壩方案的塌壩泄流過程和最大泄流量,從而能夠盡早預(yù)測塌壩產(chǎn)生的“人造洪峰”,對于橡膠壩工程的設(shè)計和梯級橡膠壩汛期控制運(yùn)用方案的編制具有一定的參考意義。

梯級橡膠壩;塌壩泄流過程;塌壩最大泄流量;塌壩時間;堰上水頭;仿真系統(tǒng)

1957年世界上首座橡膠壩在美國誕生,此后橡膠壩在我國和世界上許多國家得到了廣泛應(yīng)用。隨著同一河道上橡膠壩建設(shè)數(shù)量的增多,出現(xiàn)了梯級橡膠壩,從而汛期控制運(yùn)用問題越來越突出。實踐證明,在橡膠壩大量蓄水情況下,不合理的橡膠壩塌壩方案易造成河道“人造洪峰”現(xiàn)象,危及河道安全。因此,研究制訂科學(xué)合理的橡膠壩汛期控制運(yùn)行方案,對工程及下游河道安全度汛至關(guān)重要。

目前橡膠壩的運(yùn)行調(diào)度是研究的難點(diǎn)問題之一,得到越來越多的專家、學(xué)者的重視。例如Kahl等[1-3]對橡膠壩運(yùn)行參數(shù)的自動觀測與自動化控制進(jìn)行了研究;田忠祿[4]通過對具有決策支持系統(tǒng)的樞紐橡膠壩安全運(yùn)行管理自動化系統(tǒng)的研究,提出了利用現(xiàn)代計算機(jī)技術(shù)建立橡膠壩安全運(yùn)行管理自動化系統(tǒng);Zhang等[5-6]提出了橡膠壩建設(shè)、運(yùn)行和維護(hù)各個階段存在的問題及解決方法;徐智廷等[7]探討了梯級橡膠壩汛期調(diào)度運(yùn)用,提出了橡膠壩汛期調(diào)度運(yùn)用應(yīng)遵循的原則;李風(fēng)增[8]通過對膠南市風(fēng)河梯級橡膠壩綜合管理與運(yùn)用的研究,提出了橡膠壩運(yùn)行控制辦法;吳軍華等[9]通過對安陽市洹河梯級橡膠壩聯(lián)合運(yùn)行探討,提出了洹河梯級橡膠壩聯(lián)合運(yùn)行的要求和應(yīng)注意的問題;張世功等[10]針對2009年7月沂河兩次暴雨洪水過程,分析了沂河梯級橡膠壩調(diào)度運(yùn)用對河道洪水的影響;溫會軍[11]對遼河干流11座橡膠壩工程的聯(lián)合調(diào)度方案編制進(jìn)行了探討,分析了橡膠壩聯(lián)合調(diào)度的主要問題及技術(shù)難點(diǎn)。

在橡膠壩汛期控制運(yùn)用方案編制中,塌壩方案是主要內(nèi)容之一,即確定橡膠壩何時塌壩以及塌壩時間。而編制橡膠壩塌壩方案必須知道橡膠壩塌壩形成的泄流過程和塌壩產(chǎn)生的最大泄流量,這些問題屬于橡膠壩水力計算領(lǐng)域。在20世紀(jì)90年代中期以前,許多專家、學(xué)者開展了橡膠壩溢流等水力學(xué)問題的試驗研究,提出了橡膠壩擋水、完全塌壩后泄流量的計算方法,如1987年陳瑞[12]通過模型試驗對橡膠壩溢流的流量系數(shù)進(jìn)行了研究;1996年苑希民等[13]對溢流橡膠壩的泄流特性進(jìn)行了模型試驗; 2011年Zhang等[14]對橡膠壩塌壩泄流進(jìn)行了研究,提出了橡膠壩連續(xù)塌壩泄流過程計算方法。

梯級橡膠壩塌壩泄流計算一般采用逐時段演進(jìn)方法[14],是一個復(fù)雜而又煩瑣的過程,計算工作量大、耗時多,若用人工計算困難大。本文利用VFP語言,編制梯級橡膠壩塌壩泄流計算仿真系統(tǒng),可計算單級或梯級橡膠壩塌壩泄流過程線,同時計算塌壩產(chǎn)生的最大泄流量等數(shù)據(jù),為橡膠壩工程設(shè)計、編制橡膠壩汛期控制運(yùn)用方案提供依據(jù)。

1 梯級橡膠壩塌壩泄流計算原理

1.1 單級橡膠壩

橡膠壩塌壩泄流可按堰流公式計算。由于橡膠壩塌壩是一個連續(xù)的過程,堰流公式中的流量系數(shù)及堰上水頭隨塌壩時間而發(fā)生變化,橡膠壩塌壩后某時刻的泄流量為

式中:Q為某時刻橡膠壩的泄流量,m3/s;ε為側(cè)收縮系數(shù);σ為淹沒系數(shù);m為流量系數(shù);B為橡膠壩過流斷面凈寬,m;g為重力加速度,m/s2;H0為某時刻堰上水頭,m。

根據(jù)文獻(xiàn)[4],本文計算H0的公式為

式中:H水為某時刻橡膠壩上游水深,即水庫水深,m; H壩為某時刻橡膠壩的高度,m;v為某時刻橡膠壩上游水流行近流速,m/s;α為動能修正系數(shù)。

用式(1)(2)計算橡膠壩塌壩泄流量一般采用逐時段演進(jìn)方法,具體計算過程為:首先選定計算時段長Δt、壩袋塌落時間等基本參數(shù),然后從第一時段開始,逐時段計算壩袋塌落后壩的高度、壩上游蓄水庫容、水庫水深、堰上水頭、橡膠壩泄流量。

1.2 梯級橡膠壩

圖1為由n個橡膠壩組成的梯級橡膠壩(圖中Lk、Hk、Vk分別為第k級橡膠壩的蓄水長度、壩的高度、蓄水庫容;Lk-1,k為上級橡膠壩至下級橡膠壩設(shè)計水面之間的距離)。當(dāng)上級橡膠壩泄流時,對下級橡膠壩泄流產(chǎn)生影響,使下級橡膠壩泄流量加大,由于兩級橡膠壩之間有一定距離(Lk-1,k),上級橡膠壩的泄流到達(dá)下級橡膠壩需要一定的時間,因此,上級橡膠壩的泄流量不能簡單地疊加到下級橡膠壩泄流量,其影響可通過上級橡膠壩泄流后增加下級橡膠壩庫容(蓄水量)的方法計算。

a.上下級橡膠壩間隔距離較短。當(dāng)上下級橡膠壩間隔距離較短時(如Lk-1,k=0),上級橡膠壩的泄流很快到達(dá)下級橡膠壩,此時可不考慮上級橡膠壩泄流到達(dá)下游的延續(xù)時間,第t時段壩上游蓄水庫容為

式中:Vt、Vt-1分別為第t、第t-1時段初壩上游蓄水庫容,m3;Qt-1、Qt-1,上分別為第t-1時段本級、上級橡膠壩泄流量,m3/s。

b.上下級橡膠壩間隔距離較長。當(dāng)上下級橡膠壩間隔距離較長時,上級橡膠壩泄流到達(dá)下級橡膠壩有一定的延續(xù)時間T,由于橡膠壩塌壩泄流量是分時段進(jìn)行計算的,因此在仿真計算時,延續(xù)時間T是通過計算延續(xù)的時段數(shù)K來反映上級橡膠壩泄流對下級橡膠壩的影響,第t時段壩上游蓄水庫容可按下式計算:

式中:Qt-K,上為第t-K時段上級橡膠壩泄水流量, m3/s;K為上級橡膠壩泄水達(dá)到下游延續(xù)的時段數(shù),(取整數(shù));T為上級橡膠壩泄水到達(dá)下游橡膠壩延續(xù)的時間,s。

2 梯級橡膠壩塌壩泄流計算仿真系統(tǒng)

圖1 梯級橡膠壩示意圖

橡膠壩塌壩泄流過程計算是一個復(fù)雜而又煩瑣的過程,因此利用前述梯級橡膠壩塌壩泄流計算基本原理,編制了基于VFP語言的橡膠壩塌壩泄流計算機(jī)仿真系統(tǒng),該系統(tǒng)可根據(jù)橡膠壩的壩長、壩高、河道坡度等參數(shù),以1s時段長模擬計算橡膠壩不同塌落時間形成的流量過程,從而找到最大泄流量。該系統(tǒng)可計算最大5級橡膠壩塌壩泄流過程。

2.1 數(shù)據(jù)輸入

系統(tǒng)主要輸入3類數(shù)據(jù):①各級橡膠壩基本數(shù)據(jù),包括橡膠壩名稱、級數(shù)、凈寬、側(cè)收縮系數(shù)、正常擋水高度、正常擋水蓄水量、上游河道寬度、上游河道比降等。②橡膠壩塌壩方案數(shù)據(jù),包括上級橡膠壩泄流到達(dá)本級橡膠壩的時間、橡膠壩開始塌壩時間、橡膠壩完全塌壩時間等。③對不規(guī)則河道,還應(yīng)輸入橡膠壩擋水高度與蓄水量關(guān)系表。

2.2 塌壩泄流計算

系統(tǒng)在已知塌壩方案的基礎(chǔ)上,能進(jìn)行各級橡膠壩塌壩過程泄流量計算和塌壩形成的最大泄流量統(tǒng)計。

a.塌壩泄流過程計算。基本數(shù)據(jù)輸入完成后,運(yùn)行“梯級橡膠壩泄流計算”子程序,系統(tǒng)自動完成各級橡膠壩塌壩泄流過程計算。當(dāng)某級橡膠壩無基本數(shù)據(jù)時,系統(tǒng)僅計算其上級橡膠壩泄流過程。

b.塌壩最大泄流量統(tǒng)計。梯級橡膠壩泄流過程計算完成后,運(yùn)行“橡膠壩泄流匯總”子程序,系統(tǒng)自動完成各級橡膠壩塌壩產(chǎn)生的最大泄流量、發(fā)生時間、相應(yīng)壩高等計算結(jié)果統(tǒng)計。

2.3 計算結(jié)果查詢

系統(tǒng)可查詢橡膠壩的基本參數(shù),如各級橡膠壩塌壩泄流過程,塌壩產(chǎn)生的最大泄流量及其發(fā)生時間、相應(yīng)壩高、堰上水頭等。

2.4 計算結(jié)果輸出

各級橡膠壩塌壩泄流過程計算結(jié)果以報表的形式輸出,包括橡膠壩運(yùn)行時間、壩高、壩前蓄水庫容、堰上水頭、泄流量等。

各級橡膠壩塌壩最大泄流量統(tǒng)計結(jié)果輸出內(nèi)容包括橡膠壩名稱、蓄水量、塌壩時間、上游泄流到達(dá)時間、最大泄流量及其發(fā)生時間、相應(yīng)壩高、堰上水頭等。

以上各種計算結(jié)果還可以Excel表的形式輸出,便于統(tǒng)計分析。

3 算例計算

某河道有三級梯級橡膠壩,河道為矩形斷面,寬100 m,底坡為1/1 000。一級橡膠壩(上游)壩高4.0 m,庫容80萬m3;二級橡膠壩壩高3 m,庫容45萬m3;三級橡膠壩(下游)壩高4.0m,庫容80萬m3。一級橡膠壩距離二級橡膠壩的設(shè)計水面距離為50 m,三級橡膠壩設(shè)計水面與二級橡膠壩壩底相接。

3.1 塌壩方案

梯級橡膠壩塌壩方案見表1,其中完全塌壩時間為橡膠壩自設(shè)計擋水高度全部塌落需要的時間,塌壩開始時間為自最下游三級橡膠壩塌壩開始起算的一、二級橡膠壩開始塌壩時間(或為滯后時間)。

表1 橡膠壩塌壩方案

3.2 計算結(jié)果

啟動橡膠壩塌壩泄流計算仿真系統(tǒng),輸入橡膠壩基本數(shù)據(jù)、塌壩方案等,然后運(yùn)行橡膠壩塌壩泄流計算子程序(本例計算時段為1 s),得到梯級橡膠壩塌壩泄流過程及最大泄流量計算結(jié)果,如圖2和表2所示。

圖2 梯級橡膠壩塌壩泄流過程線

表2 塌壩最大泄流量計算結(jié)果

從圖2可以清楚地看到梯級橡膠壩各塌壩方案的泄流過程及變化趨勢,從表2中還可看到各塌壩方案的塌壩最大泄流量及發(fā)生時間,可編制橡膠壩控制運(yùn)用方案提供依據(jù)。

4 結(jié) 語

本文梯級橡膠壩塌壩泄流計算仿真系統(tǒng)具有通用性強(qiáng)、人機(jī)交互功能好、擴(kuò)展性強(qiáng)等特點(diǎn),能精確、快速地推算梯級橡膠壩不同塌壩方案的塌壩泄流過程和最大泄流量,從而能夠盡早預(yù)測塌壩產(chǎn)生的“人造洪峰”,對于橡膠壩工程的設(shè)計和梯級橡膠壩汛期控制運(yùn)用方案的編制具有一定的實用參考意義。

參考文獻(xiàn):

[1]KAHL T,RUELL S.Flashboard alternatives including rubber dams[C]//Proceedings of Waterpower,Niagara Falls.New York:[s.n.],1989:447-456.

[2]PLAUT R H,LIAPIS S I,TELIONIS D P.When the levee inflates[J].Civil Engineering,1998,68(1):62-64.

[3]CONRAD D A,MULDOON D P.Design and application of a PC/PLC utility rubberdam controlsystem[J]. Instrument Society of America.1995,6:351-362.

[4]田忠祿.具有決策支持系統(tǒng)的樞紐橡膠壩安全運(yùn)行管理自動化系統(tǒng)研究[J].中國水利水電科學(xué)研究院學(xué)報,2008,6(1):56-59.(TIAN Zhonglu.Study on the automation system for safe operation and management of rubber dams with decision supporting system[J].Journal of China Institute of Water Resources and Hydropower Research,2008,6(1):56-59.(in Chinese))

[5]ZHANG X Q,TAM P W M,ZHENG W.Construction, operation,and maintenance of rubber dams[J].Canadian Journal of Civil Engineering,2002,29(3):409-420.

[6]TAM P W M,ZHANG X Q.Management of rubber dams in Hong Kong[J].Canadian Journal of Civil Engineering, 1999,26(2):123-134.

[7]徐智廷,孫廷璽,張世功,等.沂河梯級橡膠壩汛期調(diào)度運(yùn)用原則探討[J].中國水利,2008(19):61-62.(XU Zhiting,SUN Tingxi,ZHANG Shigong,et al.Regulation of cascade rubber dams in Yihe River Basin[J].China Water Resources,2008(19):61-62.(in Chinese))

[8]李風(fēng)增.膠南市風(fēng)河梯級橡膠壩綜合管理與運(yùn)用[J].中國水利,2007(16):17.(LI Fengzeng.Integrated management and application of cascade rubber dams on Feng River of Jiaonan[J].China Water Resources,2007 (16):17.(in Chinese))

[9]吳軍華,孟紅軍.安陽市洹河梯級橡膠壩聯(lián)合運(yùn)行探討[J].河南水利與南水北調(diào),2010(4):31-32.(WU Junhua,MENG Hongjun.Integrated operation of step rubber dams on Huanhe River in Anyang City[J].Henan Water Resources&South to North Water Diversion,2010 (4):31-32.(in Chinese))

[10]張世功,孫廷璽,段偉華,等.沂河梯級橡膠壩調(diào)度運(yùn)用實況對洪水影響分析[J].治淮,2010(5):15-17. (ZHANG Shigong,SUN Tingxi,DUAN Weihua,et al. Analysis of Yihe cascade rubber dams operation and control impact on flood[J].Huaihe Management,2010 (5):15-17.(in Chinese))

[11]溫會軍.遼河干流橡膠壩聯(lián)合調(diào)度方案編制方法[J].現(xiàn)代農(nóng)業(yè)科技,2012(9):274-276.(WEN Huijun. Rubber dam joint operation scheme established method of Liaohe River[J].Modern Agricultural Sciences and Technology,2012(9):274-276.(in Chinese))

[12]陳瑞.對橡膠壩流量系數(shù)的探索[J].水利學(xué)報,1987 (6):36-39.(CHEN Rui.Research on flow coefficient of rubber dam[J].Journal of Hydraulic Engineering,1987 (6):36-39.(in Chinese))

[13]苑希民,林繼鐮.溢流橡膠壩泄流特性試驗[J].天津大學(xué)學(xué)報,1996,29(3):342-348.(YUAN Ximin,LIN Jilian.Rubber dam fall discharge characteristic test[J]. Journal of Tianjin University,1996,29(3):342-348.(in Chines))

[14]ZHANG Qinghua,DIAO Yanfang.Cascade rubber dams fall discharge calculation and analysis[J].Canadian Journal of Civil Engineering,2011,38(8):957-962.

Cascade rubber dams fall discharge calculation simulation system//

ZHANG Qinghua1,DIAO Yanfang1,CHENG Chuanmin2,LIU Wei2
(1.Water Conservancy and Civil Engineering College,Shandong Agricultural University,Tai'an 271018,China;2.Water Survey and Design Institute in Tai'an,Tai'an 271000,China)

In order to solve complex problem of cascade rubber dams fall discharge calculation,the cascade rubber dams fall discharge calculation simulation system was customized adopting VFP programming language.The results show that when parameters such as dam height,dam width,lateral contraction coefficient and complete dam fall time are input into this system,the discharge processes of cascade rubber dams under different fall schemes are calculated rapidly and the maximum discharges can be obtained as soon as possible.This system can provide basis for rubber dams design and planning the flood control and use of cascade rubber dams schemes.

cascade rubber dams;fall discharge process;maximum discharge;fall time;weir head;simulation system

TV135

:A

:1006-7647(2014)04-0053-04

10.3880/j.issn.1006-7647.2014.04.011

2013-0701 編輯:周紅梅)

山東省重大水利科研與技術(shù)推廣項目(2012);水利部公益性行業(yè)專項(201201115)

張慶華(1960—),男,山東招遠(yuǎn)人,教授,碩士,主要從事水利工程建設(shè)與管理研究。E-mail:zqh@sdau.edu.cn