黃文君 王 舸 高德利

中國石油大學(xué)（北京）石油工程教育部重點(diǎn)實(shí)驗(yàn)室

0 引言

旋轉(zhuǎn)導(dǎo)向系統(tǒng)是指在鉆柱旋轉(zhuǎn)過程中同時(shí)實(shí)現(xiàn)導(dǎo)向鉆進(jìn)功能的一種鉆井系統(tǒng)，可顯著地降低摩阻、提高鉆速、改善井眼條件并大幅度提高井眼延伸長度[1-11]，代表著現(xiàn)代導(dǎo)向鉆井技術(shù)的發(fā)展方向。

旋轉(zhuǎn)導(dǎo)向工具造斜率的準(zhǔn)確預(yù)測(cè)是導(dǎo)向工具優(yōu)化設(shè)計(jì)與井眼軌跡精確控制的前提。相關(guān)學(xué)者針對(duì)造斜率的預(yù)測(cè)問題開展了大量研究，主要經(jīng)歷了從幾何法、力學(xué)法到軌跡法的過程。幾何法主要是指“三點(diǎn)定圓法”[12-14]，力學(xué)法包括“極限曲率法”[15-16]“平衡側(cè)向力法”[17]等，軌跡法包括基于井斜趨勢(shì)角的計(jì)算方法[18-21]等。幾何法忽略了鉆具的剛度以及鉆井參數(shù)的影響，力學(xué)法忽略了鉆頭切削性能和地層各向異性的影響，軌跡法充分綜合考慮了各種因素，但軌跡結(jié)果與造斜率之間缺少明確的定量關(guān)系，不利于造斜率的準(zhǔn)確評(píng)估。

筆者在前人研究基礎(chǔ)上，提出了推靠式旋轉(zhuǎn)導(dǎo)向鉆進(jìn)“極限造斜率”和“折算造斜率”的計(jì)算方法，從“杠桿效應(yīng)”“鐘擺效應(yīng)”和“推靠效應(yīng)”角度揭示了造斜率的影響機(jī)理，分析了鉆具結(jié)構(gòu)、鉆進(jìn)參數(shù)、鉆頭與地層各向異性等參數(shù)對(duì)推靠式旋轉(zhuǎn)導(dǎo)向工具造斜率的影響規(guī)律，以期為推靠式旋轉(zhuǎn)導(dǎo)向工具造斜率預(yù)測(cè)及鉆具結(jié)構(gòu)優(yōu)化提供理論基礎(chǔ)。

1 軌跡預(yù)測(cè)模型

1.1 旋轉(zhuǎn)導(dǎo)向鉆具受力模型

推靠式旋轉(zhuǎn)導(dǎo)向鉆具可看成一系列縱橫彎曲梁的組合，其變形的控制方程如下：

式中k=1代表井斜平面，k=2代表方位平面；uk表示管柱在井斜平面和方位平面上的橫向位移，m；s表示管柱任意一點(diǎn)的弧長，m；EI表示相應(yīng)的抗彎剛度，N·m2；F表示軸向力，kN；KB表示管柱的浮力系數(shù)，KB=1-ρd/ρs（ρd表示鉆井液密度，g/cm3；ρs表示管柱的材料密度，g/cm3）；qk表示管柱線重投影到井斜平面和方位平面的分量，kN/m。

以鉆頭、偏置機(jī)構(gòu)、各跨穩(wěn)定器、管柱與井壁的可能接觸位置為節(jié)點(diǎn)，根據(jù)接觸情況，對(duì)導(dǎo)向鉆具進(jìn)行單元?jiǎng)澐?，如圖1所示。圖1中Fp表示推靠力，kN；Wb表示鉆壓，kN；Fb,k表示鉆頭側(cè)向力在井斜和方位面上的分量，kN；TS表示穩(wěn)定器接觸井壁時(shí)的支反力，kN；TC表示管柱上切點(diǎn)的支反力，kN；FC表示管柱上切點(diǎn)的軸向力，kN。

利用加權(quán)余量法，令導(dǎo)向鉆具第j跨的試撓度函數(shù)為：

式中Ck,i,j表示待定系數(shù)；N表示導(dǎo)向鉆具上節(jié)點(diǎn)個(gè)數(shù)。

鉆頭處的位移和彎矩均為0，即

上切點(diǎn)處的邊界條件如下：

若節(jié)點(diǎn)與井壁接觸，則節(jié)點(diǎn)處位移、轉(zhuǎn)角、彎矩的連續(xù)條件如下：

