基于鉆柱動力學(xué)的底部鉆具組合疲勞壽命研究*

毛良杰 甘倫科 幸雪松 李 雋

(1. 西南石油大學(xué)油氣藏地質(zhì)及開發(fā)工程國家重點實驗室 2. 中海石油(中國)有限公司北京研究中心 3. 中國石油勘探開發(fā)研究院)

0 引 言

近年來，石油與天然氣勘探開發(fā)的重點已轉(zhuǎn)向超深層/深層，超深井和深井占據(jù)主體井型，而鉆柱失效等井下復(fù)雜事故嚴重制約了此類井的鉆探成本與安全[1]。鉆柱失效形式有斷裂、磨損、變形和刺漏等[2]，一般而言，致使鉆柱失效的主要原因是鉆柱振動誘發(fā)的鉆柱疲勞[3]。其中發(fā)生在底部鉆具組合(BHA)上的多為橫向振動，且是最復(fù)雜和最具破壞力的振動，極易導(dǎo)致BHA過早疲勞[4]。因此，為了防止井下事故的發(fā)生，有必要對BHA的疲勞壽命進行預(yù)測。

累計疲勞損傷法是預(yù)測沒有初始裂紋鉆柱疲勞壽命最常用的方法。然而，由于起下鉆作業(yè)、機械加工、成形工藝和材料等因素，鉆柱不可避免地會產(chǎn)生初始裂紋缺陷[5]，傳統(tǒng)的累計疲勞損傷法不適用于含初始裂紋的鉆柱疲勞壽命預(yù)測。為此，許多學(xué)者進行了考慮鉆柱初始裂紋的疲勞壽命研究。其中，張智等[6]建立了鉆桿接頭疲勞壽命有限元計算模型，在此模型基礎(chǔ)上分析了表面質(zhì)量對疲勞壽命的影響，發(fā)現(xiàn)改善鉆桿接頭螺紋根部的表面質(zhì)量可以大幅延長其疲勞壽命。李文飛等[7]建立了復(fù)合型鉆柱疲勞壽命計算模型，通過算例說明未發(fā)生屈曲的鉆柱疲勞壽命短于發(fā)生屈曲的管段。陳勇等[8]計算了外螺紋接頭的相對疲勞壽命，結(jié)果發(fā)現(xiàn)，抗疲勞性能最佳的螺紋為帶API標準應(yīng)力減輕槽的外螺紋。O.OZGUC[9]通過有限元法對鉆桿接頭處進行建模分析，發(fā)現(xiàn)鉆桿疲勞損傷裂紋多為垂直于鉆桿軸線的平面裂紋。T.A.LEMMA等[10]采用模糊系統(tǒng)對不同鉆桿的無量綱強度因子進行建模，利用模糊系統(tǒng)和累計損傷理論預(yù)測鉆柱疲勞壽命，經(jīng)過分析發(fā)現(xiàn)，鉆柱疲勞壽命受初始裂紋尺寸測量準確性的影響極大。N.H.DAO等[11]利用三維有限元模型結(jié)合Walker模型，研究了旋轉(zhuǎn)鉆進過程中鉆桿表面裂紋的疲勞擴展，其數(shù)值模擬結(jié)果表明，軸向力對鉆桿裂紋擴展有重要作用。

