Mathematical Modeling and Numerical Simulation of Methane Production in a Hydrate Reservoir

• Published in: Industrial & Engineering Chemistry Research
• Main Authors: Gamwo, Isaac K., Liu, Yong
• Format: Article in Journal/Newspaper
• Language: English
• Published: American Chemical Society 2010
• Subjects: Ice; Methane hydrate; permafrost
• Online Access: https://oceanrep.geomar.de/id/eprint/36540/; https://oceanrep.geomar.de/id/eprint/36540/1/Gamwo.pdf; https://doi.org/10.1021/ie901452v

Methane hydrate, a potential future energy resource, is known to occur naturally in vast quantities beneath the ocean floor and in permafrost regions. It is important to evaluate how much methane is recoverable from these hydrate reserves. This article introduces the theoretical background of HydrateResSim, the National Energy Technology Laboratory (NETL) methane production simulator for hydrate-containing reservoirs, originally developed for NETL by Lawrence Berkeley National Laboratory (LBNL). It describes the mathematical model that governs the dissociation of methane hydrate by depressurization or thermal stimulation of the system, including the transport of multiple temperature-dependent components in multiple phases through a porous medium. The model equations are obtained by incorporating the multiphase Darcy’s law for gas and liquid into both the mass component balances and the energy conservation equations. Two submodels in HydrateResSim for hydrate dissociation are also considered: a kinetic model and a pure thermodynamic model. Contrary to more traditional reservoir simulations, the set of model unknowns or primary variables in HydrateResSim changes throughout the simulation as a result of the formation or dissociation of ice and hydrate phases during the simulation. The primary variable switch method (PVSM) is used to effectively track these phase changes. The equations are solved by utilizing the implicit time finite-difference method on the grid system, which can properly describe phase appearance or disappearance as well as the boundary conditions. The Newton-Raphson method is used to solve the linear equations after discretization and setup of the Jacobian matrix. We report here the application of HydrateResSim to a three-component, four-phase flow system in order to predict the methane produced from a laboratory-scale reservoir. The first results of HydrateResSim code in a peer-reviewed publication are presented in this article. The numerical solution was verified against the state-of-the art ...