Three-Phase Equilibrium Calculations of Water/Hydrocarbon/Nonhydrocarbon Systems Based on the Equation of State (EOS) in Thermal Processes

A simple and novel approach is proposed to represent the mutual solubility of water and hydrocarbon components based on equations of state at high temperatures in thermal recovery processes. Sϕreide and Whitson modifications are applied to the Peng–Robinson (PR) equation of state (EOS) so that all components, including the water component, can exist in all phases, reasonably representing gas solubility in water and water solubility in hydrocarbon phases. We propose an algorithm to assign binary interaction parameters (BIPs) for aqueous and nonaqueous phases. The water vapor pressure helps select initial K-values for stability analysis so that the aqueous phase can be split out first if present. The algorithm is tested by a wide range of variations in pressure, temperature, and composition. The results show the robustness of the algorithm and the effects of temperature and overall water mole fraction on phase behaviors in steam flooding processes.


■ INTRODUCTION
Water is present either as initial water or as injection water and hence in contact with hydrocarbon phases in equilibrium in reservoirs. During the production of heavy oils often involved with three phases including oil, gas, and water, the K-value method is adopted in commercial simulators (Eclipse Technical Description 2013; STARS 2015; Intersect 2018) to predict phase behaviors in thermal processes (steam flooding). However, this method is not thermodynamically consistent because K-values are not of thermodynamic characteristics. 1 Due to this drawback, the K-value method is competent for the performance of production where the Kvalues are weak functions of composition. 2,3 However, especially for more complicated EOR methods such as Expanding Solvent Steam-Assisted Gravity Drainage (ES− SAGD), the K-value method is not accurate enough to represent the thermodynamic properties of three-phase fluid mixtures since it is only a function of temperature and pressure. In addition, the K-value method does not consider the mutual solubility of water and hydrocarbon, which can affect calculations of viscosity, density, and other thermodynamic properties. However, as the temperature increases, the solubility of water in hydrocarbon phases and the solubility of hydrocarbon in the aqueous phase are significantly increased. 4,5 Several authors have investigated water solubility in hydrocarbon phases at high temperatures. Griswold and Kasch 6 presented the concept of water/oil solubility at high temperatures. They indicated that water solubility in 54.3°API naphtha is 16.18 mol % at 431.67°F, water solubility in 42°API kerosene is 34.97 mol % at 507.2°F, and water solubility in 29.3°API oil is 43.44 mol % at 537.8°F. Glandt and Chapman 7 showed that water solubility in oil is up to 33.3 wt %. 8−18 In commercial simulators, Henry's law is used to calculate the solubility in the aqueous phase. 19−21 The cubic equation of state (EOS), meanwhile, is used to represent the fugacity coefficients in the nonaqueous phases. Two different formulations constitute an inconsistent algorithm for phase equilibrium calculations. Also, statistical associating fluid theory (SAFT) has also been proposed to represent the phase behavior of mixtures containing the aqueous phase. But, it requires knowledge of corresponding molecular structures. Thus, it is not suitable for industrial applications. Huron and Vidal 22 presented a nonsymmetrical model on the basis of the SRK equation, but the accuracy of the model is not high enough. 23 Chapman 24 developed a sum of three Helmholtz energy terms based on the EOS called SAFT (statistic associating fluid theory). Kontogeorgis 25,26 proposed the association of compressibility factor based on cubic models called cubic plus association (CPA). Li and Firoozabadi 27 described the cross association between non-water components and water in the framework of perturbation theory. Li 28 noted that SAFT and CPA models usually exhibit good performance for standalone calculations, but they are not timeefficient when used for numerical simulations.
Another approach for equilibrium calculations for thermal compositional simulation is based on equations of state. Phase behavior representation of water-containing mixtures using the equation of state remains a challenge due to the nonideality brought by strong hydrogen bonds among water molecules. A good reproductivity of compositions in nonaqueous phases can be achieved, but the calculated composition for the aqueous phase is inconsistent with the experimental data by several orders of magnitude since slight amounts of multiple components usually exist in the water phase.
An accurate prediction of the phase behavior of mixtures using the equation of state depends on appropriately chosen binary interaction parameters, mixing rules, and the α function for pure substances. The most important modification of the original EOS for water/hydrocarbon systems is the introduction of two sets of binary interaction parameters for constant a. Whitson and Brule 29 verify the accuracy of simultaneous application of aqueous-and nonaqueous-phase interaction coefficients for mutual-solubility predictions of binaries and natural gas/water/brine mixtures, suggesting that the modification is warranted. Li and Yang 30 developed a new binary interaction parameter correlation to determine mutual solubility between CO 2 and water. The polynomial temperature-dependent BIP correlation in the aqueous phase predicts more accurately than the exponential BIP correlation. Composition-and density-dependent mixing rules have also been proposed for modifying cubic EOSs for water/hydrocarbon systems. Panagiotopoulos and Reid's 31 linear composition-dependent mixing rule has received considerable interest. Unfortunately, it violates several fundamental thermodynamic conditions. To improve the predictions of water vapor pressure and water solubility in the nonaqueous phases, Peng and Robinson proposed a modified α function for water, which can be used for 284.8 < T < 466.1 K. Sϕreide and Whitson 3232 replaced the α function with a new correlation. It can predict water vapor pressures for 288.2 < T < 598.15 K. Li and Yang 30 developed a modified α function that is able to more accurately reproduce the water vapor pressure in the full temperature range of 273.15−647.10 K.
Because of a more extensive temperature range, we adopt the α function from Li and Yang. 30 Due to the simplicity and reasonable prediction of mutual solubility between light hydrocarbons and water, we adopt BIP correlations from the Sϕreide and Whitson 32 model to match the mutual solubility of hydrocarbon and water. Mohebbinia et al. 33,34 also adopted Sϕreide and Whitson modifications to represent gas and water mutual solubility. However, the method of phase identification has not been mentioned, which is required in advance to determine the different sets of BIPs between hydrocarbon components and water in the aqueous phase and nonaqueous phases. We predict whether an aqueous phase is present by virtue of the water vapor pressure curve. If an aqueous phase may be present, the equilibrium calculation is performed first with a nearly pure water phase as a trial phase for the first stability analysis to obtain phase fractions and compositions of the hydrocarbon-rich phase and the aqueous phase.
The sequential algorithm combining phase stability analysis and phase splitting is used. 35−38 First, according to the overall composition, the stationary point method is used at the present temperature and pressure to determine how many phases exist. If the mixture is unstable, another phase is added and an ensuing flash calculation provides phase fractions and phase compositions. The heavier phase with a heavier molecular weight is applied to the second stability analysis. 39 If the heavier phase is stable, then two phases are in equilibrium; if it is not stable, another phase is added and a three-phase flash calculation is required. 40 The paper is structured along the following lines. First, we present the binary interaction parameter correlations from Sϕreide and Whitson 32 and the α function for water from Li and Yang, 30 followed by the characteristic of the water vapor pressure curve. Then, the selection of initial guesses for stability analysis is presented, followed by the gas−oil−water three-phase equilibrium algorithm. The results of our method compared with other models to represent water/hydrocarbon mutual solubility are presented. Several examples taken from the literature are tested, including the construction of phase envelopes in two-and three-phase regions. And this section also illustrates variations in phase distribution as an overall water mole fraction and temperature increase. Finally, the main results and conclusions are summarized.
Modified PR EOS for the Aqueous Phase. Phase properties are calculated by PR EOS in our work, and all components, including water, can exist in all phases. Water is not a separate phase. Li and Yang 30 developed a modified αterm to more accurately reproduce the water vapor pressure in the full temperature range of 273.15−647.10 K.
where T r is the reduced temperature for water. According to Sϕreide and Whitson, 32 two sets of binary interaction parameters (BIPs) between hydrocarbon components and water are proposed for nonaqueous phases and the aqueous phase. As a result, two different attraction terms in the EOS are calculated as a function of their respective BIPs.
Sϕreide and Whitson 32 suggested that a constant BIP for the nonaqueous phase and a temperature-dependent BIP for the aqueous phase were found adequate to match experimental mutual solubility data with reasonable accuracy. BIP constants for the nonaqueous phase are listed in Table 1, and a constant equal to 0.5 is applied for hydrocarbon/water binaries not listed in Table 1. Sϕreide and Whitson 32 proposed four correlations to calculate BIPs in the aqueous phase between hydrocarbon/ brine, CO 2 /brine, N 2 /brine, and H 2 S/brine. Recently, Li and Yang 30 proposed the BIP correlation for CO 2 /brine. This work adopts the BIP correlation for CO 2 /brine from Li and Yang, 30 while BIP correlations for hydrocarbon/brine, N 2 /brine, and H 2 S/brine are taken from Sϕreide and Whitson. 32 where T ri is the reduced temperature, c sw is the salinity of brine, and ω i is the acentric factor for each component. Compared with the original PR EOS formulation, BIP and α function modifications to PR EOS provide a more accurate representation of the true multiphase equilibrium than the original PR EOS approach. This method entails the identification of phase types upfront, so that different sets of BIPs can be applied to nonaqueous and aqueous phases.
In steam flooding processes, overall mixture compositions in reservoirs are gradually dominated by the injection steam, i.e., the water component. For stability analysis and phase-splitting calculations involving the aqueous phase, if the water component accounts for no less than 50%, the mixture is taken as waterlike, and binary interaction parameters (BIPs) for the aqueous phase are used for the evaluation of component fugacity coefficients. If the mole fraction of the water component is less than 50%, binary interaction parameters (BIPs) for the nonaqueous phase are chosen for the phase.
Vapor Pressure. In thermal settings, a key point is that the saturated steam curve can be used to guide the selection of the trial phase for stability analysis. This distinctive characteristic for hydrocarbon−water systems is sensitive to different phase states of the systems. If the reservoir pressure is lower than the water vapor pressure at its current temperature, the aqueous phase is unlikely to exist. Then, stability analysis is started with Wilson's correlation and its inverse as initial guesses. If the reservoir pressure is higher than the water vapor pressure, the aqueous phase can be identified by K 2 HO as an initial guess for stability analysis. The rationale behind the method is that the aqueous phase may only exist at pressures higher than water vapor pressure. If the first stability analysis result shows that the aqueous phase is not present, Wilson's correlation and its inverse are used as initial guesses again for the stability analysis.
We use eq 8 to calculate water vapor pressure 41 , T c is the critical temperature, and P c is the critical pressure.
Note that the validation of the above correlation corresponds to the temperature range between 273.16 and 647.14 K. As Whitson and Brulé2 9 point out, the critical temperature defines the temperature above which any gas/ liquid mixture cannot coexist, regardless of pressure. Similarly, the critical pressure defines the pressure above which liquid and vapor cannot coexist, regardless of the temperature. At the critical point, the vapor and liquid phases can no longer be distinguished, and their intensive properties are identical. So, the overall mixture contains no aqueous phase when the temperature is higher than 647.14 K. 42−45

