摘 要：為了解決泥水環(huán)流系統(tǒng)排渣困難，避免出現(xiàn)排渣口卵石滯排以及環(huán)流系統(tǒng)泵供能不足的現(xiàn)象，研究了渣土在泥漿管路中的堆積量和水力坡降變化規(guī)律。通過(guò)離散元法構(gòu)建排漿水平管段石渣輸送數(shù)值模型，研究泥漿流量、泥漿相對(duì)密度、石渣質(zhì)量流率對(duì)石渣堆積量、泥漿管路水力坡降的影響。結(jié)果表明：現(xiàn)有施工參數(shù)下大于20 mm粒徑的石渣主要在水平直管底部做推移運(yùn)動(dòng);泥漿流量的增大能夠明顯減少管路系統(tǒng)中的石渣堆積量，且能夠大大削弱泥漿相對(duì)密度和質(zhì)量流率對(duì)石渣堆積量的影響，管路水力坡降呈先減小后增大的趨勢(shì)，在流量為763.41 m3/h時(shí)取得最小值;石渣堆積量隨泥漿相對(duì)密度的增大呈下降趨勢(shì)，管路水力坡降呈增大趨勢(shì);石渣質(zhì)量流率的增大會(huì)增大管路中的石渣堆積量和管路的水力坡降。

關(guān)鍵詞：泥水盾構(gòu);環(huán)流系統(tǒng);水平直管;堆積量;管路水力坡降;質(zhì)量流率

中圖分類號(hào)：U455 文獻(xiàn)標(biāo)志碼：A

Abstract： In order to solve the problem of slagging in the slurry shield circulation system， avoid the stagnation of the gravel in the slag discharge port and the insufficient energy supply of the circulation system pump， the accumulation quantity of ballast in the mud pipe and the variation of the hydraulic gradient was studied to ensure mucking of the slurry shield circulation system. The discrete element method was used to construct the numerical model of ballast transportation in the horizontal pipe， and the influence of mud flow， mud specific gravity， mass flow rate of ballast on the accumulation quantity of ballast and the hydraulic gradient of the mud pipeline were studied. The research shows that under the current construction parameters， the ballast which is greater than 20 mm is mainly moved at the bottom of the horizontal straight pipe; the increase of mud flow can significantly reduce the accumulation of ballast in the pipeline system， and can greatly weaken the mud weight and mass flow rate， and the hydraulic gradient of the pipeline decreases first and then increases and the minimum value is 763.41 m3/h; the accumulation quantity of ballast decreases with the increase of the specific gravity of the mud， and the hydraulic gradient of the pipeline increases， and the increase of the ballast mass flow rate increases the accumulation quantity of the ballast in the pipeline and the hydraulic gradient of the pipeline.

Key words： slurry shield; circulation system; straight horizontal pipeline; accumulation quantity; hydraulic gradient; mass flow rate

1 引 言

目前泥水盾構(gòu)被廣泛應(yīng)用于我國(guó)地下空間工程以及隧道施工工程中[1-3]。我國(guó)地緣遼闊，盾構(gòu)機(jī)在地下掘進(jìn)過(guò)程中會(huì)面臨多種復(fù)雜地層，在蘭州穿黃工程中就遇到卵石地層，掘進(jìn)過(guò)程中主要遇到兩方面問(wèn)題：①盾構(gòu)機(jī)開(kāi)挖系統(tǒng)出現(xiàn)刀具過(guò)量磨損甚至刀盤面板異常磨損現(xiàn)象;②盾構(gòu)機(jī)環(huán)流系統(tǒng)排渣困難，在排渣口容易出現(xiàn)卵石滯排現(xiàn)象，石渣粒徑和數(shù)量的增大造成環(huán)流系統(tǒng)泵供能不足的現(xiàn)象。環(huán)流系統(tǒng)作為盾構(gòu)機(jī)的排渣系統(tǒng)，其工作的穩(wěn)定性對(duì)盾構(gòu)機(jī)高效掘進(jìn)和安全施工具有重要意義。

