Simulation Study on the Influence of Gas Mole Fraction and Aqueous Activity under Phase Equilibrium

This work explored the influence of gas mole fraction and activity in aqueous phase while predicting phase equilibrium conditions. In pure gas systems, such as CH4, CO2, N2 and O2, the gas mole fraction in aqueous phase as one of phase equilibrium conditions was proposed, and a simplified correlation of the gas mole fraction was established. The gas mole fraction threshold maintaining three-phase equilibrium was obtained by phase equilibrium data regression. The UNIFAC model, the predictive Soave-Redlich-Kwong equation and the Chen-Guo model were used to calculate aqueous phase activity, the fugacity of gas and hydrate phase, respectively. It showed that the predicted phase equilibrium pressures are in good agreement with published phase equilibrium experiment data, and the percentage of Absolute Average Deviation Pressures are given. The water activity, gas mole fraction in aqueous phase and the fugacity coefficient in vapor phase are discussed.


Introduction
Gas hydrate is a non-stoichiometric crystalline compound, which consists of a lattice formed by hydrogen bonds of water molecules as the host under high pressure and low temperature.Some gas molecules, such as methane, nitrogen, carbon dioxide and propane, are as the guest firmly surrounded by the crystal network formed by hydrogen bond of water molecules.The ice-like structure enables and stabilizes the existence of gas hydrates at higher temperatures and elevated pressures.The guest molecules must be of correct size to fit inside and stabilize the crystal lattice via weak van der Waals forces with the host water molecules [1,2].
Gas hydrate technology can be applied in many applications, such as gas storage and transportation [3], gas separation [4], refrigeration [5,6], etc.Meanwhile, in response to increasing carbon emissions, the hydrate-based gas separation (HBGS) has attracted the interest of researchers as an effective CO 2 capture and storage (CCS) technology [7].In the last couple of years, vast quantities of natural gas hydrate in the permafrost and deep seabed was found, which is twice as much as the amount of the other fossil fuels combined under a conservative estimate [8]; it makes natural gas hydrate as a kind of potential energy possible.However, gas hydrates can block oil and gas pipelines with high pressure and low temperature inside subsea oil and gas flow line [9].Furthermore, the methane trapped in gas hydrates is a potent greenhouse gas [10].In order to solve these problems, scholars conducted a lot of studies and found that adding thermodynamic inhibitors can effectively change the conditions of hydrate formation into higher pressure and lower temperature.
On the contrary, the formation conditions of hydrate can be changed to lower pressure and higher temperature by adding thermodynamic promoters.Regardless of whether inhibiting or promoting hydrate formation, it is necessary to predict phase equilibrium conditions for the above-mentioned applications, and it is important to use reliable and accurate predictive models for predicting hydrate phase equilibria [11].
Among the gas hydrate predictive models of three-phase equilibrium, there are two common thermodynamic prediction models for calculating the phase equilibrium conditions.One is the classical statistical thermodynamic model of van der Waals and Platteeuw [12].Some improved models designed for the phase equilibrium conditions in distilled water and aqueous solutions were proposed by Nasrifar et al. [13], Haghighi et al. [14], and Martin and Peters [15].Another hydrate model is the Chen-Guo model [16,17], based on equality of the fugacity in the hydrate and vapor phases, which assumed the activity of water to be in unity and neglected the influence of gas solubility in aqueous phase.However, the changes in water activity caused by the solubility of gases, especially acid gases (e.g., carbon dioxide and hydrogen sulfide), and the addition of inhibitors/promoters should not be ignored [18,19].Therefore, a better precision activity of water in aqueous phase is the key to improving the Chen-Guo model.Ma et al. [20,21] used the Patel-Teja equation of state (PT EoS), coupled with the Kurihara mixing rule and the one-fluid mixing rule to calculate the water activity in aqueous phase.Sun and Chen [18] combined the modified method introducing the Debye-Hückel electrostatic contribution term with the PT EoS to predict the nonideality of aqueous phase including ionic components.Liu et al. [22] built a simple correlation to calculate the activity of water in methanol-water solutions of sour gases (CH 4 /CO 2 /H 2 S/N 2 ), in which parameters were determined from the phase equilibrium data.Moreover, among the models of aqueous phase activity, the UNIQUAC model [23] and the modified UNIFAC model [24] were constantly used to calculate aqueous phase.Delavar and Haghtalab [25,26] used the Chen-Guo and UNIQUAC models, referring the Soave-Redlich-Kwong-Huron-Vidal equation of state (SRK-HV EoS) conjunction with the Henry's law, to calculate the gas hydrate formation conditions.Dehaghani and Karami [27] employed the predictive Soave-Redlich-Kwong equation of state (PSRK-EoS) along with the modified Huron-Vidal (MHV1) missing rule and UNIQUAC model to calculate fugacity and activity coefficient of water in equilibrated fluid phases.Klauda and Sander [19,28] applied the modified UNIFAC model and PSRK-EoS coupled with the classical mixing rules, and the results obtained were in good agreement with the experimental data.
However, regardless of using UNIFAC or UNIQUAC model to calculate aqueous phase activity, the gas mole fraction in aqueous phase must be obtained first.Thus, in previous literature [19,[25][26][27][28]], Henry's law was used to describe the gas solubility in aqueous phase, and the gas mole fraction in aqueous phase relied on simultaneous Henry constants, the partial molar volume at infinite dilute (presented by Heidemann and Prausnitz [29]) and vapor phase fugacity calculated by the equation of state.It should be noted that both Henry's law and the partial molar volume at infinite dilution are defined on the basis of an imaginary ideal system.Furthermore, some parameters used in calculating the Henry's law constants and the partial molar volume at infinite dilution are also obtained by regression under the assumed ideal state.As a result, when acid gases exist in the actual system, there is a deviation in the water activity and the gas mole fraction in aqueous phase.Therefore, the gas mole fraction in aqueous phase as a function of temperature and pressure is considered one of the factors that influence the phase equilibrium conditions in this work.
In order to minimize the adverse impacts of the Henry's law constants and infinite dilution partial molar volumes on the hydrate equilibrium conditions prediction, we fitted a correlation of gas mole fraction in aqueous phase according to experimental data.Moreover, the UNIFAC model [30][31][32] and the correlation proposed in this work were employed to calculate the activity coefficient of aqueous phase; PSRK [33][34][35] was used to calculate vapor phase fugacity, and the Chen-Guo model was used to calculate the fugacity of the hydrate phase.These models not only alleviate empirically fitting the intermolecular parameters required in the van der Waals and Platteeuw model but also avoid the calculation error of water activity caused by Henry's law and the infinite dilution partial molar volume.It is noteworthy that the framework proposed in this work only applied in the pure gas (such as CH 4 , CO 2 , N 2 or O 2 ) and distilled water system; the mixed gas system and the additive system will be further studied in future work.Finally, the results calculated are compared with the experimental data in literatures, and the calculated fugacity coefficient of vapor phase and water activity are given.

