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.     

Quasiamorphous state of methane hydrate

It is established experimentally that at temperatures above 0°C formation of methane hydrate starts with the appearance of thin films of this substance on the surface of water; during isothermal storage, the films permeate throughout the whole volume of the sample. The maximal content of methane hydrate is several (volume) per cent. The viscosity of the sample is much lower than that of ice and much higher than that of water. A. A. Galkin Donetsk Physicotechnical Institute, Ukrainian Academy of Sciences. Translated fromZhurnal Strukturnoi Khimii, Vol. 39, No. 1, pp. 86–91, January–February, 1998.     

Study on the recovery of hydrogen from refinery (hydrogen + methane) gas mixtures using hydrate technology

A novel technique for separating hydrogen from (H2 CH4) gas mixtures through hydrate formation/dissociation was proposed. In this work, a systematic experimental study was performed on the separation of hydrogen from (H2 CH4) feed mixtures with various hydrogen contents (mole fraction x = 40%-90%). The experimental results showed that the hydrogen content could be enriched to as high as ~94% for various feed mixtures using the proposed hydrate technology under a temperature slightly above 0℃ and a pressure below 5.0 MPa. With the addition of a small amount of suitable additives, the rate of hydrate formation could be increased significantly. Anti-agglomeration was used to disperse hydrate particles into the condensate phase. Instead of preventing hydrate growth (as in the kinetic inhibitor tests), hydrates were allowed to form, but only as small dispersed particles. Anti-agglomeration could keep hydrate particles suspended in a range of condensate types at 1℃ and 5 MPa in the water-in-oil emulsion.     

Molecular-dynamics simulations of methane hydrate dissociation

Nonequilibrium molecular-dynamics simulations have been carried out at 276.65 K and 68 bar for the dissolution of spherical methane hydrate crystallites surrounded by a liquid phase. The liquid was composed of pure water or a water-methane mixture ranging in methane composition from 50% to 100% of the corresponding theoretical maximum for the hydrate and ranged in size from about 1600 to 2200 water molecules. Four different crystallites ranging in size from 115 to 230 water molecules were used in the two-phase systems; the nanocrystals were either empty or had a methane occupation from 80% to 100% of the theoretical maximum. The crystal-liquid systems were prepared in two distinct ways, involving constrained melting of a bulk hydrate system or implantation of the crystallite into a separate liquid phase. The breakup rates were very similar for the four different crystal sizes investigated. The method of system preparation was not found to affect the eventual dissociation rates, despite a lag time of approximately 70 ps associated with relaxation of the liquid interfacial layer in the constrained melting approach. The dissolution rates were not affected substantially by methane occupation of the hydrate phase in the 80%-100% range. In contrast, empty hydrate clusters were found to break up significantly more quickly. Our simulations indicate that the diffusion of methane molecules to the surrounding liquid layer from the crystal surface appears to be the rate-controlling step in hydrate breakup. Increasing the size of the liquid phase was found to reduce the initial delay in breakup. We have compared breakup rates computed using different long-range electrostatic methods. Use of the Ewald, minimum image, and spherical cut-off techniques led to more rapid dissociation relative to the Lekner method.     

Lattice dynamical simulation of methane hydrate

INS (Inelastic Neutron Scattering) spectrum of methane hydrate was measured on MARI (a direct-geometry chopper spectrometer) at Rutherford Appleton Laboratory, UK. Compared with ice Ih, it is found that the whole spectrum of methane hydrate moves toward high-energy by about 1.5 meV. Using lattice dynamical (LD) technique, computer simulations of methane hydrate were carried out. In the simulations, four potential models (BF, TIP3P, TIP4P, MCY) were employed to calculate the phonon density of states (PDOS). ...     

Observation of hydrogen in deuterated methane hydrate by maximum entropy method with neutron powder diffraction

