Estimation of gas hydrate and free gas saturation, concentration, and distribution from seismic data

Geophysics ◽  
2002 ◽  
Vol 67 (2) ◽  
pp. 582-593 ◽  
Author(s):  
Shaoming Lu ◽  
George A. McMechan

Gas hydrates contain a major untapped source of energy and are of potential economic importance. The theoretical models to estimate gas hydrate saturation from seismic data predict significantly different acoustic/elastic properties of sediments containing gas hydrate; we do not know which to use. Thus, we develop a new approach based on empirical relations. The water‐filled porosity is calibrated (using well‐log data) to acoustic impedance twice: one calibration where gas hydrate is present and the other where free gas is present. The water‐filled porosity is used in a combination of Archie equations (with corresponding parameters for either gas hydrate or free gas) to estimate gas hydrate or free gas saturations. The method is applied to single‐channel seismic data and well logs from Ocean Drilling Program leg 164 from the Blake Ridge area off the east coast of North America. The gas hydrate above the bottom simulating reflector (BSR) is estimated to occupy ∼3–8% of the pore space (∼2–6% by volume). Free gas is interpreted to be present in three main layers beneath the BSR, with average gas saturations of 11–14%, 7–11%, and 1–5% of the pore space (6–8%, 4–6%, and 1–3% by volume), respectively. The estimated saturations of gas hydrate are very similar to those estimated from vertical seismic profile data and generally agree with those from independent, indirect estimates obtained from resistivity and chloride measurements. The estimated free gas saturations agree with measurements from a pressure core sampler. These results suggest that locally derived empirical relations between porosity and acoustic impedance can provide cost‐effective estimates of the saturation, concentration, and distribution of gas hydrate and free gas away from control wells.

Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. B125-B134 ◽  
Author(s):  
Xiujuan Wang ◽  
Myung Lee ◽  
Shiguo Wu ◽  
Shengxiong Yang

Wireline logs were acquired in eight wells during China’s first gas hydrate drilling expedition (GMGS-1) in April–June of 2007. Well logs obtained from site SH3 indicated gas hydrate was present in the depth range of 195–206 m below seafloor with a maximum pore-space gas hydrate saturation, calculated from pore water freshening, of about 26%. Assuming gas hydrate is uniformly distributed in the sediments, resistivity calculations using Archie’s equation yielded hydrate-saturation trends similar to those from chloride concentrations. However, the measured compressional (P-wave) velocities decreased sharply at the depth between 194 and 199 mbsf, dropping as low as [Formula: see text], indicating the presence of free gas in the pore space, possibly caused by the dissociation of gas hydrate during drilling. Because surface seismic data acquired prior to drilling were not influenced by the in situ gas hydrate dissociation, surface seismic data could be used to identify the cause of the low P-wave velocity observed in the well log. To determine whether the low well-log P-wave velocity was caused by in situ free gas or by gas hydrate dissociation, synthetic seismograms were generated using the measured well-log P-wave velocity along with velocities calculated assuming both gas hydrate and free gas in the pore space. Comparing the surface seismic data with various synthetic seismograms suggested that low P-wave velocities were likely caused by the dissociation of in situ gas hydrate during drilling.


Geophysics ◽  
2010 ◽  
Vol 75 (1) ◽  
pp. B35-B46 ◽  
Author(s):  
Ashwani Dev ◽  
George A. McMechan

The effect of the structural setting on gas hydrate and free-gas accumulation near the East Casey fault zone at Keathley Canyon in the Gulf of Mexico is investigated using well logs and 2D and 3D acoustic-impedance inversions. We interpret two zones from the well logs: a probable hydrate zone and a dissociating hydrate zone. We interpret hydrates as fracture fillings in the clay-dominated formation with maximum saturation of approximately 30% of the available pore space, and maximum volumetric concentration of approximately 12%. The maximum free-gas concentration below the interpreted bottom simulating reflector (BSR) is less than 4%. Resistivity measurements alone are incon-clusive; complimentary measurements, such as sonic, reduce hy-drate interpretation ambiguities. Seismic data in 3D and the derived acoustic-impedance volume are interpreted in terms of a BSR, a high [Formula: see text]-impedance hydrate zone, and free gas. We interpret two bright, low [Formula: see text]-impedance features terminating at the BSR as free-gas accumulations. The contrast in average [Formula: see text]-impedance across the fault suggests a change in lithology across the fault zone, and free-gas transport along, and to the west of, the fault. Variations in depths of interpreted free-gas zones suggest either a significant lateral variation in the hydrate stability across the fault zone, or coexistence of hydrates and free gas within the hydrate stability zone, or both. The dynamics of the tectonic environment imply nonequilibrium conditions of time-dependent temperature, pressure, or fluid transport.


Geophysics ◽  
2013 ◽  
Vol 78 (3) ◽  
pp. D169-D179 ◽  
Author(s):  
Zijian Zhang ◽  
De-hua Han ◽  
Daniel R. McConnell

