何志敏， 李九福， 劉惠康， 張玉成， 余小強

(1. 南昌市政公用投資控股有限責(zé)任公司， 江西 南昌　330000；2. 廣州地鐵設(shè)計研究院有限公司， 廣東 廣州　510000；3. 廣東省水利水電科學(xué)研究院， 廣東 廣州　510610)

?

沉管隧道管段浮運中的水流阻力性能及其在回旋區(qū)轉(zhuǎn)體的流速流場模擬分析

何志敏1， 李九福1， 劉惠康2,*， 張玉成3， 余小強1

(1. 南昌市政公用投資控股有限責(zé)任公司， 江西 南昌330000；2. 廣州地鐵設(shè)計研究院有限公司， 廣東 廣州510000；3. 廣東省水利水電科學(xué)研究院， 廣東 廣州510610)

南昌紅谷隧道沉管浮運過程存在2個關(guān)鍵的風(fēng)險控制點，分別為浮運過南昌大橋及回旋區(qū)轉(zhuǎn)體，在這2處水流流向及流速復(fù)雜，浮運施工風(fēng)險非常高。為校核浮運方案的合理性，需對管段浮運過程中水流阻力性能進行分析，以保證浮運安全?；贔luent和MIKE流體計算軟件，通過數(shù)值模擬的方法對沉管管段在2個關(guān)鍵風(fēng)險控制點中所受水流阻力進行分析。對于管段浮運過南昌大橋過程，得到了管段所受水流阻力大小及其變化情況； 對于回旋區(qū)轉(zhuǎn)體過程，先對水阻力系數(shù)Cw進行率定使其適用于本工程，再通過數(shù)值模擬得到不同水文條件下的回旋區(qū)流場，兩者結(jié)合得到管段所受水流阻力大小。以期為復(fù)雜邊界條件下管段水流阻力計算提供一套方法，計算結(jié)果為隧道管段浮運方案的制定提供參考。

紅谷沉管隧道； 管段浮運； 水流阻力； 數(shù)值模擬

0　引言

沉管法施工隧道與礦山法和盾構(gòu)法相比具有眾多優(yōu)點，在國內(nèi)外水底隧道工程中得到了廣泛應(yīng)用[1-3]。在沉管隧道發(fā)展的百年歷史中，有很多關(guān)鍵技術(shù)問題得到了解決并不斷地改進，如管段預(yù)制[3]、混凝土裂縫控制、基礎(chǔ)處理以及管段接頭處理[4]等。但是沉管隧道技術(shù)仍有一些需要繼續(xù)改進的技術(shù)問題，如在浮運過程中管段的受力及穩(wěn)定性等問題。

由于在施工過程中沉管隧道不僅受到河流、氣象以及航道等自然因素的影響，還受到管段自身的浮重力、水流阻力、浮運設(shè)備及其布置方式等人為因素的影響[5-6]，因此，成功的浮運施工離不開天時地利人和。特別在浮運路線復(fù)雜且周邊控制性建筑物多的情況下，為保證管段的成功浮運，需對隧道沉管管段浮運過程中的阻力性能及穩(wěn)定性進行研究。

目前對沉管管段浮運施工的研究主要采取數(shù)值模擬和試驗2種方法。ZHANG Dexin等[7]對沉管管段在水面浮運及沉放過程進行了數(shù)值模擬； 蔣龍皎[8]通過數(shù)值模擬的方法，模擬隧道沉管在淺水航道中浮運時的阻力情況； 周敏[9]利用模型試驗的方法，進行縮尺比例的模型實驗，然后通過相似理論把模型實驗結(jié)果換算到實際沉管，得到實際沉管的受力規(guī)律； 羅甜[10]通過模型試驗，獲得管段靜水拖曳的阻力系數(shù)以及風(fēng)浪流復(fù)雜因素影響下的浮運拖航纜力。以往的學(xué)者對管段浮運施工的研究大多集中在簡單水文條件情況下的分析，而對于周邊復(fù)雜環(huán)境引起的復(fù)雜水文條件下管段浮運施工研究較少，因此有必要針對項目的特點進行針對性的研究。

