考慮巖屑影響的大位移井機(jī)械延伸極限研究

郭驍 李思洋 周蒙恩 孟慶威 張菲菲

1.長江大學(xué)石油工程學(xué)院；2.川慶鉆探工程有限公司長慶鉆井總公司；3.大慶鉆探工程公司地質(zhì)錄井一公司

大位移井的井眼軌跡在地層中水平延伸較長距離，可以增大油氣層的裸露面積，大位移井眼連通多個儲層，使得后期的儲層改造和增產(chǎn)較為便利，增產(chǎn)效果明顯。大位移井的應(yīng)用減少了鉆井的數(shù)量，節(jié)約了開發(fā)成本，提高了效益。但是在大位移井施工中，較多的大位移井下鉆作業(yè)存在困難，所以對管柱的摩阻扭矩與大位移井極限延伸的研究是解決此類問題重要的一環(huán)。

鉆柱力學(xué)方面，國內(nèi)外學(xué)者先后建立了“軟桿”和“剛桿”模型來研究管柱在井眼中的受力情況，并逐步加以改進(jìn)。對于大位移井，巖屑堆積對井眼中摩阻扭矩的影響不可忽視。在延伸極限方面，高德利［1］提出了大位移井延伸極限的概念，閆鐵［2］、付建紅［3］、黃文君［4-5］等人對大位移井延伸極限做出修正與改進(jìn)。本文利用井眼清潔與摩阻扭矩耦合模型做出大位移井機(jī)械極限延伸預(yù)測。

1 井眼清潔與摩阻扭矩耦合模型

圖1 左表示環(huán)空中存在巖屑時的鉆柱受力情況，其中wbpRdβ 表示鉆柱單元體自重，wcRdβ表示鉆桿單元體與井壁之間的接觸力，μwcRdβ為摩擦阻力F和F+dF表示鉆柱單元體兩端的軸向力；巖屑床對鉆柱受力的具體影響如圖1 右所示，巖屑床對鉆柱，受力的影響復(fù)雜且不可忽略，其中內(nèi)部的力Fp、Frp、Fr、Fcw和Ffw表示巖屑床對鉆桿不同方向的復(fù)雜相互作用。為了模擬巖屑對鉆柱運動的確切影響，需要鉆柱在整個井眼中的位置，通過復(fù)雜的有限元模擬才能獲得。對于實際應(yīng)用，假設(shè)一個簡單的參數(shù)Fcb來表示巖屑床和鉆柱之間的復(fù)雜相互作用，并通過實驗數(shù)據(jù)與現(xiàn)場數(shù)據(jù)以數(shù)據(jù)驅(qū)動模型［6］來估算Fcb［7］。

圖1 存在巖屑床影響的鉆柱單元自由體及局部受力Fig.1 Free body and local stress of drill string unit under the influence of cuttings bed

為了考慮巖屑對井眼中扭矩和阻力的影響，通過修改軟桿模型，提出了一種新的摩阻扭矩模型。將來自巖屑床的附加力添加到標(biāo)準(zhǔn)摩阻扭矩模型中，例如，對于下鉆或滑動鉆進(jìn)的情況，將附加阻力Fcb加到鉆柱的軸向，用來表示巖屑床對鉆柱軸向力的影響。對于旋轉(zhuǎn)鉆進(jìn)的鉆柱，將阻力Fcb加至與鉆柱曲線相切并與旋轉(zhuǎn)方向相反的方向。該模型遵循傳統(tǒng)軟桿模型的假設(shè)，鉆柱被看作一條軟繩，忽略其抗彎強(qiáng)度，且鉆桿與井眼軸線一致并與井壁連續(xù)接觸。

將笛卡爾坐標(biāo)系的單位向量i，j，k轉(zhuǎn)化為弗萊納坐標(biāo)系的單位向量t，n，b，如圖2。

圖2 三維井眼中鉆柱微元自由體Fig.2 Free body of drill string infinitesimal in 3D hole

式中，φ為井斜角，°；θ為方位角，°。

通過對鉆桿元件的自由體進(jìn)行力學(xué)分析，單元體軸向力和扭矩平衡方程可以表示為

式中，F(xiàn)t為軸向力，N；Mt為扭矩，N·m；wbp為鉆柱的單位重量，kg；wc為接觸力，N；tz、nz和bz分別為單位切向量、單位法向量和單位副法向量在z方向上的分量，無量綱；κ為井眼曲率，(°)/30 m；μ為摩阻系數(shù)，無量綱；α為n-b平面中接觸力法向分量與n坐標(biāo)軸之間的夾角，°；Fcb為巖屑床導(dǎo)致的附加力，N；rp為鉆柱半徑，mm。

