李素貞, 李 翔

(1.同濟(jì)大學(xué) 建筑工程系,上海200092；2.同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海200092)

?

水平定向穿越鋼管回拖力的分析與實(shí)測(cè)

李素貞1,2, 李 翔1

(1.同濟(jì)大學(xué) 建筑工程系,上海200092；2.同濟(jì)大學(xué) 土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海200092)

為了保證水平定向穿越的順利實(shí)施,回拖力的確定至關(guān)重要.總結(jié)國(guó)內(nèi)外管道回拖力計(jì)算的規(guī)范方法,針對(duì)主要區(qū)別的泥漿阻力及管道彎曲大變形阻力兩部分,建立泥漿阻力模型和管道彎曲大變形阻力模型,形成更精細(xì)的回拖力計(jì)算公式,分析水平定向穿越回拖全過(guò)程的管道受力.對(duì)比某實(shí)際天然氣管線工程整個(gè)回拖過(guò)程的施工監(jiān)測(cè)結(jié)果表明,提出的回拖力計(jì)算方法與實(shí)際情況更吻合,管道彎曲大變形阻力占回拖力的主要部分,回拖力隨管線大轉(zhuǎn)角處拐角的增加呈指數(shù)增長(zhǎng).

水平定向穿越；回拖力；泥漿阻力；管道彎曲大變形；施工監(jiān)測(cè)

城市建設(shè)經(jīng)常伴隨著地下管線的建設(shè)、維修或更換.與傳統(tǒng)的開(kāi)挖施工比較,管道非開(kāi)挖技術(shù)對(duì)地面干擾小、綜合成本低[1].目前大約有41.6%的市政管道修復(fù)工程指定采用非開(kāi)挖技術(shù).水平定向穿越在不到50 a的發(fā)展歷史中,以環(huán)保、效率高的特點(diǎn)在世界各地得到廣泛應(yīng)用,在各種非開(kāi)挖施工方法中所占的比例達(dá)到70.2%.到2013年底,國(guó)內(nèi)水平定向鉆機(jī)的保有量已達(dá)13 030臺(tái),在城市化建設(shè)中發(fā)揮了重要作用[2].

水平定向鉆進(jìn)技術(shù)是以水平定向鉆機(jī)為主安裝管道的一種技術(shù),施工過(guò)程一般可以分為3個(gè)階段,即鉆先導(dǎo)孔、擴(kuò)孔和管道回拖[3-4].管道回拖作為穿越過(guò)程中最關(guān)鍵的一步,由于施工過(guò)程的不可預(yù)見(jiàn)性,與之相關(guān)的回拖力計(jì)算成為首要解決的問(wèn)題.穿越管道在回拖過(guò)程中的受力非常復(fù)雜,目前普遍認(rèn)為主要受到以下5種阻力[3]：1)穿越管道與孔壁間的摩擦阻力；2)管道與地表之間的摩擦阻力；3)絞盤效應(yīng)力,源于沿彎曲鉆孔軌跡拖拉管道產(chǎn)生的遞增承載壓力；4)流體阻力；5)彎曲時(shí)由管道剛度產(chǎn)生的阻力.回拖力的計(jì)算涉及土力學(xué)、工程力學(xué)、流體力學(xué)、彈性力學(xué)等方面的知識(shí),與工程實(shí)際的地質(zhì)環(huán)境、穿越軌跡等密切相關(guān)[5-6].在實(shí)際工程中,由于考慮的孔內(nèi)受力條件不同,回拖力的計(jì)算公式有多種.

