Physical and numerical modeling of horizontal gas-hydrate wells

Methane hydrate, which is usually found under deep seabed or permafrost zones, is a potential energy resource for future years. Depressurization of horizontal wells bored in methane hydrate layer is considered as one possible method for hydrate dissociation and methane extraction from the hosting soil. Since hydrate is likely to behave as a bonding material to sandy soils, supported well construction is necessary to avoid well-collapse due to the loss of the apparent cohesion during depressurization. This paper describes both physical and numerical modeling of such horizontal support wells. The experimental part involves depressurization of small well models in a large pressure cell, while the numerical part simulates the corresponding problem. While the experiment models simulate only gas saturated initial conditions, the numerical analysis simulates both gas-saturated and more realistic water-saturated conditions based on effective stress coupled flow-deformation formulation of these three phases. © 2006 Taylor & Francis Group.

Physical and numerical modeling of horizontal gas-hydrate wells

Methane hydrate, which is usually found under deep seabed or permafrost zones, is a potential energy resource for future years. Depressurization of horizontalwells bored in methane hydrate layer is considered as one possible method for hydrate dissociation and methane extraction from the hosting soil. Since hydrate is likely to behave as a bonding material to sandy soils, supported well construction is necessary to avoid wellcollapse due to the loss of the apparent cohesion during depressurization. This paper describes both physical and numerical modeling of such horizontal support wells. The experimental part involves depressurization of small well models in a large pressure cell, while the numerical part simulates the corresponding problem. While the experiment models simulate only gas saturated initial conditions, the numerical analysis simulates both gas-saturated and more realistic water-saturated conditions based on effective stress coupled flow-deformation formulation of these three phases. © 2006 Taylor & Francis Group, London.

Numerical study on hydrate bearing sediments during gas production

The fully coupled methane hydrate model developed in Cambridge was adopted in this numerical study on gas production trial at the Eastern Nankai Trough, Japan 2013. Based on the latest experimental data of hydrate soil core samples, the clay parameters at Eastern Nankai site were successfully calibrated. With updated clay parameters and site geometry, a 50 days gas production trail was numerically simulated in FLAC2D. The geomechanical behaviour of hydrate bearing sediments under 3 different depressurization strategies were explored and discussed. The results from both axisymmetrical and plane-strain models suggest, the slope of the seabed only affects mechanical properties while no significant impact on the dissociation, temperature and pore pressure. For mechanical deformation after PT recovery, there are large settlements above the perforation zone and small uplift underneath the production zone. To validate the fully coupled model, numerical simulation with finer mesh in the hydrate production zone was carried out. The simulation results suggest good agreement between our model and JOE's results on history matching of gas and water production during trial. Parameter sensitivity of gas production is also investigated and concluded the sea water salinity is a dominant factor for gas production.

Study on the Dissociation Behavior of Gas Hydrate in Porous Sediment Based on Fractal Theory

The dissociation process of gas hydrate was regarded as a gas-solid reaction without solid production layer when the temperature was above the zero centigrade. Based on the shrinking core model and the fractal theory, a fractional dimension dynamical model for gas hydrate dissociation in porous sediment was established. The new approach of evaluating the fractal dimension of the porous media was also presented. The fractional dimension dynamical model for gas hydrate dissociation was examined with the previous experimental data of methane hydrate and carbon dioxide hydrate dissociations, respectively. The calculated results indicate that the fractal dimensions of porous media acquired with this method agree well with the previous study. With the absolute average deviation (AAD) below 10%, the present model provided satisfactory predictions for the dissociation process of methane hydrate and carbon dioxide hydrate.

Control mechanisms for gas hydrate production by depressurization in different scale hydrate reservoirs

