王修武，羅 威，劉 捷，廖銳全，陳元虎

(1.長江大學(xué)，湖北 武漢 430100；2.中國石油氣舉試驗基地多相管流實驗室，湖北 武漢 430100；3.中國石化江蘇油田分公司，江蘇 揚州 225000)

0 引 言

多相管流的研究已經(jīng)有30多年的歷史，研究得到的多相管流預(yù)測方法有二十幾種，被廣泛應(yīng)用的多相流計算方法接近10種。多相管流預(yù)測方法按導(dǎo)出原理可分為半經(jīng)驗、經(jīng)驗?zāi)Ｐ蚚1-5]和機理模型[6-7]。盡管如此，多相管流還存在很多不完善和亟待解決的實際問題。早前，由于實驗條件有限(介質(zhì)、流量范圍、測量儀器精度[8]等)，能在油、氣、水三相流基礎(chǔ)上得到的多相流計算方法較少。近年來，中國著名科研院所開始側(cè)重于稠油多相流動方面的研究[9]，直接忽略了對常規(guī)黏度原油與氣、水三相流動的系統(tǒng)性研究。利用國外流行工程計算軟件PIPSIM、中國編制多相流計算軟件對中小產(chǎn)量、較高氣油比油井進行多相流預(yù)測，與實測數(shù)據(jù)對比表明，目前已有多相流計算模型預(yù)測壓力結(jié)果普遍偏高；而對中高產(chǎn)量、較高氣油比油井進行多相流預(yù)測，與油田實測數(shù)據(jù)對比表明，目前已有多相流計算模型預(yù)測壓力結(jié)果普遍偏低[10]。針對以上問題，目前缺少常規(guī)黏度原油、氣、水三相流動的系統(tǒng)性研究[11-12]。鑒于此，開展了常規(guī)黏度油、氣、水三相流動的系統(tǒng)性實驗研究，并在此基礎(chǔ)上開展了多相管流涉及到的流型、持液率、壓降3個方面的預(yù)測模型研究。

1 流型預(yù)測模型驗證及優(yōu)選

在多相管流流型圖研究方面，目前比較認(rèn)可和備受關(guān)注的多相管流流型預(yù)測方法有Zhang Hongquan[6]、L.E.Gomez[13]、Kaya[14]、J.J.Xiao[15]、Beggs-Brill[5]、Barnea[16]、Mukherjee-Brill[8]。其中，Zhang Hongquan和Kaya屬于流型預(yù)測機理模型，L.E.Gomez和J.J.Xiao對流型的研究較少，Beggs-Brill、Barnea、Mukherjee-Brill屬于流型預(yù)測，非機理模型。鑒于此，開展了油、氣、水三相流研究，實驗?zāi)M參數(shù)如表1所示。將通過實驗繪制流型圖版(高速攝像機拍攝和肉眼觀察相結(jié)合)與現(xiàn)有模型繪制流型圖版進行對比，如圖1所示。Zhang Hongquan模型和Beggs-Brill模型的實驗條件均為51 mm管徑，空氣或水為介質(zhì)，傾角為20 °。

表1 實驗?zāi)M參數(shù)選擇

圖1 不同含水及傾角多相流實驗繪制流型圖與現(xiàn)有流型圖對比

由圖1可知，流型預(yù)測方法存在不準(zhǔn)確的地方。從總體上看，Mukherjee-Brill模型預(yù)測流型效果最好，Zhang Hongquan模型在傾斜角為60～90 °時流型預(yù)測效果較好，Barnea模型水平流型預(yù)測較好。因此，優(yōu)選Mukherjee-Brill流型預(yù)測模型作為研究預(yù)測流型。

2 持液率和壓降預(yù)測模型驗證及擬合

2.1 持液率和壓降預(yù)測模型準(zhǔn)確性驗證

在多相管流持液率和壓降預(yù)測方面，首先采用實驗數(shù)據(jù)對常用的6種多相流計算方法進行壓降預(yù)測驗證(表2)，并對其中2種具有持液率計算能力的方法(Beggs-Brill、Mukherjee-Brill)進行持液率預(yù)測驗證(表3)。

2.2 不同組合條件下持液率計算誤差分析

對上述2種具有持液率計算能力的多相流方法分別進行4種不同組合條件下的壓降預(yù)測對比(表4、5)。其中，4種不同組合條件為：A為持液率、流型和壓降均預(yù)測，B為持液率給出、流型和壓降均預(yù)測，C為流型給出、持液率和壓降預(yù)測，D為流型和持液率給出、壓降預(yù)測。