綜上所述，目前關(guān)于鉆柱疲勞壽命的研究，大多是從靜態(tài)角度分析鉆柱的受力情況或者只研究某種振動形式單獨作用下鉆柱的受力情況。而事實上，鉆柱存在復(fù)雜和高度非線性的振動，振動對其疲勞壽命預(yù)測會有極大的影響，因此有必要在預(yù)測BHA疲勞壽命前獲得其動力學(xué)特性。近年來，有限元法以其精度高、適用性廣等眾多優(yōu)勢，成為研究鉆柱動力學(xué)的主要方法。而隨著研究的不斷深入，鉆柱動力學(xué)有限元模型從單一振動形式發(fā)展到多種振動形式的耦合作用，從單自由度系統(tǒng)逐步發(fā)展到多自由度系統(tǒng)，考慮的邊界條件也越來越貼近實際工況[12]。筆者充分考慮鉆柱與井壁的接觸、鉆頭與井底相互作用，建立了鉆柱動力學(xué)有限元模型，分析BHA的動力學(xué)特性；根據(jù)BHA動力學(xué)特性，考慮初始裂紋缺陷，建立了基于鉆柱動力學(xué)的底部鉆具組合疲勞壽命預(yù)測模型，并分析討論了鉆壓、轉(zhuǎn)速、初始裂紋尺寸及穩(wěn)定器安裝位置對BHA疲勞壽命的影響規(guī)律，以期為鉆柱安全提供一定的理論指導(dǎo)。

1 底部鉆具組合疲勞壽命預(yù)測模型的建立

鉆柱在充滿鉆井液的井筒中轉(zhuǎn)動，且具有超大的細長比。其受力特點如下：①上部受到大鉤的拉力和轉(zhuǎn)盤扭矩，下部受到井底與地層互作用產(chǎn)生的壓力和扭矩；②整個鉆柱都受到鉆井液帶來的黏滯阻尼力；③在與井壁發(fā)生接觸時，還會受到垂直于鉆柱軸線的沖擊力和井壁的摩擦力。針對鉆柱這些受力特點，本節(jié)基于哈密頓原理并充分考慮鉆柱與井壁的接觸、鉆頭與井底互作用，建立6自由度鉆柱動力學(xué)有限元模型，借此得到鉆柱的動力學(xué)特性；再根據(jù)Walker模型計算BHA的疲勞壽命。圖1為鉆柱系統(tǒng)的物理模型。

圖1 鉆柱系統(tǒng)物理模型Fig.1 Physical model of drill string system

1.1 鉆柱動力學(xué)模型的建立

鉆柱運動可以用單元節(jié)點位移向量來表示：

U=(u1，v1，w1，θx1，θy1，φ1，

u2，v2，w2，θx2，θy2，φ2)T

(1)

式中：U是鉆柱單元的位移矩陣；u1、v1、w1分別是鉆柱單元第1個節(jié)點X、Y、Z方向的平動位移；θx1、θy1、φ1分別是鉆柱單元第1個節(jié)點相對于X軸、Y軸和Z軸的轉(zhuǎn)角；u2、v2、w2分別是鉆柱單元第2個節(jié)點X、Y、Z方向的平動位移；θx2、θy2、φ2分別是鉆柱單元第2個節(jié)點相對于X軸、Y軸和Z軸的轉(zhuǎn)角。

拉格朗日動力方程為[13]：

(2)

單元體總動能T表達式為：

(3)

在式(3)中，前3項表示梁單元的平移動能，后1項表示梁單元的轉(zhuǎn)動動能。

單元體總勢能V表達式為：

(4)

將式(3)式(4)代入式(2)，即可得到梁單元的動力方程：

(5)

1.1.1 鉆柱與井壁接觸

假設(shè)井壁為剛性，鉆柱與井壁實質(zhì)是隨機接觸。軸向摩擦阻力Fez、切向摩擦阻力Fer、扭矩Te及彎矩Me分別為[14]：

(6)

接觸力Fr可以表示為[14]：

(7)

式中：do為井筒直徑，m；r為鉆柱半徑，m；vr是鉆柱的徑向速度，m/s；ur為鉆柱的徑向位移，m；kh為井壁的剛度，N/m；v1和v2分別為節(jié)點碰撞前、后速度，m/s。

式(7)中的負號表示在鉆柱與井壁接觸的瞬間，鉆柱受到的徑向力與鉆柱的運動方向相反。

1.1.2 鉆柱動力學(xué)模型邊界條件

