Bottom Simulating Reflector Research Articles

Characterization of gas accumulation at a venting gas hydrate system in the Shenhu area, South China Sea

Bottom-simulating reflections (BSRs) in seismic data have been widely accepted to indicate the base of the methane gas hydrate stability zone (MGHSZ), and free gas was thought to exist only below it. However, real geologic systems are far more complex. We have evaluated the results of 3D seismic, logging while drilling, in situ, and coring measurements at a venting gas hydrate system in the Shenhu area of the South China Sea. Our studies reveal that free gas has migrated upward through the thermogenic gas hydrate stability zone into the MGHSZ and become a part of the gas hydrate system. Seismic amplitude anomalies and core results suggest the presence of free gas above the base of MGHSZ at 165 mbsf and the presence of thermogenic gas hydrates below it in well SC-W01. Analyses of P-wave velocity, S-wave velocity, density, and porosity logs reveal that free gas occurs above and below the MGHSZ as well. Integrating log and core analysis with seismic interpretation suggests that the variation in seismic amplitude within the chaotic zone is associated with variable gas saturations, and a large amount of methane and thermogenic gases accumulate near the complex BSRs. We suggest that relative permeability likely plays a significant role in the free-gas distribution and the formation of gas hydrates within a venting gas hydrate system, whereas the effect of dissolved gas short migration is not ignored. Our results have important implications for understanding the accumulation and distribution of gas hydrates and free gas in the venting gas hydrate system and seeps at the seafloor.

Mud volcano possibly linked to seismogenic faults in the Kumano Basin, Nankai Trough, Japan

We investigated a mud volcano (MV) in a fault zone located at the southern edge of the Kumano Basin, the largest forearc basin along the Nankai Trough. Existing seismic sections show a truncated bottom-simulating reflection by a conduit below a topographic high, indicating the presence of an MV. New shipboard acoustic observations show that the fluid may be seeping through the seafloor, which in turn indicates that there are sufficient fluids for larger scale fluid migration in this area. Autonomous underwater vehicle-based high resolution acoustic observations and pH measurements indicate that soft sediment covers most of the MV and surrounding seafloor and that mud and small amounts of “high-backscattered materials” are sprinkled within the crater and around the MV. The MV type is different from those in the landward part of the Kumano Basin: the southern MV is smaller in size, has a steeper slope angle than the those in the northern Kumano Basin, and is located in a fault zone. The characteristics of this 14th Kumano Basin MV suggest that it is an expression of the larger scale fluid and sediment migration along the southern edge of the Kumano Basin, which may provide information regarding fluid and sediment migration along fault systems in the Nankai Trough accretionary prism.

The many double BSRs across the northern Hikurangi margin and their implications for subduction processes

The bottom simulating reflection (BSR) is widely observed along continental margins and is believed to mark the base of gas hydrate stability zone (BGHSZ). In some regions, double or multiple overlapping BSRs are observed, yet their formation mechanisms and geologic implications are not well understood. Here we present 3D seismic images from the 2018 NZ3D experiment that covers a 14×60km2 survey area on New Zealand's northern Hikurangi subduction margin. We observe double BSRs in five locations. Beneath the Tuaheni Basin in the mid-slope, a secondary BSR (BSR2) lies ∼100-360 m deeper than the primary BSR (BSR1) and its 3D geometry mimics the unconformity at the base of the basin. At three thrust ridges located 18-38 km from the deformation front, BSR2 lies ∼55-130 m below and is sub-parallel to BSR1. At another thrust ridge ∼14 km from the deformation front, BSR2 forms above the BSR1, and the two BSRs converge towards the peak of the ridge. Through 3D modeling of BGHSZ and analysis of the geometry and reflection characteristics of the double BSRs, we identify three potential mechanisms for their formation (1) rapid sedimentation, (2) tectonic uplift, (3) overpressure/heat advection caused by fluid migration. Our study demonstrates that formation of double BSRs is closely linked to subduction processes along the northern Hikurangi margin, and double BSRs may be used as indicators for areas with recent sedimentation, tectonic and/or fluid activities.

Identification and analysis of bottom simulating reflectors in the Foz do Amazonas Basin, Northern Brazil