式中Lk,j表示第j跨的長度，m；mk,j表示第j節(jié)點(diǎn)無因次系數(shù)，取決于管柱與井壁的接觸位置；rj表示第j跨處的井眼視半徑，m；yk,j表示鉆頭至當(dāng)前節(jié)點(diǎn)的狗腿角在井斜和方位平面上的投影，rad。

若單元節(jié)點(diǎn)為推靠力位置或者該點(diǎn)處管柱與井壁不接觸，則位移、轉(zhuǎn)角、彎矩、剪力的連續(xù)條件如下：

式中第4式右端，如果j節(jié)點(diǎn)為推靠位置，則Fk,p為推靠力分別投影在井斜和方位平面上的分量；如果j節(jié)點(diǎn)為懸空位置，則取0 kN。

將式（2）代入式（1）中即可得到殘值，利用子域法得殘差方程：

式中εk,j表示內(nèi)部殘值。

對(duì)于井斜平面，根據(jù)試撓度函數(shù)式（2）可知每一跨變形有5個(gè)待定參數(shù)，上切點(diǎn)至穩(wěn)定器的長度是1個(gè)待定參數(shù)，因此共有5N+6個(gè)未知參數(shù)；根據(jù)式（3）～（7）可知上下端有5個(gè)邊界方程，各節(jié)點(diǎn)有4個(gè)連續(xù)性方程，每一跨有1個(gè)殘差方程，共有5N+6個(gè)已知方程。對(duì)于方位平面而言，未知參數(shù)和已知方程個(gè)數(shù)相同，即模型可解。該模型可利用二分法粗選解的存在區(qū)間，再用牛頓迭代法進(jìn)行精細(xì)搜索至最終解。

利用求解出的系數(shù)Ck,i,j，可得到鉆頭上側(cè)向力與轉(zhuǎn)角的計(jì)算結(jié)果：

式中θb,k表示鉆頭轉(zhuǎn)角在井斜和方位平面上的分量，rad。

1.2 鉆頭—地層相互作用模型

在井眼軌跡預(yù)測(cè)中，除了考慮導(dǎo)向工具的受力和變形，還需要考慮鉆頭與地層的各向異性[22]，具體關(guān)系如圖2所示。圖2中α表示井斜角，rad；θb表示井眼軸線與鉆頭軸線的夾角，rad；β表示鉆進(jìn)方向偏角，rad；R表示鉆速，m/h；Fz表示鉆頭處的軸向力，kN；γ為地層傾角，rad。

將式（8）代入鉆頭與地層相互作用模型中，即可得到鉆速預(yù)測(cè)方程。鉆頭與地層相互作用模型可表示為如下形式：

式中Ra、R1和R2表示考慮鉆頭地層各向異性指數(shù)下的鉆速分量，m/h；Dn表示標(biāo)準(zhǔn)（法向）鉆井效率；R表示井底—地層參考坐標(biāo)系的轉(zhuǎn)換矩陣；S表示井底—鉆頭參考坐標(biāo)系的轉(zhuǎn)換矩陣；Ir1、Ir2表示地層的各向異性指數(shù)，無因次；Ib表示鉆頭的各向異性指數(shù)，無因次。

2 造斜率計(jì)算方法

2.1 造斜率的概念

前人在計(jì)算導(dǎo)向鉆具的造斜率時(shí)，認(rèn)為鉆頭側(cè)向力為零時(shí)對(duì)應(yīng)的井眼曲率為鉆具的極限造斜率，但忽略了鉆頭各向異性以及地層各向異性的影響。鉆頭與地層各向異性與鉆進(jìn)軌跡密切相關(guān)，進(jìn)而對(duì)造斜率的影響不可忽略，因此需要計(jì)入這兩因素的影響?；谇叭搜芯?，筆者提出了兩種造斜率的概念：①極限造斜率，側(cè)向鉆速為零對(duì)應(yīng)的井眼曲率。②折算造斜率，考慮折算系數(shù)修正后的造斜率。

極限造斜率是指鉆具造斜所到達(dá)的極限狀態(tài)，是基于“平衡曲率”“極限曲率”的概念發(fā)展而言。極限造斜率與井斜角、井眼尺寸、鉆具組合結(jié)構(gòu)、鉆頭各向異性、鉆壓、地層各向異性、地層走向和傾向等參數(shù)密切相關(guān)。

造斜率除了與地層、鉆具等因素相關(guān)外，還與鉆進(jìn)過程密切相關(guān)。在達(dá)到極限造斜率之前需要一個(gè)較長過程，然而在此過程中井斜角等參數(shù)也不斷變化，使得鉆進(jìn)中的造斜率一般會(huì)小于極限造斜率，因此需要利用折算造斜率來預(yù)測(cè)實(shí)際造斜率。

2.2 造斜率的計(jì)算方法

根據(jù)極限造斜率的定義，其計(jì)算模型如下：