The crystal structure of deuterated methane hydrate (structure I, space group: Pm(-)3n) was investigated by neutron powder diffraction at temperatures of 7.7-185 K. The scattering amplitude density distribution was examined by a combination of Rietveld method and maximum entropy method (MEM). The distribution of the D atoms in both D(2)O and CD(4) molecules was clarified from three-dimensional graphic images of the scattering amplitude density. The MEM results showed that there were low-density sites for the D atom of D(2)O in a particular location within the D(2)O cage at low temperatures. The MEM provided more reasonable results because of the decrease in the R factor that is attainable by this method. Accordingly, the low-density sites for the D atom of D(2)O probably exist within the D(2)O cage. This suggests that a spatial disorder of the D atom of D(2)O occurs at these sites and that hydrogen bonds between D(2)O molecules become partially weakened. With regard to the CD(4) molecules, there were high-density sites for the D atom of CD(4), and the density distribution of the C and D atoms was observed separately in the scattering amplitude density image. Consequently, the C-D bonds of CD(4) were not observed clearly because the CD(4) molecules had an orientational disorder. The D atoms of CD(4) were displaced from the line between the C and O atoms, and were located near the face center of the polygon in the cage. Accordingly, the D atoms of CD(4) were not bonded to specific O atoms. This result is consistent with the hydrophobicity of the CD(4) molecule. We also report the difference between the small and the large cages in the density distribution map and the temperature dependence of the scattering amplitude density.     

Direct measurement of methane hydrate composition along the hydrate equilibrium boundary

The composition of methane hydrate, namely n(w) for CH4.n(w)H2O, was directly measured along the hydrate equilibrium boundary under conditions of excess methane gas. Pressure and temperature conditions ranged from 1.9 to 9.7 MPa and 263 to 285 K. Within experimental error, there is no change in hydrate composition with increasing pressure along the equilibrium boundary, but n(w) may show a slight systematic decrease away from this boundary. A hydrate stoichiometry of n(w) = 5.81-6.10 H2O describes the entire range of measured values, with an average composition of CH4.5.99(+/-0.07)H2O along the equilibrium boundary. These results, consistent with previously measured values, are discussed with respect to the widely ranging values obtained by thermodynamic analysis. The relatively constant composition of methane hydrate over the geologically relevant pressure and temperature range investigated suggests that in situ methane hydrate compositions may be estimated with some confidence.     

Search for memory effects in methane hydrate: structure of water before hydrate formation and after hydrate decomposition

Neutron diffraction with HD isotope substitution has been used to study the formation and decomposition of the methane clathrate hydrate. Using this atomistic technique coupled with simultaneous gas consumption measurements, we have successfully tracked the formation of the sI methane hydrate from a water/gas mixture and then the subsequent decomposition of the hydrate from initiation to completion. These studies demonstrate that the application of neutron diffraction with simultaneous gas consumption measurements provides a powerful method for studying the clathrate hydrate crystal growth and decomposition. We have also used neutron diffraction to examine the water structure before the hydrate growth and after the hydrate decomposition. From the neutron-scattering curves and the empirical potential structure refinement analysis of the data, we find that there is no significant difference between the structure of water before the hydrate formation and the structure of water after the hydrate decomposition. Nor is there any significant change to the methane hydration shell. These results are discussed in the context of widely held views on the existence of memory effects after the hydrate decomposition.     

First determination of volume changes and enthalpies of the high-pressure decomposition reaction of the structure H methane hydrate to the cubic structure I methane hydrate and fluid methane

Pressure-temperature (P-T) conditions of the decomposition reaction of the structure H high-pressure methane hydrate to the cubic structure I methane hydrate and fluid methane were studied with a piston-cylinder apparatus at room temperature. For the first time, volume changes accompanying this reaction were determined. With the use of the Clausius-Clapeyron equation the enthalpies of the decomposition reaction of the structure H high-pressure methane hydrate to the cubic structure I methane hydrate and fluid methane have been calculated.     

Formation of methane hydrate from polydisperse ice powders

Neutron diffraction runs and gas-consumption experiments based on pressure-volume-temperature measurements are conducted to study the kinetics of methane hydrate formation from hydrogenated and deuterated ice powder samples in the temperature range of 245-270 K up to high degrees of transformation. An improved theory of the hydrate growth in a polydisperse ensemble of randomly packed ice spheres is developed to provide a quantitative interpretation of the data in terms of kinetic model parameters. This paper continues the research line of our earlier study which was limited to the monodisperse case and shorter reaction times (Staykova et al., 2003). As before, we distinguish the process of initial hydrate film spreading over the ice particle surface (stage I) and the subsequent hydrate shell growth (stage II) which includes two steps, i.e., an interfacial clathration reaction and the gas and water transport (diffusion) through the hydrate layer surrounding the shrinking ice cores. Although kinetics of hydrate formation at stage II is clearly dominated by the diffusion mechanism which becomes the limiting step at temperatures above 263 K, both steps are shown to be essential at lower temperatures. The permeation coefficient D is estimated as (1.46 +/- 0.44) x 10(-12) m2/h at 263 K with an activation energy Q(D) approximately 52.1 kJ/mol. This value is close to the energy of breaking hydrogen bonds in ice Ih and suggests that this process is the rate-limiting step in hydrate formation from ice in the slower diffusion-controlled part of the reaction.     

