High-resolution seismic velocity analysis as a tool for exploring gas hydrate systems: An example from New Zealand’s southern Hikurangi margin

2016 ◽  
Vol 4 (1) ◽  
pp. SA1-SA12 ◽  
Author(s):  
Gareth J. Crutchley ◽  
Guy Maslen ◽  
Ingo A. Pecher ◽  
Joshu J. Mountjoy
Keyword(s):  
High Resolution ◽  
Marine Sediments ◽  
Gas Hydrate ◽  
Seismic Data ◽  
Seismic Velocity ◽  
Coarse Grained ◽  
Velocity Analysis ◽  
Low Velocity ◽  
Free Gas ◽  
Velocity Models

The existence of free gas and gas hydrate in the pore spaces of marine sediments causes changes in acoustic velocities that overprint the background lithological velocities of the sediments themselves. Much previous work has determined that such velocity overprinting, if sufficiently pronounced, can be resolved with conventional velocity analysis from long-offset, multichannel seismic data. We used 2D seismic data from a gas hydrate province at the southern end of New Zealand’s Hikurangi subduction margin to describe a workflow for high-resolution velocity analysis that delivered detailed velocity models of shallow marine sediments and their coincident gas hydrate systems. The results showed examples of pronounced low-velocity zones caused by free gas ponding beneath the hydrate layer, as well as high-velocity zones related to gas hydrate deposits. For the seismic interpreter of a gas hydrate system, the velocity results represent an extra “layer” for interpretation that provides important information about the distribution of free gas and gas hydrate. By combining the velocity information from the seismic transect with geologic samples of the seafloor and an understanding of sedimentary processes, we have determined that high gas hydrate concentrations preferentially form within coarse-grained sediments at the proximal end of the Hikurangi Channel. Finer grained sediments expected elsewhere along the seismic transect might preclude the deposition of similarly high gas hydrate concentrations away from the channel.

Fluid dynamics at the base of hydrate-bearing sediments at the Vestnesa Ridge inferred from 5 years of high-resolution 4D seismic surveying

2020 ◽  
Author(s):  
Malin Waage ◽  
Stefan Bünz ◽  
Kate Waghorn ◽  
Sunny Singhorha ◽  
Pavel Serov
Keyword(s):  
High Resolution ◽  
Gas Hydrate ◽  
Seismic Data ◽  
Time Lapse ◽  
Stability Zone ◽  
Seismic Surveys ◽  
Free Gas ◽  
4D Seismic ◽  
Hydrate Bearing Sediments ◽  
Hydrate Stability

