袁蘭峰，何世明，黎學(xué)年，湯明，邵庸

（1.西南石油大學(xué)油氣藏地質(zhì)及開發(fā)工程國家重點實驗室，四川 成都 610050；2.中國石化華東石油工程公司鉆井處，江蘇 南京 210000）

考慮鉆具屈曲的環(huán)空流動仿真研究

袁蘭峰1，何世明1，黎學(xué)年2，湯明1，邵庸1

（1.西南石油大學(xué)油氣藏地質(zhì)及開發(fā)工程國家重點實驗室，四川 成都 610050；2.中國石化華東石油工程公司鉆井處，江蘇 南京 210000）

在小井眼定向井和水平井鉆井過程中，鉆具屈曲會對環(huán)空流場、壓力場產(chǎn)生較大影響。文中借助流體動力學(xué)計算軟件ANSYS FLUENT開展了考慮鉆具正弦、螺旋屈曲的環(huán)空流動仿真模擬，系統(tǒng)地研究了不同環(huán)空偏心距、不同轉(zhuǎn)速等因素對環(huán)空流場和流動壓降的影響規(guī)律。研究發(fā)現(xiàn)：同心環(huán)空中，流體流速沿徑向?qū)ΨQ分布，偏心、屈曲環(huán)空中偏離對稱分布；環(huán)空間隙越大，環(huán)空流速就越大；當(dāng)考慮鉆柱旋轉(zhuǎn)時，鉆柱旋轉(zhuǎn)對環(huán)空窄間隙處流速的影響比對環(huán)空寬間隙處的影響要大，可明顯提高環(huán)空窄間隙處的流速；同心環(huán)空、偏心環(huán)空壓降值沿井深呈線性增加，偏心度越大，壓降越?。磺h(huán)空壓降值沿井深呈非線性增加，轉(zhuǎn)速越大，壓降越小。仿真結(jié)果與實驗結(jié)果相比，具有較高的準(zhǔn)確性，對真實鉆井條件下環(huán)空壓降的計算具有一定的指導(dǎo)意義。

偏心環(huán)空；正弦屈曲；螺旋屈曲；壓降；ANSYS FLUENT軟件

0 引言

隨著油氣開采向著深部地層及深水方向發(fā)展，鉆遇窄安全密度窗口地層的情況越來越普遍，環(huán)空壓降的預(yù)測不準(zhǔn)常會造成井漏、井噴等安全事故，環(huán)空結(jié)構(gòu)作為影響環(huán)空壓降的一個重要因素，引起了重點關(guān)注。

流體在環(huán)空中的流動相較于在直管中要復(fù)雜得多，早期的油田實踐和室內(nèi)實驗對流體在環(huán)空中的流動規(guī)律已經(jīng)作了許多研究，其中對牛頓流體［1-2］、非牛頓流體［3-6］在層流、紊流流態(tài)下的流動規(guī)律主要是應(yīng)用實驗、理論解析以及數(shù)值分析方法進行了求解，并得到了許多有價值的結(jié)論。但這些研究主要是以同心、偏心環(huán)空結(jié)構(gòu)為主，對屈曲環(huán)空中的速度場、壓力場分布規(guī)律研究較少，然而在實際鉆井過程中，所鉆井多為大斜度井或水平井，鉆柱由于受到重力的作用總是靠向環(huán)空低邊發(fā)生偏心現(xiàn)象，同時當(dāng)鉆柱軸向載荷超過其臨界屈曲載荷時就會發(fā)生彎曲［7］，所以，實際環(huán)空多為偏心與屈曲環(huán)空相結(jié)合，而非單一的同心或偏心結(jié)構(gòu)。因此，傳統(tǒng)的水力學(xué)模型并不能完全、準(zhǔn)確預(yù)測實際鉆井條件下的壓降。

現(xiàn)階段國內(nèi)外關(guān)于內(nèi)管旋轉(zhuǎn)條件下屈曲環(huán)空流場、壓力場研究較少，為此，本文基于CFD仿真實驗，對考慮內(nèi)管旋轉(zhuǎn)條件下的屈曲環(huán)空流場和壓力場進行了數(shù)值模擬研究，得出了不同偏心距、不同轉(zhuǎn)速時的流場和壓力場的變化規(guī)律。

1 研究問題及模型

1.1 屈曲環(huán)空的形成

在定向井及水平井的鉆井過程中，由于鉆柱受到軸向壓力作用，一旦軸向載荷超過發(fā)生屈曲的臨界載荷，鉆柱就會發(fā)生屈曲。確定鉆柱發(fā)生螺旋屈曲時的臨界載荷方程為

