Measurement and modeling of the effective thermal conductivity for porous methane hydrate samples

The effective thermal conductivities of gas-saturated porous methane hydrates were measured by a single-sided transient plane source (TPS) technique and simulated by a generalized fractal model of porous media that based on self-similarity.The density of porous hydrate,measured by the volume of the sample in the experimental system,was used to evaluate the porosity of methane hydrate samples.The fractal model was based on Sierpinski carpet,a thermal-electrical analogy technique and one-dimensional heat flow assumption.Both the experimental and computational results show the effective thermal conductivity of methane hydrate decreases with the porosity increase.The porosity of 0.3 can reduce the thermal conductivity of the methane hydrate by 25%.By analysis of the experimental data and the simulative result,the optimized thermal conductivity of the zero-porosity methane hydrate is about 0.7 W m-1K-1.     

Molecular dynamics simulations of the thermal conductivity of methane hydrate

Nonequilibrium molecular dynamics simulations with the nonpolarizable SPC/E (Berendsen et al., J. Phys. Chem. 1987, 91, 6269) and the polarizable COS/G2 (Yu and van Gunsteren, J. Chem. Phys. 2004, 121, 9549) force fields have been employed to calculate the thermal conductivity and other associated properties of methane hydrate over a temperature range from 30 to 260 K. The calculated results are compared to experimental data over this same range. The values of the thermal conductivity calculated with the COS/G2 model are closer to the experimental values than are those calculated with the nonpolarizable SPC/E model. The calculations match the temperature trend in the experimental data at temperatures below 50 K; however, they exhibit a slight decrease in thermal conductivity at higher temperatures in comparison to an opposite trend in the experimental data. The calculated thermal conductivity values are found to be relatively insensitive to the occupancy of the cages except at low (T相似文献   

Estimation of ultra-stability of methane hydrate at 1 atm by thermal conductivity measurement

Thermal conductivity of methane hydrate was measured in hydrate dissociation self-preservation zone by means of the transient plane source (TPS) technique developed by Gustafsson. The sample was formed from 99.9% (volume ratio) methane gas with 280 ppm sodium dodecyl sulfate (SDS) solution under 6.6 MPa and 273.15 K. The methane hydrate sample was taken out of the cell and moved into a low temperature chamber when the conversion ratio of water was more than 90%. In order to measure the thermal conductivity, the sample was compacted into two columnar parts by compact tool at 268.15 K. The measurements are carried out in the temperature ranging from 263.15 K to 271.15 K at atmospheric pressure. Additionally, the relationship between thermal conductivity and time is also investigated at 263.15 K and 268.15 K, respectively. In 24 h, thermal conductivity increases only 5.45% at 268.15 K, but thermal conductivity increases 196.29% at 263.15 K. Methane hydrates exhibit only minimal decomposition at 1 atm and the temperature ranging from 263.15 K to 271.15 K. At 1 atm and 268.15 K, the total gas that evolved after 24 h was amounted to less than 0.71% of the originally stored gas, and this ultra-stability was maintained if the test was lasted for more than two hundreds hours before terminating.     

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.     

Molecular insights into the heterogeneous crystal growth of si methane hydrate

In this paper we report a successful molecular simulation study exploring the heterogeneous crystal growth of sI methane hydrate along its [001] crystallographic face. The molecular modeling of the crystal growth of methane hydrate has proven in the past to be very challenging, and a reasonable framework to overcome the difficulties related to the simulation of such systems is presented. Both the microscopic mechanisms of heterogeneous crystal growth as well as interfacial properties of methane hydrate are probed. In the presence of the appropriate crystal template, a strong tendency for water molecules to organize into cages around methane at the growing interface is observed; the interface also demonstrates a strong affinity for methane molecules. The maximum growth rate measured for a hydrate crystal is about 4 times higher than the value previously determined for ice I in a similar framework (Gulam Razul, M. S.; Hendry, J. G.; Kusalik, P. G. J. Chem. Phys. 2005, 123, 204722).     

Experimental study on geochemical characteristic of methane hydrate formed in porous media

The natural occurrence of methane hydrates in marine sediments has been intensively studied over the past decades, and geochemical charac-teristic of hydrate is one of the most attractive research fields. In this paper, we discussed the geochemical anomaly during hydrate formation in porous media. By doing so, we also investigated the temperature influence on hydrate formation under isobaric condition. It turns out that sub-cooling is an important factor to dominate hydrate formation. Larger subcooling provides more powerful driving force for hydrate formation. During the geochemical anomaly research, six kinds of ions and the total dissolved salt (TDS) were measured before and after the experiment in different porous media. The result is that all kinds of ionic concentration increased after hydrate formation which can be defined as salting out effect mainly affected by gas consumption. But the variation ratio of different ions is not equal. Ca2+ seems to be the most significantly influenced one, and its variation ratio is up to 80%. Finally, we theoretically made a model to calculate the TDS variation, the result is in good accordance with measured one, especially when gas consumption is large.     

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.     

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.     

Experimental investigation of methane hydrate decomposition by depressurizing in porous media with 3-Dimension device

In order to simulate the behavior of gas hydrate formation and decomposition, a 3-Dimension experimental device was built, consisting of a high-pressure reactor with an inner diameter of 300 mm, effective height of 100 mm, and operation pressure of 16 MPa. Eight thermal resistances were mounted in the porous media at different depthes and radiuses to detect the temperature distribution during the hydrate formation/decomposition. To collect the pressure, temperature, and flux of gas production data, the Monitor and Control Generated System (MCGS) was used. Using this device, the formation and decomposition behavior of methane hydrate in the 20～40 mesh natural sand with salinity of 3.35 wt% was examined. It was found that the front of formation or decomposition of hydrate can be judged by the temperature distribution. The amount of hydrate formation can also be evaluated by the temperature change. During the hydrate decomposition process, the temperature curves indicated that the hydrate in the top and bottom of reactor dissociated earlier than in the inner. The hydrate decomposition front gradually moved from porous media surface to inner and kept a shape of column form, with different moving speed at different surface position. The proper decomposition pressure was also determined.     

Study of the permeability characteristics of porous media with methane hydrate by pore network model

The permeability in the methane hydrate reservoir is one of the key parameters in estimating the gas production performance and the flow behavior of gas and water during dissociation. In this paper, a three-dimensional cubic pore-network model based on invasion percolation is developed to study the effect of hydrate particle formation and growth habit on the permeability. The variation of permeability in porous media with different hydrate saturation is studied by solving the network problem. The simulation results are well consistent with the experimental data. The proposed model predicts that the permeability will reduce exponentially with the increase of hydrate saturation, which is crucial in developing a deeper understanding of the mechanism of hydrate formation and dissociation in porous media.     

A binary array method for calculating thermal conductivity of porous materials

The present research on thermal conductivity of porous materials is hampered by their highly random internal structure. Until now, the researches could not describe internal structure of porous materials, but macro calculation expressions were established. In order to solve the above problem, this paper presents a method of calculation that applies a binary array. This method can describe the internal structure of porous materials, reflects its randomness, and calculates the thermal conductivity of porous materials. The technique selects large enough porous material models and divides them in small enough grids. Each grid needs a binary conversion process to determine if the grid is matrix or pore. Finally, a grid array composed of matrix and pore is established. Based on this grid array, this study sets the initial and boundary conditions, and then calculates the thermal conductivity of porous materials by a numerical calculation according to internal heat-transfer mechanism. This research adopts the described method to calculate the thermal conductivity of porous aluminum. The result is consistent with that of the experimental measurement, thereby verifying the accuracy and reliability of this method.     

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.     

Double occupancy of large cavity of diethylamin+methane sH hydrate at low pressures

Structural Chemistry - New data on the composition of hexagonal structure III (sH) diethylamine+methane double hydrate are presented. The hydrate was synthesized from methane and aqueous solutions...     

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.     

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.     

Tetrahydrofuran hydrate decomposition characteristics in porous media

Many tetrahydrofuran (THF) hydrate properties are similar to those of gas hydrates. In the present work THF hydrate dissociation in four types of porous media is studied. THF solution was cooled to 275.15 K with formation of the hydrate under ambient pressure, and then it dissociated under ambient conditions. THF hydrate dissociation experiments in each porous medium were conducted three times. Magnetic resonance imaging (MRI) was used to obtain images. Decomposition time, THF hydrate saturation and MRI mean intensity (MI) were measured and analyzed. The experimental results showed that the hydrate decomposition time in BZ-4 and BZ-3 was similar and longer than that in BZ-02. In each dissociation process, the hydrate decomposition time of the second and third cycles was shorter than that of the first cycle in BZ-4, BZ-3, and BZ-02. The relationship between THF hydrate saturation and time is almost linear.     

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.