Sorption of methane and CO2 for enhanced coalbed methane recovery and carbon dioxide sequestration

Sequestration of CO2 in deep and unmineable coal seams is one of the attractive alternatives to reduce its atmospheric concentration. Injection of CO2 in coal seams may help in enhancing the recovery of coalbed methane. An experimental study has been carried out using coal samples from three different coal seams, to evaluate the enhanced gas recovery and sequestration potential of these coals. The coals were first saturated with methane and then by depressurization some of the adsorbed methane was desorbed. After partial desorption, CO2 was injected into the coals and subsequently they were depressurized again. Desorption of methane after the injections was studied, to investigate the ability of CO2 to displace and enhance the recovery of methane from the coals. The coals exhibited varying behavior of adsorption of CO2 and release of methane. For one coal, the release of methane was enhanced by injection of CO2, suggesting preferential adsorption of CO2 and desorption of methane. For the other two coals, CO2 injection did not produce incremental methane initially, as there was initial resistance to methane release. However with continued CO2 injection, most of the remaining methane was produced. The study suggested that preferential sorption behavior of coal and enhanced gas recovery pattern could not be generalized for all coals.     

Structural and dynamical properties of methane clathrate hydrates

Equilibrium molecular dynamics (MD) simulations have been performed in both the NVT and NPT ensembles to study the structural and dynamical properties of fully occupied methane clathrate hydrates at 50, 125, and 200 K. Five atomistic potential models were used for water, ranging from fully flexible to rigid polarizable and nonpolarizable. A flexible and a rigid model were utilized for methane. The phonon densities of states were evaluated and the localized rattling modes for the methane molecules were found to couple to the acoustic phonons of the host lattice. The calculated methane density of states was found to be in reasonable agreement with available experimental data.     

Equilibrium conditions of clathrate hydrates formed from carbon dioxide and aqueous acetone solutions

Equilibrium conditions of clathrate hydrates formed from carbon dioxide and aqueous acetone solutions were experimentally measured at temperatures between 269.2 and 281.4 K and pressures up to 3.98 MPa. The acetone concentrations in solutions were investigated from 0.04 to 0.40 mass fractions. The experimental results suggested a transition in hydrate structure from structure I to another structure for acetone solutions between 0.04 and 0.12 mass fractions of acetone. The hydrate structure was suggested to be structure II which was the most stable with a 0.16 mass fraction acetone solution. For more than 0.16 mass fraction of acetone, the equilibrium conditions of the hydrate were shifted to lower temperatures as acetone concentrations increased.     

Vibrational spectra of methane clathrate hydrates from molecular dynamics simulation

Molecular dynamics simulations were performed on methane clathrate hydrates at ambient conditions. Thermal expansion results over the temperature range 60-300 K show that the unit cell volume increases with temperature in agreement with experiment. Power spectra were obtained at 273 K from velocity autocorrelation functions for selected atoms, and normal modes were assigned. The spectra were further classified according to individual atom types, allowing the assignment of contributions from methane molecules located in small and large cages within the structure I unit cell. The symmetric C-H stretch of methane in the small cages occurs at a higher frequency than for methane located in the large cages, with a peak separation of 14 cm(-1). Additionally, we determined that the symmetric C-H stretch in methane gas occurs at the same frequency as methane in the large cages. Results of molecular dynamics simulations indicate the use of power spectra obtained from the velocity autocorrelation function is a reliable method to investigate the vibrational behavior of guest molecules in clathrate hydrates.     

Nuclear magnetic resonance parameters for methane molecule trapped in clathrate hydrates

