基于相稳定性分析的油气相包络线计算新方法

doi:10.3969/j.issn.2096-1693.2026.03.016

摘要/Abstract

摘要：

油气烃类体系相包络线稳定可靠的计算是油藏流体相态研究的重点和难点，Michelsen提出的饱和压力曲线追踪算法是该领域具有里程碑意义的工作，至今仍被广泛使用。但该方法技术细节复杂，正确实施难度较大，且计算中易出现无法收敛的情况导致相图无法生成。本文将相稳定性分析和二分法相结合，提出了一种能够精确绘制油气烃类体系相图的新方法。该方法通过相稳定性分析首先计算相包络线上的最大压力点和最大温度点，据此将压力-温度平面划分为4个区域，在每个区域分别使用稳定性分析结合二分法搜索相边界，最终绘制完整的相图。采用该方法对各类油藏流体的相图进行计算并于Michelsen方法进行对比，显示该算法能够可靠稳定的计算油气烃类体系的相图。

关键词: 相图, 相包络线, 稳定性分析, 饱和压力, 二分法

Abstract:

Stable and reliable calculation of phase envelopes for hydrocarbon systems is a key and challenging aspect of reservoir fluid phase behavior research. The saturation pressure curve tracing algorithm proposed by Michelsen represents a milestone contribution in this field and is still widely used today. However, this method involves complex technical details, making it difficult to implement correctly, and is prone to non-convergence issues during computation, which may prevent the generation of a complete phase diagram. This paper presents a novel method for accurately constructing phase diagrams of oil and gas hydrocarbon systems based on phase stability analysis combined with the bisection method. In this approach, phase stability analysis is first used to compute the cricondenbar and the cricondentherm on the phase envelope. Based on these two points, the pressure-temperature plane is divided into four regions. Within each region, the phase boundary is located using a combination of phase stability analysis and the bisection method, ultimately leading to the construction of a complete phase diagram. The proposed method is applied to compute phase diagrams of various reservoir fluids, and the results are compared with those obtained using Michelsen’s method. The comparison demonstrates that the proposed algorithm is capable of calculating phase diagrams of oil and gas hydrocarbon systems reliably and stably.

Key words: phase diagram, phase envelope, stability analysis, saturation pressure, bisection method

中图分类号:

赵海宁, 吴浩楠, 贺帅铭. 基于相稳定性分析的油气相包络线计算新方法[J]. 石油科学通报, 2026, 11(3): 980-992.

ZHAO Haining, WU Haonan, HE Shuaiming. A new method for oil and gas phase envelope construction based on phase stability analysis[J]. Petroleum Science Bulletin, 2026, 11(3): 980-992.

https://sykxtb.cup.edu.cn/CN/Y2026/V11/I3/980

图/表 11

图1 油气相稳定性分析求解流程

Fig. 1 Flow chart of hydrocarbon phase stability calculation procedure

图2 典型油气烃类体系相包络线特征及分区

Fig. 2 Phase region classification of typical hydrocarbon systems

图3 压力—温度平面划分

Fig. 3 Division of pressure-temperature plane

图4 基于相稳定性分析和二分法的相图构建算法流程

Fig. 4 Flowchart of the phase envelope construction algorithm based on phase stability analysis and the bisection method

图5 区域1中的相边界计算示意图

Fig. 5 Schematic diagram of phase boundary calculation in region 1

表1 相包络线算法计算结果概述

Table 1 Summary of phase envelope algorithm calculation results

序号	油气藏流体	组分数量	P_max/MPa	T_max/℃	本文算法与Michelsen算法相对运行时间*					计算结果
序号	油气藏流体	组分数量	P_max/MPa	T_max/℃	计算P_max	计算T_max	求解包络线	局部加密	总耗时	计算结果
1	凝析气1	15	23.509	126.65	31.84%	37.55%	28.75%	1.86%	2.35	图5a
2	凝析气2	9	24.229	300.58	26.92%	12.43%	60.27%	0.38%	3.24	图5b
4	挥发油1	15	32.261	368.76	8.21%	59.91%	31.68%	0.20%	6.23	图6a
5	挥发油2	15	31.743	341.16	7.30%	64.20%	28.17%	0.34%	6.49	图6b
6	挥发油3	40	29.227	447.62	4.63%	65.54%	29.71%	0.12%	2.19	图6c
7	挥发油4	15	31.184	491.63	31.39%	24.60%	43.85%	0.16%	3.63	图6d
8	黑油1	15	15.060	512.11	3.76%	67.76%	27.86%	0.62%	1.56	图7a
9	黑油2	15	16.323	490.07	7.67%	52.24%	39.72%	0.36%	1.89	图7b
10	黑油3	15	17.304	474.72	5.43%	58.48%	35.59%	0.50%	2.84	图7c
11	黑油4	28	14.484	453.17	5.24%	63.49%	30.51%	0.76%	2.00	图7d
12	黑油5	37	11.554	502.09	7.35%	44.91%	47.19%	0.56%	1.23	图7e
13	黑油6	24	17.932	527.55	3.67%	63.42%	32.45%	0.45%	2.37	图7f
14	稠油1	40	7.539	572.18	5.03%	63.85%	30.75%	0.37%	0.77	图8