Thermodynamic Framework
To predict the phase equilibrium conditions of gas hydrate, the iso-fugacity rule is used in the three phases (vapor, water, and hydrate) and a fugacity approach is considered for both gas and water as: where f H i , f L i and f V i are the fugacity of component i including water in the hydrate, liquid and vapor phases, respectively; z, x and y represent the mole fraction of component i in the hydrate, liquid and vapor phases, respectively.The fugacity of water in hydrate phase, f H w , is expressed as: where f MT w represents the fugacity of water in the hypothetical empty hydrate lattice and is assumed equal to the saturated vapor pressure of the empty hydrate lattice [36]; −∆µ MT−  L w is the chemical potential difference calculated by the method of Holder et al. [37]; and R is the universal gas constant.

Thermodynamic Model of Vapor Phase
The PSRK group-contribution method is based on the SRK equation of state [38], which is used to calculate the fugacity of components in vapor phase, as: where P and T are the system pressure and temperature, respectively; v m represents the mole volume, which is obtained by solving the cubic equation derived from Equation (3), and the value is the same as the largest real root of the equation [35]; a and b are parameters of PSRK.The parameters a i and b i of pure component i can be calculated from the critical properties T c,i and P c,i .
where T r = T/T c ; the pure fluid parameter c 1 is taken from the study of Holderbaum and Gmehling [35].
The PSRK mixing rule is written as: where γ i stands for the activity coefficient of component i calculated by UNIFAC model; the recommended value of A 1 = −0.64663 in PSRK model.The activity coefficient γ i is a correction factor that accounts for deviations of real systems from that of an ideal solution, which can be estimated from chemical models (such as UNIFAC).Thus, the fugacity coefficient is given by: where ϕ i is the fugacity coefficient of component i; z = Pv m /RT.