This work proposed an additional approach to investigate gas hydrate occurrences in the Foz do Amazonas Basin, in the Brazilian Equatorial Margin. The automatic comparison of seafloor seismic amplitudes with those from the Bottom Simulating Reflectors (BSR) was used to exhibit the reversal of signal polarities among the seafloor (positive) and the BSR (negative), reinforcing the identification of BSR in areas where its visualization was unclear. Additionally, we used the envelope attribute to highlight the BSR in the seismic section. Subsequently, we decomposed the seismic data into different frequency bands, applied a -90 degrees phase rotation to the data and recalculated the envelope attribute for each section decomposed in frequency bands. This technique improved visualization, allowing the identification of intervals where BSR were laterally discontinuous, revealing to be valuable for mapping the gas hydrate distribution in Foz do Amazonas Basin.

A Long-Term Geothermal Observatory Across Subseafloor Gas Hydrates, IODP Hole U1364A, Cascadia Accretionary Prism

We report 4 years of temperature profiles collected from May 2014 to May 2018 in Integrated Ocean Drilling Program Hole U1364A in the frontal accretionary prism of the Cascadia subduction zone. The temperature data extend to depths of nearly 300 m below seafloor (mbsf), spanning the gas hydrate stability zone at the location and a clear bottom-simulating reflector (BSR) at ∼230 mbsf. When the hole was drilled in 2010, a pressure-monitoring Advanced CORK (ACORK) observatory was installed, sealed at the bottom by a bridge plug and cement below 302 mbsf. In May 2014, a temperature profile was collected by lowering a probe down the hole from the ROV ROPOS. From July 2016 through May 2018, temperature data were collected during a nearly two-year deployment of a 24-thermistor cable installed to 268 m below seafloor (mbsf). The cable and a seismic-tilt instrument package also deployed in 2016 were connected to the Ocean Networks Canada (ONC) NEPTUNE cabled observatory in June of 2017, after which the thermistor temperatures were logged by Ocean Networks Canada at one-minute intervals until failure of the main ethernet switch in the integrated seafloor control unit in May 2018. The thermistor array had been designed with concentrated vertical spacing around the bottom-simulating reflector and two pressure-monitoring screens at 203 and 244 mbsf, with wider thermistor spacing elsewhere to document the geothermal state up to seafloor. The 4 years of data show a generally linear temperature gradient of 0.055°C/m consistent with a heat flux of 61–64 mW/m2. The data show no indications of thermal transients. A slight departure from a linear gradient provides an approximate limit of ∼10−10 m/s for any possible slow upward advection of pore fluids. In-situ temperatures are ∼15.8°C at the BSR position, consistent with methane hydrate stability at that depth and pressure.

Distribution and development of submarine mud volcanoes on the Makran Continental Margin, offshore Pakistan

Abstract In 2018, China Geological Survey and National Institute of Oceanography of Pakistan conducted a comprehensive geological survey on the Makran Continental Margin offshore Pakistan using the R/V Haiyangdizhi10. Seismic profiles, multibeam bathymetric and water column data, and seafloor observations were obtained during the expedition. Twelve mud volcanoes and three pockmarks were discovered in the surveyed area on the Makran Continental Margin. The heights of the mud volcanoes and the depths of the pockmarks are 12–121 m and −108 to −26 m, respectively, with diameter of 500–1400 m. The mud volcanoes are distributed on top of the W-E trending ridges at 1200–2800 m water depth. Six mud volcanoes and one pockmark were active during the investigation as indicated by the gas flares mapped in the multibeam water column data. Acoustic blanking zones, faults and gas chimneys are distributed ubiquitously beneath the mud volcanoes on the seismic profiles. Their spatial correlation with mud volcanoes indicates that the faults and gas chimneys provide pathways for fluid migration. In addition, the Bottom Simulating Reflectors (BSR), interface of gas hydrate and free gas, were also discovered in multiple seismic profiles. Based on the BSR discontinuity and piercing of the gas chimneys, it is inferred that mud volcano activities might cause gas hydrate dissociation, leading to the upward movement of the BSR. Tectonic uplifting, seismic-induced shock and fluid overpressure induced by rapid sedimentation are the major driving forces for the mud volcano formation and development. Gases from deep strata and gas hydrate decomposition are suspected to provide additional forces for fluid migration. Pressure decrease in the deep subsurface and gas hydrate sealing of the faults and fracture are speculated to be important in controlling the intermittent mud volcano activities and potential pockmark formation. The new model of mud volcano activities depicted in this study might be applicable to other mud volcanoes elsewhere.

