凌燕, 張宏兵, 潘益鑫, 王強(qiáng), 尚作萍

(1.河海大學(xué)地球科學(xué)與工程學(xué)院, 江蘇 南京 210098; 2.河海大學(xué)力學(xué)與材料學(xué)院, 江蘇 南京 210098)

0 引 言

注水開發(fā)油田的油井采出液通常是原油伴隨其他流體混合物[1]。近年來(lái),中國(guó)東部油區(qū)原油中水和沙的含量較高,在脫水處理前的管道中形成了特殊的水包油乳狀液,構(gòu)成了油-水乳狀液和粒徑范圍較寬的沙在混輸管道內(nèi)的固液多相流動(dòng)[2-5]。由于沙顆粒的存在,使管道液相流動(dòng)特性參數(shù)分布特征發(fā)生改變,對(duì)原油輸運(yùn)、產(chǎn)液剖面測(cè)井方法設(shè)計(jì)產(chǎn)生影響。眾多學(xué)者對(duì)兩相流的研究始于20世紀(jì)50、60年代[6-7],20世紀(jì)80、90年代在兩相流的研究領(lǐng)域己經(jīng)取得巨大的進(jìn)步[6-9],尤其是對(duì)氣-液兩相流、氣-固兩相流、液-固兩相流等的研究。但是,由于油-水-固多相流流變特性的復(fù)雜性和特殊性,國(guó)內(nèi)外學(xué)者對(duì)其油水兩相分布特征和固態(tài)的沉降研究較少。因此,本文針對(duì)實(shí)際生產(chǎn)的應(yīng)用,建立水平彎曲管道基于混合模型[10-12](Mixture模型)進(jìn)行液相和固相的分布特征和流動(dòng)特征數(shù)值模擬分析[14-17]。

1 固液多相流瞬態(tài)數(shù)值模擬

1.1 混合模型

假設(shè)混合物為連續(xù)相和分散相,其中連續(xù)相為液體或氣體,所涉及到的物理量以c表示;分散相為顆粒、汽泡或液滴,所涉及到的物理量以p表示?；旌夏Ｐ椭械倪B續(xù)方程和動(dòng)量方程是基于連續(xù)相和分散相的混合體得出的。分散相的濃度可由該相的連續(xù)方程解得[12-16]。

(1) 獨(dú)立相方程。多相流中每一相k的連續(xù)性方程和動(dòng)量方程為

(1)

-αkpk+·[αk(τk+τTk)]+αkρkg+Mk

(2)

式中,uk為相k的平均速度;ρk為相k的平均密度;αk為相k的體積分?jǐn)?shù);Tk為相k在界面處的質(zhì)量生成率;Mk相k的動(dòng)量源項(xiàng);τk為黏性應(yīng)力項(xiàng);τTk為湍流應(yīng)力項(xiàng)

(3)

式中,ρIk為局部密度;uFk為脈動(dòng)速度,即uFk=uIk-uk。

(2) 混合連續(xù)方程

(4)

由此可得混合物的連續(xù)方程

(5)

式中,ρm為混合密度;um為混合速度,表示質(zhì)量中心速度。

(3) 混合動(dòng)量方程。由動(dòng)量方程(2)疊加所有相可得

(6)

(4) 單相體積分?jǐn)?shù)方程

(7)

若不發(fā)生相變化,則連續(xù)方程可寫為

(8)

式(8)也被作為擴(kuò)散方程(the Diffusion Equation),因此,混合模型也被稱為擴(kuò)散模型。

1.2 邊界條件

為簡(jiǎn)化起見(jiàn),針對(duì)液體(水或油)和固體顆粒兩相情況,假設(shè)水平井內(nèi)為充分發(fā)展的穩(wěn)態(tài)湍流流動(dòng),液體是不可壓縮的(密度為常數(shù)),流動(dòng)是等溫或溫度變化較小(動(dòng)力黏度為常數(shù)),則

(9)

(10)

這時(shí),水平井中混合相滿足封閉控制方程,其邊界條件為

(1) 對(duì)混合相(挾沙水或油),流速um、vm和wm在管壁上為0。

(2) 在流場(chǎng)對(duì)稱軸Y軸上,各流動(dòng)特性參數(shù)有對(duì)稱性條件式(11)。

(11)

1.3 有限差分格式

上述偏微分方程及邊界條件構(gòu)成1組封閉的、適定的偏微分方程組,可以通過(guò)數(shù)值解方式進(jìn)行求解。采用有限差分方法[2,13]對(duì)水平井圓截面進(jìn)行網(wǎng)格剖分(見(jiàn)圖1),對(duì)偏微分方程及邊界條件進(jìn)行離散化處理,并采用逐次超松弛(Successive Over Relaxation,SOR)迭代法進(jìn)行數(shù)值求解。

圖1 水平圓截面有限差分剖面圖

首先將圓截面沿x、y軸上半徑M等分,劃分圓截面為若干個(gè)網(wǎng)格,取最靠近邊界的網(wǎng)格為管壁,令

Δx=Δy=r/M

