邵陽 閆鐵 候兆凱

摘 ? ? ?要： CO2相變會導(dǎo)致天然氣黏度變化異常復(fù)雜，應(yīng)用實(shí)驗(yàn)手段測量含CO2天然氣黏度數(shù)據(jù)具有一定困難，因此需要在已有黏度實(shí)驗(yàn)數(shù)據(jù)的基礎(chǔ)上建立相關(guān)計算模型進(jìn)行預(yù)測。根據(jù)含CO2天然氣黏度的預(yù)測方法，對比計算氣體黏度的LBC、DS、LGE、Lucas等4種模型在預(yù)測CO2氣體黏度時的準(zhǔn)確性。通過計算比較發(fā)現(xiàn)，LGE黏度計算模型誤差較大，而DS方法能較準(zhǔn)確地確定含CO2天然氣的黏度。應(yīng)用DS模型對含CO2天然氣黏度進(jìn)行計算，計算結(jié)果表明：壓力升高時，氣體黏度就會增大;相同溫度、壓力條件下CO2含量越高，黏度越小。以上結(jié)論對研究超臨界CO2熱物理性質(zhì)、揭示井筒流動規(guī)律、油氣相變預(yù)測具有重要意義。

關(guān) ?鍵 ?詞：黏度模型;CO2氣體黏度;預(yù)測方法;影響因素

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

Abstract： CO2 phase change will cause extremely complicated change of natural gas viscosity. It is difficult to apply experimental methods to measure natural gas viscosity data. The accuracy of the existing CO2 natural gas viscosity calculation model is not high. Therefore， it is necessary to establish a relevant calculation model based on the existing experimental data for prediction. According to the prediction method of CO2-containing natural gas viscosity， the accuracies of four models including LBC， DS， LGE and Lucas for calculating the viscosity of gas were compared. Through calculation and comparison， it was found that the error of the LGE viscosity calculation model was large， and the DS method could determine the viscosity of the gas containing CO2 more accurately. Through analyzing the influencing factors， it was found that the gas viscosity increased when the pressure increased; the higher the CO2 content at the same temperature and pressure， the lower the viscosity. The above research is of great significance for revealing the wellbore flow law and the prediction of oil gas phase change.

Key words： Viscosity model; CO2 gas viscosity; Prediction method ; Influence factor

目前我國已發(fā)現(xiàn)的天然氣田中酸性氣田占比頗多，高含硫和CO2氣藏陸續(xù)被發(fā)現(xiàn)和開發(fā)[1-2]，研究酸性氣體的物理性質(zhì)將成為石油鉆探的熱門問題。

在酸性氣藏流體和超臨界CO2熱物理性質(zhì)中黏度均是流體的重要物理性質(zhì) [3]，其參數(shù)計算的準(zhǔn)確性為酸性氣田開發(fā)過程中的鉆井安全提供重要的理論基礎(chǔ)。目前含硫天然氣的黏度計算受到人們普遍重視[4]，而對含CO2氣體的黏度計算和預(yù)測研究相對不足。應(yīng)用實(shí)驗(yàn)手段測量CO2氣體黏度是最直接、可靠的方法[5]，由于在實(shí)際工程中，CO2氣體黏度隨溫度和壓力值變化范圍較大，難以滿足工程計算的需要，因此探尋快速、簡便、準(zhǔn)確的CO2氣體黏度計算模型十分必要。

酸性氣體中由于H2S和CO2 等非烴氣體組分的影響，氣體黏度通常較高[3，6-7]。在計算氣體黏度的主要預(yù)測方法中，經(jīng)驗(yàn)公式法一般有4種模型，即LBC模型、DS模型、LGE模型和Lucas模型[8-9]。這種方法是基于常規(guī)氣體黏度的經(jīng)驗(yàn)預(yù)測方法之上，通過擬合實(shí)驗(yàn)圖版，最后校正常規(guī)氣體黏度得到的 [3，6-7，10]，但此方法未計算實(shí)驗(yàn)值與實(shí)測值之間的誤差。因此，比較4種計算氣體黏度的預(yù)測模型，縮小計算氣體偏差因子與實(shí)驗(yàn)值之間的誤差，對有效指導(dǎo)井筒流動規(guī)律、油氣相變預(yù)測具有重要意義[11-13]。