<p>The transition from gas hydrate to gas-bearing sediments at the base of the hydrate stability zone (BHSZ) is commonly identified on seismic data as a bottom-simulating reflection (BSR). At this boundary, phase transitions driven by thermal effects, pressure alternations, and gas and water flux exist. Sedimentation, erosion, subsidence, uplift, variations in bottom water temperature or heat flow cause changes in marine gas hydrate stability leading to expansion or reduction of gas hydrate accumulations and associated free gas accumulations. Pressure build-up in gas accumulations trapped beneath the hydrate layer may eventually lead to fracturing of hydrate-bearing sediments that enables advection of fluids into the hydrate layer and potentially seabed seepage. Depletion of gas along zones of weakness creates hydraulic gradients in the free gas zone where gas is forced to migrate along the lower hydrate boundary towards these weakness zones. However, due to lack of “real time” data, the magnitude and timescales of processes at the gas hydrate – gas contact zone remains largely unknown. Here we show results of high resolution 4D seismic surveys at a prominent Arctic gas hydrate accumulation – Vestnesa ridge - capturing dynamics of the gas hydrate and free gas accumulations over 5 years. The 4D time-lapse seismic method has the potential to identify and monitor fluid movement in the subsurface over certain time intervals. Although conventional 4D seismic has a long history of application to monitor fluid changes in petroleum reservoirs, high-resolution seismic data (20-300 Hz) as a tool for 4D fluid monitoring of natural geological processes has been recently identified.<br><br>Our 4D data set consists of four high-resolution P-Cable 3D seismic surveys acquired between 2012 and 2017 in the eastern segment of Vestnesa Ridge. Vestnesa Ridge has an active fluid and gas hydrate system in a contourite drift setting near the Knipovich Ridge offshore W-Svalbard. Large gas flares, ~800 m tall rise from seafloor pockmarks (~700 m diameter) at the ridge axis. Beneath the pockmarks, gas chimneys pierce the hydrate stability zone, and a strong, widespread BSR occurs at depth of 160-180 m bsf. 4D seismic datasets reveal changes in subsurface fluid distribution near the BHSZ on Vestnesa Ridge. In particular, the amplitude along the BSR reflection appears to change across surveys. Disappearance of bright reflections suggest that gas-rich fluids have escaped the free gas zone and possibly migrated into the hydrate stability zone and contributed to a gas hydrate accumulation, or alternatively, migrated laterally along the BSR. Appearance of bright reflection might also indicate lateral migration, ongoing microbial or thermogenic gas supply or be related to other phase transitions. We document that faults, chimneys and lithology constrain these anomalies imposing yet another control on vertical and lateral gas migration and accumulation. These time-lapse differences suggest that (1) we can resolve fluid changes on a year-year timescale in this natural seepage system using high-resolution P-Cable data and (2) that fluids accumulate at, migrate to and migrate from the BHSZ over the same time scale.</p>

Planning a drilling campaign in a petroleum province using high resolution 3D seismic data – IODP proposal 909

2020 ◽  
Author(s):  
David Cox ◽  
Andrew M. W. Newton ◽  
Paul C. Knutz ◽  
Mads Huuse
Keyword(s):  
High Resolution ◽  
Gas Hydrate ◽  
Seismic Data ◽  
Large Scale ◽  
Greenland Ice Sheet ◽  
Hydrocarbon Exploration ◽  
3D Seismic ◽  
Free Gas ◽  
Shallow Subsurface ◽  
3D Seismic Data