Delineation of discontinuity using multi-channel seismic attributes: An implication for identifying fractures in gas hydrate sediments in offshore Mahanadi basin

Seismic attributes derived from post stack seismic data are important tools to delineate and characterize the subsurface geological features and discontinuities such as fracture, fault and many more. In this paper we have used two seismic profiles (MH-38A and MH-38B) of 2D multi-channel data along north-south direction with seafloor depth ranging from 1450 m to 2175 m in the deep offshore Mahanadi basin. The gas hydrate indicator, bottom simulating reflector (BSR) is easily identified but the discontinuities such as fractures are not observed directly from the seismic data. Therefore, we have computed the total four most suitable attributes which are reflection strength, instantaneous frequency, instantaneous bandwidth and instantaneous quality factor. We have compared the relative changes of the attribute values along the layers in the four derived attributes and hence, have recognized total three and seven discontinuities for MH-38A and MH-38B respectively, at near and away from the well location. Most faults are passing in NNW-SSE and N-S direction either in the gas hydrate layer or below BSR. These discontinuities are called as natural fractures which are validated by the fractures observed from resistivity image log at well site NGHP-01–19 and NGHP-01–09 for CDP ranges 1561 to 1681 in line MH-38A and 1357 to 1458 in line MH-38B respectively. The discontinuities in the gas hydrate sediments are useful for anisotropy study and to map the chimney system through which the free gas can flow from below BSR to up direction.

Shallow gas and gas hydrate occurrences on the northwest Greenland shelf margin

An extensive 3D seismic dataset was used to investigate the contemporary hydrocarbon distribution and historical fluid migration in Melville Bay offshore northwest Greenland, providing the first inventory of shallow gas and gas hydrate along this part of the Greenland margin. The shallow gas anomalies vary in seismic character and have been subdivided into four categories that represent (I) isolated shallow gas, (II) free gas trapped at the base of the gas hydrate stability zone (GHSZ), (III) gas charged glacial clinoforms and (IV) a giant mass transport deposit gas reservoir. Gas hydrate deposits have been identified across an area of 537 km2 via the identification of a discontinuous bottom simulating reflector (BSR) that marks the base of the GHSZ. The BSR has been used to estimate a geothermal gradient of 49 °C/km across the GHSZ and a heat flow of 70–90 mW/m2, providing the first publically available heat flow estimates offshore western Greenland. The contemporary hydrocarbon distribution and historical fluid migration is influenced by the underlying paleo-rift topography and multiple shelf edge glaciations since ~2.7 Ma. Continued uplift of the Melville Bay Ridge, as well as glacial-sediment redistribution and basinward margin tilting from isostatic compensation, have led to a concentration of gas within the Cenozoic stratigraphy above the ridge. Furthermore, repeated variations in subsurface conditions during glacial-interglacial cycles likely promoted fluid remigration, and possibly contributed to reservoir leakage and increased fluid migration through faults. The top of the gas hydrate occurrence at 650 m water depth is well below the hydrate-free gas phase boundary (~350 m) for the present bottom-water temperature of 1.5 °C, suggesting this hydrate province mainly adjusted to glacial-interglacial changes by expansion and dissociation at its base and is relatively inert to current levels of global warming. Glacial-related dissociation may have significantly contributed to the numerous free gas accumulations observed below the GHSZ at present day.

Multi-attribute machine learning analysis for weak BSR detection in the Pegasus Basin, Offshore New Zealand

Gas hydrates that exist in the subsurface are often difficult to detect with reflection seismic data if the seismic data lacks a strong bottom simulating reflection (BSR). In these cases, the imaging and detection of the gas hydrate stability zone (GHSZ) becomes particularly difficult, as hydrate detection relies heavily on the BSR, gas chimneys, or pockmarks on the seafloor. To address and improve upon these imaging complications, an unsupervised machine learning multi-attribute analysis is performed on 2D seismic data in the Pegasus Basin in New Zealand where the BSR is not continuously or clearly imaged. Rock physics analysis has demonstrated that the inclusion of methane gas hydrates in the pore space results in a slightly increasing amplitude at the base of the gas hydrate zone, regardless of the fluid (brine or gas) in the pore space below the hydrates. This increasing amplitude is quite weak and can be masked by noise. In the scenarios where a strong seismic impedance difference is lacking, a BSR is not typically observed in the seismic data, even though gas hydrates do exist in the subsurface. To enhance the detection of the presence of gas hydrates, a multi-attribute analysis is performed with a series of seismic attributes that can detect the minute changes in the seismic waveform due to the presence of gas hydrates. The successful attributes are those that are sensitive to attenuation, frequency, and small amplitude anomalies.

