Small Alcohols as Surfactants and Hydrate Promotors

: Many methods to produce hydrate reservoirs have been proposed in the last three decades. Thermal stimulation and injection of thermodynamic hydrate inhibitors are just two examples of methods which have seen reduced attention due to their high cost. However, different methods for producing hydrates are not evaluated thermodynamically prior to planning expensive experiments or pilot tests. This can be due to lack of a thermodynamic toolbox for the purpose. Another challenge is the lack of focus on the limitations of the hydrate phase transition itself. The interface between hydrate and liquid water is a kinetic bottle neck. Reducing pressure does not address this problem. An injection of CO 2 will lead to the formation of a new CO 2 hydrate. This hydrate formation is an efﬁcient heat source for dissociating hydrate since heating breaks the hydrogen bonds, directly addressing the problem of nano scale kinetic limitation. Adding limited amounts of N 2 increases the permeability of the injection gas. The addition of surfactant increases gas/water interface dynamics and promotes heterogeneous hydrate formation. In this work we demonstrate a residual thermodynamic scheme that allows thermodynamic analysis of different routes for hydrate formation and dissociation. We demonstrate that 20 moles per N 2 added to the CO 2 is thermodynamically feasible for generating a new hydrate into the pores. When N 2 is added, the available hydrate formation enthalpy is reduced as compared to pure CO 2 , but is still considered sufﬁcient. Up to 3 mole percent ethanol in the free pore water is also thermodynamically feasible. The addition of alcohol will not greatly disturb the ability to form new hydrate from the injection gas. Homogeneous hydrate formation from dissolved CH 4 and/or CO 2 is limited in amount and not important. However, the hydrate stability limits related to concentration of hydrate former in surrounding water are important. Mineral surfaces can act as hydrate promotors through direct adsorption, or adsorption in water that is structured by mineral surface charges. These aspects will be quantiﬁed in a follow-up paper, along with kinetic modelling based on thermodynamic modelling in this work.


Introduction
Natural gas hydrates are classes of composite structures in which water organizes to create cavities that enclathrate small non-polar molecules such as CH 4 , C 2 H 6 , C 3 H 8 , and i-C 4 H 10 .Some small slightly polar components, such as H 2 S and CO 2 , also form hydrates.The molecules that enter cavities are called guest molecules.Hydrates in nature are created from two sources.Methane released through biogenic degradation of organic material in the upper crust is the most abundant source of known hydrates worldwide.These are almost pure methane hydrates and structure I hydrates.The smallest symmetrical unit of structure I hydrate is a cubic box containing 46 water molecules, 6 large cavities (24 water molecules), and 2 small cavities (20 water molecules) that can host molecules such as CH 4 .Molecules that enter these cavities of hydrogen-bounded water molecules are called guest molecules.The size of the cubic unit cell varies with temperature.For typical temperatures Fluids 2021, 6, 345 2 of 30 above the water freezing point, a constant unit box length of 12.01 Å is generally accurate enough.
Many biogenic hydrate resources are connected to deeper sources of thermogenic hydrocarbons.These hydrocarbons are the result of long-term thermal degradation of organic matter and have a large variation in hydrocarbon composition.These resources typically contain sour gases such as H 2 S, which is formed from the degradation of organic sulfur in the absence of oxygen.Utilization of the upper hydrates may lead to an inflow of thermogenic hydrate formers through the fracture systems below.
Structure II is different from structure I in two important ways.The largest cavity has 28 water molecules and can accommodate molecules such as C 3 H 8 and i-C 4 H 10 .The small cavity in structure II is very similar to the small cavity of structure I.The second difference in structure II, as compared to structure I, is that the ratio of small cavities to large cavities in structure II is 2:1.
We will not discuss further details on the basic structures and properties of these two hydrate structures as this information is available in several books including Sloan and Koh [1] and Mokogon [2].These books also contain more information on structure H, which can clathrate guest molecules up to heptane, but is exotic in terms of natural occurrence.
In this work we limit ourselves to structure I hydrates and pure CH 4 hydrates as the energy sources.All the methods and considerations in this work also apply to hydrates from thermogenic hydrocarbons and mixed sources of biogenic and thermogenic hydrates.This work is part of a larger research project that investigates combined safe long-term storage and energy production of CO 2 .More specifically, we examine the addition of limited amounts of N 2 (roughly 20 mole%) to increase the permeability of injections of CO 2 /N 2 mixtures into CH 4 hydrate-filled sediments.As a second additive to the CO 2 /N 2 mixtures, a variety of different surfactants are evaluated experimentally and theoretically.Theoretical evaluations are mainly molecular dynamics (MD) simulations of model systems [3] used to investigate how surfactants affect the interface between liquid water and a separate phase containing hydrate forming components.This separate phase is denoted as gas later in the discussion, although it may be a liquid state hydrate in a former phase.
The primary purpose of the surfactant is to keep the liquid water/gas interface free of blocking hydrate films.Its secondary purpose is to increase the kinetics of mass transport into the liquid water of the interface.Finally, it is expected that the presence of a surfactant will increase the amount of hydrate formerly in the liquid water of the interface.
Small alcohols like methanol and ethanol have been used for many decades as thermodynamic inhibitors.Most experimental data for hydrate formation from water containing alcohols are in a concentration range in which the alcohols change the activity of water significantly.The low partial charge, relative to the size of the methyl group, results in a surfactant effect of methanol.Ethanol has greater surfactant properties due to its non-polar outer methyl group.These surfactant effects can be utilized in hydrate production using the CO 2 /CH 4 hydrate exchange method for the simultaneous release of CH 4 from the hydrate and for the safe, long-term storage of CO 2 as a hydrate.
There is, however, a symbiosis between this project and some industrial hydrate problems caused by methanol.Methanol used to keep wells free of hydrate may result in various amounts of methanol remaining during multiphase transport and/or subsequent processing.This may result in situations in which remaining methanol promotes hydrate formation.A typical example was reported by STATOIL (now EQUINOR) several years ago [4].
The main objective of this work is to examine how limited amounts of ethanol affect water.Limited amounts here refer to amounts that are typically below the amounts needed for hydrate prevention.We facilitate the surfactant properties of ethanol without reducing hydrate formation ability substantially.The focus here is how the addition of small amounts of alcohols affects thermodynamic properties responsible for phase transitions (Gibbs free energy) and associated changes in enthalpies.From a technical point of view, we seek Fluids 2021, 6, 345 3 of 30 alcohols that can easily be mechanically mixed in with the injection gas (CO 2 containing N 2 ).This excludes the high viscosity glycols and the main focus here is on methanol and ethanol.
A secondary objective is to shed light on some important factors related to design of a combined scheme for safe long-term storage of carbon dioxide, and simultaneous release of natural gas for energy.This also includes aspects related to the non-equilibrium nature of hydrates in sediments and different routes to hydrate formation.There are many experimental papers on carbon dioxide hydrate/methane hydrate swapping.It is beyond the scope of this work to discuss any of them.This also excludes our own experiments over a period of almost two decades.
Hydrate phase transition kinetics are implicit functions of thermodynamic control (Gibbs free energy changes), associated mass transport, and associated heat transport.The kinetic models will be derived and discussed in a separate follow-up paper.However, it is important in the context of this work to describe how various thermodynamic properties can be calculated.
A third objective is to demonstrate how the thermodynamic properties related to hydrate phase transitions can be calculated in a consistent way.It is critical in multiphase systems to utilize a uniform reference system.This will facilitate direct comparison of phase stability in terms of Gibbs free energies of the different co-existing phases.
To our knowledge, we are the only group that utilize residual thermodynamics for all phases, including hydrates.For that reason, there are few references to publications from other groups.There are certainly many good publications from other groups.However, within the limited theoretical focus in this work they may not fit into the main objectives discussed above.
Many of the calculated data provided in this work, like for instance Gibbs free energies, are not directly available experimentally.However, model calculations can be verified indirectly towards experimental data, for instance hydrate stability limits in pressure and temperature.Enthalpies of hydrate formation are hard to measure accurately.Additionally, I have not found any experimental data for hydrate formation from dissolved hydrate formers in water [5].To my knowledge, the method utilized in this work for calculating enthalpies [5][6][7][8] is the only general method that can be applied to any hydrate phase transition.Methods that utilize gradient along phase co-existence curves, like Clausius and Clausius-Clapeyron, are at best applicable to pure hydrate formers along pressure temperature stability limit projection.Empirical concepts that utilize a definition of hydrate fugacity are not thermodynamically consistent.Fugacity is defined on a component basis and related to a specific, component based, reference state.In residual thermodynamic formulation fugacity at reference state is pure components pressure.For symmetric excess thermodynamics the reference state is pure liquid component fugacity.Additionally, for asymmetric excess thermodynamics the reference state is infinite dilution fugacity for the component in a specific solvent.
The concept we present, and use, in this work is thermodynamically consistent.The reason for the consistency is that the equation for enthalpy is derived from Gibbs free energy using fundamental thermodynamic relationships.Practically this also implies that the corresponding entropy change for hydrate formation is consistent.It can therefore be expected that the final hydrate has the correct hydrate structure.
No similar studies on water containing ethanol have been found in the accessible literature.
The paper is organized as follows.The traditional use of alcohols for thermodynamic hydrate inhibitor is briefly discussed in Section 2. Since this is not an important focus in this work, there are no references to any publications from my research group or other international research groups.In Section 3, I give a brief overview of hydrate production philosophy.The main focus of this section is on and kinetic challenges related to the various hydrate production methods.The use of CO 2 for exchange with in situ CH 4 in hydrate is a priority in that section.Scientific methods utilized in this work are briefly described in Sec-tion 4. The non-equilibrium thermodynamic nature of hydrates in sediments is discussed in Section 5.In particular it is argued that hydrates in sediments cannot reach thermodynamic equilibrium because there are too many independent thermodynamic variables in the co-existing phases compared to constraints on these variables.These constraints are mass conservations and thermodynamic equilibrium equations.Hydrates in sediments are normally in a stationary state.The rest of the paper is devoted to thermodynamic analysis of the consequences of adding small amounts of ethanol to CO 2 and mixtures of CO 2 with N 2 .The purpose of the analysis is to investigate if the ability to form a new hydrate from CO 2 or CO 2 /N 2 is significantly reduced by adding up to 3 mole percent ethanol to the pore water.The effects of adding ethanol on heterogeneous gas/liquid formation are discussed in Section 6.The main focus in this section is on hydrate stability.A secondary focus in this section is the ability of injection gas (CO 2 or CO 2 /N 2 mixture) to form a new hydrate which releases enough heat to dissociate in situ CH 4 hydrate.A similar analysis for homogenous hydrate formation from water and dissolved hydrate former is given in Section 7. The paper is concluded with a discussion in Section 8 and conclusions in Section 9.

