趙 威
(湖北省電力勘測(cè)設(shè)計(jì)院有限公司,湖北 武漢 430040)
我國(guó)是暴雨洪水災(zāi)害較多的國(guó)家,洪水災(zāi)害是主要的自然災(zāi)害之一,蓄滯洪區(qū)是江河防洪體系的重要組成部分,是減輕災(zāi)害、保障防洪安全的有效措施。為了保證重點(diǎn)地區(qū)的防洪安全,將有條件地區(qū)開(kāi)辟為蓄滯洪區(qū),有計(jì)劃地蓄滯洪水,是兼顧區(qū)域防洪規(guī)劃實(shí)際與經(jīng)濟(jì)可持續(xù)發(fā)展的需要,也是為保全大局,而不得不犧牲局部利益的全局性考量。
江河防洪調(diào)度的一般原則為:充分利用河道的泄洪能力,適時(shí)發(fā)揮干支流水庫(kù)的攔洪錯(cuò)峰作用;當(dāng)控制站水位接近分洪水位時(shí),首先扒開(kāi)洲灘民垸,擴(kuò)大河道行洪能力;遇超過(guò)河道泄洪能力洪水時(shí),應(yīng)立即運(yùn)用蓄滯洪區(qū)分蓄超額洪水,確保重點(diǎn)堤防、地區(qū)和城市的防洪安全。以湖北省武漢市為例,作為長(zhǎng)江流域的重要城市,蓄滯洪區(qū)的建設(shè)對(duì)于保障地區(qū)防洪安全有著重要的意義。
在蓄滯洪區(qū)內(nèi)開(kāi)展電力工程建設(shè)需要明確該蓄滯洪區(qū)的分洪運(yùn)用概率,并根據(jù)工程類(lèi)別和相應(yīng)防洪要求進(jìn)行初步分析。圖1 為蓄洪區(qū)分洪。
圖1 蓄洪區(qū)分洪
一般蓄滯洪區(qū)的啟用概率低于100 年/次,例如白潭湖蓄滯洪區(qū)的啟用概率為20~30 年/次。根據(jù)《防洪標(biāo)準(zhǔn)》規(guī)定,大型火力發(fā)電設(shè)施的防洪標(biāo)準(zhǔn)為100 年一遇,由于電廠的重要程度不可忽視,一般在蓄滯洪區(qū)內(nèi)是禁止進(jìn)行大型火電工程建設(shè)的。35 kV以上變電站的防洪標(biāo)準(zhǔn)不低于50 年一遇,那么在白潭湖蓄滯洪區(qū)內(nèi)建變電站必須考慮分洪的影響,一般采取墊高方式處理。220 kV 線路的防洪標(biāo)準(zhǔn)為10~20 年,那么220 kV 線路穿越白潭湖蓄滯洪區(qū)則不用額外墊高[1]。
由于分蓄洪區(qū)影響重大,其主管單位一般為省級(jí)水利行政主管部門(mén),故蓄滯洪區(qū)內(nèi)的工程建設(shè)需要進(jìn)行洪水影響評(píng)價(jià),待評(píng)審?fù)ㄟ^(guò)后才能申請(qǐng)施工許可,行政審批流程較長(zhǎng),因此在工程可研階段就要考慮到行政審批的問(wèn)題。
另外,蓄滯洪區(qū)內(nèi)架空線路需考慮凈空高的問(wèn)題,架空線路在分蓄洪狀態(tài)下也不得影響抗洪救災(zāi)安全。線路最低弧垂高度應(yīng)在分蓄洪水位以上,且要預(yù)留搶險(xiǎn)救災(zāi)船舶凈空高度4.5 m(不含電力安全超高)。
數(shù)學(xué)模型的構(gòu)建主要是選取模型范圍和網(wǎng)格剖分,計(jì)算范圍應(yīng)根據(jù)分蓄洪區(qū)工程建設(shè)所涉及區(qū)域確定,在規(guī)劃階段也要考慮主管單位所屬的行政區(qū)域,這里以白潭湖蓄滯洪區(qū)為例。將模型分為4 個(gè)獨(dú)立的區(qū)域,主要根據(jù)各區(qū)域的分洪口門(mén)銜接而定,也需考慮擬建工程的線路布置。模型的網(wǎng)絡(luò)劃分采用一維和二維銜接的河道網(wǎng)絡(luò),外加二維地形網(wǎng)格,能夠直觀地看出主要河道的概化單元。蓄滯洪區(qū)內(nèi)的地面單元應(yīng)當(dāng)以各洼淀分洪前獨(dú)立的調(diào)度方案作為銜接的控制邊界,同時(shí)考慮周邊堤防、水庫(kù)、公路等的不規(guī)則輪廓。將水庫(kù)、碼頭等小單元的網(wǎng)格高程及糙率進(jìn)行修正,使計(jì)算區(qū)域所劃分的網(wǎng)格單元可以呈現(xiàn)出任意形狀,以便后續(xù)蓄滯洪區(qū)模擬計(jì)算更為簡(jiǎn)單。另外,根據(jù)有關(guān)管理?xiàng)l例,在非防洪建設(shè)項(xiàng)目占用蓄滯洪容量,或是影響正常蓄滯洪功能時(shí),應(yīng)當(dāng)采取補(bǔ)救措施。由于還未判定是否會(huì)影響正常蓄滯洪功能,此時(shí)應(yīng)根據(jù)洪水影響評(píng)價(jià)報(bào)告在數(shù)學(xué)模型中調(diào)整擬建工程。
一般來(lái)說(shuō)在計(jì)算蓄滯洪區(qū)前需要提高洪區(qū)的計(jì)算精度,采用非網(wǎng)格結(jié)構(gòu)對(duì)項(xiàng)目區(qū)域進(jìn)行局部加密,分析白潭湖的邊界條件,從而實(shí)現(xiàn)洪水模擬,以此為分洪運(yùn)用計(jì)算提供支撐。
根據(jù)白潭湖上游支流分布及匯流情況,可設(shè)置2個(gè)入流分區(qū)模型,而內(nèi)部堤防將白潭湖分成多個(gè)區(qū)域,需要考慮過(guò)流寬頂堰流的過(guò)水流量,數(shù)據(jù)可以根據(jù)上下游水位計(jì)算結(jié)果確定。由于白潭湖內(nèi)部具有多處陸域,還應(yīng)當(dāng)考慮地面阻水程度,分析地形條件及地表阻水特征,得出糙率。按照地形、地貌條件和地面特點(diǎn),結(jié)合過(guò)往模擬經(jīng)驗(yàn)進(jìn)行計(jì)算。一般來(lái)說(shuō),行洪道的綜合糙率為0.025~0.035;房屋建筑等阻水建筑物的糙率為0.08~0.10;葦田及其他水生植物的糙率取0.08~0.12;油田、水田、草地、菜地等區(qū)域糙率取0.035~0.040。概化處理主要分析地面流通道、不泄流邊界、內(nèi)溢堤壩、自由出流邊界等數(shù)據(jù),應(yīng)當(dāng)根據(jù)蓄滯洪區(qū)的特征進(jìn)行?;?、計(jì)算。
洪水模擬分析應(yīng)當(dāng)根據(jù)上文《防洪標(biāo)準(zhǔn)》進(jìn)行,大型火力發(fā)電的防洪標(biāo)準(zhǔn)為100 年一遇,且35 kV以上變電站的防洪標(biāo)準(zhǔn)不低于50 年一遇,選取50 年和100 年一遇兩種防御標(biāo)準(zhǔn)來(lái)確定輸電線路和輸電擬建工程的防洪標(biāo)準(zhǔn),根據(jù)入流特點(diǎn)和邊界條件差異,選用二維水流數(shù)學(xué)模型,對(duì)工程建設(shè)前后應(yīng)對(duì)50 年和100 年一遇兩種標(biāo)準(zhǔn)的實(shí)際應(yīng)用情況進(jìn)行模擬分析。得到工程前后最高洪水位與最大水深等值線數(shù)據(jù),以及洪水流速矢量數(shù)據(jù)。從數(shù)據(jù)結(jié)果可知,當(dāng)白潭湖發(fā)生50 年一遇洪水時(shí),最高洪水位23.12 m,這時(shí)白潭湖周邊分洪區(qū)將全部啟用,口門(mén)位置流速加大,分洪洪水順地勢(shì)自由滿溢;當(dāng)發(fā)生100 年一遇洪水時(shí),白潭湖最高洪水位25.28 m(長(zhǎng)江中游于1954 年8 月18日出現(xiàn)洪峰水位29.73 m),此時(shí)白潭湖周邊分洪區(qū)依然全部啟用,小關(guān)口門(mén)穩(wěn)定流速。根據(jù)擬建工程在蓄滯洪區(qū)內(nèi)壅水情況,50 年一遇標(biāo)準(zhǔn)的擬建工程總壅水高度為0.02 m;100 年一遇標(biāo)準(zhǔn)總壅水高度為0.14 m??梢缘贸鰯M建的輸電工程在工程建設(shè)模擬前后的壅水幅度較小,對(duì)白潭湖滯洪區(qū)的滯洪態(tài)勢(shì)與洪水流勢(shì)的影響程度很小[2]。
3.1.1 基本方程
蓄滯洪區(qū)的數(shù)值模擬和計(jì)算通常采用平面二維水流數(shù)學(xué)模型,基本方程如式(1)至式(3)所示。
式中:Z——水位,m;h——水深,m;u、v——垂線平均流速在x,y 方向的分量,m/s;M、N——單寬流量在x,y 方向的分量,m3/(s·m),M=hu,N=hv;n——曼寧糙率系數(shù);g——重力加速度。
3.1.2 計(jì)算方法
模型通常采用的數(shù)值方法為有限體積法。在有限體積法計(jì)算中,將基本方程對(duì)每一控制體積分,得到一組以控制體特征量平均的物理量為未知數(shù)的代數(shù)方程組,同時(shí)沿坐標(biāo)方向?qū)Ψ匠探M進(jìn)行離散,形成的離散方程與有限差分法相似。由于采用了守恒型的微分方程并對(duì)每一計(jì)算單元進(jìn)行質(zhì)量和動(dòng)量守恒形式的離散,使得微分方程包含的守恒性質(zhì)在每一個(gè)控制容積上都得到滿足,并保持各單元界面兩側(cè)相鄰控制體的計(jì)算輸運(yùn)通量相等,整個(gè)計(jì)算域上保持守恒性?;痉匠屉x散后的形式如式(4)所示。
y 方向水流運(yùn)動(dòng)方程可用與式(4)相同的方法離散。為節(jié)省篇幅,不再列出其離散結(jié)果。如式(5)至式(7)所示。
計(jì)算所需的初始條件包括:蓄滯洪區(qū)地形、分洪區(qū)內(nèi)的初始水深和流速;邊界條件包括:分洪口門(mén)處外江水位和流量。
當(dāng)分蓄洪民垸分洪運(yùn)用后,洪水對(duì)建設(shè)項(xiàng)目將產(chǎn)生一定的影響,主要的影響為在分洪運(yùn)用過(guò)程中對(duì)擬建工程基礎(chǔ)產(chǎn)生的沖刷。蓄滯洪區(qū)局部沖刷采用黏性土沖刷的計(jì)算公式。鐵道部黏土橋渡沖刷研究小組的《黏土橋渡沖刷天然資料分析報(bào)告》中,推薦以下兩個(gè)公式。
式中:hb——橋墩局部沖刷深度,m;B——橋墩計(jì)算寬度,m;Kξ——墩形系數(shù);IL——沖刷坑范圍內(nèi)黏性土液性指數(shù),該式中IL的適用范圍為0.16~1.19;e——沖刷坑范圍內(nèi)黏性土空隙比,式中e 范圍為0.6~1.5;V——橋墩前行進(jìn)垂線平均流速,m/s。當(dāng)V 小于土壤的不沖刷流速時(shí),不會(huì)發(fā)生沖刷。
由公式可以看出,影響局部沖刷的主要因素有3 個(gè):基礎(chǔ)結(jié)構(gòu)(墩形系數(shù) Kξ和基礎(chǔ)計(jì)算寬度 B)、基礎(chǔ)土質(zhì)(液性指數(shù)IL和孔隙比e)和水流條件(流速V)。
為研究擬建工程的建設(shè)對(duì)周邊區(qū)域的影響程度,確保局部計(jì)算更準(zhǔn)確,需要計(jì)算水位壅高在50 年一遇和100 年一遇時(shí)的差距,并分析兩者的流速差值。雖然水位壅高顯示擬建工程對(duì)蓄滯洪區(qū)的影響較小,但可能會(huì)對(duì)洪水流速產(chǎn)生較大的影響。通過(guò)模擬整體計(jì)算得到局部沖刷影響數(shù)據(jù),能夠反映擬建工程對(duì)整個(gè)蓄滯洪區(qū)的影響,但卻無(wú)法反映擬建工程的局部流場(chǎng)。因此,要針對(duì)擬建工程附近流場(chǎng)進(jìn)行局部加密計(jì)算,范圍以30 m×30 m 的正方形為主,底高程設(shè)為平面,在模型中概化水深、流速情況以及擬建工程的形式。本文以平均水深和平均流速計(jì)算擬建工程流場(chǎng)的變化情況。
根據(jù)實(shí)驗(yàn)結(jié)果可以得知,流場(chǎng)變化對(duì)附近洪水流場(chǎng)的影響較為明顯,水流條件與基礎(chǔ)土質(zhì)的差異會(huì)對(duì)擬建工程周?chē)鲌?chǎng)產(chǎn)生較大影響,最大流速差會(huì)產(chǎn)生在工程項(xiàng)目垂直于流速的方向上,也就是迎水面方向,最大可使流速減小56%。而且不同的基礎(chǔ)結(jié)構(gòu)也會(huì)產(chǎn)生不同影響,如圓柱群對(duì)流速影響要小于方柱群。
(1)在蓄滯洪區(qū)內(nèi)開(kāi)展電力工程建設(shè)需要明確該蓄滯洪區(qū)的分洪運(yùn)用概率,并為洪水影響評(píng)價(jià)報(bào)告的編制與行政審批預(yù)留時(shí)間。
(2)穿越蓄滯洪區(qū)的架空線路最低弧垂高度應(yīng)在分蓄洪水位以上,且要預(yù)留搶險(xiǎn)救災(zāi)船舶凈空高度4.5 m(不含電力安全超高)。
(3)在分洪運(yùn)用計(jì)算時(shí),本文僅給出了一種常見(jiàn)的水流模型與沖刷計(jì)算方式,可以根據(jù)實(shí)際情況進(jìn)行調(diào)整和優(yōu)化。
(4)蓄滯洪區(qū)內(nèi)建設(shè)電力工程還需要著重預(yù)留分洪口門(mén)位置,在分洪口門(mén)處的流速、沖刷影響較大。工程建設(shè)應(yīng)盡量遠(yuǎn)離分洪口門(mén),根據(jù)以往工程經(jīng)驗(yàn),距離保持在500 m 以上。電力工程前、后對(duì)滯洪區(qū)影響基本一致,擬建電力工程對(duì)蓄滯洪區(qū)的影響可忽略不計(jì)。
(5)電路工程處于滯洪區(qū)內(nèi),在汛期發(fā)生時(shí)應(yīng)實(shí)時(shí)保持與水利部門(mén)聯(lián)系,做好應(yīng)對(duì)洪水突發(fā)時(shí)的預(yù)案和對(duì)策,在發(fā)生水利事件時(shí)能夠迅速協(xié)商解決。