- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Journal of Geological Research

Volume 2011 (2011), Article ID 462156, 15 pages

http://dx.doi.org/10.1155/2011/462156

## Simulation of Methane Recovery from Gas Hydrates Combined with Storing Carbon Dioxide as Hydrates

Fraunhofer Institute for Environmental, Safety, and Energy Technology UMSICHT, Osterfelder Strasse 3, 46047 Oberhausen, Germany

Received 15 June 2011; Accepted 1 August 2011

Academic Editor: Michela Giustiniani

Copyright © 2011 Georg Janicki et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

In the medium term, gas hydrate reservoirs in the subsea sediment are intended as deposits for carbon dioxide (CO_{2}) from fossil fuel consumption. This idea is supported by the thermodynamics of CO_{2} and methane (CH_{4}) hydrates and the fact that CO_{2} hydrates are more stable than CH_{4} hydrates in a certain *P*-*T* range. The potential of producing methane by depressurization and/or by injecting CO_{2} is numerically studied in the frame of the SUGAR project. Simulations are performed with the commercial code STARS from CMG and the newly developed code HyReS (hydrate reservoir simulator) especially designed for hydrate processing in the subsea sediment. HyReS is a nonisothermal multiphase Darcy flow model combined with thermodynamics and rate kinetics suitable for gas hydrate calculations. Two scenarios are considered: the depressurization of an area 1,000 m in diameter and a one/two-well scenario with CO_{2} injection. Realistic rates for injection and production are estimated, and limitations of these processes are discussed.

#### 1. Introduction

Gas hydrates are ice-like solid compounds of water and gas molecules (clathrates) which are stable at low temperature and elevated pressure [1]. The water molecules build out cages by hydrogen bonds in which gas molecules are embedded. Generally, gas hydrates can contain different guest molecules in different cages, depending on their sizes and the availability of guest molecules under given thermodynamic conditions, but methane is the prevalent gas in natural gas hydrates. The exploitation of natural gas hydrate deposits that are known in various permafrost regions and submarine sediments all over the world is in the focus of several research groups because the amount of methane to be recovered could overcome future energy shortages. The greenhouse gas CO_{2} is able to build hydrates too, and these hydrates are thermodynamically more stable than methane hydrates. The possibility to destabilize methane hydrate by injecting CO_{2} as pressurized gas or in liquid form was verified in several small-scale experiments carried out by different research groups (see [2–6]). Thus, the combination of both processes offers the opportunity to open up new energy resources as well as to combat climate change by reducing CO_{2} emissions. However, the technical realization of this combination of processes has to face various challenges. Besides the technical and economic efforts for drilling in submarine sediments or in deep layers in permafrost regions, these challenges concern the reaction kinetics and transport resistances within the sediments in which methane hydrates are embedded in natural deposits.

Thus, to find the best strategy for methane recovery from a specific deposit with or without CO_{2} sequestration, a large variety of parameters describing the properties of the particular layer as well as the time- and location-dependent thermodynamic conditions have to be considered. Within the framework of the German SUGAR (SUbmarine GAs Hydrate Reservoirs) project, strategies to produce natural gas from marine methane hydrates and simultaneously store CO_{2} as hydrates are explored. Before undertaking drilling tests, numerical simulations of the local processes are necessary and helpful. For this purpose, a new scientific simulation model called UMSICHT HyReS was developed to describe the methane production from submarine hydrate layers and the exchange of methane by carbon dioxide. In addition, the commercially available simulation tool STARS (CMG Ltd., Canada) was used. In the following, the new simulation tool is described, and the results of the calculations based on particular reservoir parameters, reaction kinetics, and extraction techniques are outlined.

#### 2. Principles of Methane Extraction with CO_{2} Sequestration

##### 2.1. Exploitation Strategies and Their Limitations

All strategies to destabilize methane hydrate are aimed to change the thermodynamic and/or chemical properties within the reservoir in order to overcome the preconditions for stability. The thermodynamic equilibrium curves for methane and carbon dioxide hydrates given in Figure 1 reveal that either thermal stimulation (A-B), depressurization (A-C), or a combination of both approaches may be the best choice to destabilize a hydrate and to recover the pure gas.

The variety of measures to supply heat into a hydrate layer for thermal destabilization ranges from injecting warm water or steam into a well over injecting fuel and oxygen to especially designed on-site combustion chambers up to electromagnetic heating of a particular region around a well. Recently, Schicks et al. [9] suggested the in situ catalytic oxidation of methane to supply the heat that is needed to decompose the hydrate.

Depressurization of a hydrate-bearing region is achieved by suction through one or more wells. According to the equilibrium curve (Figure 1), lowering the pressure will shift the equilibrium to lower temperatures. However, the release of methane gas from its hydrate consumes heat (e.g., Goel [10]) so that the reaction only continues until the temperature reaches the new equilibrium point. Su et al. [11] reported on depressurization experiments that were carried out in a high-pressure reactor in which methane hydrate had been built within a layer of natural sand. It was shown that the gas production rate does not only depend on the driving force (expressed as the difference of actual and equilibrium pressure). Reaction kinetics, heat and mass transfer limitations, and multiphase flow effects are also important. In addition, the possible formation of ice that may hamper mass and heat transfer must be considered.

