高光海，仇性啟，董 輝，許俊良，任 紅

(1.中國石油大學(華東) 化學工程學院，山東 青島 266580；2.中石化勝利石油工程有限公司 鉆井工藝研究院，山東 東營 257017)

?

南海超深水取樣鉆柱縱向粘滑振動分析

高光海1，仇性啟1，董 輝1，許俊良2，任 紅2

(1.中國石油大學(華東) 化學工程學院，山東 青島 266580；2.中石化勝利石油工程有限公司 鉆井工藝研究院，山東 東營 257017)

針對南海超深水天然氣水合物取樣鉆柱縱向粘滑振動問題，推導了縱向振動微分方程，采用非線性彈簧單元模擬了鉆柱與巖層之間的粘滑摩擦。通過數值求解得到了鉆柱縱向振動頻率及振動位移，探討了鉆壓、鉆鋌長度及作業(yè)水深對鉆柱縱向振動的影響。結果表明：鉆柱低階振動頻率所對應外載激勵周期為1.5～4.0 s，接近南海波浪運動周期，鉆柱發(fā)生低階共振幾率較大；鉆壓的增加，使得鉆柱振動位移減?。汇@鋌長度及作業(yè)水深的增加使得鉆柱振動頻率減小、振動位移增加，南海環(huán)境下鉆柱更易發(fā)生低階共振；鉆柱發(fā)生1階共振時振動位移約為近頻率振動位移的8倍，共振危害程度較大。

超深水；鉆柱；縱向粘滑振動；振動頻率；振動位移

我國南海深水海域天然氣水合物鉆探取樣時，取樣鉆柱頂端與鉆井船連接，通過鉆井船中動力裝置為鉆柱提供動力；海水中無隔水管裝置，鉆柱與海水直接接觸；海底巖層中鉆柱直接與巖層接觸。鉆柱長徑比較大，鉆柱彎曲撓度較大，鉆柱結構具有幾何與接觸雙重非線性性質。鉆柱在鉆井船縱向浮沉運動、鉆柱與海底巖層之間縱向摩擦作用、鉆柱底部鉆頭與巖層的相互作用等載荷作用下，其縱向振動特性十分復雜。取樣鉆柱縱向振動劇烈程度直接關系到鉆柱內部取樣儀器能否正常定位以及鉆探取樣過程能否順利實施。

目前，國內外諸學者對海洋鉆井立管及陸地石油鉆井中鉆柱的縱向振動特性進行了相關研究，取得了豐碩成果。Y Chen[1]等通過變分迭代計算研究了海洋立管固有頻率特性；Y Kovalyshen[2]研究了深井鉆井中鉆柱粘滑振動特性；A Ghasemloonia[3]等研究了鉆柱縱向及橫向耦合非線性振動特性；Y Hu[4]等采用空間曲梁有限單元模型研究了超深井中鉆柱動力學特性；S M Sahebkar[5]等對斜直井中鉆柱在隨時間變化的軸向載荷作用下的非線性振動特性進行了分析；T G Ritto[6]等采用鉆頭-巖層非確定性模型分析了鉆柱非線性動力學特性；Y A Khulief[7]等對鉆柱粘滑自激振動特性進行了研究；R Sampaio[8]等通過建立分線性模型研究了鉆柱縱向及扭轉耦合振動特性；Y A Khulief[9]等通過拉格朗日與有限單元法相結合分析了鉆柱動力學振動特性。國內劉清友、李子豐[10-11]等建立了鉆柱縱向振動力學模型及理論計算公式。章揚烈[12]通過大量的理論分析及試驗研究，得到了鉆柱縱向振動機理及相關減振措施。王文龍[13]等研究了激勵位移、鉆井液粘度及減震器安裝位置等因素對鉆柱縱向振動軸向力的影響規(guī)律。趙宗彬、肖文生[14-15]等通過ANSYS有限元分析，研究了鉆柱縱向振動固有特性并得出了相應的結論。劉進田、閆鐵[16-17]等推導了鉆柱振動的微分方程及有限元矩陣方程，并進行了實例分析計算。孫千偉[18]等通過解析法以及相關實驗研究了管-土軸向動力作用下等效彈簧系數取值問題，確定了不同地質條件下等效彈簧系數取值范圍。

鑒于取樣鉆柱結構及工藝流程的特殊性，基于陸地鉆柱及海洋立管縱向振動特性研究現狀，建立取樣鉆柱縱向振動力學分析模型，研究鉆柱結構的縱向振動特性。

1 取樣鉆柱縱向振動微分方程

1.1 基本假設

