周 蘭 孫巧雷 王高磊 馬 越 馮 定
(長江大學(xué)機(jī)械工程學(xué)院;湖北省油氣鉆完井工具工程技術(shù)研究中心)
目前,由于常用的取換套、爆炸整形及套管補(bǔ)貼等修復(fù)技術(shù)已不適用于水平井多段破損或破損段較長的修復(fù)作業(yè),所以油田現(xiàn)場大多采用膨脹管補(bǔ)貼進(jìn)行修復(fù)[1-2]。膨脹管在水平井下入過程中,受井眼軌跡的影響會發(fā)生彎曲,且膨脹管螺紋連接處可能擠壓受損,導(dǎo)致強(qiáng)度降低甚至發(fā)生脫扣現(xiàn)象[3-5]。
針對水平井的管柱下入及通過性問題,部分學(xué)者開展了相關(guān)研究工作。馮定等[6]基于結(jié)構(gòu)力學(xué)和位移法,在充分考慮井筒摩阻、管柱自重及井眼軌跡等多種因素的影響下,推導(dǎo)出了多分層管柱整體通過性的研究方法;祝效華等[7]基于彈塑性理論,利用Abaqus非線性分析軟件,研究了套管通過水平井彎曲段時(shí)的力學(xué)行為,并總結(jié)出套管通過水平井彎曲段時(shí)的下入規(guī)律;羅有剛等[8]在考慮管柱組合結(jié)構(gòu)、流體摩阻及管柱存在多變徑部位等影響因素下,建立了沖砂管柱下入摩阻預(yù)測模型;代理震[9]針對水平井壓裂管柱封隔器的下入安全進(jìn)行了預(yù)測,并推導(dǎo)出可通過造斜段工具最大長度的計(jì)算公式。但是,目前尚無學(xué)者針對井眼彎曲現(xiàn)象嚴(yán)重及井眼曲率大的水平井開展膨脹管的通過性影響因素研究。
為此,本文基于膨脹管在水平井中的下入過程,在考慮膨脹管所受重力、浮力、摩阻、彎矩、剪力及流體作用力等的情況下[10-12],對膨脹管微元段進(jìn)行受力分析,建立膨脹管下入載荷計(jì)算模型以及井眼曲率影響下產(chǎn)生的附加軸力計(jì)算模型;在考慮膨脹管螺紋不同螺紋牙中徑和軸向載荷等不同情況下,建立膨脹管內(nèi)、外螺紋連接強(qiáng)度條件。此外,運(yùn)用所建立的模型分析了CA水平井井眼曲率對膨脹管通過性的影響,包括膨脹管通過性分析、井眼曲率對膨脹管通過性的影響分析以及極限曲率的計(jì)算與分析。所得結(jié)果可為大井眼曲率水平井下入膨脹管困難提供有效解決方案和理論依據(jù)。
基于水平井井眼軌跡,選取膨脹管為研究對象。受力分析前,將膨脹管分為若干個(gè)微元段,在綜合考慮膨脹管所受重力、軸向載荷、浮力、摩阻、彎矩、剪力和流體作用力等影響下,對膨脹管微元段進(jìn)行受力分析并建立受力載荷模型[13-18]。相比于其他下入受力模型,該模型適用于井眼彎曲現(xiàn)象嚴(yán)重和井眼曲率大的水平井,且考慮了膨脹管的剛性,更符合現(xiàn)場實(shí)際工況。建立模型前,首先做如下假設(shè):①膨脹管軸線與井眼軸線一致;②膨脹管與井壁連續(xù)性接觸,且井壁為剛性;③膨脹管微元段受力均勻;④忽略溫度效應(yīng)對膨脹管的影響;⑤膨脹管微元段的曲率為常數(shù);⑥忽略膨脹管下入過程中由于變形效應(yīng)所導(dǎo)致的內(nèi)外徑變化;⑦忽略膨脹管下入時(shí)扭矩的影響。
圖1 膨脹管微元段受力分析圖
根據(jù)受力平衡可得膨脹管微元段各作用力的關(guān)系式:
(1)
膨脹管微元段的均布接觸力為:
(2)
式中:μ為膨脹管與井壁間的摩擦因數(shù),“±”分別表示膨脹管上提、下放兩種工況;N為膨脹管微元段接觸力的合力,N/m;Nn、Nb為膨脹管微元段主法線、副法線方向所受的均布接觸力,N/m。
膨脹管單位長度浮重為:
(3)
膨脹管微元段A處所受集中力為:
(4)
式中:T(s)為膨脹管微元段A處軸向力,N;Fn(s)、Fb(s)分別為膨脹管微元段A處的主法線、副法線方向所受的剪力,N。
膨脹管微元段B處所受集中力為:
(5)
式中:dT為膨脹管微元段的軸向力增量,N;dFn、dFb分別為膨脹管微元段主、副法線上的剪力增量,N。
(6)
當(dāng)膨脹管下入到水平井的彎曲段時(shí),受井眼曲率的影響其會產(chǎn)生附加軸力。若曲率過大,則會對膨脹管產(chǎn)生一定的破壞,且膨脹管螺紋連接處可能擠壓受損,導(dǎo)致強(qiáng)度降低甚至脫扣。因此,基于材料力學(xué)的壓桿失穩(wěn)理論,建立了膨脹管的撓曲線近似微分方程,通過求解可得膨脹管曲率附加軸力的計(jì)算公式。
膨脹管撓曲線近似微分方程為:
(7)
式中:M(x)為膨脹管x處截面彎矩,N·m;Fw為膨脹管彎曲臨界力,即曲率附加軸力,N;δ為膨脹管自由端撓度,m;ω為膨脹管x處的撓度,m。
對上式簡化并求解,可得該微分方程通解為:
ω=Asin(kx)+Bcos(kx)+δ
(8)
式中:A、B、k為待定常數(shù)。
邊界條件:x=0處,ω=0;x=l處,ω=δ。
通過對通解求導(dǎo),并將邊界條件代入求解,可得膨脹管曲率附加軸力:
(9)
由于膨脹管螺紋不同位置的中徑、軸向載荷不同,導(dǎo)致螺紋牙發(fā)生破壞的危險(xiǎn)截面具有不確定性,且由于其連接失效受到軸向載荷、螺紋結(jié)構(gòu)、螺紋尺寸、螺紋各截面處復(fù)合應(yīng)力及腐蝕環(huán)境等多個(gè)因素的綜合影響,需結(jié)合相應(yīng)位置的螺紋中徑參數(shù)和軸向載荷來分別進(jìn)行螺紋牙的失效條件判定[19-20],取其最大工作應(yīng)力。因此,考慮膨脹管外螺紋存在拉應(yīng)力、內(nèi)螺紋存在壓應(yīng)力以及螺紋面上的剪應(yīng)力,基于第四強(qiáng)度理論,可分別得膨脹管內(nèi)、外螺紋連接應(yīng)力[21]。
膨脹管內(nèi)螺紋連接應(yīng)力:
(10)
膨脹管外螺紋連接應(yīng)力:
(11)
式中:σni、σwi分別為膨脹管內(nèi)、外螺紋連接應(yīng)力,MPa;[σn]、[σw]分別為膨脹管內(nèi)、外螺紋許用屈服應(yīng)力,MPa;Ti-1為第i-1圈膨脹管螺紋的軸向載荷,N;Ani、Awi分別為內(nèi)、外螺紋連接處的螺旋副接觸面積,m2;T1、T2分別為螺旋副與接觸面間的摩擦力矩,N·m;T1j為第j個(gè)螺紋牙處的摩擦力矩,N·m;Wni、Wwi分別為內(nèi)、外螺紋連接處螺旋副的接觸面抗扭截面系數(shù)。
基于以上分析,當(dāng)膨脹管下入到彎曲段時(shí),在曲率不同的情況下會產(chǎn)生不同的曲率附加軸力,從而導(dǎo)致膨脹管軸力發(fā)生變化,最終使膨脹管內(nèi)、外螺紋連接應(yīng)力變化。當(dāng)曲率越大時(shí),曲率半徑則越小,膨脹管可通過的最大長度也隨之減小,曲率附加軸力則越大,從而導(dǎo)致內(nèi)、外螺紋連接應(yīng)力隨之增大。因此,當(dāng)水平井井眼曲率增大到一定值時(shí),膨脹管內(nèi)、外螺紋連接應(yīng)力可能大于螺紋連接部位的許用屈服應(yīng)力,從而產(chǎn)生破壞導(dǎo)致失效。
為減輕膨脹管下入時(shí)井眼曲率的影響,需分析膨脹管下入到不同井深下的極限曲率。根據(jù)膨脹管螺紋連接處的許用屈服強(qiáng)度,可確定其極限連接應(yīng)力,從而可分別確定膨脹管軸力、曲率附加軸力以及極限曲率,以此可計(jì)算出膨脹管下入到各井斜角處所對應(yīng)的極限曲率。
本文以CA水平井為例,分析水平井井眼曲率對膨脹管通過性的影響。下入膨脹管管柱組合:可撈底堵×0.50 m+膨脹管×18.00 m+?105 mm膨脹錐×0.18 m +?90 mm膨脹錐連接桿×0.58 m+?73 mm工具油管。其中:膨脹前膨脹管本體外徑為108 mm,內(nèi)徑為94 mm,長度為6 m;螺紋連接處外徑為110 mm,內(nèi)徑為94 mm。通過計(jì)算并對比兩段螺紋的軸向力及內(nèi)、外螺紋連接應(yīng)力,靠近井底螺紋的第二牙受力較大,因此以下分析選取靠近井底螺紋的第二牙為研究對象,以確保兩段膨脹管螺紋的連接應(yīng)力均符合強(qiáng)度要求。
基于以上參數(shù),分析了膨脹管在CA水平井中下入的通過性、井眼曲率對膨脹管通過性的影響以及極限曲率,主要包括對膨脹管下入時(shí)的實(shí)際鉤載、最小通過軸力(下入軸力+曲率附加軸力)、內(nèi)外螺紋連接應(yīng)力及極限曲率的分析計(jì)算,主要計(jì)算步驟及分析結(jié)果如下。
為分析膨脹管在CA水平井中下入時(shí)的通過性,使用式(6)計(jì)算了膨脹管下入到不同井斜角位置時(shí)的實(shí)際鉤載,并使用式(6)和式(9)計(jì)算可得下入軸力和曲率附加軸力,將二者相加,可得膨脹管下入到各井斜角位置處,曲率為每30 m井段1°~10°范圍內(nèi)的最小通過軸力,最終選取各組數(shù)據(jù)中最小通過軸力的最大值與相對應(yīng)井斜角下入位置的實(shí)際鉤載進(jìn)行對比,結(jié)果如圖2所示。由圖2可知,膨脹管下入到各井斜角處的實(shí)際鉤載均大于最小通過軸力的最大值。因此,從鉤載角度來看,膨脹管可下入到CA水平井中各井斜角位置。
圖2 不同井斜角處最小通過軸力的最大值與實(shí)際鉤載對比曲線
為分析水平井井眼曲率對膨脹管通過性的影響,分別計(jì)算了膨脹管下入到各井斜角α位置時(shí),不同曲率下的最小通過軸力和內(nèi)、外螺紋連接應(yīng)力,計(jì)算結(jié)果如圖3和圖4所示。其中:使用式(6)和式(9)分別計(jì)算下入軸力和曲率附加軸力,二者之和即為最小通過軸力;通過式(6)計(jì)算不同曲率下的軸向載荷,將各軸向載荷代入式(10)和式(11)可分別得不同曲率下的內(nèi)、外螺紋連接應(yīng)力。
圖3分析了井斜角分別為15°、35°、55°和75°處不同曲率下的最小通過軸力。由圖3可以看出:最小通過軸力隨曲率的增大而增大;在同一曲率下,最小通過軸力的變化與井斜角的大小并無直接相關(guān)關(guān)系。這是因?yàn)榫矍试酱?,膨脹管的可通過最大長度越短,此時(shí)膨脹管的下入則需要更大的通過軸力,而水平井井眼曲率的變化并不受井斜角變化的影響,因此井斜角的變化并不會導(dǎo)致最小通過軸力發(fā)生變化。
圖3 不同曲率下的最小通過軸力曲線
圖4分析了井斜角分別為15°、35°、55°和75°處不同曲率下的膨脹管內(nèi)、外螺紋連接應(yīng)力。由圖4可以看出:膨脹管內(nèi)、外螺紋連接應(yīng)力均隨曲率的增大而增大,且同一井斜角位置處相同曲率下的外螺紋連接應(yīng)力均大于內(nèi)螺紋連接應(yīng)力;同一曲率下,內(nèi)、外螺紋連接應(yīng)力的大小并不隨井斜角的變化而變化。這是因?yàn)榍实脑龃髸?dǎo)致曲率附加軸力和最小通過軸力增大,從而使內(nèi)、外螺紋連接應(yīng)力增大。由此分析結(jié)果可知,判斷膨脹管在下入過程中是否破壞應(yīng)主要考慮井眼曲率的影響,與下入位置的井斜角大小并無關(guān)系。
圖4 不同曲率下的內(nèi)、外螺紋連接應(yīng)力曲線
圖5展示了膨脹管下入到不同井斜角位置時(shí)的極限曲率,并與膨脹管下入至CA水平井相對應(yīng)井斜角位置下的井眼曲率進(jìn)行對比。其中,膨脹管屈服強(qiáng)度取440 MPa,安全系數(shù)取1.2,由此可得許用屈服強(qiáng)度,將此作為極限連接應(yīng)力,從而通過式(10)、式(11)、式(6)以及式(9)可分別確定膨脹管軸力、曲率附加軸力以及極限曲率。
圖5 不同井斜角處的井眼曲率與極限曲率的對比曲線
計(jì)算結(jié)果表明,CA水平井各井斜角位置的井眼曲率中,存在大于極限曲率的情況。因此,膨脹管在下入過程中會受到破壞,而在CA水平井的現(xiàn)場施工過程中,起出膨脹管的脫扣及起出膨脹錐連接桿的彎曲現(xiàn)象則證實(shí)了這一結(jié)論的準(zhǔn)確性。
對結(jié)果進(jìn)一步分析可知,在計(jì)算的各井斜角位置中,井斜角為20°時(shí)的極限曲率大于井眼曲率,說明膨脹管下入到此位置時(shí),其強(qiáng)度不滿足膨脹管材料的屈服強(qiáng)度,從而會導(dǎo)致膨脹管發(fā)生破壞。
(1)針對井眼彎曲現(xiàn)象嚴(yán)重及井眼曲率大的水平井,綜合考慮膨脹管下入過程中所受重力、浮力、彎矩、剪力及流體作用力等的影響,建立了膨脹管在水平井中下入時(shí)的受力載荷模型,可用于計(jì)算膨脹管下入的軸力和鉤載,且建立了曲率附加軸力的計(jì)算模型,膨脹管下入軸力和曲率附加軸力之和即為最小通過軸力,將最小通過軸力和鉤載進(jìn)行對比即可初步判斷膨脹管的通過性。
(2)提出了計(jì)算極限曲率的方法,通過此方法可得水平井各位置的極限曲率,將其與井眼曲率對比可判斷膨脹管在下入過程中是否會受到破壞,且可初步判斷膨脹管強(qiáng)度破壞的下入位置。
(3)通過分析曲率對膨脹管通過性的影響,膨脹管內(nèi)、外螺紋連接應(yīng)力隨曲率的增大而增大,且同一井斜角位置處相同曲率下的外螺紋連接應(yīng)力均大于內(nèi)螺紋連接應(yīng)力;從鉤載和最小通過軸力最大值的對比來看,膨脹管可下入至CA水平井中的各井斜角位置處,但從井眼曲率和極限曲率的對比來看,膨脹管在下入過程中會受到破壞,分析結(jié)果與CA水平井的現(xiàn)場施工情況相符。