图6 凝析气算例

Fig. 6 Gas condensate calculation examples

图7 挥发油算例

Fig. 7 Volatile oil calculation examples

图8 黑油算例

Fig. 8 Black oil calculation examples

图9 稠油算例

Fig. 9 Heavy oil calculation example

图10 本文算法和Michelsen算法就凝析气^[21]相图计算结果对比

Fig. 10 Comparison of phase diagram calculation results of gas condensate between the proposed algorithm and Michelsen’s algorithm^[21]

参考文献 49

[1]	刘华, 韩晓, 彭光荣, 等. 惠州凹陷惠州26洼古近系油气相态特征与成因机制[J]. 中国石油大学学报(自然科学版), 2025, 49(5): 31-43.
	[LIU H, HAN X, PENG G R, et al. Characteristics and genetic mechanisms of oil and gas phases in the Paleogene of Huizhou 26 sub-sag in Huizhou Depression[J]. Journal of China University of Petroleum (Edition of Natural Science), 2025, 49(5): 31-43.]
[2]	赵星星, 李斌, 邬光辉, 等. 塔里木盆地塔中Ⅲ区奥陶系多相态油气藏成因及富集模式[J]. 天然气地球科学, 2022, 33(1): 36-48.
	[ZHAO X X, LI B, WU G H, et al. Genesis and enrichment model of Ordovician multi-phase oil and gas reservoirs in TazhongⅢ block, Tarim Basin[J]. Natural Gas Geoscience, 2022, 33(1): 36-48.]
[3]	王晓锋, 刘鹏, 王杰, 等. 烃类分子相态分异过程的化学分馏机理与油气示踪意义[J]. 西北大学学报(自然科学版), 2025, 55(3): 601-612.
	[WANG X F, LIU P, WANG J, et al. Chemical fractionation during the phase separation of hydrocarbons and its applications in tracing petroleum origins[J]. Journal of Northwest University (Natural Science Edition), 2025, 55(3): 601-612.]
[4]	黄越义, 廖玉宏, 陈承声, 等. 塔里木盆地顺南1井和顺南4井油气相态演化的数值模拟与预测[J]. 石油与天然气地质, 2023, 44(1): 138-149.
	[HUANG Y Y, LIAO Y H, CHEN C S, et al. Numerical simulation and prediction of hydrocarbon phase evolution of wells Shunnan 1 and 4, Tarim Basin[J]. Oil & Gas Geology, 2023, 44(1): 138-149.]
[5]	陈承声, 邓瑞, 张海祖, 等. 塔里木盆地轮探1井下寒武统超深层油气相态演化定量模拟[J]. 天然气地球科学, 2023, 34(1): 96-110. doi: 10.11764/j.issn.1672-1926.2022.09.008
	[CHEN C S, DENG R, ZHANG H Z, et al. Quantitative simulation of phase evolution for ultra-deep oil and gas from Lower Cambrian strata of Well Luntan-1 in the Tarim Basin[J]. Natural Gas Geoscience, 2023, 34(1): 96-110.] doi: 10.11764/j.issn.1672-1926.2022.09.008
[6]	孙志道. 油气藏流体类型判别方法[J]. 石油勘探与开发, 1996(1): 69-75, 106.
	[SUN Z D. Methods for determining the type of different oil and gas reservoirs fluid[J]. Petroleum Exploration and Development, 1996(1): 69-75, 106.]
