趙忠武 劉捷

摘 ?????要： CO₂驅(qū)作為一項具有前景的工藝技術(shù)，既能補充地層能量，又能降低原油的黏度，相比注水開發(fā)具有獨特優(yōu)勢。但由于CO₂對原油物性參數(shù)有較大影響，使得井筒中氣液兩相流動規(guī)律發(fā)生變化，不同于常規(guī)天然氣和原油體系的流動規(guī)律。在CO₂-原油物性計算模型分析的基礎(chǔ)上，建立了高含CO₂油井的節(jié)點分析模型，考慮高CO₂含量下的相態(tài)特點，在幾種不同出口壓力條件下，針對新疆油田不同CO₂含量和不同含水率情況下的停噴壓力進(jìn)行了預(yù)測。分析結(jié)果表明，停噴壓力隨含水率上升而增大，隨CO₂含量上升而減小。

關(guān) ?鍵 ?詞：二氧化碳;節(jié)點分析;狀態(tài)方程;停噴壓力

中圖分類號：TE 357????????文獻(xiàn)標(biāo)識碼：?A ??????文章編號： 1671-0460（2020）01-0125-04

Analysis on the Timing of Lifting Process for Wells With High CO₂?Content

ZHAO?Zhong-wu， LIU?Jie

（School of Petroleum Engineering， Yangtze University， Hubei Wuhan 430100，?China）

Abstract： As?a promising process，?CO₂?flooding can not only supply the formation energy， but also reduce?the?viscosity of crude oil.?It has more advantages comparing?with the development of water injection.?However， due to?the influence of CO₂?on the physical parameters of the crude oil， the gas-liquid two-phase flow pattern changes in the wellbore， which is different from the flow rules of conventional natural gas and crude oil systems.On the basis of the analysis of the physical property calculation model of CO₂-crude oil， a node analysis model for oil wells?with high CO₂?content was?established. Considering the phase characteristics under the high CO₂?content， the unflowing pressure of oil well with different CO₂?content and different water content?in Xinjiang oilfield?was?predicted under different?outlet pressure. The analysis results?showed?that?the unflowing pressure?increased?with the increase of water content， and it decreased?with the increase of CO₂?content.

Key words： Carbon dioxide; Node analysis; Equation of state; Unflowing pressure

二氧化碳驅(qū)是當(dāng)前世界范圍內(nèi)一項較為成熟、高效的采油技術(shù)，相比傳統(tǒng)水驅(qū)油具有獨特的技術(shù)優(yōu)勢^[1]，因此具有廣闊的應(yīng)用前景，也一直是我國的研究重點。國內(nèi)外大量學(xué)者對CO₂-原油體系的基本物性進(jìn)行了研究^[2-25]，結(jié)果表明，CO₂在原油中具有很好的溶解性，且CO₂在原油中的溶解會導(dǎo)致原油黏度、體積系數(shù)以及密度等物性參數(shù)發(fā)生變化，變化規(guī)律與天然氣和原油體系表現(xiàn)出來的變化規(guī)律存在明顯差異，這就使得常用多相管流計算方法無法適應(yīng)CO₂-原油-水流動過程中壓降的計算。

因此，從CO₂-原油混合物相平衡計算出發(fā)，在確定高含CO₂井筒壓力計算方法的基礎(chǔ)上，建立了高含CO₂油井的節(jié)點分析模型，并對油井的停噴壓力進(jìn)行了預(yù)測。

1 ?建立模型

1.1 ?狀態(tài)方程優(yōu)選

Phillips^[26]（1912）首次提出基于狀態(tài)方程建立黏度模型的可行性，隨后大量研究者開發(fā)出基于不同狀態(tài)方程的黏度模型。常用的狀態(tài)方程有PR^[27]、SRK^[28]、BWRS^[29]這三種，這些狀態(tài)方程半理論、半經(jīng)驗的關(guān)系式，都有一定的使用范圍，為了能夠準(zhǔn)確進(jìn)行井筒壓降的計算，首先根據(jù)CO₂-原油體系室內(nèi)實驗數(shù)據(jù)，對三種狀態(tài)方程進(jìn)行評估，選擇計算精度最高的狀態(tài)方程作為后期井筒CO₂-原油-水體系壓降計算的基礎(chǔ)。

