房建宇，李士偉

?

基于顆粒直徑大小對(duì)環(huán)形空間攜屑影響的數(shù)值模擬

房建宇1，李士偉2

（1. 東北石油大學(xué) 石油工程學(xué)院，黑龍江 大慶 163318 ; 2. 大慶油田采油六廠，黑龍江 大慶 163000）

鉆井是開發(fā)石油資源一個(gè)非常重要的環(huán)節(jié)和技術(shù)手段，井中巖屑易在井眼底部形成巖屑床，導(dǎo)致下鉆遇阻、蹩泵甚至卡鉆，因此巖屑運(yùn)移機(jī)理的研究對(duì)鉆井作業(yè)來(lái)說(shuō)至關(guān)重要。通過(guò)對(duì)現(xiàn)場(chǎng)實(shí)際情況的調(diào)研，得到井下上返顆粒的直徑范圍，進(jìn)行數(shù)值模擬。運(yùn)用SolidWorks對(duì)井下環(huán)形空間進(jìn)行建模與裝配，得到井下環(huán)形空間的三維模型。使用軟件建立井下環(huán)形空間的簡(jiǎn)化模型，通過(guò)應(yīng)用FLUENT中動(dòng)網(wǎng)格部分的UDF編寫成功地實(shí)現(xiàn)了鉆桿的旋轉(zhuǎn)。運(yùn)用FLUENT數(shù)字模擬軟件，通過(guò)對(duì)流體非定常流的數(shù)值模擬，對(duì)巖屑在環(huán)形空間的運(yùn)移規(guī)律進(jìn)行研究，模擬條件更加接近實(shí)際情況，結(jié)果更加精確，可以為鉆井井眼清潔工作提供科學(xué)參考。

計(jì)算流體力學(xué)(CFD); 環(huán)形空間; 液固兩相流; 數(shù)值模擬

井眼中的固體微粒受重力、粘滯阻力、沖擊力、浮力作用影響，巖屑濃度大于這些作用力的工作能力，巖屑就會(huì)沉降產(chǎn)生巖屑床[1]；巖屑床嚴(yán)重影響了機(jī)械鉆速，鉆進(jìn)是阻力增加，尤其是定向井、水平井時(shí)易形成井下的脫壓現(xiàn)象，導(dǎo)致鉆頭處沒有鉆壓，長(zhǎng)時(shí)間沒有進(jìn)尺，工程進(jìn)度緩慢，增加鉆井成本；巖屑床導(dǎo)致鉆柱扭矩增大，嚴(yán)重時(shí)甚至扭斷井內(nèi)鉆柱，產(chǎn)生嚴(yán)重的井下事故；巖屑床還是鉆具粘卡產(chǎn)生的主要原因，巖屑床的存在使測(cè)試工具下方受阻。

本文對(duì)井下環(huán)形空間模型進(jìn)行了優(yōu)化設(shè)計(jì)，運(yùn)用液固兩相流理論知識(shí)，并合理模擬在鉆桿旋轉(zhuǎn)的情況下不通顆粒直徑大小對(duì)井下巖屑運(yùn)移的影響。

1 計(jì)算模型和邊界條件

井下環(huán)形空間三維模擬圖，包括井筒，鉆桿。井筒長(zhǎng)度12 m，半徑10.5 cm，鉆桿長(zhǎng)度11.5 m，半徑5.5 cm在鉆桿底部開設(shè)小孔，小孔直徑0.5 cm。

首先根據(jù)井下實(shí)際工況運(yùn)用SolidWorks對(duì)環(huán)形空間進(jìn)行建模與裝配，得到環(huán)形空間的三維模型(環(huán)形空間主要由井筒和鉆桿兩部分構(gòu)成）。然后導(dǎo)入到GAMBIT軟件建立井下工況環(huán)形空間的動(dòng)態(tài)網(wǎng)格如圖1所示，繼而采用FLUENT軟件中的UDF模板對(duì)井筒內(nèi)鉆桿的流動(dòng)進(jìn)行程序的編寫，環(huán)形空間內(nèi)的流體為非定常流動(dòng)，對(duì)此情況下的流動(dòng)進(jìn)行三維仿真計(jì)算。入口邊界條件定義為速度入口，壓力為表壓；出口邊界條件定義為壓力出流邊界；壁面無(wú)滑移，默認(rèn)即可；中心線定義為軸對(duì)稱條件；材料特性：液體設(shè)置為水，固體顆粒設(shè)置為二氧化硅。壓力速度藕合方式選擇SIMPLE，壓力與動(dòng)量的離散格式和湍動(dòng)能及湍動(dòng)能耗散率默認(rèn)即可；time為非定常流動(dòng)；松弛因子為1；采用歐拉方程隱式算法。

2 不同顆粒直徑對(duì)鉆井?dāng)y巖效果影響

