呼煥苗，曹銀萍※，韋 亮，牛占山，潘 瑩，竇益華，鄭 杰

（1.西安石油大學(xué)機(jī)械工程學(xué)院，西安 710065；2.中國(guó)石油西部鉆探試油公司，新疆克拉瑪依 834000）

0 引言

連續(xù)油管（Coiled－Tubing，CT）是相對(duì)于用螺紋連接下井的常規(guī)油管而言的，是卷繞在滾筒上拉直后直接下井的長(zhǎng)油管[1]。由于連續(xù)油管能夠大幅度節(jié)省作業(yè)時(shí)間和作業(yè)成本，目前廣泛用于油田鉆井、完井、壓裂、酸化、排液、試油、采油采氣、修井等領(lǐng)域[2]。連續(xù)油管作業(yè)時(shí)，在井中起下一回都要經(jīng)歷6次拉直—彎曲塑性變形，極易發(fā)生疲勞失效[3－4]。近年來(lái)，隨著對(duì)一些高溫高壓高深井的勘探，連續(xù)油管的應(yīng)用越來(lái)越廣泛。因此，研究不同內(nèi)壓下連續(xù)油管纏繞/彎曲作用發(fā)生塑性變形進(jìn)而會(huì)發(fā)生低周疲勞具有重大的實(shí)際意義。

目前，連續(xù)油管疲勞壽命通常都是采用疲勞壽命理論模型、室內(nèi)疲勞實(shí)驗(yàn)以及現(xiàn)場(chǎng)試驗(yàn)等方法。2001年，王優(yōu)強(qiáng)等[5－6]建立了半經(jīng)驗(yàn)公式的連續(xù)油管壽命預(yù)測(cè)模型。2004年，朱小平[7]以塑性力學(xué)理論為基礎(chǔ)，進(jìn)行了連續(xù)油管在循環(huán)彎曲作用下疲勞壽命的研究。2006年，Tipton等[8]根據(jù)Miner理論建立了連續(xù)油管疲勞壽命預(yù)測(cè)模型。2008年，王海濤[9]首次結(jié)合連續(xù)油管的應(yīng)力狀態(tài)，將三參數(shù)冪函數(shù)能量法用于疲勞壽命預(yù)測(cè)。2012年，李子豐等[10]通過(guò)對(duì)預(yù)彎曲連續(xù)油管的受力分析和強(qiáng)度校核，建立了預(yù)彎曲連續(xù)油管疲勞壽命預(yù)測(cè)模型。2017年，賀海軍[11]等用ANSYS Workbench有限元分析軟件建立了滾筒處的連續(xù)油管三維有限元力學(xué)分析模型。2017年，竇益華[12]等通過(guò)分別考察不同松弛角度和鋼級(jí)的連續(xù)油管彈塑性極限彎矩以及彈性極限彎曲半徑等力學(xué)性能，得出當(dāng)松弛角度在45°～50°時(shí)，連續(xù)油管會(huì)獲得較好的力學(xué)性能。2018年，于桂杰[13]等利用人工神經(jīng)網(wǎng)絡(luò)對(duì)有表面缺陷的連續(xù)油管進(jìn)行疲勞壽命預(yù)測(cè)。Tian等[14]提出一種改進(jìn)的擴(kuò)展有限元（IXFEM）方法，從根本上克服了標(biāo)準(zhǔn)XFEM和修正XFEM的線性依賴和病態(tài)問(wèn)題。上述文獻(xiàn)雖然對(duì)連續(xù)油管疲勞壽命預(yù)測(cè)模型進(jìn)行了修正，但是沒(méi)有涉及到不同工況，因而有必要研究連續(xù)油管在內(nèi)壓作用下的疲勞壽命。

在實(shí)際應(yīng)用中，連續(xù)油管帶內(nèi)壓工作時(shí)的失效故障時(shí)常發(fā)生，尤其是壓裂工況。本文基于ABAQUS/FE－SAFE軟件，對(duì)不同內(nèi)壓下纏繞/彎曲連續(xù)油管進(jìn)行了疲勞壽命分析。

1 纏繞/彎曲連續(xù)油管有限元模型的建立

考慮到多次纏繞/彎曲過(guò)程中施加內(nèi)壓的情況下，連續(xù)油管的疲勞壽命會(huì)受到影響，所以進(jìn)行了纏繞/彎曲過(guò)程中的仿真分析。以QT－900作為研究對(duì)象，基于美國(guó)Stewart&Stevenson公司的連續(xù)油管疲勞試驗(yàn)機(jī)工作原理進(jìn)行仿真分析[15]。連續(xù)油管纏繞/彎曲有限元模型如圖1所示。圖1（a）為連續(xù)油管纏繞在滾筒上的尺寸以及連續(xù)油管、滾筒、校直模的位置，其中連續(xù)油管長(zhǎng)度為800 mm，滾筒的半徑為1 375 mm，矯直模的長(zhǎng)度為550 mm。根據(jù)上述尺寸創(chuàng)建三維有限元模型，其中，有連續(xù)油管、滾筒、滾輪和擋板，為了簡(jiǎn)化模型，將滾筒用彎曲板代替，矯直模用擋板代替，彎曲拉直用滾輪驅(qū)動(dòng)來(lái)實(shí)現(xiàn)，具體的有限元模型如圖1（b）所示。