Thermodynamic properties of methane hydrate in quartz powder

Using the experimental method of precision adiabatic calorimetry, the thermodynamic (equilibrium) properties of methane hydrate in quartz sand with a grain size of 90-100 microm have been studied in the temperature range of 260-290 K and at pressures up to 10 MPa. The equilibrium curves for the water-methane hydrate-gas and ice-methane hydrate-gas transitions, hydration number, latent heat of hydrate decomposition along the equilibrium three-phase curves, and the specific heat capacity of the hydrate have been obtained. It has been experimentally shown that the equilibrium three-phase curves of the methane hydrate in porous media are shifted to the lower temperature and high pressure with respect to the equilibrium curves of the bulk hydrate. In these experiments, we have found that the specific heat capacity of the hydrate, within the accuracy of our measurements, coincides with the heat capacity of ice. The latent heat of the hydrate dissociation for the ice-hydrate-gas transition is equal to 143 +/- 10 J/g, whereas, for the transition from hydrate to water and gas, the latent heat is 415 +/- 15 J/g. The hydration number has been evaluated in the different hydrate conditions and has been found to be equal to n = 6.16 +/- 0.06. In addition, the influence of the water saturation of the porous media and its distribution over the porous space on the measured parameters has been experimentally studied.     

Preservation phenomena of methane hydrate in pore spaces

Dissociation processes of methane hydrate synthesized with glass beads were investigated using powder X-ray diffraction and calorimetry. Methane hydrate formed with coarse glass beads dissociated quickly at 150-200 K; in this temperature range methane hydrate dissociates at atmospheric pressure. In contrast, methane hydrate formed with glass beads less than a few microns in size showed very high stability up to just below the melting point of ice, even though this temperature is well outside the zone of thermodynamic stability of the hydrate. The rate-determining steps for methane hydrate dissociation within pores are also discussed. The experimental results suggest that methane hydrate existing naturally within the pores of fine particles such as mud at low temperatures would be significantly more stable than expected thermodynamically.     

Molecular simulation of non-equilibrium methane hydrate decomposition process

We recently performed constant energy molecular dynamics simulations of the endothermic decomposition of methane hydrate in contact with water to study phenomenologically the role of mass and heat transfer in the decomposition rate [S. Alavi, J.A. Ripmeester, J. Chem. Phys. 132 (2010) 144703]. We observed that with the progress of the decomposition front temperature gradients are established between the remaining solid hydrate and the solution phases. In this work, we provide further quantitative macroscopic and molecular level analysis of the methane hydrate decomposition process with an emphasis on elucidating microscopic details and how they affect the predicted rate of methane hydrate decomposition in natural methane hydrate reservoirs. A quantitative criterion is used to characterize the decomposition of the hydrate phase at different times. Hydrate dissociation occurs in a stepwise fashion with rows of sI cages parallel to the interface decomposing simultaneously. The correlations between decomposition times of subsequent layers of the hydrate phase are discussed.     

Phase transformation of methane hydrate under high pressure

In situ Raman spectroscopy is employed to study the phase behavior of methane hydrate at high pressure. The structure 1 of methane hydrate can be maintained up to 950 MPa at 299 K. The transformation of structure I<-->structure H+water+CH4 occurs at 880 MPa and 323 K. The structure H of methane hydrate, however, decomposes to methane and water at 960 MPa and 348 K. The initiation mechanism of methane hydrate sI is also discussed.     

Electronic state of small and large cavities for methane hydrate