##### 2.2. Classification of Hydrate Deposits

In recent years, the classification of natural gas hydrate deposits in three classes according to their geologic and reservoir conditions has become accepted. According to Moridis and Collet [12], class 1 deposits comprise a hydrate layer and an underlying two-phase fluid zone with free gas and formation water. They describe class 2 deposits as a hydrate layer with an underlying mobile water zone with no free gas, and single-hydrate layers without underlying mobile fluids are characterized as class 3 accumulations. Class 1 deposits are the most attractive deposits with respect to exploitation strategies because the bottom of such a layer that is in contact with mobile gas must be in or near thermodynamic equilibrium. Therefore, the energy needed to destabilize the hydrate is minimal compared to classes 2 and 3 in which the hydrate interval may be far within the stability zone. In addition, an underlying layer of mobile water and/or gas provides a better conductivity within the sediment and supplies additional heat. This expands the reaction area where hydrates decompose. Besides these three classes, a fourth class of deposits, that is to say disperse, low-saturated hydrate accumulations in oceanic sediments, was also investigated by Moridis and Sloan [13]. They evaluated an exploitation of such deposits as noneconomic.

Structural properties of geological formations also comprise the conditions for the injection of CO_{2} and the building up and stability of CO_{2} hydrates. Uddin et al. [14] used the thermal reservoir simulator STARS to initially simulate the CO_{2} hydrate formation in four simplified geological formations that were classified with respect to their mean porosity and initial permeability. The parameters (mean porosity 0.25–0.35, permeability 10–1,000 mD) were chosen according to geological formations identified in the Mallik deposit in Canada's Mackenzie Delta which is a permafrost region. The SUGAR project is explicitly aimed at submarine gas hydrates. Hester and Brewer [15] showed over 70 different locations of subaquatic deposits all over the world that had been identified by 2009, with less than half of them being actually known deposits, the others described as inferred deposits that had been identified with indirect markers so far. Some data on the structural properties like mean porosities and intrinsic permeabilities of hydrate-bearing sediments are available for some sites that have been investigated so far (see [14, 16, 17]). Thus mean porosities reported so far range from 0.2 to 0.5 (e.g., [16, 17]), permeabilities range from 0.1 [17] to 1,000 mD [18].

##### 2.3. Simplified Reservoir Model

For simulations, a simplified reservoir model built of layers with homogeneous properties for porosity, permeability, and initial hydrate saturation was used (see Figure 2). According to the classification of Moridis and Collet [12], the reservoir corresponds to class 1 type.

In the simulations presented, the focus is on scenarios where only the hydrate layer with no-flow boundaries for mass and/or heat is considered. The initial conditions are shown in Table 2. As initial pressure, the equilibrium pressure at a given temperature was chosen in order to avoid artificial hydrate formation or decomposition during the simulation run.

#### 3. Simulation Tools

Within the frame of the SUGAR research project, two simulation programs are used:(i)the simulator STARS (Steam, Thermal and Advanced Processes Reservoir Simulator, Computer Modeling Group Ltd. (CMG), Canada), an oil and gas industry standard reservoir simulator,(ii)the simulation code UMSICHT HyReS, a newly developed in-house code.

STARS was primarily developed to model oil and gas flow through 3-dimensional sediments, but it can be set up to deal with hydrate formation and decomposition as well. Due to the complexity of subsea hydrate systems and several restrictions in using STARS regarding full thermodynamics of gas hydrate physics, a second in-house code called UMSICHT HyReS was developed within the project to achieve full flexibility.

##### 3.1. Model Concept

The underlying operation scenario for the development of a new model is methane production from subsea gas hydrate accumulations by depressurization, thermal stimulation, or injection of liquid/soluted CO_{2}. The latter can result in an exchange of methane by carbon dioxide in the hydrate lattice or in a two-step decomposition/formation process with gaseous methane and hydrate-stabilized CO_{2} in the end [5, 6, 19, 20]. Conservation laws for such a complex situation must be formulated for a gas and a liquid phase flowing through sediment partially blocked by hydrates. Assuming a stable nonmoving sediment phase, only sediment porosity will be constant, whereas the saturation of gas, liquid, and hydrates is free variables. Therefore, the flow within the sediment strongly depends on the time-dependent distribution of phase saturations within the calculated volume, and this is a fundamental feature of the UMSICHT HyReS model. The second fundamental concept is the perception of decomposition/formation kinetics within the hydrates. According to works of Haeckel et al. [21, 22], gas hydrate formation is controlled by a liquid phase formation equilibrium. Conceptually this means that hydrate formation is a two-step process of solution of the gas component in the liquid and formation of hydrate within the liquid. From a thermodynamic point of view, this is more complicated than the well-known kinetics of the *Kim-Bishnoi* type. The third remarkable feature of the new model is the distinct formulation of energy balances for gas, liquid, hydrate, and sediment phase which was made for dealing with special operation situations in thermally stimulated gas production.

##### 3.2. Flow Models