式中：Fhel為螺旋屈曲臨界載荷，N；E為彈性模量，Pa；I為轉(zhuǎn)動慣量，kg·m2；wc為定義常數(shù)；rc為鉆柱與井眼間的間距，m；θ為井斜角，（°）；ρt為管材的單位長度的質(zhì)量，kg/m；As為管柱的橫截面積，m2；g為重力加速度，m/s2。

Miska等［8-9］提出了在水平井段鉆柱發(fā)生正弦屈曲時的臨界載荷方程：

式中：Fsin為正弦屈曲臨界載荷，N。

鉆柱在軸向壓力超過其臨界載荷時，會發(fā)生正弦屈曲或螺旋屈曲，此時，由屈曲鉆柱與井壁或套管壁形成的環(huán)空稱之為屈曲環(huán)空，研究鉆柱屈曲條件下的環(huán)空壓力更切合實際情況。

1.2 CFD仿真模型及輸入?yún)?shù)

1.2.1 控制方程及邊界條件

對于不可壓縮、等溫且處于層流流態(tài)下的流體來說，其有效黏度僅僅與應(yīng)變率張量有關(guān)。液體流動模型可使用圓柱坐標(biāo)系下的連續(xù)性方程以及軸向、徑向、周向方向的動量方程進行描述。

仿真時的邊界條件包括速度入口、流量出口及固體壁面等。其中：入口軸向速度u設(shè)定為0.83 m/s，徑向速度v和切向速度w均為0；當(dāng)鉆柱旋轉(zhuǎn)時，在井壁處u=v=w=0，在鉆柱壁面處，u=w=0，v=wRi。Ri為鉆桿外徑。

1.2.2 網(wǎng)格劃分及物理模型的建立

為了研究非牛頓流體在不同環(huán)空結(jié)構(gòu)下的流場、壓力場特性，選取10 m長的環(huán)空井段，井徑為0.1524 m，鉆桿外徑為0.0889 m，環(huán)空結(jié)構(gòu)分為同心、偏心（偏心度分別為0.50，0.93）、屈曲環(huán)空（螺旋屈曲、正弦屈曲，其屈曲間距皆為5 m）3種。其中屈曲環(huán)空結(jié)構(gòu)如圖1所示，用ANSYS-Workbench中自帶的mesh模塊對網(wǎng)格進行劃分。鉆井液流體設(shè)定為屈服冪律流體，其流變方程為

式（3）中：n為流性指數(shù)，本次實驗中為 0.416；K為稠度系數(shù)，0.591 Pa·sn；τ0為屈服壓力，0.846 Pa；ρ為密度，1240 kg/m3。

圖1 正弦屈曲幾何模型及網(wǎng)格劃分

采用基于交錯網(wǎng)格的SIMPLE法對連續(xù)性方程和N-S方程進行離散化。對動量方程的對流項進行離散化處理時，為保證計算精度，采用二階格式［10］，設(shè)置速度入口邊界（0.83 m/s）、流量出口邊界、壁面邊界（其中環(huán)空外壁面始終靜止，環(huán)空內(nèi)壁面可旋轉(zhuǎn)，轉(zhuǎn)速為40，60，80 r/min）3類邊界類型。為了計算準(zhǔn)確，對靠近壁面處網(wǎng)格進行加密處理，坐標(biāo)系為柱坐標(biāo)系，速度入口采用速度分量方式，旋轉(zhuǎn)壁面邊界為無滑移壁面。

2 仿真結(jié)果與討論

2.1 內(nèi)管不旋轉(zhuǎn)時徑向流速分布

對于同心與偏心環(huán)空，取距入口5.00 m處的橫截面進行流場分布分析，對于正弦屈曲環(huán)空結(jié)構(gòu)，因其不同截面處的流速分布各不相同；故分別取距入口2.50，3.75，5.00 m的3個截面（即管柱分別位于環(huán)空高邊、環(huán)空中央、環(huán)空低邊時的3個截面）進行流速分布規(guī)律分析。同理，對于螺旋屈曲環(huán)空結(jié)構(gòu)，分別取距入口0.50，5.00，9.50 m的3個截面 （因管柱發(fā)生螺旋屈曲時始終與井壁接觸，故所取截面均為偏心環(huán)空截面）進行流速分布規(guī)律分析。鉆柱靜止時不同環(huán)空結(jié)構(gòu)下流場分布如圖2—4所示。

圖2 鉆柱靜止、管柱未發(fā)生屈曲時不同環(huán)空結(jié)構(gòu)速度云圖

