張佳沁 劉軍鋒* 邢廣俊 史圣哲 王歡歡
(1、油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室(長(zhǎng)江大學(xué)),湖北 武漢430100 2、長(zhǎng)江大學(xué)地球物理與石油資源學(xué)院,湖北 武漢430100)
同時(shí)存在兩種物質(zhì)狀態(tài)的兩相流在自然界和現(xiàn)代科學(xué)技術(shù)領(lǐng)域中是極為常見的,例如石油工程中氣水、油水兩相流。石油開發(fā)時(shí),從井底到井口,井筒和地層之間有多個(gè)射孔層或生產(chǎn)層進(jìn)行連通,地層中的油氣水在壓力差的作用下,逐漸流入到井筒內(nèi)。井筒中常見有油水、氣水、油氣兩相或油氣水三相流。有時(shí),氣水、油氣和油氣水統(tǒng)稱為氣液兩相,油水稱為液液兩相。由于井筒傾斜角度(垂直、傾斜或水平)、幾何尺寸、材質(zhì)與兩相流物性參數(shù)、流量等因素的影響,兩相流的流動(dòng)狀態(tài)(即流型)并非始終保持不變,流動(dòng)狀態(tài)發(fā)生改變使得兩相流的流型隨之轉(zhuǎn)變。
1987 年,Barnea 用空氣-水在管徑為5.1cm 管道中進(jìn)行實(shí)驗(yàn)得到流動(dòng)的流型和過度邊界,將流型分為分層流( SS)、波狀分層流(SW)、環(huán)狀流(A)、變形泡狀流(EB)、段塞流(SL)、沫狀流(CH)和分散泡狀流(DB)等,最終模型結(jié)果與實(shí)驗(yàn)數(shù)據(jù)比較良好。
Brill 用空氣-煤油、空氣-潤(rùn)滑油在管徑為38.1mm 的管徑中進(jìn)行角度對(duì)流型轉(zhuǎn)變實(shí)驗(yàn),得到以無因次準(zhǔn)數(shù)表達(dá)的氣液兩相流的流型判別式。由于在小管徑中流型劃分較多差異較少不易區(qū)分,Brill 將流型劃分為氣團(tuán)流、段塞流、環(huán)狀流和分層流這四類。
1958 年,Russell 在水-礦物油的實(shí)驗(yàn)中,觀察到了分散流、分層流和混合流這三種流型。1989 年,Arirachakaran 在直徑為25.1mm 的管道實(shí)驗(yàn)中,觀察到了分層流、混合流、環(huán)狀流、間歇流和分散流這五種流型。1996 年,Trallero 在水-礦物油兩相流實(shí)驗(yàn)中,加上對(duì)之前流型劃分的總結(jié),將水平井油水兩相流分為六種流型,分別為分層流(ST)、界面混合分層流(ST&MI)、油包水-水包油(DW/O&DO/W)、油包水(O/W)、水包油-水(DO/W&W)、水包油(W/O),該劃分在學(xué)術(shù)上被廣為認(rèn)同。
1999 年,F(xiàn)lores 在直徑為5.08cm 傾斜管道內(nèi)進(jìn)行油水兩相流實(shí)驗(yàn),將油水兩相流分為幾種流型,分別為段塞水包油(DO/WPS)、局部逆流水包油(DO/WCT)、平行流水包油(DO/WCC)、油相乳化水包油(VFDO/W)、油包水(DW/O)和水相乳化油包水(VFDW/O)。
本文以直徑20mm 的水平管油水兩相流為例,通過將物理模擬實(shí)驗(yàn)得到的各相流速和流型等資料進(jìn)行擬合,即可得到流型判別式。
(1)分層流與界面混雜分層流的判別式
油的表觀速度在0.02m/s 到0.2m/s 范圍內(nèi),水的表觀速度在0.02m/s 到0.14m/s 的范圍內(nèi)。
(2)界面混雜分層流的上邊界
油的表觀速度在0.02m/s 到0.3m/s 的范圍內(nèi),水的表觀速度在0.03m/s 到0.3m/s 的范圍內(nèi)。
(3)三層流與油包水-水包油的判別式
油的表觀速度在0.31m/s 到0.72m/s 的范圍內(nèi),水的表觀速度在0.07m/s 到0.3m/s 的范圍內(nèi)。
(4)油包水-水包油與油包水之間的判別式
油的表觀速度在0.3m/s 到0.5m/s 的范圍內(nèi),水的表觀速度在0.02m/s 到0.12m/s 的范圍內(nèi)。
(5)三層流與水包油-水之間的判別式
油的表觀速度在0.08m/s 到0.41m/s 的范圍內(nèi),水的表觀速度在0.15m/s 到0.28m/s 的范圍內(nèi)。
(6)水包油與水包油-水之間的判別式
油的表觀速度在0.09m/s 到0.51m/s 的范圍內(nèi),水的表觀速度在0.14m/s 到0.67m/s 的范圍內(nèi)。
式中,VSO、VSW分別是油相、水相的表觀速度,m/s。
從調(diào)研資料中發(fā)現(xiàn):影響流型的關(guān)鍵因素很多,主要有管道傾角、管道材質(zhì)、管道尺寸,氣液兩相流的表觀流速,流體密度、粘度、表面張力等,所以在進(jìn)行流型自動(dòng)判別時(shí),要考慮到上述因素,判別思路如下:
(1)輸入已知條件,氣液兩相的密度,液相粘度,管道材質(zhì),管道尺寸,管道傾角等;
(2)將已知的特征參數(shù)帶入設(shè)定的流型判別公式中,得到有關(guān)氣液兩相表觀流速的判別式子;
(3)輸入需要判別的氣液兩相的表觀流速;
(4)判斷輸入的氣液兩相的流速能否識(shí)別流型,若輸入數(shù)據(jù)在識(shí)別范圍內(nèi),則輸出流型;若輸入數(shù)據(jù)在識(shí)別范圍外,則不能得到結(jié)果;
(5)重新輸入數(shù)據(jù),進(jìn)行流型判別。
在不同的環(huán)境中,氣液兩相的特征參數(shù)也會(huì)發(fā)生相應(yīng)的改變。為了使此判別法有更廣泛的適用性,制作數(shù)據(jù)輸入窗口時(shí)將這些因素也當(dāng)作變量顯示出來以便能讓使用者根據(jù)實(shí)際情況調(diào)整輸入數(shù)據(jù)。
利用Excel 制作簡(jiǎn)易的數(shù)據(jù)輸入窗口,如圖1 所示,將氣液兩相的密度,液相粘度,管道的材質(zhì)、尺寸、傾角,氣液兩相的表觀速度作為輸入?yún)?shù),通過流型劃分公式調(diào)用這些數(shù)據(jù)進(jìn)行計(jì)算即可實(shí)現(xiàn)準(zhǔn)確判別流型的目的。
本文采用輸出流型圖及點(diǎn)位的方式來直觀展示流型判別結(jié)果,綜合本文上述已有資料,將流型大致劃分為分散泡狀流(Dispersed-Bubble)、泡狀流(Bubble)、光滑分層流(Stratified)、波狀分層流(Wavy)、環(huán)狀流(Annular),借助擬合所得出的流型劃分邊界條件,同時(shí)輸入具體數(shù)據(jù)來檢驗(yàn)流型能否進(jìn)行自動(dòng)判別。
3.3.1 水平井氣液兩相流型判別
在管徑為24mm,傾斜角為0 度的有機(jī)玻璃管進(jìn)行水平井氣液兩相流型判別的檢驗(yàn),利用前人真實(shí)實(shí)驗(yàn)數(shù)據(jù)及實(shí)驗(yàn)觀測(cè)到的流型進(jìn)行判別,判別結(jié)果如圖2 所示,通過比較,二者大致吻合,說明能夠進(jìn)行流型判別。
3.3.2 斜井氣液兩相流型判別
在24mm 管徑,60°傾斜角的有機(jī)玻璃管中,將預(yù)先準(zhǔn)備好的已知的多組流型數(shù)據(jù)輸入,判別結(jié)果如圖3 所示,可見判別結(jié)果與已知流型數(shù)據(jù)大致吻合。
3.3.3 流型轉(zhuǎn)換規(guī)律的分析
由流型判別圖可見,當(dāng)外部參數(shù)如管徑,傾角一定時(shí),隨著氣液兩相流速的增加,管內(nèi)的流型也隨之改變。當(dāng)氣液兩相流速都很小時(shí),管內(nèi)出現(xiàn)光滑分層流。當(dāng)氣相流速相較于液相流速快速增加時(shí),流型從光滑分層流向波狀分層流進(jìn)而向環(huán)狀流轉(zhuǎn)化。當(dāng)液相流速相較于氣相流速快速增加時(shí),管內(nèi)流型也由分層流、環(huán)狀流向分散氣泡流、泡狀流轉(zhuǎn)化。這與理論研究和實(shí)驗(yàn)現(xiàn)象都相符合。
圖2 水平管氣液兩相流型判別圖
圖3 傾斜管氣液兩相流型判別圖
4.1 本文根據(jù)實(shí)驗(yàn)所得流型數(shù)據(jù)擬合所得的流型判別公式來實(shí)現(xiàn)流型判別,從而規(guī)避了傳統(tǒng)識(shí)別方法具有較強(qiáng)的主觀性的缺點(diǎn),使判別結(jié)果較為客觀。利用Excel 編程,進(jìn)行界面交互,實(shí)現(xiàn)人機(jī)對(duì)話。在流型識(shí)別上能夠準(zhǔn)確高效的進(jìn)行,從而達(dá)到流型在線識(shí)別的目的。
4.2 本文實(shí)現(xiàn)氣液兩相流型的自動(dòng)識(shí)別,根據(jù)管道傾斜角度、幾何尺寸、材質(zhì)、兩相流物性參數(shù)、各相速度的不同可得到不同的流型,能夠?qū)Σ煌瑺顟B(tài)下的流型進(jìn)行觀察比較得到不同因素對(duì)流型轉(zhuǎn)化的影響。