A Shallow Seabed Dynamic Gas Hydrate System off SW Taiwan: Results From 3‐D Seismic, Thermal, and Fluid Migration Analyses

AbstractLarge amounts of methane, a potent greenhouse gas, are stored in hydrates beneath the seafloor. Sea level changes can trigger massive methane release into the ocean. It is not clear, however, whether surficial seafloor processes can cause comparable discharge. Previously, fluid migration was difficult to study due to a lack of spatially dense seismic and thermal observations. Here we examine a gas hydrate site at Four‐Way‐Closure Ridge off SW Taiwan using a high‐resolution 3‐D seismic cube, together with bottom‐simulating reflections (BSRs) mapped in the cube, a thermal probe data set, and 3‐D thermal modeling results. We document, on a scale of tens of meters, the interaction between surficial sedimentary processes, fluid flow, and a dynamic gas hydrate system. Fluid migrates upward through dipping permeable strata in the limb, the slope basin, and along thrust faults and ridge‐top normal faults. The seismic data also reveal several double BSRs that underlie seabed sedimentary sliding and depositional features. Abrupt changes in subsurface pressure and temperature due to the rapid seabed sedimentary processes can cause a rapid shift of the base of the gas hydrate stability zone. This shift may be either downward or upward and would result in the accumulation or dissociation of hydrate in sediments sandwiched by the double BSRs, respectively. We propose that dynamic surficial processes on the seafloor together with shallow focused fluid flow affect hydrate distribution and saturation at depth and may even result in methane expulsion into the ocean if such localized features are common along convergent plate boundaries.

Assessment of gas hydrate reservoir from inverted seismic impedance and porosity in the northern Hikurangi margin, New Zealand

Hikurangi margin, New Zealand, the shallowest subduction zone on Earth is the unique place to study the slow slip creep-like deformation and seafloor failure, and their probable link to the presence of gas hydrate, cold seeps and fluid migration. International Ocean Discovery Program expedition 372 discovered gas hydrate within thin intervals (ranging from ~5 to 25 m) of the deformed sediments in the Tuaheni landslide complex area of the northern Hikurangi margin. Seismic profiles show typical BSRs (bottom simulating reflectors) cutting across the dipping strata throughout the sections. To characterize the reservoir, here, we invert both acoustic impedance (product of velocity and density) and porosity along two perpendicular seismic profiles crossing the well using a model-based post-stack seismic inversion technique. Results show the layer structure of alternate high-impedance with low-porosity and low-impedance with high-porosity sediments. Interestingly, each layer is consisting of thin and relatively low impedance intra-layers, which clearly indicates the possible pathways for fluid and gas migration, and one of the reasons for seafloor collapse in this region. Very low impedance observed below the BSR is due to the gas-charged sediments at the dipping strata, and is the probable source of free-gas migrating upward. Next, we apply rock physics theory to know the distribution of gas hydrate at Site U1517 as well as along two perpendicular seismic profiles. Rock physics modeling using sonic velocity at well location shows that gas hydrate is distributed mainly within the depth intervals of 5–50 m, 90–115 m and 130–155 m with an average saturation of about 10–15% and the maximum concentration of about 35% of the pore space at 100 m depth. Prediction of gas hydrate saturation from the resistivity method using neutron porosity and NMR (nuclear magnetic resonance) logs in Archie's empirical formula shows less saturation compared to that from the sonic log. Whereas, estimated gas hydrate saturation from the chloride concentrations ranges from 0 to 45% of the pore space. The spatial distributions of gas hydrate along seismic profiles are predicted by rock physics theory using the inverted impedance and porosity. Saturation of gas hydrate along the seismic profiles varies from 0 to 40% with an average of ~10% of the pore space. The correlation between inverted impedance and porosity indicates that the high impedance layers are due to low porosity consolidated sediment without significant gas hydrate concentration. Gas hydrate is distributed in relatively high porosity and low impedance layers along the seismic profiles. Our result shows that the amount of gas hydrate concentration is low in this region and sediments are highly deformed/disturbed in the presence of many fluid/gas migration pathways.

A complex gas hydrate system imaged jointly using MCS and OBS data from the northern South China Sea