Thermodynamic Model of Hydrate Phase
Chen and Guo [16,17] proposed a two-step hydrate formation mechanism for gas hydrate formation.The following two processes are considered simultaneously in the nucleation process of hydrate: (1) A quasi-chemical reaction process to form basic hydrate, and (2) an adsorption process of smaller gas molecules in the linked cavities of basic hydrate, which leads to the non-stoichiometric property of hydrate.The model is expressed as: where z i denotes the mole fraction of the basic hydrate formed by gas component i, and z i = 1 for pure gas; θ j represents the fraction of the linked cavities occupied by the gas component j; α is the ratio of linked cavities and basic cavities [39], which equals 1/3 for sI hydrates and 2 for sII hydrates, respectively.
where f j denotes the fugacity of component j in vapor phase calculated by PSRK method; c j stands for the rigorous Langmuir constant, which is calculated from the Lennard-Jones potential function.
In Equation (10), f i 0 represents the fugacity of component i in vapor phase in equilibrium with the unfilled pure basic hydrate i (∑θ j = 0) [21].According to the Chen-Guo model, it can be calculated as: where A ij is the binary interaction coefficient which stands for the interplays between gas molecule i in the basic hydrate and gas molecule j in the linked cavities.A i , B i and C i are the Antoine constants, as reported by Chen and Guo [17].β is the function of water volume difference between that in the unfilled basic hydrate phase and the water phase, and the large cavity number per water molecule [20], β = 0.4242 K/bar for sI hydrates, β = 1.0224K/bar for sII hydrates.a w is the activity of water in aqueous phase, which is calculated by the UNIFAC method.For sI and sII hydrates, λ 2 = 3/23 and λ 2 = 1/17, respectively.

