杜 麥，陳小威

（1.石家莊豐水工程咨詢(xún)有限公司，河北 石家莊 050000；2.中國(guó)電建集團(tuán)華東勘測(cè)設(shè)計(jì)研究院有限公司，浙江 杭州 311122）

1 引言

近年來(lái)，由于河道管理范圍不明確，四亂等行為時(shí)有發(fā)生，影響河道行洪及人民群眾生命財(cái)產(chǎn)安全[1］，為此對(duì)潮河四道河段河道管理范圍進(jìn)行劃界對(duì)其具有重要意義。

依據(jù)《河北省河湖管理范圍復(fù)核及劃定技術(shù)指南》，無(wú)堤段河道管理范圍原則上按照歷史最高洪水位或設(shè)計(jì)洪水位確定。由于歷史最高洪水位難以確定，本次按照設(shè)計(jì)洪水位確定其管理范圍，故河道水面線(xiàn)推求是河道劃界的重要內(nèi)容和基礎(chǔ)，其推算結(jié)果將直接影響管理范圍的劃定。常規(guī)的逐段試算法，對(duì)于流態(tài)復(fù)雜的山區(qū)河道，極易出現(xiàn)不合理或無(wú)解現(xiàn)象[2］。HEC-RAS（Hydrologic Engineering Center’s River Analysis System）模型是由美國(guó)陸軍工程師團(tuán)開(kāi)發(fā)的一款模型，該模型在山區(qū)縱坡較陡的天然河道水面線(xiàn)推求中具有良好的模擬成果[3］，本文利用HEC-RAS 建立一維水力模型對(duì)潮河四道河段進(jìn)行不同流態(tài)下設(shè)計(jì)洪水位推求。

2 計(jì)算原理

HEC-RAS 可進(jìn)行一維恒定流、一維非恒定流模擬，可用于河道坡度不大于10%的天然河道的水面線(xiàn)推求，計(jì)算緩流、急流及混合流等多種流態(tài)[4］。

2.1 恒定流模擬

基于一維能量方程，通過(guò)直接步進(jìn)方法推求水面線(xiàn)。計(jì)算公式如下：

式中：Zu、Zd分別為上、下游斷面水位；au、ad分別為上、下游斷面流速系數(shù)；vu、vd分別為上、下游斷面平均流速；hf、hj分別為上下游斷面之間的沿程水頭損失和局部水頭損失。

2.2 非恒定流模擬

基于一維連續(xù)性方程和動(dòng)量方程，計(jì)算原理如下：

式中：ρw為水的密度；vi、vj為斷面流速；xi、xj為距離；p為壓力；fi為質(zhì)量力；v 為流體的運(yùn)動(dòng)粘滯系數(shù)。

3 模型構(gòu)建

以河北省豐寧滿(mǎn)族自治縣潮河四道河段為研究對(duì)象，建立一維水力模型推求設(shè)計(jì)洪水水面線(xiàn)。

3.1 研究區(qū)域概況

潮河為潮白河支流，屬北三河水系，發(fā)源于豐寧滿(mǎn)族自治縣哈拉海灣村，于天橋鎮(zhèn)的前溝門(mén)村入灤平縣境內(nèi)，在古北口鎮(zhèn)入北京市密云縣，最終匯入密云水庫(kù)。潮河流域面積6920 km2，其中豐寧縣境內(nèi)流域面積3336.01 km2。潮河于四道河村建設(shè)有大閣水文站，建站時(shí)間為1956 年，控制流域面積為1850 km2。

根據(jù)豐寧大閣氣象站的統(tǒng)計(jì)數(shù)據(jù)，流域內(nèi)多年平均氣溫6.4℃，一月份平均氣溫-11.7℃，七月份平均氣溫22.1℃；全面無(wú)霜期129 d，多年平均降水量為463 mm；多年平均凈流深60 mm；多年平均水面蒸發(fā)量為950 mm，多年平均陸面蒸發(fā)量為400 mm；年平均風(fēng)速2.1 m／s；年平均日照時(shí)數(shù)為2802 h；多年平均凍土厚度約1.5 m。

本次研究區(qū)域?yàn)槌焙铀牡篮佣危芯繀^(qū)域以上流域面積為1862 km2，該段河長(zhǎng)約4.6 km，坡降約9.7‰。該段河道為天然河道，兩岸無(wú)堤防及護(hù)岸，現(xiàn)狀主槽寬15 m～45 m，左右岸為農(nóng)田或山體，于四道河處建有大閣水文站。

3.2 參數(shù)設(shè)置

潮河四道河段為山區(qū)河道，兩岸為農(nóng)田和村莊，防洪標(biāo)準(zhǔn)取10 年一遇（P=10%），依據(jù)《豐寧滿(mǎn)族自治縣潮河防洪整治規(guī)劃》，該段10 年一遇設(shè)計(jì)洪水為689 m3／s，計(jì)算采用一維恒定流，分別采用緩流流態(tài)和急流流態(tài)推算兩種工況的設(shè)計(jì)洪水位。

（1）河道斷面

采用2018 年實(shí)測(cè)斷面，該段共計(jì)28 個(gè)大斷面，間距約30 m～220 m。

（2）河道糙率

根據(jù)現(xiàn)場(chǎng)查勘的河道平面形態(tài)、橫斷面形態(tài)及河床質(zhì)粒徑，參考《水力計(jì)算手冊(cè)》和以往工程，結(jié)合水力學(xué)參數(shù)綜合比選，河道灘地糙率采用0.05，河道主槽糙率采用0.035。

（3）起推斷面

急流工況下，通過(guò)敏感性分析，起推斷面位于樁號(hào)0+000斷面上游約1 km 處時(shí)，對(duì)該段水位影響較小，故起推斷面位于樁號(hào)0+000 斷面上游約1 km 處，則終點(diǎn)斷面為樁號(hào)4+605斷面；同理，緩流工況下，起推斷面位于樁號(hào)4+605 斷面下游約2 km 處，終點(diǎn)斷面為樁號(hào)0+000 斷面。

（4）起推水位

起推水位根據(jù)明渠均勻流計(jì)算。

4 計(jì)算成果與分析

4.1 計(jì)算成果

根據(jù)以上參數(shù)對(duì)四道河段10 年一遇洪水位進(jìn)行推求，推求成果如表1 所示。由模擬成果可知，急流工況下水位較低，流速較大，緩流工況計(jì)算的水位較高而流速較小。

表1 潮河四道河段10 年一遇水面線(xiàn)

續(xù)表1

4.2 合理性分析

該河段建有大閣水文站，于1998 年和2011 年發(fā)生較大洪水。

為驗(yàn)證一維水力模型的合理性，采用1998 年、2011 年實(shí)測(cè)洪水，利用建立的一維模型模擬驗(yàn)證，將模擬成果與實(shí)測(cè)成果進(jìn)行對(duì)比，見(jiàn)圖1 和圖2。

圖1 1998 年洪水實(shí)測(cè)水位與模擬水位對(duì)比

圖2 2011 年洪水實(shí)測(cè)水位與模擬水位對(duì)比

由圖可知，兩種工況下模擬成果與實(shí)測(cè)成果趨勢(shì)相同，模擬水位均大于實(shí)測(cè)水位，但急流工況成果更加接近實(shí)測(cè)水位，主要原因?yàn)樵摱魏拥罏樯絽^(qū)河道，坡度較陡，采用急流流態(tài)模擬更加符合實(shí)際情況[5］。

1998 年洪水條件下，模擬成果略大于實(shí)測(cè)水位，流量越大，差距越大；2011 年洪水條件下，流量較小時(shí)，兩種工況下模擬成果均與實(shí)測(cè)成果相近，流量較大時(shí)，緩流工況下模擬成果較實(shí)測(cè)水位增幅較大，急流工況下較實(shí)測(cè)值有所增加。原因主要有以下兩點(diǎn)：其一，模型采用2018 年實(shí)測(cè)斷面，相比1998 年和2011 年，河道產(chǎn)生了淤積，導(dǎo)致模擬成果與實(shí)測(cè)成果相比偏大；其二，流量較大時(shí)，洪水出槽，由于本次測(cè)量斷面偏窄，使模擬結(jié)果產(chǎn)生誤差，相比實(shí)測(cè)成果有所增大。

綜上，通過(guò)建立的一維水動(dòng)力模型模擬結(jié)果可知，與緩流工況相比，急流工況成果更加符合實(shí)際，雖在流量較大時(shí)急流工況計(jì)算成果與實(shí)際成果相比偏大，但為確保兩岸農(nóng)田及村莊安全，為河道管理留足余地，同時(shí)考慮近年河道淤積影響，采用該模型急流工況成果進(jìn)行河道劃界是合適的。

5 管理范圍劃定

根據(jù)《河北省河湖管理范圍復(fù)核及劃定技術(shù)指南》，無(wú)堤段劃定原則如下：

無(wú)堤防行洪河道，其管理范圍原則上按照歷史最高洪水位或設(shè)計(jì)洪水位確定。寬淺式、無(wú)明顯河型、出山口、過(guò)村莊等河段，以滿(mǎn)足河道行洪能力、保護(hù)對(duì)象防洪安全為原則確定，具體經(jīng)過(guò)有資質(zhì)的單位分析論證。對(duì)于已批復(fù)河道防洪整治規(guī)劃、實(shí)施方案，且河道治理工程正在實(shí)施的，按批復(fù)的管理范圍劃定。

本次潮河四道河段兩岸均無(wú)堤防，河道管理范圍按現(xiàn)狀10 年一遇淹沒(méi)范圍劃定，依據(jù)沿河道實(shí)測(cè)的1∶2000 地形圖和10 年一遇洪水位推算成果，最終得到潮河四道河段管理范圍（由于篇幅限制，成果不在計(jì)列）。

6 結(jié)論

（1）利用HEC-RAS 軟件建立潮河四道河一維水力模型，對(duì)該河段進(jìn)行急流和緩流兩種工況下洪水位推算，由模擬成果知，急流工況水位較低，流速較大，緩流工況水位較高而流速較小。

（2）將建立的一維模型通過(guò)1998 年和2011 年實(shí)測(cè)洪水進(jìn)行驗(yàn)證，發(fā)現(xiàn)山區(qū)河道推求水面線(xiàn)時(shí)，采用急流工況即從上而下推求成果更加符合實(shí)際。

（3）急流工況下模擬水位較實(shí)際水位成果偏大，但可形象的反映洪水過(guò)程特性，同時(shí)可確保兩岸農(nóng)田和村莊安全，保障河道行洪安全，為河道管理留下充足空間，考慮近年河道淤積影響，故采用該模擬成果可作為山區(qū)河道管理范圍劃定的依據(jù)。

（4）依據(jù)實(shí)測(cè)地形圖和急流工況模擬成果，最終可得到潮河四道河段管理范圍。