通過解方程得單元體上的接觸力和n-b平面中接觸力法向分量與n坐標(biāo)軸之間的夾角為

修正后的軟桿模型的算法與標(biāo)準(zhǔn)摩阻扭矩模型相似，不同之處在于，在每個單元上，都需要基于局部鉆屑床高度計算額外的阻力Fcb。從式(4)～(8)可以看出，F(xiàn)cb等效于鉆柱上的附加接觸力。

在將井眼清潔模型與摩阻扭矩模型結(jié)合應(yīng)用時，通過井眼清潔模型預(yù)測巖屑床在井眼中的整體分布。以井眼清潔模型預(yù)測的巖屑床分布作為輸入進(jìn)而改進(jìn)軟桿模型，計算出每個離散點的Fcb，最后得出鉆桿上的摩阻與扭矩值。使用表1 的測試參數(shù)對模型進(jìn)行摩阻實驗，耦合模型的關(guān)鍵是通過研究一系列實驗數(shù)據(jù)，得到Fcb和巖屑床高度的關(guān)系。

表1 測試參數(shù)Table 1 Test parameters

為了建立相關(guān)性，分別使用式(11)和(12)對巖屑床高度和巖屑床高度對扭矩的影響進(jìn)行了歸一化處理。

式中，hb為巖屑床高度，mm；Dw為井眼直徑，mm；h為無量綱巖屑床高度；Th為巖屑床引起的扭矩變化，kN·m；Tc為無巖屑時的扭矩變化，kN·m；TCR表示巖屑床引起的扭矩變化率。

經(jīng)過歸一化處理后，無量綱巖屑床高度與扭矩變化率關(guān)系如圖3 所示，可以發(fā)現(xiàn)二者大致呈指數(shù)形式。因此，在實驗數(shù)據(jù)集的基礎(chǔ)上建立了指數(shù)函數(shù)相關(guān)性，用式(13)表示。

圖3 無量綱巖屑床高度和扭矩變化率的關(guān)系圖Fig.3 Relationship between dimensionless height of cuttings bed and torque variance ratio

其中，c=20，d=2.5。

因為實驗是在模擬直井的設(shè)施上進(jìn)行的，扭矩值與井眼和鉆桿之間的接觸力、井眼和鉆桿之間的摩擦因數(shù)成正比。接觸力在實驗過程中保持不變，扭矩的變化只會導(dǎo)致摩擦因數(shù)的變化。因為參數(shù)是無量綱的，式(13)可轉(zhuǎn)換為無量綱巖屑床高度與井眼和鉆桿之間摩擦因數(shù)的關(guān)系

式中，fCR表示巖屑床引起的摩擦因數(shù)變化率；fh為無量綱巖屑床高度h下的摩擦因數(shù)；fc為相應(yīng)井眼結(jié)構(gòu)下無巖屑床時的摩擦因數(shù)。

在有巖屑的復(fù)雜井中進(jìn)行扭矩計算，式(14)可以應(yīng)用于具有給定局部巖屑床高度的摩阻扭矩模型中，巖屑床高度可從瞬態(tài)巖屑運移模型［8］得到。

與實驗測試不同，在實際鉆井中許多參數(shù)會發(fā)生變化。雖然由實驗數(shù)據(jù)得出的相關(guān)性可能不適用于所有情況，但是摩擦因數(shù)和無量綱巖屑床高度的變化趨勢是相同的。即，式(11)～(14)仍然有效，但c和d的值可能在不同的情況下有所不同。為了讓模型的應(yīng)用更加廣闊，可以采用機(jī)器學(xué)習(xí)［9］的方法，根據(jù)歷史鉆井?dāng)?shù)據(jù)對模型進(jìn)行訓(xùn)練后計算。

2 水平段機(jī)械延伸極限預(yù)測

大位移井水平段延伸極限是指在一定制約條件影響下，大位移井所能達(dá)到的極限長度。由于受到多種影響因素的制約，大位移井水平段不可能無限延伸，同時，根據(jù)相關(guān)研究，從經(jīng)濟(jì)角度來講，水平段長度也不是越長越好。大位移井水平段延伸能力影響到鉆探該類型井的施工效率以及全井成本，因此，對影響大位移井水平段延伸能力因素的分析，在大位移井鉆井工程設(shè)計中是十分重要的。制約大位移井水平段機(jī)械延伸極限的因素主要有鉆機(jī)性能約束、鉆柱屈服約束和鉆頭運動約束［5］。

2.1 鉆機(jī)性能約束

