趙蕾蕾 張玉杰

（河南省水利勘測(cè)設(shè)計(jì)研究有限公司）

1 研究區(qū)域概況

1.1 自然地理

輝縣東河屬海河流域漳衛(wèi)河水系百泉河支流，衛(wèi)河（共產(chǎn)主義渠）二級(jí)支流，發(fā)源于輝縣太行山南麓方山，流經(jīng)窯村出山，南下經(jīng)雷莊、大官莊，至輝縣市東南入百泉河，于新鄉(xiāng)縣北注入衛(wèi)河，干流長(zhǎng)12.50 km。輝縣東河流至大官莊村東與總干渠相交，交叉斷面以上集水面積17.10 km2，干流長(zhǎng)7.50 km，干流平均坡降為21.50‰，流域形狀似帶狀，地勢(shì)北高南低，由北向西南傾斜。流域內(nèi)為山丘區(qū)，山區(qū)河段為峽谷，局部基巖裸露，交叉處河床由壤土組成。輝縣東河為季節(jié)性河流，平時(shí)河床干涸，河道內(nèi)種植農(nóng)作物，兩岸樹(shù)木較少。

1.2 區(qū)域地物

總干渠運(yùn)行以來(lái)，輝縣東河下游邊界條件發(fā)生了較大變化，主要有四個(gè)變化：①由于輝縣縣城快速發(fā)展，河道兩側(cè)已發(fā)展為城市，導(dǎo)致輝縣東河下游河道束窄，行洪斷面減小。②輝縣東河暗渠交叉斷面下游約500 m 處輝常路上阻水涵洞拆除后，修建輝常路橋，橋下為3 孔7 m×6 m 的橋孔。③位于總干渠上輝縣東河暗渠下游的楊莊溝渡槽改變下游原行洪通道，在楊莊溝排水渡槽出口海漫段下游的溝道右岸，開(kāi)槽設(shè)置進(jìn)水口，經(jīng)由總干渠右岸永久占地線以外往北，在輝常路處轉(zhuǎn)向西，沿輝常路往西排水，在輝常路橋處入輝縣東河。④輝縣東河暗渠交叉斷面下游700 m處，幸福橋兩孔橋孔阻水。

1.3 風(fēng)險(xiǎn)分析

①輝縣東河暗渠交叉斷面輝縣東河水位超過(guò)設(shè)計(jì)值較大，洪水漫溢進(jìn)入總干渠，威脅到總干渠安全；②楊莊溝洪水將于輝常路處匯入輝縣東河。在遭遇洪水時(shí)，若輝縣東河水位高，楊莊溝洪水出流不暢，甚至輝縣東河洪水會(huì)倒灌楊莊溝排水渡槽，進(jìn)而威脅總干渠安全。

為了保證總干渠安全，發(fā)生洪水時(shí)，輝縣東河暗渠交叉斷面水位不能超過(guò)設(shè)計(jì)值，同時(shí)輝常路處水位不能超過(guò)楊莊溝涵洞頂高程。

文章采用Mike11 一維水動(dòng)力模型模擬在現(xiàn)狀地形、地物條件下的輝縣東河洪水潰決演進(jìn)過(guò)程，分析洪水對(duì)總干渠及輝縣東河上下游城區(qū)可能產(chǎn)生的風(fēng)險(xiǎn)，提供相應(yīng)的工程措施建議，保證總干渠的安全，也為中線工程汛期應(yīng)急搶險(xiǎn)提供決策支撐。

2 模型分析

2.1 模型原理

MIKE 11 水動(dòng)力學(xué)模型的微分方程是一維明渠非恒定流方程，即圣維南（Saint-Venant）方程組。具體形式為：

式（1）和（2）中：t 為時(shí)間坐標(biāo)；x 為河道沿程坐標(biāo)；Q 為流量；Z 為水位；A 為過(guò)水?dāng)嗝娴拿娣e；B 為水面寬度；K 為流量模數(shù)；g為重力加速度；q為旁側(cè)入流流量。

式（3）中：n為河道糙率系數(shù)；R為水力半徑。

MIKE 11采用六點(diǎn)Abbott-Ionescu 有限差分格式對(duì)圣維南方程組進(jìn)行離散，離散后的方程組用追趕法求解。該離散格式在每一個(gè)網(wǎng)格點(diǎn)并不同時(shí)計(jì)算水位和流量，而是按順序交替計(jì)算水位或流量，分別稱(chēng)為h點(diǎn)和Q點(diǎn)，如圖1所示。

圖1 六點(diǎn)有限差分格式水位點(diǎn)、流量點(diǎn)交替布置示意圖

2.2 模型構(gòu)建

2.2.1 河道斷面

利用無(wú)人機(jī)測(cè)量了輝縣東河附近地形及河道斷面，利用航測(cè)資料，每隔100 m提取一個(gè)斷面，構(gòu)建斷面文件，生成了計(jì)算范圍0+000（上游）～2+500（下游）內(nèi)共26個(gè)斷面，其中輝縣東河暗渠位于1+000 斷面，該斷面數(shù)據(jù)采用總干渠設(shè)計(jì)斷面值。輝常路橋位于1+500斷面，幸福橋位于斷面1+700斷面。

2.2.2 河道建筑物概化

