郭立兵 田福昌 曹魯贛 苑希民

摘要：為合理評價(jià)山洪溝道防沖建筑物壅水風(fēng)險(xiǎn)，基于Saint-Venant方程組和Villemonte堰流公式，構(gòu)建了防沖建筑物壅水風(fēng)險(xiǎn)評估水動(dòng)力模型。根據(jù)防沖建筑物的實(shí)際情況和設(shè)計(jì)要求，將防沖建筑物概化為薄壁堰。將該模型應(yīng)用于防沖建筑物分布密集的寧夏北武當(dāng)溝，分析評估梯級防沖坎壅水風(fēng)險(xiǎn)，并提出優(yōu)化措施建議。計(jì)算結(jié)果表明：所建模型能夠較精確地模擬山洪演進(jìn)過程和防沖建筑物壅水風(fēng)險(xiǎn)，充分體現(xiàn)了防沖建筑物的阻水效果。研究成果對溝道山洪的精細(xì)仿真模擬與風(fēng)險(xiǎn)評估具有一定參考價(jià)值。

關(guān)鍵詞：山洪模擬；防沖建筑物；壅水風(fēng)險(xiǎn)評估；Villemonte堰流公式；北武當(dāng)溝

中圖分類號：TV131.2 文獻(xiàn)標(biāo)志碼：A doi：10.3969/j.issn.1000-1379.2018.09.010

為了滿足河道的抗沖刷要求，通常會在河道中修建防沖坎等防沖建筑物。防沖建筑物減小了河道過流斷面面積，并增大局部水頭損失，導(dǎo)致河道水流情況更加復(fù)雜。山洪暴發(fā)時(shí)，防沖建筑物前水位壅高，將對沿岸防洪產(chǎn)生諸多不利影響，合理評估山洪溝道溢流防沖建筑物壅水風(fēng)險(xiǎn)是十分必要的。

近年來，相關(guān)學(xué)者針對河道—維水動(dòng)力模型的研究較多，成果頗豐。M.A.Nassar[1]利用一維SOBEK模型，模擬了灌溉改造工程對河道水位的影響；N.Pra-manik等[2]在缺少河道實(shí)測斷面數(shù)據(jù)的情況下，利用衛(wèi)星遙感影像建立了MIKE11模型，所建模型能夠較精確地模擬印度Brahmani河流的洪水演進(jìn)過程；Zhang X.等[3]提出了一種潮汐作用的表達(dá)式，并將潮汐作用表達(dá)式與Saint-Venant方程組結(jié)合，模擬了潮汐作用對洪水演進(jìn)的影響；朱茂森[4]將河道—維模型應(yīng)用于遼河流域的水質(zhì)計(jì)算，模擬了水體污染物的遷移、擴(kuò)散和衰減過程；郝紅升等[5]構(gòu)建了一維水動(dòng)力模型，分析了電站排水及水庫泄水對河流水環(huán)境的影響。雖然河道—維水動(dòng)力模型已得到較多應(yīng)用，但其在河道防沖建筑物壅水風(fēng)險(xiǎn)評估中的應(yīng)用研究比較少見。

本研究基于Saint-Venant方程組和Villemonte堰流公式，以寧夏北武當(dāng)溝為研究對象，建立山洪溝道防沖建筑物壅水風(fēng)險(xiǎn)評估的一維水動(dòng)力模型，將防沖建筑物概化為薄壁堰，重點(diǎn)分析梯級防沖坎對溝道洪水演進(jìn)過程的影響，評估防沖建筑物壅水風(fēng)險(xiǎn)，并提出減災(zāi)優(yōu)化措施建議。

1 模型原理

1.1 溝道—維水動(dòng)力模型

溝道—維水動(dòng)力模型是基于Saint-Venant方程組建立的，并利用Abbott六點(diǎn)隱式格式[6-7]求解該方程組，Saint-Venant方程組[8]如下：式中；B為溝道過流橫斷面寬度；z為溝道水深；Q為溝道總流量；q為溝道單寬旁側(cè)流量；t為時(shí)間；s為間距坐標(biāo)；g為重力加速度；A為溝道過流橫斷面面積；i為溝底坡降；R為水力半徑；C為謝才系數(shù)。

1.2 防沖建筑物概化與Villemonte堰流公式