方程（10）可采用牛頓迭代法求解井斜平面和方位平面上的造斜率κ1和κ2，可得到極限造斜率的計(jì)算公式：

式中κL表示極限造斜率的理論計(jì)算值，m－1。

折算造斜率無法從理論上直接求解，需要借助實(shí)鉆資料進(jìn)行反演計(jì)算。但是，實(shí)鉆資料只包含特定鉆具在特定條件下的軌跡數(shù)據(jù)，難以揭示各個(gè)因素對(duì)造斜率的影響規(guī)律，無法指導(dǎo)鉆具組合的優(yōu)化設(shè)計(jì)。結(jié)合理論分析以及實(shí)踐經(jīng)驗(yàn)，在極限造斜率基礎(chǔ)上引入折算系數(shù)，可得到折算造斜率的近似計(jì)算式：

式中κA表示折算造斜率，m－1；λ表示折算系數(shù)，需要利用實(shí)鉆數(shù)據(jù)反演求解。

3 實(shí)例分析

常規(guī)導(dǎo)向鉆具軌跡控制機(jī)理主要包括杠桿效應(yīng)、鐘擺效應(yīng)、滿眼效應(yīng)等，杠桿效應(yīng)具有增斜效果，鐘擺效應(yīng)具有降斜效果，滿眼效應(yīng)具有穩(wěn)斜效果。導(dǎo)向工具的造斜率取決于各種效應(yīng)的相對(duì)大小。對(duì)于推靠式旋轉(zhuǎn)導(dǎo)向系統(tǒng)而言，除了上述幾種效應(yīng)外，還包括推靠塊接觸井壁產(chǎn)生的推靠效應(yīng)。推靠效應(yīng)使得井眼軌跡沿著推靠力作用方向變化，當(dāng)推靠力在井筒截面上的投影指向上時(shí)具有增斜效果，指向下時(shí)具有降斜效果。下面結(jié)合計(jì)算實(shí)例，分析各種效應(yīng)對(duì)造斜率的影響規(guī)律。

采用一種旋轉(zhuǎn)導(dǎo)向鉆具設(shè)計(jì)方案進(jìn)行分析，如圖3所示。其具體結(jié)構(gòu)為：?215.9 mmPDC鉆頭+?189 mm旋轉(zhuǎn)導(dǎo)向鉆具+?212.7 mm導(dǎo)向穩(wěn)定器+?127 mm柔性短節(jié)+?178 mm鉆鋌+?212.7 mm上穩(wěn)定器+?178 mm鉆鋌。鉆頭各向異性指數(shù)為0.3，平均鉆壓為80 kN。

圖4-a為試驗(yàn)井鉆進(jìn)過程中推靠塊上推靠合力的測(cè)量結(jié)果，圖4-b為利用造斜率模型得到的不同折算系數(shù)下的結(jié)果，實(shí)測(cè)造斜率是根據(jù)實(shí)鉆軌跡數(shù)據(jù)并采用最小曲率插值法得到。結(jié)果表明，實(shí)測(cè)結(jié)果與計(jì)算結(jié)果比較吻合，折算系數(shù)主要處于0.65～0.85之間，平均值為0.75，驗(yàn)證了本文模型的有效性。

表1為不同推靠力和井斜角下的造斜率，其中推靠力分別為0 kN、10 kN、20 kN和30 kN。當(dāng)推靠塊全部縮進(jìn)時(shí)（Fp=0 kN），導(dǎo)向鉆具呈現(xiàn)降斜特點(diǎn)，因?yàn)殓姅[效應(yīng)大于杠桿效應(yīng)。

表1 不同推靠力作用下造斜率隨井斜角的變化規(guī)律表

隨著井斜角增大，鐘擺效應(yīng)增大，導(dǎo)致降斜率不斷增大。當(dāng)推靠塊與井壁產(chǎn)生指向上方的推靠力時(shí)，推靠效應(yīng)超過鐘擺效應(yīng)，造斜率由負(fù)值變成正值。隨著推靠力的增大，推靠效應(yīng)不斷增大，使得鉆進(jìn)造斜率也不斷增大。造斜率增大幅度與推靠力的大小近似呈正比例關(guān)系，表明推靠效應(yīng)明顯大于其他效應(yīng)，在造斜率產(chǎn)生中起到?jīng)Q定性作用。當(dāng)推靠力不為零時(shí)，隨著井斜角的增大，造斜率呈現(xiàn)先增大后降低的趨勢(shì)。對(duì)于造斜率增大階段，主要是杠桿效應(yīng)、鐘擺效應(yīng)以及推靠效應(yīng)的一種綜合結(jié)果；對(duì)于造斜率降低階段，是由于鐘擺效應(yīng)增大造成的。