本研究依托南昌紅谷沉管隧道浮運施工項目，基于RNGκ-ε紊流模型[11-12]，采用Fluent計算流體力學(xué)數(shù)值模擬軟件，對該項目浮運施工中的2個關(guān)鍵節(jié)點進行模擬分析，計算沉管結(jié)構(gòu)在這2個關(guān)鍵節(jié)點的受力情況，并對阻力系數(shù)進行率定； 同時，采用MIKE流體計算軟件，在實際河床地形的情況下，對浮運航道的流速流場進行模擬，以便于采用理論公式進行水流阻力計算。

1　工程概況

紅谷隧道工程位于南昌大橋和八一大橋之間，距南昌大橋約1.4 km，距八一大橋約2.3 km，連接紅谷灘新區(qū)與老城區(qū)。紅谷隧道采用沉管法施工，隧道由12節(jié)管段拼接而成，這些管段有的長達115 m，寬30 m，高8.3 m。管段從上游干塢浮運至隧址，先后需經(jīng)過生米、朝陽和南昌3座大橋，最終在隧址回旋區(qū)進行轉(zhuǎn)體后沉放，浮運長度8 871 m，浮運線路見圖1。管段浮運不僅距離遠，而且面臨穿越橋梁時凈空小、可能碰撞橋梁等高風(fēng)險挑戰(zhàn)，因此精準(zhǔn)控制管段浮運過程的姿態(tài)至關(guān)重要。

圖1　浮運線路圖

沉管浮運過程存在2個高風(fēng)險控制點：

1)第1個風(fēng)險控制點為管段浮運過南昌大橋。南昌大橋凈跨僅81.8 m，扣除防撞設(shè)施后凈距僅為63.4 m，拖船可操作空間小且該處實測水流流向與拖運方向不一致，存在側(cè)向水流阻力，質(zhì)量達2.8萬t的管段一旦碰撞上橋墩，其后果非常嚴(yán)重，見圖2。該處采用5艘4 000 HP港作拖輪對沉管進行拖航，拖輪布置方式見圖3。

圖2　沉管浮運過南昌大橋示意圖(風(fēng)險控制點1)(單位： m)

Fig. 2Diagram of segment passes through Nanchang Bridge (risk control point No. 1) (m)

圖3　沉管浮運過南昌大橋拖輪布置方式示意圖

Fig. 3Arrangement of tugboats when segment passes through Nanchang Bridge

2)第2個風(fēng)險控制點為在隧址回旋區(qū)進行的管段轉(zhuǎn)體?；匦齾^(qū)離東岸圍堰最近點僅25 m，拖輪在回旋區(qū)可操作的空間較小。另外東岸圍堰施工縮短了隧址河道寬度，將加大河道水流速度，且管段在回旋區(qū)轉(zhuǎn)體過程中將會出現(xiàn)縱斷面迎流的情況。該處采用5艘4 000 HP港作拖輪配合2個質(zhì)量為170 t的錨塊和1個80 t墊塊對管段進行轉(zhuǎn)體，拖航步驟見圖4。

圖4　沉管于隧址回旋區(qū)轉(zhuǎn)體示意圖(風(fēng)險控制點2)

Fig. 4Diagram of segment turning around (risk control point No. 2)

2　沉管浮運過南昌大橋過程數(shù)值模擬

對浮運過程中第1個風(fēng)險控制點采用Fluent軟件進行CFD數(shù)值模擬分析，管段模型按尺寸最大的管段考慮，即115 m×30 m×8.3 m(長×寬×高)。

2.1計算工況

依據(jù)現(xiàn)場實測水流的流速和流向，結(jié)合實際施工過程，制定出計算工況，如表1所示。

