Tuning methane content in gas hydrates via thermodynamic modeling and molecular dynamics simulation

Storage and transportation of natural gas as gas hydrate (“gas-to-solids technology”) is a promising alternative to the established liquefied natural gas (LNG) or compressed natural gas (CNG) technologies. Gas hydrates offer a relatively high gas storage capacity and mild temperature and pressure conditions for formation. Simulations based on the van der Waals–Platteeuw model and molecular dynamics (MD) are employed in this study to relate the methane gas content/occupancy in different hydrate systems with the hydrate stability conditions including temperature, pressure, and secondary clathrate stabilizing guests. Methane is chosen as a model system for natural gas. It was found that the addition of about 1% propane suffices to increase the structure II (sII) methane hydrate stability without excessively compromising methane storage capacity in hydrate. When tetrahydrofuran (THF) is used as the stabilizing agent in sII hydrate at concentration between 1% and 3%, a reasonably high methane content in hydrate can be maintained (∼85–100, v/v) without dealing with pressures more than 5 MPa and close to room temperature.     

Spectroscopic methods in gas hydrate research

Gas hydrates are crystalline structures comprising a guest molecule surrounded by a water cage, and are particularly relevant due to their natural occurrence in the deep sea and in permafrost areas. Low molecular weight molecules such as methane and carbon dioxide can be sequestered into that cage at suitable temperatures and pressures, facilitating the transition to the solid phase. While the composition and structure of gas hydrates appear to be well understood, their formation and dissociation mechanisms, along with the dynamics and kinetics associated with those processes, remain ambiguous. In order to take advantage of gas hydrates as an energy resource (e.g., methane hydrate), as a sequestration matrix in (for example) CO₂ storage, or for chemical energy conservation/storage, a more detailed molecular level understanding of their formation and dissociation processes, as well as the chemical, physical, and biological parameters that affect these processes, is required. Spectroscopic techniques appear to be most suitable for analyzing the structures of gas hydrates (sometimes in situ), thus providing access to such information across the electromagnetic spectrum. A variety of spectroscopic methods are currently used in gas hydrate research to determine the composition, structure, cage occupancy, guest molecule position, and binding/formation/dissociation mechanisms of the hydrate. To date, the most commonly applied techniques are Raman spectroscopy and solid-state nuclear magnetic resonance (NMR) spectroscopy. Diffraction methods such as neutron and X-ray diffraction are used to determine gas hydrate structures, and to study lattice expansions. Furthermore, UV-vis spectroscopic techniques and scanning electron microscopy (SEM) have assisted in structural studies of gas hydrates. Most recently, waveguide-coupled mid-infrared spectroscopy in the 3–20 μm spectral range has demonstrated its value for in situ studies on the formation and dissociation of gas hydrates. This comprehensive review summarizes the importance of spectroscopic analytical techniques to our understanding of the structure and dynamics of gas hydrate systems, and highlights selected examples that illustrate the utility of these individual methods.     

Investigation on Gas Storage in Methane Hydrate

The effect of additives (anionic surfactant sodium dodecyl sulfate (SDS), nonionic surfactantalkyl polysaccharide glycoside (APG), and liquid hydrocarbon cyclopentane (CP)) on hydrate inductiontime and formation rate, and storage capacity was studied in this work. Micelle surfactant solutions werefound to reduce hydrate induction time, increase methane hydrate formation rate and improve methanestorage capacity in hydrates. In the presence of surfactant, hydrate could form quickly in a quiescentsystem and the energy costs of hydrate formation were reduced. The critical micelle concentrations of SDS and APG water solutions were found to be 300x 10-6 and 500x 10-6 for methane hydrate formation systemrespectively. The effect of anionic surfactant (SDS) on methane storage in hydrates is more pronounced compared to a nonionic surfactant (APG). CP also reduced hydrate induction time and improved hydrateformation rate, but could not improve methane storage in hydrates.     

Hydrate phase equilibrium and structure for (methane + ethane + tetrahydrofuran + water) system