The velocity of a fluid phase *α* in the multiphase system is given by *Darcy’s* law:
Equations of the Carman-Kozeny type are widely used as permeability models for subsea sediments. In UMSICHT HyReS, formulations of *Brooks-Corey*, *van Genuchten*, and the model used in STARS are implemented. In this work, the *Brooks-Corey* model (which is also used by Sun and Mohanty [23]) is applied:
With (3), the blocking of the flow path by solid gas hydrates is included. A similar formulation for the absolute permeability is used in STARS:
Here, is the permeability under starting conditions including initial hydrate saturation, whereas in (3) is the permeability under pure conditions without hydrates. Parameterization of these equations is given later in Table 2. The parameters for (6) used in STARS are set to give equal flow as (3) in UMSICHT HyReS. With the concept of relative permeability, the flow-model includes the time-dependent progress of gas, liquid, and hydrate phase saturations and is able to describe subsea sediment flow under changing conditions during methane gas production.

##### 3.3. Kinetic Model

Decomposition and formation of gas hydrates are modelled by a linear driving force approach. For methane and carbon dioxide hydrate formation, the following gross kinetic is assumed:
As mentioned before, the kinetic model assumes that the formation process is controlled by the liquid phase concentrations of the hydrate-forming gas. Within this approach, the concentrations are liquid phase concentrations of the gas components in thermodynamic equilibrium with the hydrates present in the volume. These equilibrium concentrations are different from the solution equilibria between gas and liquid phase (solubility). The hydrate formation from a component present in the gas phase can be modelled as a two-step mechanism of gas solution (step one) and hydrate formation by (7) (step two). Decomposition is modelled by the same equations, but in opposite direction. A formulation for the thermodynamic equilibrium of solution and hydrate formation of methane and carbon dioxide is given by Tishchenko et al. [7]; kinetic constants are adapted from data given by Haeckel et al. [21, 22], and interphase areas for methane and carbon dioxide hydrate formation/decomposition are calculated according to Sun and Mohanty [23] by
Equation (8) ensures that the hydrate decomposition results in a hydrate-free environment; in case of hydrate formation in a hydrate-free environment, a small initial saturation is a prerequisite. Calculations with the STARS code are done with the well-known *Kim-Bishnoi* kinetic approach:
Here, the distance to thermodynamic equilibrium given in gas phase partial pressures is used as the driving force for hydrate formation and decomposition. In this approach, liquid phase concentrations of the hydrate forming gas have no influence on hydrate formation.

##### 3.4. Other Submodels and Property Data

Within the scope of UMSICHT HyReS, a large number of submodels for capillary pressure, property data, thermal and diffusional conductivity, and for heat and mass transfer between coexistent phases are used. Most submodels depend on one or more free solution variables (like temperature, pressure, concentration, or saturation), resulting in a strong numerical coupling of balance equations. In the following, the most relevant models are presented.

Capillary pressure is given by a modified *Brooks-Corey* approach also used in [23]:
Additionally, the *van Genuchten* equation is available:
The effective saturation is given by (5). Mass transfer of gaseous components to liquid and vice versa is modelled by a two-film theory linear approach
The mass transfer coefficient is set equal to m/s for all components, gas-liquid interphase area is calculated from
with a constant bubble diameter of m. Property data of gas phase components are taken from *NIST Standard Reference Database* [8] and are molar-weighted; for the gas phase density, the *Peng-Robinson* equation is used. Liquid phase property data are also taken from NIST and other available sources; special corrections for saline seawater are used for density (UNESCO Standard Reference Equation), thermal conductivity [24], viscosity [25], and diffusion coefficients [26]:
Overall liquid phase heat conductivity and viscosity are calculated with volumetric and molar-weighted functions:
Diffusion coefficients are calculated from *IFM-GEOMAR* data functions [26]:
and the *Boudreau* tortuosity correction [27]:
Other property data for sediment and hydrates are given in Table 1.

##### 3.5. Governing Equations System

Conservation laws are formulated for a two-phase fluid flow through a stable sediment matrix as a heterogeneous continuum model with advection-diffusion type equations for component mass conservation. In UMSICHT HyReS, there are five coexistent phases:(i)gas phase (consists of methane, water, and carbon dioxide),(ii)liquid phase (consists of water, salt, methane, and carbon dioxide),(iii)methane hydrate phase,(iv)carbon dioxide hydrate phase, (v)sediment.

The unsteady conservation conditions to describe the submarine layer flow of gas and liquid lead to a system of partial differential equations. The solution variables of the model are(i)pressure of gas and liquid phase,(ii)saturation of gas, liquid, methane hydrate, and carbon dioxide hydrate,(iii)temperature of gas, liquid, hydrate, and sediment phase,(iv)concentrations of methane, water, and carbon dioxide in the gas phase,(v)concentrations of water, salt, methane, and carbon dioxide in the liquid phase.

Based on these assumptions and settings, the conservation equations used in the UMSICHT HyReS code are given as follows.

*Mass conservation of fluid and hydrate phases*

*Mass conservation of components*

*Energy conservation for gas, liquid, hydrate, and sediment*

Fluid velocity of phase is given by *Darcy’s *law as stated in Section 3.2. The conservation laws given above are completed by two *equations-of-state: *
a *capillary pressure condition*,
respectively, and the *phase closing condition: *
respectively.

##### 3.6. Boundary Conditions

Boundary conditions (BCs) are necessary to solve the PDE system given above under specific physical constraints. Generally, two types of BC are used in UMSICHT HyReS,(i)“no-flow” conditions at boundaries impermeable for mass and/or heat flow,(ii)“flow” conditions at boundaries permeable for mass and/or heat flow.