■ INITIAL GUESS FOR STABILITY ANALYSIS
In the stability analysis, the local minimum of the tangent plane distance function (TPD) is strongly dependent on the initial guess of the trial phase composition, which can directly affect equilibrium ratios K converged at the end of stability analysis. Improper initial guesses may miss some stationary points and fail to detect phase instability. Furthermore, stability analysis provides initial K-values for phase splitting calculations. Inappropriate initial K-values may result in the failure of flash calculations. For vapor−liquid equilibrium systems, {K Wilson } can provide good initial guesses for stability analysis. 46 However, for liquid−liquid equilibrium systems, {K Wilson } becomes unreliable and may lead to failures in detecting instability. Our work adopts the stationary point location method from Michelsen 35 with an initial guess suggested by Li and Firoozabadi 47 for the second stability analysis to overcome the setback in equilibrium calculations of more than one liquid phase. 48−51 We use two types of multiple initial K-value estimates for the first stability analysis. If the aqueous phase is not present, the vapor−liquid mixtures comprised of hydrocarbon components are tested by Wilson's correlation eq 9 and its inverse as initial guesses. Otherwise, stability analysis is initiated by a nearly pure water phase as the trial phase when the aqueous phase exists in equilibrium.
where T ci , p ci , and ω i are the critical temperature, critical pressure, and acentric factor of component i. In a nearly pure aqueous phase, the fraction of water component is 99 mol %.
And the other (N c − 1) components equally share the remaining 1 mol % in the trial phase. When the aqueous phase is present in the feed composition, the first stability analysis is started with K pure (water accounts for 99% in the trial phase). If the mixture is unstable, the converged nontrivial K-values from the first stability analysis are used as the initial K-values for the two-phase splitting calculation, so that the aqueous phase can be split out of the mixture first. Also, only after the aqueous phase is distinguished from nonaqueous phases can we distribute two types of BIPs for the aqueous phase and nonaqueous phases.
In the three-phase flash calculations, the second stability analysis is performed to check whether the physical solution of the first flash is stable in the two-phase state. According to Li and Firoozabadi, 47 the heavier phase (with higher molecular weight) is selected as the test phase because the aqueous phase is mainly composed of water, and accordingly, the aqueous phase is not the heavier phase. The hydrocarbon phase is naturally chosen to be tested as the test phase. To overcome the intrinsic setback of the stationary point method, we use multiple initial estimates in the second stability analysis. Li and Firoozabadi 47 suggested that the initial composition of one component is 90% and that the other components equally share the remaining 10 mol % of the trial phase.
pure test (13) There are N c sets of initial guesses of {K j pure } because every component has the opportunity to be assigned 90 mol %.
For the second stability analysis, the initial guesses in our work are K i 1stab represents the converged nontrivial solution of the first stability analysis corresponding to the lowest TPD and K i 1flash is the K-values from the two-phase splitting calculation. The number of initial guesses in the second stability analysis is N c + 8. All of these estimates may detect phase instability.
■ ALGORITHM FOR GAS−OIL−WATER PHASE EQUILIBRIUM CALCULATION Stability analysis turns out to be a heavy burden in the whole process of multiphase equilibrium calculations because of the existence of multiple initial estimates for trial phases. We propose an efficient and reliable scheme for phase behavior calculations, taking advantage of characteristic behaviors of water/hydrocarbon mixtures.
If the current pressure is lower than the water vapor pressure, water is not present. The stability analysis is initiated with Wilson's correlation and its inverse as trial phases. The result may be single-phase, either as a liquid or as vapor, or a two-phase mixture that requires a subsequent two-phase splitting calculation. If the reservoir pressure is higher than the water vapor pressure, the aqueous phase may exist in equilibrium. Then, the stability analysis is started with a nearly pure water phase as the initial estimate. If TPD > 0, the aqueous phase does not exist. Then, the same strategy is adopted as the previous case when pressure is below the steam pressure. If the stability test result shows that the mixture is not stable, the converged nontrivial K-values from the nearly pure water trial phase are used to initiate the two-phase splitting calculation so that the aqueous phase is split out as the second phase. Thus, the identity of each phase in the splitting process is specified and different BIPs for the aqueous phase and the hydrocarbon phase are assigned for each phase. The second stability analysis is applied to the heavier phase, i.e., the hydrocarbon phase. The trial phases are related to Wilson's correlation, its inverse, its cubic, the inverse of cubic, and nearly pure phases consisting of certain components except for water. If the second stability analysis shows that the heavier phase is unstable, then the K-values corresponding to the lowest TPD in the second stability analysis are selected to initiate the three-phase splitting calculation. BIPs for the aqueous phase are assigned for the third phase, and BIPs for nonaqueous phases are assigned for the first two phases. If the result of the second stability analysis is stable, the system is a two-phase mixture. The main steps involved in the gas/oil/ water three-phase equilibrium calculations are as follows: • Step 1: Evaluate the water vapor pressure at the current temperature by eq 8 to determine whether an aqueous phase is present. If an aqueous phase is not present, go to step 3; if an aqueous phase may be present, go to step 2. • Step 2: Initiate the first stability analysis with a nearly pure water phase as the trial phase. If the result shows that the mixture is stable, then the mixture contains no aqueous phase, and so go to step 3. If the result shows that the mixture is unstable, go to step 4. • Step 3: Initiate the stability analysis with Wilson's correlation and its inverse. If two TPDs corresponding to nontrivial solutions are both positive, the mixture is single-phase; if no less than one initiate estimate leads to a negative TPD, then the converged K-value corresponding to the lowest TPD is selected to initiate a two-phase splitting calculation and obtain a liquid−vapor result. • Step 4: Perform the two-phase splitting calculation with the converged nontrivial K-values from the first stability analysis initiated with a nearly pure water phase. BIPs for the aqueous phase are assigned for the second phase (aqueous phase), and the nonaqueous phase BIPs are assigned for the first phase (hydrocarbon phase) when computing component fugacities in each phase. • Step 5: The second stability analysis is applied to the hydrocarbon phase. The initial guesses are as in eq 14. If the hydrocarbon phase is stable, the mixture is twophase. Otherwise, go to the next step. • Step 6: Perform a three-phase splitting calculation initiated with the converged nontrivial K-values corresponding to the lowest TPD from the second stability. Aqueous phase BIPs are assigned for the third phase (aqueous phase), and nonaqueous phase BIPs are assigned for the first two phases (hydrocarbon phases). One advantage of our work is that no matter whether vapor−liquid phases or liquid−liquid phases are identified by the second stability analysis, it is impossible to result in the existence of two vapor phases in equilibrium after the threephase splitting calculation, which ensures the reasonable results. If the result of the first stability analysis initiated with a nearly pure water phase as the trial phase shows that the mixture is stable, then the aqueous phase is not present. The system is only comprised of hydrocarbon phases. Wilson's correlation and its inverse are sufficient to identify whether the hydrocarbon mixture is two-phase or single-phase.
Because equal fugacity is a necessary but not sufficient condition for phase equilibrium, a poor initial guess corresponding to the lowest Gibbs free energy may lead to an erroneous flash solution. In most cases, the best initial estimates for flash calculation are K-values corresponding to the lowest TPD in the stability analysis, but it is not the whole picture. If local minima and the global minimum have close Gibbs free energy, and more than one stationary point is located corresponding to negative TPDs, then the K-values corresponding to the lowest TPD may cause erroneous flash results. If this happens, the K-values corresponding to the next lowest TPD from stability analysis are tried to start the ensuing flash calculation (Figure 1).

