黏度時(shí)變雙液漿盾構(gòu)壁后注漿過(guò)程數(shù)值模擬研究

李文廣,王東武,張晨光,陳敬軒,曹宏濤

(1.中鐵十八局集團(tuán)第五工程有限公司,天津 300450;2.浙江工業(yè)大學(xué) 土木工程學(xué)院,浙江 杭州 310023)

盾尾壁后注漿工作多數(shù)基于工程經(jīng)驗(yàn),缺乏科學(xué)、合理的理論指導(dǎo),導(dǎo)致注漿工程施工效果難以把控,盾構(gòu)隧道變形[1-2]、破損[3]和滲漏等工程問(wèn)題頻頻發(fā)生,維修成本高昂[4]。因此,深入探究各類(lèi)注漿施工參數(shù)與注漿擴(kuò)散形態(tài)的關(guān)系極具工程指導(dǎo)意義。目前,相關(guān)學(xué)者在盾構(gòu)壁后注漿漿液填充擴(kuò)散機(jī)理方面做了一定的研究工作。白云等[5]以單位時(shí)間內(nèi)形成的盾尾空隙為漿液充填的橫斷面,推導(dǎo)出盾尾空隙橫斷面內(nèi)漿液壓力分布模型;李志明等[6]進(jìn)一步將盾構(gòu)同步注漿的填充分解為環(huán)向填充和縱向填充兩個(gè)獨(dú)立的階段,沒(méi)有考慮漿液縱橫向擴(kuò)散的相互影響,得到了環(huán)向填充及縱向填充的力學(xué)模型及計(jì)算方法;李培楠等[7]則建立了盾尾同步注漿縱環(huán)向整體擴(kuò)散理論模型;胡長(zhǎng)明等[8]基于簡(jiǎn)化假設(shè)條件分析了盾尾環(huán)向流體微元空間力學(xué)特性,建立了漿液在盾尾空間中的二元流環(huán)形填充擴(kuò)散模型;茍長(zhǎng)飛[9]建立了單液漿同步注漿時(shí)的扇形充填擴(kuò)散模型;陳少波[10]依托南京軌道交通L5項(xiàng)目,分別采用環(huán)向獨(dú)立模型和整體填充擴(kuò)散模型對(duì)漿液壓力分布模式進(jìn)行計(jì)算分析,證實(shí)了整體填充擴(kuò)散模型的準(zhǔn)確性?xún)?yōu)勢(shì)。然而，目前相關(guān)學(xué)者對(duì)盾構(gòu)注漿填充擴(kuò)散的縱橫向耦合擴(kuò)散研究甚少,值得進(jìn)一步探究。

就材料而言,盾構(gòu)壁后注漿主要采用以下兩種漿液:以水泥為主的單液型漿液和以水泥-水玻璃(C-S)為主的雙液型漿液。單液漿雖然流動(dòng)性能好,利于填充擴(kuò)散,且原料便宜,但是其往往存在稀釋離析或篩濾流失等問(wèn)題。而C-S雙液漿中水玻璃類(lèi)漿液起到速凝作用,該漿液有著流失少、泵送時(shí)材料分離度小、短期和長(zhǎng)期強(qiáng)度均比較大、體積收縮低以及凝結(jié)時(shí)間短等優(yōu)勢(shì)[11-12],被越來(lái)越多地應(yīng)用于盾尾同步注漿施工。然而雙液漿凝結(jié)時(shí)間快(黏度時(shí)變特性不可忽略),在壁后注漿往往留有未充填的空隙,不能保證100%填充,會(huì)影響盾構(gòu)注漿施工安全。目前的注漿研究鮮有考慮雙液漿黏度時(shí)變性對(duì)注漿填充的影響[13]。因此,筆者深入探討雙液漿黏度時(shí)變性和施工參數(shù)對(duì)漿液縱橫向耦合擴(kuò)散的影響?；诙辔锢韴?chǎng)耦合有限元軟件COMSOL,模擬雙液漿在盾尾間隙的三維局部填充擴(kuò)散過(guò)程。通過(guò)控制變量,系統(tǒng)分析了漿液黏度時(shí)變性、入口注漿壓力和注漿孔位置對(duì)注漿填充擴(kuò)散特征的影響,以期為注漿工程采取合理的注漿策略提供參考。