表1　計算節(jié)點及計算工況一覽表

計算幾何模型見圖5。

(a) 未過橋

(b) 過橋1/2

(c) 完全過橋

2.2控制方程

本文運用數(shù)值仿真軟件Fluent中CFD方法為處理工具，選用的紊流模型為V. Yakhot等[11-12]于重整化群理論提出的RNGκ-ε紊流模型，具體方程如下。

連續(xù)性方程

ρu=0。

(1)

動量方程

。

(2)

動能方程

(3)

耗散方程

(4)

其中，

(5)

式中： ρ為流體密度； t為時間； u為速度場； F為高斯隨機力； k為動能； ε為耗散率； l為紊流長度； μ為無衰減湍流動力學(xué)黏度； μT為渦流黏滯系數(shù)； I為湍流強度； Re為雷諾數(shù)； σk、σε、Ce1、Ce2均為系數(shù)。

2.3邊界條件

考慮該處距離岸邊線尚有一定距離，且河床較為平整，故采用理想模型，認(rèn)為河床為同一標(biāo)高。邊界河流表面為自由水面，設(shè)定為開邊界，模型左右邊界設(shè)為對稱邊界，模型上游邊界為速度入口邊界，模型下游邊界為自由流出邊界。河流河床外圍邊壁、管段外壁等均為固壁，所有固壁處的節(jié)點為無滑移條件。沿水流方向取計算域500m，河寬方向取計算域240m，管段浮運干舷值取0.3m。

2.4計算結(jié)果與分析

由于篇幅原因，以下僅列出過橋1/2工況的計算結(jié)果，如圖6—9所示。

圖6　流速云圖(單位： m/s)

圖7　速度矢量圖(單位： m/s)

圖8　水壓力云圖(單位： Pa)

圖9　管段表面水壓力云圖(單位： Pa)

表2為管段浮運過南昌大橋時各工況下所受水流阻力匯總。實測南昌大橋段水流流向與航道方向存在夾角，在15～20°范圍內(nèi)，與計算工況基本一致。從數(shù)值模擬的結(jié)果可以看出，橫向水流阻力遠大于縱向水流阻力，且橫向水流阻力不僅會使管段平移，還會使管段產(chǎn)生旋轉(zhuǎn)運動。通過計算，可以校核浮運方案中拖輪所提供的橫向拖力能否滿足要求，給設(shè)計單位和施工單位在管段浮運過南昌大橋設(shè)計施工中提供參考。

表2　管段浮運過南昌大橋過程水流阻力匯總

3　沉管隧址回旋區(qū)轉(zhuǎn)體過程數(shù)值模擬

本節(jié)對浮運過程中第2個風(fēng)險控制點進行CFD數(shù)值模擬分析。河床地形見圖10。

圖10　河床地形圖(單位： m)

與第1個風(fēng)險控制點不同，回旋區(qū)由于東側(cè)圍堰、江心洲、丁壩、航道疏浚、淤積和隧道基槽開挖等因素，使得河流流向和流速十分復(fù)雜。因此，對該處而言，要想獲取符合實際情況的水流阻力，必須考慮河床地形的復(fù)雜情況。考慮到MIKE軟件具有強大的地形導(dǎo)入功能，因此，本節(jié)采用MIKE流體計算軟件計算回旋區(qū)的水流流向和流速情況，再結(jié)合理論水流阻力計算公式計算管段所受到的水流阻力作用。而按照《港口工程荷載規(guī)范》[13]，水流阻力的計算公式為

(6)

式中： Cw為水阻力系數(shù)； ρ為水的密度，取1t/m3； v為水流速度，m/s； A為迎流面積，m2； F為水流作用力，kN。

3.1水阻力系數(shù)率定