有文獻指出，隨著與鉆頭距離的增加，鉆柱橫向振動加速度幅值逐漸減小，其峰值位于鉆頭附近，由于橫向振動在傳播過程中被消耗，傳遞不到底部鉆具組合頂端。同時，為了提高計算效率，只取底部鉆具組合為研究對象，所以將上部邊界定在底部鉆具組合的頂端[15-16]，則上部邊界條件為：橫向自由度被約束，鉆柱只受軸向力和扭矩作用。

鉆頭受到與地層相互作用產(chǎn)生的扭矩Tb和激振力P，則下部邊界條件為：鉆頭鉸支在井底，在鉆頭處鉆柱的橫向位移受到約束，但繞各軸的轉(zhuǎn)動卻不受約束。

激振力和扭矩計算公式為[17]：

(8)

式中：nb為與鉆頭類型有關(guān)的系數(shù)，PDC鉆頭取1，牙輪鉆頭取3；P0為激振力幅值，N；ω為鉆頭旋轉(zhuǎn)角速度，rad/s；t為鉆頭旋轉(zhuǎn)時間，s；Nb為鉆頭轉(zhuǎn)速，r/min；h為鉆頭進尺，m；Db為鉆頭直徑，mm；μb為鉆頭與地層間的摩擦因數(shù)，無量綱；PWOB為鉆壓，N。

1.2 底部鉆具組合疲勞壽命預(yù)測

常用的疲勞壽命計算模型主要有Paris模型、Walker模型和Forman模型。其中Walker模型的參數(shù)更多，且均由現(xiàn)場實際數(shù)據(jù)計算得到，計算結(jié)果更精準，所以本文使用Walker模型作為疲勞壽命預(yù)測模型。

3軸應(yīng)力下的鉆柱受力在小變形情況下滿足第四強度理論[18]，可以應(yīng)用Mises等效應(yīng)力理論進行應(yīng)力的合成，故令：σ1max=σz+σb，σ1min=σz-σb，σ2=σr；σ3=σt，于是最大、最小Mises等效應(yīng)力為：

(9)

Tresca等效應(yīng)力為：

(10)

所以，平均應(yīng)力、應(yīng)力幅、應(yīng)力比為：

(11)

初始裂紋尺寸可以由無損探傷儀測定，考慮到鉆柱初始缺陷的影響，在疲勞壽命計算過程中，需要先計算出幾何形狀因子及應(yīng)力強度因子[19]，在此基礎(chǔ)上，臨界裂紋尺寸可根據(jù)應(yīng)力條件確定，再利用Walker模型計算鉆柱的疲勞壽命。

在使用和制造鉆柱的過程中大多產(chǎn)生表面線性裂紋，鉆柱表面線性裂紋的幾何形狀因子的計算方法為：

Fm=1+0.128(aL/Dp)-
0.288(aL/Dp)2+1.525(aL/Dp)3

(12)

式中：Fm為裂紋形狀因子，無因次；aL為裂紋半長，m；Dp為鉆柱直徑，m。

最大應(yīng)力強度因子為：

(13)

式中：Kmax為最大應(yīng)力強度因子，無因次；σmax為最大等效應(yīng)力，MPa；a為裂紋尺寸，m。

初始裂紋深度a0由探傷儀直接測得，臨界裂紋深度ac與材料的斷裂韌度和受力情況有關(guān)，其計算式為：

(14)

式中：KIC為材料斷裂韌度，MPa·m1/2；τmax為最大切應(yīng)力，MPa；μ為泊松比，無因次。

Walker模型為：

(15)

式中：N為應(yīng)力循環(huán)周次，次；c、n、m均為材料常數(shù)。

鉆柱從初始裂紋深度a0擴展到臨界裂紋深度ac的循環(huán)周次(即疲勞壽命)為：

(16)

1.3 底部鉆具組合疲勞壽命預(yù)測模型求解