1 有限元模型創(chuàng)建

COMSOL軟件應(yīng)用于非牛頓流體流動(dòng)的模擬運(yùn)算已經(jīng)相對(duì)成熟[14-16]。劉人太等[14]針對(duì)水泥-水玻璃漿液與高聚物改性水泥漿液,應(yīng)用COMSOL建立動(dòng)水條件下裂隙注漿擴(kuò)散的數(shù)值模型,研究動(dòng)水條件下裂隙注漿擴(kuò)散規(guī)律，分析不同黏度時(shí)變性、初始動(dòng)水流速與注漿速率對(duì)注漿擴(kuò)散的影響,并將數(shù)值模擬結(jié)果與模型試驗(yàn)進(jìn)行對(duì)比,驗(yàn)證了數(shù)值模擬方法的有效性;高芙蓉[15]以脲醛樹(shù)脂漿液為研究對(duì)象,在COMSOL中還原了室內(nèi)注漿模擬中漿液的黏度時(shí)變性,并研究了漿液在含水砂層中的擴(kuò)散機(jī)制,驗(yàn)證了脲醛樹(shù)脂漿液對(duì)流砂的治理效果;楊志全等[16]依托COMSOL平臺(tái),通過(guò)二次開(kāi)發(fā)得到考慮多孔介質(zhì)迂回曲折效應(yīng)的冪律流體滲透注漿的模擬程序,并以此開(kāi)展不同水灰比水泥漿液在礫石土體中滲透注漿過(guò)程的數(shù)值模擬。

筆者基于漿液和空氣兩相流體的相互關(guān)系,且考慮兩流體在擴(kuò)散全程均為層流狀態(tài),將COMSOL流體力學(xué)模塊中的層流物理場(chǎng)與兩相流物理場(chǎng)進(jìn)行耦合,對(duì)不同工況的盾尾同步注漿進(jìn)行數(shù)值模擬計(jì)算,兩相的界面追蹤計(jì)算采用相場(chǎng)方法,對(duì)漿液的空隙填充過(guò)程進(jìn)行數(shù)值模擬,漿液和空氣在整個(gè)計(jì)算區(qū)域內(nèi)的分布采用體積分?jǐn)?shù)法進(jìn)行表征。

1.1 物理場(chǎng)控制方程

1.1.1 層流接口的控制方程

模型基于Navier-Stokes方程模擬流體的質(zhì)量和動(dòng)量傳遞,為了考慮表面張力的影響,在控制方程中包含表面張力項(xiàng),因此采用的Navier-Stokes方程為

(1)

(2)

1.1.2 相場(chǎng)接口的控制方程

在相場(chǎng)接口中,兩相流動(dòng)力學(xué)由Cahn-Hilliard方程控制。該方程用于跟蹤不相溶的兩種液體的擴(kuò)散界面。擴(kuò)散界面定義為無(wú)量綱相場(chǎng)變量φ從-1到1的區(qū)域。在COMSOL中求解時(shí),Cahn-Hilliard方程分成兩個(gè)方程,即

(3)

(4)

式中:u為流體速度矢量,m/s;λ為混合能量密度,N;ε為界面厚度參數(shù),m;γ為遷移率,(m3·s)/kg;ψ為相場(chǎng)助變量。

在相場(chǎng)接口中,各種流體的體積分?jǐn)?shù)為

(5)

(6)

式中:Vf1為流體1的體積分?jǐn)?shù);Vf2為流體2的體積分?jǐn)?shù);φ為無(wú)量綱相場(chǎng)變量。

在筆者模型中,空氣定義為流體1,漿液定義為流體2。后處理時(shí),以?xún)上囿w積分?jǐn)?shù)各為0.5的等值面為相界面,追蹤漿液的實(shí)時(shí)擴(kuò)散面。

1.1.3 多物理場(chǎng)的耦合方程

多物理場(chǎng)耦合特征通過(guò)以下表達(dá)式定義混合物跨兩相界面平滑變化的密度和黏度,分別為

ρ=ρ2+(ρ1-ρ2)Vf1

(7)

