Description: 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 occurring 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 poroelasticity. The closely coordinated experimental and numerical procedures enable us to validate the proposed simplified geomechanics-to-flow coupling, and set an important precursor toward enhancing our coupled hydro-geomechanical hydrate reservoir simulator with more suitable elastoplastic constitutive models.

Description: Burial driven recycling is an important process in the natural gas hydrate (GH) systems worldwide, characterized by complex multiphysics interactions like gas migration through an evolving gas hydrate stability zone (GHSZ), competing gas-water-hydrate (i.e. fluid-fluid-solid) phase transitions, locally appearing and disappearing phases, and evolving sediment properties (e.g., permeability, reaction surface area, and capillary entry pressure). Such a recycling process is typically studied in homogeneous or layered sediments. However, there is mounting evidence that structural heterogeneity and anisotropy linked to normal and inclined fault systems or anomalous sediment layers have a strong impact on the GH dynamics. Here, we consider the impacts of such a structurally complex media on the recycling process. To capture the properties of the anomalous layers accurately, we introduce a fully mass conservative, high-order, discontinuous Galerkin (DG) finite element based numerical scheme. Moreover, to handle the rapidly switching thermodynamic phase states robustly, we cast the problem of phase transitions as a set of variational inequalities, and combine our DG discretization scheme with a semismooth Newton solver. Here, we present our new simulator, and demonstrate using synthetic geological scenarios, a) how the presence of an anomalous high-permeability layer, like a fracture or brecciated sediment, can alter the recycling process through flow-localization, and more importantly, b) how an incorrect or incomplete approximation of the properties of such a layer can lead to large errors in the overall prediction of the recycling process.

Description: Methane gas hydrates have increasingly become a topic of interest because of their potential as a future energy resource. There are significant economical and environmental risks associated with extraction from hydrate reservoirs, so a variety of multiphysics models have been developed to analyze prospective risks and benefits. These models generally have a large number of empirical parameters which are not known a priori. Traditional optimization-based parameter estimation frameworks may be ill-posed or computationally prohibitive. Bayesian inference methods have increasingly been found effective for estimating parameters in complex geophysical systems. These methods often are not viable in cases of computationally expensive models and high-dimensional parameter spaces. Recently, methods have been developed to effectively reduce the dimension of Bayesian inverse problems by identifying low-dimensional structures that are most informed by data. Active subspaces is one of the most generally applicable methods of performing this dimension reduction. In this paper, Bayesian inference of the parameters of a state-of-the-art mathematical model for methane hydrates based on experimental data from a triaxial compression test with gas hydrate-bearing sand is performed in an efficient way by utilizing active subspaces. Active subspaces are used to identify low-dimensional structure in the parameter space which is exploited by generating a cheap regression-based surrogate model and implementing a modified Markov chain Monte Carlo algorithm. Posterior densities having means that match the experimental data are approximated in a computationally efficient way.

Description: The presence of gas hydrates (GHs) increases the stiffness and strength of marine sediments. In elasto‐plastic constitutive models, it is common to consider GH saturation (Sh) as key internal variable for defining the contribution of GHs to composite soil mechanical behavior. However, the stress‐strain behavior of GH‐bearing sediments (GHBS) also depends on the microscale distribution of GH and on GH‐sediment fabrics. A thorough analysis of GHBS is difficult, because there is no unique relation between Sh and GH morphology. To improve the understanding of stress‐strain behavior of GHBS in terms of established soil models, this study summarizes results from triaxial compression tests with different Sh, pore fluids, effective confining stresses, and strain histories. Our data indicate that the mechanical behavior of GHBS strongly depends on Sh and GH morphology, and also on the strain‐induced alteration of GH‐sediment fabrics. Hardening‐softening characteristics of GHBS are strain rate‐dependent, which suggests that GH‐sediment fabrics dynamically rearrange during plastic yielding events. We hypothesize that rearrangement of GH‐sediment fabrics, through viscous deformation or transient dissociation and reformation of GHs, results in kinematic hardening, suppressed softening, and secondary strength recovery, which could potentially mitigate or counteract large‐strain failure events. For constitutive modeling approaches, we suggest that strain rate‐dependent micromechanical effects from alterations of the GH‐sediment fabrics can be lumped into a nonconstant residual friction parameter. We propose simple empirical evolution functions for the mechanical properties and calibrate the model parameters against the experimental data. Plain Language Summary Gas hydrates (GHs) are crystalline‐like solids, which are formed from natural gas molecules and water at high pressure and low temperature. GHs, and particularly methane hydrates, are naturally abundant in marine sediments. It is known that the presence of GH increases the mechanical stiffness and strength of sediments, and there is strong effort in analyzing and quantifying these effects in order to understand potential risks of sediment destabilization or slope failure. Based on our experimental results from high‐pressure geotechnical studies, we show that not only the initial amount and distribution of GH are important for the increased strength of GH‐bearing sediments but also the dynamic rearrangement of GH‐sediment fabrics during deformation characterizes the stress‐strain response and enables strength recovery after failure. We propose that different microstructural mechanisms contribute to this rearrangement and strength recovery of GH sediment. However, we consider these complicated processes in a simplified manner in an improved numerical model, which can be applied for geotechnical risk assessment on larger scales.