我國(guó)《油氣輸送管道穿越工程施工規(guī)范》[7](GB 50424-2007)提出適用于油氣輸送管道的回拖力計(jì)算公式,《給水排水管道工程施工及驗(yàn)收規(guī)范》[8](GB 50268-2008)提出針對(duì)于城鎮(zhèn)公用設(shè)施管道的回拖力計(jì)算公式.美國(guó)材料與實(shí)驗(yàn)協(xié)會(huì)法(American Society for Testing and Materials,ASTM)以絞盤法為基礎(chǔ),提出簡(jiǎn)單、實(shí)用的ASTM法[9]；Huey等[10]為定向鉆安裝鋼管提出回拖力計(jì)算方法,稱為國(guó)際管道研究協(xié)會(huì)(Pipeline Research Council International,PRCI)的計(jì)算方法.一些國(guó)內(nèi)學(xué)者提出卸荷拱土壓力計(jì)算法[5]和凈浮力計(jì)算法[5].Phillps Driscopipe給出計(jì)算聚乙烯管回拖力的計(jì)算公式[11],Polak等[12]對(duì)回拖力各組成部分建立模型,給出相應(yīng)的回拖力計(jì)算公式,在此基礎(chǔ)上,Polak等[13]提出聚乙烯管在HDD施工過(guò)程中的回拖力計(jì)算理論模型以及Cheng等[14]確定了在HDD施工過(guò)程中,管道位于16種幾何形狀下的回拖力計(jì)算理論模型.

比較國(guó)內(nèi)外規(guī)范以及相關(guān)學(xué)者的研究工作發(fā)現(xiàn),各種回拖力計(jì)算方法的主要差異在于泥漿流體阻力和管道彎曲大變形產(chǎn)生的阻力兩部分的分析.針對(duì)回拖管道的實(shí)際受力情況,筆者在前期工作中確立更精細(xì)的泥漿阻力模型和管道彎曲大變形阻力模型[15-16].本文引入這兩類阻力模型建立相應(yīng)的回拖力計(jì)算公式,分析水平定向穿越管道回拖全過(guò)程的受力.結(jié)合某實(shí)際工程水平定向穿越的監(jiān)測(cè)結(jié)果,比較各種規(guī)范方法和本文方法的計(jì)算精度.

1 油氣管道回拖力計(jì)算規(guī)范方法

1.1 《油氣輸送管道穿越工程施工規(guī)范》

《油氣輸送管道穿越工程施工規(guī)范》(GB 5040424-2007)適用于油氣輸送管道在陸上穿越人工或天然障礙的新建或擴(kuò)建工程.回拖力計(jì)算公式主要考慮了穿越管道與孔壁之間的摩擦阻力和泥漿對(duì)管道的黏阻力.

(1)

式中：Ft為計(jì)算的拉力,L2為穿越管道的長(zhǎng)度,f為摩擦系數(shù),D為管子的直徑,ρ1為泥漿的密度,δ1為管子的壁厚,kf為黏滯系數(shù).

1.2 國(guó)際管道研究協(xié)會(huì)(PRCI)的計(jì)算方法

國(guó)際管道研究協(xié)會(huì)的計(jì)算方法是由Huey等[10]在1996年為定向鉆安裝鋼管道而提出的.該方法考慮了管土之間的摩擦阻力、泥漿的黏阻力、重力和彎曲變形阻力的影響,但沒(méi)有考慮管道拖入孔道之前與地表面的摩擦阻力,并認(rèn)為管道進(jìn)入鉆孔時(shí)的回拖力為零,最大回拖力出現(xiàn)在回拖最后階段并以遞增的方法沿管道分布.該方法將整個(gè)管道分解為許多直線段和曲線段,最后的軸向拉力為每小段拉力的總和.

直線段拉力的計(jì)算公式為

(2)

Ff=μbWpLcosβ.

(3)

式中：T1為直線段拉力,Ff為孔道內(nèi)摩擦阻力,FDR為DRAG孔液的阻力，μb為管道與孔壁的摩擦系數(shù),Wp為考慮鉆孔液浮力后管道單位長(zhǎng)度的凈重,L為直線段管道長(zhǎng)度,β為管道傾角.

DRAG孔液的阻力為

FDR=πDLμmud.

(4)

式中：μmud為流體阻力系數(shù),D為管道的直徑.

彎曲段拉力的計(jì)算公式為

(5)

Ff=μbWpLarccosβ.

(6)

式中：T2為彎曲度拉力,Larc為管道弧線長(zhǎng)度. 總的軸向拉力為各段拉力之和:

(7)

式中：Ti為第i段管道上的軸向拉力,Ttot為總的軸向拉力.