目前國(guó)內(nèi)外學(xué)者針對(duì)管道內(nèi)固體顆粒輸送方面開(kāi)展了大量研究。F.Ravelet等[4]通過(guò)試驗(yàn)得到顆粒尺寸與密度對(duì)顆粒形成固定床或分散流的影響，并對(duì)經(jīng)驗(yàn)公式進(jìn)行了修正。王少勇等[5]通過(guò)試驗(yàn)研究不同水力參數(shù)對(duì)管道壓力損失的影響，確定了不同水力參數(shù)與壓力損失間的函數(shù)關(guān)系。Edelin等[6]從能量的角度研究了聚丙烯顆粒在清水中輸送的最佳濃度。周知進(jìn)等[7]研究顆粒組分特性與輸送速度之間的關(guān)系，得到顆粒體積分?jǐn)?shù)。曹斌等[8]通過(guò)高速攝影儀研究不同水力參數(shù)下粗顆粒水力學(xué)特性。隨著計(jì)算機(jī)的發(fā)展，數(shù)值模擬方法被越來(lái)越廣泛地應(yīng)用到管道輸送研究中，相比試驗(yàn)手段，其具有研究周期短、投入成本低等優(yōu)勢(shì)。夏毅敏等[9-10]研究了石渣形狀系數(shù)、等容粒徑和數(shù)量對(duì)石渣起動(dòng)速度的影響，以及盾構(gòu)機(jī)操作參數(shù)對(duì)環(huán)流系統(tǒng)壓力損失和石渣輸送速度的影響。楊端等[11]研究了泥水盾構(gòu)環(huán)流系統(tǒng)水平直管中渣土輸送特性。Gopaliya等[12]利用兩相流理論構(gòu)建了鐵礦填充料漿和沙漿的固體流動(dòng)模型，分析了鐵礦填充料漿和沙漿在水平管內(nèi)的流動(dòng)狀態(tài)。Jiang等[13]采用歐拉多相流模型研究流速和平均固體體積分?jǐn)?shù)對(duì)氮漿流動(dòng)特性的影響。

上述研究的管道輸送介質(zhì)多為清水，其漿液流變特性與泥水盾構(gòu)環(huán)流系統(tǒng)中泥漿有較大區(qū)別，且固體顆粒多為小粒徑細(xì)砂，因而所得研究成果難以適用于泥水盾構(gòu)泥漿管道的石渣輸送。本研究以蘭州某穿黃工程為例，其施工地層分為人工填土層、第四系全新統(tǒng)沖積地層、第四系下更新統(tǒng)沖積地層和第三系中新統(tǒng)泥巖夾砂巖層，施工地層中卵石含量高達(dá)81.95%，粒徑大于10 mm以及小于40 mm的卵石含量約占64.53%，粒徑大于200 mm的漂石含量較少。筆者基于固-液兩相流運(yùn)動(dòng)理論，考慮泥漿流變特性，建立排漿水平管段的石渣運(yùn)動(dòng)模型，研究了泥漿輸送速度、泥漿相對(duì)密度、石渣質(zhì)量流率對(duì)石渣堆積量、泥漿管道輸送阻力的影響。

2 水平直管石渣輸送模型

2.1 泥漿模型

泥漿流動(dòng)過(guò)程中遵循質(zhì)量和動(dòng)量守恒定律，不考慮泥漿和石渣間的熱傳遞與相變特性，其連續(xù)性方程和動(dòng)量方程如下：

2.2 石渣運(yùn)動(dòng)模型

石渣在輸送過(guò)程中受到石渣重力、石渣-石渣、石渣-管壁之間的接觸力和石渣-泥漿相互作用力，其運(yùn)動(dòng)模型如下：

2.3 幾何模型和邊界條件設(shè)置

該隧道工程石渣粒徑主要在10～40 mm范圍，平均石渣粒徑為20 mm。石渣密度為2 500～2 800 kg/m3，為便于分析，設(shè)定石渣密度為2 700 kg/m3。管道模型如圖2所示，水平管道長(zhǎng)L=5 m，直徑D為300 mm，石渣設(shè)置以一定質(zhì)量流率從管道入口進(jìn)入，質(zhì)量流率即為每秒輸送的石渣質(zhì)量。仿真邊界條件采用速度進(jìn)口、總壓出口邊界條件及壁面無(wú)滑移邊界條件，湍流模型選用Standard k-ε流場(chǎng)模型，采用SIMPLE算法求解速度、壓力耦合方程。

3 結(jié)果及分析

3.1 顆粒流態(tài)與流場(chǎng)分析