利用Newmark-β法求解鉆柱動力學(xué)模型，得到位移、速度和加速度，據(jù)此計算出最大等效應(yīng)力、應(yīng)力比和應(yīng)力幅；在此基礎(chǔ)上，利用Walker模型計算疲勞壽命，具體的求解流程如圖2所示。

圖2 底部鉆具組合疲勞壽命預(yù)測模型求解流程圖Fig.2 Flow chart for solving fatigue life prediction model of BHA

2 底部鉆具組合疲勞壽命預(yù)測模型準確性驗證

本文以表1所列參數(shù)為基礎(chǔ)，對文獻[20]報道的鉆桿疲勞壽命試驗結(jié)果進行了預(yù)測，預(yù)測結(jié)果如表2所示。對比試驗實際檢測結(jié)果、本模型預(yù)測結(jié)果及文獻[20]預(yù)測結(jié)果，可以發(fā)現(xiàn)本模型準確可靠，比從靜態(tài)角度分析受力得到的預(yù)測模型精度更高。

表1 文獻[20]中鉆桿疲勞破壞試驗基本參數(shù)Table 1 Basic parameters of drill pipe fatigue-to-failure test in Reference [20]

表2 鉆桿疲勞破壞試驗結(jié)果對比情況Table 2 Comparison of fatigue-to-failure test results of drill pipe

3 案例分析

本研究選取渤海中部一口直井(B1井)進行分析，該井鉆柱系統(tǒng)的主要參數(shù)(源于現(xiàn)場資料)如表3所示。數(shù)值模擬中，將鉆柱分為若干個單元，每個單元長3 m；模擬時間為5 s，時間步長為0.001 s。

表3 B1井底部鉆具組合主要參數(shù)Table 3 Main parameters of BHA in Well B1

3.1 動力學(xué)特性分析

圖3是BHA的橫向變形與軸向力分布圖。從圖3可以看出：BHA變形程度越靠近鉆頭越大；距離鉆頭244 m以上部分鉆柱橫向變形更小。這是因為中和點上部的鉆柱受軸向拉力，穩(wěn)定性好，不易變形；中和點下部鉆柱受軸向壓力，穩(wěn)定性差，極易發(fā)生變形，所以BHA下部在外力作用下振動更為劇烈。

圖3 BHA橫向變形與軸向力Fig.3 Lateral deformation and axial force of BHA

圖4為BHA疲勞壽命分布圖。從圖4可以看出：從距離鉆頭350 m到距離鉆頭244 m，BHA疲勞壽命逐漸縮短，在244 m處迅速縮短到一個極低的值(低于0.1×108次)；在接近244 m位置處，BHA疲勞壽命突增到6×108次左右，然后又開始逐漸縮短。這是因為距離鉆頭越近，BHA振動越劇烈，劇烈的振動加大了BHA所承受的交變應(yīng)力的幅值及頻率，所以BHA疲勞壽命逐漸縮短；在距離鉆頭244 m處，BHA橫截面積發(fā)生了變化，產(chǎn)生了應(yīng)力集中，交變應(yīng)力幅值增大，此處BHA疲勞壽命迅速縮短到一個極低的值；而距離鉆頭小于244 m后，BHA橫截面積增大，臨界裂紋尺寸減小，疲勞壽命也就延長。

圖4 BHA疲勞壽命分布(鉆壓80 kN、轉(zhuǎn)速60 r/min)Fig.4 Fatigue life distribution of BHA (WOB: 80 kN,rotary speed: 60 r/min)

通過以上分析，鉆柱振動的加劇會使得疲勞壽命迅速縮短。鉆井工程中鉆壓和轉(zhuǎn)速是最重要的兩個施工參數(shù)，合理地控制鉆壓和轉(zhuǎn)速可以有效地延長鉆柱疲勞壽命，避免鉆柱過早失效；穩(wěn)定器的存在不僅會使近鉆頭處鉆鋌保持居中，而且會抑制近鉆頭處鉆鋌的振動；鉆柱初始裂紋深度也會對鉆柱疲勞壽命產(chǎn)生巨大的影響。因此，筆者進一步研究了不同轉(zhuǎn)速、鉆壓、初始裂紋深度及穩(wěn)定器安裝位置下，BHA最大等效應(yīng)力和疲勞壽命的分布情況，得到轉(zhuǎn)速、鉆壓、初始裂紋深度及穩(wěn)定器安裝位置對BHA疲勞壽命的影響規(guī)律。