[7]	王金铎, 曾治平, 徐冰冰, 等. 准噶尔盆地沙湾凹陷二叠系上乌尔禾组流体相态及油气藏类型[J]. 岩性油气藏, 2024, 36(1): 23-31. doi: 10.12108/yxyqc.20240103
	[WANG J D, ZENG Z P, XU B B, et al. Fluid phase and hydrocarbon reservoir types of Permian Upper Urho Formation in Shawan Sag, Junggar Basin[J]. Lithologic Reservoirs, 2024, 36(1): 23-31.]
[8]	彭川, 邵才瑞, 葛祥, 等. 基于P—T相图的油气藏类型自动量化判识方法与应用[J/OL]. 天然气地球科学, 2026: 1-23. [2026-05-19].
	[PENG C, SHAO C R, GE X, et al. Automatic quantitative identification method and applications for oil and gas reservoir types based on P-T phase diagram[J/OL]. Natural Gas Geoscience, 2026: 1-23. [2026-05-19].]
[9]	汤勇, 陈焜, 胡小虎, 等. 页岩凝析气受限空间相态及开发特征研究[J]. 油气藏评价与开发, 2024, 14(3): 343-351.
	[TANG Y, CHEN K, HU X H, et al. Phase behavior and development characteristics of shale condensate gas in confined space[J]. Petroleum Reservoir Evaluation and Development, 2024, 14(3): 343-351.]
[10]	刘建仪, 杨雪, 刘勇. 低渗砂岩油藏CO₂驱相态及组分变化规律[J]. 特种油气藏, 2022, 29(6): 91-96. doi: 10.3969/j.issn.1006-6535.2022.06.011
	[LIU J Y, YANG X, LIU Y. Phase state and component change law of CO₂ flooding in low-permeability sandstone reservoirs[J]. Special Oil & Gas Reservoirs, 2022, 29(6): 91-96.]
[11]	宋传真, 赵海宁, 张慧. 氮气对塔河油田原油轻质组分抽提作用实验与评价[J]. 石油学报, 2021, 42(3): 341-349. doi: 10.7623/syxb202103007
	[SONG C Z, ZHAO H N, ZHANG H. Experiment and evaluation of extraction effect of nitrogen on light components of crude oil in Tahe oilfield[J]. Acta Petrolei Sinica, 2021, 42(3): 341-349.] doi: 10.7623/syxb202103007
[12]	戚志林, 严文德, 田杰, 等. 海上稠油注超临界多源多元热流体热采三维物理模拟实验[J]. 世界石油工业, 2026, 33(S1): 233-243.
	[QI Z L, YAN W D, TIAN J, et al. Three-dimensional physical simulation experiment of supercritical multi-source and multi-component thermal fluid injection for offshore heavy oil thermal recovery[J]. World Petroleum Industry, 2026, 33(S1): 233-243.]
[13]	汤勇, 龙科吉, 王皆明, 等. 凝析气藏型储气库多周期注采过程中流体相态变化特征[J]. 石油勘探与开发, 2021, 48(2): 337-346. doi: 10.11698/PED.2021.02.10
	[TANG Y, LONG K J, WANG J M, et al. Change of phase state during multi-cycle injection and production process of condensate gas reservoir based underground gas storage[J]. Petroleum Exploration and Development, 2021, 48(2): 337-346.] doi: 10.1016/S1876-3804(21)60027-7 URL
[14]	高德利, 姚军, 汪志明, 等. 简要论述力学与油气工程交叉研究[J]. 力学学报, 2025, 57(12): 2773-2788.
	[GAO D L, YAO J, WANG Z M, et al. A brief discussion on interdisciplinary research of mechanics and oil & gas engineering[J]. Chinese Journal of Theoretical and Applied Mechanics, 2025, 57(12): 2773-2788.]
[15]	张绪亮, 龚小卫, 李琰, 等. 超深井钻井井筒多组分流体超临界/气/液相态预测方法[J]. 科学技术与工程, 2025, 25(36): 15496-15505.
	[ZHANG X L, GONG X W, LI Y, et al. Prediction method for multi-component fluid supercritical gas-liquid phase state in ultra deep well drilling wellbore[J]. Science Technology and Engineering, 2025, 25(36): 15496-15505.]
[16]	史博会, 钱亚林, 王华青, 等. 管输天然气含水量/水露点的计算方法[J]. 油气储运, 2012, 31(3): 188-192.
	[SHI B H, QIAN Y L, WANG H Q, et al. Calculation method of water content/water dew point of natural gas[J]. Oil & Gas Storage and Transportation, 2012, 31(3): 188-192.]