xi=i·Δx,yj=r-j·Δy,

i=0,1,2,…,M,j=0,1,2,…,2M

(12)

式中,Δx和Δy分別為x和y方向網(wǎng)格長(zhǎng)度;r為半徑;xi和yj分別為網(wǎng)格的x和y坐標(biāo)值。對(duì)偏微分方程中一階、二階偏導(dǎo)數(shù)離散化,分別得

(13)

(14)

和

(15)

(16)

(17)

使用上述差分格式對(duì)控制方程和邊界條件進(jìn)行離散,對(duì)于離散化后所得到的2M×M線性方程組直接求解非常困難,這里采用逐次超松弛迭代法(SOR),其中松弛因子ω要求滿足1<ω<2。首先給定各流動(dòng)特性參數(shù)Φ的初值Φ0,然后使用SOR法依次求解,獲得各流動(dòng)特性參數(shù)Φn,重復(fù)上述過(guò)程,直到獲得混合相速度場(chǎng)和顆粒相體積分?jǐn)?shù)場(chǎng)的最終收斂解。在使用SOR法求解每個(gè)控制方程時(shí),也需要反復(fù)迭代計(jì)算,直到滿足終止迭代條件式(18),即

(18)

這里ε為迭代控制精度,可以為任意小正數(shù)。該參數(shù)取值偏大則達(dá)不到精度要求,取值偏小則收斂太慢,一般取10～10-7為宜。

圖2 模型1彎曲管道內(nèi)固液多相流中油相和水相體積分布

2 固液多相流各相體積分布特征

為進(jìn)一步了解水平井彎曲管道中各相分布特征,設(shè)計(jì)2個(gè)模型,模型直徑都為0.124 m,管道最大傾角5°。其中模型1為向下凹,進(jìn)口處油水分界線位于中線以上0.02 m;模型2為向下凹,進(jìn)口處分界線位于中線以下0.02 m,數(shù)值模擬中入口端管道軸向速度分別為0.01、0.05、0.1、0.2 m/s和0.3 m/s,相當(dāng)于流量分別為10.43、52.17、104.34、208.68 m3/d和313.02 m3/d,入口端含沙體積分?jǐn)?shù)為1.0%、2.0%和5.0%(見(jiàn)表1、表2)。

表1 特征參數(shù)的初值

表2 數(shù)值模擬參數(shù)

2.1 油相和水相體積分布特征

圖2(a)和2(b)為模型1的油相和水相的體積分布,其中入口流速為0.05 m/s,含沙2%,y=0.02 m,可以看出,在管道下降段油相體積分布區(qū)域大于上升段的油相體積,水相則相反。上述現(xiàn)象可以解釋為,在下降段油相速度明顯小于入口流速,而下降段水相速度明顯大于入口流速,為了保持總流量及各相持率的不變,致使下降段油相體積增大,而水相體積減少。同樣,在上升段油相速度大于入口流速,而下降段水相速度小于入口流速,致使下降段油相體積減少,而水相體積增大。此外,隨著流體速度的增加,如入口流速為0.2 m/s,含沙2%,y=0.02 m,油相體積分布區(qū)域變小,水相體積則相反(見(jiàn)圖3)。圖4為5%含沙的油相和水相體積分布的模擬結(jié)果,其他參數(shù)與圖2中使用的相同,可以發(fā)現(xiàn),含沙為5%的結(jié)果與2%的結(jié)果相差很小,這是由于2種含沙量的軸向速度分布相差不大所致。

圖3 模型1彎曲管道內(nèi)固液多相流中油和水相體積分布

圖4 模型1彎曲管道內(nèi)固液多相流中油相和水相體積分布

模擬模型2結(jié)果如圖5(a)和5(b)所示。在圖5模型2的油相和水相的體積分布結(jié)果中,入口流速為0.05 m/s,含沙2%,y=0.02 m,可以看出,在管道下降段油相體積分布區(qū)域大于上升段的油相體積,水相則相反,這與模型1的情況相似。對(duì)比模型2與模型1,可以發(fā)現(xiàn),模型2的油相體積要大于模型1,即隨著油水界面下降,油相體積增大,并且下部水相的軸向流速也略有降低(見(jiàn)圖5)。

圖5 模型2彎曲管道內(nèi)固液多相流中油相和水相體積分布

2.2 固相體積分布特征