3.2 初始裂紋深度對BHA疲勞壽命的影響

圖5為不同初始裂紋深度下BHA疲勞壽命。

圖5 不同初始裂紋深度下BHA疲勞壽命Fig.5 Fatigue life of BHA under different initial crack depths

從圖5a可以看出：具有更大初始裂紋深度的BHA，其疲勞壽命會更短。這是因為初始裂紋深度越大，臨界裂紋深度與初始裂紋深度之差越小，裂紋擴展到臨界裂紋深度的應(yīng)力循環(huán)次數(shù)(即疲勞壽命)越少。此外，從圖5b可以發(fā)現(xiàn)：初始裂紋深度從0增大到1 mm時，疲勞壽命從1.55×108次迅速縮減到0.4×108次；超過1 mm后隨初始裂紋深度繼續(xù)增大，疲勞壽命緩慢縮減到0.2×108次。因此，建議在生產(chǎn)和使用鉆柱的過程中，應(yīng)當經(jīng)常用無損探傷儀檢測鉆柱的裂紋深度，避免因下入初始裂紋深度過大的鉆柱而導(dǎo)致鉆柱在井下作業(yè)過程中發(fā)生疲勞斷裂。

3.3 轉(zhuǎn)速對BHA疲勞壽命的影響

圖6為不同轉(zhuǎn)速下BHA最大等效應(yīng)力和疲勞壽命分布曲線。

圖6 不同轉(zhuǎn)速下BHA各位置最大等效應(yīng)力和疲勞壽命分布曲線Fig.6 Distribution of maximum equivalent stress and fatigue life at each position of BHA at different rotary speeds

從圖6可以看出：轉(zhuǎn)速越大，BHA的最大等效應(yīng)力越大、疲勞壽命越短；在低轉(zhuǎn)速范圍內(nèi)(30～90 r/min)，轉(zhuǎn)速對BHA最大等效應(yīng)力和疲勞壽命的影響較小，在高轉(zhuǎn)速范圍內(nèi)(120～180 r/min)，轉(zhuǎn)速對BHA最大等效應(yīng)力和疲勞壽命的影響較大。這是因為在低轉(zhuǎn)速范圍內(nèi)，轉(zhuǎn)速的增加雖然會增大BHA所受離心力，但同樣也會導(dǎo)致阻尼力的增加，從而抑制鉆柱的振動。在離心力與阻尼力的綜合作用下，轉(zhuǎn)速對BHA最大等效應(yīng)力的影響并不明顯。而在高轉(zhuǎn)速范圍內(nèi)，轉(zhuǎn)速增加會大幅度增大鉆柱的離心力，過大的離心力可能加劇鉆柱與井壁間的碰撞而導(dǎo)致鉆柱的等效應(yīng)力激增(見圖6a)，所以轉(zhuǎn)速增大，鉆柱疲勞壽命明顯縮短。據(jù)此建議在正常鉆井時，可以適當選擇較高的轉(zhuǎn)速以提高機械鉆速，但在強化鉆井時，要對BHA的疲勞壽命重點關(guān)注。

3.4 鉆壓對底部鉆具組合疲勞壽命的影響