The calculations of the nuclear shielding and spin-spin coupling constants were carried out for two models of clathrate hydrates, 5(12) and 5(12)6(8), using the density functional theory three-parameter Becke-Lee-Yang-Parr method with the basis set aug-cc-pVDZ (optimization) and HuzIII-su3 (NMR parameters). Particular attention has been devoted to evaluate the influence of a geometrical arrangement, the effect of long-range interactions on the NMR shielding of methane molecule, and to predict whether (13)C and (1)H chemical shifts can distinguish between guests in two clathrate hydrates cages. The correlation of the changes in the (17)O shielding constants depend strongly on the hydrogen-bonding topology. The intermolecular hydrogen-bond transmitted (1h)J(OH) spin-spin coupling constants are substantial. The increase of their values is connected with the elongation of the intramolecular O-H bond and the shortening of the intermolecular O···H distance. These data suggests that hydrogen bonds between double donor-single acceptor (DDA)-type water molecules acting as a proton acceptor from single donor-double acceptor (DAA)-type water molecules are stronger than ones formed by DAA-type water molecules acting as an acceptor for a DDA water proton. These state-of-the-art calculations confirmed the earlier experimental findings of the cage-dependency of (13)C chemical shift of methane.     

The cages, dynamics, and structuring of incipient methane clathrate hydrates

Interest in describing clathrate hydrate formation mechanisms spans multiple fields of science and technical applications. Here, we report findings from multiple molecular dynamics simulations of spontaneous methane clathrate hydrate nucleation and growth from fully demixed and disordered two-phase fluid systems of methane and water. Across a range of thermodynamic conditions and simulation geometries and sizes, a set of seven cage types comprises approximately 95% of all cages formed in the nucleated solids. This set includes the ubiquitous 5(12) cage, the 5(12)6(n) subset (where n ranges from 2-4), and the 4(1)5(10)6(n) subset (where n also ranges from 2-4). Transformations among these cages occur via water pair insertions/removals and rotations, and may elucidate the mechanisms of solid-solid structural rearrangements observed experimentally. Some consistency is observed in the relative abundance of cages among all nucleation trajectories. 5(12) cages are always among the two most abundant cage types in the nucleated solids and are usually the most abundant cage type. In all simulations, the 5(12)6(n) cages outnumber their 4(1)5(10)6(n) counterparts with the same number of water molecules. Within these consistent features, some stochasticity is observed in certain cage ratios and in the long-range ordering of the nucleated solids. Even when comparing simulations performed at the same conditions, some trajectories yield swaths of multiple adjacent sI unit cells and long-range order over 5 nm, while others yield only isolated sI unit cells and little long-range order. The nucleated solids containing long-range order have higher 5(12)6(2)/5(12) and 5(12)6(3)/4(1)5(10)6(2) cage ratios when compared to systems that nucleate with little long-range order. The formation of multiple adjacent unit cells of sI hydrate at high driving forces suggests an alternative or addition to the prevailing hydrate nucleation hypotheses which involve formation through amorphous intermediates.     

Equilibrium conditions for clathrate hydrates formed from methane and aqueous propanol solutions

Equilibrium conditions for clathrate hydrates formed from methane and different concentrations of 1-propanol or 2-propanol aqueous solutions were experimentally determined at temperatures of 274.0–287.1 K and pressures up to 11.0 MPa. Each propanol has an inhibiting and/or promoting effect on hydrate formation depending on the propanol concentration. A structural transition from a structure I to a different hydrate structure occurred at concentrations between 3 and 5 mass% for 1-propanol and between 2 and 3 mass% for 2-propanol.     

Molecular dynamics study of the stability of methane structure H clathrate hydrates

Molecular dynamics simulations are used to study the stability of structure H (sH) methane clathrate hydrates in a 3 x 3 x 3 sH unit cell replica. Simulations are performed at experimental conditions of 300 K and 2 GPa for three methane intermolecular potentials. The five small cages of the sH unit cell are assigned methane guest occupancies of one and large cage guest occupancies of one to five are considered. Radial distribution functions, unit cell volumes, and configurational energies are studied as a function of large cage CH(4) occupancy. Free energy calculations are carried out to determine the stability of clathrates for large cage occupancies. Large cage occupancy of five is the most stable configuration for a Lennard-Jones united-atom potential and the Tse-Klein-McDonald potential parametrized for condensed methane phases and two for the most stable configuation for the Murad and Gubbins potential.     

Stability and reactivity of methane clathrate hydrates: insights from density functional theory