Thermodynamic Model of Aqueous Phase
In order to calculate the activity coefficient of components in aqueous phase, the mole fraction of each component in aqueous phase should be obtained first.Therefore, gas is also considered to be a component in aqueous phase.The pressure-corrected Henry's law is employed to calculate the mole fraction of gas component in aqueous phase exclude water as: where subscript i represents the gas component in aqueous phase; H i is the Henry's constant of component i, given by the Krichevsky-Kasarnovsky correlation [36,40]; V i ∞ is the infinite partial molar volume of the component i in water, given by Heidemann and Prausnitz [29].
With the phase equilibrium, the gas mole fraction in aqueous phase can be calculated by the correlation as: Considering the presence of the additive, whether it is a promoter that raises the phase equilibrium temperature (pressure) or an inhibitor that lowers the temperature (pressure), the components in the liquid phase should be recalculated.Furthermore, Delavar and Haghtalab [25,26] point out that the mole fraction of each component in aqueous phases can be calculated as: where i represents water and gas component; p represents the promoter (inhibitor); M i and M p are the molecular weight of component i and the promoter (inhibitor); wt% stands for the weight percentage of the promoter (inhibitor) in aqueous phase.
where i represents water, gas components and the promoter (inhibitor); n i represents the mole fraction of water, gas component and the promoter (inhibitor) in aqueous solutions of a unit mass.The activity coefficient of the components in aqueous phase is calculated by UNIFAC model [30][31][32], which consisting of the combinatorial and residual terms, as: where γ C i and γ R i stand for the combinatorial term and residual term of component i, respectively.The combinatorial term takes the different sizes and shapes of the molecules into account.
where Z = 10; j represents all components in aqueous phase; ϕ i and θ i are the volume and surface area fraction of component i, respectively; r i and q i are the volume and surface area parameters of component i, respectively.They are calculated by the van der Waals volumes R k and surface areas Q k of the individual group k using equations as follows: Processes 2019, 7, 58 6 of 15 where k is the number of group k in component i.The volume R k parameters and surface areas parameters Q k of group k are listed in Table 1.The residual term of component i in Equation ( 18) is replaced by the solution-of-groups concept [32] as: where k represents all groups in aqueous phase, including water; lnΓ k stands for the residual activity coefficient of functional group k; lnΓ k is the residual activity coefficient of group k in the reference solution containing only component i.Both lnΓ k and lnΓ (i) k are calculated as: where m and n are the summations cover different groups of all components in aqueous phase; i in Equation ( 25) is the same as component i in Equation (22); θ m and X m are the surface area fraction and the mole fraction of group m in the mixture, respectively.The group interaction parameter Ψ nm proposed by Sander et al. [30] is described as: where u nm and u mm are the adjustable group interaction parameters (energy parameters).For each group-group interaction, the two parameters have the relation of u nm = u mn .The gas-gas group interaction-energy parameters u nm and temperature range are given in Table 2.In order to properly describe the temperature dependence gas solubility, the correlation proposed by Sander et al. [30] has been used as follows: where u 0 and u 1 are temperature-independent parameters, shown in Table 3.
Table 3. Constants for the calculation of gas-water interaction-energy parameters in the temperature range 273-348K.

Calculation Procedure
The equations described above were solved using codes generated by MATLAB 2014b (MathWork, Beijing, China).The calculation procedure to obtain the phase equilibrium conditions of gas hydrates is summarized in the schematic flow diagram shown in Figure 1.
Processes 2019, 7, 58 8 of 16 Input T(K),P 0 =0(bar) Cal cul ate the m ole fraction of t he g as co mp on ent x i in aqu eou s p hase by t h e expression reported in this work Do (1), ( 2) and (3) in turn.
( ) ( ) Fit th e exp ression s of th e gas mole fraction in aqueous phase x about temperature and pressure

Results and Discussion
When using the UNIFAC model to calculate the activity of gas components, the mole gas fraction in aqueous phase needs to be obtained first.Thus, the Henryconstant Hi and the infinitely diluted partial molar volume i V ∞ are employed to calculate the gas mole fraction in aqueous phase. ∞

Figure 1.
Calculation procedure for the prediction of phase equilibrium pressures at given temperatures.
The percentage of the Average Absolute Deviation in Pressure (AADP) is calculated as: where N is the number of experimental points; P exp i and P cal i stand for the experimental and calculated pressure, respectively.

