Stability and growth of gas hydrates below the ice–hydrate–gas equilibrium line on the P–T phase diagram

Using a previously developed experimental technique, the behavior of small methane and propane hydrate samples formed from water droplets between 0.25 and 2.5 mm in size has been studied in the pressure–temperature area between the ice–hydrate–gas equilibrium line and the supercooled water–hydrate–gas metastable equilibrium line, where ice is a stable phase. The unusual persistence of the hydrates within the area bounded by these lines and the isotherms at T=253 K for methane hydrate or at T=263 K for propane hydrates was observed. This behavior has not previously been reported. For example, in the experiment carried out at 1.9 MPa and 268 K, the methane hydrates existed in a metastable state (the equilibrium pressure at 268 K is 2.17 MPa) for 2 weeks, then immediately dissociated into liquid supercooled water and gas after the pressure was isothermally decreased slightly below the supercooled water–hydrate–gas metastable equilibrium pressure. It was found that dissociation of metastable hydrate into supercooled water and gas was reversible. The lateral hydrate film growth rates of metastable methane and propane hydrates on the surface of supercooled water at a pressure below the ice–hydrate–gas equilibrium pressure were measured. The temperature range within which supercooled water formed during hydrate dissociation can exist and a role of supercooled water in hydrate self-preservation is discussed.     

Evidence of liquid water formation during methane hydrates dissociation below the ice point

Dissociation of small methane hydrate samples formed from water droplets of size 0.25-2.5 mm has been investigated below the ice melting point in the temperature range of 240-273 K, where the self-preservation effect is observed for bulk hydrates. The experiments included optical microscopy observations combined with P-T measurements of the dissociation conditions for the methane hydrates. For the first time, the formation of supercooled liquid water during the hydrate dissociation was reliably detected in the temperature range of 253-273 K. The formation of the liquid phase was visually observed. The induction time of the ice nucleation for the metastable liquid water depended from the dissociation temperature and a size of water droplets formed during the hydrate dissociation. It was found that in the temperature range of 253-273 K values of the dissociation pressure for the small hydrate samples fall on the extension of the water-hydrate-gas equilibrium curve into the metastable region where supercooled water exist. The average molar enthalpy of 51.7 kJ/mol for the dissociation of the small methane hydrate samples in the temperature range of 253-273 K was calculated using Clausius-Clapeyron equation. This value agrees with the enthalpy of dissociation of bulk methane hydrates into water and gas at temperatures above 273 K.     

Cage occupancy and dissociation enthalpy of hydrocarbon hydrates

An elaborated statistical mechanical theory on clathrate hydrates is applied to exploration of their phase equilibria and dissociation enthalpies. The experimental dissociation pressures of methane, ethane, acetylene, and propane hydrates are well recovered by the method we have proposed. We estimate water/hydrate and hydrate/guest two-phase coexisting conditions in the temperature, pressure, and composition space in addition to three-phase equilibrium conditions. It is shown that the occupancy of guest molecules and the two-phase boundaries in the phase diagram vary depending sensitively on its size. Enthalpy components arising from the host and guest interactions are separately calculated from the temperature dependence of the corresponding free energy values. This enables to evaluate the dissociation enthalpy at any stable and metastable thermodynamic state taking account of the phase transition in the coexisting phase such as melting of ice, notably that along the three-phase equilibrium line.     

An improved Clapeyron model for predicting liquid water–hydrate–liquid hydrate former phase equilibria

Determination of pressure–temperature phase diagram for the liquid water–hydrate–liquid hydrate former region is a challenge considering this diagram at these conditions is a strong function of temperature. The Clapeyron model is traditionally used for this purpose. However, the conventional Clapeyron model does not take into account the effect of pressure on the hydrate molar volume as well as the heat of dissolution of hydrate former in water and therefore cannot predict satisfactorily the gas hydrate phase diagram. In this work, the conventional Clapeyron model is extended to take into account the aforementioned factors as well as the change in the slope of the pressure–temperature phase diagram when increasing the temperature. Three hydrate formers, i.e. carbon dioxide, hydrogen sulfide and ethane, are studied. It is found that the effect of these factors on the determination of the gas hydrate phase diagram is not negligible.     

Thermodynamic stability, spectroscopic identification and cage occupation of binary CO2 clathrate hydrates

