孫　磊，李虎成

(廣東省水利水電科學(xué)研究院，廣州 510635)

?

佛岡縣四九水整治工程水文水利分析計算

孫磊，李虎成

(廣東省水利水電科學(xué)研究院，廣州 510635)

佛岡縣四九水流域整治工程上游起于石甕水庫壩下，下游止于入潖江河口。該流域面積116km2，流域內(nèi)無實測洪水資料，本次采用設(shè)計暴雨推求設(shè)計洪水，并依據(jù)其成果推求到河道整治前后的設(shè)計洪水水面線。

整治工程；設(shè)計暴雨；設(shè)計洪水；水面線

0　前　言

四九水流域(又名蘊洞水)發(fā)源于四九鎮(zhèn)中心、黃竹田村上的紅坪腦(山頂名)。在發(fā)源地的群山中先向東流，經(jīng)石甕水庫，跌落橫坑村又折向東北方向于亞鬢山南麓腳下蜿蜒經(jīng)回心村、官山材流入四九圩(四九鎮(zhèn)所在地)，再折向西南方向經(jīng)凌塘村過白石坳梳入湯塘鎮(zhèn)的上黎、良安、圍鎮(zhèn)、脈塘等村，于湯塘鎮(zhèn)的北側(cè)匯入潖江河。

四九水流域面積116km2(原黃花水是流入四九水后匯入潖江，1956 年后，因建了黃花河水庫，截流54km2，黃花水的主流于汛期經(jīng)水庫溢洪道溢出直入潖江河，因此四九水的流域面積于1956年后減去黃花水54km2，則實為116km2，河長25km，河床評價比降為3.8‰。

1　洪　水

1.1暴雨洪水特性

四九水流域是以降雨為主要補(bǔ)給的河流，中下游河道徑流及其變化主要由上游及其區(qū)間降雨情況所決定，降雨量在年內(nèi)變化決定徑流年內(nèi)變化。

四九水流域每年4月—9月為汛期，11月—翌年的3月為枯水期。全流域的大洪水都是由鋒面雨產(chǎn)生。由于暴雨量大，流域坡陡，水系又是葉脈分布，洪水容易集中，洪峰上漲速度快，具有山區(qū)洪水特點，峰型尖瘦，漲落較快，洪峰一般呈連續(xù)的多峰型。發(fā)洪時間為4月—9月，洪峰出現(xiàn)最多的時段為5月—7月，其中出現(xiàn)次數(shù)以6月最為集中。

1.2設(shè)計洪水計算

四九水流域無實測洪水資料，本次采用設(shè)計暴雨推求設(shè)計洪水。

1.2.1設(shè)計暴雨特征參數(shù)

根據(jù)廣東省水文局2003年《廣東省暴雨參數(shù)等值線圖》查得四九水各歷時點暴雨特征參數(shù)和點面換算系數(shù)，并根據(jù)Cv查取相應(yīng)于設(shè)計頻率P的kp值，流域點暴雨特征參數(shù)和點面換算系數(shù)及不同頻率的kp值見表1。

1.2.2未考慮水庫調(diào)蓄的設(shè)計洪水

根據(jù)《廣東省暴雨徑流查算圖表》，產(chǎn)流分區(qū)取“內(nèi)陸”，設(shè)計雨型分區(qū)取“北江中下游”，α/t/F關(guān)系分區(qū)取為Ⅵ，設(shè)計暴雨低區(qū)的αt～t～F關(guān)系取“暴雨低區(qū)”，m1-θ分區(qū)取大陸高區(qū)；ui-xi關(guān)系取“Ⅱ”；m-θ分區(qū)的選取，綜合考慮河道平均坡降、區(qū)域平均高程、區(qū)域土壤及植被等因素進(jìn)行選取，并與綜合單位線法計算結(jié)果進(jìn)行對比后根據(jù)區(qū)域下墊面情況做適當(dāng)調(diào)整[1]。

根據(jù)設(shè)計暴雨參數(shù)采用成果和各河段的流域地理特征參數(shù)，使用《廣東省暴雨徑流查算圖表(使用手冊)》推薦的廣東省綜合單位線法和推理公式法分別計算四九水各河段控制斷面設(shè)計洪水(未考慮四九水和沙路水因上游水庫的調(diào)洪作用)，兩種方法計算得到的成果相差小于10%，說明成果較為合理，本次采用綜合單位線計算的設(shè)計洪水成果，計算成果見表2。

表1　四九水流域暴雨參數(shù)成果

表2　設(shè)計洪水計算結(jié)果(未考慮水庫調(diào)洪)　m3/s

1.2.3水庫調(diào)蓄影響

四九水的K0+000位于石甕水庫壩下，沙路水(支流5)上游有止貝冚水庫，對四九水和沙路水的設(shè)計洪水產(chǎn)生較大影響，為考慮水庫的洪水調(diào)節(jié)作用，需分別對石甕水庫和止貝冚水庫進(jìn)行調(diào)洪演算[2]。

石甕水庫和止貝冚水庫泄水建筑物均為開敞式無閘門控制寬頂堰，堰頂高程分別為139.30m和62.75m，堰頂寬分別為20m和10m；水庫按照從正常蓄水位(即堰頂高程)起調(diào)，分別進(jìn)行調(diào)洪演算，計算結(jié)果見表3。

表3　最大下泄流量計算結(jié)果表　m3/s

1.2.4考慮水庫調(diào)蓄后的設(shè)計洪水

根據(jù)設(shè)計暴雨參數(shù)采用成果和各河段的流域地理特征參數(shù)，使用《廣東省暴雨徑流查算圖表(使用手冊)》推薦的廣東省綜合單位線法計算四九水各河段控制斷面設(shè)計洪水，其中四九水和沙路水因上游均有水庫，將區(qū)間設(shè)計洪水和水庫調(diào)洪后的下泄流量按峰峰遭遇進(jìn)行疊加，作為工程設(shè)計洪水成果，計算成果見表4。

表4　設(shè)計洪水計算結(jié)果(考慮水庫調(diào)洪)　m3/s

2　設(shè)計水面線

采用2015年6月新測斷面資料，四九水干流18.83km河長共設(shè)370個斷面，斷面間距為50m左右，在特殊斷面適當(dāng)調(diào)整。

2.1計算方法

采用mike11軟件建立四九水及各支流的一維河網(wǎng)數(shù)學(xué)模型。該模型的基本控制方程采用帶旁側(cè)入流(出流)的圣維南方程組進(jìn)行求解，即水流連續(xù)性方程

(1)

動量方程：

(2)

2.2河道糙率

根據(jù)現(xiàn)場調(diào)查，現(xiàn)狀河道河床床面不平整，底坡有凹凸?fàn)睢Ｉ嫌魏訛┑厣嫌谢牟莺偷桶墓嗄?，河床有鵝卵石、礫石間有孤石，因此按山區(qū)河流取糙率為0.040；中游下游河段有雜草和石塊，糙率取0.033。