依據(jù)新疆油田現(xiàn)場得到的混合相黏度數(shù)據(jù)資料對狀態(tài)方程進(jìn)行選擇，分別計算出不同狀態(tài)方程下對應(yīng)的混合相黏度值，然后與實驗數(shù)據(jù)進(jìn)行誤差分析，以誤差最小原則優(yōu)選適合該區(qū)塊的狀態(tài)方程。表1計算結(jié)果表明，PR狀態(tài)方程平均相對誤差1.67%，是三個狀態(tài)方程中誤差最小的，因此，節(jié)點分析過程中油氣體系物性參數(shù)的計算時狀態(tài)方程選PR方程。

對于多組分體系，PR狀態(tài)方程的形式：

式中：P?—壓力，MPa;

R?—氣體通用常數(shù)，0.008?314?MPa×m³/（kmol×K）;

T?—溫度，K;

V?—比容，m³/kmol;

X_i、X_j、n分別表示組分i和組分j的摩爾分?jǐn)?shù)以及組分個數(shù);

k_ij?—二元交互作用系數(shù);

m?—混合物。

1.2 ?考慮相態(tài)的多相管流計算模型

CO₂會對原油黏度、體積系數(shù)和密度等參數(shù)造成影響，因此必須結(jié)合相平衡和狀態(tài)方程，對多相管流模型進(jìn)行修正才能運用于高含CO₂的油井。國內(nèi)外多相管流計算方法很多，對模型的選擇是針對具體的油氣藏條件，根據(jù)經(jīng)驗和對比計算來選擇固定的計算方法，本次研究依據(jù)經(jīng)驗選用Hagedorn- Brown方法，Hagedorn-Brown^[30]導(dǎo)出的多相管流壓降公式：

式中： ?p?—壓力增量，MPa;

?H?—深度增量，m;

ρ_m?—氣液混合物密度，kg/m³;

f_m?—兩相摩阻系數(shù);

q_l?—地面產(chǎn)液量，m³/d;

M_t?—地面標(biāo)準(zhǔn)條件下每生產(chǎn)1 m³氣液混合物的總質(zhì)量，kg/m³;

d?—油管內(nèi)徑，m;

v_m?—氣液混合物速度，m/s。

為了使Hagedorn-Brown方法適合高含CO₂油井的多相流計算，相關(guān)式中的各參數(shù)要考慮井筒內(nèi)CO₂-原油體系的相態(tài)變化特點，確定如下：

（1）混合物的單位體積質(zhì)量

M_t?=1?000g ₀?+?1.205GORg_g

式中： ?GOR?—氣油比;

g_o?—油相對密度;

g_g?—CO₂相對密度。

（2）氣相密度、液相密度、氣液混合物密度

式中：?Z_g?—氣體壓縮因子;

Z_o?—原油體積系數(shù);

y_i?—組分i的摩爾分?jǐn)?shù)，無因次;

M_i?—組分i的質(zhì)量流量，kg/m³;

H_l?—持液率，無因次。

（3）氣相速度、液相速度、氣液混合物速度

式中：`T —平均溫度，K;`

P —平均壓力，MPa;

A —橫截面積，m²。

（4）氣液混合物黏度

式中：?m_m?—混合物黏度，mPa×s;

m_o?—原油黏度，mPa×s;

m_g?—氣相黏度，mPa×s。

1.3 ?流入動態(tài)模型

流入動態(tài)模型，即流入動態(tài)方程，是在一定的油氣壓力下，流體產(chǎn)量與相應(yīng)井底壓力的關(guān)系方程，流入動態(tài)計算方程最常見的有四種：Vogel方法、Jones方法、Fetkovich方法、采油指數(shù)方法。每一種方法都有各自的應(yīng)用范圍，沒有一種方法能夠適合所有的油井，某方法在一口油井中預(yù)測精度可靠，但是在另一口井中預(yù)測精度可能是最差的。本文為了便于計算，選用采油指數(shù)方程，模型數(shù)學(xué)形式為：