The hydrate phase behavior of CO₂/3-methyl-1-butanol (3M1B)/water, CO₂/tetrahydrofuran (THF)/water and CO₂/1,4-dioxane (DXN)/water was investigated using both a high-pressure equilibrium viewing cell and a kinetic pressure-temperature measurement system with a constant volume. The dissociation pressures of CO₂/3M1B/water were identical to those of pure CO₂ hydrate, indicating that CO₂ is not acting as a help gas for structure H hydrate formation with 3M1B, thus the formed hydrate is pure CO₂ structure I hydrate. The CO₂ molecules could be encaged in small cages of the structure II hydrate framework formed with both of THF and DXN. For a stoichiometric ratio of 5.56 mol% THF, we found a large shift of dissociation boundary to lower pressures and higher temperatures from the dissociation conditions of pure CO₂ hydrate. From the measurements using the kinetic pressure-temperature system, it was found that the solid binary hydrate samples formed from off-stoichiometric THF and DXN aqueous solutions are composed of pure CO₂ hydrate with a hydrate number n=7.0 and THF/CO₂ and DXN/CO₂ binary hydrates with a molar ratio of xCO₂·THF·17H₂O and xCO₂·DXN·17H₂O, respectively. The X-ray diffraction was used to identify the binary hydrate structure and Raman spectroscopy was measured to support the phase equilibrium results and to investigate the occupation of CO₂ molecules in the cages of the hydrate framework.     

Phase behavior of clathrate hydrates: a model for single and multiple gas component hydrates

Presented here is a model that accurately predicts equilibrium pressures as a function of temperature of hydrates with CH₄, C₂H₆, C₃H₈, N₂, H₂, and CO₂ and their mixtures as guests. The model parameters fit to a subset of the equilibrium pressure data for single guest hydrates allow the prediction of phase behavior in mixed guest hydrates. For single guest hydrates, our model improves upon the van der Waals and Platteeuw (vdWP) model with a percent absolute average deviation (%AAD) from all equilibrium pressure data of 5.7% compared to 15.1% for the vdWP model. Predictions of equilibrium pressures for all available mixed guest hydrates result in a 11.6%AAD with our fugacity-based model compared to 18.6% for the vdWP model. Also, our model leads to a prediction of the structure change of the methane-ethane hydrate within 5% of its known equilibrium composition in the vapor phase without any adjustment of its parameters. We have also found that at temperatures above , double occupancy of nitrogen in the large cavity of structure II hydrate is important for the prediction of accurate equilibrium pressures.     

In situ injection of THF to trigger gas hydrate crystallization: Application to the evaluation of a kinetic hydrate promoter

This paper investigates an original method to efficiently trigger gas hydrate crystallization. This method consists of an in situ injection of a small amount of THF into an aqueous phase in contact with a gas-hydrate-former phase at pressure and temperature conditions inside the hydrate metastable zone. In the presence of a CO₂–CH₄ gas mixture, our results show that the THF injection induces immediate crystallization of a first hydrate containing THF. This triggers the formation of the CO₂–CH₄ binary hydrate as proven by the pressure and temperature reached at equilibrium. This experimental method, which “cancels out” the stochasticity of the hydrate crystallization, was used to evaluate the effect of the anionic surfactant SDS at different concentrations, on the formation kinetics of the CO₂–CH₄ hydrate. The results are discussed and compared with those published in a recent article (Ricaurte et al., 2013), where THF was not injected but present in the aqueous phase from the beginning and at much higher concentrations.     

Hydrate phase equilibria for gas mixtures containing carbon dioxide: A proof-of-concept to carbon dioxide recovery from multicomponent gas stream

Three-phase equilibrium conditions (aqueous liquid-hydrate-vapor) of CO₂-N₂ binary mixtures in the temperature range of 271.75 K to 284.25 K and the pressure range of 12 to 235 bar. In addition, three-phase (aqueous liquid-hydrate-vapor) behavior for CO₂-CH₄ mixture were measured in the temperature range of 272 to 284 K at the constant pressures of 15, 20, 26, 35 and 50 bar. In high concentration of CO₂, upper quadruple points were also measured. The obtained data indicates that three-phase equilibrium temperatures become higher with increasing concentration of CO₂. For the prediction of three-phase equilibrium, the vapor and liquid phases were treated by employing the Soave-Redlich-Kwong equation of state (SRK-EOS) with the second order modified Huron-Vidal (MHV2) mixing rule and the hydrate phase with the van der Waals-Platteeuw model. The calculated results showed good agreement with experimental data. The concentration of vapor and hydrate phases was also determined experimentally. This work can be used as the basic data for selective separation process by hydrate formation. This paper was presented at The 5th International Symposium on Separation Technology-Korea and Japan held at Seoul between August 19 and 21, 1999.     

