Gibbs自由能最小化法计算二氧化碳-烃-水系统相平衡

doi:10.7623/syxb201505008

石油学报 ›› 2015, Vol. 36 ›› Issue (5): 593-599.DOI: 10.7623/syxb201505008

Gibbs自由能最小化法计算二氧化碳-烃-水系统相平衡

邸元¹, 张园¹, Wu Yu-Shu²

1. 北京大学工学院北京 100871;
2. Petroleum Engineering Department, Colorado School of Mines Golden Colorado 80401

收稿日期:2014-11-06 修回日期:2015-02-17 出版日期:2015-05-25 发布日期:2015-06-10
通讯作者: 邸元,男,1968年2月生,1985年获西安交通大学学士学位,1995年获同济大学博士学位,现为北京大学工学院副教授,主要从事多孔介质多相流数值模拟、岩土力学等方面的研究工作。Email:diyuan@mech.pku.edu.cn
作者简介:邸元,男,1968年2月生,1985年获西安交通大学学士学位,1995年获同济大学博士学位,现为北京大学工学院副教授,主要从事多孔介质多相流数值模拟、岩土力学等方面的研究工作。Email:diyuan@mech.pku.edu.cn
基金资助:
国家重大科技专项(2011ZX05014-005)和国家重点基础研究发展计划(973)项目(2011CB201004)资助。

Phase equilibrium calculation of CO₂-hydrocarbons-water system based on Gibbs free energy minimization

Di Yuan¹, Zhang Yuan¹, Wu Yu-Shu²

1. College of Engineering, Peking University, Beijing 100871, China;
2. Petroleum Engineering Department, Colorado School of Mines, Golden Colorado 80401, USA

Received:2014-11-06 Revised:2015-02-17 Online:2015-05-25 Published:2015-06-10

摘要/Abstract

摘要：

油-气-水三相的相平衡计算是二氧化碳驱物性参数分析及数值模拟研究的核心问题之一。Gibbs自由能最小化法基于相平衡的热力学理论,适应性广、计算稳定,比较适合二氧化碳-烃-水系统的相平衡计算。通过引入相稳定性变量,采用拉格朗日乘数法计算多相多组分系统Gibbs自由能的最小化问题,由Newton-Raphson迭代法求解相应的控制方程。迭代计算时,二氧化碳和烃类组分气-油平衡常数的初值采用Wilson公式进行计算。气相和油相的状态方程采用SRK方程。分别计算了凝析气、酸气和二氧化碳-烃混合物3个算例的相平衡问题,计算值同实验实测值符合良好,验证了算法的准确性。分别对二氧化碳-烃系统和二氧化碳-烃-水系统算例的相平衡进行了计算,研究了一定温度和压力条件下,系统内二氧化碳和水的含量对烃类组分平衡常数的影响。

关键词: 相平衡, Gibbs自由能最小化, 平衡常数, 多组分混合物, 二氧化碳

Abstract:

Phase equilibrium calculation of oil-gas-water system is one of the most essential problems for properties analysis and numerical simulation of carbon dioxide flooding. Based on thermodynamics of phase equilibrium, Gibbs free energy minimization is an adaptable and stable method, and is very suitable for three-phase equilibrium calculations of carbon dioxide-hydrocarbons-water system. A phase stability variable was defined, and the Lagrange multiplier method was used for Gibbs free energy minimization problem of multiphase and multicomponent systems. The corresponding governing equations were solved by Newton-Raphson method. The initial gas-oil equilibrium constants of carbon dioxide and hydrocarbon components were obtained by Wilson equation. The SRK equation of state was adopted for gas and oil phase. Several multicomponent examples were calculated, including condensate gas, sour gas and mixtures of carbon dioxide-hydrocarbons. The simulated results agree reasonably well with the experimental data reported in the literature. The examples of carbon dioxide-hydrocarbon system and carbon dioxide-hydrocarbons-water system were also calculated. The effects of contents of carbon dioxide and water on the equilibrium constants of hydrocarbons under certain conditions of temperature and pressure were presented.

Key words: phase equilibrium, Gibbs free energy minimization, equilibrium constant, multi-component mixture, carbon dioxide

中图分类号:

O642.4

邸元, 张园, Wu Yu-Shu. Gibbs自由能最小化法计算二氧化碳-烃-水系统相平衡[J]. 石油学报, 2015, 36(5): 593-599.