建立取樣鉆柱縱向振動力學模型時對一些次要因素進行相應的簡化，做如下假設：

1) 鉆柱材料均勻且各向同性，在運動和變形過程中始終處于線彈性范圍之內。

2) 鉆柱橫截面為圓環(huán)形，忽略鉆柱接頭的影響。

3) 鉆柱內部充滿海水，鉆柱外部水流為定常流動載荷。

4) 同一縱向位置處，鉆柱振動過程所受阻尼以及巖層與鉆柱之間摩擦因數皆恒定。

1.2 動力學微分方程

取樣鉆柱縱向振動模型示意圖如圖1所示。

圖1 鉆柱縱向振動模型示意

假設鉆柱材料密度為ρ，截面抗拉剛度為E(z)A(z)，阻尼系數為c(z)。假設鉆柱橫截面在縱向振動過程中始終保持平面，以u(z,t)表示z截面縱向位移，它是截面位置z與時間t的二元函數。取微元段dz，其應變量為：

根據牛頓第二定律可得鉆柱縱向振動微分方程[10-11]：

(1)

式中：E(z)為位置z處鉆柱彈性模量；A(z)為位置z處鉆柱橫截面積；ρ′為鉆井液(海水)密度；g為重力加速度；fm為鉆柱微元所受軸向粘滑摩擦力。

1.3 初始條件

1.4 邊界條件

1) 海平面處。

海平面處假定鉆柱與鉆井船連接一起，縱向方向上隨鉆井船一起做浮沉運動，即：

u(0,t)=U(t)

(2)

其中，u(0,t)為t時刻海平面處鉆柱縱向位移；U(t)為鉆井船縱向運動位移。

2) 鉆頭處。

鉆井船縱向周期性浮沉運動以及底部巖層的不均勻性使得鉆頭處所受縱向載荷不斷發(fā)生變化。根據牛頓第二定律，鉆頭與巖層互作用縱向動力學模型為：

(3)

其中，mb為鉆頭質量；uz(t)為t時刻鉆頭縱向振動位移；fz(t)為t時刻鉆頭所受到的巖石作用力；p(t)為t時刻作用到鉆頭上的鉆壓。

3) 鉆柱與海底巖層接觸處。

通過靜力學迭代計算分析，確定鉆柱系統(tǒng)與巖層接觸部位。將鉆柱與巖層之間的摩擦作用等效為非線性彈簧對鉆柱作用，將鉆柱節(jié)點與非線性彈簧單元節(jié)點耦合。通過設置彈簧單元的力學參數，來模擬鉆柱與巖層之間的粘滑摩擦作用。

2 數值分析

采用有限單元法將鉆柱連續(xù)系統(tǒng)離散為n個單元的離散系統(tǒng)，這些單元組成一個有限多自由度系統(tǒng)。

取樣鉆柱系統(tǒng)振動過程中，始終存在阻尼作用。鉆柱縱向振動系統(tǒng)所受阻尼采用工程設計中常用Rayleigh阻尼：

[C]=α[M]+β[K]

(4)

式中：[C]為阻尼矩陣；[M]為質量矩陣；[K]為剛度矩陣；α為鉆柱質量阻尼系數；β為鉆柱剛度阻尼系數。

α、β取值由鉆柱振型阻尼比及振動固有角頻率計算得出。

ξ=α/2ωi+βω2/2

(5)

式中：ξ為鉆柱振型阻尼比，一般取0.02～0.05；ωi為鉆柱第i階振動固有角頻率。

采用ANSYS靜、動力學分析模塊結合MATLAB、APDL編程模塊對取樣鉆柱縱向振動微分方程進行數值求解。選取以梁單元特性為基礎的PIPE59及PIPE16單元分別模擬海水及巖層中鉆柱，選取COMBIN39非線性彈簧單元模擬巖層與鉆柱粘滑摩擦，選取MASS21單元模擬鉆頭。

3 縱向振動頻率分析

3.1 實例分析

鉆柱結構及環(huán)境載荷相關參數如表1。

表1 鉆柱及環(huán)境相關參數

根據表1參數，求得鉆柱縱向振動頻率如圖2所示。

由圖2可知，隨著階數增加，鉆柱振動頻率逐漸增加。鉆柱前幾階振動頻率較小，外界載荷作用下，鉆柱易發(fā)生縱向共振。

圖2 鉆柱縱向振動頻率

3.2 影響因素分析

1) 鉆鋌長度。

作業(yè)水深2 000 m、巖層深度600 m時，不同鉆鋌長度下，鉆柱縱向振動頻率如表2。