The sI methane clathrate hydrate consists of methane gas molecules encapsulated as dodecahedron (5(12)CH(4)) and tetrakaidecahedron (5(12)6(2)CH(4)) water cages. The characterization of the stability of these cages is crucial to an understanding of the mechanism of their formation. In the present work, we perform calculations using density functional theory to calculate interaction energies, free energies, and reactivity indices of these cages. The contributions from polarization functions to interaction energies is more than diffuse functions from Pople basis sets, though both functions from the correlation-consistent basis sets contribute significantly to interaction energies. The interaction energies and free energies show that the formation of the 5(12)CH(4) cage (from the 5(12) cage) is more favored compared to the 5(12)6(2)CH(4) cage (from the 5(12)6(2) cage). The pressure-dependent study shows a spontaneous formation of the 5(12)CH(4) cage at 273 K (P ≥ 77 bar) and the 5(12)6(2)CH(4) cage (P = 100 bar). The reactivity of the 5(12)CH(4) cage is similar to that of the 5(12) cage, but the 5(12)6(2)CH(4) cage is more reactive than the 5(12)6(2) cage.     

Sorption of methane and CO2 for enhanced coalbed methane recovery and carbon dioxide seauestration

Sequestration of CO2 in deep and unmineable coal seams is one of the attractive alternatives to reduce its atmospheric concentration. Injection of CO2 in coal seams may help in enhancing the recovery of coalbed methane. An experimental study has been carried out using coal samples from three different coal seams, to evaluate the enhanced gas recovery and sequestration potential of these coals. The coals were first saturated with methane and then by depressurization some of the adsorbed methane was desorbed. After partial desorption, CO2 was injected into the coals and subsequently they were depressurized again. Desorption of methane after the injections was studied, to investigate the ability of CO2 to displace and enhance the recovery of methane from the coals. The coals exhibited varying behavior of adsorption of CO2 and release of methane. For one coal, the release of methane was enhanced by injection of CO2, suggesting preferential adsorption of CO2 and desorption of methane. For the other two coals, CO2 injection did not produce incremental methane initially, as there was initial resistance to methane release. However with continued CO2 injection, most of the remaining methane was produced. The study suggested that preferential sorption behavior of coal and enhanced gas recovery pattern could not be generalized for all coals.     

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.     

Raman spectra of vibrational and librational modes in methane clathrate hydrates using density functional theory

The sI type methane clathrate hydrate lattice is formed during the process of nucleation where methane gas molecules are encapsulated in the form of dodecahedron (5(12)CH(4)) and tetrakaidecahedron (5(12)6(2)CH(4)) water cages. The characterization of change in the vibrational modes which occur on the encapsulation of CH(4) in these cages plays a key role in understanding the formation of these cages and subsequent growth to form the hydrate lattice. In this present work, we have chosen the density functional theory (DFT) using the dispersion corrected B97-D functional to characterize the Raman frequency vibrational modes of CH(4) and surrounding water molecules in these cages. The symmetric and asymmetric C-H stretch in the 5(12)CH(4) cage is found to shift to higher frequency due to dispersion interaction of the encapsulated CH(4) molecule with the water molecules of the cages. However, the symmetric and asymmetric O-H stretch of water molecules in 5(12)CH(4) and 5(12)6(2)CH(4) cages are shifted towards lower frequency due to hydrogen bonding, and interactions with the encapsulated CH(4) molecules. The CH(4) bending modes in the 5(12)CH(4) and 5(12)6(2)CH(4) cages are blueshifted, though the magnitude of the shifts is lower compared to modes in the high frequency region which suggests bending modes are less affected on encapsulation of CH(4). The low frequency librational modes which are collective motion of the water molecules and CH(4) in these cages show a broad range of frequencies which suggests that these modes largely contribute to the formation of the hydrate lattice.     

Thermodynamic studies of clathrate hydrates