2 兩大阻力模型

針對(duì)各種回拖力分析方法中存在較大差異的泥漿流體阻力和管道彎曲大變形阻力,前期工作分別確立了兩大計(jì)算模型[15-16],簡(jiǎn)介如下.

2.1 泥漿阻力模型

將泥漿流動(dòng)模型簡(jiǎn)化成冪律流體在同心環(huán)形空間中的穩(wěn)定流動(dòng),得出管道外表面處泥漿的切應(yīng)力[15]：

(8)

式中：τ為管道外表面泥漿的切應(yīng)力,Κ為稠度系數(shù),v(r)為流體流速,n為流體特性指數(shù),v為管道外泥漿至管道中心的距離,Rp為管道外半徑.

單位長(zhǎng)度上管道外表面上受到的流體阻力為

fd=2πRp·τ.

(9)

式中：fd為單位長(zhǎng)度上管道外表面受到的流體阻力.

2.2 管道彎曲大變形阻力模型

在管道通過(guò)拐角處后,可以將管道在拐角處的受力情況簡(jiǎn)化成圖1所示[16].

圖1 管道拐角處的受力Fig.1 Forces on pipe in corner

AC方向上的力平衡方程為

T2cosψ=(T1+2Rμb)cosψ+2Pμb+NTμb.

(10)

垂直于AC方向上的力平衡方程為

NT=(T1+T2)sinψ.

(11)

由式(10)、(11)可得

(12)

或

ΔT=T1×C1(ψ)+P×C2(ψ).

(13)

式中：P為土體對(duì)管道施加的力,

Φ(p,φ0)=0.847 2+F(p,φ0)-2E(p,φ0),

F(p,φ0) 稱為第一類橢圓積分,E(p,φ0) 稱為第二類橢圓積分.

3 回拖力計(jì)算

將實(shí)際的管道路線進(jìn)行簡(jiǎn)化,得到如圖2所示的模型圖.

圖2 管道路線模型Fig.2 Model for geometry of bore-path

3.1 孔道外管道自重引起的阻力

(14)

或

(15)

式中：Tig為孔道外管道自重引起的阻力,ωp為管道自重,μg為管道與地表面的摩擦系數(shù),L為管道總長(zhǎng),Lk為管道第k段長(zhǎng)度.

3.2 孔道內(nèi)管道自重引起的阻力

(16)

式中：Tis為孔道內(nèi)管道自重引起的阻力,ωb為管道在鉆液中的自重,αk為管道第k段與水平面的夾角.

3.3 泥漿阻力

根據(jù)2.1節(jié)建立的泥漿阻力模型,可以得出管道在泥漿中受到的泥漿阻力：

(17)

式中：Tid為泥漿阻力所導(dǎo)致的拉力,fd為單位長(zhǎng)度上管道外表面上受到的流體阻力.

3.4 拐角處管道彎曲大變形引起的阻力

當(dāng)管道到達(dá)點(diǎn)2 時(shí),T2f由1點(diǎn)處管道的路線方向改變和彎曲剛度組成,

T2f=T12×C1(ψ1)+P1×C2(ψ1).

(18)

式中：T12=(L-L1)×C3.

當(dāng)管道到達(dá)點(diǎn)3時(shí),T3f由1、2點(diǎn)處的管道路線改變和彎曲剛度組成,T3f=T13f+T23f,

T13f=T13×C1(ψ2)+P1×C2(ψ2),

(19)

T23f=T23×C1(ψ2)+P2×C2(ψ2).

(20)

式中：T13=(L-L1-L2)×C3，T23=T3g+T2s+T2d+T13f.

同理,當(dāng)管道到達(dá)4,5,……點(diǎn)時(shí)均可迭代計(jì)算出來(lái).式(18)～(20)中,Tif為管道到達(dá)點(diǎn)i時(shí),點(diǎn)i處管道彎曲大變形引起的阻力；Tji為當(dāng)管道到達(dá)點(diǎn)i時(shí),點(diǎn)j處管道的軸向拉力；Tjif為當(dāng)管道到達(dá)點(diǎn)i時(shí),點(diǎn)j處管道彎曲大變形引起的阻力.