Much of our knowledge of gas hydrate in the Dongsha area of the South China Sea has been gained through intensive geophysical surveys and drilling. However, many factors remain unclear, such as the co-existence of shallow gas hydrate and deeper gas hydrate. This lack of clarity is partially due to the lack of a depth-domain velocity model and accurate imagery of the gas hydrate-bearing sediments. In this study, pre-stack depth migration (PSDM) is used to produce subsurface images using a depth-domain velocity model from OBS tomography as the migration velocity. A simple initial velocity model was built using the seafloor depth information derived from multi-channel seismic (MCS) data. This simple model was used to build a more accurate velocity model by using first arrival time tomography from OBS data. Three different approaches were used to assess the reliability of the velocity model. Resolution tests and uncertainties analysis were used to further evaluate the model. The final PSDM image revealed some features more clearly than they were seen on the time migration image. These features help better explain the complex gas hydrate system. The depth of a bottom-simulating reflector (BSR) was picked directly from the PSDM image. This BSR is very consistent with the drilling result. Our work highlights the importance of PSDM in the study of gas hydrate in complex settings and highlights the potential and value of using OBS refraction in building shallow velocity models.

Distribution and genesis of submarine landslides in the northeastern South China Sea

Submarine landslide is a widespread marine geohazard and an important part of the global ‘source to sink’ system. In this article, multi‐channel seismic and topographic data are utilized to investigate the distribution, calculate the geomorphologic parameters, and discuss the genesis of the submarine landslides in the northeastern South China Sea. In the study area, which can be divided into two tectonic provinces: passive continental margin province and accretionary wedge province, 101 landslides are identified, with a total size of 2,358.22 km2. Statistics show that most landslides are of small sizes (<15 km2) and developed along small slope gradient regions. There is a high correlation between the slope gradient and the ratio of headscarp height to runout distance. The size and runout distance of landslide also correlate well. However, there is a low correlation coefficient including slope gradient versus runout distance and slope gradient versus size. The distribution of submarine landslides is consistent with that of bottom simulating reflections (BSRs) and seabed activities, which may imply a possible inner correlation among them. The calculated base of gas hydrate stability zone (BGHSZ) pinches out the seafloor at about 555–556 m below the sea level, where most of the landslide headscarps are located. And the BGHSZ depths are often shallower than the BSRs bottom depths in the upper slope area. Combined with previous studies, we consider that the influence factors of submarine landslides in the study area are earthquakes, canyons, seabed fluid activities, and gas hydrate dissociation. Gas hydrate dissociation is a global influence factor in the study area. In the passive continental margin province, the submarine landslides are mainly canyon‐associated landslides. While in the accretionary wedge province, they are mainly controlled by seabed fluid activities and earthquakes.

Analysis of marine controlled source electromagnetic data for the assessment of gas hydrates in the Danube deep-sea fan, Black Sea

Marine controlled source electromagnetic (CSEM) data have been analyzed as part of a larger interdisciplinary field study to reveal the distribution and concentration of gas hydrates and free gas in two working areas (WAs) in the offshore Danube fan in the western Black Sea. The areas are located in the Bulgarian sector in about 1500 m water depth (WA1) and in the Romanian sector in about 650 m water depth (WA2). Both areas are characterized by channel levee systems and wide spread occurrences of multiple bottom simulating reflections (BSRs) suggesting the presence of gas hydrates. Electrical resistivity models have been derived from two-dimensional (2D) inversions of inline CSEM data using a seafloor-towed electric dipole-dipole system. Comparing the resistivity models with coincident reflection seismic profiles reveals insight in the sediment stratigraphy of the gas hydrate stability zone (GHSZ). Gas hydrate and free gas saturation estimates have been derived with a stochastic approach of Archie's relationship considering uncertainties in the input parameters available from drilling with the MeBo-200 seafloor rig in WA2.The resistivity models generally reflect the transition of marine to lacustrine conditions expressed by a sharp decay of pore water salinities in the top 30–40 m below seafloor caused by freshwater phases of the Black Sea due to sea level low stands in the past. In WA1, we derived saturation estimates of 10–20% within a 100 m thick layer at around 50 m depth below the channel which compares well with estimates from seismic P-wave velocities. The layer extends below the western levee with even higher saturations of 20–30%, but high gas hydrate saturations are unlikely within the fine grained, clayey sediment section, and the high resistivities may reflect different lithologies of lower permeability and porosity. The resistive layer terminates below the eastern levee where increasing resistivities at depth towards a stack of multiple BSRs indicate gas hydrate and free gas concentrations in the order of 10% to locally 30%. WA2 is characterized by a major slope failure at the landward edge of the gas hydrate stability field next to the channel. Gas hydrate saturation estimates within the slump area are close to zero within the GHSZ which is in agreement with coring results of the nearby MeBo drill sites. Elevated resistivities below the steeply upward bending BSR lead to saturation estimates less than 10% of free gas that may have accumulated.