針對防沖坎等防沖建筑物，在一維水動(dòng)力模型中，可以利用局部水頭損失的能量方程計(jì)算過流量[9]，也可以利用一些相對簡單的公式進(jìn)行計(jì)算，比如Villemonte堰流公式[10]。Villemonte堰流公式適用于在淹沒出流情況下薄壁堰的流量計(jì)算。在寧夏北武當(dāng)溝研究區(qū)段內(nèi)，溝道防沖坎采用鉛絲籠水平護(hù)腳，防沖坎頂部厚度為0.4m，防沖坎頂部厚度和坎頂水頭的關(guān)系通常滿足薄壁堰條件；防沖坎位于河道底部，坎頂高程和下游水位的關(guān)系通常滿足淹沒出流條件[11-12]。根據(jù)溝道防沖坎的實(shí)際情況和設(shè)計(jì)參數(shù)，本研究將防沖建筑物概化為薄壁堰。Villemonte堰流公式為式中：Qs為防沖建筑物頂部過流量；a為流量系數(shù)；W為堰頂寬度；Hu與Hd分別為防沖建筑物上游和下游水位；Hw為堰頂高程。

2 應(yīng)用實(shí)例

2.1 研究區(qū)概況

選擇防沖建筑物分布密集、具有代表性的寧夏北武當(dāng)溝為研究對象。北武當(dāng)溝位于寧夏石嘴山市大武口區(qū)，由歸德溝、韭菜溝匯合而形成。研究區(qū)段位于煤機(jī)鐵路專用線3#橋和大武口攔洪庫之間。研究區(qū)段內(nèi)溝道總長7.2km，最大寬度457m，最小寬度100m，平均比降0.71%，堤防全部為土堤，臨、背河坡比均為1：6。該區(qū)段防沖建筑物分布密集，共設(shè)置3座梯級防沖坎。北武當(dāng)溝地處石嘴山市，該地區(qū)為寧夏兩大暴雨中心之一，流域內(nèi)水土流失嚴(yán)重，山洪暴發(fā)頻繁（多出現(xiàn)在7-8月）。頻發(fā)的山洪極有可能導(dǎo)致防沖建筑物處發(fā)生洪水漫溢，威脅兩岸安全。針對該河段構(gòu)建防沖建筑物壅水風(fēng)險(xiǎn)評估水動(dòng)力模型，分析梯級防沖建筑物的壅水風(fēng)險(xiǎn)以及防沖建筑物對洪水演進(jìn)的影響。研究區(qū)域位置見圖1。

2.2 防沖建筑物壅水風(fēng)險(xiǎn)評估模型建立

2.2.1 河網(wǎng)概化與主要參數(shù)設(shè)置

北武當(dāng)溝模擬范圍為煤機(jī)鐵路專用線3#橋至大武口攔洪庫（中域），溝道全長7.2km。為了使模擬結(jié)果更加符合實(shí)際情況，利用溝道實(shí)測斷面數(shù)據(jù)，在防沖建筑物前后及溝道蜿蜒曲折的區(qū)域，對溝道斷面進(jìn)行內(nèi)插加密處理，加密之后壅水模型共設(shè)置124個(gè)溝道斷面，其中溝道實(shí)測斷面97個(gè)。溝道糙率是對溝道模型精度影響較大的參數(shù)，查閱相關(guān)地區(qū)資料[13]，設(shè)定北武當(dāng)溝溝道綜合糙率為0.035。綜合考慮模型穩(wěn)定及運(yùn)算效率等多種因素，設(shè)定北武當(dāng)溝一維水動(dòng)力模型計(jì)算步長為10s。

2.2.2 防沖建筑物概化處理

因北武當(dāng)溝上段較陡，不能滿足溝道抗沖刷要求，故結(jié)合城市功能區(qū)需要和沿線工程地質(zhì)情況，在比降較大、存在安全隱患的重點(diǎn)溝段設(shè)置了3座防沖坎，分段調(diào)整比降。防沖建筑物頂部厚度較小，根據(jù)防沖建筑物的實(shí)際情況和設(shè)計(jì)參數(shù)，在壅水風(fēng)險(xiǎn)評估模型中，考慮防沖建筑物的阻水作用，將防沖坎概化為薄壁堰。

根據(jù)實(shí)測斷面資料進(jìn)行整編處理，模型內(nèi)防沖建筑物參數(shù)設(shè)置見表1，1號防沖坎處溝道橫斷面見圖2。

2.2.3 邊界條件確定