圖1 連續(xù)油管纏繞/彎曲有限元模型Fig.1 Finite element model of coiled tubing winding/bending

2 纏繞/彎曲連續(xù)油管有限元分析

2.1 材料屬性設(shè)置

使用ABAQUS對(duì)連續(xù)油管的彎曲—拉直力學(xué)行為進(jìn)行有限元模擬分析，以QT－900連續(xù)油管進(jìn)行疲勞壽命分析，其材料屬性如表1所示。

表1 QT－900連續(xù)油管材料參數(shù)Tab.1 QT－900 coiled tubing material parameters

2.2 網(wǎng)格劃分

為了簡(jiǎn)化模型，將滾筒和矯直模使用殼單元進(jìn)行創(chuàng)建并設(shè)置為剛體。給連續(xù)油管的有限元體單元模型進(jìn)行材料屬性的設(shè)置。采用三維八節(jié)點(diǎn)減縮積分單元即C3D8R單元進(jìn)行網(wǎng)格劃分，單元選擇實(shí)體單元；減縮積分單元在彎曲載荷下也不易發(fā)生單元過(guò)于剛硬的問(wèn)題，還能顯著提高計(jì)算效率。同時(shí)為了提高計(jì)算的精度，對(duì)連續(xù)油管的厚度方向設(shè)置兩層六面體網(wǎng)格，其他部件均設(shè)置為剛體，網(wǎng)格類型為4結(jié)點(diǎn)的三維雙線性剛體四邊形即R3D4。

2.3 施加邊界條件和計(jì)算

采用靜態(tài)通用算法進(jìn)行有限元計(jì)算，設(shè)置兩個(gè)分析步，一是彎曲分析步，二是矯直分析步，設(shè)置的時(shí)間均為0.1 s。模型中的接觸通過(guò)相互接觸的面之間定義接觸對(duì)來(lái)實(shí)現(xiàn)，設(shè)置連續(xù)油管與滾筒面—面接觸對(duì)，連續(xù)油管與矯直模面—面接觸對(duì)以及連續(xù)油管與滾輪之間面—面接觸對(duì)，設(shè)置滾筒、矯直模和滾輪為主面，連續(xù)油管為從面。接觸屬性設(shè)置法向和切向行為，法向設(shè)置“硬”接觸，這樣可以防止連續(xù)油管穿過(guò)剛體模具；切向設(shè)置摩擦因數(shù)，模擬設(shè)置因數(shù)為0.1。這里對(duì)連續(xù)油管未設(shè)置載荷而是設(shè)置邊界條件。圖2所示為邊界條件的設(shè)定，將連續(xù)油管的一端、滾筒和矯直模固定，給滾輪設(shè)置位移，彎曲分析的時(shí)候給上滾輪設(shè)置向下的位移，矯直分析的時(shí)候給下邊的滾輪設(shè)置向上的位移，將兩個(gè)滾輪分別耦合到自己的參考點(diǎn)上。另外，給滾輪設(shè)置了轉(zhuǎn)動(dòng)慣量，一個(gè)剛體，如果在所有自由度上有一個(gè)任意自由度方向沒(méi)有被約束住的話，那么在受到載荷時(shí)，在這個(gè)自由度上必然出現(xiàn)無(wú)窮大的位移。即使受到一個(gè)很小的力，也會(huì)在這個(gè)自由度上產(chǎn)生無(wú)窮大的加速度，即產(chǎn)生無(wú)窮大的位移，所以對(duì)剛體通過(guò)設(shè)置轉(zhuǎn)動(dòng)慣量來(lái)設(shè)置質(zhì)量。

圖2 邊界條件設(shè)置Fig.2 Boundary condition setting

2.4 不同內(nèi)壓下連續(xù)油管仿真計(jì)算結(jié)果分析

為了考察不同內(nèi)壓下連續(xù)油管的疲勞壽命情況，對(duì)連續(xù)油管的管內(nèi)設(shè)置不同的壓力進(jìn)行三維有限元建模與計(jì)算，在前述力學(xué)模型計(jì)算得到連續(xù)油管彎曲和矯直的基礎(chǔ)上設(shè)置位移值，設(shè)置向上和向下的位移均為159 mm。在同樣的工況下，分別設(shè)置0～80 MPa的內(nèi)壓，分別計(jì)算得到不同內(nèi)壓作用下連續(xù)油管的彎曲—矯直應(yīng)力應(yīng)變結(jié)果的ODB文件，圖3所示為內(nèi)壓為20 MPa時(shí)的彎曲—矯直云圖。