The separation of methane and ethane through forming hydrate is a possible choice in natural gas, oil processing, or ethylene producing. The hydrate formation conditions of five groups of (methane + ethane) binary gas mixtures in the presence of 0.06 mole fraction tetrahydrofuran (THF) in water were obtained at temperatures ranging from (277.7 to 288.2) K. In most cases, the presence of THF in water can lower the hydrate formation pressure of (methane + ethane) remarkably. However, when the composition of ethane is as high as 0.832, it is more difficult to form hydrate than without THF system. Phase equilibrium model for hydrates containing THF was developed based on a two-step hydrate formation mechanism. The structure of hydrates formed from (methane + ethane + THF + water) system was also determined by Raman spectroscopy. When THF concentration in initial aqueous solution was only 0.06 mole fraction, the coexistence of structure I hydrate dominated by ethane and structure II hydrate dominated by THF in the hydrate sample was clearly demonstrated by Raman spectroscopic data. On the contrary, only structure II hydrate existed in the hydrate sample formed from (methane + ethane + THF + water) system when THF concentration in initial aqueous solution was increased to 0.10 mole fraction. It indicated that higher THF concentration inhibited the formation of structure I hydrate dominated by ethane and therefore lowered the trapping of ethane in hydrate. It implies a very promising method to increase the separation efficiency of methane and ethane.     

Study on thermal properties of TBAB–THF hydrate mixture for cold storage by DSC

In order to study the thermal properties of new type environment-friendly binary hydrate for cold storage in air-conditioning system, tests have been carried out by DSC comprehensively on the phase-change temperature and fusion heat of TBAB hydrate, THF hydrate, and TBAB–THF hydrate mixture. The results show a good trend that TBAB–THF hydrate has the superiority for more proper phase-change temperature and increased fusion heat. A broader and more developed view is that adding appropriate amount of hydrate with lower phase-change temperature to hydrate with higher one can make the hydrate mixture more suitable for cold storage (especially for 278–281 K); some hydrates with lower phase-change temperature can even make the fusion heat of mixture hydrate increased greatly. Several new environmental working pairs for binary gas hydrates have been listed to help to promote the application.     

常见客体分子对笼型水合物晶格常数的影响

Natural gas hydrates are considered as ideal alternative energy resources for the future, and the relevant basic and applied research has become more attractive in recent years. The influence of guest molecules on the hydrate crystal lattice parameters is of great significances to the understanding of hydrate structural characteristics, hydrate formation/decomposition mechanisms, and phase stability behaviors. In this study, we test a series of artificial hydrate samples containing different guest molecules (e.g. methane, ethane, propane, iso-butane, carbon dioxide, tetrahydrofuran, methane + 2, 2-dimethylbutane, and methane + methyl cyclohexane) by a low-temperature powder X-ray diffraction (PXRD). Results show that PXRD effectively elucidates structural characteristics of the natural gas hydrate samples, including crystal lattice parameters and structure types. The relationships between guest molecule sizes and crystal lattice parameters reveal that different guest molecules have different controlling behaviors on the hydrate types and crystal lattice constants. First, a positive correlation between the lattice constants and the van der Waals diameters of homologous hydrocarbon gases was observed in the single-guest-component hydrates. Small hydrocarbon homologous gases, such as methane and ethane, tended to form sI hydrates, whereas relatively larger molecules, such as propane and iso-butane, generated sⅡ hydrates. The hydrate crystal lattice constants increased with increasing guest molecule size. The types of hydrates composed of oxygen-containing guest molecules (such as CO₂ and THF) were also controlled by the van der Waals diameters. However, no positive correlation between the lattice constants and the van der Waals diameters of guest molecules in hydrocarbon hydrates was observed for CO₂ hydrate and THF hydrate, probably due to the special interactions between the guest oxygen atoms and hydrate "cages". Furthermore, the influences of the macromolecules and auxiliary small molecules on the lengths of the different crystal axes of the sH hydrates showed inverse trends. Compared to the methane + 2, 2-dimethylbutane hydrate sample, the length of the a-axis direction of the methane + methyl cyclohexane hydrate sample was slightly smaller, whereas the length of the c-axis direction was slightly longer. The crystal a-axis length of the sH hydrate sample formed with nitrogen molecules was slightly longer, whereas the c-axis was shorter than that of the methane + 2, 2-dimethylbutane hydrate sample at the same temperature.     

Phase Transition of a Structure II Cubic Clathrate Hydrate to a Tetragonal Form

The crystal structure and phase transition of cubic structure II (sII) binary clathrate hydrates of methane (CH₄) and propanol are reported from powder X‐ray diffraction measurements. The deformation of host water cages at the cubic–tetragonal phase transition of 2‐propanol+CH₄ hydrate, but not 1‐propanol+CH₄ hydrate, was observed below about 110 K. It is shown that the deformation of the host water cages of 2‐propanol+CH₄ hydrate can be explained by the restriction of the motion of 2‐propanol within the 5¹²6⁴ host water cages. This result provides a low‐temperature structure due to a temperature‐induced symmetry‐lowering transition of clathrate hydrate. This is the first example of a cubic structure of the common clathrate hydrate families at a fixed composition.     