在上述模擬基礎(chǔ)上進(jìn)一步分析彎曲管道沉沙情況。圖6提供了模型1中沙顆粒體積分布,其中入口流入速度分別為0.01、0.05、0.1 m/s和0.2 m/s,含沙2%。從圖6可以看出,沙顆?；旧隙汲恋碓趶澢艿赖牡撞?并且隨著流速增大,在下凹彎曲段沉沙量先減少、后增多并向前推進(jìn)。這一現(xiàn)象可以解釋為①沙顆粒逐漸沉淀在底部,由二維模擬結(jié)果可知,隨著流速增大沉沙量逐漸減少;②隨著流速進(jìn)一步增大,達(dá)到起揚(yáng)速度,部分沉淀的沙顆粒被帶起向前推進(jìn),而留下的空間進(jìn)一步被沉下顆粒補(bǔ)充,并形成圖6(c)、(d)中傾斜的沙平面。為此,對(duì)于大流速或流量,沉積下來(lái)的沙顆粒在水的起揚(yáng)作用下,沉沙量會(huì)增加,并且使沙顆粒的沉積面略微傾斜。油或水相中沙顆粒的起揚(yáng)速度隨沙顆粒大小的分布曲線變化,直徑為0.162 mm的沙顆粒,在水或油相中的起揚(yáng)速度在0.04 m/s左右,說(shuō)明了圖6中出現(xiàn)的現(xiàn)象。

圖6 模型1彎曲管道內(nèi)固液多相流中沙顆粒體積分布

3 結(jié) 論

(1) 同一流速和含沙量,在管道下降段油相體積分布區(qū)域大于上升段的油相體積,水相則相反。即在下降段油相速度明顯小于入口流速,而下降段水相速度明顯大于入口流速,上升段與此相反。

(2) 對(duì)于同一含沙量中,由于軸向速度分布相差不大,隨著流體速度的增加油相體積分布區(qū)域變小,水相體積則相反。

(3) 水平挾沙油水的固液多相流中,沙顆?；旧隙汲恋碓趶澢艿赖牡撞?并且隨著流速增大,在下凹彎曲段沉沙量先減少后增多并向前推進(jìn)。

(4) 對(duì)于同一沙顆粒直徑,在油相中的揚(yáng)起速度略大于水相中的揚(yáng)起速度,并且隨著顆粒直徑的增大,揚(yáng)起速度變大。

參考文獻(xiàn):

[1] 徐孝軒, 宮敬. 水平管中油氣水三相流動(dòng)研究進(jìn)展 [J]. 化工機(jī)械, 2005, 32(6): 329335, 377.

[2] 賈月梅. 流體力學(xué) [M]. 北京: 國(guó)防工業(yè)出版社, 2006.

[3] 陳杰. 水平油水管流流動(dòng)規(guī)律研究 [D]. 北京: 中國(guó)石油大學(xué)(北京), 2001.

[4] CHARLES M E. Water Layer Speeds Heavy-Crude Flow [J]. The Oil and Gas Journal, 1961, 28: 68-72.

[5] DIMITRI GORELIK, BRAUNER N. The Interface Configuration in Two-Phase Stratified Pipe Flow [J]. Int. J. Multiphase Flow, 1999, 25: 977-1007.

[6] 郭烈錦. 兩相與多相流動(dòng)力學(xué) [M]. 西安: 西安交通大學(xué)出版社, 2002.

[7] BRAUNER N, ROVINSKY J, MOALEM MARON. Determination of the Interface Curvature in Stratified Two-Phase Systems by Energy Considerations [J]. Int. J. Multiphase Flow, 1996, 22(6): 116-1185.

[8] BRAUNER N, MOALEM MARON, ROVINSKY J. A Two-Fluid Model for Stratified Flows With Curved Interfaces [J]. Int. J. Multiphase Flow, 1998, 24: 975-1004.

[9] 陳鑫, 魯傳敬, 李杰, 等. VOF和Mixture多相流模型在空泡流模型中的應(yīng)用 [C]∥第九屆全國(guó)水動(dòng)力學(xué)學(xué)術(shù)會(huì)議論文集, 2009.

[10] 謝艷芳. 多相流混合模型應(yīng)用于低含沙水流的數(shù)值模擬研究 [D]. 西安: 西安理工大學(xué), 2005.

[11] 王濤, 李紅, 馮世峰, 等. 基于Mixture模型對(duì)自吸噴灌泵自吸過(guò)程的氣液兩相流模擬 [C]∥現(xiàn)代節(jié)水高效農(nóng)業(yè)與生態(tài)灌區(qū)建設(shè)會(huì)議論文集(上), 2010.

[12] 劉亞青. 水平井氣水兩相流數(shù)值試井理論研究與應(yīng)用 [D]. 成都: 西南石油大學(xué), 2011.

[13] 張宏兵, 陳露露, 謝榮華, 等. 水平圓管固液兩相穩(wěn)態(tài)流動(dòng)特性數(shù)值模擬 [J]. 化工學(xué)報(bào), 2009, 05: 1162-1168.

[14] 邵艷宏, 陳仙江, 張紹輝, 等. 數(shù)值試井的研究現(xiàn)狀及展望 [J]. 內(nèi)蒙古石油化工, 2009, 35(21): 26-28.

[15] 宋紅偉, 郭海敏, 唐小梅. 水平井及大斜度井油水兩相流動(dòng)模型研究 [J]. 測(cè)井技術(shù), 2014, 38(4): 384-390.

[16] 周新民, 趙文謙, 李嘉. 利用相間滑移考慮固相泥沙沉降效應(yīng) [J]. 水動(dòng)力學(xué)研究與進(jìn)展(A輯), 2004, 19(5): 53-5978.