Both BC can be mixed, for example, if a boundary is permeable for heat, but impermeable for mass. If are convective and are diffusive fluxes at a boundary, the general form of these BC is for “no-flow” boundaries, and for “flow” boundaries where flow conditions at the point behind the boundary are indexed by “out.” With the general equations (26) and (27), boundary conditions for the solution vector and its derivatives can be set up for a special physical scenario. In UMSICHT HyReS, all boundaries are outer boundaries of the calculated volume and there exist no inner BC.

##### 3.7. Numerical Solution

The system of partial differential equations given in Section 3.5 has the general form in which is the vector of solution variables and both the mass matrix and the generalized diffusion coefficient matrices may depend on . The variable can take values of = 0, 1, or 2 indicating the coordinate system in the respective spatial direction (: Cartesian, : cylindrical, : spherical). Equation (28) exhibits the typical structure of a system of convection-diffusion equations. In UMSICHT-HyReS, the well-established (vertical) method of lines is applied for solving: spatial derivatives are approximated by finite differences on a structured rectangular grid in the plane, whereas temporal derivatives are left in their continuous form. The result of this semidiscretization is a large system of ordinary differential equations for the values of the solution variables in each grid point. Highly efficient solvers are available for such systems with banded Jacobian matrix, for example, the differential-algebraic integrator RADAU [28].

Below, the discretization schemes used for the first and (generalized) second-order derivatives in direction are specified. The expressions in direction are treated analogously. As will be discussed later in Section 4, the transport of the conserved quantities is strongly convection-dominated. For this reason, special care must be taken in the choice of the discretization scheme for the first-order spatial derivatives in order to avoid unphysical oscillations in the solution variables. In UMSICHT HyReS, a second-order upwind scheme combined with a flux-limiter method as described in [29] is used.

On the interface between cell and cell (indexed ), flux-limited left () and right () states of the solution variables are defined by with The first-order spatial derivatives are then approximated by The (generalized) second-order spatial derivatives are discretized by standard central differences: For the calculation of the generalized diffusion coefficients at the cell interfaces, there is a choice whether the arithmetic mean or the flux-limited left/right state of the respective solution variable is used: or

#### 4. Operation Scenarios and Simulation Results

As mentioned earlier, depressurization and the supply of heat or a combination of both methods are the most obvious strategies to destabilize methane hydrate and to recover the emerging gas. An alternative means is to change the methane concentration in the hydrate-building phase by supplying a methane-free substance to the reservoir, for example, water or carbon dioxide. If carbon dioxide is supplied and the operating conditions are chosen properly, it may even be possible to substitute methane by carbon dioxide in hydrate form.

After conforming the relevant inputs in STARS to those of UMSICHT HyReS, two scenarios were set up for simulation:(1)first scenario: pure depressurization of a reservoir with closed boundaries,(2)second scenario: depressurization combined with CO_{2} injection.

In the first scenario, the pressure in the centre of a cylindrical reservoir with a diameter of 1,000 m is reduced from equilibrium conditions by approximately 50 bars, and methane gas is produced for a period of three years (Case 1). STARS and UMSICHT HyReS are used to study the effect of a variation in model parameters on methane production. The results obtained from both simulators are compared and found in close agreement (see Figure 6).

The second scenario is a continuation of the first and was simulated with STARS. After three years of methane recovery at a bottom hole pressure of 30 bars, CO_{2} is injected at the same well for another three years at a maximum rate of 8,000 STD m^{3}/day (Case 2.1).

A variation of the second scenario is the simultaneous production of methane at one well and CO_{2} storage at a second well in a Cartesian reservoir (Case 2.2). The distance between the two wells is 500 m, and the maximum rate of injection is 8,000 STD m^{3}/day. Both reservoirs (radial and Cartesian) almost have the same volume so that the same initial amount of methane is present.

##### 4.1. First Scenario: Depressurization of a Closed Reservoir (Case 1)

The reservoir volume for the depressurization scenario is a cylindrical region with a diameter of 1,000 m and a height of 20 m. The well is situated in the centre of the region and has a diameter of 1 m. Starting from CH_{4} hydrate equilibrium conditions at about 77 bars and 10°C, the liquid pressure at the well is lowered to 30 bars. Methane hydrate is dissolved and a flow is induced in the direction to the well. Main initial values and model parameters are given in Table 2.

Figures 3, 4, and 5 show the history of pressure, temperature, and methane hydrate saturation over the radius of the reservoir. Pressure decreases slowly from 77 bars to 30 bars within the reservoir because of the restricted permeability of the hydrate formation. Temperature starts at 10°C and decreases to 1°C (the equilibrium temperature corresponding to 30 bars); in the same period of time, the methane hydrate saturation decreases due to methane production in the volume. The process will finally stop at a lower temperature in equilibrium conditions. The overall simulation time was three years (corresponds to the lowest curve in all diagrams). The curves in all diagrams correspond to times , 1, 10, and 36 months.