由表2可知，隨著鉆鋌長度的增加，鉆柱各階縱向振動頻率逐漸減小，即鉆鋌長度的增加使得鉆柱縱向振動周期增加。

2) 作業(yè)水深。

巖層深度600 m，不同作業(yè)水深時，鉆柱縱向振動頻率值如圖3所示。

表2 鉆柱縱向振動頻率

圖3 不同作業(yè)水深時鉆柱縱向振動頻率

由圖3可知，隨著作業(yè)水深的增加，鉆柱縱向振動各階振動頻率不斷減小，且隨著振動階數的增加，變化幅度越來越大。作業(yè)水深的增加，使得鉆柱縱向振動周期增加。

4 縱向振動位移響應分析

取樣鉆柱鉆探取樣過程中，鉆柱縱向振動位移大小直接關系到鉆柱內部取樣儀器的正常定位。振動位移過大將使得鉆柱內部取樣儀器不能正常定位，進而使得取樣失敗。鉆探取樣過程中，不同因素對鉆柱振動位移影響不同。以下分析不同因素對鉆柱縱向振動位移的影響。

4.1 鉆壓

作業(yè)水深2 000 m、巖層深度600 m、鉆鋌長度取80 m時，不同鉆壓作用下，鉆柱縱向振動最大位移響應曲線如圖4所示。

圖4 不同鉆壓下鉆柱縱向振動位移曲線

由圖4可知，隨著鉆壓的增加，鉆柱縱向振動位移逐漸減小，鉆柱縱向振動位移最大值所對應頻率為鉆柱1階振動頻率。相同鉆柱結構下，鉆壓的增加，使得鉆柱下部受壓區(qū)域增大，鉆柱縱向振動受到抑制，振動位移減??；鉆柱發(fā)生1階共振時，共振位移突然增大，共振位移值約為相近頻率正常振動下的8倍。

4.2 鉆鋌長度

作業(yè)水深2 000 m、巖層深度600 m、鉆壓取200 kN時，不同鉆鋌長度下，鉆柱縱向振動最大位移響應曲線如圖5所示。

圖5 不同鉆鋌長度下鉆柱縱向振動位移曲線

由圖5可知，隨著鉆鋌長度增加，鉆柱縱向振動位移逐漸增加，鉆柱縱向振動位移最大值所對應頻率為鉆柱1階振動頻率。相同鉆柱長度及鉆壓條件下，鉆鋌長度的增加，使得鉆柱下部受拉區(qū)域增加，鉆柱縱向振動增強，振動位移增大；鉆柱1階共振位移值約為相近頻率正常振動下的8倍。

4.3 作業(yè)水深

巖層深度600 m、鉆鋌長度250 m、鉆壓取200 kN時，不同作業(yè)水深時，鉆柱縱向振動最大位移響應曲線如圖6所示。

圖6 不同作業(yè)水深時鉆柱位移曲線

由圖6可知，隨著作業(yè)水深的增加，鉆柱縱向振動位移逐漸增加，鉆柱縱向振動位移最大值所對應頻率為鉆柱1階振動頻率。相同鉆壓條件下，作業(yè)水深的增加，使得鉆柱下部受拉區(qū)域逐漸增加，鉆柱縱向振動增強，振動位移增大。

5 結論

1) 推導了取樣鉆柱縱向振動微分方程。鉆柱與巖層縱向粘滑摩擦作用可通過非線性彈簧單元進行模擬。通過數值求解得出了鉆柱縱向振動頻率及振動響應位移。

2) 南海深水海域取樣鉆柱1階振動頻率所對應外載荷激勵周期在1.5～4.0 s。鉆鋌長度及作業(yè)水深的增加使得鉆柱縱向振動頻率減小，在南海環(huán)境條件下鉆柱更易發(fā)生1階共振。

3) 取樣鉆柱發(fā)生1階共振時振動位移為近頻率正常振動位移的8倍。相同條件下，鉆壓的增加使得鉆柱振動位移逐漸減小，鉆鋌長度及作業(yè)水深的增加使得鉆柱振動位移逐漸增加。

[1] Chen Y，Zhang J，Zhang H，et al.Re-examination of natural frequencies of marine risers by variational iteration method[J].Ocean Engineering，2015( 94)：132-139.

[2] Kovalyshen Y.Understanding root cause of stick-slip vibrations in deep drilling with drag bits[J].International Journal of Non-Linear Mechanics，2014( 67)：331-341.

[3] Ghasemloonia A，Rideout D G，Butt S D.Analysis of multi-mode nonlinear coupled axial-transverse drillstring vibration in vibration assisted rotary drilling[J].Journal of Petroleum Science and Engineering，2014(116)：36-49.