式中：??J —采油指數(shù)，m³/（d×MPa）;

P_r?—地層壓力，MPa;

P_wf?—井底壓力，MPa。

2 ?高含CO₂油井停噴預(yù)測

過去預(yù)測油氣藏的廢棄地層壓力最常用的方法是將同一相似氣藏的開采動態(tài)相比較，通過類比分析和經(jīng)驗公式來計算，由于計算機技術(shù)和多相流計算方法的發(fā)展，越來越多的工程師應(yīng)用更加精確地分析技術(shù)來預(yù)測停噴壓力-這就是節(jié)點分析技術(shù)（SNA）。

這種方法是繪制出不同流量和儲層靜壓下的流入、流出動態(tài)曲線，然后進(jìn)行分析。停噴壓力是指在該壓力下流入動態(tài)曲線和油管動態(tài)曲線不再相交，如圖1所示。流入動態(tài)曲線和油管動態(tài)曲線的交點即為一定井口壓力下定產(chǎn)量生產(chǎn)油井可能生產(chǎn)的最小地層壓力。

以新疆油田A區(qū)的油井為例，在油井含水10%、含CO₂分別為18.58%和51.86%的情況下，對1～4 MPa不同井口壓力時的停噴壓力進(jìn)行預(yù)測，預(yù)測過程中由于CO₂的百分含量無法直接參與計算，因此轉(zhuǎn)化為混合相中的氣液比來進(jìn)行計算，計算結(jié)果如圖2、圖3所示。

根據(jù)節(jié)點分析方法，圖2中井口壓力1～4 MPa對應(yīng)的停噴壓力分別為20、22、23、24 MPa;圖3中井口壓力1～4 MPa對應(yīng)的停噴壓力分別為16、21、23、24 MPa。綜合分析可以發(fā)現(xiàn)，隨著井口壓力的降低，停噴壓力也在降低，這主要是由于在保持產(chǎn)量一定的情況下，井底流壓會同時降低，最終使得停噴壓力變小。

二氧化碳含量越高，也就是氣液比越大，停噴壓力值越低，尤其是當(dāng)出口壓力較低時體現(xiàn)得更為明顯，因為隨著二氧化碳含量的增加，原油中溶解的二氧化碳?xì)怏w越多，原油的黏度越低，流動性能越好，在孔隙中流動阻力越小，并且大量溶解的二氧化碳?xì)怏w會使原油體積發(fā)生膨脹，膨脹能使得油井自噴能量增大。

3 ?停噴壓力的影響因素分析

3.1 ?CO₂含量對停噴壓力的影響

對比圖2和圖3可以看出，CO₂含量越高，也就是氣液比越大，停噴壓力值越低，為了能更進(jìn)一步分析CO₂含量的影響，同樣以含水10%的油井為例，用節(jié)點分析的方法確定不同CO₂含量時的停噴壓力，停噴壓力隨CO₂含量的變化趨勢如圖4所示。

從圖形中可以直觀的看出，隨著氣液比的增加，也就是二氧化碳含量的增加，初期停噴壓力保持不變，當(dāng)氣液比增大到一定程度后，油井的停噴壓力開始下降，井口壓力越大下降趨勢越平緩，井口壓力越小下降趨勢越明顯。初期由于氣液比較小，管流過程中能量的損失以滑脫為主，雖然氣液比增大，但對能量損失的影響不大，因此停噴壓力基本不變;氣液比增大到一定程度，管流過程中的能量損失就以摩擦為主，隨氣液比的增大，流動過程中的壓力損失降低，使得停噴壓力也會隨之降低。

3.2 ?含水率對停噴壓力的影響

由于CO₂在原油和水中溶解度的不同，含水率必然會影響油水混合物的物性參數(shù)，進(jìn)而影響油井的停噴壓力，如圖5-8所示，反映了不同井口壓力條件下含水率對停噴壓力的影響。

