狄勤豐,王明杰,胡以寶,趙域棟,朱衛(wèi)平,王文昌
(1.上海大學(xué)上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所,上海 200072; 2.上海大學(xué)上海市力學(xué)在能源工程中的應(yīng)用重點(diǎn)實(shí)驗室,上海 200072)
柔性短節(jié)位置對帶旋轉(zhuǎn)導(dǎo)向工具底部鉆具組合動力學(xué)特性的影響
狄勤豐1,2,王明杰1,2,胡以寶1,2,趙域棟1,2,朱衛(wèi)平1,2,王文昌1,2
(1.上海大學(xué)上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所,上海 200072; 2.上海大學(xué)上海市力學(xué)在能源工程中的應(yīng)用重點(diǎn)實(shí)驗室,上海 200072)
針對具有超長細(xì)比特征的底部鉆具組合(BHA),運(yùn)用有限元結(jié)點(diǎn)迭代方法,研究柔性短節(jié)在不同位置對帶位移工作方式旋轉(zhuǎn)導(dǎo)向工具RSBHA的三維動態(tài)變形和動態(tài)應(yīng)力的影響,并且分析在不同鉆壓和轉(zhuǎn)速情況下柔性短節(jié)位置的變化對最大動態(tài)應(yīng)力的影響。結(jié)果表明:在合理的位置加入柔性短節(jié)能夠減小RSBHA的橫向動態(tài)位移,并且使得最大動態(tài)應(yīng)力明顯下降,起到減振的效果,反之將會明顯增加動態(tài)應(yīng)力;柔性短節(jié)位置對最大動態(tài)應(yīng)力的影響遠(yuǎn)大于鉆壓和轉(zhuǎn)速,合理地確定柔性短節(jié)位置對提高井下工具的安全性和測量數(shù)據(jù)的準(zhǔn)確性有重要作用。
底部鉆具組合;柔性短節(jié);動力學(xué)特性;旋轉(zhuǎn)導(dǎo)向工具
旋轉(zhuǎn)導(dǎo)向鉆井系統(tǒng)代表了當(dāng)前井眼軌跡控制技術(shù)的最高水平,旋轉(zhuǎn)導(dǎo)向工具的出現(xiàn)使得復(fù)雜結(jié)構(gòu)井的井眼軌跡的高質(zhì)量控制成為可能。旋轉(zhuǎn)導(dǎo)向鉆井系統(tǒng)以Baker Hughes的AutoTrak系統(tǒng)、Schlumberger公司的PowerDrive系統(tǒng)、Sperry Sun公司的GeoPilot系統(tǒng)最具代表性[1-5]。帶旋轉(zhuǎn)導(dǎo)向工具的底部鉆具組合(bottom hole assembly,BHA)中帶有柔性短節(jié)(flex sub)。對鉆柱動力學(xué)特性的研究主要集中在整體鉆柱的動力學(xué)特性及安全性[6-15],也有部分文章利用轉(zhuǎn)子動力學(xué)理論研究BHA的動力學(xué)特性[16],但對帶柔性短節(jié)和旋轉(zhuǎn)導(dǎo)向工具的BHA的動力學(xué)特性分析的文獻(xiàn)很少。為了保證井下工具的安全和隨鉆測量數(shù)據(jù)的準(zhǔn)確,研究柔性短節(jié)對帶旋轉(zhuǎn)導(dǎo)向工具的BHA的動力學(xué)特性的影響十分必要。柔性短節(jié)的材料一般為鈦合金材料(彈性模量為109 GPa),長度均固定在1.2 m左右,筆者運(yùn)用有限元迭代法研究柔性短節(jié)位置對帶位移工作方式旋轉(zhuǎn)導(dǎo)向工具的BHA(這里簡稱為RSBHA)動力學(xué)特性的影響規(guī)律,為優(yōu)化RSBHA結(jié)構(gòu)、降低其動態(tài)應(yīng)力,提高井下工具安全性和測量數(shù)據(jù)的準(zhǔn)確性提供依據(jù)。
假設(shè):井眼截面為圓型;下部鉆具組合為三維彈性梁;忽略鉆柱接頭影響;鉆柱的變形為小變形。
鉆柱動力學(xué)系統(tǒng)的動力學(xué)平衡方程可表示為
為了更好地說明每個自由度的慣性質(zhì)量,將質(zhì)量矩陣M分為兩個部分:M=M1+M2。質(zhì)量矩陣M1包括3個平動及繞軸x轉(zhuǎn)動的慣性質(zhì)量,M2包括繞y、z軸轉(zhuǎn)動的慣性質(zhì)量。
剛度矩陣K可以分為線性剛度KL以及非線性剛度[9]KN,即
式中,KNA1為軸向力和彎矩耦合作用下,軸向變形對應(yīng)的非線性剛度矩陣;KNA2為軸向力和彎矩耦合作用下,彎曲變形對應(yīng)的非線性剛度矩陣;KNT為扭矩和彎矩耦合作用時的非線性剛度矩陣。
阻尼矩陣分為兩個部分[18],即
式中,CD為Rayleigh阻尼;CN為陀螺阻尼;ωi和ωj分別為系統(tǒng)的第i和第j階固有頻率;ζi和ζj為相應(yīng)于第i和第j階振型的阻尼比,由試驗確定,一般可取i=1,j=2,相應(yīng)的阻尼比為0.02~0.20。
CN是一個反對稱陣,轉(zhuǎn)速越大,其元素值也越大。雖然與Rayleigh阻尼相比,通常情況下其影響并不明顯,但當(dāng)鉆柱的瞬時渦動速度較大時,CN的影響就必須考慮。
在空間和時間上需要進(jìn)行離散,使用有限元方法求解方程(1)。對于空間的離散,本文中使用節(jié)點(diǎn)迭代法[13],其計算流程見圖1。對于時間的離散,采用Newmark方法[13]。
圖1 結(jié)點(diǎn)迭代法計算流程Fig.1 Calculation procedure of node iteration method
RSBHA的組成為:Φ215.9 mm鉆頭+類Auto-Trak旋轉(zhuǎn)導(dǎo)向工具+Φ172 mm鉆鋌×2 m+Φ213 mm穩(wěn)定器+Φ172 mm鉆鋌+Φ127 mm柔性短節(jié)×1.5 m +Φ172 mm鉆鋌+Φ213 mm穩(wěn)定器+Φ172 mm鉆鋌×50 m。
RSBHA的結(jié)構(gòu)如圖2,其中L1、L2、L3表示鉆頭、旋轉(zhuǎn)導(dǎo)向工具、穩(wěn)定器之間的距離,Ls表示柔性短節(jié)和鉆頭之間的距離。柔性短節(jié)為鈦合金材料制成,彈性模量為109 GPa,線重度為444 N/m,長度為1.2 m。
圖2 RSBHA結(jié)構(gòu)示意圖Fig.2 Sketch map of RSHBA structure
其他參數(shù):井眼直徑215.9 mm,井斜角為90°,鉆壓為100 kN,轉(zhuǎn)速為140 r/min,鉆井液密度為1.2 g/cm3。為研究柔性短節(jié)對RSBHA的動力學(xué)特性的影響,現(xiàn)僅改變Ls(6~14 m),固定其余參數(shù)(L1=0.8 m、L2=3 m、L3=10 m)。
圖3為柔性短節(jié)位于不同位置時,RSBHA從2 s到5 s的三維動態(tài)變形圖。從圖3中可以很明顯看出,雖然穩(wěn)定器的個數(shù)、位置都相同,但是有無柔性短節(jié)和柔性短節(jié)位置的改變,使RSBHA的動態(tài)變形也有明顯的不同,其中Ls=6、8、10 m和沒有柔性短節(jié)時,底部鉆具組合的橫向振動幅度較大,柔性短節(jié)的加入和合理的擺放位置能夠減少劇烈的橫向振動,防止RSBHA的運(yùn)動不穩(wěn)定所導(dǎo)致的重要工具失效和測量數(shù)據(jù)準(zhǔn)確性的問題。
圖3 柔性短節(jié)位于不同位置時RSBHA的三維動態(tài)變形Fig.3 3D dynam ic deformation of RSBHA with different flex sub's position
圖4為柔性短節(jié)位于不同位置時,從2 s到5 s RSBHA的應(yīng)力變化。
從圖4中可以看出,動態(tài)應(yīng)力達(dá)到最大時柔性短節(jié)的位置均在距鉆頭14 m處,與有無柔性短節(jié)以及柔性短節(jié)的位置沒有明顯關(guān)系,但應(yīng)力峰值與柔性短節(jié)的位置有較大關(guān)系。在靠近鉆頭3、8、14 m的位置附近均出現(xiàn)了較大的應(yīng)力。隨著柔性短節(jié)與鉆頭之間距離的增大,在靠近鉆頭處出現(xiàn)第二個應(yīng)力峰的位置在不斷遠(yuǎn)離鉆頭,逐漸靠近動態(tài)應(yīng)力最高的14m處;當(dāng)柔性短節(jié)距離鉆頭13 m時,第二個應(yīng)力峰完全與14 m處的應(yīng)力峰重合,從而使柔性短節(jié)距離鉆頭13 m時的最大的動態(tài)應(yīng)力進(jìn)一步增加,達(dá)到了94 MPa。對比無柔性短節(jié)的應(yīng)力圖可以發(fā)現(xiàn),無柔性短節(jié)的最大應(yīng)力為65 MPa,當(dāng)柔性短節(jié)位置Ls=11 m時,最大應(yīng)力為50 MPa,減小了15 MPa,所以柔性短節(jié)的加入能夠較好地降低最大應(yīng)力值,但是其位置相當(dāng)重要,如果擺放位置不合理,可能反而增加其動態(tài)應(yīng)力。當(dāng)沒有柔性短節(jié)和柔性短節(jié)位置Ls為6、8、10 m時,距離鉆頭30~60 m鉆柱的動態(tài)應(yīng)力值變化較大,主要由于此時鉆柱橫向振動較大引起;當(dāng)Ls=7、9、11、13 m時,鉆柱橫向振動較小,30~60 m鉆柱的動態(tài)應(yīng)力分布情況較為理想。圖5、6為給定轉(zhuǎn)速和鉆壓條件下RSBHA動態(tài)應(yīng)力的變化特征。
圖5對應(yīng)轉(zhuǎn)速140 r/min,鉆井液密度為1.2 g/ cm3。從圖5中可以看出,隨著鉆壓的增大,最大應(yīng)力總體增加,但是總體變化規(guī)律基本一致。在不同的鉆壓下,最大應(yīng)力均在柔性短節(jié)位置Ls為12和13 m時有較大值。因此,應(yīng)該盡量避免柔性短節(jié)位置放在12和13 m處。鉆壓的改變對最大應(yīng)力隨柔性短節(jié)位置的變化規(guī)律的影響較小。
圖6對應(yīng)鉆壓100 kN,鉆井液密度為1.2 g/ cm3。從圖6中可以看出,轉(zhuǎn)速的變化對下部鉆具組合上的最大應(yīng)力有一定的影響,但不同轉(zhuǎn)速時最大應(yīng)力隨柔性短節(jié)位置變化的規(guī)律基本一致。當(dāng)轉(zhuǎn)速為60到120 r/min時,柔性短節(jié)位置Ls≤11 m時,最大應(yīng)力較小,當(dāng)11 m<Ls<14 m時,最大應(yīng)力的值較大。當(dāng)轉(zhuǎn)速為140 r/min,柔性短節(jié)的位置Ls=10 m時,最大應(yīng)力較小,甚至小于低轉(zhuǎn)速下同樣結(jié)構(gòu)的RSBHA的最大應(yīng)力,是較為合理的柔性短節(jié)位置。
(1)有無柔性短節(jié)和柔性短節(jié)位置的改變,RSBHA的動力學(xué)特性有明顯的變化。在合理的位置加入柔性短節(jié)能夠減小RSBHA的橫向動態(tài)位移,并且使得最大動態(tài)應(yīng)力明顯下降,起到減振的效果;相反,如果柔性短節(jié)位置不合理,將會明顯增加動態(tài)應(yīng)力。
(2)鉆壓為100~180 kN時,對最大應(yīng)力隨柔性短節(jié)位置變化的規(guī)律影響較小。
(3)轉(zhuǎn)速對底部鉆具組合動態(tài)應(yīng)力有一定影響。柔性短節(jié)位置Ls≤11 m時,最大應(yīng)力較小,當(dāng)11 m<Ls<14 m時,最大應(yīng)力的值較大。
(4)柔性短節(jié)位置對最大動態(tài)應(yīng)力的影響遠(yuǎn)大于鉆壓和轉(zhuǎn)速,合理地確定柔性短節(jié)位置對提高井下工具的安全性和測量數(shù)據(jù)的準(zhǔn)確性有重要作用。
[1]狄勤豐,趙業(yè)榮.導(dǎo)向鉆具組合動力學(xué)方程建立及傳遞函數(shù)求解[J].石油學(xué)報,2000,21(4):87-92.
DIQin-feng,ZHAO Ye-rong.The establishing of dynamics equations and the transfer functions of the steering assembly of the downhole closed-loop drilling system[J].Acta Petrolei Sinica,2000,21(4):87-92.
[2]狄勤豐.旋轉(zhuǎn)導(dǎo)向井下閉環(huán)鉆井技術(shù)[M].西安:陜西科學(xué)技術(shù)出版社,1999:24-26.
[3]張紹槐.現(xiàn)代導(dǎo)向鉆井技術(shù)的新進(jìn)展及發(fā)展方向[J].石油學(xué)報,2003,24(3):82-83.
ZHANG Shao-huai.New progress and development direction ofmodern steering drilling techniques[J].Acta Petrolei Sinica,2003,24(3):82-83.
[4]閆文輝,彭勇,張紹槐.旋轉(zhuǎn)導(dǎo)向鉆井工具的研制原理[J].石油學(xué)報,2005,26(5):94-97.
YANWen-hui,PENG Yong,ZHANG Shao-huai.Mechanism of rotary steering drilling tool[J].Acta Petrolei Sinica,2005,26(5):94-97.
[5]張紹槐,狄勤豐.用旋轉(zhuǎn)導(dǎo)向鉆井系統(tǒng)鉆大位移井[J].石油學(xué)報,2000,21(1):76-80.
ZHANG Shao-huai,DI Qin-feng.Petroleum engineering drilling extended reach well with rotary steering drilling system[J].Acta Petrolei Sinica,2000,21(1):76-80.
[6]胡以寶,狄勤豐,姚建林,等.斜直井眼中鉆柱的動力學(xué)特性分析[J].石油鉆采工藝,2009,31(1):14-17.
HU Yi-bao,DIQin-feng,YAO Jian-lin,et al.Analysis of dynam ic characteristics of drillstring in straight hole[J].Oil Drilling&Production Technology,2009,31 (1):14-17.
[7]胡以寶,狄勤豐,王文昌,等.斜直井眼中轉(zhuǎn)速對鉆柱動力學(xué)特性的影響[J].工程力學(xué),2010,27(5):184-190.
HU Yi-bao,DIQin-feng,WANGWen-chang,et al.Influence of rotary table speed on the dynamic characteristics of drillsting in inclined straight hole[J].Engineering Mechanics,2010,27(5):184-190.
[8]帥健,于永南,呂英民.整體鉆柱的受力與變形分析[J].石油大學(xué)學(xué)報:自然科學(xué)版,1994,18(1):55-60.
SHUAI Jian,YU Yong-nan,LU Ying-min.Analysis of the force and deformation along the whole drillstring[J].Journal of the University of Petroleum,China(Edition of Natural Science),1994,18(1):55-60.
[9]呂英民.有限元法在鉆柱力學(xué)中的應(yīng)用[M].東營:石油大學(xué)出版社,1996.
[10]李茂生,閆相禎,高德利.鉆柱與井壁碰撞的非線性有限元分析[J].石油機(jī)械,2006,34(8):15-18.
LIMao-sheng,YAN Xiang-zhen,GAODe-li.Nonlinear finite elementanalysis on drillstring colloidingwith wellbore[J].China Petroleum Machinery,2006,34(8): 15-18.
[11]DYKSTRA M W.Nonlinear drill string dynamics[D].Tulsa:The University of Tusla,1996.
[12]BATHE K J.Finite element procedures in engineering analysis[M].New Jersey:Prentice-Hall,Inc,1982: 371-604.
[13]祝效華,劉清友,童華.三維井眼全井鉆柱系統(tǒng)動力學(xué)模型研究[J].石油學(xué)報,2008,29(2):288-295.
ZHU Xiao-hua,LIU Qing-you,TONG Hua.Research on dynamicsmodelof full hole drilling-string system with three-dimensional trajectory[J].Acta Petrolei Sinica,2008,29(2):288-295.
[14]龐東曉,劉清友,孟慶華,等.三維彎曲井眼鉆柱接觸非線性問題求解方法[J].石油學(xué)報,2009,30(1): 121-124.
PANG Dong-xiao,LIU Qing-you,MENG Qing-hua,et al.Solvingmethod for nonlinear contact problem of drill strings in 3D curved borehole[J].Acta PetroleiSinica,2009,30(1):121-124.
[15]HAKIMIH,MORADI S.Drillstring vibration analysis using differential quadraturemethod[J].Journal of Petroleum Science and Engineering,2010,70(3):235-242.
[16]ZHUWei-ping,DIQin-feng.Effect of pre-bent deflection on lateral vibration of stabilized drill collars[R].SPE 120455-PA,2011.
[17]崔俊芝,梁俊.現(xiàn)代有限元軟件方法[M].北京:國防工業(yè)出版社,1995.
Effect of flex sub's position on bottom hole assembly w ith rotary steering tool
DIQin-feng1,2,WANG Ming-jie1,2,HU Yi-bao1,2,ZHAO Yu-dong1,2,ZHUWei-ping1,2,WANGWen-chang1,2
(1.Shanghai Instituteof Applied Mathematics and Mechanics,Shanghai University,Shanghai200072,China; 2.Shanghai Key Laboratory of Mechanics in Energy and Environment Engineering,Shanghai University,Shanghai200072,China)
Devoted to the large slenderness ratio of bottom hole assembly(BHA),the finite element iterationmethod was used to analyze the effect of flex sub's position on the 3D dynamic deformation and dynamic stress,and the effectof flex sub's position on themaximal dynamic stress under differentweighton bitand speed was also analyzed.The results show that the lateral dynamic displacement and themaximum dynamic stress in RSBHA can be decreased obviously if the flex sub is put in a reasonable position,then the shock absorption effect of RSBHA can be increased.On the contrary,if the flex sub is put in an unreasonable position,the dynamic stresswill significantly increase.The flex sub's position hasmore impact on themaximum dynamic stress in RSBHA than weighton bitand speed,and how to determine the flex sub's position reasonably is important for improving the security of well tools and the accuracy ofmeasurement data.
bottom hole assembly;flex sub;dynamics characteristics;rotary steering tool
TE 21
A
10.3969/j.issn.1673-5005.2012.05.015
1673-5005(2012)05-0084-05
2012-04-09
國家自然科學(xué)基金項目(51174130;50674065);國家重大專項子課題(2008ZX05024-003);上海市科委項目(061658035);上海市重點(diǎn)學(xué)科建設(shè)項目(S30106);上海市教委科研創(chuàng)新項目(11CXY32)和上海領(lǐng)軍人才基金項目資助
狄勤豐(1963-),男(漢族),江蘇溧陽人,教授,博士生導(dǎo)師,主要從事石油工程技術(shù)及其相關(guān)力學(xué)問題研究。
(編輯 李志芬)