The Use of Alcohols for Hydrate Prevention in Industrial Settings
Problems related to the formation of hydrate from hydrocarbon mixtures and water have been important motivations for industrial hydrate research during the latest seven decades.The interest for natural gas hydrates as energy source has increased substantially during the latest four decades.Research activities on hydrates in sediments, and hydrate production technologies, have increased substantially.
Historically, small alcohols like methanol and ethanol have been used to prevent hydrate formation.Addition of alcohols to liquid water results in a reduction (more negative) of liquid water chemical potential.A consequence of this reduced liquid water chemical potential is that a higher pressure is needed in order to form hydrate.This will be discussed in more detail later.
Methanol totally dominates the market among these two alcohols.Exceptions are countries, such as Brazil, which produce ethanol at reasonable prices from a waste product in sugar production.
Thermodynamic hydrate inhibition as discussed above is based on the effect of dissolved alcohols on "bulk" water [9,10].Most efficient concentrations of alcohol are concentrations when alcohol is the solvent for water, rather than the opposite.The transition from water as solvent, over to methanol as solvent, is clearly visible as a transition in the liquid mixture dielectric constant as function of alcohol concentration [10].
There are many different ways that a hydrate can form in a multiphase system.Hydrate formation on the interface between liquid water, and a separate phase for guest molecules.Guest molecules are the hydrate forming molecules that enter the cavities created by hydrogen bonded water.Methane and natural gas are examples of phases containing guest molecules.Hydrate can only form if Gibbs free energy of the hydrate is lower than Gibbs free energy for the same amounts of water and guest molecules in the original phases.Technically a hydrate equilibrium calculation involves solving for the same water chemical potential in liquid water and hydrate water when also chemical potentials for guest molecules are the same in hydrate and in gas phase.Addition of significant amounts of alcohols to the water reduces the chemical potential of liquid water and higher pressure is needed in order to create hydrate.
However, hydrate can also form from dissolved hydrate formers in liquid water, from water dissolved in gas and towards mineral surfaces like for instance rust.As will be discussed here, all these hydrates are different in composition and stability (Gibbs free energy).
Hydrate inhibition is not a primary objective in this work.However, many of the thermodynamic calculations presented here can also provide valuable extensions in hydrate risk analysis.In particular, the multi-phase hydrate stability analysis should be incorporated.An analysis based on temperature and pressure alone is not sufficient.There are many different ways that hydrate can dissociate again even in temperature and pressure is inside hydrate formation region.
Since there are several routes to hydrate formation it is important to be able to analyze the relative stabilities of the different hydrates than are formed in presence of alcohols if the local inhibitor concentrations are insufficient for total thermodynamic inhibition of hydrate.

