TOUGH + Millstone has been developed for the analysis of coupled flow, thermal and geomechanical processes associated with the formation and/or dissociation of CH4-hydrates in geological media. It is composed of two constituent codes: (a) a significantly enhanced version of the TOUGH + HYDRATE simulator, V2.0, that accounts for all known flow, physical, thermodynamic and chemical processes associated with the behavior of hydrate-bearing systems undergoing changes and includes the most recent advances in the description of the system properties, coupled seamlessly with (b) Millstone V1.0, a new code that addresses the conceptual, computational and mathematical shortcomings of earlier codes used to describe the geomechanical response of these systems. The capabilities of TOUGH + Millstone are demonstrated in the simulation and analysis of the system flow, thermal and geomechanical behavior during gas production from a realistic complex offshore hydrate deposit. In the first paper of this series, we discuss the physics underlying the T + H hydrate simulator, the constitutive relationships describing the physical, chemical (equilibrium and kinetic) and thermal processes, the states of the CH 4+ H 2O system and the sources of critically important data, as well as the mathematical approaches used for the development of the of mass and energy balance equations and their solution. Additionally, we provide verification examples of the hydrate code against numerical results from the simulation of laboratory and field experiments.

A recent Department of Energy field test on the Alaska North Slope has increased interest in the ability to simulate systems of mixed CO2-CH4 hydrates. However, the physically realistic simulation of mixed-hydrate simulation is not yet a fully solved problem. Limited quantitative laboratory data leads to the use of various ab initio, statistical mechanical, or other mathematic representations of mixed-hydrate phase behavior. Few of these methods are suitable for inclusion in reservoir simulations, particularly for systems with large number of grid elements, 3D systems, or systems with complex geometric configurations. In this work, we present a set of fast parametric relationships describing the thermodynamic properties and phase behavior of a mixed methane-carbon dioxide hydrate system. We use well-known, off-the-shelf hydrate physical properties packages to generate a sufficiently large dataset, select the most convenient and efficient mathematical forms, and fit the data to those forms to create a physical properties package suitable for inclusion in the TOUGH+ family of codes. The mapping of the phase and thermodynamic space reveals the complexity of the mixed-hydrate system and allows understanding of the thermodynamics at a level beyond what much of the existing laboratory data and literature currently offer.

We investigate by means of numerical simulation a planned year-long field test of depressurization-induced production from a permafrost-associated hydrate reservoir on the Alaska North Slope at the site of the recently drilled Hydrate-01 Stratigraphic Test Well. The main objective of this study is to assess quantitatively the impact of temporary interruptions (well shut-ins) on the expected fluid production performance from the B1 Sand of the stratigraphic Unit B during controlled depressurization over different time scales, as well as on other relevant aspects of the system response that have the potential to significantly affect the design of the field test. We consider eight different cases of depressurization, including (a) rapid depressurization over a 60-day period to a terminal bottomhole pressure PWof 2.8 MPa and (b) a slower depressurization rate to a final PWof 0.6 MPa at the end of the year-long production test, in addition to (c) a multi-step depressurization regime and (d) a quasi-linear continuous depressurization strategy. The results of the study indicate that shut-ins obviously reduce gas release and production during and immediately after their occurrence, but their longer-term effects are strongly dependent on the depressurization regime and on the time of observation, covering the entire range of potential outcomes. Shut-ins (a) have a universally strong negative effect when quasi-linear depressurization is involved regardless of the length of the production period, and (b) have a strong positive effect in multi-step depressurization schemes that becomes apparent earlier for large initial pressure drops, but (c) can also appear to have practically no effect for slow stepwise depressurization at the end of the year-long production test. Shut-ins lead to a rapid reformation of hydrates, even to the point of disappearance of a free gas phase in the reservoir. Rapid depressurization regimes lead to early maximum rates of hydrate dissociation and gas production, while the maximum rates occur at the end of the production test for the cases of slower depressurization. Shut-ins do not appear to have a significant impact on water production, as the cessation of production is followed by higher rates production when depressurization resumes. Similarly, (a) the fraction of produced CH4originating from exsolution from the water, (b) the water-to-gas ratio, and (c) the rate of replenishment of produced water by boundary inflows do not appear significantly affected by shut-ins, the effects of which seem to be temporary in the majority of the cases. The study confirmed the superiority of multi-step depressurization methods as the most effective strategies for hydrate dissociation and gas production and showed that two observation wells (located at distances of 30 and 50 m from the production well) are appropriately positioned and both able to capture the P, T, and SGbehavior during the fluid production and shut-ins in any of the eight cases we investigated.