The production curves for methane gas obtained with UMSICHT HyReS (bold lines) and STARS (thin lines) are shown in Figure 6. To get a better impression of the curve progression, the production period was set to 30 years instead of 3 years. The maximum production rate of about 40,000 STD m^{3}/day is quickly reached right at the beginning of the depressurization. For the next five years, methane production only slightly decreases; after five years, nearly 70 million STD m^{3} have been produced at the well (UMSICHT HyReS). Thereafter, the supply of gas from the outer part of the reservoir slows down because the pressure gradient towards the well becomes increasingly smaller. Thus, although after 30 years of production, the gas saturation in the reservoir is still at a level of about 12%, the lack of the pressure gradient driving force prevents the free methane gas from being recovered.

Bearing in mind that the models used in UMSICHT HyReS and STARS are of high complexity and differ slightly, the production rate and cumulative gas amount after thirty years are in close agreement. The deviation is mainly caused by different kinetic approaches (see (7)–(9)) and various thermodynamic parameters which are -dependent (as in UMSICHT HyReS) instead of constant values (as in STARS).

Taking the above scenario as the base case, model parameters are varied in order to study the strength of their influence on methane gas production. Table 3 shows the parameters that are varied together with the resulting cumulative methane production after three years. The column in the middle represents the parameter settings of the base case. For the parameter study, the values of the base case are lowered (upper row) and increased (lower row), and the respective cumulative methane production is reported for simulations carried out with UMSICHT HyReS. As can be seen, mass transfer and kinetics are fast enough not to hinder gas recovery. The permeability of the reservoir is found to be the bottleneck of the process as it controls the overall performance of methane production.

In addition to the variation of parameters, the boundaries of the reservoir have been extended to include the over and under burdens in the simulation. Both zones are set free of hydrate and gas. As the temperature increases with the depth of the reservoir, a vertical thermal gradient of 30°C/km is assumed. The initial pressure in each layer is taken as the equilibrium pressure corresponding to the given temperature. The values for pressure and temperature listed in Table 4 are set at the uppermost layer of each zone.

As mentioned before, due to the depressurization and the endothermic decomposition of hydrates, the reservoir cools down until a new equilibrium is reached. At this state, no more gas is released from the hydrate, and the production rate decreases. The inclusion of the over and under burdens in the simulation leads to an enhanced gas production: after three years, 132 million STD m^{3} have been recovered, whereas only 50 million STD m^{3} are extracted for the base case. The reason is the supply of additional heat from the upper and lower layers. Figure 7 shows that more hydrate decomposes in the vicinity of the layer boundaries than in the middle of the hydrate layer. In order to reproduce that effect in the simulations, the grid near the layers has to be refined. Additionally, a thin impermeable layer was added just above the hydrate layer to avoid the rising of free gas from the hydrate to the upper layers.

##### 4.2. Second Scenario: Depressurization Combined with CO_{2} Injection (Case 2.1)

The same reservoir geometry and model parameters as for Case 1 are used to simulate CO_{2} storage after methane production at a single well. The results are obtained with the simulator STARS.

Within the three years of methane production at a well pressure of 30 bars, pressure and temperature in the volume decrease. This pressure decline propagates into the reservoir and causes a decomposition of methane hydrate around the well (see Figure 10). After three years, the well is closed for production and starts to function as an injection well for another three years. CO_{2} is injected at 5°C at a maximum rate of 8,000 STD m^{3}/day. As can be seen in Figures 8 and 9, pressure and temperature again start to increase close to the well due to the injection process and the exothermic formation of CO_{2} hydrate.

After the three years of depressurization, methane hydrate is decomposed completely in a radius of approximately. 10 m around the well (see Figure 10). As CO_{2} is injected at 5°C into the reservoir, methane hydrate continues decomposing within a larger area. In the region of CO_{2} hydrate formation, the - conditions are below the methane hydrate stability curve. The driving force for methane hydrate decomposition is even larger than in the case of simple depressurization. After the injection phase, the volume free of methane hydrate has a radius of approximately 60 m. In addition, while CO_{2} is injected, the gaseous methane is pushed away from the well into the (closed) reservoir, and secondary methane hydrate formation occurs because of an increase of pressure due to gas compression.

Figure 11 illustrates the development of CO_{2} hydrate saturation for the period of injection. Within the first few meters, the saturation reaches values of about 75% after three years of injection. At a distance of more than 10 m away from the well, the formation of CO_{2} hydrate is strongly limited due to the governing pressure and temperature conditions. As CO_{2} migrates into the reservoir, CO_{2} hydrate forms up to 40% within a radius of 55 m around the injection well.

Figure 12 illustrates the production and injection curves for Case 2.1 after three years.

##### 4.3. Second Scenario: Depressurization Combined with CO_{2} Injection (Case 2.2)

To simulate production of methane and storage of CO_{2} at two different wells, a Cartesian grid with an edge length of 900 m and a height of 20 m is used. The model parameters and initial conditions are the same as before (see Table 2). All simulations were performed with STARS. Starting from methane hydrate equilibrium at 72 bars (salinity is neglected) and 10°C, the bottom hole pressure in the production well is lowered to 30 bars. *Simultaneously*, CO_{2} at 10°C is injected at a maximum rate of 8,000 STD m^{3}/day at the injection well which is located 500 m apart from the production well. Both wells have a diameter of 1.0 m. The geometric configuration is shown in Figure 13.