[17]	任魁杰, 孟良, 赵海宁, 等. 黄桥高含CO₂气田产出液单井集输工艺[J]. 油气储运, 2020, 39(6): 650-655.
	[REN K J, MENG L, ZHAO H N, et al. Single well gathering and transportation process of produced fluid from high-CO₂ gasfield in Huangqiao[J]. Oil & Gas Storage and Transportation, 2020, 39(6): 650-655.]
[18]	吴玉国, 陈保东. BWRS方程在天然气物性计算中的应用[J]. 油气储运, 2003, 22(10): 16-21.
	[WU Y G, CHEN B D. The application of BWRS equation in calculating the thermo-physical properties of natural gas[J]. Oil & Gas Storage and Transportation, 2003, 22(10): 16-21.]
[19]	沈安琪, 张博源, 任思扬, 等. 基于人工神经网络优化算法的油气水三相相平衡计算[J/OL]. 地球科学, [2026-05-20].
	[SHEN A Q, ZHANG B Y, REN S Y, et al. A three-phase phase equilibrium calculation for oil, gas and water based on artificial neural network and optimization algorithm[J/OL]. Earth Science, 2026: 1-21. [2026-05-20].]
[20]	MICHELSEN M L. Calculation of phase envelopes and critical points for multicomponent mixtures[J]. Fluid Phase Equilibria, 1980, 4(1/2): 1-10. doi: 10.1016/0378-3812(80)80001-X URL
[21]	LI Y K, NGHIEM L X. The development of a general phase envelope construction algorithm for reservoir fluid studies[C]// SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 1982: SPE-11198-MS.
[22]	LINDELOFF N, MICHELSEN M L. Phase envelope calculations for hydrocarbon-water mixtures[J]. SPE Journal, 2003, 8(3): 298-303. doi: 10.2118/85971-PA URL
[23]	NIKOLAIDIS I K, ECONOMOU I G, BOULOUGOURIS G C, et al. Calculation of the phase envelope of multicomponent mixtures with the bead spring method[J]. AIChE Journal, 2016, 62(3): 868-879. doi: 10.1002/aic.v62.3 URL
[24]	NIKOLAIDIS I K, BOULOUGOURIS G C, PERISTERAS L D, et al. Efficient and robust methods for direct saturation point calculations[J]. Fluid Phase Equilibria, 2019, 500: 112242. doi: 10.1016/j.fluid.2019.112242 URL
[25]	NIKOLAIDIS I K, BOULOUGOURIS G C, PERISTERAS L D, et al. Construction of phase envelopes for binary and multicomponent mixtures with Euler-Newton predictor-corrector methods[J]. Fluid Phase Equilibria, 2020, 505: 112338. doi: 10.1016/j.fluid.2019.112338 URL
[26]	DEITERS U K, BELL I H. Calculation of phase envelopes of fluid mixtures through parametric marching[J]. AIChE Journal, 2019, 65(11): e16730. doi: 10.1002/aic.v65.11 URL
[27]	NICHITA D V. Phase envelope construction for mixtures with many components[J]. Energy & Fuels, 2008, 22(1): 488-495. doi: 10.1021/ef700251k URL
[28]	NICHITA D V. Density-based phase envelope construction[J]. Fluid Phase Equilibria, 2018, 478: 100-113. doi: 10.1016/j.fluid.2018.09.007 URL
[29]	XU L F, LI H Z. A trust-region-based phase envelope construction algorithm[J]. Geoenergy Science and Engineering, 2023, 229: 212058. doi: 10.1016/j.geoen.2023.212058 URL
[30]	DE AZEVEDO MEDEIROS F, STENBY E H, YAN W. Saturation point and phase envelope calculation for reactive systems based on the RAND formulation[J]. Chemical Engineering Science, 2022, 247: 116911. doi: 10.1016/j.ces.2021.116911 URL
[31]	MICHELSEN M L. The isothermal flash problem. Part I. Stability[J]. Fluid Phase Equilibria, 1982, 9(1): 1-19. doi: 10.1016/0378-3812(82)85001-2 URL
[32]	NICHITA D V. Phase stability testing near the stability test limit[J]. Fluid Phase Equilibria, 2016, 426: 25-36. doi: 10.1016/j.fluid.2016.01.015 URL
[33]	MICHELSEN M L. The isothermal flash problem. Part II. Phase-split calculation[J]. Fluid Phase Equilibria, 1982, 9(1): 21-40. doi: 10.1016/0378-3812(82)85002-4 URL
[34]	Baker L E, Pierce A C, Luks K D. Gibbs energy analysis of phase equilibria[C]// SPE/DOE Second Joint Symposium on Enhanced Oil Recovery, Tulsa, Oklahoma, USA, 1981.
[35]	LI Z D, FIROOZABADI A. General strategy for stability testing and phase-split calculation in two and three phases[J]. SPE Journal, 2012, 17(4): 1096-1107. doi: 10.2118/129844-PA URL
[36]	ZHAO H N, JING H B, FANG Z B, et al. Flash calculation using successive substitution accelerated by the general dominant eigenvalue method in reduced-variable space: Comparison and new insights[J]. SPE Journal, 2020, 25(6): 3332-3348. doi: 10.2118/202472-PA URL
[37]	PENG D Y, ROBINSON D B. A new two-constant equation of state[J]. Industrial & Engineering Chemistry Fundamentals, 1976, 15(1): 59-64.
[38]	WANG J, LIU H L, ZHAO H N, et al. Effects of the characterization methods of heptane plus components (C7+) on the phase behavior of volatile fluids[J]. Processes, 2024, 12(10): 2066. doi: 10.3390/pr12102066 URL
[39]	NEWLEY T M J, MERRILL R C Jr. Pseudocomponent selection for compositional simulation[J]. SPE Reservoir Engineering, 1991, 6(4): 490-496. doi: 10.2118/19638-PA URL
[40]	HEIDEMANN R A, KHALIL A M. The calculation of critical points[J]. AIChE Journal, 1980, 26(5): 769-779. doi: 10.1002/aic.v26:5 URL
[41]	LIU H L, XIE B, XIN X Z, et al. Investigation on the extent of retrograde condensation of Qianshao gas condensate reservoir using PVT experiments and compositional reservoir simulation[J]. Processes, 2024, 12(3): 503. doi: 10.3390/pr12030503 URL
[42]	KENYON D. Third SPE comparative solution project: Gas cycling of retrograde condensate reservoirs[J]. Journal of Petroleum Technology, 1987, 39(8): 981-997. doi: 10.2118/12278-PA URL
[43]	ZHAO H N, SONG C Z, ZHANG H, et al. Improved fluid characterization and phase behavior approaches for gas flooding and application on Tahe light crude oil system[J]. Journal of Petroleum Science and Engineering, 2022, 208: 109653. doi: 10.1016/j.petrol.2021.109653 URL
[44]	ZHANG J, XIA J, QIN J, et al. Insights on equation of state modeling PVT experiments for deep volatile oil reservoir[J]. Petroleum Science and Technology, 2024, 42(21): 3104-3125. doi: 10.1080/10916466.2023.2183967 URL
[45]	HEARN C L, WHITSON C H. Evaluating miscible and immiscible gas injection in the safah field, Oman[C]// SPE Reservoir Simulation Symposium, San Antonio, Texas, 1995: SPE-29115-MS.
[46]	HENDRIKS E M, VAN BERGEN A R D. Application of a reduction method to phase equilibria calculations[J]. Fluid Phase Equilibria, 1992, 74: 17-34. doi: 10.1016/0378-3812(92)85050-I URL
[47]	LANSANGAN R M, SMITH J L. Viscosity, density, and composition measurements of CO₂/west Texas oil systems[J]. SPE Reservoir Engineering, 1993, 8(3): 175-182. doi: 10.2118/21017-PA URL
[48]	JAUBERT J N, AVAULLEE L, SOUVAY J F. A crude oil data bank containing more than 5000 PVT and gas injection data[J]. Journal of Petroleum Science and Engineering, 2002, 34(1/2/3/4): 65-107. doi: 10.1016/S0920-4105(02)00153-5 URL
[49]	ZHAO H N, SONG C Z, WU D W, et al. Experimental Multiple-Contact miscibility and PVT properties with equation of state modeling for High-Pressure nitrogen injection into deep reservoirs with Tahe heavy and light crudes[J]. Fuel, 2025, 385: 134171. doi: 10.1016/j.fuel.2024.134171 URL

选择文件类型/文献管理软件名称

选择包含的内容