3.5 總的回拖力

在點(diǎn)i處的回拖力為

(21)

式中：Ti為管道到達(dá)點(diǎn)i時(shí),管道回拖所需要的拉力.

4 實(shí)例計(jì)算

結(jié)合某實(shí)際管道工程的水平定向穿越施工監(jiān)測(cè),比較各種回拖力計(jì)算方法的結(jié)果.穿越處地質(zhì)為淤泥質(zhì)粉質(zhì)黏土、淤泥質(zhì)黏土以及流砂土.穿越管道為高頻直縫電阻焊鋼管(SS400)D325x8,最大穿越設(shè)計(jì)深度為21.4 m,管道穿越長(zhǎng)度為520 m,選用的鉆機(jī)型號(hào)為DDW1200型拖車式定向鉆機(jī).在穿越管道前焊接同種材質(zhì)和大小的4 m測(cè)試管道,在測(cè)試管道圓周上下、左右對(duì)稱布置FBG光纖傳感器.如圖3所示,測(cè)得整個(gè)穿越過(guò)程中的應(yīng)變,并在管道內(nèi)布置陀螺儀,實(shí)時(shí)測(cè)試管道在回拖過(guò)程中的位置信息.整個(gè)施工監(jiān)測(cè)方案如圖4所示.測(cè)試過(guò)程與結(jié)果詳見(jiàn)文獻(xiàn)[17].

圖3 光纖傳感器的布置Fig.3 Sensor placement

圖4 HDD施工監(jiān)測(cè)方案[17]Fig.4 Scheme for HDD construction monitoring

圖5 管道截面圖[17]Fig.5 Cross section of pipe

如圖5所示，通過(guò)管道上傳感器布置的位置及應(yīng)變數(shù)據(jù)，得出測(cè)試管道的軸力和彎矩的計(jì)算公式如下：

(22)

和

(23)

式中：N為測(cè)試管道的軸力,M為測(cè)試管道的彎矩,hA、hC分別為點(diǎn)A、C距離中性軸的距離,σA、σC分別為點(diǎn)A、C處的應(yīng)力,R為管道的外半徑,r為管道的內(nèi)半徑,S為管道的橫截面積.

取管道與地表的摩擦系數(shù)為0.3,管道在孔內(nèi)的摩擦系數(shù)為0.25,泥漿密度為1.2 t/m3,鋼管的彈性模量為2.05×105MPa.將實(shí)際的三維管道穿越路線進(jìn)行平面擬合,并將其投影至平面,對(duì)投影點(diǎn)進(jìn)行最小二乘法線段擬合,如圖6所示.同時(shí)得出關(guān)鍵點(diǎn)的拐角,如表1所示.1、7分別為管道路線入口和管道路線與水平面的交點(diǎn),2～6分別為擬合后的管線拐點(diǎn),x為水平方向坐標(biāo)值,y為豎直方向坐標(biāo)值.

將上述關(guān)于管道路線、土體和管材等條件分別代入GB 50424-2007、PRCI方法以及本文方法的回拖力計(jì)算公式中,得到的結(jié)果如圖7所示.圖中,T為回拖力,i為管線的各個(gè)關(guān)鍵點(diǎn),i=1,2,3,4,5,6,7.在管道回拖力的計(jì)算上,PRCI方法和GB 50424-2007在計(jì)算結(jié)果上相差無(wú)幾；本文方法的計(jì)算結(jié)果與PRCI方法和GB 50424-2007在整體趨勢(shì)上保持一致,數(shù)值上存在一個(gè)較穩(wěn)定的差值.實(shí)測(cè)值是保證管道順利安全穿越,鉆機(jī)所保持的回拖力,理論上應(yīng)為設(shè)計(jì)值的最大值乘以安全系數(shù),該工況的實(shí)測(cè)值遠(yuǎn)大于計(jì)算最大值,偏于保守.

圖6 管道路線擬合圖Fig.6 Fitting diagram of pipe route

關(guān)鍵點(diǎn)α/(°)23．4136．8948．1654．8464．80