μ=μ2+(μ1-μ2)Vf1

(8)

式中:μ為流體黏度,Pa·s;下角標(biāo)1,2分別表示流體1和流體2;其他符號(hào)含義同上。

1.2 幾何創(chuàng)建、邊界條件及基本參數(shù)

參考杭州下沙隧道工程盾構(gòu)施工段的相關(guān)參數(shù),構(gòu)建局部注漿模型?；诨鸀橹钡暮?jiǎn)化思想,將弧形的局部盾尾間隙創(chuàng)建為7 000 mm×2 000 mm×200 mm的長(zhǎng)方體空腔,并在前側(cè)面的居中位置設(shè)定直徑為80 mm的同步注漿管口,長(zhǎng)方體空腔的兩側(cè)設(shè)為無(wú)壓出流邊界。其他壁面則設(shè)為無(wú)滑移壁(圖1)。采用自由四面體網(wǎng)格進(jìn)行劃分,網(wǎng)格劃分示意圖見(jiàn)圖2。

圖1 三維幾何建模Fig.1 Three-dimensional geometric modeling

圖2 網(wǎng)格劃分(單位:mm)Fig.2 Mesh division (unit: mm)

模型計(jì)算初始狀態(tài)設(shè)定局部盾尾間隙內(nèi)填充相為空氣,而注漿管道內(nèi)的初始填充相為漿液,因此初始的兩相界面位于注漿管體與待填充腔體的銜接面上。在恒定注漿入口壓力條件下,漿液逐漸擴(kuò)散并填充盾尾間隙,完成局部間隙的填充。

創(chuàng)建如圖3所示的3條觀測(cè)壓力場(chǎng)與速度場(chǎng)的觀測(cè)線(xiàn),所有觀測(cè)線(xiàn)均處在盾尾間隙水平對(duì)稱(chēng)面上。其中,觀測(cè)線(xiàn)1通過(guò)注漿孔圓心且平行于側(cè)面出口;觀測(cè)線(xiàn)2平行于盾尾間隙前側(cè)壁面且與其相距100 mm;觀測(cè)線(xiàn)3平行于右側(cè)出口且與其相距100 mm。觀測(cè)線(xiàn)1將表征注漿入口正前方漿液注入方向的壓力場(chǎng)分布與速度場(chǎng)分布;觀測(cè)線(xiàn)2將表征與注漿方向垂直的橫向壓力場(chǎng)分布與速度場(chǎng)分布;觀測(cè)線(xiàn)3的觀測(cè)結(jié)果則反映了出口邊界附近的壓力場(chǎng)與速度場(chǎng)的分布情況。

圖3 觀測(cè)線(xiàn)布置Fig.3 Observation wire cloth

1.3 計(jì)算工況

為研究漿液黏度時(shí)變性、入口注漿壓力以及注漿孔位置對(duì)局部填充注漿的影響,共設(shè)定6種計(jì)算工況以供后期對(duì)比分析?？紤]到漿液在配制完成后還需一段時(shí)間從調(diào)配位置運(yùn)輸?shù)蕉芪瞇17],設(shè)其黏度μ=0.018 64(t+50)2.066,單位為Pa·s;工況1,2,3采用不同的入口注漿壓力,分別為100,200,300 kPa;工況1,5,6將盾尾間隙局部模型與重力的夾角設(shè)定為90°,60°,30°,從而模擬不同位置注漿孔的漿液擴(kuò)散特征;工況4不考慮漿液的黏度時(shí)變性,認(rèn)為漿液黏度值始終為t=0 s時(shí)的黏度值60 Pa·s,將其作為工況1的對(duì)照組,探究漿液黏度時(shí)變性對(duì)填充擴(kuò)散的影響。綜合以上各類(lèi)參數(shù)的設(shè)定,通過(guò)控制變量的研究方法得出不同因素對(duì)漿液填充的影響。

表1 算例參數(shù)Table 1 Example parameter setting

2 漿液擴(kuò)散典型工況分析

2.1 漿液擴(kuò)散形態(tài)分析