模型上游控制邊界條件為煤機(jī)鐵路專用線3#橋斷面設(shè)計(jì)洪水過程。北武當(dāng)溝并無水文測站長期監(jiān)測資料，考慮該區(qū)域暴雨洪水特性，根據(jù)區(qū)域暴雨洪水資料[14]，利用設(shè)計(jì)暴雨推求北武當(dāng)溝設(shè)計(jì)洪水，采用地區(qū)均衡法得到3#橋斷面100a—遇設(shè)計(jì)洪水過程，見圖3。

模型下游控制邊界為大武口攔洪庫入庫斷面（見圖4），控制出流條件為該斷面水位一流量關(guān)系。采用曼寧公式推求入庫斷面水位一流量關(guān)系，斷面比降為0.0071、糙率為0.035，推求結(jié)果見圖5。

2.3 山洪計(jì)算結(jié)果與風(fēng)險(xiǎn)評估建議

2.3.1 溝道洪水計(jì)算結(jié)果

根據(jù)所建北武當(dāng)溝溝道—維水動(dòng)力模型，考慮梯級防沖坎壅水作用，模擬北武當(dāng)溝100a—遇洪水演進(jìn)過程，推算各斷面水位、流量等，提取河道計(jì)算斷面最高水位，連接各斷面最高水位，繪制河道最高水面線，見圖6。

由圖6可知，在100a—遇洪水情況下，北武當(dāng)溝洪水最高水位低于兩岸堤頂高程，但3座防沖坎壅水作用尤為明顯，最高壅水水位接近堤頂高程。如遇超百年稀遇洪水，則應(yīng)特別注意洪水漫溢淹沒影響，潛在風(fēng)險(xiǎn)較大。

2.3.2 防沖建筑物壅水風(fēng)險(xiǎn)分析

由北武當(dāng)溝溢流防沖建筑物壅水風(fēng)險(xiǎn)評估模型計(jì)算結(jié)果可知，防沖建筑物很大程度上影響了溝道行洪，導(dǎo)致溝道水位壅高。經(jīng)統(tǒng)計(jì)，北武當(dāng)溝100a—遇洪水防沖建筑物壅水水位見表2。

當(dāng)北武當(dāng)溝遭遇100a—遇洪水時(shí)，1號防沖坎最高水位1106.75m，距左堤頂2.23m，距右堤頂2.56m，溝道最大水深3.95m，坎前最大壅水高度2.75m；2號防沖坎最高水位1105.70m，距左堤頂0.23m，距右堤頂0.20m，河道最大水深3.98m，坎前最大壅水高度2.78m；3號防沖坎最高水位1103.67m，距左堤頂0.95m，距右堤頂0.93m，溝道最大水深4.07m，坎前最大壅水高度2.87m。1、2、3號防沖坎頂部厚度分別是最大坎頂水頭的0.15倍、0.14倍、0.14倍，頂部厚度小于0.67倍坎頂水頭，防沖坎頂部厚度和坎頂水頭的關(guān)系滿足薄壁堰條件，驗(yàn)證了將防沖建筑物作為薄壁堰處理的合理性。分析可知，防沖建筑物壅水效果較為明顯，壅水高度較大，導(dǎo)致溝道行洪能力減弱，將對兩岸堤防安全和沿岸防洪等產(chǎn)生不利影響。雖然坎前最高水位未超過堤頂高程，但防沖建筑物處水位與堤頂接近，屬防汛重點(diǎn)關(guān)注位置。

2.3.3 溝道防沖建筑物優(yōu)化措施建議

由山洪計(jì)算結(jié)果和防沖建筑物壅水風(fēng)險(xiǎn)分析結(jié)果可知，防沖建筑物減小了溝道過流斷面面積，壅水抬升作用較為明顯，導(dǎo)致溝道行洪能力減弱。根據(jù)溝道防沖建筑物的水位、流量等，從降低防沖建筑物壅水風(fēng)險(xiǎn)角度提出以下幾點(diǎn)優(yōu)化措施建議：①在山區(qū)溝道滿足抗沖刷要求的條件下，應(yīng)盡可能避免或減少防沖坎等防沖建筑物，包括防沖建筑物的數(shù)量和規(guī)模；②為了減弱防沖坎的阻水效果和減小銜接處上下游水位差，防沖坎坎頂石料宜采用斜面石；③因防沖坎壅水水位接近堤頂高程，故為了避免堤防在長期高水位作用下發(fā)生流土、管涌和潰堤等險(xiǎn)情，對防沖建筑物處的兩岸堤防應(yīng)進(jìn)行加高、加固處理，并在迎水坡面設(shè)置混凝土防滲墻、防滲斜板及防滲趾墻；④為了避免超百年稀遇洪水造成洪水漫溢，可在防沖建筑物的某一側(cè)或兩側(cè)修建適當(dāng)?shù)膶?dǎo)洪堤，導(dǎo)洪堤結(jié)構(gòu)可采用漿砌石重力式直墻。3結(jié)語