Small Alcohols as Surfactants in Hydrate Production
Countries such as China and Japan are very actively working towards full scale hydrate production schemes.
A kinetically critical element in hydrate dissociation is the transport of guest molecules across a thin (roughly 1.2 nm) interface between a hydrate and its surroundings [8,9,[11][12][13].
During the hydrate dissociation process, there will be liquid water facing the dissociating hydrate.The interface between liquid water and the hydrate will reproduce itself continuously.This is determined by the molecular physics of water.Partial charges on the surface of a hydrate are fairly fixed except for small translational and rotational movements relative to a minimum energy situation.Liquid water molecules, on the other hand, are highly dynamic.However, liquid water molecules have to relate to the hydrate water molecules in order to optimize entropy.
Methanol has surfactant properties because of the large methyl group.The low partial charge on the methyl group results in a small charge per atomistic group surface.The concentration of methanol, on the surface of water that faces the gas, will therefore be higher than the average concentration of methanol in "bulk" water.Carbon dioxide has a significant quadrupole moment.Solubility of carbon dioxide in water is therefore higher than solubility of small hydrocarbons in water.
Efficient hydrate dissociation depends on a favorable Gibbs free energy change.Temperature and pressure are two of all the independent thermodynamic variables.Concentrations of all molecules in all the co-existing are the other independent thermodynamic variables.If temperature and pressure are outside hydrate formation limits then gas and liquid water is more stable than the hydrate and Gibbs free energy change is favorable for hydrate dissociation.This is, however, only one of the thermodynamic conditions that has to be fulfilled.The first law of thermodynamics also requires that enough the enthalpy of hydrate dissociation can be supplied.
Reducing pressure to a condition outside hydrate stability is considered as a low-cost method for producing natural gas from a hydrate.The challenge is that the natural supply of heat may not be sufficient.Pressure reduction results in a temperature gradient that can supply heat.Local geological depth involves a geothermal gradient.Both of these possible sources of heat are low temperature sources and may not be efficient for breaking hydrogen bonds.The duration of early stages pilot tests [14,15] were extremely short.The results from these pilot studies are far too short to provide significant learning about hydrate production.The first offshore test [16] ended after 6 days.Officially, the reasons were massive sand production and reservoir freezing.The second test [17] was supposed to last for 6 months but the sediments froze after 24 days [17][18][19].
The enthalpy needed to dissociate CH 4 hydrate is roughly 58 kJ/mole CH 4 at 273.15 K [5][6][7][8][9].Heat is transported fast through water phases [20].Injection of steam or hot water is expensive.Heat will also be lost to minerals.However, it still addresses the breaking of hydrogen bonds.Injection of alcohols or salts is also very efficient in breaking hydrogen bonds.
Natural gas hydrates in sediments are unable to touch mineral surfaces due to extremely low chemical potential of the structured adsorbed water [21][22][23][24][25][26].Another way to look at this fact is to consider distribution of atomistic charges in mineral surfaces.Distribution of atomistic charges on mineral surfaces are noy compatible with average partial charges in liquid water.The result is rather extreme water densities in the first adsorbed water layers on mineral surfaces.Many theoretical studies are available in the literature.A more limited number of experimental studies are also available.Water densities in the order of 3 times liquid water density is quite common.The paper by Geissbühler et al. [27] is a representative experimental report.Geo-mechanical models based on hydrate "cementing" pores are therefore physically wrong.Practically there will be minimum distances between mineral surfaces and hydrate.This minimum distance is controlled by the mineral/water structures discussed above, and also hydrate/liquid water interfaces.Molecular diffusion and induced hydrodynamic flow increase the distance between mineral surfaces and hydrate.
Fracture systems lead to hydrodynamic flow in offshore hydrates.Fractures that connect the seafloor to hydrate filled sediments lead to hydrate dissociation.The reason for this dissociation is the seawater concentration of guest molecules, which are almost at infinite dilution in the seawater.Limits for hydrate stability, in terms of guest concentration in surrounding liquid water, are discussed in the section describing homogeneous hydrate formation.Fracture systems below the hydrate can lead to inflow of new hydrocarbons, and can lead to subsequent hydrate formation when the hydrocarbons enter hydrate forming conditions.Hydrates in sediments are not thermodynamically stable (see Section 5).The system of a hydrate, liquid water, and gas will establish a situation that can be close to stationary.The dynamics of this stationary situation varies from one hydrate reservoir to another.The stationary situation is not constant.A net hydrate dissociation, which is higher than the creation of new hydrate formation, will over time lead to reduced hydrate saturation.Even in the permafrost in Alaska the average hydrate saturation is around 75% and it is rare to find hydrate deposits that exceed 85% of hydrate pore filling volume.
Injection of carbon dioxide is efficient because there is available free pore water that can create carbon dioxide hydrate.Enthalpy of hydrate formation for CO 2 hydrate is roughly −68 kJ/mole CO 2 at 273.15 K 5-9 .Addition of up to 30 mole% N 2 to the CO 2 increases this value to around −66 kJ/mole guest.The net difference between what is needed to dissociate the in situ CH 4 hydrate is still sufficient.
The first question is whether addition of limited amounts of alcohol will significantly disturb the hydrate formation ability (Gibbs free energy).The second question is whether the released enthalpy from formation of new hydrate from CO 2 or CO 2 /N 2 mixture is sufficient to dissociate the in situ CH 4 hydrate.

Methodology
The primary scientific tool in this work is classical thermodynamics.There are elements of statistical mechanics in the historical development of the theory for chemical potential of water in hydrate, and interactions between water and guest in the hydrate lattice.As will be discussed in Section 6, our model for chemical potential of water in hydrate differs from older models.Using molecular dynamics simulations, we were able to incorporate effects of guest movements on water lattice destabilization.For example, carbon dioxide movements interfere with water movements in the hydrate lattice.The interference affects specific water movement frequencies and leads to water lattice destabilization.The result is a difference of roughly 1 kJ/mole between a rigid hydrate water lattice and a dynamic lattice.Ethane in a large cavity also disturbs water movements but less so than carbon dioxide.For methane as a guest, the difference between a dynamic water lattice and a rigid lattice is almost zero.

Non-Equilibrium Nature of Hydrates in Sediments
The early hydrate equilibrium experiments started around 1940.With CH 4 gas, liquid water, and a hydrate phase the researchers knew from Gibbs phase rule that only one independent thermodynamic variable could be fixed.Despite this fact, these hydrate equilibrium curves are also used as equilibrium curves when both temperature and pressure are fixed.As will be discussed later, hydrate stability also depends on a minimum concentration of guest molecules in the surrounding liquid water.Hydrates can form from liquid water solutions of guest molecules.The concentration range for this hydrate formation is between liquid water solubility concentration and minimum concentration for hydrate stability.Hydrates formed along this route have different composition and density.They are therefore thermodynamic definition.
Hydrates formed from different phases are different because composition, density, and Gibbs free energy are different.The reason for this is that there is a mathematical imbalance between the number of independent thermodynamic variables and the associated constraint (mass conservation and equilibrium equations).If the thermal (same temperatures in all phases) equilibrium and the mechanical equilibrium (same pressures in all phases) is fulfilled, we can return to the simple system of CH 4 , liquid water, and a hydrate.There are six remaining mole-fractions in the three phases.The constraints are the conservation of mole-fractions in each of the three phases.In addition, there are 4 independent equations of equal chemical potential for CH 4 and water in the three phases.With 7 equations in 6 variables there is no unique mathematical solution.The first and second laws of thermodynamics reduce to the following form when temperatures and pressures are the same for all phases: in which the line below G denotes extensive property (unit Joule), j is a phase index, and N is the number of moles in each phase, as indicated by counter j.Total number of co-existing phases is p. i is a component index and x is mole-fraction.The line above x in the chemical potential dependency denotes the vector of all mole-fractions in the specific phase.This also includes mineral surfaces since they serve as catalysts for hydrate nucleation while at the same time not being able to "host" the hydrate on the mineral surface.It also includes all different hydrate phases and in principle also adsorbed phases on hydrate surfaces.The abbreviation min denotes minimum.
Number of degrees of freedom is the number of thermodynamic independent variables that must be specified for a system to be able to reach equilibrium.I prefer to use discrete counting rather than the compact Gibbs phase rule formula.One reason for this is that mixtures of guest molecules in a gas phase increases the complexity.The result can be that the system is brought even more out of the possibility to read the thermodynamic equilibrium.Many of these aspects will not be captured by Gibbs phase rule.As an example, consider an equimolar mixture of carbon dioxide and methane.Carbon dioxide is closer to condensation on water than methane.Carbon dioxide has a stronger attraction to the liquid water surface.In more rigorous physics, one example of a two-dimensional adsorption theory is described by Kvamme [37].If the system is closed, then the result of selective adsorption is that many different hydrates can form.Carbon dioxide will have a higher mole-fraction in the adsorbed layer on liquid water.This also implies a higher mole-fraction than methane of the liquid side of the interface.The first hydrates that form will be dominated by carbon dioxide.Gradually the carbon dioxide content of the gas decreases and the final hydrates that form will be pure methane hydrate.Similar physical differences apply to other ways that hydrate can form, such as hydrate formation towards solid surfaces.
In summary, there is no way to establish thermodynamic equilibrium for hydrate in sediments.A theoretical philosophy about phases that can be totally consumed and then disappear from the balance will generally not be realistic for real hydrate reservoirs.As mentioned above, offshore hydrates are normally in a stationary balance of hydrate dissociation and formation of new hydrates from upcoming gas.On top of this, there are active phases that are normally not counted for, including mineral adsorption and hydrate nucleation towards mineral surfaces In view of the above, the thermodynamic equilibrium is not generally possible for hydrates in sediments.A similar line of arguments can be established for hydrates in industrial processing units and in pipeline transport of gas containing dissolved water.Multiphase transport of gas and water is another example.
Due to general non-equilibrium situation, I will use the term stability limits rather than equilibrium curves.For heterogeneous hydrate formation on gas/liquid water interface the set natural set of independent thermodynamic variables are temperature and pressure.For homogeneous hydrate formation from liquid water and dissolved gas components, the independent variables are mole-fractions of the guest molecules in water.

Heterogeneous Hydrate Formation in Systems Containing Alcohols
20 mole% N 2 seems feasible based on earlier studies in open literature.In order to limit the number of figures, I investigate hydrate free energy for pure CH 4 , pure CO 2 , and an injection mixture with 20 mole% N 2 in CO 2 .Most natural gas hydrates are formed from biogenic sources and almost pure CH 4 hydrates.There are, however, no limits in the analysis here.The same calculations can be applied to thermogenic gas and mixed structure I and structure II hydrates.
From an environmental point of view, and also in terms of surfactant properties, ethanol is more interesting than methanol.I will therefore mainly focus on ethanol.In Section 6.1, I focus on the free energy changes needed to form hydrates and in Section 6.2 I examine changes in released enthalpy during hydrate formation.
The main goal is to establish knowledge on the stability of CH 4 and CO 2 hydrates in real thermodynamic variables (Gibbs free energy) and the ability of hydrate formation from injection gas (CO 2 and CO 2 /N 2 mixtures) to deliver enough heat to dissociate in situ CH 4 hydrate.
The thermodynamic models are described in Section 6.1 after some examples for model verification purposes.The model for water activity coefficients used in Equation ( 2) is given in Table 1.Note that this is a simplified model.For NaCl it only applies to calculations of water activity.For situations in which salinity in the water is important for interactions with mineral reactions then a totally different model is needed.

Hydrate Phase Transition and Free Energy Changes
Ethanol and methanol only have one hydroxyl group.It can be of interest to compare hydrates formed from liquid water containing glycols and CH 4 with hydrates formed from water containing ethanol or methanol.MEG is mono-ethylene glycol and TEG is tri-ethylene glycol.In Figure 1a, I plot hydrate stability pressure limits as function of temperature and concentration of MEG or TEG in the water.Gibbs free energies for the formed hydrates are plotted in Figure 1b.
Similar plots for water containing ethanol and methanol are plotted in Figure 2a,b.

Hydrate Phase Transition and Free Energy Changes
Ethanol and methanol only have one hydroxyl group.It can be of interest to compare hydrates formed from liquid water containing glycols and CH4 with hydrates formed from water containing ethanol or methanol.MEG is mono-ethylene glycol and TEG is triethylene glycol.In Figure 1a, I plot hydrate stability pressure limits as function of temperature and concentration of MEG or TEG in the water.Gibbs free energies for the formed hydrates are plotted in Figure 1b.
Similar plots for water containing ethanol and methanol are plotted in Figure 2a  Stability of hydrates formed from water containing glycols, and hydrates formed from water containing ethanol and methanol are very similar.Fugacity formulations dominate hydrate stability limit calculations in open literature on hydrates.Academic and commercial hydrate software is also based on the so-called reference scheme from around 1970.It might be useful to provide the residual thermodynamic model behind the curves in Figures 1 and 2. Liquid water chemical potential is  [40] and solid curves are calculated) and methanol (squares are experimental data [41] and dashed curves are calculated).Blue circles and solid blue curves are for 15 wt% ethanol.Green squares and dashed curve are for 10 wt% methanol.Dashed blue is for 15 wt% methanol.Red squares and dashed are for 20 wt% methanol.Solid black curve is calculated data for pure water.Crosses are experimental data from Sabil et al. [42] and plusses are experimental data from Tumba et al. [43].(b) Gibbs free energy for hydrates formed from CH 4 and water containing various concentrations of ethanol and methanol.Same notations as in Figure 2a.
Stability of hydrates formed from water containing glycols, and hydrates formed from water containing ethanol and methanol are very similar.Fugacity formulations dominate hydrate stability limit calculations in open literature on hydrates.Academic and commercial hydrate software is also based on the so-called reference scheme from around 1970.It might be useful to provide the residual thermodynamic model behind the curves in Figures 1 and 2. Liquid water chemical potential is formulated by symmetric excess thermodynamics as: in which p is an indicator of liquid phase and p is water in this case.Subscript i denotes component index and i is water in this case.Superscript mix denote mixture and superscript pure denotes pure component.x is liquid mole-fraction and the arrow on x denote vector of all mole-fractions in the liquid mixture.γ (p) i is the activity coefficient for component i. Asymptotic limit is unity when mole-fraction of component i approaches unity.Activity correction is in excess thermodynamics.The whole quantity is, however, in residual thermodynamics since the pure water term (first term on right hand side) is derived from Molecular Dynamics (MD) simulations for ice.Consistent chemical potential for pure liquid water is derived from enthalpy of ice dissociation at 0 • C and extended in the liquid region using experimental heat capacities.The second term on the right-hand side is the ideal mixing term due to entropy of mixing.The final term corrects for the deviations from ideal mixing.A simple model for water activity coefficients in mixtures with some alcohols is given in Table 1.These are based on 1 bar pressure.Correction from 1 bar to actual pressure is included through a Poynting correction for pure liquid water in Equation (2).A Poynting correction for water in empty hydrate is also included in Equation (3) below.
Chemical potential for water in hydrate is given by: is the chemical potential for water in an empty clathrate.Number of cavities is ν, with subscripts k for large and small cavities, respectively.For structure I, which is the main focus in this overview on thermodynamics, ν large = 3/24 and ν small = 1/24.Within the scope of this work, I will assume that only one guest molecule can enter a cavity.The harmonic oscillator approach model can then be expressed as: Chemical potential for molecule type i in cavity type k.I will assume that small and large cavities are at equilibrium so that: For a system at equilibrium the chemical potential for a guest molecule in a cavity is equal to the chemical potential for the same molecule in the equilibrium phase.∆g ki is the free energy change for inclusion for guest molecule i in a cavity of type k.
The most classical example is a hydrate former phase (gas, liquid, or supercritical) in contact with liquid (or ice) water which forms a hydrate at the interface.For these three phases there are 12 independent thermodynamic variables, 3 mass conservation equations, and 8 conditions of equilibrium.The difference between the number of independent variables minus mass conservations and equilibrium conditions leaves only one independent thermodynamic variable that can be fixed for equilibrium to be achieved.For this particular case of equilibrium Equation ( 4) can be written as: For the heterogeneous hydrate formation, the phase index p in Equation ( 7) is gas.
is the fugacity coefficient for CH 4 .In this work I utilize the Soave-Redlich-Kwong (SRK) equation of state [14].The first term on the right-hand side is the ideal gas chemical potential as function of temperature and density.Ideal gas chemical potential is trivial to calculate using the momentum space of the canonical partition function.The methane model is a sphere and the methane ideal gas chemical potential is easy to calculate.The analytical expression can be found in any textbook on physical chemistry.Guest molecules with rotational degrees of freedom also require the moments of inertia for the rotational degrees of freedom ideal gas chemical potential.In this work, all models are rigid and there are no ideal gas contributions from intramolecular degrees of freedom.The second term is the ideal gas mixing terms that comes from the entropy of mixing ideal gas components at constant pressure.
Since methane is pure in the first systems illustrated in Figures 1 and 2 mole-fraction is unity.Equation ( 7) enters into Equation ( 6).The free energy change of inclusion for methane in small and large cavity is given elsewhere [9,13,26,[44][45][46].With Equation (5) in Equation ( 2) and solving Equation ( 2) equal to Equation ( 3) for the same chemical potential of water in liquid water as in the hydrate.Chemical potential for water in empty clathrate is also given by Kvamme and Tanaka [46].I have solved the equations iteratively for pressures, with temperature as the chosen fixed variable.The chemical potential for the gas changes significantly with pressure.I therefore also plot the composition of the formed hydrate as well as chemical potential of CH 4 .
The filling fractions are trivially available from the semi grand canonical ensemble used for derivation of (3).This part is the same for the van der Waal and Platteeuw [47] model and our model [46]: θ ki is the filling fraction of component i in cavity type k.Also: x H i,small = θ small,i ν small 1 + θ large,i ν large + θ small,i ν small (10) where ν is the fraction of cavity per water for the actual cavity type, as indicated by subscripts.The corresponding mole-fraction water is then given by: and the associated hydrate free energy is then: Additionally, note that the results in Figures 1a and 2a are verification of the model system since these figures are solutions to: The molecular weight of ethanol and methanol is very different.The weight% alcohol used in industry is not convenient for comparison of thermodynamic effects since these are on a molar basis.
The changes in chemical potential for CH 4 as guest varies substantially along the stability limits of temperature and pressure, and also as function of increasing ethanol content.The composition of the hydrates also varies substantially with amount of ethanol added to the water.
Addition of ethanol as a surfactant to CO 2 /N 2 system should rarely exceed 2 mole% but I also include 3 mole% as a possible upper limit and 1.5 mole% as minimum addition of ethanol.CO 2 hydrate is generally a more stable hydrate than CH 4 hydrate.This is illustrated in Figure 3b below.The density of CO 2 is also higher than CH 4 for the same range of conditions.There is a CO 2 phase transition to a higher density that results in a (relative) steep change in hydrate formation pressure as function of temperature.The stability limits of the hydrate are not substantially shifted for ethanol concentrations up to 3 mole% in water.In Figure 4a, we plot mole-fractions CO 2 in the formed hydrates as function of ethanol mole-fractions up to 3 mole%.Based on Figure 3, it is considered feasible to add up to 3 mole% ethanol to the water phase.Thermodynamic stability of the formed hydrate is not affected substantially for ethanol concentrations below this limit.The changes in hydrate composition and chemical   Based on Figure 3, it is considered feasible to add up to 3 mole% ethanol to the water phase.Thermodynamic stability of the formed hydrate is not affected substantially for ethanol concentrations below this limit.The changes in hydrate composition and chemical potential of CO2 in hydrate are also small, as is illustrated in Figure 5. Changes in temper- Based on Figure 3, it is considered feasible to add up to 3 mole% ethanol to the water phase.Thermodynamic stability of the formed hydrate is not affected substantially for ethanol concentrations below this limit.The changes in hydrate composition and chemical potential of CO 2 in hydrate are also small, as is illustrated in Figure 5. Changes in temperature and pressure hydrate formation conditions are also limited.Technically the best way is to add the ethanol to the injection gas.This will likely keep the surfactant effect for longer times before ethanol distributes and dissolve into groundwater.
The same calculations presented in Figures 3 and 4 are also conducted for 20 mole% N 2 in CO 2 /N 2 mixture.Results are presented in Figures 5 and 6.There is no experimental evidence that CO2 enters the small cavity of structure I at liquid water conditions.Hydrate formation mechanisms in the ice range is different and can likely lead to trapping of CO2 in small cavity.There are some references to work conducted in the group of Werner Kuhs [48][49][50] that provides evidence of small cavity filling of CO2 for temperatures far below zero and atmospheric pressures.We have conducted theoretical studies on different levels.These range from quantum mechanics to molecular dynamics simulations in the classical regime.All conclusions are the same.There is no free energy contribution to hydrate stability from CO2 filling in small cavities at liquid water temperatures.At some conditions of full filling of large cavities with CO2, a limited portion of CO2 in small cavities is not enough to result in destabilization of the hydrate as a whole.However, that does not mean that CO2 "willingly" enters a small cavity.According to Equation ( 5), the canonical partition function for CO2 in small cavity is zero at liquid water conditions.
Another challenge in the experiments that are conducted at atmospheric pressure and 1 bar [48][49][50][51] is that the conversion between CH4 hydrate and CO2 hydrate is totally different since there is no liquid side interface.Specifically, a solid-state mechanism has been proposed [51].A solid-state mechanism is not in accordance with kinetic rates we have observed during 2 decades of CH4/CO2 exchange experiments in liquid water regime.As mentioned in the introduction, this is a theoretical paper and we are therefore not reporting from our experiment or any other experiments on this specific issue.
For CO2/N2 mixtures at liquid water state conditions the hydrate formation from injection gas will facilitate from N2 filling of the small cavities.This can be seen from a different stability limit window in Figure 5a.However, the free energies of the formed hydrates are less stable than hydrates created from pure CO2, as can be seen by comparing Figures 3b and 5b.This is due to the change in CO2 chemical potential though dilution with N2.The reduction in stability is, however, limited and stability of the mixed hydrates formed from up to 3 mole% ethanol in water is significantly more stable than the CH4 hydrate.Thermodynamically it is feasible to use ethanol in CO2/N2 with 20 mole% N2, and an amount of ethanol that would result in maximum 3 mole% ethanol in the pore water.Higher mole-fractions of ethanol in water have not been investigated since 3 mole% ethanol in water is equivalent to 7.33 weight%.

Enthalpies of Hydrate Formation
A consistent way to calculate is given by the general thermodynamic rela- There is no experimental evidence that CO 2 enters the small cavity of structure I at liquid water conditions.Hydrate formation mechanisms in the ice range is different and can likely lead to trapping of CO 2 in small cavity.There are some references to work conducted in the group of Werner Kuhs [48][49][50] that provides evidence of small cavity filling of CO 2 for temperatures far below zero and atmospheric pressures.We have conducted theoretical studies on different levels.These range from quantum mechanics to molecular dynamics simulations in the classical regime.All conclusions are the same.There is no free energy contribution to hydrate stability from CO 2 filling in small cavities at liquid water temperatures.At some conditions of full filling of large cavities with CO 2 , a limited portion of CO 2 in small cavities is not enough to result in destabilization of the hydrate as a whole.However, that does not mean that CO 2 "willingly" enters a small cavity.According to Equation ( 5), the canonical partition function for CO 2 in small cavity is zero at liquid water conditions.
Another challenge in the experiments that are conducted at atmospheric pressure and 1 bar [48][49][50][51] is that the conversion between CH 4 hydrate and CO 2 hydrate is totally different since there is no liquid side interface.Specifically, a solid-state mechanism has been proposed [51].A solid-state mechanism is not in accordance with kinetic rates we have observed during 2 decades of CH 4 /CO 2 exchange experiments in liquid water regime.As mentioned in the introduction, this is a theoretical paper and we are therefore not reporting from our experiment or any other experiments on this specific issue.
For CO 2 /N 2 mixtures at liquid water state conditions the hydrate formation from injection gas will facilitate from N 2 filling of the small cavities.This can be seen from a different stability limit window in Figure 5a.However, the free energies of the formed hydrates are less stable than hydrates created from pure CO 2 , as can be seen by comparing Figures 3b and 5b.This is due to the change in CO 2 chemical potential though dilution with N 2 .The reduction in stability is, however, limited and stability of the mixed hydrates formed from up to 3 mole% ethanol in water is significantly more stable than the CH 4 hydrate.
Composition of the hydrate is plotted in Figure 6a.Chemical potentials of the guests in the hydrate are plotted in Figure 6b.The values are reasonable in view of the influence of CO 2 on chemical potential for N 2 through the density of the mixture and the fugacity coefficient.
Thermodynamically it is feasible to use ethanol in CO 2 /N 2 with 20 mole% N 2 , and an amount of ethanol that would result in maximum 3 mole% ethanol in the pore water.Higher mole-fractions of ethanol in water have not been investigated since 3 mole% ethanol in water is equivalent to 7.33 weight%.

Enthalpies of Hydrate Formation
A consistent way to calculate ∆H Total is given by the general thermodynamic relationship: The route to enthalpy changes, from the residual thermodynamic scheme described in this work, is fairly straightforward utilizing Equation (3) for hydrate water.Liquid water chemical potential from Equation (2) as basis for corresponding enthalpy of liquid water is trivial and available in thermodynamic textbooks.
I now use H 1 to denote heterogeneous hydrate formation from liquid water (or ice) and methane gas.The same proportionality as in Equation (14) applies to relationship between chemical potential and partial molar enthalpy: in which ∆H is the partial molar enthalpy change for water during the phase transition.
Differentiation of Equation (3) gives the following results for partial molar enthalpy of water in hydrate [6,[12][13][14][15]: and from Equation (2) for liquid water: and the enthalpy of hydrate formation is then: The last term is the partial molar enthalpy of the guest molecules in the gas (or liquid) hydrate former phase.This is trivially given by the Equation of state using the same type of relationship as Equation ( 15) with chemical potential for methane from Equation (7).The analytical expression for the SRK [52] version of this term is available from engineering handbooks and does not need space here.The partial molar enthalpy of each guest molecule in the hydrate, H H 1 i needs some more attention.The ideal gas enthalpy depends only on temperature and will be the same for the guest molecule inside the cavities and in the gas mixture.As such we only need the residual contributions for these two terms.Entropy related properties like chemical potentials and free energies are sensitive to possible interactions between water movements and guest movements.Average interaction energy is less sensitive to these effects.I have therefore used a Monte Carlo approach [53,54] with fixed water lattice to evaluate average interaction energies between guest and water for the various cavity types in structure I.I also sample the efficient volume for the guest in the cavity in order to calculate a compressibility factor for the guest molecule inside the two types of cavities.
where V ki is the sampled molar volume of guest molecule i in cavity type k. x ) ∂T (21) In Figure 7 below, I plot enthalpies of hydrate formation for pure CO 2 and water containing various concentrations of ethanol up to mole-fraction 0.03.These are plotted as per mole hydrate in (a) as well as per guest in (b) in order to illustrate that the temperature dependency of the coordination number (moles water per mole guest in the formed hydrate) changes significantly when adding N 2 .As an example: for the pure CO 2 hydrate formed form water containing 3 mole% ethanol the coordination number changes from around 7.0 for the lowest temperature range and down to 6.5 after the CO 2 phase transition temperature.The coordination number is actually the numbers in Figure 8b divided by corresponding numbers in Figure 7a.Similar comparison for CO 2 with 20 mole% N 2 is illustrated is Figure 7c,d.As can be seen from Figure 6a, there is a significant change in the nitrogen filling of the hydrate with temperature.The coordination number changes substantially with temperature as compared to the pure CO 2 case.However, as can be seen by comparing Figure 7b,d, the available enthalpy that can be used to dissociate in situ CH 4 hydrate is significantly decreased by the addition of N 2 .
Fluids 2021, 6, x FOR PEER REVIEW 17 of 31 by corresponding numbers in Figure 7a.Similar comparison for CO2 with 20 mole% N2 is illustrated is Figure 7c,d.As can be seen from Figure 6a, there is a significant change in the nitrogen filling of the hydrate with temperature.The coordination number changes substantially with temperature as compared to the pure CO2 case.However, as can be seen by comparing Figure 7b,d, the available enthalpy that can be used to dissociate in situ CH4 hydrate is significantly decreased by the addition of N2.
How much ethanol that will dissolve into "bulk" pore water during the average residence time of the injection gas in the pore is unknown.As such, it makes sense to also compare enthalpies per guest for methane as function of mole% ethanol in water as well.This will enable us to illustrate the range between maximum and minimum available formation enthalpy that can be used for dissociation of the in situ CH4 hydrate.In Figure 8, I plot the same curves as in Figure 7 but now also with CH4 hydrate formation enthalpy included.In the worst case scenario, the available net enthalpy of hydrate formation from the CO2/N2 mixture is roughly 6 kJ/mole more than what is needed for dissociation of the in situ CH4 hydrate.These are all values at pressure temperature phase stability boundaries.Heat transport is at least 2-3 times faster than mass transport through water systems [20].The hydrate formation enthalpy for hydrates formed from CO2/N2 mixtures containing up to 20 mole% N2 is therefore sufficient to dissociate in situ CH4 hydrate for up to 3 mole% ethanol in the pore water.As illustrated in the previous section, the free energy of the hydrates formed from CO2/N2 mixture is thermodynamically more stable than CH4 hydrates.This is illustrated by comparing Figures 2b and 5b.The difference is in the order of 2 kJ/mole hydrate.Formation of a new hydrate from the injection gas (CO2/N2/surfactant) will result in higher salinity of the remaining liquid pore water.Eventually, the resulting changes in liquid water chemical potential will lead to dissociation of hydrate in the pore.Additionally, since the stability of CH4 hydrate is lower than stability of CO2 hydrate then the CH4 hydrate will dissociate first.

Homogeneous Hydrate Formation in Systems Containing Alcohols
CH4 from dissociating hydrates is expected to distribute as gas bubbles escaping upwards in the reservoir due to buoyancy.Some methane will dissolve rapidly in water and can generate hydrate cores for further growth elsewhere.However, it is also important to investigate hydrate stability limits towards concentration of CH4 in surrounding water.This is also a route to hydrate dissociation even for other production schemes.Pumping out water during pressure reduction in a hydrate reservoir that has no free gas leads to change in water surrounding the hydrate.If the surrounding water contains less CH4 than hydrate stability limit concentration then the hydrate will dissociate.
Of course, we need to know the enthalpy of homogeneous hydrate formation since the released enthalpy of hydrate formation will generate local temperature gradients.
Homogeneous CH4 hydrate formation is discussed in Section 7.1 and homogeneous CO2 hydrate formation is discussed in Section 7.2.How much ethanol that will dissolve into "bulk" pore water during the average residence time of the injection gas in the pore is unknown.As such, it makes sense to also compare enthalpies per guest for methane as function of mole% ethanol in water as well.This will enable us to illustrate the range between maximum and minimum available formation enthalpy that can be used for dissociation of the in situ CH 4 hydrate.
In Figure 8, I plot the same curves as in Figure 7 but now also with CH 4 hydrate formation enthalpy included.In the worst case scenario, the available net enthalpy of hydrate formation from the CO 2 /N 2 mixture is roughly 6 kJ/mole more than what is needed for dissociation of the in situ CH 4 hydrate.These are all values at pressure temperature phase stability boundaries.
Heat transport is at least 2-3 times faster than mass transport through water systems [20].The hydrate formation enthalpy for hydrates formed from CO 2 /N 2 mixtures containing up to 20 mole% N 2 is therefore sufficient to dissociate in situ CH 4 hydrate for up to 3 mole% ethanol in the pore water.As illustrated in the previous section, the free energy of the hydrates formed from CO 2 /N 2 mixture is thermodynamically more stable than CH 4 hydrates.This is illustrated by comparing Figures 2b and 5b.The difference is in the order of 2 kJ/mole hydrate.Formation of a new hydrate from the injection gas (CO 2 /N 2 /surfactant) will result in higher salinity of the remaining liquid pore water.Eventually, the resulting changes in liquid water chemical potential will lead to dissociation of hydrate in the pore.Additionally, since the stability of CH 4 hydrate is lower than stability of CO 2 hydrate then the CH 4 hydrate will dissociate first.

Homogeneous Hydrate Formation in Systems Containing Alcohols
CH 4 from dissociating hydrates is expected to distribute as gas bubbles escaping upwards in the reservoir due to buoyancy.Some methane will dissolve rapidly in water and can generate hydrate cores for further growth elsewhere.However, it is also important to investigate hydrate stability limits towards concentration of CH 4 in surrounding water.This is also a route to hydrate dissociation even for other production schemes.Pumping out water during pressure reduction in a hydrate reservoir that has no free gas leads to change in water surrounding the hydrate.If the surrounding water contains less CH 4 than hydrate stability limit concentration then the hydrate will dissociate.
Of course, we need to know the enthalpy of homogeneous hydrate formation since the released enthalpy of hydrate formation will generate local temperature gradients.
Homogeneous CH 4 hydrate formation is discussed in Section 7.1 and homogeneous CO 2 hydrate formation is discussed in Section 7.2.

Methane Hydrate Formation from Water Solution
Solubility of CO 2 in water is significantly larger than solubility of CH 4 in water.Both of these hydrate formers can form hydrate homogeneously in water.They can also form hydrate heterogeneously since dissolved hydrate formers in water can adsorb on a hydrate film separating liquid water from the guest phase.Density of CH 4 hydrate is lower that water density and will reside on top of liquid water.A hydrate film which has been formed heterogeneously, as discussed in Section 6, might already exist.
Generally, enthalpies of hydrate formation from guest molecules dissolved in water are smaller in absolute value than enthalpies of hydrate formation from gas and liquid water.The reason is that contributions to enthalpy changes for guest molecules are smaller.Partial molar volumes of these guest molecules in water are significantly smaller than the molar volumes for the same molecules in gas phase.The guest interactions with surrounding water are also closer to that of interactions between guest and water in a closed hydrate cavity [5,53,54].
A residual thermodynamic formulation for methane dissolved in water was developed by Kvamme [5] and briefly the resulting equation is: where T R is temperature divided by critical temperature for methane.The maximum usable temperature for Equation ( 23) is 325 K, which is of course far above most hydrate forming condition.The ideal gas term on right hand side is calculated using infinite dilution partial molar volume for CH 4 as function of temperature.Pressure dependency is negligible since water is almost incompressible.The ideal gas term is therefore calculated using the analytical solution for translational movement of CH 4 in momentum space since CH 4 is approximated as a sphere.
The activity coefficient for methane in water is based on the asymmetric excess convention.In this convention, the activity coefficient for CH 4 in water approaches unity when mole-fraction CH 4 goes to zero.A fitted equation, and corresponding parameters, is given in Table 2.
. The critical temperature for CH 4 is 190.6 K.
The lower summation 1,2 indicates starting from 1 and counting in steps of 2. Parameters for Equation (24).The density of CH 4 at ideal gas reference state is almost constant.With constant density the ideal gas chemical potential can be approximated by the following fit: The homogeneous hydrate formation region is limited by solubility of methane in water, and minimum mole-fraction of CH 4 in water needed to keep hydrate stable.Solubility in liquid water is calculated by solving gas chemical potential, Equation (7), being equal to dissolved CH 4 chemical potential for CH 4 according to Equation ( 22).Minimum mole-fraction for keeping the hydrate stable is calculated by inserting Equation (22) into Equation ( 4) and then assuming equilibrium (or rather, the stability limit) between hydrate CH 4 and CH 4 in solution.For defined temperature and pressure, we can then solve for equal water chemical potentials in liquid, Equation (2), and hydrate, Equation (3).In Figure 9, I plot CH 4 solubility and CH 4 hydrate stability limit as function of mole-fraction CH 4 in surrounding water.These are plotted for various mole-fractions ethanol in water.Solubility plots are shown for the same range of conditions of homogeneous hydrate formation.The density of CH4 at ideal gas reference state is almost constant.With constant density the ideal gas chemical potential can be approximated by the following fit: The homogeneous hydrate formation region is limited by solubility of methane in water, and minimum mole-fraction of CH4 in water needed to keep hydrate stable.Solubility in liquid water is calculated by solving gas chemical potential, Equation (7), being equal to dissolved CH4 chemical potential for CH4 according to Equation (22).
Minimum mole-fraction for keeping the hydrate stable is calculated by inserting Equation (22) into Equation ( 4) and then assuming equilibrium (or rather, the stability limit) between hydrate CH4 and CH4 in solution.For defined temperature and pressure, we can then solve for equal water chemical potentials in liquid, Equation (2), and hydrate, Equation (3).In Figure 9, I plot CH4 solubility and CH4 hydrate stability limit as function of mole-fraction CH4 in surrounding water.These are plotted for various mole-fractions ethanol in water.Solubility plots are shown for the same range of conditions of homogeneous hydrate formation.(a) Solubility for 0 ethanol, (b) hydrate stability limit for 0 ethanol, (c) solubility for 0.015 mole-fraction ethanol in water, (d) hydrate stability limit for 0.015 mole-fraction ethanol in water, (e) solubility for 0.02 mole-fraction ethanol in water, (f) hydrate stability limit for 0.02 mole-fraction ethanol in water, (g) solubility for 0.03 mole-fraction ethanol in water, and (f) hydrate stability limit for 0.03 mole-fraction ethanol in water.
These 3D plots are not easy to read even with grids on.However, in the context of this work they are mainly intended as illustrations of the concentration window for possible hydrate grow from dissolved CH4 in water.As mentioned earlier, the hydrates formed from solution of CH4 are different from the heterogeneous hydrate discussed in Section 6.This can be seen by comparing Figure 10 below and gas CH4 chemical potentials (a) Solubility for 0 ethanol, (b) hydrate stability limit for 0 ethanol, (c) solubility for 0.015 mole-fraction ethanol in water, (d) hydrate stability limit for 0.015 mole-fraction ethanol in water, (e) solubility for 0.02 mole-fraction ethanol in water, (f) hydrate stability limit for 0.02 mole-fraction ethanol in water, (g) solubility for 0.03 mole-fraction ethanol in water, and (h) hydrate stability limit for 0.03 mole-fraction ethanol in water.These 3D plots are not easy to read even with grids on.However, in the context of this work they are mainly intended as illustrations of the concentration window for possible hydrate grow from dissolved CH 4 in water.As mentioned earlier, the hydrates formed from solution of CH 4 are different from the heterogeneous hydrate discussed in Section 6.This can be seen by comparing Figure 10 below and gas CH 4 chemical potentials from Figure 3b.This will result in different cavity partition functions, e.g., Equation (4).Compositions of hydrate will therefore be different and then also Gibbs free energies are different.The differences between heterogeneous and homogenous hydrate can be seen by comparing Figure 10 and gas CH4 chemical potentials from Figure 3b.Since the gas chemical potential and the liquid water chemical potential of CH4 are different in a non-equilibrium system then cavity partition functions, such as in Equation ( 3), are also different.As a consequence, Gibbs free energies for heterogeneous and homogeneous hydrates are different.By thermodynamic definition hydrates with different composition and different Gibbs free energy are a different phase.In Figure 11, I plot free energy for the homogeneously formed hydrates for 0 ethanol and 3 mole percent ethanol in water.Hydrate can form in the whole concentration range between solubility and hydrate stability limits.Each mole-fraction guest in water will result in a unique hydrate composition.Additionally, note from the plots in Figures 9-11 that pressure dependency in the hydrate plots is very small since these are phase transitions that only involve condensed phases.The differences between heterogeneous and homogenous hydrate can be seen by comparing Figure 10 and gas CH 4 chemical potentials from Figure 3b.Since the gas chemical potential and the liquid water chemical potential of CH 4 are different in a nonequilibrium system then cavity partition functions, such as in Equation ( 3), are also different.As a consequence, Gibbs free energies for heterogeneous and homogeneous hydrates are different.By thermodynamic definition hydrates with different composition and different Gibbs free energy are a different phase.In Figure 11, I plot free energy for the homogeneously formed hydrates for 0 ethanol and 3 mole percent ethanol in water.Hydrate can form in the whole concentration range between solubility and hydrate stability limits.Each mole-fraction guest in water will result in a unique hydrate composition.Additionally, note from the plots in Figures 9-11 that pressure dependency in the hydrate plots is very small since these are phase transitions that only involve condensed phases.The differences between heterogeneous and homogenous hydrate can be seen by comparing Figure 10 and gas CH4 chemical potentials from Figure 3b.Since the gas chemical potential and the liquid water chemical potential of CH4 are different in a non-equilibrium system then cavity partition functions, such as in Equation ( 3), are also different.As a consequence, Gibbs free energies for heterogeneous and homogeneous hydrates are different.By thermodynamic definition hydrates with different composition and different Gibbs free energy are a different phase.In Figure 11, I plot free energy for the homogeneously formed hydrates for 0 ethanol and 3 mole percent ethanol in water.Hydrate can form in the whole concentration range between solubility and hydrate stability limits.Each mole-fraction guest in water will result in a unique hydrate composition.Additionally, note from the plots in Figures 9-11 that pressure dependency in the hydrate plots is very small since these are phase transitions that only involve condensed phases.A very important practical implication of the hydrate stability limits towards outside water methane concentration is the role it plays in destabilizing offshore hydrate.Most offshore hydrate reservoirs worldwide are in a stationary balance between formation of new hydrate from upcoming gas, and dissociation towards incoming seawater.Fracture systems in offshore hydrate frequently connect hydrate filled sections to the seafloor.The CH 4 concentration in this seawater is close to infinite dilution.The chemical potential for CH 4 at infinite dilution in water is lower than chemical potential for CH 4 in hydrate.As a result, the hydrate will dissociate due to a more beneficial situation for CH 4 in the outside water.The result of this mechanism for hydrate dissociation is observed worldwide as fluxes of methane from natural gas hydrate deposits.The flux rates can range from invisible bubbles to massive gas fluxes that bubble through the water column.Another example is all the hydrate mounds that sit on top of fractures that connect from seafloor and downwards towards conventional hydrocarbon reservoirs.Hydrate forms from upcoming gas and seawater, due to the temperature and pressure.The formed hydrate dissociates from the top due to chemical potential gradients of CH 4 between hydrate and seawater.The hydrate formation/dissociation couplings are kinetically limited by the hydrate dissociation and these mounds can seem stationary.Dissociation fluxes can be monitored by measuring local CH 4 concentrations in surrounding water.Since the CH 4 is an attractive source for several biological ecosystems the biological activity on these mounds is normally high.A complete dynamic modelling of CH 4 fluxes from these mounds will therefore include dynamics related the thermodynamics of hydrate dissociation towards seawater, as well as dynamics of biological consumption of CH 4 .
Homogeneous hydrate formation will happen from injection gas when CO 2 /N 2 / surfactants are used for hydrate production.This can be significant but is likely not very important.This route to hydrate formation is likely more important for CH 4 .Release of CH 4 from a mole-fraction of roughly 0.14 in a hydrate will result in CH 4 bubbles.However, it is also likely that the water surrounding the dissociating hydrate will be super-saturated with CH 4 .Possible CH 4 reformation nucleation can therefore potentially happen heterogeneously as well as homogeneously.Examples for enthalpies of hydrate formation for two conditions of temperature and pressure are plotted in Figure 12 below.A very important practical implication of the hydrate stability limits towards outside water methane concentration is the role it plays in destabilizing offshore hydrate.Most offshore hydrate reservoirs worldwide are in a stationary balance between formation of new hydrate from upcoming gas, and dissociation towards incoming seawater.Fracture systems in offshore hydrate frequently connect hydrate filled sections to the seafloor.The CH4 concentration in this seawater is close to infinite dilution.The chemical potential for CH4 at infinite dilution in water is lower than chemical potential for CH4 in hydrate.As a result, the hydrate will dissociate due to a more beneficial situation for CH4 in the outside water.The result of this mechanism for hydrate dissociation is observed worldwide as fluxes of methane from natural gas hydrate deposits.The flux rates can range from invisible bubbles to massive gas fluxes that bubble through the water column.Another example is all the hydrate mounds that sit on top of fractures that connect from seafloor and downwards towards conventional hydrocarbon reservoirs.Hydrate forms from upcoming gas and seawater, due to the temperature and pressure.The formed hydrate dissociates from the top due to chemical potential gradients of CH4 between hydrate and seawater.The hydrate formation/dissociation couplings are kinetically limited by the hydrate dissociation and these mounds can seem stationary.Dissociation fluxes can be monitored by measuring local CH4 concentrations in surrounding water.Since the CH4 is an attractive source for several biological ecosystems the biological activity on these mounds is normally high.A complete dynamic modelling of CH4 fluxes from these mounds will therefore include dynamics related the thermodynamics of hydrate dissociation towards seawater, as well as dynamics of biological consumption of CH4.
Homogeneous hydrate formation will happen from injection gas when CO2/N2/surfactants are used for hydrate production.This can be significant but is likely not very important.This route to hydrate formation is likely more important for CH4.Release of CH4 from a mole-fraction of roughly 0.14 in a hydrate will result in CH4 bubbles.However, it is also likely that the water surrounding the dissociating hydrate will be super-saturated with CH4.Possible CH4 reformation nucleation can therefore potentially happen heterogeneously as well as homogeneously.Examples for enthalpies of hydrate formation for two conditions of temperature and pressure are plotted in Figure 12 below.

Carbon Dioxide Hydrate Formation from Water Solution
As discussed in the previous section, hydrate mounds on the seafloor can be connected to hydrate deposits.These mounds can, however, also be connected to conventional hydrocarbons sources that leak towards seafloor through fractures.When these hydrocarbons enter the seafloor at hydrate forming conditions, they create a dynamically unstable situation.These hydrate mounds are dynamically controlled by biological activity and also by chemical potential for CH 4 in surrounding seawater.As a result of the biological activity also CO 2 will be generated.CO 2 hydrate formation from released CO 2 will often be slow enough to dissolve in water.Seawater enriched with CO 2 is heavier than regular seawater and will sink.In either case, CO 2 hydrate can form locally from the dissolved CO 2 in water.Although the main focus herein is on hydrate production using CO 2 /N 2 mixtures modelling here, and the follow-up paper on kinetic modelling will also be useful for other groups involved in studies of hydrocarbon fluxes and the resulting changes in the ocean carbon cycle.
The solubility of N 2 in water is extremely low and of no practical importance within the scope of this work.With 20 mole% N 2 in the CO 2 /N 2 mixture the chemical potential for CO 2 "bulk" gas can be calculated using SRK [14].On the other hand, the properties of CO 2 and N 2 are extremely different in terms of selective adsorption on liquid water surface [32].It might therefore be a fair approximation to calculate solubility based on almost pure CO 2 in the adsorbed water layer.
The density of CO 2 as dissolved in water will correspond to the partial molar volume of CO 2 at infinite dilution.The density of CO 2 at infinite dilution is simply the inverse of the partial molar volume.Chemical potential for ideal gas is a trivial analytical function of the center of mass and molecular weight and the moments of inertia for rotation.It is the Boltzmann factor integrals of the momentum space in the canonical partition function.The integrals end up as analytical simple expressions that can be found in any textbook on physical chemistry or statistical mechanics.The infinite dilution ideal gas chemical potential is not very sensitive to pressure.The following fitted approximation to only temperature dependency is considered as adequate: where T 0,R is 273.15K divided by the actual temperature.Equation (26) does not apply to temperatures above 303 K due to the limited range of temperatures for which infinite partial molar volumes are used and for temperatures above 273.15K.Many methods for production of hydrate reservoirs have been proposed during the latest three decades.Thermal stimulation, or injection of thermodynamic hydrate inhibitors are just two examples of methods that have more or less lost attention due to the high cost.One problem is, however, that different methods for producing hydrates are not evaluated thermodynamically prior to planning expensive experiments or even more expensive pilot tests.This can be due to lack of a thermodynamic toolbox for the purpose.Another challenge is the lack of focus on the limitations of the hydrate phase transition itself.The interface between hydrate and liquid water is a kinetic bottle neck.Reducing pressure does not address this problem.Injection of CO 2 will lead to formation of a new CO 2 hydrate.This hydrate formation is an efficient heat source for dissociating hydrate since heating will break the water hydrogen bonds and directly address the important nano scale kinetic limitation.Adding limited amounts of N 2 increases permeability for the injection gas.The addition of surfactant increases gas/water interface dynamics and promote heterogeneous hydrate formation.In this work, we demonstrate a residual thermodynamic scheme that opens up for thermodynamic analysis of different routes to hydrate formation and dissociation.It is demonstrated that 20 mole% N 2 added to the CO 2 is thermodynamically feasible for generating a new hydrate in the pores.The available hydrate formation enthalpy from CO 2 /N 2 hydrate formation is also feasible for dissociating in situ CH 4 hydrate.Up to 3 mole% ethanol in the free pore water is also thermodynamically feasible.The addition of alcohol will not significantly disturb the ability to form new hydrate from the injection gas.Homogeneous hydrate formation from dissolved CH 4 and/or CO 2 is limited in amount and not important.However, the hydrate stability limits related to concentration of hydrate former in surrounding water is important.Mineral surfaces can act as hydrate promotors through direct adsorption, or adsorption in water which is structured by the mineral surface charges.These aspects will be quantified in a follow-up paper, along with kinetic modelling based on thermodynamic modelling in this work.
For CO 2 in water, I utilize the following function: where T R is reduced temperature and defined as actual T in Kelvin divided by critical temperature for CO 2 (304.35K).The lower summation 1, 2 indicates starting from 1 and counting in steps of 2. Parameters are given in Table 3 below.The arrow on top of x denotes the vector of all mole-fractions in the actual phase.The chemical potential for CO 2 that applies to Equation ( 23) for an equilibrium case is then given as: In Figure 13a, I plot solubility of CO 2 in pure water.Hydrate stability limits for CO 2 as function of mole-fractions CO 2 in surrounding water is plotted in Figure 13b.Similar results for water containing 3 mole% ethanol are plotted in Figure 14.In order to get a picture of the mole-fraction range for homogeneous hydrate growth from CO2 dissolved in water, I plot solubility in water minus hydrate stability limit in water as a function of mole-fraction CO2 in water from solubility mole-fraction and down to hydrate stability limit mole-fraction.These are plotted in Figure 15 as function of solubility and temperature.Solubility axis reflects a pressure range from 90 bar to 300 bar.In order to get a picture of the mole-fraction range for homogeneous hydrate growth from CO2 dissolved in water, I plot solubility in water minus hydrate stability limit in water as a function of mole-fraction CO2 in water from solubility mole-fraction and down to hydrate stability limit mole-fraction.These are plotted in Figure 15 as function of solubility and temperature.Solubility axis reflects a pressure range from 90 bar to 300 bar.In order to get a picture of the mole-fraction range for homogeneous hydrate growth from CO 2 dissolved in water, I plot solubility in water minus hydrate stability limit in water as a function of mole-fraction CO 2 in water from solubility mole-fraction and down to hydrate stability limit mole-fraction.These are plotted in Figure 15 as function of solubility and temperature.Solubility axis reflects a pressure range from 90 bar to 300 bar.
The stability of hydrate formed homogeneously from dissolved CO 2 in water is lower than heterogeneously formed CO 2 hydrate, as illustrated in Figure 16.The stability of hydrate formed homogeneously from dissolved CO2 in water is lower than heterogeneously formed CO2 hydrate, as illustrated in Figure 16.The change in the enthalpy model for CO2 is the same that was done for CH4 (but not shown here).It simply involves replacing the partial molar enthalpy from a gas phase as in the heterogeneous case to partial molar enthalpy for guest molecules dissolved in water.Specifically for CO2:  The stability of hydrate formed homogeneously from dissolved CO2 in water is lower than heterogeneously formed CO2 hydrate, as illustrated in Figure 16.The change in the enthalpy model for CO2 is the same that was done for CH4 (but not shown here).It simply involves replacing the partial molar enthalpy from a gas phase as in the heterogeneous case to partial molar enthalpy for guest molecules dissolved in water.Specifically for CO2: The change in the enthalpy model for CO 2 is the same that was done for CH 4 (but not shown here).It simply involves replacing the partial molar enthalpy from a gas phase as in the heterogeneous case to partial molar enthalpy for guest molecules dissolved in water.Specifically for CO 2 : Examples for two conditions of temperatures and pressures are plotted in Figure 17.The enthalpies of hydrate formation are small due to smaller changes for guest molecules as compared to changes for a gas molecule that gets trapped into a higher density.Practically this homogeneous hydrate formation may not be critically important for the hydrate production planning.The hydrate stability limits towards CO2 in surrounding water may, however, be more important in terms of long terms stability of the stored CO2.Nevertheless-storing CO2 as solid hydrate in a reservoir with proven sealing is a safe option.The in situ CH4 hydrate may have been trapped for millions of years due to sealing formations (clay and shale).

Discussion
Interest for using CO2 in a concept for combined safe long-term storage of CO2, and simultaneous release of CH4 is increasing worldwide.Associated thermodynamic properties and knowledge of conversion mechanisms are critical.I have many years of experience on the water/gas interface dynamics and how it affects conditions for hydrate nucleation on liquid water side of the interface.All the papers in the PhD theses [49][50][51][52] and other references here [5,7,53] support the physical picture that heterogeneous hydrate nucleation is facilitated by water structures on the liquid water side of the interface.Recent studies on surfactant stimulated water/carbon dioxide interfaces underscore the importance of surfactants for hydrate nucleation [3,10,55,56].There is a fundamental difference between heterogeneous hydrate nucleation in liquid water region, and at conditions far into the ice region.CH4 hydrates made from powdered ice and CH4, and then exposed to CO2 at ice conditions do not have a liquid phase.There will for sure be some sort of interface between ice and hydrate, but conversion from CH4 hydrate to CO2 hydrate does not have much flexibility.It is understandable that experimental cells in atomistic imaging equipment cannot handle high pressures.It is therefore no surprise that hydrate exchange The enthalpies of hydrate formation are small due to smaller changes for guest molecules as compared to changes for a gas molecule that gets trapped into a higher density.Practically this homogeneous hydrate formation may not be critically important for the hydrate production planning.The hydrate stability limits towards CO 2 in surrounding water may, however, be more important in terms of long terms stability of the stored CO 2 .Nevertheless-storing CO 2 as solid hydrate in a reservoir with proven sealing is a safe option.The in situ CH 4 hydrate may have been trapped for millions of years due to sealing formations (clay and shale).

Discussion
Interest for using CO 2 in a concept for combined safe long-term storage of CO 2 , and simultaneous release of CH 4 is increasing worldwide.Associated thermodynamic properties and knowledge of conversion mechanisms are critical.I have many years of experience on the water/gas interface dynamics and how it affects conditions for hydrate nucleation on liquid water side of the interface.All the papers in the PhD theses [49][50][51][52] and other references here [5,7,53] support the physical picture that heterogeneous hydrate nucleation is facilitated by water structures on the liquid water side of the interface.Recent studies on surfactant stimulated water/carbon dioxide interfaces underscore the importance of surfactants for hydrate nucleation [3,10,55,56].There is a fundamental difference between heterogeneous hydrate nucleation in liquid water region, and at conditions far into the ice region.CH 4 hydrates made from powdered ice and CH 4 , and then exposed to CO 2 at ice conditions do not have a liquid phase.There will for sure be some sort of interface between ice and hydrate, but conversion from CH 4 hydrate to CO 2 hydrate does not have much flexibility.It is understandable that experimental cells in atomistic imaging equipment cannot handle high pressures.It is therefore no surprise that hydrate exchange experiments at atmospheric conditions and very low temperatures [51] reveal a solid-state exchange mechanism.There is no liquid water and a totally different ice/hydrate interface as compared to the dynamic and broad (relative) hydrate/liquid water interface, and also a dynamic interface between liquid water and gas.I have been quite purposeful in not referring to experimental data, including our own data.Kvamme et al. [11] and Tegze et al. [12] are referred to in this work for interface details and modelling aspects.However, these references also contain information from experiments.A solid-state mechanism is a totally impossible explanation of our experiments over 15 years.These experiments were all conducted at liquid water conditions.An important element of Baig's study [57] was to shed light on conversion mechanisms.The results revealed a fast conversion mechanism based on creation of new CO 2 hydrate as long as free water was available.After the stage when

Figure 1 .
Figure 1.(a) Temperature pressure stability limits for CH4 hydrate formed from CH4 and water containing MEG (circles are experimental data [38] and solid curves are calculated) or TEG (squares are experimental data [39] and dashed curves are calculated).Blue curves and circles are for 10 wt% MEG in water.Red circles and curve are for 30 wt% MEG. Green squares and dashed curve are for 10 wt% TEG.Red squares and dashed curve are for 20 wt% TEG and blue squares and dashed curve are for 40 wt% TEG.(b) Gibbs free energy for hydrates formed along the stability curves in Figure 1a.Same line notations as in Figure 1a.

Figure 1 .Figure 2 .
Figure 1.(a) Temperature pressure stability limits for CH 4 hydrate formed from CH 4 and water containing MEG (circles are experimental data [38] and solid curves are calculated) or TEG (squares are experimental data [39] and dashed curves are calculated).Blue curves and circles are for 10 wt% MEG in water.Red circles and curve are for 30 wt% MEG. Green squares and dashed curve are for 10 wt% TEG.Red squares and dashed curve are for 20 wt% TEG and blue squares and dashed curve are for 40 wt% TEG.(b) Gibbs free energy for hydrates formed along the stability curves in Figure 1a.Same line notations as in Figure 1a.Fluids 2021, 6, x FOR PEER REVIEW 10 of 31

Figure 2 .
Figure 2. (a) Temperature pressure stability limits for CH 4 hydrate formed from CH 4 and water containing ethanol (circles are experimental data[40] and solid curves are calculated) and methanol (squares are experimental data[41] and dashed curves are calculated).Blue circles and solid blue curves are for 15 wt% ethanol.Green squares and dashed curve are for 10 wt% methanol.Dashed blue is for 15 wt% methanol.Red squares and dashed are for 20 wt% methanol.Solid black curve is calculated data for pure water.Crosses are experimental data from Sabil et al.[42] and plusses are experimental data from Tumba et al.[43].(b) Gibbs free energy for hydrates formed from CH 4 and water containing various concentrations of ethanol and methanol.Same notations as in Figure2a.

FluidsFigure 3 .
Figure 3. (a) Temperature pressure stability limits for CO2 hydrate formed from CO2 and water containing ethanol.Black curve is for 0 mole-fraction ethanol, green is for 0.015 mole-fraction ethanol, blue is for 0.02 mole-fraction ethanol, and red is for 0.03 mole-fraction ethanol.(b) Gibbs free energy for hydrates formed from CO2 and water containing various concentrations of ethanol.The color codes are the same as in Figure 3a.

Figure 4 .
Figure 4. (a) Mole-fraction CO2 in CO2 hydrate formed from CO2 and water containing ethanol.Black curve is for 0 molefraction ethanol, green curve is for 0.015 mole-fraction ethanol, blue curve is for 0.02 mole-fraction ethanol, and red curve is for 0.03 mole-fraction ethanol.(b) Chemical potential for CO2 in CO2 hydrate formed from CO2 and water containing various concentrations of ethanol.Color codes are the same as in Figure 4a.

Figure 3 .
Figure 3. (a) Temperature pressure stability limits for CO 2 hydrate formed from CO 2 and water containing ethanol.Black curve is for 0 mole-fraction ethanol, green is for 0.015 mole-fraction ethanol, blue is for 0.02 mole-fraction ethanol, and red is for 0.03 mole-fraction ethanol.(b) Gibbs free energy for hydrates formed from CO 2 and water containing various concentrations of ethanol.The color codes are the same as in Figure 3a.

FluidsFigure 3 .
Figure 3. (a) Temperature pressure stability limits for CO2 hydrate formed from CO2 and water containing ethanol.Black curve is for 0 mole-fraction ethanol, green is for 0.015 mole-fraction ethanol, blue is for 0.02 mole-fraction ethanol, and red is for 0.03 mole-fraction ethanol.(b) Gibbs free energy for hydrates formed from CO2 and water containing various concentrations of ethanol.The color codes are the same as in Figure 3a.

Figure 4 .
Figure 4. (a) Mole-fraction CO2 in CO2 hydrate formed from CO2 and water containing ethanol.Black curve is for 0 molefraction ethanol, green curve is for 0.015 mole-fraction ethanol, blue curve is for 0.02 mole-fraction ethanol, and red curve is for 0.03 mole-fraction ethanol.(b) Chemical potential for CO2 in CO2 hydrate formed from CO2 and water containing various concentrations of ethanol.Color codes are the same as in Figure 4a.

Figure 4 .
Figure 4. (a) Mole-fraction CO 2 in CO 2 hydrate formed from CO 2 and water containing ethanol.Black curve is for 0 mole-fraction ethanol, green curve is for 0.015 mole-fraction ethanol, blue curve is for 0.02 mole-fraction ethanol, and red curve is for 0.03 mole-fraction ethanol.(b) Chemical potential for CO 2 in CO 2 hydrate formed from CO 2 and water containing various concentrations of ethanol.Color codes are the same as in Figure 4a.

Fluids 2021, 6 ,Figure 5 .
Figure 5. (a) Temperature pressure stability limits for mixed CO2/N2 hydrate formed from a CO2/N2 mixture with 20 mole percent N2 and water containing ethanol.Black curve is for 0 mole-fraction ethanol, green curve is for 0.015 mole-fraction ethanol, blue curve is for 0.02 mole-fraction ethanol, and red curve is for 0.03 mole-fraction ethanol.(b) Gibbs free energy for hydrates formed from a CO2/N2 mixture with 20 mole percent N2 and water containing various concentrations of ethanol.The color notation is the same as in Figure 5a.

Figure 5 .
Figure 5. (a) Temperature pressure stability limits for mixed CO 2 /N 2 hydrate formed from a CO 2 /N 2 mixture with 20 mole percent N 2 and water containing ethanol.Black curve is for 0 mole-fraction ethanol, green curve is for 0.015 mole-fraction ethanol, blue curve is for 0.02 mole-fraction ethanol, and red curve is for 0.03 mole-fraction ethanol.(b) Gibbs free energy for hydrates formed from a CO 2 /N 2 mixture with 20 mole percent N 2 and water containing various concentrations of ethanol.The color notation is the same as in Figure 5a.

Fluids 2021, 6 ,Figure 6 .
Figure 6.(a) Hydrate composition in mixed CO2/N2 hydrate formed from a CO2/N2 mixture containing 20 mole percent N2 and water containing ethanol.Dashed curves are N2 and solid curves are CO2.Black curves are for 0 mole-fraction ethanol, green curves are for 0.015 mole-fraction ethanol, blue curves are for 0.02 mole-fraction ethanol, and red curves are for 0.03 mole-fraction ethanol.(b) Chemical potentials for CO2 (solid) and N2 (dashed curves) in hydrates formed from a CO2/N2 mixture with 20 mole percent N2 and water containing various concentrations of ethanol.Color codes are the same as in Figure 6a.

Figure 6 .
Figure 6.(a) Hydrate composition in mixed CO 2 /N 2 hydrate formed from a CO 2 /N 2 mixture containing 20 mole percent N 2 and water containing ethanol.Dashed curves are N 2 and solid curves are CO 2 .Black curves are for 0 mole-fraction ethanol, green curves are for 0.015 mole-fraction ethanol, blue curves are for 0.02 mole-fraction ethanol, and red curves are for 0.03 mole-fraction ethanol.(b) Chemical potentials for CO 2 (solid) and N 2 (dashed curves) in hydrates formed from a CO 2 /N 2 mixture with 20 mole percent N 2 and water containing various concentrations of ethanol.Color codes are the same as in Figure 6a.

Figure 7 .
Figure 7. (a) Enthalpy of hydrate formation per mole hydrate as function of temperature and mole-fraction ethanol in water for pure CO2 hydrate.The lowest curve is for no ethanol in water, then 0.015 mole-fraction ethanol in water, then 0.02 mole-fraction ethanol in water, and upper curve 0.03 mole-fraction ethanol in water.(b) Enthalpy of hydrate formation per mole guest in hydrate as function of temperature and mole-fraction ethanol in water for pure CO2 hydrate.The lowest curve is for no ethanol in water, then 0.015 mole-fraction ethanol in water, then 0.02 mole-fraction ethanol in water, and upper curve 0.03 mole-fraction ethanol in water.(c) Enthalpy of hydrate formation per mole hydrate as function of temperature and mole-fraction ethanol in water for hydrate formed from 20 mole percent N2 in CO2.The lowest curve is for no ethanol in water, then 0.015 mole-fraction ethanol in water, then 0.02 mole-fraction ethanol in water, and upper curve 0.03 mole-fraction ethanol in water.(d) Enthalpy of hydrate formation per mole guest in hydrate as function of temperature and mole-fraction ethanol in water for hydrate formed from 20 mole percent N2 in CO2.The lowest curve is for no ethanol in water, then 0.015 mole-fraction ethanol in water, then 0.02 mole-fraction ethanol in water, and upper curve 0.03 mole-fraction ethanol in water.

Figure 7 .
Figure 7. (a) Enthalpy of hydrate formation per mole hydrate as function of temperature and mole-fraction ethanol in water for pure CO 2 hydrate.The lowest curve is for no ethanol in water, then 0.015 mole-fraction ethanol in water, then 0.02 mole-fraction ethanol in water, and upper curve 0.03 mole-fraction ethanol in water.(b) Enthalpy of hydrate formation per mole guest in hydrate as function of temperature and mole-fraction ethanol in water for pure CO 2 hydrate.The lowest curve is for no ethanol in water, then 0.015 mole-fraction ethanol in water, then 0.02 mole-fraction ethanol in water, and upper curve 0.03 mole-fraction ethanol in water.(c) Enthalpy of hydrate formation per mole hydrate as function of temperature and mole-fraction ethanol in water for hydrate formed from 20 mole percent N 2 in CO 2 .The lowest curve is for no ethanol in water, then 0.015 mole-fraction ethanol in water, then 0.02 mole-fraction ethanol in water, and upper curve 0.03 mole-fraction ethanol in water.(d) Enthalpy of hydrate formation per mole guest in hydrate as function of temperature and mole-fraction ethanol in water for hydrate formed from 20 mole percent N 2 in CO 2 .The lowest curve is for no ethanol in water, then 0.015 mole-fraction ethanol in water, then 0.02 mole-fraction ethanol in water, and upper curve 0.03 mole-fraction ethanol in water.

Figure 8 .
Figure 8.(a) Enthalpy of hydrate formation per mole hydrate and (b) enthalpy of hydrate formation per mole guest.Both figures are plotted as function of ethanol concentration in water.Dashed curves are for CH4 hydrates formed from water with various mole-fractions ethanol in water.Solid curves are for hydrates formed from CO2 with 20 mole% N2.Lowest curves are for no ethanol, then 0.015 mole-fraction ethanol in water, then 0.02 mole-fraction ethanol in water, and upper curves for 0.03 mole-fraction ethanol in water.

Figure 8 .
Figure 8.(a) Enthalpy of hydrate formation per mole hydrate and (b) enthalpy of hydrate formation per mole guest.Both figures are plotted as function of ethanol concentration in water.Dashed curves are for CH 4 hydrates formed from water with various mole-fractions ethanol in water.Solid curves are for hydrates formed from CO 2 with 20 mole% N 2 .Lowest curves are for no ethanol, then 0.015 mole-fraction ethanol in water, then 0.02 mole-fraction ethanol in water, and upper curves for 0.03 mole-fraction ethanol in water.

Figure 9 .
Figure 9. Liquid water solubility and hydrate stability limit mole-fractions of hydrate for various mole-fractions of ethanol.(a)Solubility for 0 ethanol, (b) hydrate stability limit for 0 ethanol, (c) solubility for 0.015 mole-fraction ethanol in water, (d) hydrate stability limit for 0.015 mole-fraction ethanol in water, (e) solubility for 0.02 mole-fraction ethanol in water, (f) hydrate stability limit for 0.02 mole-fraction ethanol in water, (g) solubility for 0.03 mole-fraction ethanol in water, and (h) hydrate stability limit for 0.03 mole-fraction ethanol in water.

Fluids 2021, 6 ,Figure 10 .
Figure 10.Chemical potential of CH4 entering hydrate as function of temperature pressure and mole-fractions CH4 (see Figure 10) for: (a) 0 ethanol and (b) 3 mole percent ethanol in water.

Figure 11 .
Figure 11.Gibbs free energy for homogeneously formed CH4 hydrate as function of temperature pressure and mole-fractions CH4 (see Figure 10) for: (a) 0 ethanol and (b) 3 mole percent ethanol in water.

Figure 10 .
Figure 10.Chemical potential of CH 4 entering hydrate as function of temperature pressure and mole-fractions CH 4 (see Figure 10) for: (a) 0 ethanol and (b) 3 mole percent ethanol in water.

Fluids 2021, 6 ,Figure 10 .
Figure 10.Chemical potential of CH4 entering hydrate as function of temperature pressure and mole-fractions CH4 (see Figure 10) for: (a) 0 ethanol and (b) 3 mole percent ethanol in water.

Figure 11 .
Figure 11.Gibbs free energy for homogeneously formed CH4 hydrate as function of temperature pressure and mole-fractions CH4 (see Figure 10) for: (a) 0 ethanol and (b) 3 mole percent ethanol in water.Figure 11.Gibbs free energy for homogeneously formed CH 4 hydrate as function of temperature pressure and mole-fractions CH 4 (see Figure 10) for: (a) 0 ethanol and (b) 3 mole percent ethanol in water.

Figure 11 .
Figure 11.Gibbs free energy for homogeneously formed CH4 hydrate as function of temperature pressure and mole-fractions CH4 (see Figure 10) for: (a) 0 ethanol and (b) 3 mole percent ethanol in water.Figure 11.Gibbs free energy for homogeneously formed CH 4 hydrate as function of temperature pressure and mole-fractions CH 4 (see Figure 10) for: (a) 0 ethanol and (b) 3 mole percent ethanol in water.

Figure 13 .Figure 14 .
Figure 13.Pure water solubility for CO2 and hydrate stability limits for CO2.(a) Solubility of CO2 in water for the range of hydrate forming conditions.(b) Hydrate stability limit concentrations for CO2.

Figure 13 .
Figure 13.Pure water solubility for CO 2 and hydrate stability limits for CO 2 .(a) Solubility of CO 2 in water for the range of hydrate forming conditions.(b) Hydrate stability limit concentrations for CO 2 .

Fluids 2021, 6 ,Figure 13 .Figure 14 .
Figure 13.Pure water solubility for CO2 and hydrate stability limits for CO2.(a) Solubility of CO2 in water for the range of hydrate forming conditions.(b) Hydrate stability limit concentrations for CO2.

Figure 14 .
Figure 14.Solubility of CO 2 in water containing 3 mole percent ethanol and hydrate stability limits for CO 2 .(a) Solubility of CO 2 in water for the range of hydrate forming conditions.(b) Hydrate stability limit concentrations for CO 2 .

Figure 15 .
Figure15.Hydrate growth potential in terms of mole-fraction, from solubility down to hydrate stability limit (z-axis), as function of solubility mole-fraction (x-axis) and temperature (y-axis) and pressure.For each temperature on the temperature axis, solubilities and growth potentials are calculated for pressures from 90 bar to 300 bar This is reflected in the increased solubility along the x-axis.(a) Pure water and (b) 3 mole% ethanol in water.

Figure 16 .
Figure 16.Hydrate Gibbs free energy for hydrate formed heterogeneously from CO2 and liquid water (solid) and homogeneously formed hydrate from dissolved CO2 in water (dashed).(a) Pure water.(b) Water containing 3 mole percent ethanol.

Figure 15 .
Figure15.Hydrate growth potential in terms of mole-fraction, from solubility down to hydrate stability limit (z-axis), as function of solubility mole-fraction (x-axis) and temperature (y-axis) and pressure.For each temperature on the temperature axis, solubilities and growth potentials are calculated for pressures from 90 bar to 300 bar This is reflected in the increased solubility along the x-axis.(a) Pure water and (b) 3 mole% ethanol in water.

Fluids 2021, 6 ,Figure 15 .
Figure15.Hydrate growth potential in terms of mole-fraction, from solubility down to hydrate stability limit (z-axis), as function of solubility mole-fraction (x-axis) and temperature (y-axis) and pressure.For each temperature on the temperature axis, solubilities and growth potentials are calculated for pressures from 90 bar to 300 bar This is reflected in the increased solubility along the x-axis.(a) Pure water and (b) 3 mole% ethanol in water.

Figure 16 .
Figure 16.Hydrate Gibbs free energy for hydrate formed heterogeneously from CO2 and liquid water (solid) and homogeneously formed hydrate from dissolved CO2 in water (dashed).(a) Pure water.(b) Water containing 3 mole percent ethanol.

Figure 16 .
Figure 16.Hydrate Gibbs free energy for hydrate formed heterogeneously from CO 2 and liquid water (solid) and homogeneously formed hydrate from dissolved CO 2 in water (dashed).(a) Pure water.(b) Water containing 3 mole percent ethanol.

Fluids 2021, 6 Figure 17 .
Figure 17.Enthalpy of hydrate formation from dissolved CO2 in water.(a) Temperature 274 K and pressure 100 bar.Solid curve is for pure water and dashed curve is for water containing 3 mole percent ethanol.(b) Temperature 280 K and pressure 180 bar.Solid curve is for pure water and dashed curve is for water containing 3 mole percent ethanol.

Figure 17 .
Figure 17.Enthalpy of hydrate formation from dissolved CO 2 in water.(a) Temperature 274 K and pressure 100 bar.Solid curve is for pure water and dashed curve is for water containing 3 mole percent ethanol.(b) Temperature 280 K and pressure 180 bar.Solid curve is for pure water and dashed curve is for water containing 3 mole percent ethanol.

Table 2 .
Activity model for CH 4 dissolved in water.