當(dāng)盾尾間隙局部模型與重力的夾角為90°時(shí),在不同漿液黏度及入口注漿壓力條件下,漿液的擴(kuò)散形態(tài)具有相似性。因此,選擇典型的工況3進(jìn)行漿液的擴(kuò)散形態(tài)分析,取6個(gè)代表時(shí)刻的漿液體積分?jǐn)?shù)剖面圖反映漿液的整個(gè)擴(kuò)散過(guò)程。圖4為工況3漿液擴(kuò)散形態(tài)圖,黑色輪廓線(xiàn)代表兩相的分界面,即漿液的擴(kuò)散邊界。由于注漿條件、幾何尺寸及邊界條件的對(duì)稱(chēng)性,漿液的擴(kuò)散形態(tài)整體呈現(xiàn)出對(duì)稱(chēng)性。擴(kuò)散初期,漿液自注漿入口開(kāi)始以扇形形態(tài)逐漸擴(kuò)散,當(dāng)漿液擴(kuò)散面抵觸到后壁時(shí),漿液擴(kuò)散面開(kāi)始分成左右對(duì)稱(chēng)的兩個(gè)界面向兩側(cè)開(kāi)口繼續(xù)擴(kuò)散,直至填充整個(gè)盾尾間隙,完成局部填充。

圖4 漿液擴(kuò)散形態(tài)(工況3)Fig.4 Slurry diffusion pattern (case 3)

2.2 速度場(chǎng)與壓力場(chǎng)變化規(guī)律分析

2.2.1 速度場(chǎng)特性分析

以工況3測(cè)線(xiàn)上的速度場(chǎng)變化為例(圖5),說(shuō)明注漿過(guò)程中流體速度場(chǎng)隨時(shí)間的發(fā)展規(guī)律。考慮到漿液的黏度時(shí)變性,即漿液黏度隨時(shí)間推進(jìn)不斷增大。在入口壓力恒定的條件下,入口處注漿流速隨著時(shí)間推進(jìn)銳減,觀測(cè)線(xiàn)1與2中其他觀測(cè)點(diǎn)的速度也隨時(shí)間的推進(jìn)而減小。此外,隨著漿液擴(kuò)散邊界的發(fā)展,相同點(diǎn)后側(cè)的漿液更為廣泛,且黏度更大,從而提高了對(duì)該點(diǎn)的漿液推進(jìn)阻力,這也是導(dǎo)致速度場(chǎng)值隨著時(shí)間推進(jìn)而逐漸減小的原因。觀測(cè)線(xiàn)1與觀測(cè)線(xiàn)2沿著y軸和x軸的速度場(chǎng)值分布基本相似,速度下降梯度在靠近注漿處最大,隨著遠(yuǎn)離注漿口的位置逐漸減小。

圖5 工況3速度場(chǎng)Fig.5 Velocity field of case 3

2.2.2 壓力場(chǎng)特性分析

以工況3的3條測(cè)線(xiàn)上的壓力分布為例(圖6),說(shuō)明注漿過(guò)程中壓力場(chǎng)隨時(shí)間的變化情況。由圖6可知:觀測(cè)線(xiàn)1與2的壓力峰值均位于注漿口處,流體壓力以注漿口為中心向四周遞減。這是因?yàn)殡S著注漿范圍增大,流體壓力因內(nèi)部剪切和邊界摩擦沿程損失?？紤]到漿液的黏度時(shí)變性,流速隨著黏度增大而減小(圖5)。根據(jù)流體力學(xué)中壓力衰減與流速關(guān)系,流體流速越小,流體壓降梯度越小[18]。因此,觀測(cè)線(xiàn)1遠(yuǎn)離注漿口的壓力值隨時(shí)間推進(jìn)略有提高。觀測(cè)線(xiàn)2的壓力值呈對(duì)稱(chēng)分布,流體壓力隨距離的變化規(guī)律和時(shí)間的推進(jìn)規(guī)律與觀測(cè)線(xiàn)1基本相同,即流體壓力隨著距注漿口距離的增大而不斷減小。隨著時(shí)間的推進(jìn),流速降低,因此沿程壓力損失逐漸減小,遠(yuǎn)離注漿口的壓力略有提高。在觀測(cè)線(xiàn)2的兩端,由于出流邊界的壓力始終為0,導(dǎo)致壓力場(chǎng)最終均趨于0值。觀測(cè)線(xiàn)3表征了沿著出口邊界附近的壓力分布規(guī)律。由于兩側(cè)出口附近的區(qū)域與注漿入口相距較遠(yuǎn),且緊挨著出口邊界,因此這些區(qū)域表征為相對(duì)均勻的氣體低壓場(chǎng)。通過(guò)對(duì)觀測(cè)線(xiàn)3不同時(shí)刻壓力場(chǎng)的對(duì)比分析,可以得出流體壓力隨時(shí)間的推進(jìn)(即隨漿液的擴(kuò)展)而提高。在t=160 s時(shí),由于漿液擴(kuò)展到了觀測(cè)線(xiàn)3的位置,在該時(shí)刻的局部區(qū)段壓力(漿液壓力)明顯增大。