Di Yuan, Zhang Yuan, Wu Yu-Shu. Phase equilibrium calculation of CO₂-hydrocarbons-water system based on Gibbs free energy minimization[J]. Acta Petrolei Sinica, 2015, 36(5): 593-599.

参考文献

[1] Brock W R, Bryan L A.Summary results of CO₂ EOR field tests, 1972-1987[R].SPE 18977, 1989.
[2] Reid B G, David S S.State of the industry in CO₂ floods[R].SPE 38849, 1997.
[3] Lambert M R, Marino S D.Implementing CO₂ floods:no more delays[R].SPE 35187, 1996.
[4] 邸元, 吴玉树, 鞠斌山.裂缝性油藏二氧化碳驱的多组分数值计算模型[J].西南石油大学学报:自然科学版, 2015, 37(2):93-100.
Di Yuan, Wu Yushu, Ju Binshan.Model formulation for compositional numerical simulation of CO₂ flooding in fractured reservoirs[J].Journal of Southwest Petroleum University:Science & Technology Edition, 2015, 37(2):93-100.
[5] Sanderson R V, Chien H H Y.Simultaneous chemical and phase equilibrium calculation[J].Industrial and Engineering Chemistry Process Design and Development, 1973, 12(1):81-85.
[6] Almehaideb R A, Abdulkarim M A, Al-Khanbashi A S.Improved K-value correlation for UAE crude oil components at low pressure using PVT laboratory data[J], Fuel, 2001, 80(1):117-124.
[7] Michelsen M L.The isothermal flash problem.Part I.Stability[J].Fluid Phase Equilibria, 1982, 9(1):1-19.
[8] Michelsen M L.The isothermal flash problem.Part II.Phase-split calculation[J].Fluid Phase Equilibria, 1982, 9(1):21-40.
[9] Haugen K B, Sun L, Firoozabadi A.Three-phase equilibrium calculations for compositional simulation[R].SPE 106045, 2007.
[10] Pan H, Firoozabadi A.Complex multiphase equilibrium calculations by direct minimization of Gibbs free energy by use of simulated annealing[J].SPE Reservoir Evaluation and Engineering, 1998, 1(1):36-42.
[11] Nichita D V, Gomez S, Luna E.Multiphase equilibria calculation by direct minimization of Gibbs free energy with a global optimization method[J].Computers & Chemical Engineering, 2002, 26(12):1703-1724.
[12] Ballard A L, Sloan Jr.E D.The next generation of hydrate prediction:Part III.Gibbs energy minimization formalism[J].Fluid Phase Equilibria, 2004, 218(1):15-31.
[13] Yarborough L.Vapor-liquid equilibrium data for multicomponent mixtures containing hydrocarbon and nonhydrocarbon components[J].Journal of Chemical & Engineering Data, 1972, 17(2):129-133.
[14] Alvarado, F J.Studies on phase equilibrium calculations from equations of state[D].College Station:Texas A & M University, 1993.
[15] 郑德馨, 刘芙蓉.多组分气体分离[M].西安:西安交通大学出版社, 1988, 6-11.
Zheng Dexin, Liu Furong.Multi-component gas separation[M].Xi'an:Xi'an Jiaotong University Press, 1988, 6-11.
[16] 郭天民.多元汽-液平衡和精馏[M].2版.北京:石油工业出版社, 2002:11-14.
Guo Tianmin.Multiple vapor-liquid equilibrium and distillation[M].2nd ed.Beijing:Petroleum Industry Press, 2002:11-14.
[17] Gupta A K.Steady state simulation of chemical process[D].Calgary:University of Calgary, 1988.
[18] Wilson G M.A modified Redlich-Kwong equation of state, application to general physical data calculations:proceedings of 65th National AIChE Meeting, Cleveland, USA, 1969[C].AIChE:Cleveland, 1969.
[19] Rahman S.Experimental K-value prediction[R].SPE 125413, 2009.
[20] Al-Jawad M S, Hassan O F.Comprehensive model for flash calculations of heavy oils using the Soave-Redlich-Kwong equation of state[R].SPE 150083, 2012.
[21] Soave G.Equilibrium constants from a modified Redlich-Kwong equation of state[J].Chemical Engineering Science, 1972, 27(6):1197-1203.
[22] 刘光启, 马连湘, 刘杰.化学化工物性数据手册[M].北京:化学工业出版社, 2002.
Liu Guangqi, Ma Lianxiang, Liu Jie.Chemical property data manual[M].Beijing:Chemical Industry Press, 2002.
[23] Habiballah W A, Startzman R A, Barrufet M A.Use of neural networks for prediction of vapor/liquid equilibrium K-values for light-hydrocarbon mixtures[J].SPE Reservoir Engineering, 1996, 11(2):121-126.
[24] Mehra R K, Heidemann R A, Aziz K.Computation of multiphase equilibrium for compositional simulation[J].Society of Petroleum Engineers Journal, 1982, 22(1):61-68.
[25] Turek E A, Metcalfs R S, Yarborough L, et al.Phase equilibria in CO₂-multicomponent hydrocarbon systems:experimental data and an improved prediction technique[J].Society of Petroleum Engineers Journal, 1984, 24(3):308-324.
[26] Metcalfe R S, Yarborough L.The effect of phase equilibria on the CO₂ displacement mechanism[J].Society of Petroleum Engineers Journal, 1979, 19(4):242-252.