整治后護(hù)腳采用較低的格賓石籠擋土墻，河床進(jìn)行清淤，邊坡采用人工草皮護(hù)坡，阻礙行洪的低矮灌木清除，上游將鵝卵石和塊石用來護(hù)腳。四九水整治后糙率中下游取0.030，上游取0.035。

2.3起推水位

一維河網(wǎng)模型上游開邊界采用流量控制，包括石翁水庫和止貝冚水庫的下泄流量；下邊界條件采用水位控制。

四九水流域整治工程的下游斷面為四九水與潖江交匯處的河口斷面，起推水位取河口斷面的潖江水位。根據(jù)《清遠(yuǎn)市洪情概貌》(梁榮燦編著)，結(jié)合潖江河佛岡縣城段洪水水面線成果及近期大洪水調(diào)查成果，采用極值倍比法推求河口斷面水位成果見表5。

表5　河口斷面起推水位計算結(jié)果　m

2.4水面線成果

根據(jù)以上計算參數(shù)和計算方法，計算四九水干流各頻率現(xiàn)狀河道水面線，現(xiàn)狀部分?jǐn)嗝嫠怀^了兩岸地面高程，部分橋涵已漫水，阻水嚴(yán)重；整治后的水面線根據(jù)現(xiàn)狀河道水面線成果及各斷面的兩岸現(xiàn)狀，以盡量不占地為原則，按照原堤防河寬布置河流堤線，進(jìn)行護(hù)岸和河道的清淤疏浚，整治前、后河道水面線成果見表6。

表6　四九水干流整治前后水面線計算結(jié)果　m

續(xù)表6　四九水干流整治前后水面線計算結(jié)果　m

從表6可以看出，四九水通過整治后，比整治前水面線有所降低，平均降低0.42m，最大降低0.72m，整治后使四九水行洪能力加大明顯。

[1]王利，曹秀菊，王玉坤 .勞動湖綜合治理工程水文分析[J].黑龍江水利科技，2002，30(01)：38-39.

[2]韋翠娥 .蓮塘河流域治理工程水文水利分析計算[J].吉林水利，2015(02)：41-44.

1007-7596(2016)05-0049-03

2016-03-15

孫磊(1981-)，男，河南蘭考人，工程師，從事水利水電工程規(guī)劃與設(shè)計工作；李虎成(1982-)，男，河南開封人，高級工程師，從事水力學(xué)及河流動力學(xué)方面研究工作。

TV882

B