微裂縫內液固兩相流動特性研究

楊 倩，陳彥超

?

微裂縫內液固兩相流動特性研究

楊 倩，陳彥超

（東北石油大學， 黑龍江 大慶 163318）

應用歐拉-拉格朗日模型，對微裂縫中液-固兩相流動過程進行數(shù)值模擬。在數(shù)學模型中，液相采用連續(xù)相模型，顆粒當做離散分散在連續(xù)相中，并考慮了相間耦合作用。對于物理模型，建立了二維直微裂縫模型，研究裂縫的液相粘度、液相速度和顆粒密度等因素對裂縫中兩相流動的影響。模擬得到了微通道內顆粒瞬時分布狀態(tài)、顆粒運動軌跡、顆粒停留時間、顆粒軸向速度分布、液相軸向速度分布等參數(shù)的變化規(guī)律，定性的揭示了儲層裂縫內液固兩相流動過程。

微裂縫；計算流體力學；離散元方法；數(shù)值模擬

石油主要儲存巖石的儲集空間中，它是儲集油氣的場所。儲集空間包括粒間孔隙、粒內孔隙、裂縫、溶洞等各種類型的孔、洞、縫，就其形態(tài)和分布而論是十分復雜的孔喉網(wǎng)絡。儲集石油的沉積巖中含有豐富的粘粒礦物，砂巖顆粒以及石英晶體的細小顆粒[1,2]。石油運移時也會對這些微細顆粒有攜帶的作用，而固體微粒的存在可以影響原油運移的速度和路徑。在開采疏松砂巖油藏時，由于地層微粒隨石油的流動過程中會發(fā)生運移、沉積造成地層損害，使產油量下降，影響油田的正常生產[3]。因此，微粒遷移問題的研究對于疏松砂巖油藏的開發(fā)具有非常重要的意義。

如今，微尺度流動已經(jīng)得到了學者們的廣泛關注[4-8]，形成了一個新的學科分支。微觀尺度可能會出現(xiàn)與宏觀流動不同的規(guī)律，并且許多在宏觀流動中被忽略的因素，在微流動中可能會成為主要的影響因素。也就是說，當研究對象是微觀時，流動的尺度效應就會凸顯出來。近年來，大量的對微尺度流動特性的研究表明：常規(guī)尺度下的流動關聯(lián)式因為微通道的許多固有特點而失效，導致基于宏觀經(jīng)驗的傳統(tǒng)流體力學不再完全適用于微尺度現(xiàn)象，因此，需要對原有的一些規(guī)律和理論進行補充和修正[9]。

1 模型的建立

1.1 數(shù)學模型

利用Euler-Lagrange方法里的離散元（DEM）模型并結合計算流體動力學（CFD）的相關知識，對裂縫中的液-固兩相流動過程進行數(shù)值模擬。FLUENT 軟件可以在拉氏坐標下模擬流場中離散的第二相。由球形顆粒（代表油滴和砂粒）構成的離散分布在連續(xù)相中，計算這些顆粒的軌道以及由顆粒引起的熱量或質量傳遞。液體-顆粒流動遵循質量、動量和牛頓守恒定律［10,11］。

歐拉-拉格朗日模型的關鍵是計算顆粒的運動軌跡。在兩相流動的過程中，每個運動的顆粒都會受到流體以及相鄰顆粒對它的作用。因此，可將顆粒的運動分解為隨流體運動的懸浮過程和與其他運動顆粒的碰撞過程。

采用計算流體力學和離散單元方法對微裂縫內液固兩相流動過程進行了數(shù)值模擬。采用-ε模型模擬液體流動，離散單元方法模擬顆粒流動，Huilin- Gidaspow 曳力模型模擬相間作用力。考慮液體速度、粘度以及固體顆粒密度等因素對流動的影響。

1.2 物理模型

如圖1所示為裂縫網(wǎng)格劃分圖 ，依據(jù)此圖，本論文模擬了水平方向微裂縫內的液固兩相流動。具體參數(shù)如下：裂縫的長為1 000 μm，寬200 μm。左側為射流口，出口設為壓力出口。

圖1 裂縫網(wǎng)格圖

液相選擇液態(tài)水為參數(shù)，密度為1 000/m3。入口速度分別為0.1、0.3和0.5 mm/s。砂粒密度選取1 550、2 000和2 500 kg/m3，直徑為2 μm。油滴的密度為800 kg/m3，直徑為2 μm。流動為非定常流、松弛因子為1、考慮重力影響（9.81 m/s2）。