圖7為不同鉆壓下BHA各位置最大等效應(yīng)力和疲勞壽命分布曲線。從圖7可以看出：中和點以下部分，鉆壓越大，BHA承受的最大等效應(yīng)力越大，疲勞壽命越短。這是因為鉆壓的增加會增加處于壓縮狀態(tài)鉆柱的彎曲應(yīng)力，同時鉆壓增加又會降低處于拉伸狀態(tài)鉆柱的拉伸應(yīng)力。因此在大鉆壓下，中和點以下鉆柱的最大等效應(yīng)力較大。中和點以上部分則相反，疲勞壽命的變化情況與之對應(yīng)。由此可知鉆壓對BHA疲勞壽命影響顯著，在滿足高機械鉆速要求的鉆壓范圍內(nèi)，建議選用較低鉆壓。

圖7 不同鉆壓下BHA各位置最大等效應(yīng)力和疲勞壽命分布曲線Fig.7 Distribution of maximum equivalent stress and fatigue life at each position of BHA under different weights on bit

3.5 穩(wěn)定器安裝位置對BHA疲勞壽命的影響

圖8為不同穩(wěn)定器安裝位置下BHA近鉆頭處最大等效應(yīng)力和疲勞壽命分布曲線。從圖8可以發(fā)現(xiàn)：安裝穩(wěn)定器與不安裝穩(wěn)定器相比，BHA近鉆頭處最大等效應(yīng)力降低，疲勞壽命延長；穩(wěn)定器安裝在距離鉆頭9 m處時，延長疲勞壽命的效果最佳。這是因為安裝穩(wěn)定器后，BHA被約束，橫向變形減小(見圖9)；橫向變形減小使得BHA振動減弱，最大等效應(yīng)力減小，疲勞壽命延長。

圖8 不同穩(wěn)定器安裝位置下BHA近鉆頭處最大等效應(yīng)力和疲勞壽命分布曲線Fig.8 Distribution of maximum equivalent stress and fatigue life of BHA near the bit at different installation positions of stabilizer

圖9 不同穩(wěn)定器安裝位置下BHA近鉆頭處橫向變形Fig.9 Lateral deformation of BHA near the bit at different installation positions of stabilizer

此外，在圖9還可以發(fā)現(xiàn)，穩(wěn)定器安裝在距離鉆頭9 m時，穩(wěn)定器到鉆頭這一段BHA的橫向變形最小，所以穩(wěn)定器安裝在距離鉆頭9 m時，延長BHA疲勞壽命的效果最佳。在實際鉆井過程中，可以通過安裝穩(wěn)定器來延長近鉆頭處BHA的疲勞壽命，但需綜合考慮鉆具特性、鉆柱振動以及穩(wěn)定器特性來確定合適的穩(wěn)定器安裝位置。

4 結(jié) 論

(1)BHA橫向振動比縱向振動更激烈，橫向振動是導(dǎo)致BHA疲勞的主要原因；越靠近鉆頭，BHA變形程度越大，橫向振動越劇烈，疲勞壽命也越短；鉆柱橫截面積變化使得鉆柱產(chǎn)生應(yīng)力集中，劇烈振動和應(yīng)力集中的共同作用會極大縮短鉆柱疲勞壽命。

(2)初始裂紋深度對BHA疲勞壽命影響明顯，具有更大初始裂紋深度的BHA，其疲勞壽命會更短；從無初始裂紋增加到裂紋深度為1 mm，BHA疲勞壽命迅速縮短，隨后初始裂紋深度繼續(xù)增加，BHA疲勞壽命緩慢縮短。

(3)鉆壓對BHA疲勞壽命影響顯著，隨著鉆壓的提高，中和點以下部分BHA疲勞壽命縮短；轉(zhuǎn)速越大BHA疲勞壽命越短，在低轉(zhuǎn)速范圍內(nèi)轉(zhuǎn)速對BHA疲勞壽命影響較小，在高轉(zhuǎn)速范圍內(nèi)轉(zhuǎn)速對BHA疲勞壽命影響較大；安裝穩(wěn)定器可以延長近鉆頭處BHA疲勞壽命，在本文案例中，穩(wěn)定器安裝在距離鉆頭9m處，穩(wěn)定器到鉆頭的橫向變形最小，延長疲勞壽命的效果最佳。