由于石渣粒徑較大，因此大部分石渣無(wú)法在泥漿浮力和湍流作用下克服重力以懸浮狀態(tài)進(jìn)行移動(dòng)。石渣以一定速度進(jìn)入，泥漿與石渣發(fā)生動(dòng)能交換，推動(dòng)石渣向前運(yùn)動(dòng)。石渣群開(kāi)始做加速運(yùn)動(dòng)，石渣群的運(yùn)動(dòng)軌跡類似拋物線，總體趨勢(shì)呈現(xiàn)順泥漿流動(dòng)方向與重力方向運(yùn)動(dòng)。石渣顆粒從懸移運(yùn)動(dòng)狀態(tài)變成完全推移運(yùn)動(dòng)狀態(tài)移動(dòng)的距離稱為沉降距離。由圖3可知，沉降距離隨泥漿流量Q的增大而增大，當(dāng)石渣開(kāi)始做推移運(yùn)動(dòng)時(shí)，大量石渣與管壁接觸，受到管壁摩擦阻力影響，石渣群整體運(yùn)動(dòng)速度開(kāi)始降低。由于入口石渣顆粒的質(zhì)量流率一定，且管道足夠長(zhǎng)，顆粒群在達(dá)到沉降距離之前的整體運(yùn)動(dòng)速度是高于達(dá)到沉降距離之后的，因此石渣在管道底部會(huì)發(fā)生一定堆積，且泥漿流量越大，堆積量越少。越靠近管道底部的石渣受到泥漿流動(dòng)的影響越小，造成石渣的速度分布不均勻，如圖4所示。管道底部石渣的運(yùn)動(dòng)速度較低，呈推移運(yùn)動(dòng);堆積層上部區(qū)域石渣運(yùn)動(dòng)速度高，呈翻滾式跳躍運(yùn)動(dòng)。石渣堆積也造成泥漿流速分布不均，堆積層較厚的地方流速高，堆積層較薄的地方流速低，如圖5所示。

3.2 泥漿流量的影響

定義泥漿密度與清水密度的比值為泥漿相對(duì)密度s，圖6為不同泥漿相對(duì)密度下石渣的堆積量隨泥漿流量增大的變化曲線?？芍?dāng)泥漿相對(duì)密度一定時(shí)，石渣堆積量隨泥漿流量的增大呈下降趨勢(shì)。

泥漿流量對(duì)石渣堆積量影響明顯，當(dāng)泥漿流量為508.94 m3/h、漿液相對(duì)密度為1.1時(shí)，石渣堆積量為114.8 kg，石渣體積為管道總體積的12%;排漿管的泥漿流量增大至763.41 m3/h時(shí)，石渣的堆積量急劇減小，系統(tǒng)殘留的石渣質(zhì)量?jī)H為37.06 kg;進(jìn)一步增大漿液流量后，石渣堆積量降幅變得平緩。

由圖7可知，當(dāng)泥漿相對(duì)密度、石渣質(zhì)量流率一定時(shí)，水力坡降隨泥漿流量的增大呈先降低后增大的趨勢(shì)，水力坡降均在泥漿流量為763.41 m3/h時(shí)取得最小值。當(dāng)泥漿流量超過(guò)763.41 m3/h后不同泥漿相對(duì)密度下水力坡降增幅幾乎相等。這是由于流量超過(guò)763.41 m3/h后，漿液流速成為了管道水力坡降的主要影響因素，此時(shí)石渣相對(duì)密度對(duì)管道水力坡降的影響程度會(huì)被削弱，因此在大流量工況下石渣相對(duì)密度的增大對(duì)管道水力坡降的影響變得不明顯。

3.3 泥漿相對(duì)密度的影響

圖8為不同泥漿流量下石渣堆積量隨漿液相對(duì)密度變化的規(guī)律。由圖8可知，當(dāng)泥漿流量為508.94 m3/h時(shí)，管道石渣堆積量隨泥漿相對(duì)密度的增大呈衰減趨勢(shì)。隨著泥漿流量的增大，泥漿相對(duì)密度對(duì)石渣堆積量的影響減弱。

3.4 石渣質(zhì)量流率的影響

不同泥漿流量下石渣堆積量隨石渣質(zhì)量流率變化規(guī)律見(jiàn)圖10。由圖10可知，當(dāng)泥漿流量一定時(shí)，石渣堆積量隨石渣質(zhì)量流率的增大呈增大趨勢(shì)，增長(zhǎng)率隨泥漿流量的增大而減小。當(dāng)流量為508.94 m3/h時(shí)，石渣質(zhì)量流率對(duì)石渣堆積量影響較大，隨著泥漿流量的增大，石渣質(zhì)量流率對(duì)石渣堆積量的影響越來(lái)越小。當(dāng)石渣質(zhì)量流率從13 kg/s增大到26 kg/s，泥漿流量為763.41、1 017.88、1 272.35 m3/h時(shí)對(duì)應(yīng)的石渣堆積量分別增加了11.45、7.14、4.60 kg。

4 結(jié) 語(yǔ)

以蘭州某穿黃河工程為例，基于固-液兩相流理論，考慮泥漿流變特性，建立排漿水平管段卵石運(yùn)動(dòng)模型，通過(guò)有限體積法和離散單元法研究了泥漿流量、泥漿相對(duì)密度、石渣質(zhì)量流率對(duì)石渣堆積量、泥漿管道輸送阻力的影響，得到如下結(jié)論：

（1）當(dāng)管路泥漿流量為508.94～1 272.35 m3/h，泥漿相對(duì)密度為1.10～1.30，質(zhì)量流率為13.36～26.71 kg/s時(shí)，大于20 mm粒徑的石渣以底部推移運(yùn)動(dòng)為主。