1 ?黏度計算及其校正模型

1.1 ?經(jīng)驗(yàn)公式

目前，確定酸性天然氣黏度的方法主要是經(jīng)驗(yàn)公式法，即對已有的黏度的經(jīng)驗(yàn)公式進(jìn)行必要的校正，得到適用于含CO2或H2S組分的酸性天然氣黏度計算方法。

1.1.1 ?LBC模型

2 ?計算結(jié)果對比分析

利用實(shí)驗(yàn)手段測量CO2氣體黏度數(shù)據(jù)有一定困難，筆者在已有實(shí)驗(yàn)數(shù)據(jù)基礎(chǔ)上構(gòu)建模型來關(guān)聯(lián)和預(yù)測這些物性數(shù)據(jù)。本文數(shù)據(jù)來自文獻(xiàn)[21]提供的氣體組成及實(shí)驗(yàn)值，測試溫度80 ℃，壓力范圍5～44 MPa，文獻(xiàn)中流體組成與特征參數(shù)如表1所示。為了優(yōu)選出能夠適用于含CO2天然氣的黏度計算方法，分別采用以上經(jīng)驗(yàn)公式對上述樣品進(jìn)行計算，并以模型平均相對誤差為判斷標(biāo)準(zhǔn)對壓縮因子計算方法進(jìn)行優(yōu)選。

從圖1中可以看出，DS方法計算的平均相對誤差最小。因此，筆者在計算天然氣黏度時采用的是DS方法，并使用楊繼盛方法進(jìn)行校正。

圖2和圖3所示為含不同比例二氧化碳混合氣體的黏度隨溫度、壓力的變化曲線。隨壓力增加，混合氣體黏度逐漸增大。相同壓力、溫度條件下，CO2含量越高，黏度增大。

3 ?超臨界CO2氣體黏度性質(zhì)

在CO2氣體黏度的研究基礎(chǔ)上，討論超臨界CO2氣體黏度的物理性質(zhì)。CO2的溫度和壓力大于臨界點(diǎn)溫度和壓力時達(dá)到超臨界狀態(tài)，此時的CO2就是超臨界CO2。超臨界CO2的擴(kuò)散性和動力黏度均符合理想要求，是非常安全且潔凈的理想介質(zhì)[22-23]。

這里分析不同壓力條件下CO2黏度隨溫度的變化情況。由圖4可以看出，隨著溫度的增加，CO2黏度逐漸減小;在相同溫度條件下，壓力越大，氣體黏度也越大。從圖4中也可以看到，在臨界點(diǎn)附近，變化最為劇烈，也就是說在臨界點(diǎn)附近，CO2氣體黏度受溫度和壓力的影響最大。

通過對CO2黏度性質(zhì)的總結(jié)分析，處于超臨界區(qū)域內(nèi)的CO2的物理性質(zhì)受溫度和壓力的影響非常大。因此，這類氣體進(jìn)入井筒后不能按照一般流體進(jìn)行分析，避免造成井涌、井噴事故的發(fā)生，以免鉆井設(shè)備損壞、人員傷亡和環(huán)境污染。

4 ?結(jié) 論

1）基于經(jīng)驗(yàn)法的LGE模型和LUCAS模型雖然可以應(yīng)用在CO2氣體黏度的預(yù)測中，但誤差較大。

2） DS經(jīng)驗(yàn)公式法預(yù)測CO2氣體黏度的準(zhǔn)確性要高于其他的預(yù)測模型，若不滿足DS模型的計算條件，也可適當(dāng)考慮LBC計算模型。

3） 溫度對氣體壓縮系數(shù)影響不大，可以忽略。

