An All-At-Once Newton Strategy for Marine Methane Hydrate Reservoir Models

The migration of methane through the gas hydrate stability zone (GHSZ) in the marine subsurface is characterized by highly dynamic reactive transport processes coupled to thermodynamic phase transitions between solid gas hydrates, free methane gas, and dissolved methane in the aqueous phase. The mar...

Full description

Bibliographic Details
Published in:Energies
Main Authors: Shubhangi Gupta, Barbara Wohlmuth, Matthias Haeckel
Format: Article in Journal/Newspaper
Language:English
Published: MDPI AG 2020
Subjects:
ncp
T
Online Access:https://doi.org/10.3390/en13020503
https://doaj.org/article/1cc7a08cd64944ab97027b0daa49714b
Description
Summary:The migration of methane through the gas hydrate stability zone (GHSZ) in the marine subsurface is characterized by highly dynamic reactive transport processes coupled to thermodynamic phase transitions between solid gas hydrates, free methane gas, and dissolved methane in the aqueous phase. The marine subsurface is essentially a water-saturated porous medium where the thermodynamic instability of the hydrate phase can cause free gas pockets to appear and disappear locally, causing the model to degenerate. This poses serious convergence issues for the general-purpose nonlinear solvers (e.g., standard Newton), and often leads to extremely small time-step sizes. The convergence problem is particularly severe when the rate of hydrate phase change is much lower than the rate of gas dissolution. In order to overcome this numerical challenge, we have developed an all-at-once Newton scheme tailored to our gas hydrate model, which can handle rate-based hydrate phase change coupled with equilibrium gas dissolution in a mathematically consistent and robust manner.