曹 侃,吳德安,丁 堅,趙明志
(河海大學(xué)海岸災(zāi)害及防護教育部重點實驗室,南京210098)
蕪申運河處于長三角高等級航道網(wǎng)的建設(shè)范圍內(nèi),具有極其重要的地位,是長三角地區(qū)整體發(fā)展的組成部分,安徽省主要骨干航道架之一。它的建成對發(fā)展皖南經(jīng)濟及沿河開發(fā)有重要作用,同時大大縮短蕪申航線里程,提高了航行的安全性,有利內(nèi)陸防洪排澇。主要支流青弋江、水陽江為山區(qū)性河流,來水暴漲暴落,洪水含沙量較大,航道不穩(wěn)定。同時受長江來水影響,水沙條件復(fù)雜。三汊河河段(圖1)工程前航道曲率半徑過小,航道距離護堤過近,受青弋江來水影響,航道流態(tài)不穩(wěn)、橫流流速過大,對行船安全形成威脅,影響航道的正常使用;水流交匯處沖刷形成深坑,影響岸堤穩(wěn)定,威脅附近居民生命財產(chǎn)安全。
為了船舶的安全航行,改善該段的航道現(xiàn)狀,穩(wěn)定護堤,較為可行的是進行工程方案整治。根據(jù)安徽省交通規(guī)劃設(shè)計研究院提供的初步航道設(shè)計方案,建立物理模型。通過試驗發(fā)現(xiàn)雖然解決了航道水深不足和彎道曲率半徑過小的問題,但在中大水流情況下橫流對航行安全的影響仍然存在。在以往的航道整治工作中,數(shù)值模擬技術(shù)得到了廣泛的應(yīng)用[1-3],為解決上述問題,運用Delft 3D 水動力模塊建立數(shù)值模型,根據(jù)實測水文資料,找出最不利的水流組合[4]進行模擬計算,調(diào)整優(yōu)化初步設(shè)計工程方案,改善流態(tài),提出整治措施及推薦方案,對數(shù)學(xué)模型優(yōu)化工程方案在物理模型上進行驗證和進一步優(yōu)化[5-9]。
流場的模擬采用了荷蘭Delft 水利機構(gòu)的Delft 3D 模型,主要使用Delft 3D-flow 模塊。Delft-3D-flow 模型是一個三維的水動力輸運模型,該模型采用曲線正交網(wǎng)格。水動力模塊建立在Navier-Stokes 方程的基礎(chǔ)上,采用交替方向法(ADI)對該坐標(biāo)系下的控制方程組進行離散求解。應(yīng)用沿水深平均的兩維水流數(shù)學(xué)模型來計算,采用大小模型嵌套的技術(shù)進行數(shù)值模擬。其特點是由多個模塊靈活組成,模塊間聯(lián)系性強,采用貼體正交曲線網(wǎng)格,計算穩(wěn)定,精度高。
(1)正交曲線坐標(biāo)系下的水深平均的連續(xù)方程。
(2)動量方程。
水平ξ 方向上
水平η 方向上
式中:u,v,ω 分別表示在正交曲線坐標(biāo)系下ξ,η,σ3 個方向上的速度分量,其中ω 是定義在運動的σ 空間的豎向速度,在σ 坐標(biāo)系統(tǒng)中由以下的連續(xù)方程求得
式中:Gξξ,Gηη為曲線坐標(biāo)系轉(zhuǎn)換為直角坐標(biāo)系的轉(zhuǎn)換系數(shù);Fξ,F(xiàn)η分別為ξ 和η 方向的紊動棟梁通量;Mξ,Mη分別表示ξ 和η 兩方向上動量的源或匯,包括建筑物引起的外力、波浪切應(yīng)力,以及引排水產(chǎn)生的外力;ρ0為水體密度;vv為豎向渦動系數(shù);f 是柯氏力系數(shù),取決于地理緯度和地球自轉(zhuǎn)的角速度Ω,f 可以用下式表示:f=2Ωsinφ,φ 為北緯緯度;Pξ,Pη分別表示ξ 和η 兩方向上的水壓力梯度
式中:Patm包括浮體建筑物引起的壓力在內(nèi)的自由面壓力。
根據(jù)實際控制區(qū)域分別順應(yīng)河流兩個方向生成橫向167 個,縱向92 個的正交矩形貼體網(wǎng)格,有效網(wǎng)格控制在10~15 m,滿足計算精度要求。
為充分論證蕪申運河青弋江段航道整治工程的建設(shè)規(guī)模,擬定合適的工程規(guī)模布置方案,需在清水河與青弋江交匯的三河口處開展水文測驗,以便為航道整治物理模型的驗證以及工程布置階段的試驗提供技術(shù)支撐。各斷面具體位置如表1 所示。
運用軟件Flow 模塊計算。根據(jù)地形圖中提取的水深樣本文件生成網(wǎng)格節(jié)點上的水深文件。它是利用在網(wǎng)格節(jié)點附近的樣本點通過三角形插值在網(wǎng)格節(jié)點上生成水深點。流場模擬時間12 h,根據(jù)生成網(wǎng)格密度和計算的穩(wěn)定性,時間步長不宜取得過大,本計算過程中,時間步長選為0.1 min。底部糙率選用曼寧系數(shù)0.02,邊界為無滑移條件。每隔30 min 輸出一次計算結(jié)果。
分別計算了用流量和水位控制的兩種情況,結(jié)果表明2 種情況流場、流速基本一致,接近實際測驗情況。為了更加準(zhǔn)確的對比模擬工程前后流場變化,3 個開口邊界分別通過流量和對應(yīng)水位控制,斷面1 通過流量(1 210 m3/s),斷面2、3 通過對應(yīng)流量的水位(9.7 m)控制。
表1 斷面位置Tab.1 Cross-sectional position
根據(jù)驗證結(jié)果,流態(tài)上數(shù)值模型計算結(jié)果和物理模型基本吻合;流速數(shù)值模擬計算結(jié)果對比實測值略小。航道流速最大值平均誤差在10%以內(nèi),符合規(guī)范要求,可應(yīng)用于本工程計算。
根據(jù)設(shè)計部門提供數(shù)據(jù),航道底高程-1.5 m,邊坡1:5,做出如圖9 所示的設(shè)計航道地形。在此基礎(chǔ)上計算流場、流速。
分析計算結(jié)果表明:圖9 中①處橫向流速過大(最大1.2 m/s)依然是影響通航安全的主要因素,為了使船舶的航行穩(wěn)定安全,根據(jù)設(shè)計要求航道橫向流速要控制在0.3 m/s 左右。主要問題存在于青弋江來流過于集中且相對設(shè)計航道接近垂直,為了解決此問題,提供典型四種計算結(jié)果供研究比對。
通過用數(shù)值模型和物理模型的比證分析,航道的大部分滿足通航要求。為了取證方案改善效果,截取了計算網(wǎng)格中的關(guān)鍵部分(圖20)流速進行比對。原設(shè)計航道在圖9 中①處的流速達到了1.2 m/s,斷面平均流速0.69 m/s。
方案一在計算中挖掉大部分礙流洲灘,以此來驗證其對流速流向的影響。實驗證明,對于不同水流條件下,水流流態(tài)和主流流向都作相應(yīng)調(diào)整,但從流路上看,在高水位大流量情況下,口門交匯區(qū)西北側(cè)橫向流速偏大(0.8 m/s),高水位大流速的情況下,水流流向發(fā)生較大改變。
方案二主要思想是退建右岸護堤,研究右岸大堤對水流的影響。經(jīng)計算,退堤后水流導(dǎo)向明顯,但是水流過于集中,在設(shè)計航道處不能完全滿足通航要求,而且對原河道分流比影響較大,青弋江下游流量明顯減少,影響河道穩(wěn)定。
方案三是在方案二的基礎(chǔ)之上加以修改,因為圖9 中②處凸起岸堤對水流影響較大,為了達到最佳效果,對其進行退建。經(jīng)計算,此方案效果較好,在設(shè)計航道處的流速已經(jīng)降到0.4 m/s,平均流速0.24 m/s,在增大了有效泄洪和通航面積面積的同時對河道的分流影響較小,此方案有效的解決了航道隱患,最大程度保證了通航船舶安全。但此方案工程量較大且涉及到堤防的重建工作,可考慮長期規(guī)劃。
方案四從經(jīng)濟和實用性出發(fā),口門擴挖,分散水流使流速減小,在大流量時,設(shè)計航道處流速較大,原航道適當(dāng)加寬,水流條件良好,可以使用原航道通航;枯水季節(jié)小流量時設(shè)計航道滿足通航要求,根據(jù)季節(jié)流量的大小變化,適當(dāng)調(diào)整通航路線。
目前航道整治工程的實施多為以一個水道或兩個水道作為一個工程單元,而水道一般是以歷史上的2個自然節(jié)點間的河段定義的。節(jié)點在沖積性河流河床演變中扮演著重要的角色,起著十分重要的作用。上游河段的整治一般只會出現(xiàn)兩種情況,或強化節(jié)點的挑流作用或弱化節(jié)點的挑流作用,基本不改變節(jié)點的挑、導(dǎo)流作用難以做到。節(jié)點作用的改變,作為下一水道的輸入水流條件將隨之發(fā)生改變,從而影響下游水道的河床演變。
以一個水道作為一個航道整治工程的實施對象應(yīng)該不存在大的問題,而我們研究的對象是一個節(jié)點上的三段水道,且河岸的形態(tài)為彎曲半徑適度的弧線形,青弋江水道較短且曲折,挑流作用往往占優(yōu),限制河流向一側(cè)或兩側(cè)發(fā)展而且這種挑流作用在年內(nèi)的不同水期和年際間都在不斷地變化,也由于彎曲半徑的不盡合理,水流動力軸線不斷地遷徙,表現(xiàn)出不穩(wěn)定的一面。故此,解決這個問題的辦法就是盡可能以一個較長的河段作為航道整治的單元研究對象,而物理模型試驗研究時,河段的長度是有限的,我們認(rèn)為更應(yīng)該關(guān)注的是洲灘的合理布局。
由上,以2 個主導(dǎo)河岸型節(jié)點之間的河段作為工程河段我們認(rèn)為是合適的。我們在研究它的時候要預(yù)判一個灘段對它將來是否有利。在預(yù)測的基礎(chǔ)上,適時采取工程措施,對造床流量下的河勢進行控制,以穩(wěn)定或塑造航道整治的目標(biāo)河型。
[1]李琳琳,余錫平.裁彎對上下游流態(tài)影響的三維數(shù)學(xué)模型研究[J].水力發(fā)電學(xué)報,2010,29(5):183-189.LI L L,YU X P. Study on effects of river cutoffs by a 3-D numerical model[J]. Journal of Hydroelectric Engineering,2010,29(5):183-189.
[2]鄭金海,王義剛.黃浦江航道改善方案的三維水流數(shù)學(xué)模型計算[J].水運工程,336(1):34-37.ZHENG J H,WANG Y G. Application of Three-Dimensional Flow Numerical Model to Channel Improvement Project in the Huangpu River[J]. Port & Waterway Engineering,336(1):34-37.
[3]余新明,談廣鳴,趙連軍,等.天然分汊河道平面二維水流泥沙數(shù)值模擬研究[J].四川大學(xué)學(xué)報:工程科學(xué)版,2007,39(1):33-37.YU X M,TAN G M,ZHAO L J,et al.Planar 2-D Flow and Sediment Numerical Modeling of Branching River[J].Journal of Sichuan University:Engineering Science Edition,2007,39(1):33-37.
[4]姚仕明,張超,王龍,等.分汊河道水流運動特性研究[J].水利發(fā)電學(xué)報,2006,25(3):49-57.YAO S M,ZHANG C,WANG L,et al.Study on the characteristics of flow movement in branching river[J]. Journal of Hydroelectric Engineering,2006,25(3):49-57.
[5]張蔚,嚴(yán)以新,諸裕良,等.人工采沙及航道整治對珠江三角洲水流動力條件的影響[J].水利學(xué)報,2008,39(9):1 098-1 104.ZHANG W,YAN Y X,ZHU Y L,et al. Impact of sand excavation and waterway regulation on hydrodynamics of Pearl River networks[J]. Journal of Hydraulic Engineering,2008,39(9):1 098-1 104.
[6]柴曉玲,余啟輝,要威,等.復(fù)雜河網(wǎng)地區(qū)航道整治工程對防洪的影響分析[J].人民長江,2011(10):1-6.CHAI X L,YU Q H,YAO W, et al. Research on influence of navigation improvement project on complex river network area[J].Yangtze River,2011(10):1-6.
[7]楊昌道.沙潁河航道安徽段航道整治工程航道尺度研究[J].水運工程,2011(7):131-135.YANG C D.On channel dimensions of Shayinghe Anhui section channel regulation engineering[J]. Port& Waterway Engineering,2011(7):131-135.
[8]蔡大富.長江中游牯牛沙水道航道整治一期工程工程可行性研究報告[R].武漢:長江航道規(guī)劃設(shè)計研究院,2008.
[9]李旺生.長江中下游航道整治技術(shù)問題的幾點思考[J].水道港口,2007,28(6):419-423.LI W S. Some thoughts about technical problems on regulation of waterway in the middle and lower reach of the Yangtze River[J]. Journal of Waterway and Harbor, 2007,28(6):419-423.