為了合理評估溝道防沖建筑物壅水風(fēng)險(xiǎn)，基于Saint-Venant方程組和Villemonte堰流公式，建立了適用于模擬防沖建筑物雍水過程的一維水動(dòng)力模型。針對溝道防沖建筑物的過流特點(diǎn)，將防沖建筑物概化為薄壁堰。模型應(yīng)用于防沖建筑物密集的寧夏北武當(dāng)溝防沖坎壅水分析，結(jié)果表明：在100a—遇洪水情況下，防沖建筑物壅水效果明顯，壅水高度較大，防沖建筑物處水位距堤頂較近，導(dǎo)致溝道行洪能力減弱，設(shè)置防沖建筑物的溝段屬防汛重點(diǎn)關(guān)注位置。從降低防沖建筑物塞水風(fēng)險(xiǎn)角度提出了幾點(diǎn)優(yōu)化措施建議，所建模型適用于評估山洪溝道防沖建筑物的壅水風(fēng)險(xiǎn)，也可用于溝道山洪演進(jìn)過程的精細(xì)仿真模擬。

參考文獻(xiàn)：

[1]NASSARm A.One Dimensional Hydrodynamicmodel Simu-lating Water Stage in Open Channels（WS-1）[J].Interna-tional Journal ofmodeling Simulation&ScientificComputing，2010，1（2）：303-316.

[2]PRAMANIK N，PANDA R K，SEN D.One DimensionalHydrodynamicmodeling of River Flow Using DEM ExtractedRiver Cross-Sections[J].Water Resourcesmanagement，2010，24（5）：835-852.

[3]ZRANG X，BAO W，QU S，et al.One Dimensional Hydro-dynamicmodel Accounting for Tidal Effect[J].HydrologyResearch，2012，43（1-2）：113-122.

[4]朱茂森.基于MIKE11的遼河流域一維水質(zhì)模型[J].水資源保護(hù)，2013，29（3）：6-9.

[5]郝紅升，李克鋒，莊春義.關(guān)于河道—維非恒定流水溫預(yù)測模型的研究[J].四川大學(xué)學(xué)報(bào)（自然科學(xué)版），2005，42（6）：1189-1193.

[6]ABBOTT M B.Computational Hydraulics：Elements of theTheory of Free Surface Flows[M].London：Pitman，1979：31-46.

[7]ABBOTT M B，IONESCU F.On the Numerical Computationof Nearly Horizontal Flows[J].Journal of Hydraulic Re-search，1967，5（2）：97-117.

[8]苑希民，薛文宇，馮國娜，等.潰堤洪水分析的一、二維水動(dòng)力耦合模型及應(yīng)用[J].水利水電科技進(jìn)展， 2016，36（4）：53-58.

[9]衣秀勇.DHI MIKE FLOOD洪水模擬技術(shù)應(yīng)用與研究[M].北京：中國水利水電出版社，2014：55-61.

[10]TULLIS B P，YOUNG J C，CHANDLERm A.Head-Dis-charge Relationships for Submerged Labyrinth Weirs[J].Journal of Hydraulic Engineering，2007，133（3）：248-254.

[11]彭儒武，李保棟，王青，等.矩形薄壁堰貼壁堰流試驗(yàn)研究[J].山東農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2002，33（2）：197-202.

[12]高學(xué)平，張效先.水力學(xué)[M].北京：中國建筑工業(yè)出版社，2006：236-242.

[13]寧夏回族自治區(qū)水利廳.石嘴山市大峰溝治理工程初步設(shè)計(jì)報(bào)告[R].銀川：寧夏水利廳，2014：43-65.

[14]寧夏回族自治區(qū)水文總站.寧夏回族自治區(qū)暴雨洪水圖集[M].銀川：寧夏回族自治區(qū)水文總站，1986：35-56.