（2）泥漿流量的增大能大大降低環(huán)流系統(tǒng)中石渣的堆積量，隨著流量的增大環(huán)流系統(tǒng)水力坡降呈先減小后增大趨勢(shì)，在泥漿流量為763.41 m3/h時(shí)取得最小值。

（3）泥漿相對(duì)密度的增大會(huì)減少環(huán)流系統(tǒng)中石渣的堆積量，當(dāng)漿液流量大于763.41 m3/h時(shí)，該作用會(huì)被明顯削弱。當(dāng)泥漿流量為508.94 m3/h時(shí)，水力坡降隨泥漿相對(duì)密度的增大呈減小趨勢(shì)，而當(dāng)泥漿流量大于763.41 m3/h時(shí)，水力坡降隨泥漿相對(duì)密度的增大呈增長(zhǎng)趨勢(shì)。

（4）漿液流量較低時(shí)，石渣質(zhì)量流率的增大會(huì)明顯增大環(huán)流系統(tǒng)中石渣的堆積量，隨著漿液流量的增大該作用會(huì)被削弱。水力坡降隨著石渣質(zhì)量流率的提升呈近似線性增長(zhǎng)趨勢(shì)。

現(xiàn)有研究中石渣粒徑單一且為球形顆粒，可以進(jìn)一步對(duì)非球形石渣和不同粒徑石渣展開(kāi)研究，另外本研究管道局限于水平直管，可以進(jìn)一步對(duì)復(fù)雜管路結(jié)構(gòu)進(jìn)行探究。

參考文獻(xiàn)：

[1] MIN F， ZHU W， LIN C， et al. Opening the Excavation Chamber of the Large-Diameter Size Slurry Shield： a Case Study in Nanjing Yangtze River Tunnel in China[J].Tunnelling and Underground Space Technology，2015，46：18-27.

[2] LIU C， ZHANG Z， REGUEIRO R A. Pile and Pile Group Response to Tunnelling Using a Large Diameter Slurry Shield-Case Study in Shanghai[J]. Computers and Geotechnics，2014，59：21-43.

[3] LIN C， ZHANG Z， WU S， et al.Key Techniques and Important Issues for Slurry Shield Under-Passing Embankments： a Case Study of Hangzhou Qiantang River Tunnel[J]. Tunnelling and Underground Space Technology，2013，38：306-325.

[4] RAVELET F， BAKIR F， KHELLADI S， et al. Experimental Study of Hydraulic Transport of Large Particles in Horizontal Pipes[J]. Experimental Thermal and Fluid Science （EXP THERM FLUID SCI）， 2013， 45（2）：187-197.

[5] 王少勇，吳愛(ài)祥，尹升華，等.膏體料漿管道輸送壓力損失的影響因素[J].北京科技大學(xué)學(xué)報(bào)，2015，37（1）：7-12.

[6] EDELIN D， CZUJKO P C， CASTELAIN C， et al. Experimental Determination of the Energy Optimum for the Transport of Floating Particles in Pipes[J]. Experimental Thermal & Fluid Science，2015，68：634-643.

[7] 周知進(jìn)，劉愛(ài)軍，夏毅敏，等.顆粒組分特性對(duì)揚(yáng)礦硬管輸送速度的影響[J].中南大學(xué)學(xué)報(bào)（自然科學(xué)版），2011，42（9）：2692-2697.

[8] 曹斌，徐心一，夏建新.水力輸送管道中粗顆粒運(yùn)動(dòng)狀態(tài)變化及其判別[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2016，24（4）：672-679.

[9] 夏毅敏，姚菁，吳遁，等.泥水盾構(gòu)水平直管內(nèi)石渣起動(dòng)速度研究[J].隧道建設(shè)，2018，38（3）：392-398.

[10] 夏毅敏，王洋，吳遁，等.泥水盾構(gòu)環(huán)流系統(tǒng)管道輸送特性[J].中南大學(xué)學(xué)報(bào)：自然科學(xué)版，2017，48（11）：2889.

[11] YANG D， XIA Y， WU D， et al. Numerical Investigation of Pipeline Transport Characteristics of Slurry Shield Under Gravel Stratum[J]. Tunnelling and Underground Space Technology， 2018，71：223-230.

[12] GOPALIYA M， KAUSHAL D. Modeling of Sand-Water Slurry Flow Through Horizontal Pipe Using CFD[J]. Journal of Hydrology and Hydromechanics， 2016， 64（2）：261-272.

[13] JIANG Y Y， ZHANG P. Numerical Investigation of Slush Nitrogen Flow in a Horizontal Pipe[J]. Chemical Engineering Science， 2012， 73：169-180.

【責(zé)任編輯 張華巖】