Results and Discussion
When using the UNIFAC model to calculate the activity of gas components, the mole gas fraction in aqueous phase needs to be obtained first.Thus, the Henryconstant H i and the infinitely diluted partial molar volume V i ∞ are employed to calculate the gas mole fraction in aqueous phase.Heidemann and Prausnitz [29] provided a correlation for solving V i ∞ as follows: where c 11 represents the cohesive energy for water, which was evaluated at each temperature from thermodynamic properties tabulated; h 0 is the molar enthalpy at the given temperature but at zero pressure, and v s w is the molar volume of the saturated liquid [29].However, it is not convenient to get the infinitely diluted partial molar volume of gas in the actual system by the parameters regressed from the assumed ideal state system, especially in the industrial application.Moreover, the parameters are not available for a liquid of unknown components.Therefore, as one of the factors affecting the formation of gas hydrate, the gas mole fraction in aqueous phase cannot be obtained accurately.
As seen from Equation ( 14), when gas type was given, the gas mole fraction in aqueous phase is only a function of temperature and pressure in phase equilibrium.Furthermore, the gas mole fraction in aqueous phase should be a fixed value in the three-phase equilibrium, which is related to the phase equilibrium temperature and pressure.Therefore, when the equilibrium temperature and pressure are given, the gas mole fraction in aqueous phase can be found by interval search using the framework mentioned in this work.Remarkably, in the process of numerical calculation, due to the unreasonable search step length and the inevitable error of experimental data, a set of phase equilibrium temperature and pressure may correspond to multiple gas mole fraction values.In this case, the average of these values can be taken as the gas mole fraction in aqueous phase under the phase equilibrium.As a result, the correlation of the gas mole fraction in aqueous phase is fitted by temperature and pressure, which is determined by the experimental data and defined as: where a, b and c are constants for gas component i + water system, given in Table 4.Although the correlation of the gas mole fraction in aqueous phase is simple and is multivariate linear form, its precision and calculation accuracy are satisfactory.The number of data points used in fitting Equation ( 31) and the R-Square are given in Table 4.
As shown in Figure 2, there are two kinds of tendencies for the fugacity coefficient in methane system with the increase of temperature.First, when the temperature is lower than 286.5 K, the fugacity coefficient decreases with the increase of temperature, which has the same tendency with the fugacity coefficient in the carbon dioxide system, although the reduction rate is smaller.On the other hand, when the temperature is higher than 286.5 K, the fugacity coefficient increases exponentially with the temperature increment, which can also be seen in nitrogen and oxygen systems, as displayed in Figure 2. It should be pointed out that the fugacity coefficient decreases with the temperature increment in carbon dioxide system.This indicates that carbon dioxide is more likely to yield to pressure when the temperature is below the critical temperature.Moreover, the exponential growth of the fugacity coefficient in the N 2 and O 2 systems is mainly related to the temperature increment.
Processes 2019, 7, 58 10 of 16 systems, as displayed in Figure 2. It should be pointed out that the fugacity coefficient decreases with the temperature increment in carbon dioxide system.This indicates that carbon dioxide is more likely to yield to pressure when the temperature is below the critical temperature.Moreover, the exponential growth of the fugacity coefficient in the N2 and O2 systems is mainly related to the temperature increment.The mole fraction of CH4, CO2, N2 and O2 in aqueous phase under phase equilibrium condition is shown in Figure 3-5.In this work, the gas mole fraction was considered as one of the factors affecting the phase equilibrium.It represents the ratio of the number of gas molecules maintaining the three-phase equilibrium to the number of all molecules in aqueous phase.The mole fraction of CH 4 , CO 2 , N 2 and O 2 in aqueous phase under phase equilibrium condition is shown in Figures 3-5.In this work, the gas mole fraction was considered as one of the factors affecting the phase equilibrium.It represents the ratio of the number of gas molecules maintaining the three-phase equilibrium to the number of all molecules in aqueous phase.
In Figures 3-5, under phase equilibrium condition, the gas mole fraction in aqueous phase decrease with temperature increment, and all the changing ranges are less than 1 × 10 −3 .Therefore, there may exist a threshold value for the gas mole fraction in aqueous phase.In other words the hydrate will form when the gas mole fraction in aqueous phase reaches a certain threshold value.Furthermore, for methane hydrate, the results in this work are in agreement with the views of Walsh et al. [47] and Guo and Rodger [48].Walsh et al. suggested that the threshold value of gas mole fraction triggering hydrate formation calculated by the molecular dynamics (MD) simulation was 1.5 × 10 −3 .The threshold value is a reasonable explanation for reducing the temperature or increasing the pressure, which could effectively promote the formation of hydrate.This is because lowering the temperature or increasing the pressure will enhance the gas dissolution, which in turn causes the gas mole fraction in aqueous phase exceeding the threshold value, then hydrate forms.The mole fraction of CH4, CO2, N2 and O2 in aqueous phase under phase equilibrium condition is shown in Figure 3-5.In this work, the gas mole fraction was considered as one of the factors affecting the phase equilibrium.It represents the ratio of the number of gas molecules maintaining the three-phase equilibrium to the number of all molecules in aqueous phase.In particular, although carbon dioxide has a large solubility in water, the mole fraction of carbon dioxide calculated in this work is not very large under phase equilibrium, as shown in Figure 4.It implied that the total number of carbon dioxide gas molecules in aqueous phase for maintaining the three-phase equilibrium is not large, and it may be much smaller than the sum of the gas molecules dissolved in the water.This is mainly because part of carbon dioxide gas molecules dissolved in water turn into carbonic acid, thus reducing the amount of carbon dioxide gas molecules existing in water.Meanwhile, the pH of aqueous phase will be changed, and affecting the activity of the water and phase equilibrium conditions.The gas molar fraction in aqueous phase (×10 In Figure 3-5, under phase equilibrium condition, the gas mole fraction in aqueous phase decrease with temperature increment, and all the changing ranges are less than 1 × 10 −3 .Therefore, there may exist a threshold value for the gas mole fraction in aqueous phase.In other words the hydrate will form when the gas mole fraction in aqueous phase reaches a certain threshold value.Furthermore, for methane hydrate, the results in this work are in agreement with the views of Walsh et al. [47] and Guo and Rodger [48].Walsh et al. suggested that the threshold value of gas mole fraction triggering hydrate formation calculated by the molecular dynamics (MD) simulation was 1.5 × 10 −3 .The threshold value is a reasonable explanation for reducing the temperature or increasing the pressure, which could effectively promote the formation of hydrate.This is because lowering the temperature or increasing the pressure will enhance the gas dissolution, which in turn causes the gas mole fraction in aqueous phase exceeding the threshold value, then hydrate forms.

The water activity in aqueous phase
In particular, although carbon dioxide has a large solubility in water, the mole fraction of carbon dioxide calculated in this work is not very large under phase equilibrium, as shown in Figure 4.It implied that the total number of carbon dioxide gas molecules in aqueous phase for maintaining the In Figure 3, water activity increases with the increase of temperature, in the methane system, reaching a maximum value of 0.985 when the temperature is about 278 K, and then decreases rapidly with the increase of temperature.The general variation trend of the water activity in carbon dioxide system is similar to that of methane system, as seen in Figure 4.In addition, the water activity of the CO 2 system reached its maximum value at about 279.5 K.However, the maximum water activity in the CO 2 system is only about 0.5658, which is probably because of the effect of the carbonic acid.Nevertheless, the activity of water in aqueous phase decreases almost linearly with temperature increase in nitrogen and oxygen systems, as shown in Figure 5.
Figures 6 and 7 show the experimental and predicted phase equilibrium conditions for the single gas hydrate systems.The temperature range, pressure range and AADP (%) are listed in Table 5.
three-phase equilibrium is not large, and it may be much smaller than the sum of the gas molecules dissolved in the water.This is mainly because part of carbon dioxide gas molecules dissolved in water turn into carbonic acid, thus reducing the amount of carbon dioxide gas molecules existing in water.Meanwhile, the pH of aqueous phase will be changed, and affecting the activity of the water and phase equilibrium conditions.
In Figure 3, water activity increases with the increase of temperature, in the methane system, reaching a maximum value of 0.985 when the temperature is about 278 K, and then decreases rapidly with the increase of temperature.The general variation trend of the water activity in carbon dioxide system is similar to that of methane system, as seen in Figure 4.In addition, the water activity of the CO2 system reached its maximum value at about 279.5 K.However, the maximum water activity in the CO2 system is only about 0.5658, which is probably because of the effect of the carbonic acid.Nevertheless, the activity of water in aqueous phase decreases almost linearly with temperature increase in nitrogen and oxygen systems, as shown in Figure 5.
Figure 6 shows the experimental and predicted phase equilibrium pressures for CH 4 and CO 2 .It can be seen the predicted results for all the gas systems are in excellent agreement with the experimental data.It should be noted that the type of carbon dioxide hydrate structure was set to sI, and, because the carbon dioxide gas molecule is too big to be encaged in the linked cavities, the filling rate of the gas molecules in the linked cavities, θ j , was set to 0, as described by Chen and Guo [16].Exp.pure N 2 [43] Exp.pure N 2 [44] Exp.pure N 2 [45] Calc.N 2 Exp.pure O 2 [43] Exp.pure O 2 [44] Exp.pure O 2 [46] Calc.O 2 Pressure (bar) Temperature (K) Figure 7. Experimental and predicted phase equilibrium conditions for N 2 /O 2 + water systems.van Cleeff and Diepen [43], squares , ; Mohammadi et al. [44], triangles , ; Duc et al. [45], stars ; van Cleeff and Diepen [46], circles .
The experimental and predicted phase equilibrium pressures for N 2 and O 2 are displayed in Figure 7.The predicted phase equilibrium pressures are in good agreement with the experiment.It is especially noteworthy that, when calculating oxygen and nitrogen hydrate, the hydrate structure was set to sII, which was based on the ideas proposed by Chen and Guo [16].This is because the gas molecules of N 2 and O 2 are small and have a high filling rate in the connected cavities.
The gas mole fraction in Figure 8 was obtained by inverse phase equilibrium data using the framework proposed in this work.The gas mole fraction threshold value for maintaining the three-phase equilibrium state is different to the critical gas concentration.The gas mole fraction threshold value calculated in this work does not contradict the critical gas concentration proposed by Zhang et al. [49].They pointed out that there is a critical gas concentration in aqueous phase that can spontaneously nucleate in the induction period, and the critical gas concentration is calculated by the total amount of carbon dioxide consumed in vapor phase until hydrate nucleation.
Processes 2019, 7, 58 13 of 16 Figure 6 shows the experimental and predicted phase equilibrium pressures for CH4 and CO2.It can be seen the predicted results for all the gas systems are in excellent agreement with the experimental data.It should be noted that the type of carbon dioxide hydrate structure was set to sI, and, because the carbon dioxide gas molecule is too big to be encaged in the linked cavities, the filling rate of the gas molecules in the linked cavities, θj, was set to 0, as described by Chen and Guo [16].
The experimental and predicted phase equilibrium pressures for N2 and O2 are displayed in Figure 7.The predicted phase equilibrium pressures are in good agreement with the experiment.It is especially noteworthy that, when calculating oxygen and nitrogen hydrate, the hydrate structure was set to sII, which was based on the ideas proposed by Chen and Guo [16].This is because the gas molecules of N2 and O2 are small and have a high filling rate in the connected cavities.
The gas mole fraction in Figure 8 was obtained by inverse phase equilibrium data using the framework proposed in this work.The gas mole fraction threshold value for maintaining the threephase equilibrium state is different to the critical gas concentration.The gas mole fraction threshold value calculated in this work does not contradict the critical gas concentration proposed by Zhang et al. [49].They pointed out that there is a critical gas concentration in aqueous phase that can spontaneously nucleate in the induction period, and the critical gas concentration is calculated by the total amount of carbon dioxide consumed in vapor phase until hydrate nucleation.However, the phase equilibrium data cited in this work were recorded at the end of decomposition rather than in the preliminary stage of nucleation.Gas molecules entrapped in hydrate cage cannot be released totally, which was owing to memory effect [50].Moreover, there theoretically exists a concentration difference as a force in mass transfer during hydrate nucleation and decomposition.Therefore, the gas mole fraction threshold value calculated in this work is less However, the phase equilibrium data cited in this work were recorded at the end of decomposition rather than in the preliminary stage of nucleation.Gas molecules entrapped in hydrate cage cannot be released totally, which was owing to memory effect [50].Moreover, there theoretically exists a concentration difference as a force in mass transfer during hydrate nucleation and decomposition.Therefore, the gas mole fraction threshold value calculated in this work is less than the critical gas concentration.
In Figure 8, the threshold value of gas mole fraction achieves a maximum of 5.61 × 10 −3 at 0 • C. A possible reason is that part of carbon dioxide molecules in aqueous phase react with water to form carbonic acid.When the temperature is above 0 • C and below 0 • C, the pressure increment and the temperature decrement become a dominant factor that results in more stability for the carbonic acid and less solubility of carbon dioxide, respectively.However, this analysis should be proved by further study.
Furthermore, since the correlation of gas mole fraction fitted in this work is a multivariate linear form, the trend of the gas mole fractions in Figures 4 and 8 are different.Therefore, a large number of accurate and reliable experiment data can effectively improve the prediction accuracy of the model in this work.

Conclusions
In this work, the Chen-Guo model coupled with the PSRK method were employed to predict phase equilibrium conditions of CH 4 , CO 2 , N 2 or O 2 in pure water systems.The gas mole fraction in aqueous phase is one of the factors that affect the phase equilibrium of gas hydrate proposed in this work.The gas mole fraction threshold value maintaining the three-phase equilibrium was obtained by reversed phase equilibrium data.Meanwhile, in order to obtain the water activity in aqueous phase, the correlation of the gas mole fraction threshold value in aqueous phase was fitted though UNIFAC model.The calculated water activity can effectively improve the accuracy of the prediction results, and the predicted results of this work are in good agreement with the experimental data reported in the references.

Save
experiment data T(K) and P(bar); x 0 = 0 (1):Calculate the activity coefficient of component i, including water, through UNIFAC model (2):Calculate the fugacity coefficient of component i in the gas phase f i V by the PSRK model (3) :Calcu la te the fugacity of component i in the hydrate phase f i H by the Chen-Guo model, with the specified structure, f i V and the activity of water a w .

Figure 1 .
Figure 1.Calculation procedure for the prediction of phase equilibrium pressures at given temperatures.

Figure 2 .
Figure 2. The fugacity coefficient of CH 4 , CO 2 , N 2 and O 2 in vapor phase.

Figure 3 .
Figure 3.The mole fraction of CH4 and water activity in aqueous phase.

Figure 4 . 3 )
Figure 4.The mole fraction of CO2 and water activity in aqueous phase.