圖3 鉆柱靜止、管柱發(fā)生正弦屈曲時環(huán)空結(jié)構(gòu)速度云圖

由圖2可知：在同心環(huán)空中，流速沿徑向?qū)ΨQ分布，最大流速位于環(huán)空中心區(qū)域，離壁面越近，流速越小。這是因為：靠近壁面時，流體要受到壁面摩擦力的作用，從而產(chǎn)生速度梯度；偏心環(huán)空時，流速偏離對稱分布，且偏心度越大，偏離越明顯，隨著偏心度的增加，窄環(huán)空間隙處速度減小，寬環(huán)空間隙處速度增大，這是由于流體在窄環(huán)空間隙處的流動受到了限制，進而大量流向?qū)挱h(huán)空間隙處的結(jié)果。對于定向井及水平井，窄環(huán)空間隙處的低流速可能會造成巖屑顆粒的沉積，最終導(dǎo)致巖屑床的形成。

圖4 鉆柱靜止、管柱發(fā)生螺旋屈曲時環(huán)空結(jié)構(gòu)速度云圖

由圖3可知，正弦屈曲環(huán)空是一個偏心度不斷變化的環(huán)空結(jié)構(gòu)，在一個屈曲間距內(nèi)，其環(huán)空經(jīng)歷了偏心—同心—偏心的變化過程，環(huán)空中流速偏離對稱分布，環(huán)空間隙越大，流速越大。

圖4為螺旋屈曲環(huán)空結(jié)構(gòu)速度云圖。由圖可知，螺旋屈曲環(huán)空窄間隙處存在一個流動靜止區(qū)，在該區(qū)域軸向流速接近為0。這與完全偏心環(huán)空時觀察到的情況一致，最大流速出現(xiàn)在環(huán)空間隙最大處，螺旋環(huán)空中管柱屈曲位置在不停地變化，相應(yīng)的，其最大流速區(qū)與流動靜止區(qū)也在不停變化。這將極大地影響流體在環(huán)空中的流動，使得流體在環(huán)空中流動時產(chǎn)生更多的壓耗，故其最大流速區(qū)流速呈現(xiàn)逐漸減小的趨勢。

2.2 內(nèi)管旋轉(zhuǎn)時徑向流速分布

圖5—7為鉆柱旋轉(zhuǎn)（轉(zhuǎn)速為40 r/min）時不同環(huán)空結(jié)構(gòu)下徑向流速分布（因轉(zhuǎn)速為60，80 r/min時的流場分布規(guī)律與40 r/min時基本一致，這里不再給出）。將圖5—7分別與圖2—4進行對比可知，鉆柱旋轉(zhuǎn)對同心環(huán)空流場分布影響很小，幾乎可忽略不計，對偏心、屈曲環(huán)空結(jié)構(gòu)窄間隙處的流速卻有較大的影響，轉(zhuǎn)速越大，窄環(huán)空間隙處的流速越大。這是由于屈服冪律流體具有剪切稀釋性，鉆柱旋轉(zhuǎn)使流體黏度減小，相應(yīng)的壁面剪切力減小，故鉆柱旋轉(zhuǎn)有利于增大窄環(huán)空間隙處的流速。

圖8為鉆柱靜止時，不同環(huán)空結(jié)構(gòu)下的壓降分布。從圖中可以看出：同心、偏心環(huán)空結(jié)構(gòu)下的壓降隨井深的增加呈線性增大的趨勢，且偏心度越大，壓降越小，即偏心度與壓降成反比；模擬正弦屈曲環(huán)空壓降居于偏心度為0.5與完全偏心時的壓降之間，沿井深呈非線性增加，這是因為正弦屈曲環(huán)空時其等效偏心度一直在0.5～1.0不斷變化的原因；螺旋屈曲環(huán)空結(jié)構(gòu)下的壓降在0～5.0 m的井深范圍內(nèi)要低于完全偏心環(huán)空結(jié)構(gòu)下的壓降值，5.0～10.0 m后，其壓降值大于完全偏心環(huán)空下的壓降值，這是因為鉆柱屈曲產(chǎn)生額外壓耗的緣故，在所討論的10 m井深范圍內(nèi)，螺旋屈曲環(huán)空與同心環(huán)空壓降值最大達到6 kPa，且隨著井深的增加，其壓降值還會繼續(xù)增大。

圖5 鉆柱旋轉(zhuǎn)、管柱未發(fā)生屈曲時不同環(huán)空結(jié)構(gòu)速度云圖