The methane hydrate was formed in a pressure vessel 38 mm in id and 500 mm in length. Experimental works on gas production from the hydrate-bearing core by depressurization to 0.1, 0.93, and 1.93 MPa have been carried out. The hydrate reservoir simulator TOUGH-Fx/Hydrate was used to simulate the experimental gas production behavior, and the intrinsic hydration dissociation constant (K-0) fitted for the experimental data was on the order of 104 mol m(-2) Pa-1 s(-1), which was one order lower than that of the bulk hydrate dissociation. The sensitivity analyses based on the simulator have been carried out, and the results suggested that the hydrate dissociation kinetics had a great effect on the gas production behavior for the laboratory-scale hydrate-bearing core. However for a field-scale hydrate reservoir, the flow ability dominated the gas production behavior and the effect of hydrate dissociation kinetics on the gas production behavior could be neglected.

Study on the kinetics of hydrate formation in a bubble column

Gas hydrate formation experiments were performed using methane in the presence of tetrahydrofuran (THF) in aqueous solution in a transparent bubble column in which a single pipe or a sintered plate was used to produce bubbles. The mole fraction of THF in aqueous solution was fixed at 6%. The hydrate formation kinetic behaviors on the surface of the rising bubble, the mechanical stability of hydrate shell formed on the surface of the bubble, the interactions among the bubbles with hydrate shell were observed and investigated morphologically. The rise velocities of individual bubbles with hydrate shells of different thickness and the consumption rates of methane gas were measured. A kinetic model was developed to correlate the experimentally measured gas consumption rate data. It was found that the hydrate formation rate on the surface of the moving bubble was high, but the formed hydrate shell was not very easy to be broken up. The bubbles with hydrate shells tended to agglomerate rather than merge into bigger bubble. This kind of characteristic of hydrate shell hindered the further formation of hydrate and led to the lower consumption rate of methane. The consumption rate of methane was found to increase with the decrease of temperature or increase of pressure. The increase of gas flux led to a linear increase in consumption rate of methane. It was demonstrated that the developed kinetic model could be used to correlate the consumption rate satisfyingly.

Analysis for temperature and pressure fields in process of hydrate dissociation by depressurization

An improved axisymmetric mathematic modeling is proposed for the process of hydrate dissociation by depressurization around vertical well. To reckon in the effect of latent heat of gas hydrate at the decomposition front, the energy balance equation is employed. The semi-analytic solutions for temperature and pressure fields are obtained by using Boltzmann-transformation. The location of decomposition front is determined by solving initial value problem for system of ordinary differential equations. The distributions of pressure and temperature along horizontal radiate in the reservoir are calculated. The numeric results indicate that the moving speed of decomposition front is sensitively dependent on the well pressure and the sediment permeability. Copyright (C) 2010 John Wiley & Sons, Ltd.

Diversity, abundance and distribution of amoA-encoding archaea in deep-sea methane seep sediments of the Okhotsk Sea

The authors would like to thank Jin Sun, Jian Sun, Liangliang Kong, Nianshuang Wang, Chunhui Wang, Linbao Zhang and Ying Zhang for their assistance in the project. This work was supported by China Ocean Mineral Resources R&D Association grants DYXM-115-02-2-20 and DYXM-115-02-2-6, Hi-Tech Research and Development Program of China grant 2007AA091903, China National Natural Science Foundation grant 40576069, National Basic Research Program of China grant 2009CB219506 and the Fundamental Research Funds for the Central Universities of China grant 09CX05005A. M. G. K. was funded by incentive funds provided by the UofL-EVPR office and the US National Science Foundation (EF-0412129).

Diverse and Novel nifH and nifH-Like Gene Sequences in the Deep-Sea Methane Seep Sediments of the Okhotsk Sea

Diverse nifH and nifH-like gene sequences were obtained from the deep-sea surface sediments of the methane hydrate-bearing Okhotsk Sea. Some sequences formed novel families of the NifH or NifH-like proteins, of currently unresolved bacterial or archaeal origin. Comparison with other marine environments indicates environmental specificity of some of the sequences, either unique to the methane seep sediments of the Okhotsk Sea or to the general deep-sea methane seep sedimentary environments.