■ RESULTS AND DISCUSSION
Mutual Solubility. To quantify the mutual solubility of water and hydrocarbons, we examine water solubility in the oleic phase and hydrocarbon solubility in the aqueous phase at different temperatures in the three-phase region. Compared with Henry's law, hydrocarbon solubility in the aqueous phase is examined through a synthetic quaternary mixture taken from Mohebbinia et al., 33 which consists of CO 2 , methane, normal hexadecane, and water. The fluid properties are listed in Table  2, and BIPs for nonaqueous phases are provided in Table 3. BIPs for the aqueous phase are calculated according to eqs 4 and 5. Figure 2 shows that methane solubility in the aqueous phase increases to 1% with temperature on the order of 10 −2 . Under Henry's law, there is no gas solubility in the aqueous phase. Figure 3 illustrates that the CO 2 solubility in the aqueous phase calculated through our method is significantly higher, compared with Henry's law and original PR EOS formulation. We evaluate water solubility in the oleic phase through our method in comparison with Henry's law. This example uses a synthetic oil mixture taken from Luo and Barrufet, 52 which consists of 0.25 water, 0.15 pseudocomponent 1 (PC1), 0.15 pseudocomponent 2 (PC2), 0.2 pseudocomponent 3 (PC3), and 0.25 pseudocomponent 4 (PC4). 53 The component characteristics are described in Table 10. Since the BIPs for the nonaqueous phases used in Luo and Barrufet 52 are not given, we use the BIPs between water and pseudocomponents for the nonaqueous phases from    Zhu and Okuno, 54 described in Table 11. The pressure ranges are selected where the gas/oil/water three phases exist at the given temperatures. Simulation results are shown in Figure 4.
With temperature elevated, the solubility of water in the oleic phase increases up to 0.6 at 480K so that the solubility is nonnegligible in a thermal simulation process. There is an excellent agreement of water solubility in the oleic phase between Henry's law and our method at moderate temperatures and low pressures. However, it turns out that the gas/oil/water three-phase flash using Henry's law leads to a lower water solubility at elevated temperatures than our work, and then it predicts a lower oleic production rate in a thermal simulator. So Henry's law is suitable for low temperatures but not for high temperatures in thermal recovery processes. Power of our Algorithm. Two examples are presented to validate two-phase splitting calculation and three-phase splitting calculation in comparison with results from the dissertation of Varavie 4 and CMG-WinProp. The algorithm adopted in the dissertation of Varavie is a free-water threephase flash, under the assumption that water is a separate phase. The routine in CMG-WinProp employs Henry's law to represent the aqueous phase. The two-phase equilibrium calculation is carried out at a pressure of 34.5 bar and a temperature of 333.15 K. The input data are listed in Tables 4   and 5, and the results of the two-phase equilibrium are shown in Table 6. The agreement of two-phase equilibrium calculation among our works and other methods is excellent. There is no hydrocarbon solubility in the aqueous phase since the feed mixture contains no hydrocarbon gas at the test temperature. It is shown that water solubility calculated by our method is slightly higher than those from two other methods, which validates a more accurate representation of water solubility in hydrocarbon-rich phases provided by our approach. The pressure and the temperature of the threephase splitting calculation are performed at 13.79 bar and 366.5 K. Tables 7, 8, and 9 present the input data and the three-phase splitting results. It can be observed that the concentration of the light hydrocarbon component in the aqueous phase calculated by our method is higher than that by the two other approaches. This is further validation of involving gas solubility in the aqueous phase in our work. Furthermore, the concentration of the water component in the oleic phase calculated by our work is much higher than those obtained with the other three-phase flash methods; a similar conclusion can be drawn for the vapor phase.
Effects of Temperature and Overall Water Mole Fraction. We test our algorithm across a large parameter space to validate the robust calculations of phase compositions in difficult cases. This case uses the mixture taken from Zhu and Okuno. 54,55 The mixture consists of 50% pseudocomponent heavy crude oil and 50% water. The pure-component properties are listed in Table 10. The binary interaction parameters (BIPs) for nonaqueous phases are shown in Table  11, and the binary interaction parameters (BIPs) for aqueous phases are assigned through eq 4. Figure 5 shows the phase diagram of the water/heavy oil mixture. There are no convergence problems even in the vicinity of phase boundaries. Figure 6 indicates that phase fractions are sensitive to temperature. At 30 bar, the vapor phase fraction increases from 0 to 54%, and the aqueous phase decreases from 49% to 0 with the temperature elevated in the three-phase region. The reason behind this phenomenon is correlated with the variations in component distributions in each phase. As shown in Figures 7  and 8, the water concentration increases in both the oleic phase and the vapor phase as the temperature increases. Therefore, the aqueous phase fraction falls drastically. The aqueous phase is dominated by 99.9% water component, and so its composition is not shown here. The concentration of PC1 in the oleic phase and the vapor phase reduces significantly with the increasing temperature, but the other three components change slightly in both phases.
Another set of equilibrium calculations are performed on the water/heavy oil mixture to demonstrate the effect of overall water mole fraction in feed composition on phase distributions in the three-phase region. The flash calculations are performed at 15 bar, with variations in temperatures and overall water mole fractions. With increasing temperature, the oleic phase fraction decreases slightly; the vapor phase fraction increases, but the aqueous phase fraction drops. In the heating process, the water component in the aqueous phase and the light hydrocarbons in the oleic phase become volatile, and phase transitions occur with the elevated temperature. When the temperature is less than 390 K, the oleic phase transits into the vapor phase, and the aqueous phase seldom changes. The variations in both oleic and vapor phase fractions are less than 10% because the overall mole fraction of volatile hydrocarbon   components in the feed mixture is small. However, if the temperature is over 390 K, the vapor phase fraction increases significantly mainly because of the noticeable decrease in the amount of the aqueous phase. The phase amounts change dramatically when the temperature is more than 390 K since the overall water mole fraction is larger than that of the light hydrocarbon component in the feed composition. As is well known, the boiling temperature of water is higher than that of the light hydrocarbon component, and that is why a significant decrease in the aqueous phase mole fraction occurs after the temperature reaches 390 K. We test our method on the basis of a compositional range by combining different fractions of water with the above hydrocarbon mixture taken from Zhu and Okuno. 54,55 The overall water mole fraction in the system is increased from 40 to 60%, to evaluate the effect of water mole fraction on phase distributions. The increase in overall water mole fraction, as is the case in steam flooding, shifts the vapor phase fraction curve toward higher temperatures. That is because water is introduced into the feed system in the form of vapor. With the addition of water, the oleic phase fraction drops, but the aqueous phase fraction mounts at a certain temperature ( Figure 9).     Figure 5. PT diagram for the water/heavy crude oil mixture. V is the vapor phase, L is the oleic phase, and W is the aqueous phase.

■ CONCLUSIONS
A simple but robust algorithm was developed based on the modified PR EOS to describe the gas/oil/water three-phase behaviors of water-containing systems. The water vapor curve is used to predict whether or not an aqueous phase may be present, and accordingly, one of the two sets of initial estimates is chosen to initiate the first stability analysis. We employ the initial estimates of the second stability analysis from Li and Firoozabadi. 47 The developed algorithm is evaluated on a wide range of variations in temperature, pressure, and composition to validate the power of the algorithm.
A standard is stipulated that if the water component accounts for no less than 50%, the mixture is taken as waterlike and binary interaction parameters (BIPs) for the aqueous phase are used for evaluation of component fugacity coefficients; otherwise, BIPs for the nonaqueous phase are selected.
The mutual solubility of water and hydrocarbons caused by high temperatures is shown.
The robustness of the proposed algorithm is shown by several examples. Convergence is obtained without problems even close to the phase boundaries.
The effects of temperature and overall water mole fraction on phase distributions are displayed. As the temperature increases, the vapor phase fraction mounts; with the addition of water, the amount of oleic phase decreases. Specific variations are dependent on the feed composition and current conditions.