Influence of N2 on Formation Conditions and Guest Distribution of Mixed CO2 + CH4 Gas Hydrates

In this contribution, a method based on a solid solution theory of clathrate hydrate for multiple cage occupancy, host lattice relaxation, and guest-guest interactions is presented to estimate hydrate formation conditions of binary and ternary gas mixtures. We performed molecular modeling of the structure, guest distribution, and hydrate formation conditions for the CO2 + CH4 and CO2 + CH4 + N2 gas hydrates. In all considered systems with and without N2, at high and medium content of CO2 in the gas phase, we found that CO2 was more favorable in occupying clathrate hydrate cavities than CH4 or N2. The addition of N2 to the gas phase increased the ratio concentration of CO2 in comparison with the concentration of CH4 in clathrate hydrates and made gas replacement more effective. The mole fraction of CO2 in the CO2 + CH4 + N2 gas hydrate rapidly increased with the growth of its content in the gas phase, and the formation pressure of the CO2 + CH4 + N2 gas hydrate rose in comparison to the formation pressure of the CO2 + CH4 gas hydrate. The obtained results agreed with the known experimental data for simple CH4 and CO2 gas hydrates and the mixed CO2 + CH4 gas hydrate.


Introduction
Due to concerns of an increasing global warming effect, the capture from industrial flue gas and long-term storage of carbon dioxide are among the most important challenges facing the world scientific community today.
Various strategies [1] to sequester carbon dioxide have been proposed, but a technology for large-scale and safely stored CO 2 has not been completely developed.
Currently, the leading approach to this problem involves the injection of CO 2 into depleted deep underground natural gas reservoirs [2]. In seismically active zones such as Japan, the use of this method of CO 2 storage may lead to gas leakage due to geological perturbations such as earthquakes or fractures. Another approach offered by Ohgaki and Inoue [3] is the sequestration of CO 2 as solid hydrates through the formation of CO 2 clathrate hydrate. Recently, the injection of CO 2 into porous sediments at a depth of several hundred meters below the deep ocean floor has been proposed as an alternative long-term sequestration option that would be resistant to geophysical perturbations [4]. Such deposition of CO 2 prevents the transport of the CO 2 back to the surface due to the formation of a CO 2 clathrate hydrate capping layer that reduces the migration. The authors estimate this storage strategy could remain intact for millions of years. It has also been shown that CO 2 hydrates have anomalously low dissociation rates at atmospheric pressure. This self-preservation effect takes place in the temperature range 245-271 K, which could be significant in practice for the CO 2 storage in the form of clathrate hydrates [5].
Clathrate hydrates are nonstoichiometric inclusion compounds consisting of water (host) molecules forming a crystalline framework in which cavities (cages) guest molecules can be included. In nature, three types of gas hydrates are common: Hydrates of cubic structure (CS)-I, cubic structure CS-II, and hexagonal structure (HS)-III [6]. The hydrate structure is determined primarily by the size of guest molecules. Thus, large guest molecules such as propane and isobutane as well as small guests, in particular oxygen, nitrogen, and hydrogen, form the cubic structure CS-II. Guest molecules of intermediate size such as methane, xenon, and carbon dioxide form hydrates of the cubic structure CS-I. The presence of two types of guest molecules is necessary to form the hexagonal structure HS-III: A very large molecule (e.g., 2,2-dimethylbutane) and a small molecule (e.g., methane). The hydrate structures differ by size and number of cavities in their unit cells.
Large amounts of natural gas hydrates are composed mainly of methane in the form of solid hydrates stored on continental margins and in permafrost regions [6]. Technologies involving the simultaneous production of raw hydrocarbon and greenhouse-gas sequestration is promising. Carbon dioxide sequestration in deep-sea sediments or permafrost regions can be performed simultaneously with natural gas recovery by swapping hydrocarbon molecules in hydrate cages for carbon dioxide molecules, thus providing a mechanism of hydrocarbons production and greenhouse gas sequestration [7][8][9]. The replacement of CH 4 hydrates by CO 2 hydrates has been studied [7] for recovering CH 4 gas. When a mixture of water with gas or liquid CO 2 itself is put under certain pressure, a solid CO 2 hydrate can be formed at much milder P-T conditions than a CH 4 hydrate can [10]. Thus, the swapping process between two gaseous guests is considered to be a promising approach to long-term storage of CO 2 . When the CH 4 hydrate is put under a certain pressure [8,9] of a CO 2 /N 2 gas mixture, a decomposition of CH 4 hydrates and a solid mixed hydrate containing CO 2 can appear with recovered CH 4 gas. The direct use of a CO 2 + N 2 gas mixture (20 mol% CO 2 and 80 mol% N 2 to reproduce flue gas from a power plant) instead of pure CO 2 greatly enhances the overall CH 4 recovery rate in complex marine systems and reduces the costs [11] of CO 2 separation from flue gas. A great number of experimental and theoretical studies concerning the stability and composition of gas hydrates formed from gas mixtures has been published in the last decade [6,[12][13][14][15][16][17][18][19][20][21][22]. In particular, a comparison of numerous experimental data on phase equilibria in a water-methane-carbon dioxide and a water-nitrogen-carbon dioxide systems were presented in papers [23][24][25][26][27][28]. The stability of CO 2 /N 2 or CO 2 /CH 4 mixed hydrates were studied in various conditions, and it was shown that the three-phase hydrate-water-vapor equilibrium curves were shifted to higher pressures at all considered temperatures with decreases [25][26][27][28] in CO 2 concentration in the vapor phase. The statistical thermodynamic theory of van der Waals and Platteeuw [29] was used for modeling the hydrate phase containing nitrogen [23][24][25]. The main assumptions were made in the original van der Waals-Platteeuw model: It stated that each cavity can contain at most one gas molecule. However, Kuhs and coworkers [30,31] found the first direct evidence for multiple occupancy of the cages in nitrogen hydrate. These results have been confirmed by molecular dynamics simulations of CS-II nitrogen hydrate with varying cage occupancies and at different conditions [32,33]. For the more correct prediction of hydrate phase equilibria, it is necessary to consider the possibility of multiple occupancy of the cages in the gas hydrate containing nitrogen.
The aim of this paper is the investigation of the possibility of recovering methane from methane hydrates using either CO 2 or a CO 2 /N 2 gas mixture. The hydrate phase was treated with the solid solution theory of clathrate hydrate for multiple cage occupancy, host lattice relaxation, and guest-guest interactions [34][35][36][37][38]. With this goal, we determined the dependencies of the compositions of the gas hydrates formed from methane + carbon dioxide, nitrogen + carbon dioxide binary gas mixtures, and methane + carbon dioxide + nitrogen ternary gas mixtures, as well as of the formation conditions of these hydrates in dependence on temperature and pressure for different compositions of the gas phase.

Implemented Theoretical Models
In order to accurately estimate the thermodynamic properties of clathrate hydrates, we developed an approach based on the solid solution theory of clathrate hydrate for multiple cage occupancy, host lattice relaxation, and guest-guest interactions [34][35][36][37][38]. The method is based on only one of several assumptions of the original van der Waals-Platteeuw theory [29]: The free energy of clathrate hydrate does not depend on the arrangement of guest molecules in cavities at fixed values of filling degrees for each definite type of cavity. In this approach, the lattice dynamics method that takes quantum effects into account is used and the crystalline host lattice is considered to be nonrigid, is able to change volume depending on the type of guest molecules, and is permitted to describe first-order phase transitions.
The mathematical formalism of the present model for the general case and in the case of clathrate hydrates with two types of cavities and one type of guest was described in our previous studies [34,35]. In contrast with our previous work, in the present work we formulated our approach for the hydrate having two types of cavities, large (L) and small (S), and with the possibility of single occupancy of small and large cavities by a, b, and c type guests; single occupancy of large cavities by a, b type guests; and multiple occupancy of large cavities by c type guests.
In the mentioned approximation, free energy of the clathrate hydrates could be presented as where: (a) For binary clathrate hydrates with cavities of two types including two types of guest molecules and the possibility of single occupancy for type a guests and single occupancy for type b guests of both the small and large cavities, i l = 1, l = a, b, t = S, L; F 1 is the part of free energy at a given degree of filling of the guest molecules in the S and L cavities; {y} = y a S , y a L , y b S , y b L , y l t = N l t /N t are filling degrees for cavities of the tth type (t = S, L) by guest molecules of the lth type (l = a, b); N t is the number of cavities of the tth type; and N l t is the number of guest molecules of the lth type contained in cavities of the tth type. (b) For ternary clathrate hydrates with cavities of two types including three types of guest molecules and the possibility of single occupancy for type a guests and single occupancy for type b guests of both the small and large cavities, as well as the possibility of single occupancy for type c guests of the small and multiple occupancy large cavities, l = a, b, c, i a = 1, i b = 1, i c = 1, 2, t = S, L; F 1 is the part of free energy at a given degree of filling of the guest molecules in the S and L cavities; {y} = y a S , y a L , y b S , y b L , y c S , y c L , y 2c L , y i l l t = N i l l t /N t is the filling degree for a cavity of the tth type (t = S, L) by guest molecules of the lth type (l = a, b, c); N t is the number of cavities of the tth type; N i l l t is the number of guest molecules of the lth type contained in cavities of the tth type.
In our cases, the binary clathrate hydrates were i l = 1, l = CO 2 , CH 4 , t = S, L, and the ternary clathrate hydrates were i CO 2 = 1, i CH 4 = 1, i N 2 = 1, 2, l = CO 2 , CH 4 , N 2 , t = S, L. In the models, it was considered that the molecules of CO 2 , CH 4 , and N 2 could single occupy both the small and large cavities, whereas molecules of N 2 could also double occupy both the large cavities.
For a given arrangement {y} of the guest molecules in the cavities the free energy F 1 (V, T, {y}) of the crystal could be calculated within the framework of a lattice dynamics approach as where U is the potential energy, and F vib is the vibrational contribution where ω j → q is the jth eigenfrequency of crystal vibration, and → q is the wave vector. Free energy was computed for several values of volume, and it had a minimum corresponding to the equilibrium structure at zero pressure.
The equation of state was found by numerical differentiation of the free energy with respect to volume: Then we found the chemical potentials µ i l l t of guest molecules in the hydrate by numerical differentiation of the free energy with respect to the number of guest molecules: If the free energy F is known, then the Gibbs free energy is As an expression in terms of the chemical potentials of the host and guest molecules, it can be found that where ν t = N t N Q , and N Q is the number of water molecules. The P-T line of monovariant equilibrium of different hydrates and ices could be found from the equality condition of the chemical potentials of water molecules in hydrates and in ice or in the liquid phase: Analogously, the equality of the chemical potentials of guest molecules in the hydrate and gas phases could be written as The chemical potential of guest molecules in the gas phase were calculated using the following equations for a non-ideal gas mixture with a Lennard-Jones interaction between molecules [20]: where ε i and σ i are Lennard-Jones parameters. Interaction parameters between molecules of different types are defined by the combination rules ε ij = √ ε i ε j and σ ij = σ i + σ j 2 :x i = N g i /N g is the mole fraction of the ith component in the gas mixture,N g is the number of guest molecules in the gas phase; N g i is the number of guest molecules of the lth type in the gas phase; and m i is the molar mass of the ith component. The first term in Equation (12) corresponded to the chemical potential of the ideal gas, and the second two corrections appeared for real gases.
The chemical potential of liquid water µ wA q was taken from the model proposed earlier [39] and was given by where T 0 is the initial temperature, and P 0 is the initial pressure. The following constants were also used, as defined in Reference [39]: The Gibbs energy of formation g w 0 L pure , the molar enthalpy of water h wL pure , and the water volume v wL pure , all constants that refer to pure water. Instead of using the empirical value of g w 0 L pure , we calculated this parameter directly by using the lattice dynamic method [38]. The constants were recalculated because those given in Reference [39] were related to standard conditions. We calculated the Gibbs energy and enthalpy at T = 273.15 K (i.e., the temperature of the equilibrium between the ice and water phases, but with pressure remaining at 1 atmosphere). It was connected with the differences between experimental Gibbs energy and that estimated within the model's interaction potential. This parameter could be evaluated from the chemical potential of water, µ wA q , which should be equal to the µ ice 0 of hexagonal ice calculated at the ice I h melting point at standard pressure and temperature. The degrees of cage filling were found from Equations (5) and (11): For determination of the hydrate composition as a function of the gas phase composition, the following relations for fractions of the filled large (C i l l L ) and small (C i l l S ) cages were used, and mole fractions in guest subsystems of l type guest molecules included in the hydrate phase (x h l ) were used:

Simulations Details
The unit cells were chosen as the simulation cell of CS-I (46 water molecules forming 6 large and 2 small cages) and CS-II hydrates (136 water molecules forming 8 large and 16 small cages). Large cages as well as small ones could be filled by one carbon dioxide or methane molecule. The possibility of double filling of large or small cages by these guest molecules was not considered due to the comparatively large size of these molecules. For modeling of ice I h , the simulation supercell containing 32 unit cells (i.e., 128 water molecules) was used. Coulomb interactions were calculated by the Ewald method. The protons were placed according to Bernal-Fowler rules [40], and the water molecules were oriented such that the total dipole moments of the simulation cells of ice and the hydrates were zero with a precision of better than 0.1% of the magnitude of the dipole moment of a single water molecule. The interaction of water-water molecules in hydrates and in ice were described by the modified SPC/E (Simple Point Charge/Extended) potential [41]: where the Lennard-Jones parameters are σ O = 3.1556 Å, ε O = 0.65063 kJ/mole. The usage of other water-water interactions models could influence the absolute values of chemical potentials and slightly influence the ice-gas-hydrate equilibria. Although the chosen water-water interaction model was simple, it was in agreement with known experimental data and was earlier successively applied to consider properties of a number of other simple and mixed hydrates such as CH 4 , C 2 H 6 , Xe, CH 4 + C 2 H 6 , and it allowed us to establish the correctness of our model [34]. Charges on hydrogen atoms were q H = + 0.4238|e| and on oxygen atoms q O = −0.8476|e|. This parameters selection allowed us to reach good agreement with the experimental data [34,35]. For description of the interactions of guest molecules between each other and with water molecules, the Lennard-Jones potential was used with the parameters σ = 3.73 Å, ε = 1.2305 kJ/mole for methane molecules [42], σ = 4.00 Å, ε = 1.5801 (1) kJ/mole for carbon dioxide molecules [43], and σ = 3.6154 Å, ε = 0.844 kJ/mole for nitrogen molecules [44].

Gas-Hydrate Phase Equilibria
Gas-hydrate divariant equilibria are described by Equation (11). This equation represents the conditions of equality of chemical potentials of guest molecules in a hydrate with the gas phase of the same kind of molecules in dependence on pressure and temperature. The comparison of degrees of filing for binary (50% CH 4 and 50% CO 2 in the gas phase) and ternary (15% CH 4 and 15% CO 2 and 70% N 2 in the gas phase) mixed hydrates at the temperature T = 277 K are presented in Figure 1a,b. One can see (Figure 1a) that occupation of both small and large cavities by carbon dioxide molecules was more preferable than by methane. The difference in the degrees of filling was the result of a slightly larger size of CO 2 molecules and stronger interaction with water molecules. After the addition of nitrogen into the gas phase, the tendency was the same (Figure 1b), but in this case cavities' occupation by N2 molecules could be concurred by CO 2 and CH 4 molecules. In spite of more than a two times higher concentration of nitrogen in the gas phase, methane and carbon dioxide molecules more rapidly occupied the large cavities. Thus, at the pressure 10 MPa, only 19% of large cavities were filled by N 2 , whereas 24% and 57% were filled by CH 4 and CO 2 , respectively. A different situation was observed for small cavities filling. In this case, N 2 molecules became preferable and could concur with larger molecules of CH 4 and CO 2 . Therefore, at the same pressure, 10 MPa, the cages were filled by N 2 , CH 4 , and CO 2 in the amounts of 37%, 19%, and 31%, respectively. A tendency of increasing of the cavities number occupied by nitrogen molecules with the pressure increasing was also observable. In both binary and ternary hydrates, one could see the noticeable growth of ratios         Figure 3. The arrows show the equilibrium formation points for hydrates. The intriguing result was that with pressure increasing, the CO2 fraction in hydrate decreased while the CH4 fraction grew. Such behavior  Another interesting finding was that the CO2 fraction decrease rate in the hydrate phase was small for high CO2 concentrations in the gas phase (Figure 2a), and it increased with CO2 concentration decreases. At the temperature T = 273 K (Figure 2a,b), in the pressure interval 1 to 10 MPa for the gas mixture of 90% CO2 and 10% CH4, the change of or was about 0.022 whereas for the gas mixture of 50% CO2 + 50% CH4 it became about 0.082, almost four times larger.
After the addition of nitrogen to the binary carbon dioxide + methane mixture (Figure 3), the amount of these gases stored in the hydrate phase decreased, but not drastically. At the formation pressure (1.8 MPa at T = 273 K) for the ternary gas mixture of 27% CO2 + 3% CH4 + 70% N2, (Figure  3a), the relative fraction of CO2 in the hydrate was found to be 0.847 instead of 0.96 for the binary mixture 90% CO2 + 10% CH4. Methane with its low content in the gas phase (3%) had a fractional content of 0.035. In comparison to the hydrate phase formed from the ternary gas mixture 15% CO2 + 15% CH4 + 70% N2 (Figure 3b), the situation changed notably. At first, methane occupied a notable part of hydrate cavities and the mole fraction of methane reached 0.229 at the formation pressure. Nitrogen content in hydrate became 0.163, almost 40% higher than in the previous case. Methane and nitrogen could replace carbon dioxide in cavities, but CO2 molecules still occupied more than 60% of cavities (mole fraction is 0.614). The increase of relative gas fraction contents in hydrates with pressure was almost equal for different gas phase compositions (0.06 for 27% CO2 and 0.08 for 15% CO2 in gas mixtures).
In all considered systems with and without N2, at high and medium content of CO2 in the gas phase, we found that CO2 was more favorable in occupying clathrate hydrate cavities than CH4 or N2. Moreover, the addition of N2 to the gas phase increased the ratio. For mixtures of 50% CO2 + 50% CH4 and 15% CO2 + 15% CH4 + 70% N2, this ratio increased by 1.5%, and for mixtures of 90% CO2 + 10% CH4 and 27% CO2 + 3% CH4 + 70% N2, by about 1%. Therefore, the addition of N2 made gas replacement more effective.
Another interesting finding was that the CO 2 fraction decrease rate in the hydrate phase was small for high CO 2 concentrations in the gas phase (Figure 2a), and it increased with CO 2 concentration decreases. At the temperature T = 273 K (Figure 2a,b), in the pressure interval 1 to 10 MPa for the gas mixture of 90% CO 2 and 10% CH 4 , the change of x CH 4 or x CO 2 was about 0.022 whereas for the gas mixture of 50% CO 2 + 50% CH 4 it became about 0.082, almost four times larger.
After the addition of nitrogen to the binary carbon dioxide + methane mixture (Figure 3), the amount of these gases stored in the hydrate phase decreased, but not drastically. At the formation pressure (1.8 MPa at T = 273 K) for the ternary gas mixture of 27% CO 2 + 3% CH 4 + 70% N 2 , (Figure 3a), the relative fraction of CO 2 in the hydrate was found to be 0.847 instead of 0.96 for the binary mixture 90% CO 2 + 10% CH 4 . Methane with its low content in the gas phase (3%) had a fractional content of 0.035. In comparison to the hydrate phase formed from the ternary gas mixture 15% CO 2 + 15% CH 4 + 70% N 2 (Figure 3b), the situation changed notably. At first, methane occupied a notable part of hydrate cavities and the mole fraction of methane reached 0.229 at the formation pressure. Nitrogen content in hydrate became 0.163, almost 40% higher than in the previous case. Methane and nitrogen could replace carbon dioxide in cavities, but CO 2 molecules still occupied more than 60% of cavities (mole fraction is 0.614). The increase of relative gas fraction contents in hydrates with pressure was almost equal for different gas phase compositions (0.06 for 27% CO 2 and 0.08 for 15% CO 2 in gas mixtures).
In all considered systems with and without N 2 , at high and medium content of CO 2 in the gas phase, we found that CO 2 was more favorable in occupying clathrate hydrate cavities than CH 4 or N 2 . Moreover, the addition of N 2 to the gas phase increased the x CO 2 x CH 4 ratio. For mixtures of 50% CO 2 + 50% CH 4 and 15% CO 2 + 15% CH 4 + 70% N 2 , this ratio increased by 1.5%, and for mixtures of 90% CO 2 + 10% CH 4 and 27% CO 2 + 3% CH 4 + 70% N 2 , by about 1%. Therefore, the addition of N 2 made gas replacement more effective.

3.2.Gas-Hydrate-Ice (Water)Phase Equilibria
We conducted a calculation of P-T diagrams for gas-hydrate-ice (water) phase equilibria (described by Equations (10) and (11), carried out earlier) for one-component hydrates of methane and reproduced the experimental data with good accuracy [45,46].
The modeling was performed within the molecular model framework described above: Here such calculations were conducted for carbon dioxide hydrates. The obtained lines of the phase equilibria also were in reasonable agreement with the experimental data. In Figure 4, the calculated curves of the ice-gas-hydrate phase equilibria are presented for the considered one-component hydrates of carbon dioxide and methane as well as the available experimental data [6] for comparison. The modeling was performed within the molecular model framework described above: Here such calculations were conducted for carbon dioxide hydrates. The obtained lines of the phase equilibria also were in reasonable agreement with the experimental data. In Figure 4, the calculated curves of the ice-gas-hydrate phase equilibria are presented for the considered one-component hydrates of carbon dioxide and methane as well as the available experimental data [6] for comparison.
The calculation of hydrate formation pressure as well as the CO2 fraction in the hydrate phase in dependence on the gas phase composition were performed for binary CO2 + CH4 mixtures at temperatures of 273 K and 277 K (Figure 5a,b). The temperatures were chosen in order to describe gas equilibria for both gas-hydrate-ice and gas-hydrate-water. It is notable that carbon dioxide molecules filled cages better than methane. For example, at the The calculation of hydrate formation pressure as well as the CO 2 fraction in the hydrate phase in dependence on the gas phase composition were performed for binary CO 2 + CH 4 mixtures at temperatures of 273 K and 277 K (Figure 5a,b). The temperatures were chosen in order to describe gas equilibria for both gas-hydrate-ice and gas-hydrate-water. The modeling was performed within the molecular model framework described above: Here such calculations were conducted for carbon dioxide hydrates. The obtained lines of the phase equilibria also were in reasonable agreement with the experimental data. In Figure 4, the calculated curves of the ice-gas-hydrate phase equilibria are presented for the considered one-component hydrates of carbon dioxide and methane as well as the available experimental data [6] for comparison.
The calculation of hydrate formation pressure as well as the CO2 fraction in the hydrate phase in dependence on the gas phase composition were performed for binary CO2 + CH4 mixtures at temperatures of 273 K and 277 K (Figure 5a,b). The temperatures were chosen in order to describe gas equilibria for both gas-hydrate-ice and gas-hydrate-water.  [47] of mole fraction carbon dioxide in the hydrate phase; crosses = experimental data [47].
It is notable that carbon dioxide molecules filled cages better than methane. For example, at the temperature 273 К, the hydrate formation pressure and equimolar composition of the gas phase (50%  [47] of mole fraction carbon dioxide in the hydrate phase; crosses = experimental data [47]. It is notable that carbon dioxide molecules filled cages better than methane. For example, at the temperature 273 K, the hydrate formation pressure and equimolar composition of the gas phase (50% methane, 50% carbon dioxide), the fraction of CO 2 molecules in hydrate reached 73% and the fraction of CH 4 about 27%. At the temperature 277 K, the fractions of gas molecules in hydrate were 70% and 30%, respectively. The ratio of occupancies by CO 2 and CH 4 was 2.7:1 at 273 K and 2.3:1 at 277 K. These calculation results agreed well with the experimental data [47]. It has to be noted that with increasing temperature, the fraction of CO 2 in hydrate decreased. That could be connected with the increase of pressure, which was necessary for gas hydrate formation. In this case, the methane molecules could concur with the carbon dioxide molecules to occupy mainly small hydrate cavities.
It was connected with the more suitable size of hydrate small cavities for methane molecules in spite of their weaker interactions with surrounding water molecules. In Table 1, the data for binary hydrates formation conditions at several temperatures and gas phase compositions are presented. Calculations were performed for water-gas-hydrate equilibria (T = 277 K) and for ice-gas-hydrate equilibria (T = 273 K, 258 K). Table 1. P-T-x equilibria of ice (water)-gas-hydrate systems for the CO 2 + CH 4 gas mixture.  At relatively low pressures, the solubility of considered gases in ice and water was neglected in our calculations. As could be expected, the equilibrium pressure in systems of gas-hydrate-ice (water) increased with temperature and decreased with a rising amount of carbon dioxide in the gas phase. Analysis of the data also showed nontrivial increasing of the carbon dioxide mole fraction in hydrate with decreasing temperature. That could be connected with the lowering formation pressure while temperature decreased. We could conclude that at low pressure, methane was less favorable than carbon dioxide in cavities occupation.
The formation pressure of binary hydrates rose with increasing methane content in the gas phase and increasing temperature. With temperatures of about 277 K, which corresponds to water temperatures near the bottom of oceans, methane hydrates could form and exist in thermodynamic equilibrium with water and gas at the pressure 4.2 MPa, corresponding to 420 meters in depth (depths of continental slope), whereas for the one-component hydrate of carbon dioxide the formation pressure was lower, 1.2 MPa. At lower temperatures (273 K, 258 K), the decrease of the hydrate formation pressure with the addition of carbon dioxide to methane became not so significant. Thus, with methane content diminishing from 100% to 0%, the change in the formation pressure at T = 273 K was about 1.5 MPa and at T = 258 K, it was about 0.8 MPa.
With an increase of methane content in the gas phase, the hydrate formation pressure gradually rose. The formation pressure of double hydrates of methane and carbon dioxide appeared to be a linear function of the methane content in the hydrate. This was remarkable having in mind the significant difference in interaction strengths between the guest molecules (carbon dioxide-carbon dioxide, methane-methane).
The calculations of the formed hydrates composition at different temperatures showed that for replacement of methane in hydrate by carbon dioxide, the low temperature was preferable.
To understand the influence of additional nitrogen into the carbon dioxide gas phase on hydrate formation conditions and compositions, the dependencies of the formation pressure and mole fraction of CO 2 in hydrate on the gas phase composition for temperatures of T = 272 K (equilibria of gas-hydrate-ice) and T = 274 K (equilibria of gas-hydrate-water), and in the range for CO 2 mole fractions from 0.0 to 1.0 in the gas phase (Figure 6a,b), were found. With an increase of methane content in the gas phase, the hydrate formation pressure gradually rose. The formation pressure of double hydrates of methane and carbon dioxide appeared to be a linear function of the methane content in the hydrate. This was remarkable having in mind the significant difference in interaction strengths between the guest molecules (carbon dioxide-carbon dioxide, methane-methane).
The calculations of the formed hydrates composition at different temperatures showed that for replacement of methane in hydrate by carbon dioxide, the low temperature was preferable.
To understand the influence of additional nitrogen into the carbon dioxide gas phase on hydrate formation conditions and compositions, the dependencies of the formation pressure and mole fraction of CO2 in hydrate on the gas phase composition for temperatures of T = 272 K (equilibria of gas-hydrate-ice) and T = 274 K (equilibria of gas-hydrate-water), and in the range for CO2 mole fractions from 0.0 to 1.0 in the gas phase (Figure 6a,b), were found.
The calculations showed that for all gas phase compositions, the fraction of CO2 of more than 0.035 was higher in the hydrate phase relative to the gas phase, and the hydrate structure CS-I appeared to be more stable than structure CS-II. Whereas for all gas phase compositions the fraction of CO2 less than 0.035 was higher in the hydrate phase relative to the gas phase, the hydrate structure CS-II appeared to be more stable than structure CS-I. That could be connected with the larger largeto-small cavities ratio in the CS-I structure that was more suitable for CO2 molecules. Even at 20% CO2 in the gas phase, the CO2 fraction in hydrate reached 0.75 at T = 272 K and 0.70 at T = 274 K. As one can see from Figure 6b, experimental data were described with a reasonable accuracy, at least for relatively low temperatures.
The absence of phase transition CS-I-CS-II in our results at T = 274 K could be connected with the roughness of the used approximation, in which we did not take into account the solubility of gases on water. It is not very important at low pressure, but is significant at high pressure, when N2 and CO2 solubility rapidly increase. On the other hand, one can neglect the solubility of gases on the ice at formation conditions. Another system we considered was CO2 + CH4 +N2 mixtures, which can form hydrates with water or ice. In the gas mixture, the N2 mole fraction was fixed at 0.7 and the relative content of CO2 and CH4 varied from 0.00 CO2 to 0.30 CO2 in the gas phase. At low Figure 6. Gas-hydrate-ice (water)phase equilibria for carbon dioxide and nitrogen gas mixtures for cubic structure (CS)-I and CS-II hydrates in comparison with the experimental data, skew crosses [25] and crosses [27], The calculations showed that for all gas phase compositions, the fraction of CO 2 of more than 0.035 was higher in the hydrate phase relative to the gas phase, and the hydrate structure CS-I appeared to be more stable than structure CS-II. Whereas for all gas phase compositions the fraction of CO 2 less than 0.035 was higher in the hydrate phase relative to the gas phase, the hydrate structure CS-II appeared to be more stable than structure CS-I. That could be connected with the larger large-to-small cavities ratio in the CS-I structure that was more suitable for CO 2 molecules. Even at 20% CO 2 in the gas phase, the CO 2 fraction in hydrate reached 0.75 at T = 272 K and 0.70 at T = 274 K. As one can see from Figure 6b, experimental data were described with a reasonable accuracy, at least for relatively low temperatures.
The absence of phase transition CS-I-CS-II in our results at T = 274 K could be connected with the roughness of the used approximation, in which we did not take into account the solubility of gases on water. It is not very important at low pressure, but is significant at high pressure, when N 2 and CO 2 solubility rapidly increase. On the other hand, one can neglect the solubility of gases on the ice at formation conditions. Another system we considered was CO 2 + CH 4 +N 2 mixtures, which can form hydrates with water or ice. In the gas mixture, the N 2 mole fraction was fixed at 0.7 and the relative content of CO 2 and CH 4 varied from 0.00 CO 2 to 0.30 CO 2 in the gas phase. At low concentration of CO 2 in gas phase pressures of methane replacement by CO 2 and N 2 are high and it is not suitable for technology. In this case large cages of hydrate are can be partly filled by nitrogen, but CH 4 recovery degree would be smaller.
In Figure 7a,b, the calculated formation pressures and mole fractions of guests in hydrates of ternary mixtures CO 2 + CH 4 + N 2 are presented.
is not suitable for technology. In this case large cages of hydrate are can be partly filled by nitrogen, but CH4 recovery degree would be smaller.
In Figure 7a,b, the calculated formation pressures and mole fractions of guests in hydrates of ternary mixtures CO2 + CH4 + N2 are presented. The mole fraction of CO2 in hydrate increased rapidly with the growth of its content in the gas phase. The fraction of CO2 at equilibria both with water and with ice, even at small CO2 concentrations, was at least three times higher in the hydrate phase than in the gas mixture ( Figure  7a,b). The behavior of N2 and CH4 guests at equilibria with ice and water were quite different. At equilibrium with ice, CH4 occupied many more cavities than N2 up to a 0.15 mole fraction of CH4 in the gas phase. At higher temperatures, at equilibrium with water and thus for higher formation pressures, the nitrogen became more suitable for occupation of hydrate small cavities.
The solubility of carbon dioxide in water was comparatively high and reached one mole/liter at the hydrate formation pressure and temperature (1.24 MPa at 273 K). At these conditions, the methane solubility was not higher than 0.03 mole/liter. Thus, we could conclude that even at small CO2 concentrations in the gas phase, it could be the excess in the reaction mixture that promotes the methane displacement from hydrate.
In Table 2, the calculated P-T-x equilibria conditions of gas-hydrate-ice (water) systems at 70% in the gas phase are presented.  The mole fraction of CO 2 in hydrate increased rapidly with the growth of its content in the gas phase. The fraction of CO 2 at equilibria both with water and with ice, even at small CO 2 concentrations, was at least three times higher in the hydrate phase than in the gas mixture (Figure 7a,b). The behavior of N 2 and CH 4 guests at equilibria with ice and water were quite different. At equilibrium with ice, CH 4 occupied many more cavities than N 2 up to a 0.15 mole fraction of CH 4 in the gas phase. At higher temperatures, at equilibrium with water and thus for higher formation pressures, the nitrogen became more suitable for occupation of hydrate small cavities.
The solubility of carbon dioxide in water was comparatively high and reached one mole/liter at the hydrate formation pressure and temperature (1.24 MPa at 273 K). At these conditions, the methane solubility was not higher than 0.03 mole/liter. Thus, we could conclude that even at small CO 2 concentrations in the gas phase, it could be the excess in the reaction mixture that promotes the methane displacement from hydrate.
In Table 2, the calculated P-T-x equilibria conditions of gas-hydrate-ice (water) systems at 70% in the gas phase are presented. Table 2. P-T-x equilibria for CO 2 + CH 4 + N 2 +H 2 O in gas-hydrate-ice (water) systems at 70% N 2 in the gas phase.  is the mole fraction of CH 4 in the gas phase; g S 273 h is a type of hydrate structure.
If we assume that the formation pressure rose rapidly with the temperature and with lowering CO 2 content in the gas phase, the obtained data showed very significant increases in the formation pressure after transition from ice to liquid water.

Conclusions
In this work, the method based on the solid solution theory of clathrate hydrate [34][35][36][37][38] was presented to investigate the effects of the influence of nitrogen on the equilibrium pressure and on the hydrate composition of clathrate hydrates formed from methane + carbon dioxide and nitrogen + carbon dioxide binary gas and methane + carbon dioxide + nitrogen ternary gas mixtures. The comparison of cavities' filling degrees for binary CO 2 + CH 4 with ternary CO 2 + CH 4 + N 2 mixed hydrates showed that carbon dioxide molecules occupied both small and large cavities more preferably than methane, but in the case of ternary mixed hydrate, cavities occupation by N 2 , molecules could be concurred by CO 2 and CH 4 molecules.
In all considered systems with and without N 2 , at high and medium CO 2 content in the gas phase, we found that CO 2 was more favorable in occupying clathrate hydrate cavities than CH 4 or N 2 . The addition of N 2 into the gas phase increased the ratio concentration of CO 2 in comparison with the concentration of CH 4 in clathrate hydrates and made gas replacement more effective. The calculation results of the CO 2 + CH 4 hydrates confirmed that for all gas phase compositions, the fraction of CO 2 was higher in the hydrate phase relative to the gas phase, and the hydrate structure CS-I appeared to be more stable than structure CS-II. The mole fraction of CO 2 in the CO 2 + CH 4 + N 2 gas hydrate increased rapidly with the growth of its content in the gas, and the formation pressure of the CO 2 + CH 4 + N 2 gas hydrate rose in comparison with the formation pressure of the CO 2 + CH 4 gas hydrate. Our calculated data were compared with the experimental data [6,25,47], and it was shown that the used theory generally over predicted the experimental data.