由圖2、圖3可以看出：在鉆桿旋轉(zhuǎn)時(shí)，鉆井液在鉆桿旋轉(zhuǎn)的離心力作用下，使得鉆井液圍繞鉆桿做周向運(yùn)動(dòng)（也稱作剪切運(yùn)動(dòng)），同時(shí)鉆井液在井內(nèi)壓差的作用下，沿著環(huán)形空間中軸線的方向做層流運(yùn)動(dòng)，鉆井液的周向運(yùn)動(dòng)和軸向運(yùn)動(dòng)共同構(gòu)成了鉆桿旋轉(zhuǎn)時(shí)鉆井液在環(huán)空中的螺旋流流動(dòng)。

圖2 不同顆粒直徑的顆粒環(huán)空體積濃度圖

鉆井液作螺旋流流動(dòng)時(shí),鉆井液的軸向流動(dòng)主要對(duì)巖屑在井內(nèi)產(chǎn)生向上的升力，而鉆井液的周向運(yùn)動(dòng)對(duì)巖屑起到保持平衡的懸浮力的作用，在上升力和懸浮力等的作用下使得具有沉積趨勢(shì)的顆粒，或側(cè)壁已經(jīng)沉積的顆粒，或者井底沉積的顆粒重新進(jìn)入到環(huán)空高速流動(dòng)區(qū)域內(nèi)。由以上可以看出，在鉆桿旋轉(zhuǎn)的情況下[2]，巖屑能夠獲取更大的動(dòng)能從而運(yùn)移出井口。

晚清以降，三千年未有之大變局警醒了部分國(guó)人，亡國(guó)滅種的危機(jī)逼迫這些“先知先覺”者告別相對(duì)穩(wěn)定的時(shí)間意識(shí)和生活狀態(tài)。風(fēng)行一時(shí)的《天演論》正在重塑他們的世界觀，逐漸形成一種奠基于進(jìn)化論的時(shí)間觀。龔自珍《己亥雜詩(shī)》之四四曰：“何敢自矜醫(yī)國(guó)手，藥方只販古時(shí)丹?！钡诤m看來(lái)，中國(guó)傳統(tǒng)的“古時(shí)丹”已根本無(wú)法應(yīng)對(duì)“今”之?。骸敖袢瘴釃?guó)之急需……以吾所見言之，有三術(shù)焉，皆起死之神丹也：一曰歸納的理論，二曰歷史的眼光，三曰進(jìn)化的觀念。”[注]胡適：《留學(xué)日記》卷3，胡適著、季羨林主編：《胡適全集》第27卷，合肥：安徽教育出版社，2003年，第261頁(yè)。

圖3 顆粒環(huán)空橫截面矢量圖

Fig.3 Vector profile in particle annular cross section

同心環(huán)空中液固兩相流動(dòng)紊流的軸向速度分布情況如圖4所示：由圖可以看出環(huán)形空間軸向速度呈軸對(duì)稱分布。軸向速度在環(huán)形空間中心存在一個(gè)基本穩(wěn)定、類似于流核的區(qū)域,速度趨于均勻，這是由于在紊流的條件下，由于顆粒與顆粒之間的，以及液體和顆粒之間的的互相摻混和碰撞，使得固體顆粒的運(yùn)動(dòng)方向雜亂無(wú)章，因此，導(dǎo)致軸向橫截面的速度分布趨于均勻。然后向管壁快速降低,形成速梯區(qū)，這是由于在近壁面顆粒受到壁面的阻力的影響所以速度逐漸降低。由圖還可以看出隨著顆粒直徑的增加，顆粒的軸向速度逐漸增大，到8 mm速度不再增加，說(shuō)明存在一個(gè)臨界顆粒直徑范圍使得顆粒在環(huán)空中軸向速度達(dá)到最大，此項(xiàng)目仍在研究中[3,4]。

圖4 不同顆粒直徑的顆粒軸向速度圖

由圖5可以看出隨著顆粒直徑的增加，顆粒環(huán)空體積分?jǐn)?shù)先增加后減少，當(dāng)顆粒直徑較小時(shí)，顆粒由于慣性作用顆粒直徑越大軸向速度越大，液體能攜帶更多的顆粒進(jìn)入高速流核區(qū)，使得顆粒環(huán)空體積分?jǐn)?shù)增大。

圖5 不同顆粒直徑的顆粒環(huán)空體積分?jǐn)?shù)圖