4）利用DS模型對含CO2天然氣各種物性參數(shù)的影響因素進(jìn)行了分析，結(jié)果表明：CO2含量和溫度越高，偏差系數(shù)就越小;相同含量條件下，溫度越高，黏度越小;相同溫度、壓力條件下，CO2含量越高，黏度越小。

5）CO2黏度隨溫度增加而逐漸減小，相同溫度條件下，壓力越大，黏度也越大。處于超臨界區(qū)域內(nèi)的CO2物理性質(zhì)受溫度和壓力的影響非常大。此類氣體進(jìn)入井筒后不能按照一般流體進(jìn)行分析，避免安全事故的發(fā)生。

參考文獻(xiàn)：

[1]郭肖，王彭.含水對普光酸性氣田流體物性的影響[J].天然氣地球科學(xué)，2017，28（7）：1054-1058.

[2]王正東.各國由酸性氣體生產(chǎn)硫磺產(chǎn)量概況和酸性氣田開采中的一些問題[J].石油與天然氣化工，1980（4）：38-48.

[3]于忠.超臨界酸性天然氣密度黏度變化規(guī)律實(shí)驗(yàn)研究[D].青島：中國石油大學(xué)，2011.

[4]王彭.高含硫有水氣藏水侵動態(tài)分析[D].成都：西南石油大學(xué)，2018.

[5]楊曉鴻，張順喜，朱薇玲.天然氣黏度計算新模型[J].天然氣工業(yè)，2016，36（12）：113-118.

[6]楊學(xué)鋒，林永茂.酸性氣藏氣體黏度預(yù)測方法對比研究[J].特種氣藏，2005，12（5）：42-45.

[7]楊學(xué)鋒.高含硫氣藏特殊流體相態(tài)及硫沉積對氣藏儲層傷害研究[D].成都：西南石油大學(xué)，2006.

[8]歷勇.天然氣動力黏度的簡便計算方程[J].油氣田地面工程，2017，36（9）：9-11.

[9]郭肖，王彭.含水酸性氣體黏度計算方法[J]. 油氣藏評價與開發(fā)，2017，7（6）：46-51.

[10]楊繼盛.采氣工藝基礎(chǔ)[ M].北京：石油工業(yè)出版社，1992.

[11]馮金禹，田龍，劉嘉康，李世昌.氣液兩相流在臥式螺旋管內(nèi)摩阻特性實(shí)驗(yàn)研究與計算模型的建立[J]. 當(dāng)代化工，2019，48（7）：1452-1455.

[12]甘慶明，雷宇，伍振華，等.黏度對垂直管氣液兩相流壓降的影響[J].科學(xué)技術(shù)與工程，2019，19（19）：134-142.

[13]馮金禹，閆鐵，孫士慧，等.超臨界二氧化碳鉆井技術(shù)的研究進(jìn)展[J].現(xiàn)代化工，2018，38（6）：11-14.

[14]LOHRENZ J， BRAY B G， CLARK C R. Calculating viscosities of reservoir fluids from their compositions[J]. Journal of Petroleum Technology， 1964， 16（10）： 1171-1176.

[15]YANG X H， ZHU W L. A theoretical model for the density and temperature dependent viscosity of hydrocarbon gases[J]. Petroleum ?Science and Technology， 2016， 34（8）： 765-770.

[16]JARRAHIAN A， AGHEL B， HEIDARYAN E. ?On the viscosity of natural gas [J]. Fuel， 2015， 150（6）： 609-618.

[17]蘇云河，湯勇，肖云，等.CO2含量對火山巖氣藏開發(fā)指標(biāo)的影響[J].天然氣工業(yè)，2011，31（8）：69-72.

[18]白執(zhí)松，羅光熹.石油及天然氣物性預(yù)測[M].北京：石油工業(yè)出版社，1995.

[19] 楊繼盛.計算含H2S和CO2酸性天然氣高壓黏度的新方法[J].天然氣工業(yè)，1986，6（4）：108-110.

[20]吳晗，吳曉東，張慶生，等.普光氣田高含硫天然氣黏度計算模型優(yōu)選與評價[J].石油天然氣學(xué)報，2011，33（7）：157-160.