圖6 鉆柱旋轉(zhuǎn)、管柱發(fā)生正弦屈曲時環(huán)空結(jié)構(gòu)速度云圖

圖7 鉆柱旋轉(zhuǎn)、管柱發(fā)生螺旋屈曲時環(huán)空結(jié)構(gòu)速度云圖

圖8 鉆柱靜止時不同環(huán)空結(jié)構(gòu)下的壓降分布

3 仿真結(jié)果驗證與應(yīng)用

由于現(xiàn)階段關(guān)于屈曲環(huán)空壓降的室內(nèi)實驗以及現(xiàn)場數(shù)據(jù)幾乎沒有，為了驗證CFD仿真試驗的準(zhǔn)確性，應(yīng)用湯明在俄克拉荷馬大學(xué)所獲得的關(guān)于同心、偏心環(huán)空結(jié)構(gòu)下壓降的實驗測量數(shù)據(jù)［11］與CFD仿真結(jié)果進行了對比，如圖9所示。仿真結(jié)果與實驗結(jié)果吻合度較高，同理，用CFD仿真屈曲環(huán)空結(jié)構(gòu)下的壓降也具有一定的準(zhǔn)確性。

圖9 仿真結(jié)果與實驗結(jié)果對比

鉆井過程中，鉆柱發(fā)生正弦屈曲時，偏心率會發(fā)生周期性變化，如果以同心環(huán)空來計算環(huán)空壓降，則壓降過大；如果以完全偏心環(huán)空計算，則壓降過小。為了避免因壓降計算過大或者過小而引起深部地層的井噴、井漏等安全問題，必須準(zhǔn)確計算鉆柱發(fā)生正弦屈曲時的壓降值。此時可將正弦屈曲環(huán)空與同心環(huán)空壓降比和偏心環(huán)空與同心環(huán)空壓降比相結(jié)合，利用Tang等［12］關(guān)于偏心環(huán)空與同心環(huán)空壓降比擬合式（3）與本文的仿真壓降數(shù)據(jù)相結(jié)合，將正弦屈曲環(huán)空近似等效為某一偏心度條件下的偏心環(huán)空（偏心度恒定）進行處理，從而找到一個適合正弦屈曲的偏心率，以達到真實條件下環(huán)空壓降的精確計算。

4 結(jié)論

1）同心環(huán)空中，流體流速沿徑向?qū)ΨQ分布，偏心、屈曲環(huán)空中偏離對稱分布，最大流速出現(xiàn)在環(huán)空間隙最大處，最小流速出現(xiàn)在環(huán)空間隙最小處；鉆柱旋轉(zhuǎn)對同心環(huán)空流速的影響不大，轉(zhuǎn)速越大，偏心、屈曲環(huán)空窄間隙處流速越大，有利于鉆井過程中巖屑的運移。

2）同心環(huán)空、偏心環(huán)空壓降值沿井深線性增加，偏心度越大，壓降值越?。磺h(huán)空壓降值沿井深非線性增加，隨鉆柱旋轉(zhuǎn)速度增加，環(huán)空壓降減小。

3）CFD仿真結(jié)果與實驗結(jié)果相對比，具有較高的準(zhǔn)確性，將正弦屈曲環(huán)空等效為壓降相等的偏心環(huán)空近似處理，可以更方便的計算環(huán)空壓降。

4）屈曲環(huán)空結(jié)構(gòu)下的壓降變化趨勢與常規(guī)的同心、偏心結(jié)構(gòu)有所不同，其壓降值沿井深方向呈非線性增加，發(fā)生周期性變化，考慮到鉆遇窄安全密度窗口地層時易引起的井噴、井漏等井下安全問題，精確計算屈曲環(huán)空結(jié)構(gòu)下的壓降值十分重要，本文為真實鉆井條件下環(huán)空壓降的計算提供了一定的指導(dǎo)。

［1］REDBERGER P J，CHARLES M E.Axial laminar flow in a circular pipe containing a fixed eccentric core［J］.The Canadian Journal of Chemical Engineering，1962，40（4）：148-151.

［2］HEYDA J F.A green′s function solution for the case of laminar incompressible flow between non-concentric circular cylinders［J］. Journal of the Franklin Institute，1959，267（1）：25-34.

［3］鐘鑫，孔慶祥.卡森流體軸向同心環(huán)空流動規(guī)律研究［J］.石油鉆采工藝，1987，9（5）：27-36

［4］劉永建，陳家瑯.偏心環(huán)空中冪律流體軸向?qū)恿鲏航档挠嬎悖跩］.石油鉆采工藝，1989，2（5）：9-15.