圖5-8中停噴壓力曲線都是呈上升趨勢的，并且上升的幅度越來越大，這說明不管在哪種出口壓力條件下，停噴壓力都隨含水率的上升而增大，而且含水率越高，停噴壓力增大的也越快。

由于二氧化碳在水中的溶解度遠(yuǎn)遠(yuǎn)小于在原油中的溶解度，所以含水率上升后，流體中溶解的二氧化碳含量降低，氣體膨脹能減小，流動過程中的各種壓降損失增大，就造成了停噴壓力隨含水率增加的趨勢;圖中還有一個共同點就是高二氧化碳含量對應(yīng)著低停噴壓力曲線，這是由于二氧化碳含量增加，原油中溶解的量也會增加，溶解的二氧化碳一方面增大了原油體積，儲存了彈性能量，驅(qū)替原油流動，另一方面降低了原油黏度和流動過程中的壓降損失、摩阻損失，使得自噴時間延長，停噴壓力降低。圖5-8停噴壓力曲線間距逐漸變小，因為出口壓力越低，生產(chǎn)壓差相對就變大，壓差增大了，就會使得自噴時間延長，停噴壓力值降低。

4 ?結(jié)論

（1）依據(jù)PR狀態(tài)方程而建立的用來預(yù)測停噴壓力的節(jié)點分析模型準(zhǔn)確性較強，通過模型分析的停噴壓力值對新疆油田A區(qū)實際生產(chǎn)具有較強的指導(dǎo)意義;

（2）在含水率相同的條件下，停噴壓力值隨二氧化碳含量的增加而降低，油井自噴生產(chǎn)周期變長;在二氧化碳含量相同的條件下，停噴壓力值隨含水率的上升而增大，含水率上升，油井自噴生產(chǎn)周期變短;

（3）綜合考慮二氧化碳含量和含水率因素，發(fā)現(xiàn)二氧化碳含量在對停噴壓力值預(yù)測時起主導(dǎo)作用，也就是說二氧化碳對自噴的有利作用大于含水率帶來的不利影響，表明了注氣生產(chǎn)相比于注水開發(fā)的獨特優(yōu)勢，也說明了二氧化碳驅(qū)油具有適用范圍廣，應(yīng)用前景大的優(yōu)點。

參考文獻(xiàn)：

[1]伍曉妮 .二氧化碳采油技術(shù)研究[D]. 東北石油大學(xué)， 2015.

[2]Welker J.. R， Dunlop D. D. Physical properties of carbonated oils [J].?Journal of Petroleum Technology， 1963， 15（8）：873-876.

[3]Lohrenz J.， Bray B. G.， Aime M.， Clark C. R. Calculating viscosity of reservoir ?fluids from their compositions[J]. Journal of Petroleum Technology， 1964， 16（10）：?1171-1176.

[4]Simon R.， Graue D. J. Generalized correlations for predicting solubility， swelling and viscosity behavior of CO₂-crude oil systems[J].Journal of Petroleum Technology， 1965， 17（1）：102-106.

[5]Mulliken C. A.， Sandler S. I.The prediction of CO₂?solubility and swelling factors for enhanced oil recovery[J]. Industrial?and Engineering Chemistry Process Design and Development， 1980，19（4）：709-711.

[6]Miller J. S.， Jones R. A. A laboratory study to determine physical characteristics of heavy oil after CO₂?saturation[C]. SPE 9789， 1981.

[7]Miller J. S.， Jones R. A. A laboratory study to determine physical characteristics of heavyoil after CO₂?saturation[R]. Washington D. C.： U. S. DOE， 1984.

[8]Mehrotra A. K.， Svrcek W. Y. Correlations for properties of bitumen saturated with CO₂， CH4 and N₂， and experiments with combustion gas mixtures[J].Journal of Canadian Petroleum Technology， 1982， 21 （6）：95-104.

[9]?Emanuel A. S. A mixing rule method for calculating the viscosity of CO₂?and low-gravity reservoir oil[C]. SPE 14229，1985.