The distributions of pressure, temperature, and hydrate saturations along the -axis of the grid in the middle layer in -direction are illustrated in Figures 14, 15, 16, and 17, respectively, for a period of six years. As far as the production well (left hand side) is concerned, the situation is the same as for Case 1. The pressure decrease (Figure 14) causes a decomposition of CH_{4} hydrate and therefore a temperature decline (Figure 15). A front propagates into the reservoir, and gas is produced from hydrates as long as the system is not at equilibrium conditions and a pressure gradient exists. The higher pressure at the injection well (right hand side) leads to the exothermic formation of CO_{2} hydrate and thus to an increase in temperature around the well (see Figure 15).

The decomposition of methane hydrate takes place at both wells (see Figure 16): at the production well because of the depressurization and at the injection well because of heat supply by the injection of CO_{2} and the heat release due to CO_{2} hydrate formation. As discussed for Case 2.1, the methane hydrate saturation declines due to depressurization within a smaller radius (10 to 20 m) than due to CO_{2} injection (approximately 100 m).

After six years, a maximum peak in CO_{2} hydrate saturation of about 75% is reached near the CO_{2} injection well. Due to unfavourable conditions for CO_{2} hydrate formation near the injection, the saturation decreases significantly within a few meters apart from the well. At a distance of about 20 m, hydrate has only been formed up to a saturation of 45%. The radius within which CO_{2} hydrate is formed is up to 150 m (see Figure 17).

An effect which can be seen in Figure 17 is clarified in Figure 18. Due to the fact that pressure is lowered at the production well, CO_{2} is sucked towards it. Therefore, more CO_{2} hydrate forms between the wells than to the right of the injection well.

In addition to the production rate for Case 2.2 (bold lines), Figure 19 shows the results for a depressurization case without CO_{2} injection (dashed lines; simulated with the same grid). After 1.5 years, the production rate with CO_{2} injection slightly exceeds the results for the case without CO_{2} injection. The reason might be that methane is pushed towards the production well while CO_{2} is injected. CO_{2} is also sucked towards the production well, but within the simulated period, the amount reaching the production well is almost zero.

#### 5. Conclusions and Outlook

First, a pure depressurization scenario for a reservoir with 1,000 m in diameter was simulated both with STARS and UMSICHT HyReS. For the given reservoir settings (see Table 2), permeability turns out to be the most important parameter for methane gas production, whereas decomposition kinetics and gas-liquid mass transfer are of lower relevance. A production rate of nearly 40,000 STD m^{3}/day can be achieved over a period of 5 years; over 10 years, the production rate decreases to half of its initial value. The reported rates can be enhanced by means of extended simulation boundaries (including over and under burdens). The results predicted by both simulators are in close agreement because the differences of the models for decomposition kinetics and gas-liquid mass transfer are of lower importance in the cases under consideration. As shown above (see Section 3.3-3.4), the enhanced hydrate formation/decomposition model used in UMSICHT HyReS depends on the concentration of the dissolved gas in the liquid phase and therefore strongly on gas-liquid mass transfer effects. Further simulation cases will be carried out with reservoir data which make mass transfer characteristics more relevant, for example, restricted gas-liquid interphase area. In such cases, model differences between STARS and UMSICHT HyReS will become more obvious.

The second simulation case aims at the industrially relevant CO_{2} sequestration problem and is set up as a one-well scenario with two consecutive steps (CH_{4} hydrate decomposition with methane production/CO_{2} hydrate formation) and as a two-well scenario with simultaneous methane production and CO_{2} sequestration. These simulations were carried out with STARS in a 2-D radial and 3-D Cartesian volume, respectively. In continuation of the pure depressurization scenario, CO_{2} is injected with a constant rate of 8,000 STD m^{3}/day in a partially depleted radial reservoir for three years. Within this period, CH_{4} hydrate has decomposed further due to unfavourable - conditions during the injection of CO_{2}. Simultaneously, CO_{2} hydrate has formed around the well. At the end of the operation period, substitution of methane hydrate by carbon dioxide hydrate has occurred only within a radius of about 60 m around the well. The saturation of CH_{4} hydrate has decreased to zero, and CO_{2} hydrate saturation has increased up to 40% and 75% just next to the well, respectively.

In the case of simultaneous CH_{4} production and CO_{2} storage in a two-well setting with a Cartesian grid, the same effects can be observed as mentioned before. Since the productivity depends on the pressure gradient driving force towards the production well, a local increase of pressure due to injection of CO_{2} leads to an enhanced gas rate (compared to the case without CO_{2} injection). A maximum rate of more than 40,000 STD m^{3}/day at the beginning of the production decreases to 35,000 STD m^{3}/day after a period of six years. Within this period, the produced amount of injected CO_{2} is almost zero.

The performed simulations show that depressurization with or without simultaneous CO_{2} injection is a practicable way to produce methane from subsea hydrate fields with assured rates over many years. Thermal stimulation is an alternative from a thermodynamic point of view, but there is a lack of appropriate technology to heat up subsea sediment in a wider range around a production well. Some promising ideas (as methane-driven point heaters, e.g.) will be investigated in subsequent work.

#### Symbols

: | Volumetric interphase area (m^{2}/m^{3}) |

: | van Leer limiter (1) |

: | Molar concentration (mole/m^{3}) |