表2 6種方法壓降預(yù)測平均誤差統(tǒng)計

表3 2種方法持液率預(yù)測平均誤差統(tǒng)計

由表4、5可知，不同組合方法下，當(dāng)給出持液率時，與其他組合條件下相比，壓降預(yù)測準(zhǔn)確率大大提高。這也說明，準(zhǔn)確預(yù)測持液率能大大提高傾斜管流壓降預(yù)測的精度。

表4 Beggs-Brill方法不同組合條件壓降預(yù)測平均誤差統(tǒng)計

表5 Mukherjee-Brill方法不同組合條件壓降預(yù)測平均誤差統(tǒng)計

同時，對Beggs-Brill方法和Mukherjee-Brill方法進行不同組合條件下持液率預(yù)測檢驗，結(jié)果如表6、7所示。

表6 Beggs-Brill方法不同組合條件持液率預(yù)測平均誤差統(tǒng)計

由表6、7可知，Beggs-Brill方法流型預(yù)測與實際有較大差別，傾斜段預(yù)測存在較大偏差，因此，給出流型時，持液率計算誤差差別較大，而Mukherjee-Brill則基本相同，進一步驗證該模型流型圖判斷的準(zhǔn)確性。

表7 Mukherjee-Brill方法不同組合條件持液率預(yù)測平均誤差統(tǒng)計

2.3 持液率預(yù)測模型擬合

鑒于持液率的預(yù)測準(zhǔn)確性對壓降預(yù)測影響較大，準(zhǔn)確預(yù)測持液率能夠大幅度提高壓降預(yù)測準(zhǔn)確性?；诔Ｒ?guī)黏度油氣水三相多相管流系統(tǒng)性實驗數(shù)據(jù)(40、60、75 mm 3種管徑、不同傾角、不同含水率條件)對Mukherjee-Brill[14]的持液率計算結(jié)果進行了重新擬合。由表7可知，水平狀態(tài)和傾斜狀態(tài)計算誤差差別較大，因此，對水平和傾斜2種條件分別進行了擬合(表8)。

3 新預(yù)測模型驗證和誤差分析

分別從實驗數(shù)據(jù)、現(xiàn)場實測數(shù)據(jù)對新預(yù)測模型進行驗證。由于新模型是在系統(tǒng)性實驗數(shù)據(jù)基礎(chǔ)上擬合得到，因此，首先需要通過實驗數(shù)據(jù)對新模型進行檢驗，驗證其與實驗數(shù)據(jù)的符合情況，即驗證擬合過程的正確性，再通過現(xiàn)場實測數(shù)據(jù)對新預(yù)測模型進行檢驗。

3.1 實驗數(shù)據(jù)對新模型的檢驗

持液率和壓降預(yù)測結(jié)果與實驗數(shù)據(jù)之間平均誤差如表9所示。

表8 水平狀態(tài)、斜井狀態(tài)持液率擬合結(jié)果

表9 新模型預(yù)測結(jié)果與實驗數(shù)據(jù)平均誤差統(tǒng)計

將表9與表2、3、5、7對比可知，持液率預(yù)測和壓降預(yù)測精度均有較大提升，在傾斜狀態(tài)(包括垂直狀態(tài))可以滿足工程計算預(yù)測需要，證明擬合方法過程可行；由于水平管流壓降中摩阻項占主導(dǎo)，持液率重新擬合后壓降預(yù)測精度提升較小，而傾斜管流壓降中重力項占主導(dǎo)，持液率預(yù)測準(zhǔn)確后該項預(yù)測精度提升幅度較大。

3.2 現(xiàn)場實測數(shù)據(jù)對新模型檢驗

以新疆瑪湖油田多口生產(chǎn)井的實測流壓數(shù)據(jù)為例，對現(xiàn)有常用的6種井筒多相管流計算方法和新模型進行壓力預(yù)測驗證。

3.2.1 實測數(shù)據(jù)

統(tǒng)計瑪湖地區(qū)7井次(5口井，測試井套管管徑為139.7 mm，油管平均傾斜角為0 °)的測試數(shù)據(jù)(表10)。

3.2.2 預(yù)測結(jié)果誤差分析

通過實測數(shù)據(jù)對7種多相管流壓力預(yù)測方法(包含新模型)進行檢驗，計算結(jié)果平均誤差統(tǒng)計如表11所示。