鉆機(jī)的承載能力包括提升能力和扭矩輸出能力。提升能力指在鉆井作業(yè)中，起鉆、下鉆、鉆進(jìn)等工況下井口大鉤載荷必須在鉆機(jī)的額定載荷范圍內(nèi)；扭矩輸出能力是指在鉆柱旋轉(zhuǎn)的工況下，井口扭矩不能超過鉆機(jī)能夠輸出的最大扭矩。為保證鉆機(jī)在水平段鉆進(jìn)過程中能夠正常工作，需要準(zhǔn)確預(yù)測不同工況下的井口大鉤載荷和井口扭矩。

使用井眼清潔與摩阻扭矩耦合模型，考慮了巖屑對井眼中扭矩和阻力的影響，更加精準(zhǔn)地模擬大位移井鉆柱的真實受力。通過式(4)～(8)計算出井口的軸向力和扭矩。

2.2 鉆柱屈服約束

根據(jù)井眼清潔與摩阻扭矩耦合模型計算出沿井深分布的軸向力和扭矩，然后計算出相應(yīng)的軸向應(yīng)力和切應(yīng)力，根據(jù)第四強(qiáng)度理論計算鉆柱微元上的等效應(yīng)力 σ，為確保鉆柱不發(fā)生屈服失效，鉆柱等效應(yīng)力應(yīng)小于某值，其約束條件為

式中，σ為鉆柱等效應(yīng)力，Pa；[σ]為許用等效應(yīng)力，P a；σs為鉆柱的屈服強(qiáng)度，Pa；n4是為安全系數(shù)，無因次。

2.3 鉆頭運動約束

鉆柱在滑動上提與旋轉(zhuǎn)上提的過程中，有可能發(fā)生卡鉆。在鉆頭處添加一個上提阻力來進(jìn)行卡鉆情況模擬。所以鉆頭位置的軸向力與扭矩必須等于模擬卡鉆阻力與扭矩，而且鉆頭處鉆壓與扭矩必須大于巖石的門限鉆壓與門限扭矩。

3 實例分析

3.1 實例井簡況

某井是一口水垂比大于5 的大位移井，井身結(jié)構(gòu)為四開井，如圖4 所示，井眼軌跡數(shù)據(jù)如表2 所示。該井平臺井架提升能力為3 629.75 kN，頂驅(qū)提升能力為7 352.91 kN，頂驅(qū)最大連續(xù)輸出扭矩為88.81 kN·m，旋轉(zhuǎn)速度為100 r/min。

圖4 井身結(jié)構(gòu)Fig.4 Casing program

表2 井眼軌跡數(shù)據(jù)Table 2 Data of well trajectory

3.2 摩阻扭矩

摩阻扭矩計算參數(shù)如下：裸眼摩阻系數(shù)0.3，套管摩阻系數(shù)0.25，鉆井液密度1.03 g/cm3，鉆具活動速度0.3 m/s，旋轉(zhuǎn)速度100 r/min，鉆進(jìn)鉆壓44 kN，鉆頭扭矩2 711 N·m。假定水平段巖屑床高度是40 mm 和80 mm，使用井眼清潔與摩阻扭矩耦合模型計算不同工況下有效軸向力與扭矩分布。

圖5 為使用耦合模型分別計算水平段巖屑床高度40 mm 和80 mm 不同工況下的鉆柱有效軸向力分布，結(jié)果表明，起鉆工況下，鉆柱受拉影響最大，但在巖屑床高度40 mm 的情況下未達(dá)到鉆柱拉力極限，當(dāng)巖屑床高度提升至80 mm 時，鉆柱在起鉆工況下達(dá)到拉力極限，發(fā)生斷裂；下鉆及滑動鉆進(jìn)工況下，鉆柱受壓，在巖屑床高度80 mm 情況下，鉆柱在井口～400 m 時達(dá)到螺旋屈曲臨界值，發(fā)生螺旋屈曲，軸向壓力陡增。圖6 為使用耦合模型分別計算巖屑床高度40 mm 和80 mm 不同工況下的鉆柱扭矩分布，結(jié)果表明，在巖屑床高度80 mm 時，旋轉(zhuǎn)鉆進(jìn)與倒劃眼均達(dá)到扭矩極限，且在旋轉(zhuǎn)鉆進(jìn)工況下，鉆柱在井口～993 m 達(dá)到扭矩極限。

圖5 不同工況下有效軸向力隨鉆柱分布Fig.5 Distribution of effective axial force along drill string under different working conditions

圖6 不同工況下扭矩隨鉆柱分布Fig.6 Distribution of torque along drill string under different working conditions