Gas hydrate formation from hydrogen/carbon dioxide and nitrogen/carbon dioxide gas mixtures

Gas hydrates from CO₂/N₂ and CO₂/H₂ gas mixtures were formed in a semi-batch stirred vessel at constant pressure and temperature of 273.7 K. These mixtures are of interest to CO₂ separation and recovery from flue gas and fuel gas, respectively. During hydrate formation the gas uptake was determined and the composition changes in the gas phase were obtained by gas chromatography. The rate of hydrate growth from CO₂/H₂ mixtures was found to be the fastest. In both mixtures CO₂ was found to be preferentially incorporated into the hydrate phase. The observed fractionation effect is desirable and provides the basis for CO₂ capture from flue gas or fuel gas mixtures. The separation from fuel gas is also a source of H₂. The impact of tetrahydrofuran (THF) on hydrate formation from the CO₂/N₂ mixture was also observed. THF is known to substantially reduce the equilibrium formation conditions enabling hydrate formation at much lower pressures. THF was found to reduce the induction time and the rate of hydrate growth.     

Pressure variation due to heat shock of CO2hydrate desserts

CO₂ hydrate desserts are carbonated frozen desserts in which the CO₂ is trapped in a crystalline water‐carbon dioxide structure called a CO₂ clathrate hydrate. The CO₂ concentration of the dessert enables strong perception of carbonation, but CO₂ hydrate dissociation during heat shock can cause high package pressures during storage and distribution. In this work, a model is developed for package pressure as a function of temperature, CO₂ content, package volume, dessert mass, and recipe. The model is validated by comparison with an experimental measurement of the pressure and mass of a CO₂ hydrate dessert subjected to heat shock. It is shown that during heat shock a sealed package can reach pressures greater than the ice‐CO₂ hydrate equilibrium pressure. At pressures above the ice‐CO₂ hydrate equilibrium pressure, the fraction of water crystallized in the dessert can be increased, potentially mitigating heat shock damage. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Phase and kinetic behavior of the mixed methane and carbon dioxide hydrates

Large amounts of CH₄ are stored as hydrates on continental margins and permafrost regions. If the CH₄ hydrates could be converted into CO₂ hydrate, they would serve double duty as CH₄ sources and CO₂ storage sites in the deep ocean sediments. As preliminary investigations, both the phase behavior of CH₄ hydrates and kinetic behavior of CO₂ hydrate were measured at versatile conditions that can simulate actual marine sediments. When measuring three-phase equilibria (H-LW-V) containing CH₄ hydrate, we also closely examined pore and electrolyte effects of clay and NaCl on hydrate formation. These two effects inhibited hydrate nucleation and thus made the hydrate equilibrium line shift to a higher pressure region. In addition, the kinetic data of CO₂ hydrate in the mixtures containing clay and NaCl were determined at 2.0 MPa and 274.15 K. Clay mineral accelerated an initial formation rate of CO₂ hydrate by inducing nucleation as initiator, but total amount of formed CO₂, of course, decreased due to the capillary effect of clay pores. Also, the addition of NaCl in sample mixtures made both initial formation rate and total amount of CO₂ consumption decrease.     

Kinetic simulation of methane hydrate formation and dissociation in porous media

