李文明,鐘 瑾,朱 莎
(1.江蘇省淮沭新河管理處,江蘇 淮安 223001,2.淮陰工學(xué)院,江蘇 淮安 223001)
船閘是水利樞紐工程中重要通航設(shè)施[1],能保障水上交通安全穩(wěn)定運行。 研究水閘樞紐工程通航水流狀態(tài),可為船閘調(diào)度、維護(hù)加固等提供理論依據(jù)。 麥建清等[2]、徐剛等[3]為研究船閘結(jié)構(gòu)水力特征,采用滲流場模擬計算方法,對船閘樞紐水力特性開展計算分析,探討了船閘內(nèi)流速、泥沙淤積以及水流流態(tài)、流線等特征,對推動船閘工程的高標(biāo)準(zhǔn)運營提供了數(shù)據(jù)支撐。 韓瑩瑩等[4]、沈立群等[5]為研究船閘樞紐的最優(yōu)化調(diào)度運營方案,設(shè)計開展水工模型試驗,根據(jù)模型試驗結(jié)果,評價了閘體流速、水位以及壓強等水力參數(shù),對提高水工建筑物建設(shè)水平具有重要意義。 導(dǎo)流堤是船閘樞紐中重要構(gòu)件,楊峰等[6]為優(yōu)化導(dǎo)流堤平面布置,采用對比分析方法,對導(dǎo)流堤結(jié)構(gòu)體型進(jìn)行優(yōu)化,并結(jié)合滲流場特征、結(jié)構(gòu)靜力學(xué)特征等,分析了導(dǎo)流堤設(shè)計對船閘運營的影響。
為了研究沭新北船閘工程的流量調(diào)度、導(dǎo)流堤體型特征,本文建立船閘水工模型,分析不同流量調(diào)度方案、不同導(dǎo)流堤體型下船閘航道水力特性,并提出閘體運營參考建議。
淮沭新河是分淮入沂、淮水北調(diào)的重要人工河道,其聯(lián)通了連云港、淮安等地區(qū),對地區(qū)用水緊張、農(nóng)業(yè)灌溉效率較低等具有重要緩解作用。
薔薇樞紐又稱沭新河渠尾控制樞紐,通過人為影響河道水位,使淮沭新河、沭新干渠等能夠保證輸水勢能,形成一定的水位落差效應(yīng),方便下游灌區(qū)輸水灌溉、通航以及水力發(fā)電,薔薇河地涵區(qū)域分布見圖1。 該地涵工程共有3 孔,單孔凈寬4.5m,設(shè)計過閘流量65m3/s,規(guī)劃年輸水超過2 000×104m3,惠及下游農(nóng)田超過3.333×104hm2。
圖1 薔薇河地涵區(qū)域分布
沭新北船閘是薔薇河地涵樞紐工程重要組成部分,位于薔薇河地涵工程下游1 000m,具有800t 級通航作用,同時承擔(dān)著區(qū)域排澇防洪。 設(shè)計最大排澇流量30m3/s,設(shè)計為六級航道,是淮沭新河引調(diào)水工程的重要交通線。
根據(jù)對沭新北船閘工程調(diào)研得知,其上下游最大水位分別為7.8、7.5m,通航運行水位分別為6、3.5m,而通航水位最大可分別調(diào)整為7.5、6m,設(shè)計通航標(biāo)準(zhǔn)為800t 船隊。 從沭新北船閘工程結(jié)構(gòu)分析,全閘室長度135m,凈寬度12.5m,閘首凈寬10m,閘室體型設(shè)計采用漿砌石承重結(jié)構(gòu),底板厚0.6m,所在高程-0.5m,閘頂高程7.5m,全閘身漿砌石砌塊在多次運營監(jiān)測及除險加固中仿真分析,均表明閘體自身結(jié)構(gòu)安全可靠度好,最大沉降位移不超過5mm。 船閘啟閉動力裝置來源于4 臺推桿啟閉機,最大轉(zhuǎn)速880 轉(zhuǎn)/min,上下游閘門頂高程均為8m,底高程分別為3、-0.5m,閘門為平面鋼閘門,寬度6.06m,可承受最大時均壓強超過50kPa,能夠滿足峰值水位下?lián)跛蜷]水。
基于沭新北船閘工程運營監(jiān)測分析,通航水流條件常遇非穩(wěn)定水流、流線集中以及局部航段泥沙淤積嚴(yán)重的現(xiàn)象,這與沭新北船閘樞紐工程結(jié)構(gòu)體型布置有關(guān)。 如上下游航道進(jìn)流調(diào)度方式以及導(dǎo)流堤體型結(jié)構(gòu)等,均會改變沭新北船閘航道水流狀態(tài),從而影響通航條件。 為此,沭新北船閘管理部門考慮在不干擾通航條件現(xiàn)狀的基礎(chǔ)上,先期對影響航道水流特征的顯著條件開展水工模型試驗研究,為樞紐工程優(yōu)化平面布置提供依據(jù)。
在考慮沭新北船閘工程實際情況的前提下,按照水工模型制作標(biāo)準(zhǔn),將船閘樞紐工程包含的節(jié)制閘、船閘、上下游航道以及導(dǎo)流堤等組成部分,均制作出相應(yīng)室內(nèi)模型。 該模型研究影響范圍包括上游航道4km、下游3. 5km,河段總長9km,河面寬度1.5km,幾何比尺1/80,研究試驗段主軸長125m,進(jìn)出水口長5m,模型總面積超過1 500m2。
為貼合船閘樞紐工程實際模型,采用水泥砂漿作為河道硬化層,表層凸起厚度不超過5mm,局部泥沙淤積地采用預(yù)制面板加密處理。 模型中,不論是閘身或是航道,其混凝±剛性材料糙率均需滿足0.007 5。 圖2 為所建船閘樞紐工程模型示意,供、退水管路均采用管徑20cm 的剛性管道輸水,不考慮其輸水耗散,且認(rèn)為管壁為不透水邊界。
圖2 船閘樞紐工程模型
根據(jù)水工模型試驗原理,本模型流量比尺為48 600,流速與時間比尺均為8.5,且模型水流為紊流特征[7],雷諾數(shù)、水深等水力參數(shù)分別不低于1 000、2cm。 同時,鑒于船閘航道泥沙淤積影響,故需考慮渾水河床受沖刷影響,其試驗?zāi)嗌沉奖3植怀^2mm,且為天然砂、塑料砂以及滑石粉等材料的混合體,容重2.1g/cm3,沉降相似以及起動相似比尺分別為2.35、3.2,顆粒細(xì)沙的中值粒徑0.05mm,顆粒沙的各項物理參數(shù)均與淮沭新河道泥沙沉降屬性一致。 模型的進(jìn)出水口均設(shè)置有擋水閘門,為剛性平面結(jié)構(gòu)體型,寬、高分別為2.6、3.2m,圖3 為該船閘實際模型的細(xì)部構(gòu)造特征。
圖3 船閘水工模型細(xì)部構(gòu)造
為全面準(zhǔn)確獲取船閘運營期內(nèi)航道水流狀態(tài),采用水位測試儀、流速標(biāo)記儀等測量設(shè)備開展數(shù)據(jù)監(jiān)測,見圖4。
各測量設(shè)備均布置在水流上、中、下3 個層次深度內(nèi),設(shè)備傳感器誤差不超過0.1%。 從船閘及航道運營監(jiān)測考慮,所有數(shù)據(jù)分析需覆蓋全樞紐工程,故在航道、閘室內(nèi)布置有各監(jiān)測斷面點,圖5(a)為閘室、下游航道上分布的監(jiān)測斷面點,閘室內(nèi)間距0.5m,航道上間距1m,共有A1-A10、B1-B15 斷面監(jiān)測點。
在開始水工模型試驗前,對航道相關(guān)數(shù)據(jù)先期開展了在模型監(jiān)測數(shù)據(jù)與實際河段數(shù)據(jù)對比。圖5(b)為上游航道流量320m3/s 時,下游河段4~8m 內(nèi)流速數(shù)據(jù)對比,最大誤差0.018m/s,流速分布以及趨勢特征具有一致性,滿足試驗相似性要求
本模型試驗中,需要探討航道水流影響因素,設(shè)計試驗組包括流量調(diào)度組與導(dǎo)流堤結(jié)構(gòu)體型優(yōu)化組。 前一組參數(shù)為不同方式下的流量組合,設(shè)計通航流量為300、600、900、1 200 和1 500m3/s 共5 組,均為敞泄流量組合方式。 后一組影響因素為導(dǎo)流堤長度參數(shù),并只考慮導(dǎo)流堤身向上游延伸長度,分別設(shè)計延伸長度為25、40、55、70、85 和100m 共6 種方案。 其中,延伸長度100m 后,導(dǎo)流堤已超過閘身與河道彎角處。 基于兩組因素下的模型試驗探討,分析航道、閘室內(nèi)水流水力參數(shù),為工程建設(shè)提供數(shù)據(jù)支撐。
基于不同流量調(diào)度方式下水力模型試驗,獲得各方案下航道回流特征,見圖6。
圖6 航道回流特征
分析圖6(a)可知,不論調(diào)度流量工況為何值,只要保證航道口門區(qū)仍具有相應(yīng)穩(wěn)定水流,則回流區(qū)所在航段仍位于口門區(qū)與導(dǎo)流堤左側(cè),該區(qū)域為航道回流出現(xiàn)的集中區(qū)。 從回流區(qū)分布尺寸來看,調(diào)度流量愈大,則回流區(qū)分布面積愈廣,回流區(qū)的長、寬參數(shù)均增大。 在圖6(a)中,調(diào)度流量為300m3/s 時,回流區(qū)分布的最大長、寬分別為82、24.5m;而調(diào)度流量為600、1 200m3/s時,相應(yīng)的回流區(qū)最大長度分別達(dá)85、95m;而調(diào)度流量超過1 200m3/s 以后,為1 500m3/s 工況時,回流區(qū)長、寬較之調(diào)度流量300m3/s 下分別增大52.4%、140.8%。 調(diào)度流量增加,河道水流流速提高,水流黏滯性影響河道主流發(fā)展,從而提高了回流區(qū)分布范圍,回流強度、回流漩渦均會顯著提升,特別是在調(diào)度流量超過1 200m3/s后,影響更敏感[8-9]。
從圖6(b)航道回流區(qū)最大回流速度對比可知,流速與調(diào)度流量為正相關(guān)關(guān)系。 在調(diào)度流量為300~900 m3/s 下,回流速度的增長仍較緩慢,流速分布為0.1~0.24m/s;而在1 200m3/s 工況后,流速陡增,調(diào)度流量1 200~1 500 m3/s 兩者間增幅可達(dá)77.9%,調(diào)度流量梯次變化300m3/s,則回流速度增幅可達(dá)67.3%。
綜合對比分析認(rèn)為,調(diào)度流量會改變河道主流特征,從而影響航道回流區(qū)分布擾動,引起回流、斜流等非穩(wěn)定水流,不利于航道水流紊動,限制了船閘通航水流條件。
在考慮導(dǎo)流堤延伸長度的前提下,開展導(dǎo)流堤結(jié)構(gòu)體型優(yōu)化對比試驗,圖7 為航道B1-B15各斷面上流速分布變化特征。 從圖7 流速對比可知,在導(dǎo)流堤不同延伸長度下,航道斷面流速分布變化具有明顯差異性,特別是在峰值流速斷面等航段處。 當(dāng)延伸長度為25~55m 時,航道各斷面流速整體呈遞增-遞減-穩(wěn)定變化,峰值流速位于斷面B5;而延伸長度為70~100m 時,相應(yīng)斷面流速為遞增-小幅遞減-再次遞增變幅,二次遞增段分布于B11-B16,一次遞減段降幅分布為B7-B11。 對比可知,當(dāng)導(dǎo)流堤延伸長度過長,不利于控制航道下游末端流速穩(wěn)定,促發(fā)水流裹挾過多水利勢能,從而造成航道水流紊動[10]。
圖7 航道B1-B15 各斷面上流速分布變化
分析流速水平可知,導(dǎo)流堤延伸長度愈大,航道斷面流速水平愈低。 在導(dǎo)流堤延伸長度25m方案內(nèi),其峰值流速為1.56m/s,航道全斷面平均流速0.83m/s;而延伸長度為55、85、100m 時,相應(yīng)的峰值流速較之前者分別減少34.5%、85.5%,而斷面平均流速分別為0.55、0.12m/s。 從斷面流速的降幅對比來看,在延伸長度70~100m 方案后,其流速降幅實質(zhì)上較小,方案間各斷面測點處流速較為接近,B1-B15 各測點斷面流速最大變幅為12.5%,屬于延伸長度85~100m 方案間B7 測點。 綜合對比可知,延伸導(dǎo)流堤長度,可減小非穩(wěn)定滲流,約束水流活動,確保航道水流平穩(wěn)。
航道內(nèi)水流狀態(tài)不僅需要考慮流速特征,還需考慮航道泥沙淤積變化,確保航道水流安全,圖8 為航道內(nèi)各測點處泥沙淤積厚度分布變化特征。 由圖8 可看出,航道斷面上泥沙淤積變化具有一致性,淤積峰值厚度位于斷面測點B12。在閘室-航道區(qū)段內(nèi),泥沙淤積仍處于低水平;而在航道斷面B7 后,泥沙淤積較顯著,部分測點淤積厚度超過30cm。 總體上,泥沙淤積曲線呈穩(wěn)定-遞增-遞減變化。 分析認(rèn)為,改變導(dǎo)流堤延伸長度,不會影響航段泥沙淤積分布趨勢,只會改變不同斷面處泥沙淤積厚度,而泥沙淤積分布特征主要受航段流量影響。
圖8 航道B1-B15 各斷面上泥沙淤積分布變化
分析泥沙淤積與導(dǎo)流堤延伸長度關(guān)系可知,兩者為負(fù)相關(guān)變化。 在延伸長度25m 時,航道全斷面上泥沙淤積厚度分布7.5~35.5cm;而延伸長度40、70m 時,峰值厚度較之前者分別減少28.1%、56.7%;總體隨導(dǎo)流堤延伸長度每梯次變幅15m,其淤積厚度峰值平均可下降4.2cm,降幅為15.9%。 此外,從泥沙淤積厚度變幅也可看出,導(dǎo)流堤延伸長度在70m 以內(nèi)時,泥沙淤積厚度受削弱明顯高于延伸長度70~100m 內(nèi),即導(dǎo)流堤延伸長度為25~70m 內(nèi)時,導(dǎo)流堤長度能夠適配于泥沙淤積狀態(tài)[11],有利于航段水沙懸浮運動。
1)回流區(qū)分布于航道口門區(qū)與導(dǎo)流堤左側(cè),調(diào)度流量改變,不影響回流區(qū)分布集中區(qū)。 調(diào)度流量愈高,回流區(qū)分布長、寬愈大,回流強度、回流漩渦等均提高。 回流速度與調(diào)度流量呈正相關(guān),調(diào)度流量每梯次變化300m3/s,則最大回流速度可增長67.3%。
2)導(dǎo)流堤延伸長度不同,則航道斷面流速分布變化具有顯著差異。 在延伸長度25~55m 與70~100m 時,流速分別呈遞增-遞減-穩(wěn)定、遞增-小幅遞減-再次遞增變化。 導(dǎo)流堤延伸長度與航道斷面流速為負(fù)相關(guān),且延伸長度70~100m 后,流速降幅較小。
3)導(dǎo)流堤延伸長度不同,但航道斷面泥沙淤積分布變化均為穩(wěn)定-遞增-遞減。 泥沙淤積厚度與導(dǎo)流堤延伸長度為負(fù)相關(guān)變化,隨延伸長度每梯次變幅15m,其淤積厚度峰值平均下降15.9%,尤以延伸長度25~70m 內(nèi)限制泥沙沉降更顯著。