Description: 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.

Description: The exchange of CH4 by CO2 in gas hydrates is of interest for the production of natural gas from methane hydrate with net zero climate gas balance, and for managing risks that are related to sediment destabilization and mobilization after gas-hydrate dissociation. Several experimental studies on the dynamics and efficiency of the process exist, but the results seem to be partly inconsistent. We used confocal Raman spectroscopy to map an area of several tens to hundreds µm of a CH4 hydrate sample during its exposure to liquid and gaseous CO2. On this scale, we could identify and follow different processes in the sample that occur in parallel. Next to guest-molecule exchange, gas-hydrate dissociation also contributes to the release of CH4. During our examination period, about 50% of the CO2 was bound by exchange for CH4 molecules, while the other half was bound by new formation of CO2 hydrates. We evaluated single gas-hydrate grains with confirmed gas exchange and applied a diffusion equation to quantify the process. Obtained diffusion coefficients are in the range of 10−13–10−18 m2/s. We propose to use this analytical diffusion equation for a simple and robust modeling of CH4 production by guest-molecule exchange and to combine it with an additional term for gas-hydrate dissociation.

Description: Seafloor heat flow measurements are utilized to determine the geothermal regime of the Danube deep-sea fan in the western Black Sea and are presented in the larger context of regional gas hydrate occurrences. Heat flow data were collected across paleo-channels in water depths of 550–1460 m. Heat flow across levees ranges from 25 to 30 mW m−2 but is up to 65 mW m−2 on channel floors. Gravity coring reveals sediment layers typical of the western Black Sea, consisting of three late Pleistocene to Holocene units, notably red clay within the lowermost unit cored. Heat flow derived from the bottom-simulating reflector (BSR), assumed to represent the base of the gas hydrate stability zone (GHSZ), deviates from seafloor measurements. These discrepancies are linked either to fast sedimentation or slumping and associated variations in sediment physical properties. Topographic effects account of up to 50% of heat flow deviations from average values. Combined with climate-induced variations in seafloor temperature and sea-level since the last glacial maximum large uncertainties in the prediction of the base of the GHSZ remain. A regional representative heat flow value is ~30 mW m−2 for the study region but deviations from this value may be up to 100%.

Description: The "guest exchange"of methane (CH4) by carbon dioxide (CO2) in naturally occurring gas hydrates is seen as a possibility to concurrently produce CH4 and sequester CO2. Presently, process evaluation is based on CH4-CO2 exchange yields of small-or medium-scale laboratory experiments, mostly neglecting mass and heat transfer processes. This work investigates process efficiencies in two large-scale experiments (210 L sample volume) using fully water-saturated, natural reservoir conditions and a gas hydrate saturation of 50%. After injecting 50 kg of heated CO2 discontinuously (E1) and continuously (E2) and a subsequent soaking period, the reservoir was depressurized discontinuously. It was monitored using electrical resistivity, temperature and pressure sensors, and fluid flow and gas composition measurements. Phase and component inventories were analyzed based on mass and volume balances. The total CH4 production during CO2 injection was only 5% of the initial CH4 inventory. Prior to CO2 breakthrough, the produced CH4 roughly equaled dissolved CH4 in the produced pore water, which balanced the volume of the injected CO2. After CO2 breakthrough, CH4 ratios in the released CO2 quickly dropped to 2.0-0.5 vol %. The total CO2 retention was the highest just before the CO2 breakthrough and higher in E1 where discontinuous injection improved the distribution of injected CO2 and subsequent mixed hydrate formation. The processes were improved by the succession of CO2 injection by controlled degassing at stability limits below that of the pure CH4 hydrate, particularly in experiment E2. Here, a more heterogeneous distribution of liquid CO2 and larger availability of free water led to smaller initial degassing of liquid CO2. This allowed for quick re-formation of mixed gas hydrates and CH4 ratios of 50% in the produced gases. The experiments demonstrate the importance of fluid migration patterns, heat transport, sample inhomogeneity, and secondary gas hydrate formation in water-saturated sediments.