Figure 4 .
Figure 4.The mole fraction of CO 2 and water activity in aqueous phase.

Figure 4 .
Figure 4.The mole fraction of CO2 and water activity in aqueous phase.

Figure 5 .
Figure 5.The mole fraction of N2/O2 and water activity in aqueous phase.

Figure 5 .
Figure 5.The mole fraction of N 2 /O 2 and water activity in aqueous phase.

Figure 8 .
Figure8.The experimental data and the mole fraction of carbon dioxide in aqueous phase for CO2 + water systems under the phase equilibrium.The experimental data were reported by Ma et al.[42].

3 )Figure 8 .
Figure8.The experimental data and the mole fraction of carbon dioxide in aqueous phase for CO 2 + water systems under the phase equilibrium.The experimental data were reported by Ma et al.[42].

Table 1 .
The UNIFAC group volume and surface-area parameters.

Table 2 .
Gas-gas group interaction-energy parameters u nm and temperature range.
a 56, 57, 58 and 60 are the group numbers of CO 2 , CH 4 , O 2 and N 2 in the UNIFAC group parameter list, respectively.

Table 4 .
Parameters for the correlation of the gas mole fraction in aqueous phase.

Table 5 .
The phase equilibrium pressure and temperature range of experimental data and the Average Absolute Deviation in Pressure (AADP) for predicted results.