2 計算結果和分析

2.1 微裂縫內液固兩相流動行為分析

圖2為Fluent模擬的流體速度為0.1 mm/s，粘度為0.001 kg/m.s算例不同時刻油滴及砂粒顆粒分布圖，左側為入口處，兩種顆粒的初始射入速度為0m/s，液體速度為0.1 mm/s。圖中可以清楚看出油和砂粒在微裂縫中的分布情況。兩種顆粒隨液體流動而運移，最后由壓力出口流出。兩種顆粒均勻混合，未表現(xiàn)出明顯的重力作用。

圖2 不同時刻顆粒分布圖

Fig.2 Particle distribution of different time

圖3為流體、油滴及砂粒的速度矢量圖，由圖可知流體和顆粒運動方向及速度大小。液體的流動是比較有序的，在壁面附近流體出現(xiàn)回流，沿壁面向左運動，在模型中有多處呈旋渦流動狀態(tài)；對于油滴顆粒來說，顆粒速度較小，在邊界處的顆粒碰擊壁面速度較大；而固相砂粒的速度相對較大，而且表現(xiàn)出強烈的無序性。

圖3 三相速度矢量圖

圖4分別表示油滴和砂粒運動軌跡圖，由圖可以看出，顆粒的初始堆積位置不一致，使得不同顆粒的運動軌跡明顯不同，但總體趨勢是向出口運動。距離射流口較近的位置，顆粒的運動范圍較大，這主要是因為入口射流為顆粒提供更大的能量，顆粒更活躍，致使運動范圍擴大。

圖4 油滴-砂粒運動軌跡

2.2 液相粘度對顆粒流動特性的影響

圖5表示在液體速度均為0.1 mm/s條件下，平行于入口距離300 μm直線上液體軸向速度和切向速度沿軸向的分布。隨著粘度的增加，液體在中心線處運動劇烈，靠近壁面處速度逐漸降低。顆粒的徑向速度受粘度影響較大，并且與粘度大小不成正比，在模型的上部出現(xiàn)負速度，說明液體受重力作用向下運動，形成裂縫內的循環(huán)流動。而在接近壁面處變化較小，趨于穩(wěn)定。

同正：撲鼻的香料、櫻桃香氣，十分有爆發(fā)力。入口酸度有力，帶來十分不錯的清爽感。酒體中等，單寧柔滑。余味中可可、巧克力香氣明顯。

圖5 不同粘度下液相軸向和切向速度

圖6分別是，油滴和砂粒的軸向速度圖。由圖中可以看出，兩種顆粒速度變化受流體粘度影響較小，并且沿裂縫中心處成軸對稱分布。油滴速度略大于砂粒速度，兩種顆粒初始速度約為零，隨著流體的沖擊作用，速度增大到流體速度大小的1/4，由此可見，流體的流動可以有效的攜帶顆粒運動。

圖6 不同粘度下油滴-砂粒軸向速度

2.3 液相速度對顆粒流動特性的影響

在同一粘度的條件下，改變液體的速度觀察顆粒運動特性，圖7分別是油滴和砂粒的軸向速度圖。由圖中可以看出，隨著流體速度增大，兩種顆粒速度都明顯增大，而油滴的運動更為劇烈，最大速度約是砂粒的2.5倍。砂粒運動相對平緩，這是由于砂粒密度比油滴大，所以受流體作用力相對較小。

圖7 不同速度下油滴-砂粒軸向速度

2.4 固相砂粒密度對流動的影響

在模擬流動中密度大的顆粒在流動中更容易發(fā)生沉降，顆粒的停留時間增大。在砂粒密度為2 500 kg/m3時，微裂縫中顆粒整體呈向下的旋轉流動，說明固體顆粒的密度對整個流場有很大的作用力影響。

圖8 不同密度下砂粒軸向速度圖

圖9 不同密度下液體軸向速度圖

圖8-圖9分別是以砂粒密度為變量的條件下固液兩相瞬時軸向速度圖，由圖可知，密度越大，砂粒脈動越劇烈，而密度對液體流動狀態(tài)的影響較小，隨著固相砂粒密度的增大，液體速度的波動范圍有減小的趨勢。

3 結 論

