甲烷水合物分子间势能的量子化学研究

用Hartree-Fock SCF和密度泛函(BLYP,B3LYP,MPW1PW91)方法对以结构-Ⅰ为单元的甲烷水合物进行了分子间势能的理论研究.该结构单元为正十二面体,其中包括20个水分子,甲烷分子在其中心.采用从头算HF/6-31G(d,p)对甲烷分子进行几何优化,采用ST2模型对水分子作几何优化.水-水间氢键势能Ehb(l)和水-甲烷间范德华势能Evdw(l)作为边长l的函数进行计算,计算时固定水和甲烷分子的几何形状.所有计算中均使用6-31G(d,p)基组.基组重叠误差(BSSE)经校正其上限和下限为水-水氢键能加以确定.由B3LYP经基组重叠误差(BSSE)校正得到的O—O距离为RO—O=0.280 nm,C—O距离RC—O=0.392 nm,比其他方法更接近实验值的0.282和0.395 nm.结果表明,在天然气水合物结构-Ⅰ中水-水分子对的氢键能(30~36 kJ/mol)大于水的二聚体(H2O)2氢键能(-22.6±2.9)kJ/mol,亦大于六角形冰的(-21.7±0.5)kJ/mol,十二面体结构为一稳定单元.以上分子间相互作用势能的结果为得出Lennard-Jones和Kihara势能参数提供了坚实的基础,此参数对分子动力学模拟天然气水合物是非常有用的.     

A computational study of the reaction of methane with methyl radical

The activation energy and optimized transition-state geometry for the abstraction of a hydrogen atom from methane by methyl radical have been calculated by the semiempirical methods MINDO /3 and MNDO . These results are compared with other semiempirical and ab initio results. The MINDO /3 method, based upon accuracy of the computed energy of activation, appears to be the computational method of greatest reliability. A method of locating the transition state on semiempirical surfaces is demonstrated.     

High-pressure transitions in methane hydrate

Three recent studies of the high-pressure transformations of methane hydrate [Chem. Phys. Lett. 325 (2000) 490; Proc. Natl. Acad. Sci. 97 (2000) 13484; Nature 410 (2001) 661] have reported apparently different behaviours. Detailed comparison of our X-ray and neutron diffraction data with those obtained in earlier work shows that there is in fact consistent behaviour on isothermal compression at room temperature with a transition from the clathrate I structure hydrate to a hexagonal hydrate with an unknown structure at 0.9 GPa.     

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.     

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). ...     

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.     

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.     

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.     

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.     

Secondary nucleation in the formation of methane crystal hydrate

The kinetic data on crystallization and a morphological analysis of a layer of CH₄ · 6H₂O hydrate crystals formed on the surface of water as a result of methane absorption showed that secondary nucleation occurred during hydrate crystallization. The mutual arrangement of crystals in the layer revealed photographically in situ was evidence that part of nuclei produced on the surface of previously formed crystals went away from the surface into solution and grew there independently of “mother” crystals, although the probability of such transfer into an immobile solution remained low. In view of this, a model of crystal growth generating secondary crystals was developed.     

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.     

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.     

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.     

Effect of cooling rate on methane hydrate formation in media

In order to study gas hydrate in media, formation of methane hydrate in three different media including loess, fine and coarse sands were studied. Five cooling rates were applied to form the methane hydrate. The nucleation time of the formation of methane hydrate with each cooling rate were measured for comparison. The experimental results show that the cooling rate is a significant factor affecting nucleation of methane hydrate and gas conversion. Under the same initial conditions, the faster the cooling rate, the shorter the nucleation time and the lower the methane gas conversion rate. The media also affect the formation process of methane hydrate within it. In loess, the gas conversion rate is lowest; in coarse sand, the gas conversion rate is the greatest; and in fine sand, it is in between. According to the study, it is found that the smaller the particle size of the media, the harder the methane hydrate forms within it.     

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.