Methane concentration and delta deuterium of methane from the NGRIP ice core

The causes of past changes in the global methane cycle and especially the role of marine methane hydrate (clathrate) destabilization events are a matter of debate. Here we present evidence from the North Greenland Ice Core Project ice core based on the hydrogen isotopic composition of methane [dD(CH4)] that clathrates did not cause atmospheric methane concentration to rise at the onset of Dansgaard-Oeschger (DO) events 7 and 8. Box modeling supports boreal wetland emissions as the most likely explanation for the interstadial increase. Moreover, our data show that dD(CH4) dropped 500 years before the onset of DO 8, with CH4 concentration rising only slightly. This can be explained by an early climate response of boreal wetlands, which carry the strongly depleted isotopic signature of high-latitude precipitation at that time.

Oxygen isotopes and hydrocarbon composition of gas hydrate and hydrate water from ODP Leg 165 sites

Ocean Drilling Program (ODP) Leg 164 recovered a number of large solid gas hydrate from Sites 994, 996, and 997 on the Blake Ridge. Sites 994 and 997 samples, either nodular or thick massive pieces, were subjected to laboratory analysis and measurements to determine the structure, molecular and isotopic composition, thermal conductivity, and equilibrium dissociation conditions. X-ray computed tomography (CT) imagery, X-ray diffraction, nuclear magnetic resonance (NMR), and Raman spectroscopy have revealed that the gas hydrates recovered from the Blake Ridge are nearly 100% methane gas hydrate of Structure I, cubic with a lattice constant of a = 11.95 ± 0.05 angström, and a molar ratio of water to gas (hydration number) of 6.2. The d18O of water is 2.67 per mil to 3.51 per mil SMOW, which is 3.5-4.0 heavier than the ambient interstitial waters. The d13C and dD of methane are -66 per mil to -70 per mil and -201 per mil to -206 per mil, respectively, suggesting that the methane was generated through bacterial CO2 reduction. Thermal conductivity values of the Blake Ridge hydrates range from 0.3 to 0.5 W/(m K). Equilibrium dissociation experiments indicate that the three-phase equilibrium for the specimen is 3.27 MPa at 274.7 K. This is almost identical to that of synthetic pure methane hydrate in freshwater.

(Table 1) Chemical composition of interstitial waters from the area of gas hydrate generation, Sea of Okhotsk

Distributions of halogens (Cl, Br and I) in interstitial waters from sediments containing methane hydrate and in water of the hydrate itself are presented. High concentrations of halogens do not occur in interstitial waters from sediments that contain gas hydrates. The main reason for their low concentrations is the poverty of organic matter in sediments.

(Table 3) Noble gas composition of dissociated gas hydrate specimens from ODP Leg 164 sites

Fractionation of the noble gases should occur during formation of a Structure I gas hydrate from water and CH4 such that CH4 hydrate is greatly enriched in Xenon. Noble gas concentrations and fractionation factors (F[4He], F[22Ne], F[86Kr], and F[132Xe] as well as R/Ra) were determined for eight gas hydrate specimens collected on Leg 164 to evaluate this theoretical possibility and to assess whether sufficient quantities of Xe are hosted in oceanic CH4 hydrate to account for Xe "missing" from the atmosphere. The simplest explanation for our results is that samples contain mixtures of air and two end-member gases. One of the end-member gases is depleted in Ne, but significantly enriched in Kr and Xe, as anticipated if the source of this gas involves fractionation during Structure I gas hydrate formation. However, although oceanic CH4 hydrate may be greatly enriched in Xe, simple mass balance calculations indicate that oceanic CH4 hydrate probably represents only a minor reservoir of terrestrial Xe. Noble gas analyses may play an important role in understanding the dynamics of gas hydrate reservoirs, but significantly more work is needed than presented here.