對儲層裂微縫進行系統(tǒng)結構的分析，應用Gambit建模軟件完成對裂縫二維模型的建立。在以液固兩相流的基本理論為基礎，采用歐拉-拉格朗日雙流體模型，建立數(shù)學模型，應用Fluent軟件對微裂縫內固液兩相流動狀態(tài)進行數(shù)值模擬。結果表明Fluent 可以在拉氏坐標下模擬流場中離散的第二相，由球形顆粒（代表油滴和砂粒）構成的第二相分布在連續(xù)相中，并且可以計算這些顆粒的運動軌跡以及由顆粒引起的能量傳遞。通過單一變量方法模擬案比較不同條件下的流動狀態(tài)，發(fā)現(xiàn)當液體入口速度相對較大時、顆粒運動速度也越大，但并不成正比；顆粒的停留時間減小。而液體粘度對流動的影響較大，當液體粘度增大時，顆粒運移明顯，表現(xiàn)為顆粒停留計算域時間更短。固相密度能明顯的影響顆粒的運動軌跡，密度大的顆粒在流動中更容易發(fā)生沉降。隨著運動顆粒密度的增大，顆粒的瞬時軸向速度減小，顆粒的停留時間增大。顆粒密度的大小對液相流動狀態(tài)的影響表現(xiàn)為：隨著固相顆粒密度的增大，液相速度的波動范圍有減小的趨勢。

［1］紀友亮.油氣儲層地質學[M].第2版.東營：中國石油大學出版社,2009:6-7.

［2］ Tannenbaum, E., Huiainga, B., Kaplan, I.. Role of minerals in thermal alteration of organic matter. II—a material balance[J]. American Association of Petroleum Geologists Bulletin, 1986,70:1156-1165.

［3］王尤富.疏松砂巖地層微細粒運移特征實驗研究[J].特種油氣藏,2005,12(5):91-92.

［4］Gao P, Peterson S L. Favre-Marinet M. Seale Effects on Hydrodynamics and Heat Transfer in Two-dimensiolal Mini and Microchannels [J]. Int. J. Ther. Sei, 2002, 41 (11): 1017-1027.

［5］Pfund D A, Shekarriz A, Popescu A, et al. Pressure Drops Measurements in Microchannels[J]. Proceedings of MEMS, ASME DSC. 1998 (66): 193-198.

［6］Herwig H, Hausner O. Critical View on“New Results in Micro-fluid Mechanics”an Example[J]. Int. J. Heat and Mass Transfer, 2003(46): 935-937.?

［7］Guo Z Y, Li Z X. Size Effect on Microscale Single-phase Flow and Heat Transfer[J]. Int. J. Heat and Mass Transfer, 2003 (46): 149-159.

［8］劉君.微通道流動特性的數(shù)值分析[D].哈爾濱工業(yè)大學,2010.?

［9］聶海寬,張金川,張培先,等.福特沃斯盆地Bamet頁巖氣藏特征及啟動［J］.地球科學進展，2009,28（2）：87-93.

［10］Gidaspow, D., Multiphase flow and fluidization: continuum and kinetic theory descriptions[M]. Academic Press, San Diego, 1994.

［11］王嘉駿，顧雪萍，楊富軍，馮連芳，雙流體模型中曳力及恢復系數(shù)對氣固流動的影響[J].高校化學工程學報，20，164-168，2006。

Study on Flow Characteristics of Liquid-solid Two-phase Flow in Microfractures

YANG Qian，CHEN Yan-chao

（Northeast Petroleum University, Heilongjiang Daqing 163318，China）

Numerical simulation of liquid-solid two-phase flow in micro-fractures was carried out by Euler-Lagrange model. In the mathematical model, the liquid phase used the continuous phase model, the particles phase was treated as discrete phase dispersed in a continuous phase, and the coupling effect was considered. For the physical model, a two-dimensional direct micro-fracture model was established to study the effect of liquid viscosity, liquid velocity and particle density on the two-phase flow in fracture. The instantaneous distribution of particles in micro-fracture was obtained as well as the particle trajectory and particle retention time, particle axial velocity distribution, the variation law of the axial velocity distribution in liquid phase. The process of liquid-solid two-phase flow in micro-fracture reservoir was revealed qualitatively.

Micro-fracture; Computational fluid dynamics; Discrete element method; Numerical simulation

TE 122

A

1671-0460（2016）06-1194-04

2016-04-15

楊倩，女，黑龍江大慶人，在讀碩士研究生，就讀于東北石油大學石油與天然氣專業(yè)，研究方向：復雜流體力學。E-mail：1036573319@qq.com。