圖6 工況3壓力場(chǎng)Fig.6 Pressure field of case 3

3 不同參數(shù)對(duì)漿液擴(kuò)散規(guī)律的影響

3.1 不同黏度模型漿液擴(kuò)散規(guī)律

為研究雙液漿黏度時(shí)變性對(duì)漿液擴(kuò)散的影響,對(duì)工況1,4的模擬結(jié)果進(jìn)行對(duì)比分析,結(jié)果如圖7所示。在整個(gè)填充擴(kuò)散過(guò)程中,相同時(shí)刻下不考慮漿液黏度時(shí)變性的漿液擴(kuò)散面始終比考慮漿液黏度時(shí)變性的漿液擴(kuò)散面發(fā)展更為廣泛。

圖7 工況1,4漿液擴(kuò)散面(t=30 s)Fig.7 Slurry diffusion surface of case 1 and 4 (t=30 s)

進(jìn)一步分析模擬結(jié)果中盾尾間隙觀測(cè)線(xiàn)數(shù)據(jù),探究黏度時(shí)變性對(duì)速度場(chǎng)和壓力場(chǎng)的影響作用。工況1,4的觀測(cè)線(xiàn)對(duì)應(yīng)的速度分布如圖8所示,壓力

圖8 工況1,4速度場(chǎng)(t=30 s)Fig.8 Velocity field of case 1 and 4 (t=30 s)

分布如圖9所示。因?yàn)殡p液漿黏度的增加使得漿液流動(dòng)過(guò)程中受到的摩擦阻力和內(nèi)部剪切增大,造成流體動(dòng)能在傳遞過(guò)程中的損耗提高,最終導(dǎo)致漿液壓力場(chǎng)和速度場(chǎng)的減小[19]。因此,黏度與速度場(chǎng)、壓力場(chǎng)呈負(fù)相關(guān)。在考慮漿液黏度時(shí)變性的工況中，隨著時(shí)間推進(jìn)漿液黏度增加,故在相同時(shí)刻下工況4的速度場(chǎng)與壓力場(chǎng)比工況6的速度場(chǎng)與壓力場(chǎng)更小。因此,在考慮漿液黏度時(shí)變性的工況中,漿液的整體擴(kuò)散過(guò)程先快后慢。由此看來(lái),雙液漿的黏度時(shí)變性對(duì)漿液的填充擴(kuò)散發(fā)展速度影響顯著,在數(shù)值模擬和理論計(jì)算中均不可忽略。進(jìn)一步對(duì)漿液擴(kuò)散面在y軸方向上的推進(jìn)距離d進(jìn)行量化分析,結(jié)果如圖10所示,經(jīng)數(shù)據(jù)擬合得到工況1,4的d—t曲線(xiàn)的函數(shù)方程。考慮漿液為常黏度工況的d隨t增長(zhǎng)迅速;考慮漿液黏度時(shí)變性工況的漿液因黏度的不斷提高,擴(kuò)散阻力不斷增大,使得d隨t的發(fā)展由快轉(zhuǎn)慢,這正是實(shí)際工程中注漿不充分的主要原因。

圖9 工況1,4壓力場(chǎng)(t=30 s)Fig.9 Pressure field of case 1 and 4 (t=30 s)

圖10 工況1,4的擴(kuò)散距離擬合曲線(xiàn)Fig.10 Diffusion distance fitting curve of case 1 and 4