[4] Hu Y，Di Q，Zhu W，et al.Dynamic characteristics analysis of drillstring in the ultra-deep well with spatial curved beam finite element[J].Journal of Petroleum Science and Engineering，2012(82)：166-173.

[5] Sahebkar S M，Ghazavi M R，Khadem S E，et al.Nonlinear vibration analysis of an axially moving drillstring system with time dependent axial load and axial velocity in inclined well[J].Mechanism and Machine Theory，2011，46(5)：743-760.

[6] Ritto T G，Soize C，Sampaio R.Non-linear dynamics of a drill-string with uncertain model of the bit-rock interaction[J].International Journal of Non-Linear Mechanics，2009，44(8)：865-876.

[7] Khulief Y A，Al-Sulaiman F A，Bashmal S.Vibration analysis of drillstrings with self-excited stick-slip oscillations[J].Journal of sound and vibration，2007，299(3)：540-558.

[8] Sampaio R，Piovan M T，Lozano G V.Coupled axial/torsional vibrations of drill-strings by means of non-linear model[J].Mechanics Research Communications，2007，34(5)：497-502.

[9] Khulief Y A，Al-Naser H.Finite element dynamic analysis of drillstrings[J].Finite elements in analysis and design，2005，41(13)：1270-1288.

[10] 劉清友，馬德坤，湯小文.鉆柱縱向振動模型的建立及求解方法[J].西南石油大學學報(自然科學版)，1998，20(4)：55-58.

[11] 李子豐，李志剛.鉆柱縱向振動分析[J].天然氣工業(yè)，2004，24(6)：70-73.

[12] 章揚烈.鉆柱運動學與動力學[M].北京：石油工業(yè)出版社，2001.

[13] 王文龍，胡群愛，劉化偉，等.鉆柱縱向振動分析與應用[J].振動與沖擊，2011，30(6)：229-233.

[14] 趙宗彬，仇性啟，許俊良.深海天然氣水合物鉆探取心鉆柱振動模態(tài)分析[J].天然氣工業(yè)，2011，31(1)：73-76.

[15] 肖文生，王現鋒，裴艷麗，等.基于ANSYS的鉆柱縱向振動固有頻率分析[J].鉆采工藝，2011，34(5)：93-95.

[16] 劉進田.下部鉆柱振動分析及計算[J].石油礦場機械，2013，42(3)：50-54.

[17] 閆鐵，王雪剛，李杉，等.鉆柱軸向與橫向耦合振動的有限元分析[J].石油礦場機械，2012，41(3)：39-42.

[18] 孫千偉，劉威，李杰.管-土軸向動力相互作用等效彈簧系數取值[J].同濟大學學報(自然科學版)，2012，40(8)：1123-1128.

Analysis on Longitudinal Stick-slip Vibration of Ultra-deep Water Drillstring

GAO Guanghai1，QIU Xingqi1，DONG Hui1，XU Junliang2，REN Hong2

(1.CollegeofChemicalEngineering，ChinaUniversityofPetroleum(EastChina)，Qingdao266580，China；2.DrillingTechnologyResearchInstitute，SinopecOilfieldServiceShengliCorporation，Dongying257017，China)

Based on the longitudinal stick-slip vibration of the drillstring in South China Sea，the differential equation of the longitudinal vibration was deduced.The stick-slip friction of the drillstring and rock was simulated through nonlinear spring unit.The vibration frequency and displacement were obtained by using numerical simulation method.The effects of drilling pressure，the length of the drill collar and the depth of water were also discussed in the paper.According to the results，the period of the external load excitation which corresponding to the low-level vibration frequencies of the drillstring is 1.5～4.0 s.The period is close to the period of the wave in south sea of China.The probability of the resonance of the drillstring is large.The vibration displacements decrease owing to the increase of the drilling pressure.The vibration frequencies decrease and the vibration displacements increase owing to the increase of the length of the drill collar and water depth.The resonance displacement is nearly 8 times to the displacement of normal vibration.

ultra-deep water；drill string；longitudinal stick-slip vibration；vibration frequency；vibration displacement

2016-05-18

國家高技術研究發(fā)展計劃(863計劃)課題“天然氣水合物鉆探取芯工程樣機及配套技術”(2013AA092602)；中央高?；究蒲袠I(yè)務費專項資金資助項目(16CX06015A)

高光海(1989-)，男，山東禹城人，博士研究生，從事海洋油氣管柱力學研究，E-mail：ggh8912@126.com。

1001-3482(2016)11-0008-06

TE952

A

10.3969/j.issn.1001-3482.2016.11.002