由表11可知，新預(yù)測模型平均相對誤差最小，為5.08%，其次為Mukherjee-Brill[14]方法，為8.53%?，F(xiàn)場實測數(shù)據(jù)驗證表明，以實驗數(shù)據(jù)為基礎(chǔ)建立的多相管流預(yù)測模型可降低多相管流壓降預(yù)測誤差，新預(yù)測模型可用于現(xiàn)場多相管流壓力分布預(yù)測。

4 結(jié) 論

(1) 在多相管流流型劃分方面，結(jié)合實驗對流型圖版對比發(fā)現(xiàn)，最適合的流型判斷方法為Mukherjee-Brill方法。

(2) 實驗數(shù)據(jù)檢驗表明，現(xiàn)有6種多相流壓降預(yù)測方法誤差均較大，準(zhǔn)確預(yù)測持液率能極大提高傾斜管流壓降預(yù)測的準(zhǔn)確性，進一步說明Mukherjee-Brill方法流型預(yù)測的準(zhǔn)確性。

(3) 對Mukherjee-Brill持液率方法進行了重新擬合，檢驗表明，新模型擬合效果較好，且其傾斜管流壓降預(yù)測誤差明顯降低?，F(xiàn)場實測數(shù)據(jù)驗證表明，新預(yù)測模型較其他模型預(yù)測精度至少提高3個百分點，預(yù)測更加可靠。

(4) 研究表明，傾斜管流壓降中重力項占主導(dǎo)，持液率預(yù)測準(zhǔn)確后，壓降預(yù)測精度提升幅度較大，而由于水平管流壓降中摩阻項占主導(dǎo)，持液率重新擬合后壓降預(yù)測精度提升較小。

表10 7井次(5口井)測試數(shù)據(jù)統(tǒng)計

表11 壓力計算平均相對誤差統(tǒng)計(平均7井次)

[1] DUNS H J，ROS N C J.Vertical flow of gas and liquid mixtures in wells[C].WPC10132，1963：451-466.

[2] HAGEDORN A R，BROWN K E.Experimental study of pressure gradients occurring during continuous two-phase flow in small-diameter vertical conduits[J].Journal of Petroleum Technology，1965，17(4)：475-484.

[3] ORKISZEWSKI J.Predicting two-phase pressure drops in vertical pipe[J].Journal of Petroleum Technology，1967，19(6)：829-838.

[4] AZIZ K，GOVIER G W.Pressure drop in wells producing oil and gas[J].Journal of Canadian Petroleum Technology，1972，11(3)：38-48.

[5] BEGGS D H，BRILL J P.An experimental study of two-phase flow in inclined pipes[J].Journal of Petroleum Technology，1973，25(5)：607-617.

[6] WANG S，PEREYRA E，SARICA C，et al.A mechanistic slug liquid holdup model for wide ranges of liquid viscosity and pipe inclination angle[C].OTC-24046，2014：329-336.

[7] WANG S，ZHANG H Q，SARICA C，et al.A mechanistic slug-liquid-holdup model for different oil viscosities and pipe-inclination angles[J].SPE Production & Operations，2014，29(4)：329-336.

[8] 李穎川.采油工程[M].北京：石油工業(yè)出版社，2002：41-43.

[9] NOSSEN J，LAWRENCE C J.Pressure drop in laminar slug flow with heavy oil[C].BHR-A006，2012：91-104.

[10] WANG S.Experiments and model development for high-viscosity oil/water/gas horizontal and upward vertical pipe flows[D].Oklahoma：University of Tulsa，2012.

[11] 張榮軍，陳軍斌，任月玲，等.水平井筒中流型識別的數(shù)學(xué)方法[J].特種油氣藏，2009，16(2)：67-69.

[12] 丁景辰.致密氣藏同步回轉(zhuǎn)增壓排水采氣工藝[J].特種油氣藏，2017，24(3)：145-149.

[13] GOMEZ L，SHOHAM O，SCHMIDT Z，et al.A unified mechanistic model for steady-state two-phase flow in wellbores and pipelines[C].SPE56520，1999：1-14.

[14] KAYA A S，SARICA C，BRILL J P.Mechanistic modeling of two-phase flow in deviated wells[J].SPE Production & Facilities，1999，16(3)：156-165.

[15] XIAO J J，SHOHAM O，BRILL J P.A comprehensive mechanistic model for two-phase flow in pipelines[C].SPE20631，1990：167-180.

[16] BARNEA D.A unified model for predicting flow-pattern transitions for the whole range of pipe inclinations[J].International Journal of Multiphase Flow，1987，13(1)：1-12.