圖3 內(nèi)壓20 MPa連續(xù)油管有限元計(jì)算結(jié)果Fig.3 Finite element calculation results of coiled tubing with internal pressure of 20 MPa

通過(guò)仿真計(jì)算結(jié)果分析可知，連續(xù)油管在滾筒上纏繞/彎曲是一個(gè)復(fù)雜的塑形變形過(guò)程。連續(xù)油管完全纏繞/彎曲狀態(tài)下，除彎曲中性層附近的部分，大部分截面已經(jīng)產(chǎn)生了塑性變形。因此，矯直狀態(tài)下會(huì)有明顯的殘余應(yīng)力，殘余應(yīng)力較彎曲狀態(tài)減小。

3 連續(xù)油管低周疲勞壽命預(yù)測(cè)

3.1 Fe－safe疲勞壽命計(jì)算

在有限元分析結(jié)果的基礎(chǔ)上，利用Fe－safe軟件進(jìn)行連續(xù)油管的疲勞壽命分析。將有限元計(jì)算的ODB結(jié)果導(dǎo)入Fe－safe，選取兩個(gè)分析步當(dāng)中最后一幀進(jìn)行計(jì)算，計(jì)算過(guò)程中單位設(shè)置一致，材料參數(shù)設(shè)置和前述有限元模擬當(dāng)中設(shè)置也一致，然后選取有限元模型中連續(xù)油管集合進(jìn)行分析計(jì)算。設(shè)置表面粗糙度為0.6＜Ra＜1.6，最后選取粗糙度算法，選取默認(rèn)的Normal Strain：－Morrow算法進(jìn)行計(jì)算。

3.2 連續(xù)油管疲勞壽命預(yù)測(cè)分析

為了考慮壓裂工況下的連續(xù)油管疲勞壽命，將不同內(nèi)壓下的有限元計(jì)算結(jié)果按照上述計(jì)算過(guò)程進(jìn)行疲勞壽命的計(jì)算，得到0～80 MPa壓力下的疲勞壽命結(jié)果，圖4所示分別為20 MPa和70 MPa內(nèi)壓下的計(jì)算結(jié)果云圖。

圖4 20 MPa和70 MPa內(nèi)壓下連續(xù)油管疲勞壽命云圖Fig.4 Fatigue life cloud of coiled tubing at 20 MPa and 70 MPa internal pressure

為了直觀地得到連續(xù)油管疲勞壽命隨內(nèi)壓變化的趨勢(shì)，依據(jù)仿真結(jié)果繪制了不同內(nèi)壓下連續(xù)油管疲勞壽命曲線圖形，如圖5所示?？梢灾庇^地看到，隨著內(nèi)壓的不斷增大，連續(xù)油管循環(huán)次數(shù)呈現(xiàn)非線性下降的趨勢(shì)，沒(méi)有內(nèi)壓的時(shí)候約1 000次，當(dāng)內(nèi)壓增加到60 MPa時(shí)，循環(huán)次數(shù)約為96次，而當(dāng)內(nèi)壓為80 MPa時(shí)，循環(huán)次數(shù)為50次左右。

圖5 不同內(nèi)壓下連續(xù)油管疲勞壽命曲線Fig.5 Fatigue life curves of coiled tubing under different internal pressures

4 結(jié)束語(yǔ)

考慮到連續(xù)油管彎折作用下會(huì)有殘余應(yīng)力，在壓裂過(guò)程中多次彎折作用下會(huì)造成塑性變形，進(jìn)而影響其疲勞壽命。所以采用有限元方法進(jìn)行了壓裂工況下連續(xù)油管纏繞/彎曲低周疲勞壽命分析，得到如下結(jié)論。

（1）通過(guò)對(duì)不同內(nèi)壓下連續(xù)油管纏繞/彎曲有限元仿真分析表明，完全纏繞在滾筒上的連續(xù)油管截面大部分進(jìn)入了塑性區(qū)，只有中性軸附近的部分還處于彈性階段。

（2）當(dāng)內(nèi)壓增大到10 MPa，循環(huán)次數(shù)降到855次；當(dāng)內(nèi)壓增大到80 MPa，循環(huán)次數(shù)降到50次，與現(xiàn)有的實(shí)驗(yàn)結(jié)果吻合較好，所以有限元仿真方法提供了一種便捷的疲勞壽命分析的思路。

（3）在壓裂工況下，隨著內(nèi)壓的增大，連續(xù)油管的疲勞壽命是呈現(xiàn)非線性下降的趨勢(shì)，為了保證壓裂安全，應(yīng)該合理控制壓裂內(nèi)壓。為連續(xù)油管的現(xiàn)場(chǎng)操作提供安全保障，具有重要實(shí)際意義。