The effect of the water/methane interface on methane hydrate cages: the potential of mean force and cage lifetimes

Molecular dynamics simulations were used to determine the influence of a methane-water interface on the position and stability of methane hydrate cages. A potential of mean force was calculated as a function of the separation of a methane hydrate cage and a methane-water interface. The hydrate cages are found to be strongly repelled from the methane gas into the water phase. At low enough temperatures, however, the most favorable location for the hydrate cage is at the interface on the water side. Cage lifetime simulations were performed in bulk water and near a methane-water interface. The methane-water interface increases the cage lifetime by almost a factor of 2 compared to cage lifetimes of cages in bulk water. The potential of mean force and the cage lifetime results give additional explanations for the proposed nucleation of gas hydrates at gas-water interfaces.     

Molecular dynamics study of structure H clathrate hydrates of methane and large guest molecules

Methane storage in structure H (sH) clathrate hydrates is attractive due to the relatively higher stability of sH as compared to structure I methane hydrate. The additional stability is gained without losing a significant amount of gas storage density as happens in the case of structure II (sII) methane clathrate. Our previous work has showed that the selection of a specific large molecule guest substance (LMGS) as the sH hydrate former is critical in obtaining the optimum conditions for crystallization kinetics, hydrate stability, and methane content. In this work, molecular dynamics simulations are employed to provide further insight regarding the dependence of methane occupancy on the type of the LMGS and pressure. Moreover, the preference of methane molecules to occupy the small (5(12)) or medium (4(3)5(6)6(3)) cages and the minimum cage occupancy required to maintain sH clathrate mechanical stability are examined. We found that thermodynamically, methane occupancy depends on pressure but not on the nature of the LMGS. The experimentally observed differences in methane occupancy for different LMGS may be attributed to the differences in crystallization kinetics and/or the nonequilibrium conditions during the formation. It is also predicted that full methane occupancies in both small and medium clathrate cages are preferred at higher pressures but these cages are not fully occupied at lower pressures. It was found that both small and medium cages are equally favored for occupancy by methane guests and at the same methane content, the system suffers a free energy penalty if only one type of cage is occupied. The simulations confirm the instability of the hydrate when the small and medium cages are empty. Hydrate decomposition was observed when less than 40% of the small and medium cages are occupied.     

用分子动力学模拟水合物储氢

用分子动力学(MD)模拟方法系统研究了结构Ⅰ型(SⅠ)和结构H型(SH)氢气水合物中氢气的占据情况并确定了氢气水合物的稳定结构: SⅠ水合物氢气分子数小胞中为2, 大胞中为3; SH水合物氢气分子数小胞中为2, 中胞中为2, 大胞中为11. 分析了稳定情况下水合物各胞腔内氢气分子之间的径向分布函数(RDF), 得出了氢气分子在各胞腔内的稳定位置. 由稳定位置得到了稳定结构下氢气水合物的储氢质量分数: SⅠ为5.085%, SH为6.467%. 与实验对比得出结论: SH水合物稳定结构下的储氢能力最强.     

Compressibility of Gas Hydrates

Experimental data on the pressure dependence of unit cell parameters for the gas hydrates of ethane (cubic structure I, pressure range 0–2 GPa), xenon (cubic structure I, pressure range 0–1.5 GPa) and the double hydrate of tetrahydrofuran+xenon (cubic structure II, pressure range 0–3 GPa) are presented. Approximation of the data using the cubic Birch–Murnaghan equation, P=1.5B₀[(V₀/V)^7/3?(V₀/V)^5/3], gave the following results: for ethane hydrate V₀=1781 ų, B₀=11.2 GPa; for xenon hydrate V₀=1726 ų, B₀=9.3 GPa; for the double hydrate of tetrahydrofuran+xenon V₀=5323 ų, B₀=8.8 GPa. In the last case, the approximation was performed within the pressure range 0–1.5 GPa; it is impossible to describe the results within a broader pressure range using the cubic Birch–Murnaghan equation. At the maximum pressure of the existence of the double hydrate of tetrahydrofuran+xenon (3.1 GPa), the unit cell volume was 86 % of the unit cell volume at zero pressure. Analysis of the experimental data obtained by us and data available from the literature showed that 1) the bulk modulus of gas hydrates with classical polyhedral structures, in most cases, are close to each other and 2) the bulk modulus is mainly determined by the elasticity of the hydrogen‐bonded water framework. Variable filling of the cavities with guest molecules also has a substantial effect on the bulk modulus. On the basis of the obtained results, we concluded that the bulk modulus of gas hydrates with classical polyhedral structures and existing at pressures up to 1.5 GPa was equal to (9±2) GPa. In cases when data on the equations of state for the hydrates were unavailable, the indicated values may be recommended as the most probable ones.     

