李會，楊翠巧，牛彥輝

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

滏陽新河左堤常村潰口洪水風(fēng)險分析及應(yīng)對研究

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

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

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

洪水演進(jìn)模擬；洪災(zāi)影響；洪水風(fēng)險圖；滏陽新河

1 引言

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

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

2 研究目標(biāo)及基礎(chǔ)資料

2.1 研究目標(biāo)

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

2.2 分析范圍

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

2.3 基礎(chǔ)資料

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

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

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

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

（5）社會經(jīng)濟(jì)。根據(jù)各縣市2013年經(jīng)濟(jì)統(tǒng)計(jì)年鑒，以鄉(xiāng)鎮(zhèn)為單元統(tǒng)計(jì)了區(qū)域社會經(jīng)濟(jì)指標(biāo)，包括面積、地區(qū)生產(chǎn)總值、常住人口、耕地面積等指標(biāo)。

圖1 滏陽新河設(shè)計(jì)洪水過程線

圖2 保護(hù)區(qū)排瀝河道天平溝設(shè)計(jì)瀝水過程線

圖3 獻(xiàn)縣泛區(qū)水位過程線

3 洪水演進(jìn)分析

3.1 洪水分析模型及構(gòu)建思路

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

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

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

3.2 網(wǎng)格剖分

保護(hù)區(qū)二維計(jì)算采用MIKE21 FM不規(guī)則網(wǎng)格剖分，并將主要阻水建筑物加入網(wǎng)格生成器中參與網(wǎng)格生成，保證構(gòu)筑物落在不規(guī)則三角網(wǎng)的一個邊上。本次共剖分網(wǎng)格101 117個，計(jì)算節(jié)點(diǎn)51 155個，主要分為3個剖分區(qū)域。

（1）滏陽河。滏陽河為保護(hù)區(qū)內(nèi)的主要排澇河道，為提高模型計(jì)算的準(zhǔn)確性和合理性，將河道內(nèi)作局部網(wǎng)格加密處理，最大網(wǎng)格面積控制在0.02 km2以內(nèi)。

（2）保護(hù)區(qū)東側(cè)。該區(qū)域緊鄰滏陽新河左堤，受洪水威脅較大，將此區(qū)域網(wǎng)格適當(dāng)加密，最大網(wǎng)格面積控制在0.05 km2以內(nèi)。

（3）保護(hù)區(qū)西側(cè)。該區(qū)域距離滏陽新河較遠(yuǎn)，且保護(hù)區(qū)西高東低，經(jīng)模型初步分析模擬后，受洪水威脅較小，故適當(dāng)加大網(wǎng)格面積，最大網(wǎng)格面積控制在0.1 km2以內(nèi)。

網(wǎng)格剖分結(jié)果，如圖4所示。

圖4 網(wǎng)格剖分結(jié)果

3.3 參數(shù)選取

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

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

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

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

3.4 邊界條件確定

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

3.5 洪水模擬成果

滏陽新河50 a一遇標(biāo)準(zhǔn)洪水時，模擬結(jié)果顯示，在第5天左右，潰口處河道水位漲至潰決水位，堤防潰決，2 h后潰口流量達(dá)到最大為1 621.2 m3/s，之后潰口流量逐漸變小，至第12天逐漸變?yōu)? m3/s，后因退水較慢，潰口出流變?yōu)樨?fù)值，即保護(hù)區(qū)內(nèi)少量洪水返入滏陽新河，最終潰口出流逐漸趨于0 m3/s。洪水主要在滏陽新河左堤與滏陽河之間的夾道內(nèi)行洪，沿滏陽河方向向下游演進(jìn)，直至滹沱河南大堤。部分洪水漫過滏陽河堤防，向西北方向演進(jìn)，直至天平溝右堤。

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

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

圖5 常村潰口流量過程線

3.6 洪災(zāi)影響分析

根據(jù)洪水風(fēng)險計(jì)算得到的淹沒范圍、淹沒水深等要素，結(jié)合保護(hù)區(qū)內(nèi)各鄉(xiāng)鎮(zhèn)的社會經(jīng)濟(jì)情況，綜合分析評估洪水影響及損失情況。

50 a一遇標(biāo)準(zhǔn)時，保護(hù)區(qū)受淹面積為229.34 km2，淹沒農(nóng)田面積21 437.45 hm2，淹沒房屋面積1 288.34萬m2，受影響人口總數(shù)10.55萬人，受影響GDP為181 583.04萬元。

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

3.7 洪水風(fēng)險圖繪制

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

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

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

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

4 應(yīng)對措施

洪水進(jìn)入保護(hù)區(qū)后，主要通過滏陽河退入獻(xiàn)縣泛區(qū)。因滏陽河過流能力較低，且受滹沱河南大堤的阻擋，大量洪水滯留在保護(hù)區(qū)，致使洪災(zāi)損失較大。為加速退水，建議在滏陽河入獻(xiàn)縣泛區(qū)位置相機(jī)扒開滹沱河南大堤。

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

TV122+.4

A

1004-7328（2017）04-0030-04

10.3969/j.issn.1004-7328.2017.04.010

2017—03—12

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