The vast amount of hydrocarbon gas deposited in the earth's crust as gas hydrates has significant implications for future energy supply and global climate. A 3-D simulator for methane hydrate formation and dissociation in porous media is developed for designing and interpreting laboratory and field hydrate experiments. Four components (hydrate, methane, water and salt) and five phases (hydrate, gas, aqueous-phase, ice and salt precipitate) are considered in the simulator. The intrinsic kinetics of hydrate formation or dissociation is considered using the Kim-Bishnoi model. Water freezing and ice melting are tracked with primary variable switch method (PVSM) by assuming equilibrium phase transition. Mass transport, including two-phase flow and molecular diffusions, and heat transfer involved in formation or dissociation of hydrates are included in the governing equations, which are discretized with finite volume difference method and are solved in a fully implicit manner. The developed simulator is used here to study the formation and the dissociation of hydrates in laboratory-scale core samples. In hydrate formation from the system of gas and ice (G+I) and in hydrate dissociation systems where ice appears, the equilibrium between aqueous-phase and ice (A-I) is found to have a “blocking” effect on heat transfer when salt is absent from the system. Increase of initial temperature (at constant outlet pressure), introduction of salt component into the system, decrease of outlet pressure, and increase of boundary heat transfer coefficient can lead to faster hydrate dissociation.     

Experimental determination of methane hydrate dissociation curve up to 55 MPa by using a small amount of surfactant as hydrate promoter

Methane hydrate equilibrium has been studied upon continuous heating of the water-hydrate-gas system within the temperature range of 275-300 K. This temperature range corresponds to equilibrium pressures of 3.15-55 MPa. The hydrate formation/dissociation experiments were carried out in a high-pressure reactor under isochoric conditions and with no agitation. A small amount of surfactant (0.02 wt% sodium dodecyl sulfate, SDS) was added to water to promote hydrate formation. It was demonstrated that SDS did not have any influence on the gas hydrate equilibrium, but increased drastically both the hydrate formation rate and the amount of water converted into hydrate, when compared with the experiments without surfactant. To understand and clarify the influence of SDS on hydrate formation, macroscopic observations of hydrate growth were carried out using gas propane as hydrate former in a fully transparent reactor. We observed that 10^-3 wt% SDS (230 times less than the Critical Micellar Concentration of SDS) were sufficient to prevent hydrate particles from agglomerating and forming a rigid hydrate film at the liquid-gas interface. In the presence of SDS, hydrates grew mainly on the reactor walls as a porous structure, which sucked the solution due to capillary forces. Hydrates grew with a high rate until about 97 wt% of the water present in the reactor was transformed into hydrate.Our data on methane hydrate equilibrium both confirm already published literature data and complement them within the pressure range of 20-55 MPa.     

Metastable state of gas hydrate during decomposition: A novel phenomenon

Natural gas hydrates are solid compounds with cage-like structures formed by gas and water. An intriguing phenomenon that gas hydrates can dissociate at a low rate below the ice freezing point has been viewed as the metastability of hydrate. The mechanisms of hydrate metastability have been widely studied, and many mechanisms were proposed involving the self-preservation effect, supercooled water-gas-hydrate metastable equilibrium, and supersaturated liquid–gas-hydrate system etc. The metastable state of hydrate could be of crucial significance in the kinetics of hydrate formation and decomposition, heat and mass transfer during gas production processes, and the application of hydrate-based technique involving desalination, energy storage and transportation, and gas separation and sequestration. Few researches have systematically considered this phenomenon, and its mechanism remains unclear.In this work, various mechanisms and hypothesis explaining the metastable state of gas hydrates were introduced and discussed. Further studies are still required to reveal the intrinsic nature of this metastable state of gas hydrate,and this work could give some implications on the existing theory and current status of relevant efforts.     

An investigation into the laboratory method for the evaluation of the performance of kinetic hydrate inhibitors using superheated gas hydrates

Kinetic hydrate inhibitors (KHIs) are used to prevent gas hydrate formation in gas and oilfield operations. Recently, a new KHI test method was reported in which hydrates are formed and re-melted just above the equilibrium temperature, before the fluids are re-cooled and the performance of the chemical as a KHI is determined. The method, which we have called the superheated hydrate test method, is claimed to be more reliable for KHI ranking in small equipment, giving less scattering in the hold time data due to avoiding the stochastic nature of the first hydrate formation. We have independently investigated this superheated hydrate test method in steel and sapphire autoclave tests using a gas mixture forming Structure II hydrates and a liquid hydrocarbon phase, which was necessary for satisfactory results. Our results indicate that hold times are shorter than using non-superheated hydrate test methods, but they are more reproducible with less scattering. The reduced scattering occurs in isothermal or slow ramping experiments even when the hydrates are melted at more than 10 °C above the equilibrium temperature (T_eq). However, if a rapid cooling method is used, the improved reproducibility is retained when melting hydrate at 2.4 °C above T_eq but lost when warming to 8.4 °C above T_eq. Using the ramping test method, most, but not all the KHIs tested agreed with the same performance ranking obtained using traditional non-superheated hydrate test methods. This may be related to the variation in the dissociation temperature of gas hydrates with different KHIs and different KHI inhibition mechanisms. Results also varied between different size autoclave equipments.     