: | Molar heat capacity (J/(mole K)) |

: | Generalized diffusion coefficient (specific) |

: | Molar total energy (J/mole) |

: | Gravity constant (m/s^{2}) |

: | Molar enthalpy (J/mole) |

: | Molar flux (mol/(m^{2}*s)) |

: | Dissociation constant (m/s) |

: | Relative permeability (1) |

: | Absolute/intrinsic permeability (m^{2}) |

: | Length of hydrate layer (m) |

: | Molar mass (mole) |

: | Molar flow per unit volume (mole/(m^{3}s)) |

: | Pressure (Pa) |

: | Capillary pressure (Pa) |

: | Heat flow per unit volume (W/(m^{3}s)) |

: | Volumetric molar conversion rate (mole/(m^{3}s)) |

: | Saturation = fraction of pore volume (m^{3}/m^{3}) |

: | Salinity = mass of salt/mass of pure water (kg/kg) |

: | Absolute temperature (K) |

: | Time (s) |

: | Velocity (volume flow per unit area), solution vector (m/s) |

: | Volume (m^{3}) |

: | Axial coordinate (m) |

: | Molar fraction (mol/mol) |

: | Vertical coordinate (m) |

β: | Mass transfer coefficient (m/s) |

: | Diffusion coefficient (m^{2}/s) |

: | Volume fraction per unit volume (m^{3}/m^{3}) |

: | Porosity (1) |

: | Permeability model parameter in STARS code (1) |

: | Dynamic viscosity (Pa*s) |

: | Brooks-Corey permeability model parameter (1) |

: | Heat conductivity (W/(m*K)) |

: | Hydrate number (1) |

: | Celsius temperature (deg C) |

: | Specific density (kg/m^{3}) |

: | Molar density (mole/m^{3}) |

σ: | Brooks-Corey permeability model parameter (1). |

#### Indices

*: | In equilibrium condition |

0: | Intrinsic/entry condition |

α: | Phase α |

abs: | Absolute |

: | Capillary |

: | Carbon dioxide |

CH: | Carbon dioxide hydrate |

: | Diffusive |

eff: | Effective |

: | Gas phase |

: | Hydrate phase |

: | Component |

: | Liquid phase |

: | Methane |

MH: | Methane hydrate |

rel: | Relative |

res: | Resident |

: | Sediment phase |

SW: | Sea water |

: | Pure water |

: | In axial direction |

: | In vertical direction. |

#### Acknowledgments

The support of our research activities by the German Federal Ministry of Economics and Technology and the German Federal Ministry of Education and Research within the framework of the SUGAR project is gratefully acknowledged.

#### References