［5］VIEIRA NETOJ L，MARTINS A L，ATAIDE C H，et al.The effect of the inner cylinder rotation on the fluid dynamics of non-Newtonian fluids in concentric and eccentric annuli［J］.Brazilian Journal of Chemical Engineering，2014，31（4）：829-838.

［6］張景富，李邦達.冪律液體偏心環(huán)空螺旋流層流壓降計算［J］.石油鉆采工藝，1992，5（5）：13-29.

［7］高德利.油氣井管柱力學(xué)與工程［M］.東營：中國石油大學(xué)出版社，2006：35-38.

［8］MISKA S，CUNHA J C S.Helical buckling of long weightless strings subjected to axial and torsional loads［C］.New York：American Society of Mechanical Engineers，1995.

［9］MISKA S，QIU W，VOLK L，et al.An improved analysis of Axial Force along coiled tubing in inclined/horizontal wellbores［R］.SPE 37056，1996.

［10］張晉凱，李根生，劉廣軍，等.旋轉(zhuǎn)鉆柱內(nèi)赫巴流體流動模型與流場特性研究［J］.石油機械，2013，41（5）：15-18.

［11］湯明.波動壓力與裂縫性地層極限水平段長度研究［D］.成都：西南石油大學(xué)，2016.

［12］TANG M，AHAMED R，SRIVASTAV R，et al.Simplified surge pressure model for yield power law fluid in eccentric annuli［J］. Journal of Petroleum Science&Engineering，2016，145（38）：346-356.

（編輯 孫薇）

Simulation of annular flow considering buckling drilling tools

YUAN Lanfeng1,HE Shiming1,LI Xuenian2,TANG Ming1,SHAO Yong1
(1.State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation,Southwest Petroleum University,Chengdu 610050, China;2.Drilling Department,Huadong Petroleum Engineering Company,SINOPEC,Nanjing 210000,China)

During the process of slim-hole directional well and horizontal well drilling,the buckling of drilling tools can affect not only the mechanical properties of the drilling tool assembly,but also the annulus flow field and pressure field greatly.In this paper, the simulations of annular flow about the sinusoidal and helical buckling were carried out by ANSYS Fluent software,and the effects of different annular eccentricity and rotational speed on annular flow velocity field and flow pressure drop were studied systematically.The results show that the flow velocity is distributed symmetrically in the radial direction in the concentric annulus, however,the flow velocity distribution is deviated from the radial direction in the eccentric and buckling annulus;without considering rotation,the larger the annular clearance is,the higher the velocity is;the influence of drill string rotation on the flow velocity at the narrow gap is larger than that at the wide gap;the concentric annulus and eccentric annulus pressure drops increase linearly along the depth of the well,and the larger the eccentricity is,the smaller the pressure drop is;as the depth of well increases,the pressure drop of buckled annulus will also increase nonlinearly,the faster the speed of pipe rotation is,the smaller the pressure drop will be. The simulation results agree well with the experimental results,which has some guidance to the calculation of annular pressure drop under real drilling conditions.

eccentric annulus;sinusoidal buckling;helical buckling;pressure drop;ANSYS FLUENT software

TE347

A

國家自然科學(xué)基金項目“基于流-固-化耦合理論的欠平衡鉆水平井井壁失穩(wěn)機理研究”（51474186）、“基于多場耦合理論的頁巖氣水平井井壁失穩(wěn)機理研究”（51574202）；國家科技重大專項專題“欠平衡鉆井工藝技術(shù)研究”（2008ZX05045-03-07HZ）、“特殊結(jié)構(gòu)井井眼穩(wěn)定及摩阻控制技術(shù)”（2011ZX05045-03-01WX）

10.6056/dkyqt201705026

2017-02-20；改回日期：2017-06-22。

袁蘭峰，男，1992年生，主要從事油氣井流體力學(xué)與井控相關(guān)研究。E-mail：yuanlanfeng2015@sina.com。

何世明，男，1966年生，教授，博士生導(dǎo)師，主要從事油氣井工程力學(xué)（井壁穩(wěn)定、井控與欠平衡鉆井）和鉆完井新技術(shù)等方面的研究及教學(xué)工作。E-mail：hesming@sina.com。

袁蘭峰，何世明，黎學(xué)年,等.考慮鉆具屈曲的環(huán)空流動仿真研究［J］.斷塊油氣田，2017，24（5）：714-718.

YUAN Lanfeng，HE Shiming，LI Xuenian，et al.Simulation of annular flow considering buckling drilling tools［J］.Fault-Block Oil&Gas Field，2017，24（5）：714-718.