對于水流阻力計算公式中的水阻力系數(shù)Cw，本文采用Fluent數(shù)值模擬進行率定。數(shù)值模擬考慮沉管管段與水流方向的不同角度，分為夾角 0°、22.5°、45°、67.5°和90°，計算最低設(shè)計浮運水位(12.5m)及易出現(xiàn)洪水位(20m)考慮，計算流速取0.2、0.4、0.6m/s。圖11為流速0.6m/s、計算水位20m情況下的水壓力結(jié)果。

(a) 0°

(b) 22.5°

(c) 45°

(d) 67.5°

(e) 90°

Fig. 11Nephogramsofwaterpressuresunderdifferentanglesbetweensegmentandflowdirection(Pa)

由圖11可以看出，管段與水流方向夾角越大，管段上水壓力最大值越大，90°情況下最大水壓力為505 Pa，約為0°情況下最大水壓力的2倍。

提取管段上水流阻力大小，各工況水流阻力大小如表3所示。

表3　隧道管段水流阻力匯總

把數(shù)值模擬結(jié)果與水流阻力計算公式進行率定，采用最小二乘法進行擬合，得到適用于管段的水阻力系數(shù)Cw=1.82。

3.2回旋區(qū)流向和流速數(shù)值模擬分析

計算采用MIKE軟件和二次開發(fā)程序在考慮實際邊界條件的情況下以三維方式模擬非恒定流，分析流速流場。

當(dāng)模型建立后，通過實測數(shù)據(jù)對模型參數(shù)進行率定。實測數(shù)據(jù)采用2015年6月3日的實測流速資料，模型上游流量和下游水位給定為2015年6月3日實測值5 780 m3/s和16.6 m，并考慮隧址下游丁壩缺口。通過將回旋區(qū)數(shù)值模擬結(jié)果與測量點實測結(jié)果比較，見圖12。結(jié)果顯示，隧址各測站的計算值與實測值無論是流速還是流向均吻合良好，采用該模型能較好地預(yù)測隧址處的水流流速和流向。不同工況下流場計算結(jié)果見圖13—15。

藍色箭頭表示實測流速； 紅色箭頭表示數(shù)值模擬流速。

圖12回旋區(qū)流速、流向?qū)崪y值與計算值比較(單位： m/s)

Fig. 12Comparison between measured flow velocity and flow direction and calculated flow velocity and flow direction at segment turning area (m/s)

圖13　流速、流向彩圖(無沉管)(單位： m/s)

Fig. 13Nephogram of flow velocities and flow directions of the river (m/s)

圖14　流速、流向彩圖(縱斷面迎流)(單位： m/s)

Fig. 14Nephogram of flow velocities and flow directions during vertical section incident flowing (m/s)

圖15　流速、流向彩圖(橫斷面迎流)(單位： m/s)

Fig. 15Nephogram of flow velocities and flow directions during cross-section incident flowing (m/s)

數(shù)值模擬結(jié)果顯示，水流流速和流向在隧址回旋區(qū)十分復(fù)雜。在江心洲和圍堰的影響下，河道寬度變窄，流速迅速增加，并且在圍堰外邊緣存在旋渦現(xiàn)象； 西汊由于下游老丁壩阻隔水流，水流經(jīng)江心洲基槽通往東汊主航道，并且在越過基槽后，河床變窄，流速進一步增大； 而由于江心洲的擋水作用，江心洲后方水流流速較慢。

3.3水流阻力計算

采用MIKE軟件建立的河流河道模型能夠很好地預(yù)測隧址回旋區(qū)的水流流速和流向，結(jié)合采用經(jīng)數(shù)值模擬軟件率定后的水流阻力計算公式能夠較好地計算沉管隧道實際受到的水流阻力大小。以上游流量5 970 m3/s、水位16.5 m的工況為例，隧址處最大的流速為0.84 m/s，對應(yīng)的管段橫斷面迎流水流阻力為176 kN，縱斷面迎流水流阻力為674 kN。

4　結(jié)論與討論