3.2 不同入口注漿壓力下的漿液擴(kuò)散規(guī)律

如圖11所示,對(duì)比入口注漿壓力分別為100,200,300 kPa時(shí)漿液的填充擴(kuò)散結(jié)果。通過(guò)對(duì)比分析可得：在相同時(shí)刻下,隨著注漿壓力的提高,漿液填充擴(kuò)散發(fā)展越充分。

圖11 工況1,2,3漿液擴(kuò)散面(t=30 s)Fig.11 Dispersion surface of slurry in case 1, 2and 3 (t=30 s)

進(jìn)一步分析漿液在填充盾尾間隙時(shí)速度場(chǎng)和壓力場(chǎng)的不同,研究不同入口注漿壓力對(duì)漿液擴(kuò)散的影響。取工況1,2,3在t=30 s這一典型時(shí)刻的觀測(cè)線(xiàn)數(shù)據(jù)進(jìn)行分析,結(jié)果如圖12,13所示。由圖12，13可以得出:不同入口壓力并不改變盾尾間隙內(nèi)全局的速度場(chǎng)和壓力場(chǎng)變化趨勢(shì),只是影響速度和壓力的數(shù)值大小,且入口注漿壓力與速度、壓力均呈正相關(guān)。以觀測(cè)線(xiàn)2為例,隨著注漿壓力的增大,觀測(cè)線(xiàn)2的速度場(chǎng)隨之提高,但速度場(chǎng)分布始終呈現(xiàn)對(duì)稱(chēng)的倒V型走勢(shì)。注漿壓力由100 kPa增至200 kPa,對(duì)應(yīng)測(cè)線(xiàn)2上速度峰值增量為0.49 m/s。注漿壓力由200 kPa增至300 kPa,對(duì)應(yīng)測(cè)線(xiàn)2上速度峰值增量為0.57 m/s。隨著注漿壓力的增大,在相同的注漿壓力增幅下,速度場(chǎng)增量略有提高。隨著注漿壓力的增大,觀測(cè)線(xiàn)2的壓力峰值分別為2 911,7 029,12 466 Pa?？梢?jiàn)壓力場(chǎng)與注漿壓力也呈正相關(guān),且隨著注漿壓力的增大,在相同的注漿壓力增幅下,壓力場(chǎng)增幅略有提高。通過(guò)對(duì)比分析所有觀測(cè)線(xiàn)的壓力場(chǎng)、速度場(chǎng),可得入口注漿壓力與研究域內(nèi)的速度、壓力均呈正相關(guān)。值得注意的是:由于工況1注漿壓力較小,漿液發(fā)展相對(duì)較慢,在t=30 s時(shí),盾尾間隙內(nèi)僅在注漿口附近存在少量漿液,而觀測(cè)線(xiàn)3距離注漿口較遠(yuǎn)且緊挨著零壓出口邊界,因此圖13(c)的100 kPa的壓力曲線(xiàn)表征為相對(duì)均勻的氣體低壓場(chǎng)。

圖12 工況1,2,3速度場(chǎng)(t=30 s)Fig.12 Velocity field of case 1, 2 and 3 (t=30 s)

圖13 工況1,2,3壓力場(chǎng)(t=30 s)Fig.13 Pressure field of case 1, 2 and 3 (t=30 s)

3.3 不同位置注漿孔漿液擴(kuò)散規(guī)律

通過(guò)改變局部填充模型與重力的夾角,模擬不同位置注漿孔的局部填充擴(kuò)散過(guò)程。當(dāng)盾尾間隙局部模型與重力的夾角小于90°時(shí),漿液擴(kuò)散形態(tài)因重力作用體現(xiàn)出明顯的不對(duì)稱(chēng)性,如圖14所示,在工況6條件下,漿液右側(cè)擴(kuò)散面填充速度明顯大于左側(cè)擴(kuò)散面。起初隨著時(shí)間的推進(jìn),漿液迅速填充右側(cè)盾尾間隙,隨后漿液黏度的提高使得其可注性不斷減小,漿體擴(kuò)散面在注漿后期趨于穩(wěn)定,因而左側(cè)的空隙始終沒(méi)有得到漿液填充,該模擬結(jié)果詮釋了實(shí)際注漿工程中盾尾間隙填充不滿(mǎn)工況的發(fā)展過(guò)程。

