LBL logo Lawrence Berkeley National Laboratory » Earth Sciences Division
Questions & Comments | Privacy & Security Notice

TOUGH+ Software


TOUGH+ is the result of a re-engineering effort in which the capabilities of TOUGH2 have been recast in a modular design that adheres to the tenets of object-oriented programming as implemented in Fortran 95. New capabilities have been added. I particular, the range of thermodynamic properties of water has been extended to include ice for the modeling of freezing and thawing in permafrost regions, and - specifically - for the simulation of processes in hydrate-bearing sediments, which involve. formation and dissociation of ice-like hydrates with associated changes in porosity, permeability, and two-phase flow characteristics. The first module released in the new framework is TOUGH+Hydrate for the simulation of these processes in hydrate-bearing formations. TOUGH+ maintains backwards compatibility with TOUGH2 input files.

Features & Capabilities

TOUGH+ solves mass and energy balance equations that describe fluid and heat flow in general multiphase, multicomponent, porous-media systems. TOUGH+ can model the following processes and phenomena in geologic systems:

  1. Conduction
  2. Advection/convection
  3. Radiation
  4. Latent heat related to phase changes
  5. Gas and inhibitor dissolution
  6. The transport of water-soluble gases and solutes (such as salts and alcohols, if included in the properties model), accounting for advection and molecular diffusion

When using the HYDRATE properties module:

A deliberate effort was made to keep the simplifying assumptions involved in the development of the underlying physical, thermodynamic and mathematical model to a minimum. The assumptions used include:

A detailed description of the constitutive equations used in TOUGH+ can be found in the TOUGH+Hydrate Users Guide.

TOUGH+ maintains backward compatibility with TOUGH2 input files and input conventions, while also providing new user features such as dynamic memory allocation and modern memory management, enhanced meshing capabilities, namelist-based input designations, improved labeling of input and output files, and the option of coupling to commercial geomechanical simulators.

TOUGH+ has been validated and verified through comparison to different analytical and numerical solutions, laboratory results, and data from field-scale experiments.


Gas Production from Oceanic Class 2 Hydrate Deposits

George J. Moridis and Matthew T. Reagan

Production from gas hydratesClass 2 hydrate deposits comprise two zones: a hydrate-bearing layer overlying a zone of mobile water. The objective of this study is to evaluate the production potential of marine Class 2 methane hydrate accumulations by means of numerical simulation, and to develop appropriate strategies for gas production. The simulated system is based the “Tigershark” area, located in the Gulf of Mexico. We use a 2-D radial mesh representing the hydrate deposit, underlying aquifer, and impermeable shale boundaries. Production simulations were performed using TOUGH+Hydrate. Depressurization-induced dissociation was selected as the most promising gas production strategy because of its simplicity and its technical and economic effectiveness.

For all strategies, production at the well is characterized by distinct cycles (see figure), with the end of each cycle marked by cavitation (rapid pressure drop at the well) caused by either secondary hydrate blockages near the well or by the increasing participation of lower-density gas. This is remedied by warm seawater injection and a reduction in the total mass production rate. Gas production continuously increases during each cycle, with a corresponding reduction in the water production.

The most effective production strategy used a combination of electrical heating to initiate production, warm seawater circulation during continued production, and the application of a novel well design that combines warm seawater injection with a long production interval to access the remaining hydrate and avoid blockages and cavitation (see the related publication). This approach is clearly superior to conventional production techniques for production under the Tigershark reference conditions, and demonstrates that gas can be produced at high rates over long periods from Class 2 oceanic hydrate deposits using conventional technology.

The figure shows rates of (a) hydrate-originating CH4 release in the reservoir (red) and (b) CH4 production at the well (blue) during production from a Class 2 oceanic hydrate deposit. The stages of production and the average production rate (Qavg, in green) over the simulation period are also shown.


Gas Production from Oceanic Class 3 Hydrate Deposits

George J. Moridis and Matthew T. Reagan

Gas Production from Oceanic Class 3 Hydrate Deposits