[10]?Sankur V.， Creek J. L.， Di Julio S. S.， Emanuel A. S. A laboratory study of Wilmingtontar zone CO₂?injection project [J]. SPE Reservoir Engineering， 1986， 1 （1）：95-104.

[11]Chung F. T. H.， Jones R. A.， Burchfield T. E. Recovery of viscous oil under high pressure by CO₂?displacement： a laboratory study [C]. SPE 17588， 1988.

[12]Chung F. T. H.， Jones R. A.， Nguyen H. T. Measurements and correlations of the physical properties of CO₂/heavy-crude-oil mixtures [J]. SPE Reservoir Engineering， 1988， 3 （3）： 822-828.

[13]Quail B.， Hill G. A.， Jha K. N. Correlations of viscosity， density and gas ?solubility forSaskatchewan heavy oils[C]. PETSOC-SS-87-6，1987.

[14]?Quail B.， Hill G. A.， Jha K. N. Correlations of viscosity， gas solubility， and ?density for Saskatchewan heavy oils[J]. IndustrialEngineering Chemistry Research， 1988， 27 （3）?：519-523.

[15]?Dyer S. B.， Farouq Ali S. M. The potential of the immiscible carbon dioxide ?flooding process for the recovery of heavy oil[C]. PETSOC-SS-89-27，1989.

[16]Sayegh S.G.， Rao D. N.， Kokal S.， Najman J. Phase behaviour and physical properties of Lindbergh heavy oil/CO₂?mixtures [J].Journal of Canadian Petroleum Technology， 1990， 29 （6）： 31-39.

[17]De Ruiter R. A.， Nash L. J.， Singletary M. S. Solubility and displacement behavior of viscous crude with CO₂?and hydrocarbon gases [J]. SPE Reservoir Engineering， 1994，9（2）： 101-106.

[18]?Emera M. K.， Sarma H. K. A genetic algorithm-based model to predict CO₂-oil physicalproperties for dead and live oil [C]. PETSOC-2006-197，2006.

[19]Emera M. K.， Sarma H. K. A genetic algorithm-based model to predict CO₂-oil physical properties for dead and live oil [J].Journal of Canadian Petroleum Technology， 2008， 47（2）： 52-61.

[20]Al-Jarba M.， Al-Anazi B. D. A comparison study ofthe CO₂-oil physical properties literature correlations accuracy using visual basic modeling technique[J]. Nafta Scientific Journal， 2009， 60 （5）：287-291.

[21]?Mehrotra A. K.， Svrcek W. Y. Measurement and correlation of viscosity， density andgas solubility for Marguerite Lake bitumen saturated with carbon dioxide[J].AOSTARJournal of Research， 1984，1（1）：51-62.

[22]郝永卯，薄啟煒，陳月明.CO₂驅(qū)油實驗研究[J].石油勘探與開發(fā)，2005，32（2）：110-112.

[23]陶磊，王勇，李兆敏等.CO₂降粘劑改進(jìn)超稠油物性研究[J].陜西科技大學(xué)學(xué)報，2008，26（6）：25-29.

[24]周配配.特稠油CO₂超臨界降粘技術(shù)研究[D].東營：中國石油大學(xué)（華東），2010.

[25]董燕.CO₂超臨界稠油降粘輸送技術(shù)研究[D].東營：中國石油大學(xué)（華東），2011.

[26]Phillips P. The viscosity of carbon dioxide[J].Proc Roy Soc.， 1912 （87A）;48-51.

[27]Peng D Y and Robinson D B.A new two-constant equation of state[J].Ind.Eng.Chem.Fundam.1976（15）：59-63.

[28]Soave G.Equilibriumcomstants from a modified redlich-kwong equation of state[J].Chem.Eng.Sci.，1972，27（6）：1197-1203.

[29]?Starling K E， Han M S.Thermo data refined for LPG（Part 14：Mixture）[J]. Hydrocarbon Processing， 1972 51（5）;129.

[30]Hagedorn A R，Brown K E.Experimental study of pressure gradients occurring during continuous two phase flow in small diameter vertical conduits[J].J Pet Tech.，1965（3）：125-134.