本文采用Fluent和MIKE數(shù)值模擬軟件，模擬了南昌紅谷隧道浮運過程中的2個關(guān)鍵風(fēng)險控制點，計算分析了沉管管段在不同施工階段、不同水流流向和流速情況下所受水流阻力大小，得到了以下結(jié)論：

1)由于橋墩及防撞箱對側(cè)向水流有一定阻擋作用，因此隨著沉管管段浮運過南昌大橋，其作用在沉管上的水流阻力先變小后變大，并且流速對水流阻力影響很大。

2)依據(jù)實測的水流流向，經(jīng)計算，對于沉管管段浮運過南昌大橋時，管段受橫向水流阻力較大，流速大時施工存在一定風(fēng)險。

3)采用數(shù)值模擬手段，對水流阻力計算公式中的水流阻力系數(shù)Cw進行率定，使之適用于紅谷隧道浮運工程。

4)采用MIKE計算軟件對整個浮運航道進行模擬，考慮實際疏浚后的河床地形，計算了不同水文條件下河道的流場，數(shù)值模擬結(jié)果表明，受河床地形影響，回旋區(qū)的水流流速十分復(fù)雜。

5)依據(jù)MIKE軟件的流速流向計算結(jié)果，結(jié)合采用Fluent數(shù)值模擬率定的規(guī)范經(jīng)驗計算公式，為管段在復(fù)雜周邊環(huán)境中水阻流力計算提供一套方法。

沉管管段在周邊存在控制性建筑或流場復(fù)雜的情況下浮運往往存在很高的風(fēng)險，獲取管段所受水流阻力大小對指導(dǎo)浮運施工具有十分重要的意義。由于筆者精力有限，未能對更多的水文工況進行分析，對施工過程的模擬也未進一步細分，接下來還需進行更全面的研究。

[1]孫鈞. 海底隧道工程設(shè)計施工若干關(guān)鍵技術(shù)的商榷[J]. 巖石力學(xué)與工程學(xué)報,2006,25(8): 1513-1521. (SUN Jun. Discussion on some key technical issues for design and construction of undersea tunnels[J]. Chinese Journal of Rock Mechanics and Engineering,2006,25(8): 1513-1521.(in Chinese))

[2]王夢恕. 水下交通隧道發(fā)展現(xiàn)狀與技術(shù)難題： 兼論“臺灣海峽海底鐵路隧道建設(shè)方案”[J]. 巖石力學(xué)與工程學(xué)報,2008,27(11): 2161-2172.(WANG Mengshu. Current developments and technical issues of underwater traffic tunnel: Discussion on construction scheme of Taiwan Strait undersea railway tunnel[J].Chinese Journal of Rock Mechanics and Engineering,2008,27(11): 2161-2172.(in Chinese))

[3]王秋林,崔玉國,陳旺. 南昌紅谷沉管隧道預(yù)制關(guān)鍵設(shè)備選型及配套技術(shù)[J]. 隧道建設(shè),2015,35(增刊2): 1345-1350. (WANG Qiulin, CUI Yuguo, CHEN Wang. Type selection and matching of key equipment for tube element prefabrication of Honggu Immersed Tunnel in Nanchang[J].Tunnel Construction,2015, 35(S2): 1345-1350.(in Chinese))

[4]孫鈞. 論跨江越海建設(shè)隧道的技術(shù)優(yōu)勢與問題[J]. 隧道建設(shè),2013,33(5): 337-342. (SUN Jun. Study of technological advantages and problem in construction of tunnels crossing rivers and seas[J]. Tunnel Construction,2013,33(5): 337-342.(in Chinese))

[5]呂衛(wèi)清,應(yīng)宗權(quán),蘇林王,等. 沉管管節(jié)浮運過程中波浪附加阻力的水動力學(xué)分析[J]. 水運工程,2013(增刊1): 86-91. (LYU Weiqing, YING Zongquan, SU Linwang, et al. Hydrodynamic analysis of added resistance in waves of immersed tunnel elements during floating transportation[J].Port & Water Engineering, 2013(S1): 86-91.(in Chinese))

