Description: Highlights • Code comparisons build confidence in simulators to model interdependent processes. • International hydrate reservoir simulators are compared over five complex problems. • Geomechanical processes significantly impact response of gas hydrate reservoirs. • Simulators yielded comparable results, however many differences are noted. • Equivalent constitutive models are required to achieve agreement across simulators. Geologic reservoirs containing gas hydrate occur beneath permafrost environments and within marine continental slope sediments, representing a potentially vast natural gas source. Numerical simulators provide scientists and engineers with tools for understanding how production efficiency depends on the numerous, interdependent (coupled) processes associated with potential production strategies for these gas hydrate reservoirs. Confidence in the modeling and forecasting abilities of these gas hydrate reservoir simulators (GHRSs) grows with successful comparisons against laboratory and field test results, but such results are rare, particularly in natural settings. The hydrate community recognized another approach to building confidence in the GHRS: comparing simulation results between independently developed and executed computer codes on structured problems specifically tailored to the interdependent processes relevant for gas hydrate-bearing systems. The United States Department of Energy, National Energy Technology Laboratory, (DOE/NETL), sponsored the first international gas hydrate code comparison study, IGHCCS1, in the early 2000s. IGHCCS1 focused on coupled thermal and hydrologic processes associated with producing gas hydrates from geologic reservoirs via depressurization and thermal stimulation. Subsequently, GHRSs have advanced to model more complex production technologies and incorporate geomechanical processes into the existing framework of coupled thermal and hydrologic modeling. This paper contributes to the validation of these recent GHRS developments by providing results from a second GHRS code comparison study, IGHCCS2, also sponsored by DOE/NETL. IGHCCS2 includes participants from an international collection of universities, research institutes, industry, national laboratories, and national geologic surveys. Study participants developed a series of five benchmark problems principally involving gas hydrate processes with geomechanical components. The five problems range from simple geometries with analytical solutions to a representation of the world's first offshore production test of methane hydrates, which was conducted with the depressurization method off the coast of Japan. To identify strengths and limitations in the various GHRSs, study participants submitted solutions for the benchmark problems and discussed differing results via teleconferences. The GHRSs evolved over the course of IGHCCS2 as researchers modified their simulators to reflect new insights, lessons learned, and suggested performance enhancements. The five benchmark problems, final sample solutions, and lessons learned that are presented here document the study outcomes and serve as a reference guide for developing and testing gas hydrate reservoir simulators.

Description: Highlights • Sedimentation-driven gas hydrate recycling is cyclic in nature with time scales set by reactive multi-phase transport. • Each cycle can be divided into three distinct phases: 1) gas accumulation phase, 2) gas breakthrough phase and 3) uninhibited hydrate build-up phase. • In the presence of sufficient accumulated gas, convex deposition of hydrate acts like a mechanical nozzle for the ascending gas flow. Gas hydrate recycling is an important process in natural hydrate systems worldwide and frequently leads to the high gas hydrate saturations found close to the base of the gas hydrate stability zone (GHSZ). However, to date it remains enigmatic how, and under which conditions, free gas invades back into the GHSZ. Here we use a 1D compositional multi-phase flow model that accounts for sedimentation to investigate the dominant mechanisms that control free gas flow into the GHSZ using a wide-range of parameters i.e. hydrate formation kinetics, sediment permeability, and capillary pressure. In the first part of this study, we investigate free gas invasion into the GHSZ without any sedimentation, and analyse the dynamics of hydrate formation in the vicinity of the base of GHSZ. This helps establish plausible initial conditions for the main part of the study, namely, hydrate recycling due to rapid and continuous sedimentation. For the case study, we apply our numerical model to the Green Canyon Site 955 in the Gulf of Mexico, where the reported high hydrate saturations are likely a result of hydrate recycling driven by rapid sedimentation. In the model, an initial hydrate layer forms due to the invasion of a specified volume of rising free gas. This hydrate layer is consistent with the local pressure, temperature and salinity state. This hydrate layer is then thermally de-stabilised by sedimentation resulting in free gas formation and hydrate recycling. A key finding of our study is that gas hydrate recycling is a cyclic process which can be divided into three phases of 1) gas hydrate melting and free gas nozzling through the hydrate layer, 2) formation of a new gas hydrate layer as the old layer vanishes, and 3) fast uninhibited grow of a new hydrate layer. High hydrate saturations of about 80% can be attained purely through physical, burial-driven recycling of gas hydrates, without any additional gas input from other sources. Hydrate recycling is, therefore, highly dynamic with its own inherent cyclicity rather than a gradual process paced by the rate of sediment deposition.

