滏陽新河左堤常村潰口洪水風險分析及應對研究

李會，楊翠巧，牛彥輝

（1.河北浩川工程咨詢有限公司，河北石家莊050000；2.邯鄲市水利水電勘測設計研究院，河北邯鄲056000）

滏陽新河左堤常村潰口洪水風險分析及應對研究

李會1，楊翠巧2，牛彥輝1

（1.河北浩川工程咨詢有限公司，河北石家莊050000；2.邯鄲市水利水電勘測設計研究院，河北邯鄲056000）

考慮地形、堤防工程、道路、橋涵、房屋、河流對洪水演進的影響，建立了滏陽新河MIKE11模型、堤外保護區MIKE21模型及一二維耦合MIKE FLOOD模型，分析了50 a一遇、“63·8”校核洪水常村段堤防潰決后的洪水演進過程，結合淹沒區社會經濟數據進行了洪災影響分析，繪制了該區域洪水風險圖，提出減小洪災的應對措施，為防汛指揮及決策提供科學依據。

洪水演進模擬；洪災影響；洪水風險圖；滏陽新河

1 引言

滏陽新河是河北省為避免滏陽河流域發生類似1963年型洪水而開挖的人工排洪河道，起自艾辛莊樞紐，入獻縣泛區與滏陽河、滹沱河匯合，河道長度132.4 km，設計標準為50 a一遇，校核標準為“63·8”洪水。左、右堤均為2級提防，左堤長127.39 km，右堤長128.8 km，兩堤間距平均為1 500 m，現狀堤頂寬6～8 m，邊坡1∶3。

滏陽新河左堤常村險工位于衡水市武強縣豆村鄉常村南，寬度約240 m，堤防樁號107+660—107+ 900，堤身以壤土為主，局部夾砂壤土，外堤腳緊鄰夾道排水渠，堤頂與渠底高差達9 m，堤身相對變高，抗滑能力降低，且堤防外坡面不平整，局部有沖溝，行洪期間極易形成潰口。堤外保護區為滹滏區間洪澇災害最為頻繁區域，涉及衡水市武強縣、武邑縣，是河北省防洪除澇的重點區域。

2 研究目標及基礎資料

2.1 研究目標

主要目標是建立滏陽新河MIKE11模型、堤外保護區MIKE21模型及一二維耦合MIKE FLOOD模型，模擬滏陽新河50 a一遇、“63·8”校核洪水常村段堤防潰決后的洪水演進過程，評估淹沒區的洪水影響及損失，繪制洪水風險圖并提出減小洪災的應對措施。

2.2 分析范圍

分析范圍為滏陽新河左堤保護區，保護區為滏陽新河左堤以西區域，南起寧晉泊北圍堤，北至滹沱河南大堤，西至30 m等高線，涉及面積3 973 km2。

2.3 基礎資料

（1）地形資料。主要包括滏陽新河河道縱橫斷面資料、分析范圍內1∶10 000地形圖、常村潰口附近1∶2 000地形圖。

（2）水文資料。滏陽新河50 a一遇、“63·8”校核洪水流量過程，如圖1所示；保護區內主要排瀝河道滏陽河及其支溝汪洋溝、邵村溝、胡堂排干、小西河、龍治河、天平溝、留楚排干20 a一遇瀝水過程，以較大支溝天平溝為例，如圖2所示；獻縣泛區水位過程線，如圖3所示。

（3）防洪工程。收集了保護區防洪工程分布情況，測量了滏陽新河、滏陽河、龍治河、天平溝等河道堤防及滹沱河南大堤的高程資料。

（4）交通道路。測量并收集整理了分析范圍內阻水交通干線如大廣高速、黃石高速、青銀高速、邢衡高速、京九鐵路、石德鐵路、G307、G308的路面高程、過水橋涵尺寸。

（5）社會經濟。根據各縣市2013年經濟統計年鑒，以鄉鎮為單元統計了區域社會經濟指標，包括面積、地區生產總值、常住人口、耕地面積等指標。

圖1 滏陽新河設計洪水過程線

圖2 保護區排瀝河道天平溝設計瀝水過程線

圖3 獻縣泛區水位過程線

3 洪水演進分析

3.1 洪水分析模型及構建思路

洪水分析模型選用丹麥水資源及水環境研究所開發的MIKE系列模型——MIKE11、MIKE21、MIKE FLOOD。該軟件目前廣泛應用于世界各國有關水力工程的分析計算，計算的可靠性和準確性可以滿足本次洪水分析的要求。

將滏陽新河河道入流過程加入MIKE11模型，建立滏陽新河一維河道模型；利用MIKE21模型建立保護區二維非恒定流洪水演進模型，區內各排澇河道澇水概化為點源加入到各河道入滏陽河口處；通過滏陽新河左堤常村潰口，利用MIKE FLOOD將一維河道模型與二維保護區模型進行耦合，實現河道洪水傳播與保護區洪水演進的同步交換與實時模擬。

為充分考慮保護區內線狀構筑物對洪水演進的影響，對高于地面0.5 m以上的線狀構筑物，如道路、堤防，以堤壩（Dikes）的形式添加至模型中，并將每個線性地物在橋涵或不連續堤防位置斷開。

3.2 網格剖分

保護區二維計算采用MIKE21 FM不規則網格剖分，并將主要阻水建筑物加入網格生成器中參與網格生成，保證構筑物落在不規則三角網的一個邊上。本次共剖分網格101 117個，計算節點51 155個，主要分為3個剖分區域。