圖7 各種回拖力計(jì)算公式的結(jié)果對(duì)比Fig.7 Contrast of calculation results of different pullback force formulas

圖8 回拖力各組成部分的比例Fig.8 Proportion of each component of pullback force

圖9 回拖力與拐角的關(guān)系曲線Fig.9 Relation of pullback force and transit angle

改變拐角最大的4點(diǎn)處拐角α,研究管道回拖至點(diǎn)7處時(shí)T的變化情況.如圖9所示,回拖力隨著拐角的變化呈指數(shù)增加.以該穿越工程為例,選用的鉆機(jī)可提供的最大拉力為1 200 kN,理論上以回拖力為限制條件,拐角不應(yīng)超過(guò)20.6°.同理,該方法可以為相似穿越工程的路線設(shè)計(jì)和HDD鉆機(jī)的工作機(jī)制設(shè)計(jì)提供依據(jù).

5 結(jié) 論

(1) 本文提出的回拖力計(jì)算公式充分反映了管道在回拖各個(gè)階段的受力情況,比GB 50424-2007、PRCI的方法接近實(shí)際情況,且與GB 50424-2007、PRCI計(jì)算結(jié)果的穩(wěn)定差值反映出所考慮的泥漿阻力和管道彎曲大變形阻力.

(2) 縱觀管道的整個(gè)回拖過(guò)程可知,在回拖前期階段,影響最大的是管道與地表面的摩擦力；在后期階段中,影響最大的是管道彎曲大變形阻力；泥漿阻力和管道與孔壁的摩擦力相對(duì)影響較小.在整個(gè)回拖過(guò)程中,可以通過(guò)調(diào)整管道拐角來(lái)有效減小施工所需要的回拖力.

(3) 回拖力隨管線大轉(zhuǎn)角處拐角的增加呈指數(shù)增長(zhǎng).在實(shí)際施工中,可以預(yù)先根據(jù)采用的鉆機(jī)型號(hào),利用拐角與回拖力的關(guān)系圖來(lái)判斷穿越路線中的拐角是否合理,從而使得整個(gè)施工過(guò)程能夠安全進(jìn)行.

[1] 王軍,吳君侯,章宇,等.上海地區(qū)深基坑工程的環(huán)境安全問(wèn)題[J].巖土工程學(xué)報(bào), 2007, 28(B11):1328-1331. WANG Jun, WU Jun-hou, ZHANG Yu, et al. Environmental safety of deep excavation engineering in Shanghai [J]. Chinese Journal of Geotechnical Engineering, 2007, 28(B11): 1328-1331.

[2] 顏純文. 2013年非開(kāi)挖行業(yè)統(tǒng)計(jì)分析[C]∥2014年非開(kāi)挖技術(shù)會(huì)議論文集. 廣東：[s.n.], 2014: 1-11. YAN. Chun-wen. Review and analysis of trenchless technology in 2013 [C]∥International Trenchless Technology Conference. Guangdong: [s.n.], 2014: 1-11.

[3] 馬保松.非開(kāi)挖工程學(xué)[M].北京:人民交通出版社, 2008: 532.

[4] 顏純文, STEIN D. 非開(kāi)挖地下管線施工技術(shù)及其應(yīng)用[M].北京:震出版社, 1992: 74.

[5] 安金龍.水平定向鉆穿越回拖力的計(jì)算方法及其分析[J].石油工程建設(shè), 2008, 34(1):21-26. AN Jin-long. Calculation and analysis of pull-back force of pipeline crossing by horizontal directional drilling [J]. Petroleum Engineering Construction, 2008, 34(1):21-26.

[6] 劉大鵬,尤曉偉.土力學(xué)[M].北京:清華大學(xué)出版社,北京交通大學(xué)出版社,2005: 58.

[7] GB 50424-2007,油氣輸送管道穿越工程施工規(guī)范[S].北京：中國(guó)計(jì)劃出版社,2008:17.

[8] GB 50268-2008,給水排水管道工程施工及驗(yàn)收規(guī)范[S].北京：中國(guó)建筑工業(yè)出版社,2008: 99.

