Vapor-liquid two-phase and vapor-liquid-solid three-phase behavior in three ternary hydrocarbon systems containing methane

1976 ◽  
Vol 21 (2) ◽  
pp. 220-222 ◽  
Author(s):  
W. F. O'Reilly ◽  
T. E. Blumer ◽  
K. D. Luks ◽  
J. P. Kohn
2020 ◽  
Vol 55 (1) ◽  
pp. 32-37
Author(s):  
A. Yu. Vorob’ev ◽  
V. A. Nebol’sin ◽  
N. Swaikat ◽  
V. A. Yuriev

SPE Journal ◽  
2011 ◽  
Vol 16 (04) ◽  
pp. 751-767 ◽  
Author(s):  
R.. Okuno ◽  
R.T.. T. Johns ◽  
K.. Sepehrnoori

Summary CO2 floods at temperatures typically below 120°F can involve complex phase behavior, where a third CO2-rich liquid (L2) phase coexists with the oleic (L1) and gaseous (V) phases. Results of slimtube measurements in the literature show that an oil displacement by CO2 can achieve high displacement efficiency of more than 90% when three hydrocarbon phases coexist during the displacement. However, the mechanism for the high-displacement efficiency is uncertain because the complex interaction of phase behavior with flow during the displacement is not fully understood. In this paper, we present the first detailed study of three-phase behavior predictions and displacement efficiency for low-temperature CO2 floods. Four-component EOS models are initially used to investigate systematically the effects of pressure, temperature, and oil properties on development of three-phase regions and displacement efficiency. Multicomponent oil displacements by CO2 are then considered. We use a compositional reservoir simulator capable of robust three-phase equilibrium calculations. Results show that high displacement efficiency of low-temperature CO2 floods is a consequence of both condensing and vaporizing behavior. The L2 phase serves as a buffer between the immiscible V and L1 phases within the three-phase region. Components in the L1 phase first transfer efficiently to the L2 phase near a lower critical endpoint (LCEP). These oil components then transfer to the V phase near an upper critical endpoint (UCEP) at the trailing edge of the three-phase region. The CEPs are defined where two of the three coexisting phases merge in the presence of the other immiscible phase. Unlike two-phase displacements, condensation and vaporization of intermediate components occur simultaneously within the three-phase region. The simultaneous condensing/vaporizing behavior involving the CEPs is also confirmed for simulations of several west Texas oil displacements. Quaternary fluid models can predict qualitatively the complex displacements because four is the minimum number of components to develop CEP behavior in composition space at a fixed temperature and pressure.


SPE Journal ◽  
2014 ◽  
Vol 19 (05) ◽  
pp. 943-955 ◽  
Author(s):  
Saeid Khorsandi ◽  
Kaveh Ahmadi ◽  
Russell T. Johns

Summary Minimum miscibility pressure (MMP) is one of the most important parameters in the design of a successful gasflooding process. The most-reliable methods to calculate the MMP are based on slimtube experiments, 1D slimtube simulations, mixing-cell calculations, and the analytical methods known as the method of characteristics (MOC). The calculation of MMP by use of MOC is the fastest method because it relies solely on finding the key tie lines in the displacement path. The MOC method for MMP estimation in its current form assumes that the composition path is a series of shocks from one key tie line to the next. For some oils, however, these key tie lines do not control miscibility, and the MMP calculated by use of the key-tie line approach can be significantly in error. The error can be as high as 5,000 psia for heavier oils or CO2 displacements at low temperature in which three-phase hydrocarbon regions can exist (L1–L2–V). At higher pressures, the two- or three-phase region can split (or bifurcate) into two separate two-phase regions (L1–L2 and L1–V regions). Thus, for the MMP calculation from MOC to be correct, we must calculate the entire composition path for this complex phase behavior, instead of relying on the shock assumption from one key tie line to the next. In this paper, the MOC-composition route is developed completely for the bifurcating phase-behavior displacement for pure CO2 injection by use of a simplified pseudoternary system that is analogous to the complex phase behavior observed for several real displacements with CO2. We develop the MOC analytical solutions by honoring all constraints required for a unique solution—velocity, mass balance, entropy, and solution continuity. The results show that a combination of shocks and rarefaction waves exists along the nontie-line path, unlike previous MOC solutions reported to date. We show that by considering the entire composition path, not just the key tie lines, the calculated MMP agrees with the mixing-cell method. We also show that, in this complex ternary displacement, the displacement mechanism has features of a both condensing and vaporizing (C/V) drive, which was thought to be possible only for gasfloods with four or more components. For pure CO2 injection, the solution also becomes discontinuous for oils that lie on the tie line envelope curve. Finally, we show that shock paths within the two-phase region are generally curved in composition space and that there is no MMP for some oil compositions considered in the displacements by CO2. Recovery can be large even though the MMP is not reached.


Author(s):  
Ruixue Li ◽  
Huazhou Andy Li

Multiphase isenthalpic flash calculations are often required in compositional simulations of steam-based enhanced oil recovery methods. These flash calculations are challenging in the narrow-boiling regions and in the determination of the correct number of existing phases. Based on the free-water assumption that the aqueous phase is pure water, a robust and efficient algorithm is proposed to perform isenthalpic three-phase flash calculations in this work. Multiphase equilibria can be considered by this algorithm, including single-phase equilibria, two-phase equilibria, and three-phase vapor-liquid-aqueous equilibria. Isenthalpic flash is a type of flash calculation conducted at given pressure and enthalpy for a feed mixture. In the proposed algorithm, assuming the feed is stable, the temperature is first determined by solving the energy conservation equation. Then the stability test on the feed mixture is conducted at the calculated temperature and the given pressure. If the mixture is found unstable, the two-phase and three-phase vapor-liquid-aqueous isenthalpic flash calculations can be simultaneously initiated without resorting to stability tests. To achieve simultaneous flashes, the outer loop is used to update the temperature by solving the energy conservation equation. The inner loop is used to obtain phase fractions and compositions by performing a three-phase free-water isothermal flash. Note that a two-phase isothermal flash will be initiated if an open feasible region in the phase fractions appears in any iteration during the three-phase isothermal flash or any of the ultimately calculated phase fractions from the three-phase flash do not belong to [0,1]. Negative flash is allowed in the three-phase free-water isothermal flash. A number of example calculations for water/hydrocarbon mixtures are carried out to test the robustness of the proposed algorithm. At low to medium pressures, a good agreement can be achieved between the results obtained by this algorithm and those obtained by the conventional algorithm. This algorithm performs well for the narrow-boiling regions, for example, the three-phase vapor-liquid-aqueous equilibrium region encountered for the water/hydrocarbon mixtures. During the iteration, the new algorithm can readily handle the appearance and disappearance of phases in the inner loop as temperature updates in the outer loop. The number of stability tests involved in the new algorithm is significantly reduced, helping to boost its computational efficiency.


Sign in / Sign up

Export Citation Format

Share Document