The TOUGH+Millstone simulator has been developed for the analysis of coupled flow, thermal and geomechanical processes associated with the formation and/or dissociation of CH 4 hydrates in geological media. It is composed of two constituent codes: (a) a significantly enhanced version of the TOUGH+HYDRATE simulator, v2.0, that accounts for all known flow, physical, thermodynamic and chemical processes associated with the behavior of hydrate-bearing systems undergoing changes and includes the most recent advances in the description of the system properties, coupled seamlessly with (b) Millstone v1.0, a new code that addresses the conceptual, computational and mathematical shortcomings of earlier codes used to describe the geomechanical response of these systems. The capabilities of the TOUGH+Millstone code are demonstrated in the simulation and analysis of the system flow, thermal and geomechanical behavior during gas production from a realistic complex offshore hydrate deposit. In the second part of this series, we describe the Millstone geomechanical simulator. The hydrate-dependent, rate-based poromechanical formulation is presented and solved using a finite element discretization. A novel multimesh coupling scheme is introduced, wherein interpolators are automatically built to transfer data between the finite difference discretization of TOUGH+ and the finite element discretization of Millstone. We provide verification examples against analytic solutions for poroelasticity and a simplified demonstration problem for mechanically induced phase change in a hydrate sediment.

The TOUGH+Millstone simulator has been developed for the analysis of coupled flow, thermal and geomechanical processes associated with the formation and/or dissociation of CH 4-hydrates in geological media. It is composed of two constituent codes: (a) a significantly enhanced version of the TOUGH+Hydrate simulator, v2.0, that accounts for all known flow, physical, thermodynamic and chemical processes associated with the evolution of hydrate-bearing systems and includes the most recent physical properties relationships, coupled seamlessly with (b) Millstone v1.0, a new code that addresses the conceptual, computational and mathematical shortcomings of earlier codes used to describe the geomechanical response of these systems. The capabilities of the TOUGH+Millstone code are demonstrated in the simulation and analysis of the system flow, thermal, and geomechanical behavior during gas production from a realistic complex offshore hydrate deposit. In the third paper of this series, we apply the simulators described in parts 1 and 2 to a problem of gas production from a complex, multilayered system of hydrate-bearing sediments in an oceanic environment. We perform flow simulations of constant-pressure production via a vertical well and compare those results to a coupled flow-geomechanical simulation of the same process. The results demonstrate the importance of fully coupled geomechanics when modeling the evolution of reservoir properties during production.

We investigated numerical simulation strategies for a long-term test of depressurization-induced gas production from the B1 Sand of Unit B at the Hydrate-01 Stratigraphic Test Well. The main objective of this study was to estimate fluid production rates (with emphasis on water production) under a variety of conditions and production scenarios and contribute new insights to the design and management of the field test. In the first part of the study, we investigated the system response to a three-step depressurization process using two limiting sets of flow properties─the expected maximum and minimum intrinsic and effective permeabilities─for the very heterogeneous reservoir. In the second part, we investigated the effect of the production interval length and placement within the formation relative to the boundaries of the hydrate-bearing unit. The best performing well configuration was used in the third part of the study, which used the most representative subsurface flow properties to investigate the effect of the depressurization strategy on the production performance. The best overall performance (largest gas production with modest water production and a strong response at the observation wells) was obtained with a 10 m-long well situated 3 m below the top of the formation and a three-step depressurization scheme at 15-day intervals to a terminal bottomhole pressure of 2.8 MPa. The overall production performance was enhanced by a faster rate of depressurization. Estimated water production rates in all cases were limited and easily manageable. None of the tested well configurations or depressurization strategies significantly reduced water production without also severely reducing gas production. In all the investigated cases, 95% of the long-term fraction of produced water was replenished by inflows from the boundaries and could not be reduced. These substantial water inflows are an unavoidable feature of HU-B and cannot be easily mitigated by a hydraulic control.

The objective of this study is to quantify, by means of numerical simulation, the response of the complex system of gas hydrate accumulations at Site NGHP-02-09, Krishna-Godavari Basin, Indian Ocean, to different production conditions, and to determine the technical feasibility of gas production through depressurization-induced dissociation. The study assesses the suitability of the site for a long-term production test involving a single vertical well, and the long-term potential of the deposit under full-field production using a system of multiple vertical wells. We simulate gas and water flow, estimate the production performance of the accumulation and separately investigate the corresponding geomechanical response of the system. Results indicate that production from Site NGHP-02-09 under the conditions of a long-term field test involving a single vertical well is technically feasible and can yield high gas production rates. However, an inability to fully isolate the water bearing zones results in production that is largely from dissolved gas rather than hydrate dissociation and is thus burdened by excessive water production. Given the estimated physical properties of the reservoir system, Site NGHP-02-09 does not appear to be a promising location for a single-well field test of gas production, but may be a promising production target for full-field production operations using a multi-well system in which exterior wells can mitigate water inflows to allow interior wells to more effectively depressurize the formation and capture methane from gas hydrate dissociation. Geomechanical issues need to be carefully considered as significant displacements are possible, which can be challenging to well construction and stability.