圖14 漿液擴(kuò)散形態(tài)(工況6)Fig.14 Slurry diffusion pattern (case 6)

如圖15所示,對(duì)比不同重力夾角工況下模型內(nèi)漿液擴(kuò)散的發(fā)展形態(tài)可知:在相同時(shí)刻下,隨著夾角的減小,漿液擴(kuò)散形態(tài)的不對(duì)稱(chēng)性就越為明顯。進(jìn)一步觀測(cè)模型內(nèi)速度場(chǎng)和壓力場(chǎng)的不同,探究不同位置注漿孔的漿液擴(kuò)散規(guī)律,結(jié)果如圖16,17所示。觀測(cè)線(xiàn)2的速度場(chǎng)隨著盾尾間隙與重力夾角的減小,速度數(shù)值出現(xiàn)左低右高的非對(duì)稱(chēng)性,且非對(duì)稱(chēng)性與重力夾角呈負(fù)相關(guān)。分析認(rèn)為:隨著重力夾角的減小,重力在盾尾間隙平面的分力隨之增大,從而使得右側(cè)漿液下滑,擴(kuò)散加速發(fā)展;左側(cè)漿液上爬,擴(kuò)散減緩發(fā)展。由圖17可以看出與速度場(chǎng)對(duì)應(yīng)的壓力場(chǎng)的不對(duì)稱(chēng)性也與夾角成負(fù)相關(guān)。該模擬結(jié)果說(shuō)明：盾尾不同位置注漿孔呈現(xiàn)不同的漿液擴(kuò)散形態(tài),在實(shí)際施工時(shí),應(yīng)采取不同的控制參數(shù)。

圖15 工況1,5,6漿液擴(kuò)散面(t=30 s)Fig.15 Dispersion surface of slurry in case 1, 5 and 6(t=30 s)

圖16 工況1,5,6觀測(cè)線(xiàn)2速度場(chǎng)(t=30 s)Fig.16 Velocity field of observation line 2 in case 1,5 and 6 (t=30 s)

圖17 工況1,5,6觀測(cè)線(xiàn)2壓力場(chǎng)(t=30 s)Fig.17 Pressure field of observation line 2 in case 1,5 and 6 (t=30 s)

4 結(jié) 論

通過(guò)控制變量法,組合對(duì)比了不同工況下,盾尾壁后注漿局部模型內(nèi)雙液漿擴(kuò)散面的發(fā)展過(guò)程以及速度場(chǎng)與壓力場(chǎng)的變化規(guī)律,探究了盾尾壁后漿液的填充擴(kuò)散機(jī)制,得出了以下結(jié)論:1) 漿液擴(kuò)散速度與漿液的黏度呈負(fù)相關(guān),而與入口注漿壓力呈正相關(guān);2) 由于雙液漿具有黏度時(shí)變性,其擴(kuò)散進(jìn)程表現(xiàn)出顯著的先快后慢的趨勢(shì),導(dǎo)致施工過(guò)程中盾尾間隙填充不滿(mǎn),因此在實(shí)際工程中應(yīng)適當(dāng)增加注漿壓力,保證盾尾間隙填充率;3) 在盾構(gòu)壁后注漿過(guò)程中,壓力場(chǎng)和速度場(chǎng)的峰值出現(xiàn)在注漿口附近,且兩者數(shù)值的大小基本都與注漿口距離呈負(fù)相關(guān);4) 在盾尾不同位置的注漿孔中,頂部和底部注漿孔的漿液擴(kuò)散形態(tài)基本對(duì)稱(chēng),側(cè)邊注漿孔的漿液同時(shí)受壓力和重力作用,擴(kuò)散形態(tài)不對(duì)稱(chēng),絕大多數(shù)漿液沿重力方向向下擴(kuò)散并堆積,上側(cè)區(qū)域填充不

滿(mǎn),因此在實(shí)際工程中可以適當(dāng)增加頂部注漿孔的注漿量或上側(cè)注漿孔的排布密度,以保證注漿效果。