Density and phase equilibrium for ice and structure I hydrates using the Gibbs–Helmholtz constrained equation of state

A new, rigorous framework centered around the multi-scale GHC equation of state is presented for predicting bulk density and phase equilibrium for light gas–water mixtures at conditions where hexagonal ice and structure I hydrate phases can exist. The novel aspects of this new framework include (1) the use of internal energies of departure for ice and empty hydrate respectively to determine densities, (2) contributions to the standard state fugacity of water in ice and empty hydrate from lattice structure, (3) computation of these structural contributions to standard state fugacity from compressibility factors and EOS parameters alone, and (4) the direct calculation of gas occupancy from phase equilibrium. Numerical results for densities and equilibrium for systems involving ice and/or gas hydrates predicted by this GHC-based framework are compared to predictions of other equations of state, density correlations, and experimental data where available. Results show that this new GHC-based EOS framework accurately predicts the densities of hexagonal water ice and structure I gas hydrates as well as phase equilibrium for methane–water and CO₂–water mixtures.     

Prediction of carbon dioxide gas hydrate formation conditions in aqueous electrolyte solutions

A methodology for predicting the incipient equilibrium conditions for carbon dioxide gas hydrates in the presence of electrolytes such as NaCl, KCl and CaCl₂ is presented. The method utilizes the statistical thermodynamics model of van der Waals and Platteeuw (1959) to describe the solid hydrate phase. Three different models were examined for the representation of the liquid phase: Chen and Evans (1986), Zuo and Guo (1991), and Aasberg-Petersen et al. (1991). It was found that the model of Zuo and Guo (1991) gave the best results for predicting incipient CO₂ gas hydrate conditions in aqueous single salt solutions. The model was then extended for prediction of CO₂ gas hydrates in mixed salts solutions. The predictions agree very well with experimental data.     

Dynamic measurements of hydrate based gas separation in cooled silica gel

Hydrate based gas separation is a promising method for carbon dioxide capture. The purpose of this study is to analyze hydrates formation and dissociation characters when gas mixture flows through cooled silica gel. The additives mixture (THF/SDS) was used to saturate the silica gel partly, and gas mixture (CO₂/H₂) was injected into it to form hydrates. Magnetic resonance imaging (MRI) images were obtained using fast spin echo multi-slice pulse sequence. Hydrates saturations were calculated quantitatively using MRI data. The experimental results showed that the optimal initial solution saturation was 34.2% in this investigation. The gas component was analyzed to assess the separation efficiency. For hydrates dissociation processes at 1 atmospheric pressure, CO₂ concentrations increased obviously. Half of the six cycles showed that more than 85.00 mol% CO₂ contained in the capture gas, and the lowest CO₂ concentration was 64.83 mol%. Hydrate blockages appeared frequently, which restricted the contact of gas and solution and caused the incomplete transformations of residual solution to hydrates. It was a key restricted factor for hydrate based CO₂ capture.     

Prediction of hydrate formation conditions based on the vdWP‐type models at high pressures

Prediction of phase boundaries of gas hydrates has been done for several decades based on the vdWP (van der Waals and Platteeuw) hydrate equation and the classical thermodynamic equations for describing the water fugacities in water or ice phase. This procedure gives a reasonable prediction at low pressures, but when the pressure increases, above 10⁵ kPa, it shows a significant error. In the conventional vdWP‐type models it has been assumed that the volume difference between the empty hydrate lattice and pure liquid water is independent of the system pressure and temperature. In this work, different approaches for describing the volume dependency of pure liquid water and the empty hydrate lattice on the system pressure have been used to predict the hydrate equilibria based on the vdWP‐type model. Also, an expression is introduced to estimate the volume of methane hydrate lattice as a function of pressure and temperature. Finally, this method is extended to other hydrate formers, that is, ethane, carbon dioxide, xenon, and nitrogen. The predicted values are in good agreement with the experimental data both for L_w–H–V and L_w–H–L_hf phase boundaries.