- E. D. Sloan, “Clathrate hydrates: the other common solid water phase,”
*Industrial and Engineering Chemistry Research*, vol. 39, no. 9, pp. 3123–3129, 2000. View at Google Scholar · View at Scopus - M. Ota, Y. Abe, M. Watanabe, R. L. Smith, and H. Inomata, “Methane recovery from methane hydrate using pressurized CO
_{2},”*Fluid Phase Equilibria*, vol. 228-229, pp. 553–559, 2005. View at Publisher · View at Google Scholar · View at Scopus - M. Ota, K. Morohashi, Y. Abe, M. Watanabe, R. Lee Smith, and H. Inomata, “Replacement of CH
_{4}in the hydrate by use of liquid CO_{2},”*Energy Conversion and Management*, vol. 46, no. 11-12, pp. 1680–1691, 2005. View at Publisher · View at Google Scholar · View at Scopus - X. Zhou, S. Fan, D. Liang, and J. Du, “Replacement of methane from quartz sand-bearing hydrate with carbon dioxide-in-water emulsion,”
*Energy and Fuels*, vol. 22, no. 3, pp. 1759–1764, 2008. View at Publisher · View at Google Scholar · View at Scopus - G. Ersland, J. Husebø, A. Graue, B. A. Baldwin, J. Howard, and J. Stevens, “Measuring gas hydrate formation and exchange with CO
_{2}in Bentheim sandstone using MRI tomography,”*Chemical Engineering Journal*, vol. 158, no. 1, pp. 25–31, 2010. View at Publisher · View at Google Scholar · View at Scopus - H. Lee, Y. Seo, Y. T. Seo, I. L. Moudrakovski, and J. A. Ripmeester, “Recovering methane from solid methane hydrate with carbon dioxide,”
*Angewandte Chemie*, vol. 42, no. 41, pp. 5048–5051, 2003. View at Publisher · View at Google Scholar · View at Scopus - P. Tishchenko, C. Hensen, K. Wallmann, and C. S. Wong, “Calculation of the stability and solubility of methane hydrate in seawater,”
*Chemical Geology*, vol. 219, no. 1–4, pp. 37–52, 2005. View at Publisher · View at Google Scholar · View at Scopus - NIST Chemistry Webbook, http://webbook.nist.gov/chemistry/.
- J. M. Schicks, E. Spangenberg, R. Giese, B. Steinhauer, J. Klump, and M. Luzi, “New approaches for the production of hydrocarbons from hydrate bearing sediments,”
*Energies*, vol. 4, no. 1, pp. 151–172, 2011. View at Publisher · View at Google Scholar - N. Goel, “In situ methane hydrate dissociation with carbon dioxide sequestration: current knowledge and issues,”
*Journal of Petroleum Science and Engineering*, vol. 51, no. 3-4, pp. 169–184, 2006. View at Publisher · View at Google Scholar · View at Scopus - K. Su, C. Sun, X. Yang, G. Cgen, and S. Fan, “Experimental investigation of methane hydrate decomposition by depressurizing in porous media with 3-Dimension device,”
*Journal of Natural Gas Chemistry*, vol. 19, pp. 210–216, 2010. View at Google Scholar - G. R. Moridis and T. S. Collet, “Strategies for gas production from hydrate accumulations under various geological and reservoir conditions,” in
*Proceedings Tough Symposium*, 2003. - G. J. Moridis and E. D. Sloan, “Gas production potential of disperse low-saturation hydrate accumulations in oceanic sediments,” Tech. Rep. LBNL-5268, Lawrence Berkeley National Laboratory, Berkely, Calif, USA, 2006. View at Google Scholar
- M. Uddin, D. Coombe, and F. Wright, “Modeling of CO
^{2}-hydrate formation in geological reservoirs by injection of^{2}Gas,”*Journal of Energy Resources Technology, Transactions of the ASME*, vol. 130, no. 3, Article ID 032502, pp. 032502-1–032502-11, 2008. View at Publisher · View at Google Scholar · View at Scopus - K. C. Hester and P. G. Brewer, “Clathrate hydrates in nature,”
*Annual review of marine science*, vol. 1, pp. 303–327, 2009. View at Google Scholar - S. Hustoft, S. Bünz, J. Mienert, and S. Chand, “Gas hydrate reservoir and active methane-venting province in sediments on < 20Ma young oceanic crust in the Fram Strait, offshore NW-Svalbard,”
*Earth and Planetary Science Letters*, vol. 284, no. 1-2, pp. 12–24, 2009. View at Publisher · View at Google Scholar · View at Scopus - Z. Chen, W. Bai, W. Xu, and Z. Jin, “An analysis on stability and deposition zones of natural gas hydrate in Dongsha Region, North of South China Sea,”
*Journal of Thermodynamics*, vol. 2010, Article ID 185639, 6 pages, 2010. View at Publisher · View at Google Scholar - B. J. Anderson, M. Kurihara, M. D. White et al., “Regional long-term production modeling from a single well test, Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope,”
*Marine and Petroleum Geology*, vol. 28, no. 2, pp. 493–501, 2011. View at Publisher · View at Google Scholar - B. A. Baldwin, J. Stevens, J. J. Howard et al., “Using magnetic resonance imaging to monitor CH
_{4}hydrate formation and spontaneous conversion of CH_{4}hydrate to CO_{2}hydrate in porous media,”*Magnetic Resonance Imaging*, vol. 27, no. 5, pp. 720–726, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. W. Jung, D. N. Espinoza, and J. C. Santamarina, “Properties and phenomena relevant to CH
_{4}-CO_{2}replacement in hydrate-bearing sediments,”*Journal of Geophysical Research B*, vol. 115, no. 10, Article ID B10102, 2010. View at Publisher · View at Google Scholar · View at Scopus - M. Haeckel, E. Suess, K. Wallmann, and D. Rickert, “Rising methane gas bubbles form massive hydrate layers at the seafloor,”
*Geochimica et Cosmochimica Acta*, vol. 68, no. 21, pp. 4335–4345, 2004. View at Publisher · View at Google Scholar · View at Scopus - M. Haeckel, A. Reitz, and I. Klauke, “Methane budget of a large gas hydrate province offshore Georgia, Black Sea,” in
*Proceedings of the 6th International Conference on Gas Hydrates (ICGH '08)*, Vancouver, British Columbia, Canada, 2008. - X. Sun and K. K. Mohanty, “Kinetic simulation of methane hydrate formation and dissociation in porous media,”
*Chemical Engineering Science*, vol. 61, no. 11, pp. 3476–3495, 2006. View at Publisher · View at Google Scholar · View at Scopus - D. R. Caldwell, “Thermal conductivity of sea water,”
*Deep-Sea Research and Oceanographic Abstracts*, vol. 21, no. 2, pp. 131–137, 1974. View at Google Scholar · View at Scopus - D. J. Kukulka, B. Gebhart, and J. C. Mollendorf, “Thermodynamic and transport properties of pure and saline water,”
*Advances in Heat Transfer*, vol. 18, pp. 325–363, 1987. View at Publisher · View at Google Scholar · View at Scopus - M. Haeckel, K. Wallmann et al., “Main equations for gas hydrate modeling,”
*SUGAR Internal Communication*, 2010. View at Google Scholar - B. P. Boudreau,
*Diagenetic Models and Their Implementation: Modelling Transport and Reaction in Aquatic Sediments*, Springer, Berlin, Germany, 1997. - S. V. Pennington and M. Berzins, “New NAG library software for first-order partial differential equations,”
*ACM Transactions on Mathematical Software*, vol. 20, no. 1, pp. 63–99, 1994. View at Publisher · View at Google Scholar · View at Scopus - E. Hairer and G. Wanner,
*Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic Problems*, Springer Series in Computational Mathematics 14, Springer, 2nd edition, 1996.