如圖7 所示，在滑動上提工況下，隨著巖屑床的堆積，即巖屑床高度的提高，鉤載隨井深增大。如圖8 所示，在滑動鉆進(jìn)工況下，隨著巖屑床高度的提高，鉤載隨井深減小，但當(dāng)井深增加到一定深度后，鉤載迅速減小，造成此現(xiàn)象的原因是鉆柱發(fā)生屈曲產(chǎn)生的高摩擦阻力，當(dāng)鉤載為0 時的井深為延伸極限。通過對比分析，滑動鉆進(jìn)是滑動工況下更為限制延伸極限的工況，以巖屑床高度為40 mm 為標(biāo)準(zhǔn)計算，機(jī)械延伸極限為6 156 m。

圖7 滑動上提大鉤載荷與井深關(guān)系Fig.7 Relationship between hook load and well depth during slide lifting

圖8 滑動鉆進(jìn)大鉤載荷與井深關(guān)系Fig.8 Relationship between hook load and well depth during slide drilling

如圖9 所示，在旋轉(zhuǎn)上提工況下，隨著巖屑床的堆積，即巖屑床高度的提高，井口扭矩隨井深增大。如圖10 所示，在旋轉(zhuǎn)鉆進(jìn)工況下，隨著巖屑床高度的提高，井口扭矩隨井深增大，當(dāng)井口扭矩達(dá)到額定扭矩時的井深為延伸極限。通過對比分析，旋轉(zhuǎn)鉆進(jìn)是旋轉(zhuǎn)工況下更為限制延伸極限的工況，以巖屑床為40 mm 為標(biāo)準(zhǔn)計算，機(jī)械延伸極限為9 035 m。

圖9 旋轉(zhuǎn)上提井口扭矩與井深關(guān)系Fig.9 Relationship between wellhead torque and well depth during rotary lifting

圖10 旋轉(zhuǎn)鉆進(jìn)井口扭矩與井深關(guān)系Fig.10 Relationship between wellhead torque and well depth during rotary drilling

3.3 大位移井機(jī)械延伸極限預(yù)測

對于滑動工況來說，機(jī)械延伸極限指起鉆和滑動鉆進(jìn)延伸極限的最小值；對于旋轉(zhuǎn)工況來說，機(jī)械延伸極限指倒劃眼和旋轉(zhuǎn)鉆進(jìn)延伸極限的最小值。

在相同鉆機(jī)性能、僅改變巖屑床高度的情況下，使用上述鉆進(jìn)參數(shù)計算不同鉆進(jìn)方式的機(jī)械延伸極限，如圖11 所示。結(jié)果表明，隨著巖屑床高度的增加，機(jī)械延伸極限相應(yīng)降低，而且隨著巖屑堆積到一定程度后，對延伸極限的影響顯著。巖屑床高度的增加對不同鉆進(jìn)方式的機(jī)械延伸極限的變化趨勢相同，但滑動鉆進(jìn)的延伸極限明顯低于旋轉(zhuǎn)鉆進(jìn)。

圖11 不同巖屑床高度對機(jī)械延伸極限的影響Fig.11 Influence of the height of cuttings bed on the mechanical extension limit

4 結(jié)論

(1)考慮巖屑床對井筒中摩阻扭矩的影響，通過對前人的軟桿模型進(jìn)行修正，建立井眼清潔與摩阻扭矩耦合模型。提出大位移井水平段機(jī)械延伸極限的約束條件，并建立預(yù)測模型。

(2)使用井眼清潔與摩阻扭矩耦合模型進(jìn)行實鉆井摩阻扭矩分析，分別計算水平段巖屑床高度40 mm 和80 mm 下鉆柱有效軸向力和扭矩分布，校核鉆柱強(qiáng)度。巖屑堆積到一定程度后，鉆柱屈服約束制約大位移井機(jī)械延伸極限；計算不同水平段巖屑床高度與不同工況下大鉤載荷與井口扭矩，巖屑堆積對鉤載與井口扭矩影響顯著；分析不同巖屑床高度對機(jī)械延伸極限的影響，在無法提升鉆機(jī)性能情況下，提高井眼清潔效率可有效提高延伸極限。

(3)井眼清潔與摩阻扭矩耦合模型仍需要根據(jù)歷史鉆井?dāng)?shù)據(jù)通過機(jī)器學(xué)習(xí)進(jìn)行訓(xùn)練；大位移井延伸極限也受到鉆井水力學(xué)［10］的限制，而且井眼清潔會緩解巖屑堆積。本文僅討論巖屑床對機(jī)械延伸極限的影響，實際使用應(yīng)綜合分析。