Controlled-source electromagnetic and seismic delineation of sub-seafloor fluid flow structures in a gas hydrate province, offshore Norway

Deep sea pockmarks underlain by chimney-like or pipe structures that contain methane hydrate are abundant along the Norwegian continental margin. In such hydrate provinces the interaction between hydrate formation and fluid flow has significance for benthic ecosystems and possibly climate change. The Nyegga region, situated on the western Norwegian continental slope, is characterized by an extensive pockmark field known to accommodate substantial methane gas hydrate deposits. The aim of this study is to detect and delineate both the gas hydrate and free gas reservoirs at one of Nyegga's pockmarks. In 2012, a marine controlled-source electromagnetic (CSEM) survey was performed at a pockmark in this region, where high-resolution three-dimensional seismic data were previously collected in 2006. Two-dimensional CSEM inversions were computed using the data acquired by ocean bottom electrical field receivers. Our results, derived from unconstrained and seismically constrained CSEM inversions, suggest the presence of two distinctive resistivity anomalies beneath the pockmark: a shallow vertical anomaly at the underlying pipe structure, likely due to gas hydrate accumulation, and a laterally extensive anomaly attributed to a free gas zone below the base of the gas hydrate stability zone. This work contributes to a robust characterization of gas hydrate deposits within sub-seafloor fluid flow pipe structures

Comparison of 2-D and 3-D full waveform inversion imaging using wide-angle seismic data from the Deep Galicia Margin

Full waveform inversion (FWI) is a data-fitting technique capable of generating high-resolution velocity models with a resolution down to half the seismic wavelength. FWI is applied typically to densely sampled seismic data. In this study, we applied FWI to 3-D wide-angle seismic data acquired using sparsely spaced ocean bottom seismometers (OBSs) from the Deep Galicia Margin west of Iberia. Our data set samples the S-reflector, a low-angle detachment present in this area. Here we highlight differences between 2-D, 2.5-D and 3-D-FWI performances using a real sparsely spaced data set. We performed 3-D FWI in the time domain and compared the results with 2-D and 2.5-D FWI results from a profile through the 3-D model. When overlaid on multichannel seismic images, the 3-D FWI results constrain better the complex faulting within the pre- and syn-rift sediments and crystalline crust compared to the 2-D result. Furthermore, we estimate variable serpentinization of the upper mantle below the S-reflector along the profile using 3-D FWI, reaching a maximum of 45 per cent. Differences in the data residuals of the 2-D, 2.5-D and 3-D inversions suggest that 2-D inversion can be prone to overfitting when using a sparse data set. To validate our results, we performed tests to recover the anomalies introduced by the inversions in the final models using synthetic data sets. Based on our comparison of the velocity models, we conclude that the use of 3-D data can partially mitigate the problem of receiver sparsity in FWI

Elastic and electrical properties and permeability of serpentinites from Atlantis Massif, Mid-Atlantic Ridge

Serpentinized peridotites co-exist with mafic rocks in a variety of marine environments including subduction zones, continental rifts and mid-ocean ridges. Remote geophysical methods are crucial to distinguish between them and improve the understanding of the tectonic, magmatic and metamorphic history of the oceanic crust. But, serpentinite peridotites exhibit a wide range of physical properties that complicate such a distinction. We analyzed the ultrasonic P- and S-wave velocities (Vp, Vs) and their respective attenuation (Qp−1, Qs−1), electrical resistivity and permeability of four serpentinized peridotite samples from the southern wall of the Atlantis Massif, Mid-Atlantic Ridge, collected during International Ocean Discovery Program (IODP) Expedition 357. The measurements were taken over a range of loading-unloading stress paths (5 - 45 MPa), using ∼1.7 cm length, 5 cm diameter samples horizontally extracted from the original cores drilled on the seafloor. The measured parameters showed variable degrees of stress dependence, but followed similar trends. Vp, Vs, resistivity and permeability show good inter-correlations, while relationships that included Qp−1 and Qs−1 are less clear. Resistivity showed high contrast between highly serpentinized ultramafic matrix (> 50 Ω m) and mechanically/geochemically altered (magmatic/hydrothermal-driven alteration) domains (< 20 Ω m). This information together with the elastic constants (Vp/Vs ratio and bulk moduli) of the samples allowed us to infer useful information about the degree of serpentinization and the alteration state of the rock, contrasted by petrographic analysis. This study shows the potential of combining seismic techniques and controlled source electromagnetic surveys for understanding tectono-magmatic processes and fluid pathways in hydrothermal systems

Laboratory observations of frequency-dependent ultrasonic P-wave velocity and attenuation during methane hydrate formation in Berea sandstone