Gibbs自由能最小化法计算二氧化碳-烃-水系统相平衡

Phase equilibrium calculation of CO₂-hydrocarbons-water system based on Gibbs free energy minimization

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

[1]	史晓东, 孙灵辉, 战剑飞, 李斌会, 韩雪, 陆慧敏, 蒋成刚. 松辽盆地北部致密油水平井二氧化碳吞吐技术及其应用[J]. 石油学报, 2022, 43(7): 998-1006.
[2]	孙小辉, 孙宝江, 王志远, 李航宇, Pan Lehua, 高永海, 廖友强. 酸性天然气与水体系内超临界-气-液-水合物多相平衡模型[J]. 石油学报, 2021, 42(9): 1212-1223.
[3]	宋传真, 赵海宁, 张慧. 氮气对塔河油田原油轻质组分抽提作用实验与评价[J]. 石油学报, 2021, 42(3): 341-349.
[4]	韩海水, 李实, 姚小琪, 陈兴隆, 张可, 秦积舜. 基于摩尔密度的原油-CO₂体系膨胀能力预测方法[J]. 石油学报, 2018, 39(4): 456-462.
[5]	赵跃军, 宋考平, 范广娟, 皮彦夫, 刘丽. 酯类化合物降低原油与二氧化碳体系最小混相压力实验[J]. 石油学报, 2017, 38(9): 1066-1072.
[6]	贺振国, 李根生, 王海柱, 田守嶒, 沈忠厚, 邵尚奇, 刘新云. 超临界二氧化碳射流中的颗粒跟随及影响因素[J]. 石油学报, 2017, 38(4): 475-484.
[7]	侯磊, 孙宝江, 蒋廷学, Geng Xueyu, 王志远, 李庆超. 支撑剂在超临界二氧化碳中的跟随性计算[J]. 石油学报, 2016, 37(8): 1061-1068.
[8]	韩海水, 李实, 陈兴隆, 秦积舜, 曾保全. CO₂对原油烃组分膨胀效应的主控因素[J]. 石油学报, 2016, 37(3): 392-398.
[9]	赵青, 李玉星, 李顺丽. 超临界二氧化碳管道杂质对节流温降的影响[J]. 石油学报, 2016, 37(1): 111-116.
[10]	郭建春, 曾冀, 张然, 周长林. 井筒注二氧化碳双重非稳态耦合模型[J]. 石油学报, 2015, 36(8): 976-982.
[11]	李元生, 李相方, 藤赛男, 徐大融. 考虑孔喉与流体分子相互作用及毛管力的相平衡计算[J]. 石油学报, 2015, 36(4): 511-515.
[12]	郭建春, 曾冀. 超临界二氧化碳压裂井筒非稳态温度-压力耦合模型[J]. 石油学报, 2015, 36(2): 203-209.
[13]	孙小辉, 孙宝江, 王志远. 超临界CO₂压裂裂缝温度场模型[J]. 石油学报, 2015, 36(12): 1586-1592.
[14]	张运军, 沈德煌, 高永荣, 李军辉, 李杰, 刘岩. 二氧化碳气体辅助SAGD物理模拟实验[J]. 石油学报, 2014, 35(6): 1147-1152.
[15]	程宇雄, 李根生, 沈忠厚, 王海柱, 李敬彬. 超临界二氧化碳射流冲击压力参数影响规律[J]. 石油学报, 2014, 35(4): 765-770.

Gibbs自由能最小化法计算二氧化碳-烃-水系统相平衡

Phase equilibrium calculation of CO2-hydrocarbons-water system based on Gibbs free energy minimization

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics

本文评价

Phase equilibrium calculation of CO₂-hydrocarbons-water system based on Gibbs free energy minimization