[21]卞小強(qiáng)，杜志敏.高含CO2天然氣相變及其物性參數(shù)實(shí)驗(yàn)測試[J].新疆石油地質(zhì)，2013（1）：63-65.

[22]AL-ADWANI F A， LANGLINAIS J， HUGHES R G. Modeling of an under- balanced-drilling operation using supercritical carbon dioxide[J]. SPE drilling & completion， 2009， 24（04）： 599-610.

[23]CABEZA L F， DE GRACIA A， FERN?NDEZ A I， et al. Supercritical CO2 as heat transfer fluid： A review[J]. Applied thermal engineering， 2017， 125： 799-810.

[24]朱自強(qiáng).超臨界流體技術(shù)：原理和應(yīng)用[M].北京：化學(xué)工業(yè)出版社，2000.

[25]ZHANG X， HEINONEN S， Lev?nen E. Applications of supercritical carbon dioxide in materials processing and synthesis[J]. Rsc Advances， 2014， 4（105）： 61137-61152.

[26]馮金禹，閆鐵，孫士慧，等.超臨界二氧化碳鉆井技術(shù)的研究進(jìn)展[J].現(xiàn)代化工，2018，38（6）：11-14.

[27]沈忠厚，王海柱，李根生.超臨界CO2連續(xù)油管鉆井可行性分析[J].石油勘探與開發(fā)，2010，37（6）：743-747.

[28]李昊.超臨界條件下井筒環(huán)空多相流動規(guī)律研究[D].青島：中國石油大學(xué)，2010.

[29]REDLICH O， KWONG J N S. An equation of state. Fugacities of gaseous solutions[J]. Chem. Rev， 1949， 44： 233-244.

[30]SOAVE G. Equilibrium constants from a modified Redlich-Kwong equation of state[J]. Chemical engineering science， 1972， 27（6）： 1197-1203.

[31]PENG D Y， ROBINSON D B. A new two-constant equation of state[J]. Industrial & Engineering Chemistry Fundamentals， 1976， 15（1）： 59-64.

[32]KNOTT R F， ANDERSON R N， ACRIVOS A. An experimental study of heat transfer to nitrogen-oil mixtures[J]. Industrial & Engineering Chemistry， 1959， 51（11）： 1369-1372.

[33]DRANCHUK P M， PURVIS R A， ROBINSON D B. Computer calculation of natural gas compressibility factors using the Standing and Katz correlation[C]. Annual Technical Meeting， Petroleum Society of Canada， 1973.

[34]KHOZE A N， DUNAEV S V， SPARIN V A. Heat and mass transfer in rising two-phase flows in rectangular channels[J]. Heat Transfer- Sov Res， 1976， 8（3） ： 154-161.

[35]TAITEL Y， BARNEA D. ?Counter current gas-liquid vertical flow， model for flow pattern and pressure drop[J]. International Journal of Multiphase Flow， 1983， 9（6）： 637-647.

[36]TENGESDAL J O， KAYA A S， SARICA C. Flow-pattern transition and hydrodynamic modeling of churn flow[J]. SPE Journal， 1999， 4（4）： 342-348.

[37]DITTUS F W， BOELTER L M K. Heat transfer in automobile radiators of the tubular type[J]. International Communications in Heat and Mass Transfer， 1985， 12（1）： 3-22.

[38]CHILTON T H， COLBURN A P. Mass transfer （absorption） coefficients prediction from data on heat transfer and fluid friction[J]. Industrial & engineering chemistry， 1934， 26（11）： 1183-1187.

[39]KISELEV S B， ELY J F. Generalized crossover description of the thermodynamic and transport properties in pure fluids[J]. Fluid Phase Equilibria ，2004， 252（1）：57-65.

[40]VIJAY M M， AGGOUR M A， SIMS G E. A correlation of mean heat- transfer coefficients for two-phase two-component flow in a vertical tube[C]. International Heat Transfer Conference Digital Library， Begel House Inc， 1982.