[6]劉慧. 洲頭咀沉管隧道施工穩(wěn)定性與可操縱性研究[D].廣州： 華南理工大學(xué),2012. (LIU Hui. Study of stability and maneuverability of Zhoutouzui immersed tube tunnel[D].Guangzhou: South China University of Technology,2012.(in Chinese))

[7]ZHANG Dexin, WANG Xingquan. Experiments of hydrodynamics and stability of immersed tube tunnel on transportation and immersing[J]. Journal of Hydrodynamics, 2001, 13(2): 121-126.

[8]蔣龍皎. 隧道沉管在淺水航道中浮運時的阻力性能與浮態(tài)研究[D].廣州：華南理工大學(xué),2011. (JIANG Longjiao. Study of the resistance and floatation performance of the immersed tube tunnel towing in the shallow water[D].Guangzhou: South China University of Technology,2011.(in Chinese))

[9]周敏. 沉管隧道浮運水動力學(xué)的模型試驗研究[D].北京： 北京交通大學(xué),2014.(ZHOU Min. Experimental research on the floating transportation hydrodynamics of the immersed tube tunnel element[D].Beijing: Beijing Jiaotong University,2014. (in Chinese))

[10]羅甜. 沉管管節(jié)浮運動力特性研究[D].武漢： 武漢理工大學(xué),2013.(LUO Tian.Research on motion characteristics in floating transportation of immersed tunnel element[D].Wuhan: Wuhan University of Technology,2013.(in Chinese))

[11]Yakhot V, Orszag S A. Renormalization group analysis of turbulent I: Basic theory[J]. Journal of Scientific Computing, 1986,1(1): 3-51.

[12]Yakhot V, Orszag S A,Thangam S,et al. Development of turbulent models for shear flows by a double expansion technique[J].Physics of Fluids A Fluid Dynamics, 1992,4(7): 1510-1520.

[13]港口工程荷載規(guī)范： JTS 144-1—2010[S]. 北京： 人民交通出版社，2010. (Load code for harbour engineering: JTS 144-1—2010[S]. Beijing: China Communications Press, 2010. (in Chinese))

Numerical Analysis of Characteristics of Water Resistance during Floating Transportation of Immersed Tunnel Segment and Flow Field of Segment Turning Area

HE Zhimin1, LI Jiufu1, LIU Huikang2, *, ZHANG Yucheng3, YU Xiaoqiang1

(1.NanchangMunicipalPublicGroup,Nanchang330000,Jiangxi,China; 2.GuangzhouMetroDesign&ResearchInstituteCo.,Ltd.,Guangzhou510000,Guangdong,China; 3.GuangdongResearchInstituteofWaterResourcesandHydropower,Guangzhou510610,Guangdong,China)

For segment floating transportation, Nanchang Bridge and segment turning area are the key points for Honggu Immersed Tunnel in Nanchang. The characteristics of water resistance during floating transportation of the segment are analyzed by Fluent and MIKE fluid dynamics software, so as to guarantee the safe floating transportation. The characteristics of water resistance when the segment passes through Nanchang Bridge are obtained. The characteristics of water resistance at segment turning area are obtained by analyzing water resistance coefficientCwand the flow field under different hydrological conditions. The results can provide reference for water resistance calculation during floating transportation of immersed tunnel segment.

Honggu Immersed Tunnel; segment floating transportation; water resistance; numerical simalation

2016-05-19;

2016-08-16

何志敏(1964—)，男，福建光澤人，1984年畢業(yè)于江蘇工學(xué)院，修造專業(yè)，本科，工程師，現(xiàn)從事工程項目管理工作。E-mail: hezhimin88888@163.com。*通訊作者： 劉惠康， E-mail: 616075123@qq.com。

10.3973/j.issn.1672-741X.2016.09.001

U 455.46

A

1672-741X(2016)09-1023-07