Hydrate-bearing sands and shallow nodular hydrate are potential energy resources and geohazards, and they both need to be better understood and identified. Therefore, it is useful to develop methodologies for modeling and simulating elastic constants of these hydrate-bearing sediments. A gas-hydrate rock-physics model based on the effective medium theory was successfully applied to dry rock, water-saturated rock, and hydrate-bearing rock. The model was used to investigate the seismic interpretation capability of hydrate-bearing sediments in the Gulf of Mexico by computing elastic constants, also known as seismic attributes, in terms of seismic interpretation, including the normal incident reflectivity (NI), Poisson’s ratio (PR), P-wave velocity ([Formula: see text]), S-wave velocity ([Formula: see text]), and density. The study of the model was concerned with the formation of gas hydrate, and, therefore, hydrate-bearing sediments were divided into hydrate-bearing sands, hydrate-bearing sands with free gas in the pore space, and shallow nodular hydrate. Although relations of hydrate saturation versus [Formula: see text] and [Formula: see text] are different between structures I and II gas hydrates, highly concentrated hydrate-bearing sands may be interpreted on poststack seismic amplitude sections because of the high NI present. The computations of elastic constant implied that hydrate-bearing sands with free gas could be detected with the crossplot of NI and PR from prestack amplitude analysis, and density may be a good hydrate indicator for shallow nodular hydrate, if it can be accurately estimated by seismic methods.


2020 ◽  
Author(s):  
Malin Waage ◽  
Stefan Bünz ◽  
Kate Waghorn ◽  
Sunny Singhorha ◽  
Pavel Serov

<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>


Geophysics ◽  
2000 ◽  
Vol 65 (2) ◽  
pp. 565-573 ◽  
Author(s):  
Christine Ecker ◽  
Jack Dvorkin ◽  
Amos M. Nur

Marine seismic data and well‐log measurements at the Blake Ridge offshore South Carolina show that prominent seismic bottom‐simulating reflectors (BSRs) are caused by sediment layers with gas hydrate overlying sediments with free gas. We apply a theoretical rock‐physics model to 2-D Blake Ridge marine seismic data to determine gas‐hydrate and free‐gas saturation. High‐porosity marine sediment is modeled as a granular system where the elastic wave velocities are linked to porosity; effective pressure; mineralogy; elastic properties of the pore‐filling material; and water, gas, and gas‐hydrate saturation of the pore space. To apply this model to seismic data, we first obtain interval velocity using stacking velocity analysis. Next, all input parameters to the rock‐physics model, except porosity and water, gas, and gas hydrate saturation, are estimated from geologic information. To estimate porosity and saturation from interval velocity, we first assume that the entire sediment does not contain gas hydrate or free gas. Then we use the rock‐physics model to calculate porosity directly from the interval velocity. Such porosity profiles appear to have anomalies where gas hydrate and free gas are present (as compared to typical profiles expected and obtained in sediment without gas hydrate or gas). Porosity is underestimated in the hydrate region and is overestimated in the free‐gas region. We calculate the porosity residuals by subtracting a typical porosity profile (without gas hydrate and gas) from that with anomalies. Next we use the rock‐physics model to eliminate these anomalies by introducing gas‐hydrate or gas saturation. As a result, we obtain the desired 2-D saturation map. The maximum gas‐hydrate saturation thus obtained is between 13% and 18% of the pore space (depending on the version of the model used). These saturation values are consistent with those measured in the Blake Ridge wells (away from the seismic line), which are about 12%. Free‐gas saturation varies between 1% and 2%. The saturation estimates are extremely sensitive to the input velocity values. Therefore, accurate velocity determination is crucial for correct reservoir characterization.


Geophysics ◽  
2010 ◽  
Vol 75 (2) ◽  
pp. C1-C6 ◽  
Author(s):  
Maheswar Ojha ◽  
Kalachand Sain ◽  
Timothy A. Minshull

We estimate the saturations of gas hydrate and free gas based on measurements of seismic-reflection amplitude variation with offset (AVO) for a bottom-simulating reflector coupled with rock-physics modeling. When we apply the approach to data from a seismic line in the Makran accretionary prism in the Arabian Sea, the results reveal lateral variations of gas-hydrate and free-gas saturations of 4–29% and 1–7.5%, respectively, depending on the rock-physics model used to relate seismic velocity to saturation. Our approach is simple and easy to implement.


2021 ◽  
pp. 1-60
Author(s):  
Darrell A. Terry ◽  
Camelia C. Knapp

The presence of marine gas hydrates is routinely inferred based on the identification of bottom simulating reflectors (BSRs) in common depth-point (CDP) seismic images. Additional seismic studies such as amplitude variation with offset (AVO) analysis can be applied for corroboration. Though confirmation is needed by drilling and sampling, seismic analysis has proven to be a cost-effective approach to identify the presence of marine gas hydrates. Single channel far offset seismic images are investigated for what appears to be a more reliable and cost-effective indicator for the presence of bottom simulating reflectors than traditional CDP processing or AVO analysis. A non-traditional approach to processing seismic data is taken to be more relevant to imaging the gas/gas hydrate contact. Instead of applying the traditional CDP seismic processing workflows from the oil industry, we more carefully review the significant amount of information existing in the data to explore how the character of the data changes as offset angle increases. Three cases from different environments are selected for detailed analysis. These include 1) stratigraphy running parallel with the ocean bottom; 2) a potential bottom simulating reflector, running parallel to the ocean bottom, and cross-cutting dipping reflections, and 3) a suspected thermal intrusion without a recognizable bottom simulating reflector. This investigation considers recently collected multi-channel seismic data from the deep waters of the central Aleutian Basin beneath the Bering Sea, the pre-processing of the data sets, and the methodology for processing and display to generate single channel seismic images. Descriptions are provided for the single channel near and far offset seismic images for the example cases. Results indicate that BSRs related to marine gas hydrates, and originating due to the presence of free gas, are more easily and uniquely identifiable from single channel displays of far offset seismic images than from traditional CDP displays.


Sign in / Sign up

Export Citation Format

Share Document