High‐Performance of Gas Hydrates in Confined Nanospace for Reversible CH4/CO2 Storage

The molecular exchange of CH₄ for CO₂ in gas hydrates grown in confined nanospace has been evaluated for the first time using activated carbons as a host structure. The nano‐confinement effects taking place inside the carbon cavities and the exceptional physicochemical properties of the carbon structure allows us to accelerate the formation and decomposition process of the gas hydrates from the conventional timescale of hours/days in artificial bulk systems to minutes in confined nanospace. The CH₄/CO₂ exchange process is fully reversible with high efficiency at practical temperature and pressure conditions. Furthermore, these activated carbons can be envisaged as promising materials for long‐distance natural gas and CO₂ transportation because of the combination of a high storage capacity, a high reversibility, and most important, with extremely fast kinetics for gas hydrate formation and release.     

Kinetics and stability of CH4-CO2 mixed gas hydrates during formation and long-term storage.

The formation of CH4-CO2 mixed gas hydrates was observed by measuring the change of vapor-phase composition using gas chromatography and Raman spectroscopy. Preferential consumption of carbon dioxide molecules was found during hydrate formation, which agreed well with thermodynamic calculations. Both Raman spectroscopic analysis and the thermodynamic calculation indicated that the kinetics of this mixed gas hydrate system was controlled by the competition of both molecules to be enclathrated into the hydrate cages. However, the methane molecules were preferentially crystallized in the early stages of hydrate formation when the initial methane concentration was much less than that of carbon dioxide. According to the Roman spectra, pure methane hydrates first formed under this condition. This unique phenomenon suggested that methane molecules play important roles in the hydrate formation process. These mixed gas hydrates were stored at atmospheric pressure and 190 K for over two months to examine the stability of the encaged gases. During storage, CO2 was preferentially released. According to our thermodynamic analysis, this CO2 release was due to the instability of CO2 in the hydrate structure under the storage conditions.     

Effect of Guest–Host Hydrogen Bonding on the Structures and Properties of Clathrate Hydrates

To provide improved understanding of guest–host interactions in clathrate hydrates, we present some correlations between guest chemical structures and observations on the corresponding hydrate properties. From these correlations it is clear that directional interactions such as hydrogen bonding between guest and host are likely, although these have been ignored to greater or lesser degrees because there has been no direct structural evidence for such interactions. For the first time, single‐crystal X‐ray crystallography has been used to detect guest–host hydrogen bonding in structure II (sII) and structure H (sH) clathrate hydrates. The clathrates studied are the tert‐butylamine (tBA) sII clathrate with H₂S/Xe help gases and the pinacolone + H₂S binary sH clathrate. X‐ray structural analysis shows that the tBA nitrogen atom lies at a distance of 2.64 Å from the closest clathrate hydrate water oxygen atom, whereas the pinacolone oxygen atom is determined to lie at a distance of 2.96 Å from the closest water oxygen atom. These distances are compatible with guest–water hydrogen bonding. Results of molecular dynamics simulations on these systems are consistent with the X‐ray crystallographic observations. The tBA guest shows long‐lived guest–host hydrogen bonding with the nitrogen atom tethered to a water HO group that rotates towards the cage center to face the guest nitrogen atom. Pinacolone forms thermally activated guest–host hydrogen bonds with the lattice water molecules; these have been studied for temperatures in the range of 100–250 K. Guest–host hydrogen bonding leads to the formation of Bjerrum L‐defects in the clathrate water lattice between two adjacent water molecules, and these are implicated in the stabilities of the hydrate lattices, the water dynamics, and the dielectric properties. The reported stable hydrogen‐bonded guest–host structures also tend to blur the longstanding distinction between true clathrates and semiclathrates.     

Gas hydrate fast nucleation from melting ice and quiescent growth along vertical heat transfer tube