<p>A drilling hazard assessment has been completed for a large area of the NW Greenland-Baffin Bay continental shelf. This assessment was in relation to International Ocean Discovery Program (IODP) proposal 909 that aims to drill several sites across the shelf in an attempt to better understand the evolution and variability of the northern Greenland Ice Sheet. The assessment utilised high quality and extensive 3D seismic data that were acquired during recent hydrocarbon exploration interest in the area – a fact that highlights the risk of drilling in a petroleum province and therefore, the importance of this assessment with regards to safety.</p><p>Scattered seismic anomalies are observed within the Cenozoic sedimentary succession covering the rift basins of the Melville Bay region. These features, potentially representing the presence of free gas or gas-rich fluids, vary in nature from isolated anomalies, fault flags, stacked fluid flow features and canyons; all of which pose a significant drilling risk and were actively avoided during site selection. In areas above the Melville Bay Ridge – a feature that dominates the structure of this area – free gas is also observed trapped beneath extensive gas hydrate deposits, identified via a spectacularly imaged bottom simulating reflector marking the base of the gas hydrate stability zone. The location of the hydrate deposits, and the free gas beneath, are likely controlled by a complicated migration history, due to large scale rift-related faulting and migration along sandy aquifer horizons. In other areas, gas is interpreted to have reached the shallow subsurface due to secondary leakage from a deeper gas reservoir on the ridge crest.</p><p>It is clear that hydrocarbon related hazards within this area are varied and abundant, making it a more challenging location to select sites for an IODP drilling campaign. However, due to the extensive coverage and high resolution (up to 11 m vertical resolution (45 Hz at 2.0 km/s velocity) of the 3D seismic data available, as well as the use of recently acquired ultra-high resolution site survey lines, these features can be accurately imaged and confidently mapped. This allowed for the development of a detailed understanding of the character and distribution of fluids within the shallow subsurface, and the use of this knowledge to select site localities that maximise the potential for drilling to be completed safely and successfully if proposal 909 were to be executed.</p>

scholarly journals Gas hydrate and free gas detection using seismic quality factor estimates from high-resolution P-cable 3D seismic data

2016 ◽  
Vol 4 (1) ◽  
pp. SA39-SA54 ◽  
Author(s):  
Sunny Singhroha ◽  
Stefan Bünz ◽  
Andreia Plaza-Faverola ◽  
Shyam Chand
Keyword(s):  
High Resolution ◽  
Frequency Shift ◽  
Gas Hydrates ◽  
Gas Hydrate ◽  
Seismic Data ◽  
Seismic Attenuation ◽  
3D Seismic ◽  
Free Gas ◽  
Centroid Frequency ◽  
3D Seismic Data

We have estimated the seismic attenuation in gas hydrate and free-gas-bearing sediments from high-resolution P-cable 3D seismic data from the Vestnesa Ridge on the Arctic continental margin of Svalbard. P-cable data have a broad bandwidth (20–300 Hz), which is extremely advantageous in estimating seismic attenuation in a medium. The seismic quality factor (Q), the inverse of seismic attenuation, is estimated from the seismic data set using the centroid frequency shift and spectral ratio (SR) methods. The centroid frequency shift method establishes a relationship between the change in the centroid frequency of an amplitude spectrum and the Q value of a medium. The SR method estimates the Q value of a medium by studying the differential decay of different frequencies. The broad bandwidth and short offset characteristics of the P-cable data set are useful to continuously map the Q for different layers throughout the 3D seismic volume. The centroid frequency shift method is found to be relatively more stable than the SR method. Q values estimated using these two methods are in concordance with each other. The Q data document attenuation anomalies in the layers in the gas hydrate stability zone above the bottom-simulating reflection (BSR) and in the free gas zone below. Changes in the attenuation anomalies correlate with small-scale fault systems in the Vestnesa Ridge suggesting a strong structural control on the distribution of free gas and gas hydrates in the region. We argued that high and spatially limited Q anomalies in the layer above the BSR indicate the presence of gas hydrates in marine sediments in this setting. Hence, our workflow to analyze Q using high-resolution P-cable 3D seismic data with a large bandwidth could be a potential technique to detect and directly map the distribution of gas hydrates in marine sediments.

High-resolution seismic velocity analysis by sign-based weighted semblance

Geophysics ◽  
2021 ◽  
pp. 1-35
Author(s):  
M. Javad Khoshnavaz
Keyword(s):  
Seismic Velocity ◽  
Resolution Enhancement ◽  
Synthetic Data ◽  
Correlation Coefficients ◽  
Velocity Model ◽  
Seismic Attenuation ◽  
Velocity Analysis ◽  
Velocity Models ◽  
Amplitude Versus Offset ◽  
Amplitude Versus

Building an accurate velocity model plays a vital role in routine seismic imaging workflows. Normal-moveout-based seismic velocity analysis is a popular method to make the velocity models. However, traditional velocity analysis methodologies are not generally capable of handling amplitude variations across moveout curves, specifically polarity reversals caused by amplitude-versus-offset anomalies. I present a normal-moveout-based velocity analysis approach that circumvents this shortcoming by modifying the conventional semblance function to include polarity and amplitude correction terms computed using correlation coefficients of seismic traces in the velocity analysis scanning window with a reference trace. Thus, the proposed workflow is suitable for any class of amplitude-versus-offset effects. The approach is demonstrated to four synthetic data examples of different conditions and a field data consisting a common-midpoint gather. Lateral resolution enhancement using the proposed workflow is evaluated by comparison between the results from the workflow and the results obtained by the application of conventional semblance and three semblance-based velocity analysis algorithms developed to circumvent the challenges associated with amplitude variations across moveout curves, caused by seismic attenuation and class II amplitude-versus-offset anomalies. According to the obtained results, the proposed workflow is superior to all the presented workflows in handling such anomalies.

Poststack velocity analysis by separation and imaging of seismic diffractions

Geophysics ◽  
2007 ◽  
Vol 72 (6) ◽  
pp. U89-U94 ◽  
Author(s):  
Sergey Fomel ◽  
Evgeny Landa ◽  
M. Turhan Taner
Keyword(s):  
High Resolution ◽  
Seismic Data ◽  
Field Data ◽  
Small Scale ◽  
Velocity Analysis ◽  
High Resolution Imaging ◽  
Diffracted Waves ◽  
Resolution Imaging ◽  
Seismic Reflections ◽  
Optimal Focusing

Small geologic features manifest themselves in seismic data in the form of diffracted waves, which are fundamentally different from seismic reflections. Using two field-data examples and one synthetic example, we demonstrate the possibility of separating seismic diffractions in the data and imaging them with optimally chosen migration velocities. Our criteria for separating reflection and diffraction events are the smoothness and continuity of local event slopes that correspond to reflection events. For optimal focusing, we develop the local varimax measure. The objectives of this work are velocity analysis implemented in the poststack domain and high-resolution imaging of small-scale heterogeneities. Our examples demonstrate the effectiveness of the proposed method for high-resolution imaging of such geologic features as faults, channels, and salt boundaries.

scholarly journals High-resolution coherency functionals for velocity analysis: An application for subbasalt seismic exploration

Geophysics ◽  
2013 ◽  
Vol 78 (5) ◽  
pp. U53-U63 ◽  
Author(s):  
Andrea Tognarelli ◽  
Eusebio Stucchi ◽  
Alessia Ravasio ◽  
Alfredo Mazzotti
Keyword(s):  
High Resolution ◽  
Seismic Data ◽  
Field Data ◽  
Signal To Noise Ratio ◽  
Synthetic Seismograms ◽  
Seismic Exploration ◽  
Velocity Analysis ◽  
Signal To Noise ◽  
Geologic Setting ◽  
Analytic Signals

We tested the properties of three different coherency functionals for the velocity analysis of seismic data relative to subbasalt exploration. We evaluated the performance of the standard semblance algorithm and two high-resolution coherency functionals based on the use of analytic signals and of the covariance estimation along hyperbolic traveltime trajectories. Approximate knowledge of the wavelet was exploited to design appropriate filters that matched the primary reflections, thereby further improving the ability of the functionals to highlight the events of interest. The tests were carried out on two synthetic seismograms computed on models reproducing the geologic setting of basaltic intrusions and on common midpoint gathers from a 3D survey. Synthetic and field data had a very low signal-to-noise ratio, strong multiple contamination, and weak primary subbasalt signals. The results revealed that high-resolution coherency functionals were more suitable than semblance algorithms to detect primary signals and to distinguish them from multiples and other interfering events. This early discrimination between primaries and multiples could help to target specific signal enhancement and demultiple operations.

Bringing slow slip processes into focus

2020 ◽  
Author(s):  
Rebecca Bell
Keyword(s):  
High Resolution ◽  
Subduction Zones ◽  
Seismic Velocity ◽  
Numerical Models ◽  
Fluid Pressure ◽  
Geophysical Methods ◽  
Rock Properties ◽  
Slow Slip ◽  
Velocity Models ◽  
The North

<p>The discovery of slow slip events (SSEs) at subduction margins in the last two decades has changed our understanding of how stress is released at subduction zones. Fault slip is now viewed as a continuum of different slip modes between regular earthquakes and aseismic creep, and an appreciation of seismic hazard can only be realised by understanding the full spectrum of slip. SSEs may have the potential to trigger destructive earthquakes and tsunami on faults nearby, but whether this is possible and why SSEs occur at all are two of the most important questions in earthquake seismology today. Laboratory and numerical models suggest that slow slip can be spontaneously generated under conditions of very low effective stresses, facilitated by high pore fluid pressure, but it has also been suggested that variations in frictional behaviour, potentially caused by very heterogeneous fault zone lithology, may be required to promote slow slip.</p><p>Testing these hypotheses is difficult as it requires resolving rock properties at a high resolution many km below the seabed sometimes in km’s of water, where drilling is technically challenging and expensive. Traditional geophysical methods like travel-time tomography cannot provide fine-scale enough velocity models to probe the rock properties in fault zones specifically. In the last decade, however, computational power has improved to the point where 3D full-waveform inversion (FWI) methods make it possible to use the full wavefield rather than just travel times to produce seismic velocity models with a resolution an order of magnitude better than conventional models. Although the hydrocarbon industry have demonstrated many successful examples of 3D FWI the method requires extremely high density arrays of instruments, very different to the 2D transect data collection style which is still commonly employed at subduction zones.</p><p> The north Hikurangi subduction zone, New Zealand is special, as it hosts the world’s most well characterised shallow SSEs (<2 km to 15 km below the seabed).  This makes it an ideal location to collect 3D data optimally for FWI to resolve rock properties in the slow slip zone. In 2017-2018 an unprecedentedly large 3D experiment including 3D multi-channel seismic reflection, 99 ocean bottom seismometers and 194 onshore seismometers was conducted along the north Hikurangi margin in an 100 km x 15 km area, with an average 2 km instrument spacing. In addition, IODP Expeditions 372 and 375 collected logging-while drilling and core data, and deployed two bore-hole observatories to target slow slip in the same area. In this presentation I will introduce you to this world class 3D dataset and preliminary results, which will enable high resolution 3D models of physical properties to be made to bring slow slip processes into focus.  </p>

An image of the Columbia Plateau from inversion of high‐resolution seismic data

Geophysics ◽  
1994 ◽  
Vol 59 (8) ◽  
pp. 1278-1289 ◽  
Author(s):  
William J. Lutter ◽  
Rufus D. Catchings ◽  
Craig M. Jarchow
Keyword(s):  
High Resolution ◽  
Seismic Data ◽  
Velocity Structure ◽  
Interface Velocity ◽  
Seismic Profile ◽  
Fault Zones ◽  
Reliable Estimate ◽  
Low Velocity ◽  
Depth Of Penetration ◽  
Image Blurring

We use a method of traveltime inversion of high‐resolution seismic data to provide the first reliable images of internal details of the Columbia River Basalt Group (CRBG), the subsurface basalt/sediment interface, and the deeper sediment/basement interface. Velocity structure within the basalts, delineated on the order of 1 km horizontally and 0.2 km vertically, is constrained to within ±0.1 km/s for most of the seismic profile. Over 5000 observed traveltimes fit our model with an rms error of 0.018 s. The maximum depth of penetration of the basalt diving waves (truncated by underlying low‐velocity sediments) provides a reliable estimate of the depth to the base of the basalt, which agrees with well‐log measurements to within 0.05 km (165 ft). We use image blurring, calculated from the resolution matrix, to estimate the aspect ratio of imaged velocity anomaly widths to true widths for velocity features within the basalt. From our calculations of image blurring, we interpret low velocity zones (LVZ) within the basalts at Boylston Mountain and the Whiskey Dick anticline to have widths of 4.5 and 3 km, respectively, within the upper 1.5 km of the model. At greater depth, the widths of these imaged LVZs thin to approximately 2 km or less. We interpret these linear, subparallel, low‐velocity zones imaged adjacent to anticlines of the Yakima Fold Belt to be brecciated fault zones. These fault zones dip to the south at angles between 15 to 45 degrees.

Sign in / Sign up

Export Citation Format

Share Document