[9] ASTM F1962-99, Standard guide for use of maxi-horizontal directional drilling for placement of polyethylene pipe or conduit under obstacles, including rivers crossings[S]. West Conshohocken: ASTM International, 2005.

[10] HUEY D P, HAIR J D, MCLEOD K B. Installation loading and stress analysis involved with pipelines installed by horizontal directional drilling [R]. Chicago： North American Society for Trenchless Technology,1996.

[11] BAUMERT M E, ALLOUCHE E N. Methods for estimating pipe pullback loads for horizontal directional drilling (HDD) crossings [J]. Journal of Infrastructure Systems, 2002, 8(1): 12-19.

[12] POLAK M A, LASHEEN A. Mechanical modelling for pipes in horizontal directional drilling [J]. Tunneling and Underground Space Technology, 2001,16(supple.1): 47-55.

[13] POLAK M A, CHU D. Pulling loads for polyethylene pipes in horizontal directional drilling: theoretical modeling and parametric study [J]. Journal of Infrastructure Systems, 2005, 11(2): 142-150.

[14] CHENG E, POLAK M A. Theoretical model for calculating pulling loads for pipes in horizontal drectional drilling [J]. Tunneling and Underground Space Technology, 2007, 22(5): 633-643.

[15] 蔡亮學(xué),何利民,呂宇玲,等.水平定向鉆管道穿越孔底泥漿的力學(xué)特性[J]. 油氣儲(chǔ)運(yùn), 2011, 30(1): 25-29. CAI Liang-xue, HE Li-min, LV Yu-ling, et al. Hole-bosom slurry mechanical properties of horizontal directional drilling in pipeline crossing project [J]. Oil and Gas Storage and Transportation, 2011, 30(1): 25-29.

[16] 李翔,李素貞.管道穿越大角度借轉(zhuǎn)附加應(yīng)力分析[R].上海：同濟(jì)大學(xué),2015. LI Xiang, LI Su-zhen. The additional stress analysis of pipeline in the big angle corner [R]. Shanghai: Tongji University, 2015.

[17] LI S. Construction monitoring of a municipal gas pipeline during horizontal directional drilling [J]. Journal of Pipeline Systems Engineering and Practice, 2013,4(4): 04013005.

Analysis and measurement on pullback force of steel pipeline during horizontal directional crossing

LI Su-zhen1,2, LI Xiang1

(1.DepartmentofBuildingEngineering,TongjiUniversity,Shanghai200092,China；2.StateKeyLaboratoryofDisasterResearchinCivilEngineering,TongjiUniversity,Shanghai200092,China)

The determination of pullback force on the pipeline during installation is critical to the successful implement of pipeline crossing. The domestic and international standards regarding the calculation of pullback force were summarized and compared. The main difference among these regulations lies in the analysis of the resistances due to mud flow and pipeline large deformation. The new method was proposed based on the evaluation on the accuracy of the two resistance models in order to calculate the pullback force of a steel pipeline during horizontal directional drilling (HDD) construction. The analysis on the mechanical behavior of pipeline was conducted. A practical engineering project on the installation of a gas pipeline was employed. The proposed method accorded with the actual situation. The resistance caused by pipeline large deformation is a dominant factor contributed to the total pullback force. The pullback force increased exponentially with the increase of the transit angle due to pipeline deformation.

horizontal directional crossing; pullback force; mud flow resistance; pipeline large deformation; construction monitoring

2015-10-20. 浙江大學(xué)學(xué)報(bào)(工學(xué)版)網(wǎng)址： www.journals.zju.edu.cn/eng

土木工程防災(zāi)國(guó)家重點(diǎn)實(shí)驗(yàn)室自主研究課題基金資助項(xiàng)目(SLDRCE14-B-19);霍英東教育基金會(huì)高等院校青年教師基金資助項(xiàng)目(142004).

李素貞(1978—),女,副教授，從事生命線工程的研究.E-mail: lszh@#edu.cn

10.3785/j.issn.1008-973X.2016.04.016

TU 990

A

1008-973X(2016)04-0714-06