Description: Abrupt fluid emissions from shallow marine sediments pose a threat to seafloor installations like wind farms and offshore cables. Quantifying such fluid emissions and linking pockmarks, the seafloor manifestations of fluid escape, to flow in the sub-seafloor remains notoriously difficult due to an incomplete understanding of the underlying physical processes. Here, using a compositional multi-phase flow model, we test plausible gas sources for pockmarks in the south-eastern North Sea, which recent observations suggest have formed in response to major storms. We find that the presence of free gas in the subsurface effectively damps storm wave-induced pressure changes due to its high compressibility, so that the mobilization of pre-existing gas pockets is unlikely. Rather, our results point to spontaneous appearance of a free gas phase via storm-induced gas exsolution from pore fluids. This mechanism is primarily driven by the pressure-sensitivity of gas solubility. We show that in highly permeable sediments containing gas-rich pore fluids, wave-induced pressure changes result in the appearance of a persistent gas phase. This suggests that seafloor fluid escape structures are not always proxies for overpressured shallow gas and that periodic seafloor pressure changes can induce persistent free gas phase to spontaneously appear.

Description: Offshore meteoric groundwater (OMG) has long been hypothesised to be a driver of seafloor geomorphic processes in continental margins worldwide. Testing this hypothesis has been challenging because of our limited understanding of the distribution and rate of OMG flow and seepage, and their efficacy as erosive/destabilising agents. Here we carry out numerical simulations of groundwater flow and slope stability using conceptual models and evolving stratigraphy - for passive siliciclastic and carbonate margin cases – to assess whether OMG and its evolution during a late Quaternary glacial cycle can generate the pore pressures required to trigger mechanical instabilities on the seafloor. Conceptual model results show that mechanical instabilities by OMG flow are most likely to occur in the outer shelf to upper slope, at or shortly before the Last Glacial Maximum sea level lowstand. Models with evolving stratigraphy show that OMG flow is a key driver of pore pressure development and instability in the carbonate margin case. In the siliciclastic margin case, OMG flow plays a secondary role in preconditioning the slope to failure. The higher degree of spatial/stratigraphic heterogeneity of carbonate margins, lower shear strengths of their sediments, and limited generation of overpressures by sediment loading may explain the higher susceptibility of carbonate margins, in comparison to siliciclastic margins, to mechanical instability by OMG flow. OMG likely played a more significant role in carbonate margin geomorphology (e.g. Bahamas, Maldives) than currently thought. Key Points Offshore meteoric groundwater (OMG) flow can drive mechanical instabilities in the outer shelf to upper slope Such instabilities occur at, or shortly after, the Last Glacial Maximum sea level lowstand Carbonate margins are more susceptible to mechanical instability by OMG than siliciclastic margins

Description: Subsurface flows, particularly hyporheic exchange fluxes, driven by streambed topography, permeability, channel gradient and dynamic flow conditions provide prominent ecological services such as nitrate removal from streams and aquifers. Stream flow dynamics cause strongly nonlinear and often episodic contributions of nutrient concentrations in river-aquifer systems. Using a fully coupled transient flow and reactive transport model, we investigated the denitrification potential of hyporheic zones during peak-flow events. The effects of streambed permeability, channel gradient and bedform amplitude on the spatio-temporal distribution of nitrate and dissolved organic carbon in streambeds and the associated denitrification potential were explored. Distinct peak-flow events with different intensity, duration and hydrograph shape were selected to represent a wide range of peak-flow scenarios. Our results indicated that the specific hydrodynamic characteristics of individual flow events largely determine the average positive or negative nitrate removal capacity of hyporheic zones, however the magnitude of this capacity is controlled by geomorphological settings (i.e. channel slope, streambed permeability and bedform amplitude). Specifically, events with longer duration and higher intensity were shown to promote higher nitrate removal efficiency with higher magnitude of removal efficiency in the scenarios with higher slope and permeability values. These results are essential for better assessment of the subsurface nitrate removal capacity under the influence of flow dynamics and particularly peak-flow events in order to provide tailored solutions for effective restoration of interconnected river-aquifer systems.

