Matthias Haeckel

Shubhangi Gupta, Barbara Wohlmuth, Matthias Haeckel

Marine gas hydrate systems are characterized by highly dynamic transport-reaction processes in an essentially water-saturated porous medium that are coupled to thermodynamic phase transitions between solid gas hydrates, free gas and dissolved methane in the aqueous phase. These phase transitions are highly nonlinear and strongly coupled, and cause the mathematical model to rapidly switch the phase states and pose serious convergence issues for the classical Newton's method. One of the common methods of dealing with such phase transitions is the primary variable switching (PVS) method where the choice of the primary variables is adapted locally to the phase state **outside** the Newton loop. In order to ensure that the phase states are determined accurately, the PVS strategy requires an additional iterative loop, which can get quite expensive for highly nonlinear problems. For methane hydrate reservoir models, the PVS method shows poor convergence behaviour and often leads to extremely small time step sizes. In order to overcome this issue, we have developed a nonlinear complementary constraints method (NCP) for handling phase transitions, and implemented it within a non-smooth Newton's linearization scheme using an active-set strategy. Here, we present our numerical scheme and show its robustness through field scale applications based on the highly dynamic geological setting of the Black Sea.

Shubhangi Gupta, Christian Deusner, Matthias Haeckel et al.

Natural gas hydrates are considered a potential resource for gas production on industrial scales. Gas hydrates contribute to the strength and stiffness of the hydrate-bearing sediments. During gas production, the geomechanical stability of the sediment is compromised. Due to the potential geotechnical risks and process management issues, the mechanical behavior of the gas hydrate-bearing sediments needs to be carefully considered. In this study, we describe a coupling concept that simplifies the mathematical description of the complex interactions occuring during gas production by isolating the effects of sediment deformation and hydrate phase changes. Central to this coupling concept is the assumption that the soil grains form the load-bearing solid skeleton, while the gas hydrate enhances the mechanical properties of this skeleton. We focus on testing this coupling concept in capturing the overall impact of geomechanics on gas production behavior though numerical simulation of a high-pressure isotropic compression experiment combined with methane hydrate formation and dissociation. We consider a linear-elastic stress-strain relationship because it is uniquely defined and easy to calibrate. Since, in reality, the geomechanical response of the hydrate bearing sediment is typically inelastic and is characterized by a significant shear-volumetric coupling, we control the experiment very carefully in order to keep the sample deformations small and well within the assumptions of poro-elasticity. The closely co-ordinated experimental and numerical procedures enable us to validate the proposed simplified geomechanics-to-flow coupling, and set an important precursor towards enhancing our coupled hydro-geomechanical hydrate reservoir simulator with more suitable elasto-plastic constitutive models.