Class 3 hydrate accumulations are composed of a single zone of methane hydrate bounded by confining strata with no underlying zone of mobile fluids. Early production attempts using thermal stimulation suggested these formations as poor targets for gas production. A new strategy, simulated using TOUGH+HYDRATE, uses constant-pressure production to overcome the very low initial effective permeability of the deposit and the lack of a zone of mobile fluids.

The simulated gas production occurs in a cyclical pattern. Each cycle consists of a long stage of increasing gas production, followed by sharp production decline, followed by a recovery. This pattern is repeated until the cessation of production. Peak production regularly exceeds 5 ST m3/s of CH4 (15 MMSCFD), and the average gas production during the 6000-day simulation period is 2.61 ST m3/s (8.10 MMSCFD). Recovery from each production decline occurs naturally, without heating or fluid injection.

The figure shows the evolution of the spatial distribution of hydrate saturation over time, and provides an explanation for the cyclical pattern. The precipitous drops in production are caused by the evolution of traveling barriers of secondary hydrate around the well. As gas flows to the well, a cylindrical sheath of secondary hydrate forms next to the wellbore (Panel A). Because the inner radius of the barrier is exposed to intense depressurization, it dissociates, but additional hydrate forms on its outer radius, resulting in the appearance of a traveling barrier (B). A second barrier forms (C) and moves outward through the same process (D). The outer barrier is now shielded from steep pressure gradients, and dissociates (D through F). New barriers appear and disappear (G though K) over time. Each of the temporary drops in production occur when a new inner secondary hydrate barrier is formed.

This simulation demonstrates that large volumes of natural gas can be produced at high rates over long periods from Class 3 oceanic hydrate deposits using straightforward depressurization-induced dissociation. Hence, these deposits may be viable candidates for future exploitation and development.

The figure shows the evolution of spatial distribution of hydrate saturation, SH, during gas production from a Class 3 oceanic hydrate deposit, showing the traveling barriers of secondary hydrate.


Oceanic Gas Hydrate Instability and Dissociation Under Climate Change Scenarios

Matthew T. Reagan and George J. Moridis

Schematic of gas hydrate phase diagramThe amount of methane hydrate currently residing in the deep ocean and along continental margins is estimated to be from 500 billion to over 3,000 billion tons. The dissociation of accumulated oceanic gas hydrate deposits and the release of large quantities of methane—a powerful greenhouse gas—in a short period of time could have dramatic climatic consequences, leading to positive climate feedbacks. Using TOUGH+Hydrate, we simulated the response of several types of oceanic gas hydrate accumulations to temperature changes at the seafloor and assessed the potential for methane release into the ecosystem.

The figures shows a schematic of the gas hydrate stability zone in the seafloor environment and its relationship to ocean temperature changes.
The figure presents a general schematic of the hydrate phase boundary (blue line) and gas hydrate stability zone (GHSZ) for oceanic hydrates (shaded area). An increase in water temperature at the seafloor (a shift from temperature profile (1) to profile (2) lowers the position of the top of the GHSZ (A) and raises the position of the bottom of the GHSZ (B). Such a shift could induce hydrate dissociation and lead to the release of methane into the ocean and atmosphere.

We modeled oceanic hydrate deposits, represented as 1-D columns, under a range of depths and temperatures using the TOUGH+HYDRATE code. We simulated temperature increases of 1 °C to 5 °C over a 100 yr period, representing possible changes in seafloor temperature as predicted by climate simulators, and determined the fate of the hydrate deposits and tracked the transport of released methane to the seafloor.

These results indicate that while many deep oceanic hydrate deposits are stable under the influence of significant temperature variations, shallow deposits, such as those found in arctic regions or in the Gulf of Mexico, can undergo rapid dissociation and release large quantities of carbon. In arctic regions in particular, temperature changes are expected to be more pronounced, increasing the risk of hydrate destabilization even further. These simulations will provide a source term to global climate models, yielding a prediction of the possible effects of oceanic hydrate decomposition on global climate.


Awards & Recognition

Licensing & Download


TOUGH+Hydrate Users Guide