Description: Highlights • A new numerical model for permafrost in alpine regions. • Importance of lateral fluxes in mountain permafrost modeling. • Influence of unsaturated conditions on freezing processes. • Development of mountain permafrost during warming scenario. Abstract Alpine permafrost environments are highly vulnerable and sensitive to changes in regional and global climate trends. Thawing and degradation of permafrost has numerous adverse environmental, economic, and societal impacts. Mathematical modeling and numerical simulations provide powerful tools for predicting the degree of degradation and evolution of subsurface permafrost as a result of global warming. A particularly significant characteristic of alpine environments is the high variability in their surface geometry which drives large lateral thermal and fluid fluxes along topographic gradients. The combination of these topography-driven fluxes and unsaturated ground makes alpine systems markedly different from Arctic permafrost environments and general geotechnical ground freezing applications, and therefore, alpine permafrost demands its own specialized modeling approaches. In this work, we present a multi-physics permafrost model tailored to subsurface processes of alpine regions. In particular, we resolve the ice–water phase transitions, unsaturated conditions, and capillary actions, and account for the impact of the evolving pore space through freezing and thawing processes. Moreover, the approach is multi-dimensional, and therefore, inherently resolves the topography-driven horizontal fluxes. Through numerical case studies based on the elevation profiles of the Zugspitze (DE) and the Matterhorn (CH), we show the strong influence of lateral fluxes in 2D on active layer dynamics and the distribution of permafrost.

Description: Highlights • Well-defined periodic states are embedded within the steady-state hydrate dynamics. • Periodic states lead to cyclic formation and dissociation of massive hydrate layers. • Periodic states are fully self-sustaining even in the absence of external triggers. • Spontaneous gas migration & pressure release occur in supposedly unperturbed systems. • Existence of periodic states implies an irreducible uncertainty in hydrate dynamics. Abstract Gas hydrates are one of the largest marine carbon reservoirs on Earth. The conventional understanding of hydrate dynamics assumes that the system, in the absence of external triggers, converges to a steady-state over geological time-scales, achieving fixed concentrations of gas hydrate and free gas phase. However, using a high-fidelity numerical model and consistently resolving phase states across multiple fluid-fluid and fluid-solid phase boundaries, we have identified well-defined periodic states embedded within hydrate system dynamics. These states lead to cyclic formation and dissolution of massive hydrate layers that is self-sustaining for the majority of natural marine settings. This previously unresolved characteristic could manifest as spontaneous gas migration and pressure release in, supposedly, unperturbed systems. Our findings show that the gas hydrate systems are not bound to have unique steady-state solutions. Instead, existence of periodic states introduces an irreducible, but, quantifiable uncertainty in gas hydrate dynamics which adds significant error bars to global gas hydrate inventory estimates.

Description: Highlights • This study simulates the sedimentation-driven development of multiple stacked BSRs in the Danube paleo-delta, Black Sea. • Formation of multiple BSRs in the Black Sea is controlled by the sequence of sedimentation events of the levees induced by sea-level changes. • Kinetics of phase transitions plays a key role in the coexistence, location, and timing of the multiple BSRs. • Development of multiple stacked BSRs is possible only under a narrow range of parameters, unique for the Danube delta setting. Abstract The gas hydrate stability zone (GHSZ) is defined by pressure-temperature-salinity (pTS) constraints of natural gas hydrate (GH) system. It refers to a depth interval which usually extends several hundred meters into the sediment column at sufficient water depths. The lower boundary of the GHSZ often coincides in seismic reflection data with a bottom simulating reflector (BSR), which indicates the transition between the underlying free gas and the overlying no-free gas zone at the thermodynamic stability boundary. The GHSZ in geological systems is dynamic and can shift in response to sedimentation processes and/or changes in environmental conditions such as bottom water temperatures, hydrostatic pressure, and water salinity. The appearance of multiple BSRs has been interpreted as remnants of former GHSZ shifts which have persisted over geological timescales. In this study, we numerically simulate the sedimentation-driven development of multiple stacked BSRs in the Danube deep-sea fan in the Black Sea. We show that in this dynamic sediment depositional regime sufficient amounts of residual gas remain trapped in the former GHSZ, given sufficiently high initial gas hydrate saturations, so that paleo-BSRs could persist over long time scales (similar to 300 kyr). In particular, the formation and persistence of multiple BSRs in the Danube Delta is controlled by the sequence of sedimentation events of the levees induced by sea-level change. The kinetics of methane phase transitions between gas hydrate, dissolved methane, and free gas plays a key role in the coexistence, location and timing of the multiple BSRs. Thus, For a given permeability, distinct multiple BSRs appear only for a narrow range of GH formation (10(-14) 〈 k(f) [mol/m(2) Pa s] 〈= 10(-12)) and dissociation rates (10(-16) 〈 k(d) [mol/m(2) Pa s] 〈 10(-14)).