It is well known that methane hydrate is aggregates of small and large hydrogen bonded water cavities (composed of 12 pentagonal faces of 20 water molecules, and 12 pentagonal and two hexagonal faces of 24 water molecules, respectively) where one methane molecule is encaged. We calculated the methane molecule in vacuum, the small and large cavities by ab initio MO method to clarify the electronic state. The proton of methane in the cavities is shown to form the weak hydrogen bond (O...H[bond]C) between methane and four water molecules, and the H-bond lengths and energies in the small and large cavities were estimated as (0.293 nm, 6.8 kJ/mol) and (0.309 nm, 5.2 kJ/mol), respectively. The calculated values of symmetric C[bond]H stretching frequencies and (13)C-NMR chemical shieldings of the methane in the two cluster cavities show good agreement with the experimental ones observed by Sum et al. and Ripmeester and coworker, respectively.     

Molecular dynamics study of methane hydrate formation at a water/methane interface

We present molecular dynamics simulation results of a liquid water/methane interface, with and without an oligomer of poly(methylaminoethylmethacrylate), PMAEMA. PMAEMA is an active component of a commercial low dosage hydrate inhibitor (LDHI). Simulations were performed in the constant NPT ensemble at temperatures of 220, 235, 240, 245, and 250 K and a pressure of 300 bar. The simulations show the onset of methane hydrate growth within 30 ns for temperatures below 245 K in the methane/water systems; at 240 K there is an induction period of ca. 20 ns, but at lower temperatures growth commences immediately. The simulations were analyzed to calculate hydrate content, the propensity for hydrogen bond formation, and how these were affected by both temperature and the presence of the LDHI. As expected, both the hydrogen bond number and hydrate content decreased with increasing temperature, though little difference was observed between the lowest two temperatures considered. In the presence of PMAEMA, the temperature below which sustained hydrate growth occurred was observed to decrease. Some of the implications for the role of PMAEMA in LDHIs are discussed.     

Thermal conductivity of methane hydrate from experiment and molecular simulation

A single-sided transient plane source technique has been used to determine the thermal conductivity and thermal diffusivity of a compacted methane hydrate sample over the temperature range of 261.5-277.4 K and at gas-phase pressures ranging from 3.8 to 14.2 MPa. The average thermal conductivity, 0.68 +/- 0.01 W/(m K), and thermal diffusivity, 2.04 x 10(-7) +/- 0.04 x 10(-7) m2/s, values are, respectively, higher and lower than previously reported values. Equilibrium molecular dynamics (MD) simulations of methane hydrate have also been performed in the NPT ensemble to estimate the thermal conductivity for methane compositions ranging from 80 to 100% of the maximum theoretical occupation, at 276 K and at pressures ranging from 0.1 to 100 MPa. Calculations were performed with three rigid potential models for water, namely, SPC/E, TIP4P-Ew, and TIP4P-FQ, the last of which includes the effects of polarizability. The thermal conductivities predicted from MD simulations were in reasonable agreement with experimental results, ranging from about 0.52 to 0.77 W/(m K) for the different potential models with the polarizable water model giving the best agreement with experiments. The MD simulation method was validated by comparing calculated and experimental thermal conductivity values for ice and liquid water. The simulations were in reasonable agreement with experimental data. The simulations predict a slight increase in the thermal conductivity with decreasing methane occupation of the hydrate cages. The thermal conductivity was found to be essentially independent of pressure in both simulations and experiments. Our experimental and simulation thermal conductivity results provide data to help predict gas hydrate stability in sediments for the purposes of production or estimating methane release into the environment due to gradual warming.     

Mechanical property of artificial methane hydrate under triaxial compression

Methane production from hydrate reservoir may induce seabed slide and deformation of the hydrate-bearing strata. The research on mechanical properties of methane hydrate is considered to be important for developing an efficient methane exploitation technology. In this paper, a triaxial test system containing a pressure crystal device was developed with the conditions to stabilize the hydrate. A series of triaxial shear tests were carried out on artificial methane hydrate specimen. In addition, mechanical characteristics of methane hydrate were studied with the strain rates of 0.1 and 1.0 mm/min, respectively, under the conditions of different temperatures (T = -5, -10, and -20 ℃) and confining pressures (P = 0, 5, 10, 15, and 20 MPa). The preliminary results show that when the confining pressure was less than 10 MPa, the increase of confining pressure leaded to the enhancement of shear strength. Furthermore, the decreasing temperature and the increasing strain rate both caused the increase in shear strength.