空氣鉆井偏心環(huán)空氣固兩相流動數(shù)值模擬

2011-09-20 08:40:44赫娟張清兵劉康任強燕青海油田采油三廠

油氣田地面工程 2011年5期

赫娟張清兵劉康任強燕青海油田采油三廠

赫娟張清兵劉康任強燕青海油田采油三廠

基于歐拉多相流模型，建立空氣鉆井偏心環(huán)空氣固兩相流動數(shù)學模型，借助FLUENT軟件對計算模型進行求解；得到了穩(wěn)定流動后井筒內氣固兩相流動速度場、壓力場及兩相濃度分布規(guī)律。分析表明，整個流場沿著兩圓中心連線呈對稱分布，寬、窄間隙處軸向速度差異較大；且受邊界層效應的影響，靠近內、外兩圓壁面處流速最小。出口處巖屑速度小于空氣速度；而軸向越靠近出口，寬間隙區(qū)域空氣速度越大，窄間隙區(qū)域空氣速度越小。

空氣鉆井；環(huán)空流動；氣固兩相流；數(shù)值模擬；速度

空氣鉆井是以空氣作循環(huán)介質攜帶巖屑的鉆井技術，作為欠平衡鉆井技術的一個分支，因其具有鉆速快、成本低，對儲層損害小等特點備受人們的關注[1-2]。但氣體鉆井的整體技術還有待于進一步完善，尤其環(huán)空氣固兩相流場分布的控制還需要進一步的研究[3]，以求有效衡量空氣的攜巖能力和沖蝕效應。因此，本文從計算流體動力學出發(fā)，基于歐拉多相流模型，建立空氣鉆井偏心環(huán)空氣固兩相流動數(shù)學模型，借助FLUENT軟件對計算模型進行求解。

1 模型的建立

1.1 物理模型

取12m環(huán)空井段進行模擬分析，環(huán)空外徑為316.5mm，環(huán)空內徑為127mm，套管與鉆桿的絕對粗糙度均為0.04678mm。進口為速度入口邊界條件，出口為自由出口邊界條件[4]。采用結構化六面體網(wǎng)格對模型進行單元劃分。圖1為偏心度e= 0.5的環(huán)空偏心流場計算網(wǎng)格示意圖。

圖1 計算區(qū)域網(wǎng)格模型

1.2 數(shù)學模型

任何流體的流動均滿足連續(xù)性方程和動量守恒方程，假設流體的定常、等溫、湍流模型采用標準的K-ε雙方程模型。故描述三維流動的控制方程為：

（1）湍流脈動動能方程（K方程）。

湍流黏度

式中GK為平均速度梯度引起的湍動能k的產(chǎn)生項；Gb為浮力引起的湍動能k的產(chǎn)生項；YM為可壓縮湍流中脈動擴張貢獻；模型中C1ε=1.44，C3ε=1，Cμ=0.09，為經(jīng)驗常數(shù)。

（3）連續(xù)性方程。

式中αq為第q相流體在單元中的體積分數(shù)，mpq為第p項到第q項的輸運質量。

（4）動量方程。

（5）能量方程。

2 數(shù)值模擬

偏心度e為內外管圓心間距與內外管半徑差之比，即e=d/(R-r)。模擬的巖屑密度為2600kg/ m3，比熱為1100J/kg·K，導熱系數(shù)為0.3W/m·K。氣固兩相混合后，從下部入口進入，空氣流速25m/s，巖屑流速15m/s，巖屑體積分數(shù)0.02，顆粒直徑0.001m。圖2、圖3分別給出了出口及軸向截面巖屑和空氣的速度場分布云圖。分析可知：整個流場沿著兩圓中心連線呈對稱分布，寬、窄間隙處軸向速度差異較大；且受邊界層效應的影響，靠近內、外兩圓壁面處流速最小；出口處巖屑速度小于空氣速度；而軸向越靠近出口，寬間隙區(qū)域空氣速度越大，窄間隙區(qū)域空氣速度越小。