Gas hydrates, or clathrate hydrates, are ice-likecrystal, composed of host lattice (cavities) formed byhydrogen-bonded water molecules, and other guestmolecules called guest molecules. The guest mole-cules act with host lattice in weak van der Waals force…     

On the thermodynamic stability and structural transition of clathrate hydrates

Gas mixtures of methane and ethane form structure II clathrate hydrates despite the fact that each of pure methane and pure ethane gases forms the structure I hydrate. Optimization of the interaction potential parameters for methane and ethane is attempted so as to reproduce the dissociation pressures of each simple hydrate containing either methane or ethane alone. An account for the structural transitions between type I and type II hydrates upon changing the mole fraction of the gas mixture is given on the basis of the van der Waals and Platteeuw theory with these optimized potentials. Cage occupancies of the two kinds of hydrates are also calculated as functions of the mole fraction at the dissociation pressure and at a fixed pressure well above the dissociation pressure.     

Intensification of methane and hydrogen storage in clathrate hydrate and future prospect

Gas hydrate is a new technology for energy gas (methane/hydrogen) storage due to its large capacity of gas storage and safe. But industrial application of hydrate storage process was hindered by some problems. For methane, the main problems are low formation rate and storage capacity, which can be solved by strengthening mass and heat transfer, such as adding additives, stirring, bubbling, etc. One kind of additives can change the equilibrium curve to reduce the formation pressure of methane hydrate, and the other kind of additives is surfactant, which can form micelle with water and increase the interface of water-gas. Dry water has the similar effects on the methane hydrate as surfactant. Additionally, stirring, bubbling, and spraying can increase formation rate and storage capacity due to mass transfer strengthened. Inserting internal or external heat exchange also can improve formation rate because of good heat transfer. For hydrogen, the main difficulties are very high pressure for hydrate formed. Tetrahydrofuran (THF), tetrabutylammonium bromide (TBAB) and tetrabutylammonium fluoride (TBAF) have been proved to be able to decrease the hydrogen hydrate formation pressure significantly.     

Gas hydrate nucleation and cage formation at a water/methane interface

Nucleation of gas hydrates remains a poorly understood phenomenon, despite its importance as a critical step in understanding the performance and mode of action of low dosage hydrate inhibitors. We present here a detailed analysis of the structural and mechanistic processes by which gas hydrates nucleate in a molecular dynamics simulation of dissolved methane at a methane/water interface. It was found that hydrate initially nucleates into a phase consistent with none of the common bulk crystal structures, but containing structural units of all of them. The process of water cage formation has been found to correlate strongly with the collective arrangement of methane molecules.     

Antifreezes Act as Catalysts for Methane Hydrate Formation from Ice

Contrary to the thermodynamic inhibiting effect of methanol on methane hydrate formation from aqueous phases, hydrate forms quickly at high yield by exposing frozen water–methanol mixtures with methanol concentrations ranging from 0.6–10 wt % to methane gas at pressures from 125 bars at 253 K. Formation rates are some two orders of magnitude greater than those obtained for samples without methanol and conversion of ice is essentially complete. Ammonia has a similar catalytic effect when used in concentrations of 0.3–2.7 wt %. The structure I methane hydrate formed in this manner was characterized by powder X‐ray diffraction and Raman spectroscopy. Steps in the possible mechanism of action of methanol were studied with molecular dynamics simulations of the Ih (0001) basal plane exposed to methanol and methane gas. Simulations show that methanol from a surface aqueous layer slowly migrates into the ice lattice. Methane gas is preferentially adsorbed into the aqueous methanol surface layer. Possible consequences of the catalytic methane hydrate formation on hydrate plug formation in gas pipelines, on large scale energy‐efficient gas hydrate formation, and in planetary science are discussed.     

Micro-Raman investigations of mixed gas hydrates

We report laser Raman spectroscopic measurements on mixed hydrates (clathrates), with guest molecules tetrahydrofuran (THF) and methane (CH₄), at ambient pressure and at temperatures from 175 to 280 K. Gas hydrates were synthesized with different concentrations of THF ranging from 5.88 to 1.46 mol%. In all cases THF molecules occupied the large cages of sII hydrate. The present studies demonstrate formation of sII clathrates with CH₄ molecules occupying unfilled cages for concentrations of THF ranging from 5.88 to 2.95 mol%. The Raman spectral signature of hydrates with 1.46 mol% THF are distinctly different; hydrate growth was non-uniform and structural transformation occurred from sII to sI prior to clathrate melting.