Summary: 3-D X-ray micro-CT (XCT) is a non-destructive 3-D imaging method, increasingly used for a wide range of applications in Earth Science. An optimal XCT image-processing workflow is derived here for accurate quantification of porosity and absolute permeability of heterogeneous sandstone samples using an assessment of key image acquisition and processing parameters: image resolution, segmentation method, representative elementary volume (REV) size and fluid-simulation method. XCT image-based calculations obtained for heterogeneous sandstones are compared to two homogeneous standards (Berea sandstone and a sphere pack), as well as to the results from physical laboratory measurements. An optimal XCT methodology obtains porosity and permeability results within ±2 per cent and vary by one order of magnitude around the direct physical measurements, respectively, achieved by incorporating the clay fraction and cement matrix (porous, impermeable components) to the pore-phase for porosity calculations and into the solid-phase for permeability calculations. Two stokes-flow finite element modelling (FEM) simulation methods, using a voxelized grid (Avizo) and tetrahedral mesh (Comsol) produce comparable results, and similarly show that a lower resolution scan (∼5 µm) is unable to resolve the smallest intergranular pores, causing an underestimation of porosity by ∼3.5 per cent. Downsampling the image-resolution post-segmentation (numerical coarsening) and pore network modelling both allow achieving of a REV size, whilst significantly reducing fluid simulation memory requirements. For the heterogeneous sandstones, REV size for permeability (≥1 mm3) is larger than for porosity (≥0.5 mm3) due to tortuosity of the fluid paths. This highlights that porosity should not be used as a reference REV for permeability calculations. The findings suggest that distinct image processing workflows for porosity and permeability would significantly enhance the accurate quantification of the two properties from XCT.

A better understanding of the effect of methane hydrate morphology and saturation on elastic wave velocity of hydrate‐bearing sediments is needed for improved seafloor hydrate resource and geohazard assessment. We conducted X‐ray synchrotron time‐lapse 4‐D imaging of methane hydrate evolution in Leighton Buzzard sand and compared the results to analogous hydrate formation and dissociation experiments in Berea sandstone, on which we measured ultrasonic P and S wave velocities and electrical resistivity. The imaging experiment showed that initially hydrate envelops gas bubbles and methane escapes from these bubbles via rupture of hydrate shells, leading to smaller bubbles. This process leads to a transition from pore‐floating to pore‐bridging hydrate morphology. Finally, pore‐bridging hydrate coalesces with that from adjacent pores creating an interpore hydrate framework that interlocks the sand grains. We also observed isolated pockets of gas within hydrate. We observed distinct changes in gradient of P and S wave velocities increase with hydrate saturation. Informed by a theoretical model of idealized hydrate morphology and its influence on elastic wave velocity, we were able to link velocity changes to hydrate morphology progression from initial pore‐floating, then pore‐bridging, to an interpore hydrate framework. The latter observation is the first evidence of this type of hydrate morphology and its measurable effect on velocity. We found anomalously low S wave velocity compared to the effective medium model, probably caused by the presence of a water film between hydrate and mineral grains.

Methane is a potent greenhouse gas and large‐scale rapid release of methane from hydrate may have contributed to past abrupt climate change inferred from the geological record. The discovery in 2008 of over 250 plumes of methane gas escaping from the seabed of the West Svalbard continental margin at ~400 m water depth (mwd) suggests that hydrate is dissociating in the present‐day Arctic. Here we model the dynamic response of hydrate‐bearing sediments over a period of 2300 years and investigate ocean warming as a possible cause for present‐day and likely future dissociation of hydrate, within 350–800 mwd, west of Svalbard. Future temperatures are given by two climate models, HadGEM2 and CCSM4, and scenarios, Representative Concentration Pathways (RCPs) 8.5 and 2.6. Our results suggest that over the next three centuries 5.3–29 Gg yr⁻¹ of methane may be released to the Arctic Ocean on the West Svalbard margin.

The majority of presently exploitable marine methane hydrate reservoirs are likely to host hydrate in disseminated form in coarse grain sediments. For hydrate concentrations below 25–40%, disseminated or pore‐filling hydrate does not increase elastic frame moduli, thus making impotent traditional seismic velocity‐based methods. Here, we present a theoretical model to calculate frequency‐dependent P and S wave velocity and attenuation of an effective porous medium composed of solid mineral grains, methane hydrate, methane gas, and water. The model considers elastic wave energy losses caused by local viscous flow both (i) between fluid inclusions in hydrate and pores and (ii) between different aspect ratio pores (created when hydrate grows); the inertial motion of the frame with respect to the pore fluid (Biot's type fluid flow); and gas bubble damping. The sole presence of pore‐filling hydrate in the sediment reduces the available porosity and intrinsic permeability of the sediment affecting Biot's type attenuation at high frequencies. Our model shows that attenuation maxima due to fluid inclusions in hydrate are possible over the entire frequency range of interest to exploration seismology (1–10⁶ Hz), depending on the aspect ratio of the inclusions, whereas maxima due to different aspect ratio pores occur only at sonic to ultrasound frequencies (10⁴–10⁶ Hz). This frequency response imposes further constraints on possible hydrate saturations able to reproduce broadband elastic measurements of velocity and attenuation. Our results provide a physical basis for detecting the presence and amount of pore‐filling hydrate in seafloor sediments using conventional seismic surveys.

Methane hydrate saturation estimates from remote geophysical data and borehole logs are needed to assess the role of hydrates in climate change, continental slope stability, and energy resource potential. Here we present laboratory hydrate formation/dissociation experiments in which we determined the methane hydrate content independently from pore pressure and temperature and from electrical resistivity. Using these laboratory experiments, we demonstrate that hydrate formation does not take up all the methane gas or water even if the system is under two phase water‐hydrate stability conditions and gas is well distributed in the sample. The experiment started with methane gas and water saturations of 16.5% and 83.5%, respectively; during the experiment, hydrate saturation proceeded up to 26% along with 12% gas and 62% water remaining in the system. The coexistence of hydrate and gas is one possible explanation for discrepancies between estimates of hydrate saturation from electrical and acoustic methods. We suggest that an important mechanism for this coexistence is the formation of a hydrate film enveloping methane gas bubbles, trapping the remaining gas inside.