（1）滏陽河。滏陽河為保護區內的主要排澇河道，為提高模型計算的準確性和合理性，將河道內作局部網格加密處理，最大網格面積控制在0.02 km2以內。

（2）保護區東側。該區域緊鄰滏陽新河左堤，受洪水威脅較大，將此區域網格適當加密，最大網格面積控制在0.05 km2以內。

（3）保護區西側。該區域距離滏陽新河較遠，且保護區西高東低，經模型初步分析模擬后，受洪水威脅較小，故適當加大網格面積，最大網格面積控制在0.1 km2以內。

網格剖分結果，如圖4所示。

圖4 網格剖分結果

3.3 參數選取

因滏陽新河無實測歷史洪水資料，無法進行參數率定與驗證。本次依據《子牙河流域防洪規劃報告》以及以往經驗成果綜合考慮，滏陽新河主槽糙率選取0.022 5，灘地糙率選取0.05。保護區內無實測潰堤洪水資料，糙率值按地形、地貌條件和地面特點[1]及類似項目進行取值，確定村莊糙率為0.1、林地糙率為0.065、耕地糙率為0.055。

洪水模擬時段設置為2016年3月1日至4月15日，共計45 d，主時間步長取15 s，共計259 200個步長單位。

選擇Smagorinsky系數方法設定渦粘系數，程序缺省值為0.28。

常村險工決口方式采用潰決方式，過程設定為漸潰，潰口發展時間設定為3 h，潰口深度選取潰至堤防兩側地面較高者，潰決時機為滏陽新河50 a一遇設計水位。

3.4 邊界條件確定

模型上邊界條件為滏陽新河50 a一遇、“63·8”校核洪水流量過程及保護區各排瀝河道20 a一遇瀝水過程，下邊界為獻縣泛區水位過程。

3.5 洪水模擬成果

滏陽新河50 a一遇標準洪水時，模擬結果顯示，在第5天左右，潰口處河道水位漲至潰決水位，堤防潰決，2 h后潰口流量達到最大為1 621.2 m3/s，之后潰口流量逐漸變小，至第12天逐漸變為0 m3/s，后因退水較慢，潰口出流變為負值，即保護區內少量洪水返入滏陽新河，最終潰口出流逐漸趨于0 m3/s。洪水主要在滏陽新河左堤與滏陽河之間的夾道內行洪，沿滏陽河方向向下游演進，直至滹沱河南大堤。部分洪水漫過滏陽河堤防，向西北方向演進，直至天平溝右堤。

滏陽新河“63·8”校核洪水時，在模型模擬的第3天左右，潰口處河道水位漲至潰決水位，堤防潰決，65 h后潰口流量達到最大為2 271.0 m3/s，后期亦存在保護區內洪水返入滏陽新河現象。洪水演進過程起初與50 a一遇情形基本一致，部分洪水在夾道內行洪，部分洪水漫過滏陽河堤防，向西北方向演進，至天平溝位置受到堤防阻擋，當洪水達到一定高度時漫過天平溝堤防，向北演進，直至滹沱河南大堤。

常村潰口流量過程線，如圖5所示。

圖5 常村潰口流量過程線

3.6 洪災影響分析

根據洪水風險計算得到的淹沒范圍、淹沒水深等要素，結合保護區內各鄉鎮的社會經濟情況，綜合分析評估洪水影響及損失情況。

50 a一遇標準時，保護區受淹面積為229.34 km2，淹沒農田面積21 437.45 hm2，淹沒房屋面積1 288.34萬m2，受影響人口總數10.55萬人，受影響GDP為181 583.04萬元。

“63·8”校核洪水時，保護區受淹面積為641.77km2，淹沒農田面積59 985.51 hm2，淹沒房屋面積3 511.63萬m2，受影響人口總數29.18萬人，受影響GDP為531 682.36萬元。

3.7 洪水風險圖繪制

洪水風險圖可直觀反映某個區域發生洪水后可能淹沒的范圍、水深等洪水要素[2]。本次選用基于ArcGIS平臺的中國水利水電科學研究院洪水風險圖繪制系統，繪制滏陽新河左堤防洪保護區的洪水

圖6 滏陽新河“63·8”校核洪水常村潰口洪水風險（淹沒水深）

風險圖[3]。風險圖顯示內容主要包括：①地理底圖，包括居民地、主要河湖、渠道、主要交通道路、鄉鎮及以上行政界等信息；②水利工程要素，主要包括控制站、堤防、水閘、險工險段信息等；③洪水風險專題信息，包括淹沒水深、洪水流速、到達時間信息等。

將上述內容放入繪制系統，并增加圖名、圖例、潰口流量過程線、淹沒損失情況及方案說明等內容，繪制洪水風險圖。其中，滏陽新河“63·8”校核洪水常村潰口洪水風險（淹沒水深），如圖6所示。

4 應對措施

洪水進入保護區后，主要通過滏陽河退入獻縣泛區。因滏陽河過流能力較低，且受滹沱河南大堤的阻擋，大量洪水滯留在保護區，致使洪災損失較大。為加速退水，建議在滏陽河入獻縣泛區位置相機扒開滹沱河南大堤。

本次按滹沱河南大堤扒口寬度200 m，結合滏

TV122+.4

A

1004-7328（2017）04-0030-04

10.3969/j.issn.1004-7328.2017.04.010

2017—03—12

李會（1987—），女，工程師，主要從事水文規劃、防洪影響評價工作。