Automated design of a convolutional neural network with multi-scale filters for cost-efficient seismic data classification

Geoscientists mainly identify subsurface geologic features using exploration-derived seismic data. Classification or segmentation of 2D/3D seismic images commonly relies on conventional deep learning methods for image recognition. However, complex reflections of seismic waves tend to form high-dimensional and multi-scale signals, making traditional convolutional neural networks (CNNs) computationally costly. Here we propose a highly efficient and resource-saving CNN architecture (SeismicPatchNet) with topological modules and multi-scale-feature fusion units for classifying seismic data, which was discovered by an automated data-driven search strategy. The storage volume of the architecture parameters (0.73 M) is only ~2.7 MB, ~0.5% of the well-known VGG-16 architecture. SeismicPatchNet predicts nearly 18 times faster than ResNet-50 and shows an overwhelming advantage in identifying Bottom Simulating Reflection (BSR), an indicator of marine gas-hydrate resources. Saliency mapping demonstrated that our architecture captured key features well. These results suggest the prospect of end-to-end interpretation of multiple seismic datasets at extremely low computational cost.

Crustal processes sustain Arctic abiotic gas hydrate and fluid flow systems

The Svyatogor Ridge and surroundings, located on the sediment-covered western flank of the Northern Knipovich Ridge, host extensive gas hydrate and related fluid flow systems. The fluid flow system here manifests in the upper sedimentary sequence as gas hydrates and free gas, indicated by bottom simulating reflections (BSRs) and amplitude anomalies. Using 2D seismic lines and bathymetric data, we map tectonic features such as faults, crustal highs, and indicators of fluid flow processes. Results indicate a strong correlation between crustal faults, crustal highs and fluid accumulations in the overlying sediments, as well as an increase in geothermal gradient over crustal faults. We conclude here that gas generated during the serpentinization of exhumed mantle rocks drive the extensive occurrence of gas hydrate and fluid flow systems in the region and transform faults act as an additional major pathway for fluid circulation.

Downward shift of gas hydrate stability zone due to seafloor erosion in the eastern Dongsha Island, South China Sea

New acquired and reprocessed three-dimensional (3D) seismic data were used to delineate the distribution and characterization of bottom simulating reflections (BSRs) in the Chaoshan Sag, in the eastern part of Dongsha Island, South China Sea. Three submarine canyons with different scales were interpreted from the 3D seismic data, displaying three stages of canyon development and are related with the variation of BSR. Abundant faults were identified from the coherence and ant-tracing attributions extracted from 3D seismic data, which provide the evidence for fluid migration from deeper sediments to the gas hydrate stability zone (GHSZ). The uplift of Dongsha Island created a large number of faults and leads to the increased seafloor erosion. The erosion caused the cooling of the seafloor sediments and deepening of the base of the gas hydrate stability zone, which is attributed to the presence of paleo-BSR and BSR downward shift in the study area. Hence, methane gas may be released during the BSR resetting and gas hydrate dissociation related with seafloor erosion.

Dynamic accumulation of gas hydrates associated with the channel-levee system in the Shenhu area, northern South China Sea