Description: Abrupt fluid emissions from shallow marine sediments pose a threat to seafloor installations like wind farms and offshore cables. Quantifying such fluid emissions and linking pockmarks, the seafloor manifestations of fluid escape, to flow in the sub-seafloor remains notoriously difficult due to an incomplete understanding of the underlying physical processes. Here, using a compositional multi-phase flow model, we test plausible gas sources for pockmarks in the south-eastern North Sea, which recent observations suggest have formed in response to major storms. We find that the mobilization of pre-existing gas pockets is unlikely because free gas, due to its high compressibility, damps the propagation of storm-induced pressure changes deeper into the subsurface. Rather, our results point to spontaneous appearance of a free gas phase via storm-induced gas exsolution from pore fluids. This mechanism is primarily driven by the pressure-sensitivity of gas solubility, and the appearance of free gas is largely confined to sediments in the vicinity of the seafloor. We show that in highly permeable sediments containing gas-rich pore fluids, wave-induced pressure changes result in the appearance of a persistent gas phase. This suggests that seafloor fluid escape structures are not always proxies for overpressured shallow gas and that periodic seafloor pressure changes can induce persistent free gas phase to spontaneously appear. Key Points - Storm-induced pressure changes can lead to spontaneous appearance of free gas phase near the seafloor - This process is driven by pressure-sensitive phase instabilities - This mechanism could help explain elusive gas sources in recently observed pockmarks in the North Sea

Description: Focused fluid flow through sub-seafloor pipes and chimneys, and their seafloor manifestations as pockmarks, are ubiquitous. However, the dynamics of flow localization and evolution of fluid escape structures remain poorly understood. Models based on geomechanical mechanisms like hydro-fracturing and porosity wave propagation offer some useful insights into fluid flow and escape dynamics, but face limitations in capturing features like mobilized granular matter, especially in the upper sediment layers where the link between fracture and pockmark is not always clear. Here, we propose a mathematical model based on the multiphase theory of porous media, where changes in subsurface and seafloor morphology are resolved through seepage-induced erosion, fluidization, transport, and re-deposition of granular material. Through simulation of an idealized scenario of gas escape from overpressured shallow gas reservoir, we demonstrate that our model can capture flow localization and formation of pipes, chimneys, and pockmarks. Our simulations show (1) formation of conical focused-flow conduits with a brecciated core and annular gas channels; (2) pockmarks of W and ring shapes; and (3) pulsed release of gas. Sediment erodibility and flow anisotropy control the morphology of focused fluid flow and escape structures, while permeability shows negligible impact. While the geological setting for this study is theoretical, we show that our results have real-world analogs.

Description: Burial driven recycling is an important process in the natural gas hydrate (GH) systems worldwide, characterized by complex multiphysics interactions like gas migration through an evolving gas hydrate stability zone (GHSZ), competing gas-water-hydrate (i.e., fluid-fluid-solid) phase transitions, locally appearing and disappearing phases, and evolving sediment properties (e.g., permeability, reaction surface area, and capillary entry pressure). Such a recycling process is typically studied in homogeneous or layered sediments. However, there is mounting evidence that structural heterogeneity and anisotropy linked to normal and inclined fault systems or anomalous sediment layers have a strong impact on the GH dynamics. Here, we consider the impacts of such a structurally complex media on the recycling process. To capture the properties of the anomalous layers accurately, we introduce a fully mass conservative, high-order, discontinuous Galerkin (DG) finite element based numerical scheme. Moreover, to handle the rapidly switching thermodynamic phase states robustly, we cast the problem of phase transitions as a set of variational inequalities, and combine our DG discretization scheme with a semi-smooth Newton solver. Here, we present our new simulator, and demonstrate using synthetic geological scenarios, (a) how the presence of an anomalous high-permeability layer, like a fracture or brecciated sediment, can alter the recycling process through flow-localization, and more importantly, (b) how an incorrect or incomplete approximation of the properties of such a layer can lead to large errors in the overall prediction of the recycling process. Key Points Structural heterogeneity linked to inclined fault systems or anomalous sediment layers have a strong impact on the gas hydrate dynamics The presence of anomalous high-permeability layers within gas hydrate stability zone alters the recycling process through flow-localization The presented discontinuous Galerkin scheme is able to accurately capture the gas hydrate recycling processes through strongly anisotropic materials