河道上共有兩個(gè)建筑物，分別為改建后的輝常路橋（3 孔7 m×6 m 的橋孔）和幸福橋（共一孔半，其中半孔不過(guò)水，單孔為城門(mén)洞型10 m×3 m）。模型設(shè)置上將兩座橋梁概化為涵洞與寬頂堰的組合建筑物，水位低的時(shí)候，水流過(guò)涵，水位過(guò)高時(shí)，變?yōu)楹吹挠袎毫?，同時(shí)出現(xiàn)堰流，模型會(huì)根據(jù)水位自動(dòng)調(diào)整采用堰流公式或者有壓流公式進(jìn)行計(jì)算。

考慮到洪水實(shí)際演進(jìn)情況，在輝常路橋和幸福橋處洪水會(huì)大量集中漫溢，所以在模型中，這兩處均設(shè)置潰口，進(jìn)行潰口潰決模擬，當(dāng)水位超過(guò)橋面高程時(shí)，洪水沿橋面向周?chē)?，隨周?chē)貏?shì)繼續(xù)演進(jìn)。

2.2.3 邊界條件

模擬采用輝縣東河遭遇5、10、20、50 和200 a 一遇共5 個(gè)洪水量級(jí)。模型共設(shè)置3 個(gè)邊界條件，上邊界0+000 處為5 個(gè)頻率對(duì)應(yīng)的洪峰流量，下邊界2+500 處為該斷面的水位流量關(guān)系，最后一個(gè)為輝常路橋處楊莊溝渡槽入流輝縣東河，該入流作為點(diǎn)源在輝常路橋處加入。

2.2.4 參數(shù)選取與模型驗(yàn)證

一維模型中涉及到的參數(shù)主要為河道糙率，根據(jù)查勘情況，結(jié)合設(shè)計(jì)階段輝縣東河糙率值，初步選取糙率值，通過(guò)模型計(jì)算不斷調(diào)試，最終確定輝縣東河主槽糙率為0.03、灘地糙率為0.06。

模型驗(yàn)證與總干渠設(shè)計(jì)標(biāo)準(zhǔn)和校核標(biāo)準(zhǔn)一致，選取50 a和200 a 一遇兩個(gè)工況，同時(shí)河道沿程不漫溢，河道建筑物也與設(shè)計(jì)階段一樣，通過(guò)模型計(jì)算出輝縣東河沿程水面線，并選取輝縣東河暗渠處洪水位與總干渠設(shè)計(jì)值對(duì)比分析，對(duì)比結(jié)果見(jiàn)表1。結(jié)果顯示，模型計(jì)算結(jié)果均高于總干渠設(shè)計(jì)值，差值較小，均不超過(guò)0.20 m，精度滿足要求，說(shuō)明模型是可靠的。

表1 模型計(jì)算結(jié)果與總干渠設(shè)計(jì)結(jié)果對(duì)比情況表

2.2.5 洪水方案

采用最不利洪水組合，輝縣東河上游洪水頻率對(duì)應(yīng)相同頻率楊莊溝涵洞入流組合，計(jì)算中有5、10、20、50 和200 a 一遇共5個(gè)頻率洪水組合。

3 模擬結(jié)果

通過(guò)模型計(jì)算，得到了現(xiàn)狀地形及行洪邊界條件下5 中工況輝縣東河沿程水面線。經(jīng)過(guò)計(jì)算結(jié)果分析，①總干渠交叉斷面處50 a 和200 a 一遇水位均低于總干渠設(shè)計(jì)階段水位成果，對(duì)總干渠安全。②只有在5 a 一遇時(shí)輝常路橋處水位低于楊莊溝涵洞頂高程106.00 m，輝縣東河洪水不會(huì)倒灌入楊莊溝，洪水頻率高于5 a 一遇就可能倒灌楊莊溝，威脅總干渠的安全。所以現(xiàn)狀邊界條件對(duì)總干渠存在很大的風(fēng)險(xiǎn)。

從計(jì)算結(jié)果來(lái)看，幸福橋處水位落差接近1 m，明顯是幸福橋的阻水作用推高了上游水位，如果將幸福橋拆除，輝常路橋處水位就會(huì)降低很多，也許可以滿足總干渠防洪要求。為了尋求更合理的方案，排出總干渠的防洪安全隱患，將上述模型調(diào)整，拆除幸福橋，模擬計(jì)算5個(gè)工況下水面線。

從調(diào)整后結(jié)果分析，幸福橋拆除后，輝常路橋處的水位降低明顯，200 a 一遇洪水時(shí)，輝常路橋水位為105.38 m，低于楊莊溝涵洞頂高程106.00 m，不會(huì)威脅到楊莊溝的安全，同時(shí)總干渠處水位也明顯降低，50 a 和200 a 一遇洪水時(shí)，總干渠處水位分別低于設(shè)計(jì)結(jié)果1.74 m和1.81 m，總干渠也是安全的。

4 結(jié)論

①為了分析輝縣東河暗渠下游邊界條件改變后，洪水對(duì)總干渠的風(fēng)險(xiǎn)，文章基于Mike11 軟件，以輝縣東河為研究對(duì)象，建立輝縣東河一維河道潰決水動(dòng)力模型，模擬計(jì)算了輝縣東河在現(xiàn)狀邊界條件下的洪水潰決演進(jìn)過(guò)程，得到了5種工況下輝縣東河沿程水面線。②通過(guò)模擬計(jì)算結(jié)果來(lái)看，現(xiàn)狀邊界條件改變后，洪水對(duì)總干渠的安全有很大威脅，風(fēng)險(xiǎn)很大。經(jīng)過(guò)調(diào)整下游邊界，模擬得出，拆除幸福橋以后，總干渠的洪水風(fēng)險(xiǎn)基本消除。為了保證總干渠安全，建議拆除幸福橋。