Abstract Research on the formation and distribution of submarine channel systems and associated gas-bearing fluids is of great significance for gas hydrate exploration. Disseminated gas hydrates with high saturation up to 65% were recovered from a submarine ridge, equivalent to the levee of the channel–levee system in the Shenhu area, northern South China Sea. Sedimentary deposits in the submarine ridge were dominated by fine-grained silt and clay-rich silt; gas hydrates with relatively high saturation preferentially accumulated in coarser sediments with less clay content. Although abundant foraminifera fossils may have increased reservoir pore space, their presence was not a necessary condition for high-saturation hydrates. Higher levels of pyrite appeared in the reservoirs corresponding to high-saturation hydrates, which suggests that the reducing environment caused by sufficient methane provided adequate gas to form higher-saturation hydrates. Because of the migration of the channel–levee system, different channels formed their respective depositional systems composed of channel-filling, buried channel-filling, erosion grooves, and slumped turbidities. Relatively coarse-grained deposits were identified in the channel fillings and levees, and the accumulation of hydrates was affected by the lithological features of the sediments and their spatial coupling with the gas hydrate stability zone (GHSZ). GHSZ modeling based on in situ measurements indicated that erosion and sedimentation, as well as variations of the geothermal gradient, resulted in the upward/downward migration of bottom simulating reflectors (BSRs). On the erosion flank of the channel, the strata thinned, and rapid erosion was likely to destroy the shallower BSR, causing gas hydrate decomposition and methane release, and may have caused turbidite slumping and seepage, whereas the strata thickened on the deposition flank of the channel. The BSR in the channel–levee system would gradually move toward the new GHSZ, eventually forming a new BSR; parts of the BSR that formed under the original P–T conditions have remained, and double BSRs occurred in the seismic profile. The thermal fluid that moved upward through a gas chimney may also have caused the migration of the GHSZ, resulting in the emergence of double BSRs. During the lateral migration of the channel and the vertical migration of the gas-bearing fluid, there was a dynamic adjustment relationship between the GHSZ and the erosion–deposition process of the channel, resulting in the dynamic accumulation of hydrates in the Shenhu area. A model to demonstrate the relationship between channel migration and variation of the BSR was established, which is of great significance for understanding the formation and accumulation mechanisms of gas hydrates.

A rock magnetic perspective of gas hydrate occurrences in a high-energy depositional system in the Krishna-Godavari basin, Bay of Bengal

We conducted a detailed rock magnetic study complemented by sedimentological and mineralogical methods on a 177.2-m-long sediment core of Hole NGHP-01-14A to constrain the influence of high-energy depositional environment on the magnetic mineral diagenesis and formation of gas hydrates in the Krishna-Godavari (K-G) basin, Bay of Bengal. Five sediment magnetic zones were identified based on the downcore variation in rock magnetic parameters. The magnetic mineralogies comprised of detrital titanomagnetite, diagenetically formed magnetic iron sulfides and their mixture. A distinct magnetically enhanced zone (Z-3) close to a bottom simulating reflector (BSR) is dominated by higher magnetite content, sediment bulk density, and grain size and probably represents several sand-rich sediment layers formed as a result of short-duration intense sedimentation events and got rapidly buried and provided conducive environment for the formation of gas hydrate deposits at the studied site. The concurrence of gas hydrates and high gas saturations within the sand-rich layers (Z-3) suggests that the high-energy depositional environment had a larger control over hydrate formation at Hole NGHP-01-14A. Two magnetic iron sulfide bearing sediment intervals in Z-1b and Z-1c below sulfate-methane transition zone (SMTZ) suggest that their formation is controlled by microbially mediated diagenetic reactions fuelled by presence of gas hydrates and probably represents the fossil gas hydrate zones. We propose that the rock magnetic variations in the studied sediment core at Hole NGHP-01-14A is controlled by both differential loading of detrital magnetic minerals as well as hydrate-induced late diagenetic magnetic iron sulfide formation.

Distributions of gas hydrate and free gas accumulations associated with upward fluid flow in the Sanriku-Oki forearc basin, northeast Japan

Abstract We applied automated velocity analysis to a 3D seismic data volume to obtain a high-resolution seismic velocity model that we used to investigate the influence of fluid behavior on the subsurface distribution of gas in the Sanriku-Oki forearc basin, northeast Japan. We identified free gas accumulations as zones of low P-wave velocity separated from overlying high-velocity gas hydrates by the clear seismic boundary of the bottom-simulating reflection (BSR) between them. We then used conductive modeling to map upward heat flow in our study area from the depth of the BSR derived from the velocity model. The estimated heat flow demonstrates that upward fluid flux has considerable influences on the distributions of both gas hydrate and free gas. The areas of high heat flow (representing high fluid flux) correspond to underlying permeable geological features such as gas chimneys, faults, and the edges of porous slumps, suggesting that these features provide fluid migration pathways from Eocene to Oligocene source rocks to free gas and hydrate accumulations in overlying sedimentary rocks. The techniques we have employed can contribute to the high-resolution mapping of gas hydrate and free gas accumulations in similar deep-water reservoirs off the eastern continental margin of Japan and at tectonically similar plate subduction margins in other regions.