Thermodynamic studies of clathrate hydrates, mainly of structures I and II, are considered in this review which is based on 147 references. There are two main subjects. The first is the host lattice energy and the guest-host interaction energy, both of these quantities being related to the enthalpy of dissociation and composition of the hydrates. The second subject concerns orientational ordering phenomena occurring in both host and guest, as reflected in the low temperature heat capacity. The classical theoretical treatment of clathrate formation has been reconsidered on the basis of recent experimental results. Particular emphasis has been given to orientational ordering since this topic is undoubtedly central to clarifying the nature of clathrate hydrates.Ausgehend von 147 Literaturangaben wurden in diesem Review thermodynamische Untersuchungen von Klathrathydraten hauptsächlich der Struktur I und II betrachtet. Es gibt zwei Hauptaugenmerke. Als erstes die Wirtsgitterenergie und die Gast-Wirt-Wechselwirkungsenergie, beide bezogen auf die Dissoziationsenthalpie und die Bildungsenthalpie der Hydrate. Das zweite Hauptaugenmerk betrifft Orientierungs-Konditionierungserscheinungen sowohl in Wirt als auch Gast, wie in den Wärmekapazitäten bei niedrigen Temperaturen widergespiegelt wird. Auf der Basis jüngster experimenteller Ergebnisse wurde die klassische theoretische Betrachtung über die Bildung von Klathraten überprüft. Der Orientierung-Konditionierung wurde besonderer Nachdruck verliehen, da dies zweifellos eine entscheidende Rolle bei der Klärung der Natur der Klathrathydrate spielt. 147 I II. . «» « — », . «» « », . . , .

---

Contribution No. 155 from the Chemical Thermodynamics Laboratory.     

Phase equilibria of clathrate hydrates of methane + carbon dioxide: New experimental data and predictions

In this communication, we report experimental dissociation conditions for region clathrate hydrates of methane + carbon dioxide in gas–liquid water–hydrate (G–Lw–H) equilibrium. The temperature and pressure conditions are in the range of (279.1–289.9) K and (2.96–13.06) MPa, respectively. Concentrations of carbon dioxide in the feed gas are also varied. An isochoric pressure-search method was used to perform the measurements. The dissociation data generated in this work along with the literature data are compared with the predictions of a thermodynamic model and a previously reported empirical equation. A discussion is made on the deviations between the experimental and predicted data.     

Polymorphism in Br2 clathrate hydrates

The structure and composition of bromine clathrate hydrate has been controversial for more than 170 years due to the large variation of its observed stoichiometries. Several different crystal structures were proposed before 1997 when Udachin et al. (Udachin, K. A.; Enright, G. D.; Ratcliffe, C. I.; Ripmeester, J. A. J. Am. Chem. Soc. 1997, 119, 11481) concluded that Br2 forms only the tetragonal structure (TS-I). We show polymorphism in Br2 clathrate hydrates by identifying two distinct crystal structures through optical microscopy and resonant Raman spectroscopy on single crystals. After growing TS-I crystals from a liquid bromine-water solution, upon dropping the temperature slightly below -7 degrees C, new crystals of cubic morphology form. The new crystals, which have a limited thermal stability range, are assigned to the CS-II structure. The two structures are clearly distinguished by the resonant Raman spectra of the enclathrated Br2, which show long overtone progressions and allow the extraction of accurate vibrational parameters: omega(e) = 321.2 +/- 0.1 cm(-1) and omega(e)x(e) = 0.82 +/- 0.05 cm(-1) in TS-I and omega(e) = 317.5 +/- 0.1 cm(-1) and omega(e)x(e) = 0.70 +/- 0.1 cm(-1) in CS-II. On the basis of structural analysis, the discovery of the CS-II crystals implies stability of a large class of bromine hydrate structures and, therefore, polymorphism.     

Thermal expansion of structure-H clathrate hydrates

The temperature dependence of the unit cell parameters of two newly identified hexagonal structure clathrate hydrates of hexamethylethane (HME) and 2,2-dimethylbutane (DMB) have been measured by X-ray powder diffraction. The thermal expansion of the two distinct crystallographic axes was found to be inequivalent. However, the coefficients of cubic expansion are comparable to that in the cubic structure I and II hydrates. The larger thermal expansivity in the clathrate hydrates relative to ice is attributed to the weakening of the host lattice due to the internal pressure generated by the rattling motions of the encaged guests.Dedicated to Dr D. W. Davidson in honor of his great contributions to the sciences of inclusion phenomena.