Knowledge of the effect of methane hydrate saturation and morphology on elastic wave attenuation could help reduce ambiguity in seafloor hydrate content estimates. These are needed for seafloor resource and geohazard assessment, as well as to improve predictions of greenhouse gas fluxes into the water column. At low hydrate saturations, measuring attenuation can be particularly useful as the seismic velocity of hydrate-bearing sediments is relatively insensitive to hydrate content. Here, we present laboratory ultrasonic (448–782 kHz) measurements of P-wave velocity and attenuation for successive cycles of methane hydrate formation (maximum hydrate saturation of 26 per cent) in Berea sandstone. We observed systematic and repeatable changes in the velocity and attenuation frequency spectra with hydrate saturation. Attenuation generally increases with hydrate saturation, and with measurement frequency at hydrate saturations below 6 per cent. For hydrate saturations greater than 6 per cent, attenuation decreases with frequency. The results support earlier experimental observations of frequency-dependent attenuation peaks at specific hydrate saturations. We used an effective medium rock-physics model which considers attenuation from gas bubble resonance, inertial fluid flow and squirt flow from both fluid inclusions in hydrate and different aspect ratio pores created during hydrate formation. Using this model, we linked the measured attenuation spectral changes to a decrease in coexisting methane gas bubble radius, and creation of different aspect ratio pores during hydrate formation

High-resolution resistivity imaging of marine gas hydrate structures by combined inversion of CSEM towed and ocean-bottom receiver data

We present high-resolution resistivity imaging of gas hydrate pipe-like structures, as derived from marine controlled-source electromagnetic (CSEM) inversions that combine towed and ocean-bottom electric field receiver data, acquired from the Nyegga region, offshore Norway. 2.5-D CSEM inversions applied to the towed receiver data detected four new prominent vertical resistive features that are likely gas hydrate structures, located in proximity to a major gas hydrate pipe-like structure, known as the CNE03 pockmark. The resistivity model resulting from the CSEM data inversion resolved the CNE03 hydrate structure in high resolution, as inferred by comparison to seismically constrained inversions. Our results indicate that shallow gas hydrate vertical features can be delineated effectively by inverting both ocean-bottom and towed receiver CSEM data simultaneously. The approach applied here can be utilized to map and monitor seafloor mineralization, freshwater reservoirs, CO2 sequestration sites and near-surface geothermal systems

Laboratory insights into the effect of sediment-hosted methane hydrate morphology on elastic wave velocity from time-lapse 4D synchrotron X-ray computed tomography.

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 4D 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 velocity, 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 inter‐pore 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 velocity 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 inter‐pore 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

Hydrophone component of Ocean bottom recording sgy-files of seismic refraction and wide angle data from profile 1172 and 2222 of expedition MGL1307 and POSEIDON cruise POS453 with links to data files

Hyperextension of continental crust at the Deep Galicia rifted margin in the North Atlantic has been accommodated by the rotation of continental fault blocks, which are underlain by the S reflector, an interpreted detachment fault, along which exhumed and serpentinized mantle peridotite is observed. West of these features, the enigmatic Peridotite Ridge has been inferred to delimit the western extent of the continent‐ocean transition. An outstanding question at this margin is where oceanic crust begins, with little existing data to constrain this boundary and a lack of clear seafloor spreading magnetic anomalies. Here we present results from a 160 km long wide‐angle seismic profile (Western Extension 1). Travel time tomography models of the crustal compressional velocity structure reveal highly thinned and rotated crustal blocks separated from the underlying mantle by the S reflector. The S reflector correlates with the 6.0–7.0 km s−1 velocity contours, corresponding to peridotite serpentinization of 60–30%, respectively. West of the Peridotite Ridge, shallow and sparse Moho reflections indicate the earliest formation of an anomalously thin oceanic crustal layer, which increases in thickness from ~0.5 km at ~20 km west of the Peridotite Ridge to ~1.5 km, 35 km further west. P wave velocities increase smoothly and rapidly below top basement, to a depth of 2.8–3.5 km, with an average velocity gradient of 1.0 s−1. Below this, velocities slowly increase toward typical mantle velocities. Such a downward increase into mantle velocities is interpreted as decreasing serpentinization of mantle rock with depth

The effect of heterogeneities in hydrate saturation on gas production from natural systems

Understanding the rate and time evolution of gas release from natural gas hydrate systems is important when evaluating the potential of gas hydrate as a future energy source, or the impact of gas from hydrate on climate. The release of gas from hydrate is heavily influenced by a number of factors, many of which vary through the hydrate system. The fundamental heterogeneity of natural gas hydrate systems is often poorly represented in models. Here we simulate depressurisation-induced gas production from a single vertical well in 34 models with heterogeneous 2D distributions of hydrate that include layered, columnar or random configurations and comparable models with homogenous saturation distributions. We found that the temporal evolution of gas production rate follows a consistent trend for all models, but at any time the gas production rate across the models varied by up to ±35% in the first year of production, and by up to ±25% thereafter. The primary control on the gas production rate is the overall amount of hydrate in the system, but local variations in hydrate saturation cause significant fluctuations in the time evolution of production. These hydrate variations can cause changes in the gas flow path through the system and associated drops in gas production rate continuing for multiple years. Overall, our results suggest that small levels of heterogeneity in hydrate systems can cause variations in the gas production rate similar in scale to much larger variations in homogenous systems. Our work provides an error margin for previously modelled gas production rates, and a note of caution for potential commercial development of gas hydrate