馬靜輝
(河北省水利水電勘測(cè)設(shè)計(jì)研究院,天津 300250)
七里莊西溝排水倒虹吸下游防洪工程為南水北調(diào)京石段防洪影響處理工程(保定市段)項(xiàng)目之一。
七里莊西溝為兩溝合一型河溝,溝溝合并前主溝50年一遇洪水洪峰流量47m3/s,倒虹吸工程建成后50年一遇洪水洪峰流量53m3/s。海漫段底高程56.0m,溝底地面高程56~57m,河溝兩側(cè)地面高程61.5m。倒虹吸下游40m為北易水灌渠,渠頂高程61.6~62.0m,渠底高程60.8m。總干渠左右堤頂高程均為64.24m。
北易水灌渠本段基本與總干渠平行,七里莊西溝穿易水干渠交叉建筑物為單孔1.1m×1.3m(寬×高)漿砌石拱涵結(jié)構(gòu)。
為保證洪水順利下泄,按50年一遇洪水標(biāo)準(zhǔn)重建穿北易水灌渠涵洞。涵洞采用2孔一聯(lián)鋼筋混凝土結(jié)構(gòu),單孔斷面尺寸2.5m×2.5m,長(zhǎng)20m。北易水灌渠涵洞軸線(xiàn)兩側(cè)各10m范圍采用0.3m厚漿砌石防護(hù)。涵洞進(jìn)出口各20m范圍河道,采用0.3m厚漿砌石進(jìn)行全斷面襯砌。
七里莊西溝排水倒虹吸下游防洪工程是擴(kuò)建穿北易水灌渠涵洞,當(dāng)發(fā)生50年一遇洪水標(biāo)準(zhǔn)時(shí),保證洪水順利下泄,保護(hù)北易水灌渠的安全。防洪工程實(shí)施時(shí),地方反映北易水灌渠現(xiàn)已無(wú)灌溉功能,現(xiàn)狀處于基本廢棄狀態(tài),因此建議取消防洪工程。
為保障南水北調(diào)總干渠安全,在防洪工程不實(shí)施的情況下,按現(xiàn)狀條件對(duì)倒虹吸進(jìn)出口水位重新進(jìn)行驗(yàn)算、復(fù)核。
為分析七里莊西溝防洪影響工程對(duì)七里莊西溝行洪及南水北調(diào)中線(xiàn)總干渠安全造成影響,根據(jù)區(qū)域地形等相關(guān)條件,采用二維水動(dòng)力DHI MIKE21模型進(jìn)行計(jì)算,模擬現(xiàn)狀條件下七里莊西溝的行洪過(guò)程。
3.1.1 二維模型建立
本次建立模型范圍主要為七里莊西溝、七里莊溝附近地區(qū),邊界由南水北調(diào)中線(xiàn)總干渠、高易鐵路、供電公司宅基、七里莊村南側(cè)田間路等組成。南部以高易鐵路為界,西起七里莊西溝以西250m,東至易縣供電公司以南,長(zhǎng)度1200m;北部以南水北調(diào)中線(xiàn)總干渠為界,西起七里莊西溝以西250m,東至易縣供電公司,長(zhǎng)度1000m;西部以田間路為界,長(zhǎng)度200m;東部以易縣供電公司和九江石材廠房基為界,長(zhǎng)400m。模型范圍總面積0.3km2。
3.1.2 地形剖分及糙率采用
3.1.2.1 地形剖分
對(duì)模擬區(qū)域進(jìn)行非結(jié)構(gòu)三角形網(wǎng)格剖分,網(wǎng)格高程代表網(wǎng)格范圍的地形平均高程,而方程解則是網(wǎng)格中心的特征值。三角形網(wǎng)格邊長(zhǎng)不超過(guò)25m,最大網(wǎng)格面積不超過(guò)200m2;模型區(qū)域共剖分了4747個(gè)網(wǎng)格單元。
3.1.2.2 糙率
結(jié)合實(shí)際并參照以往相關(guān)設(shè)計(jì)文件采用糙率情況,河槽糙率取0.025,農(nóng)田、道路糙率取0.05。
3.1.3 內(nèi)部邊界條件
計(jì)算區(qū)域內(nèi)高出地面的堤防及村莊房屋等對(duì)洪水起到阻擋作用,在計(jì)算區(qū)內(nèi)部形成阻水邊界。模型中堤防圍埝按兩種類(lèi)型處理,一是可潰決堤防,當(dāng)擋水位與堤頂相平時(shí),堤防自動(dòng)潰決;二是非潰決堤防,當(dāng)水位高于堤頂時(shí),洪水漫流。
本次模擬中設(shè)置北易水灌渠堤防兩處,分別位于北易水灌渠兩岸,均屬于非潰決堤防。
3.1.4 計(jì)算成果
通過(guò)分析計(jì)算,七里莊西溝洪水受北易水灌渠堤防頂托作用,漫出河槽向東南方向演進(jìn),一部分洪水漫過(guò)廢棄的高易鐵路下泄,一部分繼續(xù)向東通過(guò)七里莊溝鐵路橋下泄。七里莊西溝排水倒虹吸出口50年一遇設(shè)計(jì)水位61.56m,200年一遇設(shè)計(jì)水位61.70m。
3.2.1 計(jì)算方法
倒虹吸進(jìn)口水位,采用SL482—2011《灌溉與排水渠系建筑物設(shè)計(jì)規(guī)范》附錄C公式,根據(jù)出口水位進(jìn)行計(jì)算。
倒虹吸管過(guò)流能力采用SL482—2011《灌溉與排水渠系建筑物設(shè)計(jì)規(guī)范》公式計(jì)算。
式中 Q為管道的泄流量;μ為流量系數(shù);ω為管道出口的橫斷面積;ΔZ2為進(jìn)口漸變段末端斷面至管道出口斷面的水面落差;g為重力加速度;ω2為進(jìn)口漸變段末端的過(guò)水面積;Ci為管身計(jì)算段水流的舍齊系數(shù);Ri為管身計(jì)算段水流的水力半徑;R為斷面水力半徑;Li為管身計(jì)算段水流的管長(zhǎng);ωi為管身計(jì)算段水流的橫斷面積;n為邊壁糙率;v為管道的斷面平均流速;v2為進(jìn)口漸變段末斷面的平均流速;ξi為某一局部水頭損失系數(shù)。
3.2.2 進(jìn)口水位復(fù)核
根據(jù)倒虹吸出口水位及結(jié)構(gòu)尺寸,計(jì)算倒虹吸進(jìn)口50年一遇設(shè)計(jì)洪水位62.67m,200年一遇校核洪水位63.79m。
為防止河道洪水漫溢進(jìn)入渠道,總干渠防洪標(biāo)準(zhǔn)與相鄰排水建筑物的洪水標(biāo)準(zhǔn)相同。七里莊西溝流域面積2.55km2,小于20km2,屬于小河溝渠段,左排及總干渠渠道按50年一遇洪水設(shè)計(jì),200年一遇洪水校核。
小型河溝范圍內(nèi)的總干渠渠道防洪堤頂高程應(yīng)滿(mǎn)足以下要求,則:
式中 Z設(shè)為渠堤外設(shè)計(jì)洪水位,即50年一遇洪水位;Z校為渠堤外校核洪水位,即200年一遇洪水位。
根據(jù)倒虹吸進(jìn)出口50年一遇設(shè)計(jì)水位和200年一遇校核水位計(jì)算結(jié)果,按照小型河溝范圍內(nèi)的渠道防洪堤堤頂高程確定原則,設(shè)計(jì)水位加1m,校核水位加0.5m,取兩者較大者。
根據(jù)計(jì)算,七里莊西溝排水倒虹吸左堤堤頂高程64.29m,右堤堤頂高程62.56m,此處總干渠左右堤頂高程均為64.24m。
(1)取消七里莊西溝排水倒虹吸下游防洪工程后,考慮現(xiàn)狀條件下七里莊西溝的行洪過(guò)程,北易水灌渠底部單孔1.1m×1.3m涵洞行洪通暢,總干渠右堤滿(mǎn)足安全加高要求,左堤200年一遇設(shè)計(jì)水位63.79m,加上安全超高0.5m后64.29m,僅僅與現(xiàn)狀堤頂高程64.24m高出5cm,洪水不漫堤,基本對(duì)總干渠安全影響不大。
(2)采用二維水動(dòng)力DHI MIKE21模型進(jìn)行計(jì)算時(shí),內(nèi)部邊界條件計(jì)算區(qū)域內(nèi),模擬中設(shè)置北易水灌渠堤防屬于非潰決堤防。北易水灌渠堤防高程由西往東逐漸降低,由62.0m逐漸降到61.0m,由于此區(qū)域內(nèi)200年一遇洪水位為61.7m,北易水灌渠堤防為土質(zhì)結(jié)構(gòu),洪水漫頂過(guò)流后,灌渠堤防會(huì)逐漸產(chǎn)生潰決,隨著潰決加大,倒虹吸出口水位降低,滿(mǎn)足防洪超高要求。