表2為井斜角為30°條件下不同推靠力和鉆壓下的造斜率。當(dāng)推靠塊全部縮進(jìn)時(shí)（Fp=0 kN），隨著鉆壓的增大，杠桿效應(yīng)增大，使得造斜率略微增大；當(dāng)具有增斜效應(yīng)的推靠力存在時(shí)，隨著鉆壓增大，造斜率反而降低。這是因?yàn)殡S著鉆壓增大，推靠力與鉆壓之比減小，推靠力影響鉆進(jìn)軌跡的相對(duì)作用減小，即推靠效應(yīng)降低，因而造斜率呈現(xiàn)緩降趨勢(shì)。

表2 不同推靠力作用下造斜率隨鉆壓的變化規(guī)律表

表3為推靠力為20 kN條件下不同柔性短節(jié)長度和鉆壓下的造斜率。隨著柔性短節(jié)長度的增大，造斜率也不斷增大。由于柔性短節(jié)的剛度明顯小于相鄰管柱的剛度，因此可降低柔性短節(jié)上部鉆具剛度對(duì)下部鉆具造斜能力的影響。隨著柔性短節(jié)長度增大，上部鉆具剛度的影響作用減小，更有利于推靠效應(yīng)的發(fā)揮。當(dāng)柔性短節(jié)長度較長時(shí)（大于2 m），再提高柔性短節(jié)長度時(shí)上部鉆具剛度的影響變?nèi)?，此時(shí)造斜率增大的幅度將明顯降低。

表3 不同柔性短節(jié)長度下造斜率隨鉆壓的變化規(guī)律表

表4為鉆頭各向異性指數(shù)（Ib）對(duì)造斜率的影響。隨著Ib的增大，鉆頭的切削能力不斷提高。當(dāng)Ib=1時(shí)為各向同性鉆頭，側(cè)向切削能力與軸向切削能力相同；當(dāng)Ib=0時(shí)鉆頭只具有軸向切削能力，無側(cè)向切削能力。推靠效應(yīng)的機(jī)理是提高鉆頭的側(cè)向力，進(jìn)而提高軌跡變化能力。軌跡變化能力與鉆頭側(cè)向力的關(guān)系主要取決于鉆頭側(cè)向切削能力，即Ib。因此，隨著Ib的增大，推靠效應(yīng)增大，使得造斜率也不斷提高。

表4 造斜率隨鉆頭各向異性指數(shù)的變化規(guī)律表

4 結(jié)論

1）極限造斜率是一種理想結(jié)果，折算造斜率一般小于極限造斜率，它們之間存在一個(gè)近似的比例關(guān)系，該比例系數(shù)與鉆進(jìn)過程密切相關(guān)，需要借助實(shí)鉆數(shù)據(jù)反演來確定該系數(shù)。

2）推靠式旋轉(zhuǎn)導(dǎo)向鉆具的造斜率主要取決于杠桿效應(yīng)、鐘擺效應(yīng)、推靠效應(yīng)三者的綜合結(jié)果。推靠效應(yīng)對(duì)鉆具造斜率的影響最顯著，因此鉆具造斜率與推靠力的相關(guān)性最高。杠桿效應(yīng)和鐘擺效應(yīng)起到次要作用，因此井斜角和鉆壓對(duì)造斜率的影響相對(duì)較小。

3）對(duì)于常規(guī)滑動(dòng)導(dǎo)向鉆具而言，隨著鉆壓增大，杠桿效應(yīng)增大，造斜率也增大；對(duì)于推靠式旋轉(zhuǎn)導(dǎo)向鉆具而言，隨著鉆壓增大，杠桿效應(yīng)所占比例增大，推靠效應(yīng)所占比例減小，導(dǎo)致造斜率反而減小。推靠效應(yīng)的存在，使得推靠式旋轉(zhuǎn)導(dǎo)向鉆具相對(duì)于常規(guī)滑動(dòng)導(dǎo)向鉆具而言，其造斜機(jī)理更加復(fù)雜，甚至出現(xiàn)相反的影響規(guī)律。

4）造斜率與鉆具結(jié)構(gòu)參數(shù)（鉆柱尺寸和長度、柔性短節(jié)等）、推靠力（大小、方向）、鉆頭性能（各向異性指數(shù)等）、鉆進(jìn)參數(shù)（鉆壓等）等因素密切相關(guān)。為了更加精確地控制井眼軌跡，一方面在設(shè)計(jì)階段要合理地優(yōu)選鉆具組合和匹配的鉆頭，另一方面在實(shí)鉆過程中要合理地控制鉆井參數(shù)和推靠力大小。