當(dāng)顆粒直徑增大到一定極限由于自身重力作用，顆粒下滑速度逐漸增大，環(huán)空軸向速度不再增加，顆粒重新返回至井底或形成巖屑床，環(huán)空顆粒體積分?jǐn)?shù)降低[5-7]。

3 結(jié)論及建議

本文通過(guò)對(duì)井下環(huán)形空間的液固流動(dòng)進(jìn)行數(shù)值模擬，并且運(yùn)用液固兩相流理論基礎(chǔ)，使用SolidWorks軟件建立井下三維環(huán)形空間模型。然后利用UDF模擬環(huán)形空間內(nèi)鉆桿的轉(zhuǎn)動(dòng)，并用FLUENT對(duì)流場(chǎng)進(jìn)行數(shù)值模擬，進(jìn)一步分析了不同顆粒直徑對(duì)環(huán)空攜巖效果的影響。得出以下結(jié)論：鉆桿旋轉(zhuǎn)有利于巖屑的運(yùn)移。隨著顆粒直徑的增加，顆粒的軸向速度先增大后趨于平緩，說(shuō)明存在臨界直徑，使得環(huán)空返速最大，環(huán)空顆粒體積分?jǐn)?shù)達(dá)到峰值，更有利于巖屑的運(yùn)移，這對(duì)于現(xiàn)場(chǎng)施工應(yīng)用具有一定的參考意義。

［1］ 宋巍. 水平井巖屑運(yùn)移規(guī)律研究[D]. 成都：西南石油大學(xué)，2013.

［2］ 陳浩，李冰，史禹. 大井斜大位移井鉆井液攜巖問題分析及對(duì)策[J]. 勝利油田職工大學(xué)學(xué)報(bào)，2009, 23(3)：54-56.

［3］ 袁照永. 水平井環(huán)空泡沫攜巖流動(dòng)規(guī)律研究[D]. 青島：中國(guó)石油大學(xué)，2009.

［4］ 鹿傳世．傾斜井眼內(nèi)巖屑運(yùn)移的數(shù)值計(jì)算與清巖工具設(shè)計(jì)流場(chǎng)計(jì)算[D]．東營(yíng)：中國(guó)石油大學(xué)（華東），2008．

［5］ 王延民，孟英峰，李皋，顏建華，李永杰，林鐵軍，莒拉提. 氣體鉆井出水條件下環(huán)空攜巖機(jī)理研究[J]. 重慶科技學(xué)院學(xué)報(bào)(自然科學(xué)版)，2008(4)：1-3.

［6］ G. Espinosa-Paredes, O. Cazarez-Candia. Two-Region Average Model for Cuttings Transport in Horizontal Wellbores I: Transport Equationsroleum Science and Technology[G] . 2011 (13) .

［7］ Chang-Heon Song; Ki-Beom Kwon; Min-Gi Cho; Joo-Young Oh; Dae-Young Shin;Jung-Woo Cho.Development of lab-scale rock drill apparatus for testing performance of a drill bit[C]. International Journal of Precision Engineering and Manufacturing，2015-07.

Numerical Simulation on Influence of Particle Diameter to Cuttings Transportation in Annular Space

FANG Jian-yu1，LI Shi-wei2

(1. School of Petroleum Engineering, Northeast Petroleum University, Heilongjiang Daqing 163318, China 2. Oil Production Plant No.6, Daqing Oilfield Co.,Ltd., Heilongjiang Daqing 163000, China)

Drilling is a very important link and technical means in the development of oil resources. Cuttings in the well are easy to form the cuttings bed at the bottom of the borehole, which will cause slacking off, pump chocking up and even pipe-sticking. So study of cuttings migration mechanism is crucial for drilling operations. Based on research of the actual situation, the diameter range of upward-moving particle in the borehole was obtained, and then numerical simulation was carried out.SolidWorks was used to model and assemble the downhole annular space, the 3D model of the downhole annular space was obtained. The software was used to establish a simplified model of the downhole annular space. Through application of the UDF writing in the moving grid of FLUENT, the rotation of the drill rod was achieved successfully. Using the numerical simulation software FLUENT, through the numerical simulation of fluid unsteady flow, the migration rule of cuttings in the annular space was studied. The simulation conditions are more close to the actual situation. The result is more accurate, and it can provide scientific reference for drilling hole cleaning.

computational fluid dynamics (CFD); annular space; liquid-solid two-phase flow; numerical simulation

TQ 018

A

1671-0460（2016）06-1276-03

2016-03-08

房建宇（1992-），男，黑龍江省大慶市人，碩士學(xué)位，2010年畢業(yè)于東北石油大學(xué)油氣儲(chǔ)運(yùn)專業(yè)，研究方向：復(fù)雜流體力學(xué)。E-mail：1907017622@qq.com。