Appraisal of gas-hydrate/free-gas from VP/VS ratio in the Makran accretionary prism

Existence of gas-hydrate in the marine sediments elevates both the P- and S-wave seismic velocities, whereas even a small amount of underlying free-gas decreases the P-wave velocity considerably and the S-wave velocity remains almost unaffected. Study of both P- and S-wave seismic velocities or their ratio (V_P/V_S) for the hydrate-bearing sediment provides more information than that obtained by the P- or S-wave velocity alone for the quantitative assessment of gas-hydrate. We estimate the P- and S-wave seismic velocities across a BSR (interface between gas-hydrate and free-gas bearing sediments) using the travel time inversion followed by a constrained AVA modeling of multi channel seismic (MCS) data at two locations in the Makran accretionary prism. Using this V_P/V_S ratio, we then quantify the amount of gas-hydrate and free-gas based on two rock-physics models. The result shows an estimate of 12–14.5% gas-hydrate and 4.5–5.5% free-gas of the pore volume based on first model, and 13–20% gas-hydrate and 3–3.5% free-gas of the pore volume based on the second model, respectively.     

Rock physics modeling for assessing gas hydrate and free gas: a case study in the Cascadia accretionary prism

We investigate the estimation of gas hydrate and free gas concentration using various rock physics models in the Cascadia accretionary prism, which is one of the most intensively studied regions of natural gas hydrate occurrences. Surface seismic reflection data is the most useful and cost-effective in deriving seismic velocity, and hence estimating gas hydrate and free gas across a BSR with depth, if a proper background (without gas hydrate and free gas) velocity is chosen. We have used effective medium theory of Helgerud et al. (EMTH) and, a combination of self-consistent approximation and differential effective medium (SCA-DEM) theory coupled with smoothing approximation for crystalline aggregate. Using the SCA-DEM (non-load-bearing) and EMTH (load-bearing) modeling, we calculate the average saturations of gas hydrate as 17 and 19%, respectively within ~100 m thick sedimentary column using velocity, derived from the surface seismic data. The saturations of gas hydrate are estimated as 15 and 18% using the SCA-DEM, and 20 and 25% using EMTH from the logging-while-drilling and wire-line sonic velocities, respectively. Estimations of gas hydrate from Poisson’s ratio are in average 50% for EMTH and 10% for SCA-DEM theory. We obtain the maximum saturation of free gas as 1–2% by employing the SCA-DEM theory either to seismic or sonic velocities, whereas the free-gas saturation varies between 0.1 and 0.4% for EMTH model. The gas hydrate saturation estimated from the sonic velocity and the free gas saturation derived from both the seismic and sonic velocities using the SCA-DEM modeling match quite well with those determined from the pressure core data in the study region.     

Amplitude blanking related to the pore-filling of gas hydrate in sediments

Seismic indicators of gas-hydrate-bearing sediments include elevated interval velocities and amplitude reduction of seismic reflections owing to the presence of gas hydrate in the sediment's pore spaces. However, large amplitude blanking with relatively low interval velocities observed at the Blake Ridge has been enigmatic because realistic seismic models were absent to explain the observation. This study proposes models in which the gas hydrate concentrations vary in proportion to the porosity. Where gas hydrate concentrations are greater in more porous media, a significant amplitude blanking can be achieved with relatively low interval velocity. Depending on the amount of gas hydrate concentration in the pore space, reflection amplitudes from hydrate-bearing sediments can be much less, less or greater than those from corresponding non-hydrate-bearing sediments. This revised version was published online in November 2006 with corrections to the Cover Date.     

Seismic velocities and quantification of gas-hydrates from AVA modeling in the western continental margin of India

The most commonly used marker for the investigation of gas-hydrates is the bottom simulating reflector (BSR), which is caused by gas-hydrate laden sediment underlain by either brine or gas-saturated sediment. A BSR has been identified by seismic experiment in the Kerala-Konkan Basin of the western continental margin of India. Here we perform AVA modeling of seismic reflection data from a BSR to investigate the seismic velocities for quantitative assessment of gas-hydrates and to understand the origin of the BSR. The result reveals a P-wave velocity of 2.245 km/s and an S-wave velocity of 0.895 km/s for the sediments above the BSR. This corresponds to a Poisson ratio of 0.406 and hydrates saturation of ∼30% in the study area. The comparison of estimated P-wave velocity (1.77 km/s) above the hydrates-bearing sediment to that (1.78 km/s) below the BSR implies that the origin of the BSR is mainly due to gas-hydrates, as the presence (even in small quantities) of free-gas reduces the P-wave velocity considerably.     

Tying seismic data to geologic information from core data: an example from ODP Leg 177

The integration of seismic data with core data should provide ground-truth to a structural interpretation of seismic data. The main difficulty in such an integration effort is the correct translation of physical property measurements on cores to a form which can be used in seismostratigraphic interpretation. In the absence of down-hole well data and check-shots, required knowledge of the velocity structure at the drilling locations can be obtained directly from measurements of the physical properties of core samples. This involves upscaling of the data from physical properties of cores to the sample interval used in the seismic data. In the present study, three of the seven drill-sites of ODP (Ocean Drilling Program) Leg 177 in 1997/1998, located on the Agulhas Ridge in the south-eastern Atlantic (sites 1088–1090), were connected with eight seismic profiles. Physical properties data measured on the cores from the various holes at each site were combined to create a single continuous log and used to construct synthetic seismograms. The synthetics generally show a good agreement with real seismic data in terms of amplitude and waveform. Some reflections in these generated traces may have a time-shift due to sections with incomplete or spurious P-wave velocity measurements in the ODP datasets. The main reflectors identified in the real seismic data correspond to hiatuses or periods of reduced sedimentation rates, and correlate well with density variations. One particular hiatus, clearly observable in the real seismic data, was not unequivocally identifiable in the various types of core data, and tying core data to seismic data can confirm its existence in the core data, showing the benefit of including seismic data in an interpretation of core log data. On the other hand, core data provide a calibration tool for the geological timescale of seismic data and information about the lithology, needed in the interpretation of seismic data.     

南海陆坡天然气水合物饱和度估计

基于双相介质理论和热弹性理论,建立了沉积层纵波速度与天然气水合物饱和度、弹性性质及地层孔隙度之间的关系。通过对比饱和水的理论P波速度与实际P波速度,可以得到天然气水合物饱和度。根据ODP184航次的电阻率、声波速度、密度等测井资料以及地质资料,初步推断南海陆坡存在天然气水合物。根据声波测井的纵波速度估算出南海1146和1148井天然气水合物饱和度分别为孔隙空间的25％～30％和10％～20％,1148井个别沉积层天然气水合物饱和度可达40％～50％。沉积层的纵波速度与饱和水速度差值越大,天然气水合物饱和度越高。     

Volcanic stratigraphy of DSDP/ODP Hole 395A: An interpretation using well-logging data

Deep Sea Drilling Project/Ocean Drilling Program Hole 395A was drilled approximately 500 m deep into young oceanic crust west of the Mid-Atlantic Ridge. Core recovery is very poor in this hole and therefore continuous downhole measurements are important to understand the drilled lithology. Geophysical downhole measurements were carried out during several cruises. A new set of logs was recorded during Leg 174B in summer 1997. The new logging data show a significant improvement in data quality compared to older measurements from Leg 109. The lithostratigraphy established from cores gives only limited information because of the poor core recovery in this hole. The gaps in the core lithostratigraphy are filled by reconstructing a synthetic lithological profile using the standard well-logging data. Three types of lava morphologies, massive basalts, altered lava flows, and pillow basalts, may be distinguished using the logs because the lava morphologies show differences in their physical properties due to differences in fracturing and alteration. The synthetic lithological profile gives a more detailed and precise vertical definition of single layers than the core profile. The integration of further logging and core data enables a detailed reconstruction of the accretion history at the drill site. Cyclic, upward decreasing trends in the resistivity logs were already observed during earlier cruises and were referred to magmatic cycles. Similar trends occur in the density log and, inversely, in the total gamma ray log. The trends reflect gradual changes in fracturing, porosity, permeability, and alteration and cover depth intervals of several tens of meters. Boundaries between cycles are interpreted to correspond to periods of volcanic quiescence. Two types of boundaries may be identified. Boundaries correlating with reversals in the magnetic field and/or changes in the geochemical composition of the basalts are interpreted as long pauses. Basalts separated by these boundaries were probably fed by separate magma reservoirs. Boundaries identified only by changes in alteration but not in geochemistry are interpreted to represent shorter pauses. They separate basalts that were probably fed by the same magma chamber. This revised version was published online in November 2006 with corrections to the Cover Date.     

Effective medium modeling to assess gas hydrate and free-gas evident from the velocity structure in the Makran accretionary prism, offshore Pakistan

Seismic properties of sediments are strongly influenced by pore fluids. Stiffness of unconsolidated marine sediment increases with the presence of gas hydrate and decreases with the presence of gas. A strong bottom-simulating reflector (BSR) observed on a seismic profile in the Makran accretionary prism, offshore Pakistan, indicates the presence of gas hydrate and free-gas across the BSR. Elastic properties of gas depend largely on pressure and temperature. We, therefore, first determine the elastic modulus of gas at pressure and temperature calculated at the BSR depth in the study region. The interval velocities derived from the seismic data are interpreted by the effective medium theory, which is a combination of self-consistent approximation and differential effective medium theories, together with a smoothing approximation, for assessment of gas hydrate and free-gas. The results show the saturations of gas hydrate and free-gas as 22 and 2.4% of pore space, respectively, across the BSR.     

Unique problems associated with seismic analysis of partially gas-saturated unconsolidated sediments

Gas hydrate stability conditions restrict the occurrence of gas hydrate to unconsolidated and high water-content sediments at shallow depths. Because of these host sediments properties, seismic and well log data acquired for the detection of free gas and associated gas hydrate-bearing sediments often require nonconventional analysis. For example, a conventional method of identifying free gas using the compressional/shear-wave velocity (V_p/V_s) ratio at the logging frequency will not work, unless the free-gas saturations are more than about 40%. The P-wave velocity dispersion of partially gas-saturated sediments causes a problem in interpreting well log velocities and seismic data. Using the White, J.E. [1975. Computed seismic speeds and attenuation in rocks with partial gas saturation. Geophysics 40, 224–232] model for partially gas-saturated sediments, the difference between well log and seismic velocities can be reconciled. The inclusion of P-wave velocity dispersion in interpreting well log data is, therefore, essential to identify free gas and to tie surface seismic data to synthetic seismograms.     

Downhole well log and core montages from the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope

The BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well was an integral part of an ongoing project to determine the future energy resource potential of gas hydrates on the Alaska North Slope. As part of this effort, the Mount Elbert well included an advanced downhole geophysical logging program. Because gas hydrate is unstable at ground surface pressure and temperature conditions, a major emphasis was placed on the downhole-logging program to determine the occurrence of gas hydrates and the in-situ physical properties of the sediments. In support of this effort, well-log and core data montages have been compiled which include downhole log and core-data obtained from the gas-hydrate-bearing sedimentary section in the Mount Elbert well. Also shown are numerous reservoir parameters, including gas-hydrate saturation and sediment porosity log traces calculated from available downhole well log and core data.     

High-resolution well-log derived dielectric properties of gas-hydrate-bearing sediments, Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope

A dielectric logging tool, electromagnetic propagation tool (EPT), was deployed in 2007 in the BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well (Mount Elbert Well), North Slope, Alaska. The measured dielectric properties in the Mount Elbert well, combined with density log measurements, result in a vertical high-resolution (cm-scale) estimate of gas hydrate saturation. Two hydrate-bearing sand reservoirs about 20 m thick were identified using the EPT log and exhibited gas-hydrate saturation estimates ranging from 45% to 85%. In hydrate-bearing zones where variation of hole size and oil-based mud invasion are minimal, EPT-based gas hydrate saturation estimates on average agree well with lower vertical resolution estimates from the nuclear magnetic resonance logs; however, saturation and porosity estimates based on EPT logs are not reliable in intervals with substantial variations in borehole diameter and oil-based invasion.EPT log interpretation reveals many thin-bedded layers at various depths, both above and below the thick continuous hydrate occurrences, which range from 30-cm to about 1-m thick. Such thin layers are not indicated in other well logs, or from the visual observation of core, with the exception of the image log recorded by the oil-base microimager. We also observe that EPT dielectric measurements can be used to accurately detect fine-scale changes in lithology and pore fluid properties of hydrate-bearing sediments where variation of hole size is minimal. EPT measurements may thus provide high-resolution in-situ hydrate saturation estimates for comparison and calibration with laboratory analysis.     

Structure and sediment distribution in the western Bering Sea

Eleven seismic reflection profiles across Shirshov Ridge and the adjacent deep-water sedimentary basins (Komandorsky and Aleutian Basins) are presented to illustrate the sediment distribution in the western Bering Sea. A prominent seismic reflecting horizon, Reflector P (Middle—Late Miocene in age), is observed throughout both the Aleutian and Komandorsky Basins at an approximate subbottom depth of 1 km. This reflector is also present, in places, on the flanks and along the crest of Shirshov Ridge. The thickness of sediments beneath Reflector P is significantly different within the two abyssal basins. In the Aleutian Basin, the total subbottom depth to acoustic basement (basalt?) is about 4 km, while in the Komandorsky Basin the depth is about 2 km.Shirshov Ridge, a Cenozoic volcanic feature that separates the Aleutian and Komandorsky Basins, is an asymmetric bathymetric ridge characterized by thick sediments along its eastern flank and steep scarps on its western side. The southern portion of the ridge has more structural relief that includes several deep, sediment-filled basins along its summit.Velocity data from sonobuoy measurements indicate that acoustic basement in the Komandorsky Basin has an average compressional wave velocity of 5.90 km/sec. This value is considerably larger than the velocities measured for acoustic basement in the northwestern Aleutian Basin (about 5.00 km/sec) and in the central Aleutian Basin (5.40–5.57 km/sec). In the northwestern Aleutian Basin, the low-velocity acoustic basement may be volcaniclastic sediments or other indurated sediments that are overlying true basaltic basement. A refracting horizon with similar velocities (4.6–5.0 km/sec) as acoustic basement dips steeply beneath the Siberian continental margin, reaching a maximum subbottom depth of about 8 km. The thick welt of sediment at the base of the Siberian margin may be the result of sediment loading or tectonic depression prior to Late Cenozoic time.     

基于OBS和MCS数据的水合物地层速度特征

天然气水合物是一种新型的清洁能源, 南海北部神狐海域的地质条件有利于水合物的形成和储藏。传统的多道地震(MCS)数据难以得到精确的速度信息, 并且只能从时间域上判断地质体纵向分布。海底地震仪(OBS)是一种常用的主动源地震仪器, 可以接收到更清晰的气枪信号。相比于MCS, OBS剖面上的折射震相可以揭示较深部的地层速度信息。文章结合MCS和OBS的优势, 识别水平叠加剖面上的反射层位, 并得到初始模型; 将OBS剖面和水平叠加剖面拼合, 从而判断OBS剖面上反射震相所对应层位; 拾取OBS台站上的反射和折射震相, 使用RayInvr软件正演模拟得到水合物存在区域的二维速度模型, 解决了MCS中较为困难的时深转换问题。最终模型显示了水合物、游离气区域的埋深、厚度和速度, 以及似海底反射(BSR)下方更深部界面的深度和速度特征。     

A seismic velocity inversion on Bjørnøya — The western Barents Shelf

Seismic velocity measurements on the exposed Triassic to Cambro-Silurian sedimentary sequence on Bjørnøya show a velocity inversion in the stratigraphic section with a minimum velocity (3.8 km/s) in the Devonian fluvial sandstones of the Røedvika Formation, which also has the highest porosity (～ 17%). The low porosity in the younger formations is due to silicification and dolomitization probably by chemical alterations of the pore water rather than an overburden effect. Sonobuoy measurements in the vicinity of Bjørnøya show seismic velocities at the sea floor comparable with velocities of the exposed formations on the island itself.     

The ODP color digital imaging system: Color logs of quaternary sediments from the Santa Barbara basin,site 893

Digital color logs of cores from leg 146, holes 893A and 893B, have been taken from images captured during January 1993, within days after the cores were split and described. The images were captured and color analyses performed on the Ocean Drilling Program (ODP) color digital imaging system, which was assembled from relatively inexpensive, off‐the‐shelf components. The images were used to calculate sedimentation rates by fitting chronological data from hole 893A to void‐corrected depths determined by eliminating all voids mapped from the images as >1 cm in length measured downcore. Color measurements were made at intervals between 0.22 and 1.0 mm in length, and then Commission Internationale de l'Eclairage (CIE) 1931 chromaticity values were computed. Results plotted within CIE chro‐maticity space lie in close proximity to the Munsell 5Y hue plane, confirming that the instrumental color analysis technique has produced results consistent with those of the human observers who described colors visually. Significant periodicities occurring at 12, 17, 31, and 90 years correlate with sunspot activity cycles, suggesting that color variations may reflect global climatic forcing functions. Linear correlation between color data sets from holes 893A and 893B suggests that as much as 1.2 m of material present at or near the top of hole 893B may not have been recovered from hole 893A, and that there is a 70‐cm depth discrepancy in the opposite direction at 51.5 m below seafloor (mbsf). The article recommends that the 1931 CIE chromaticity system be used routinely for describing colors of geological materials because it readily accommodates the mathematical manipulations required for statistical and time‐series analyses and avoids the subjectivity and other weaknesses inherent in the Munsell Color System.     

Pre- and post-drill comparison of the Mount Elbert gas hydrate prospect, Alaska North Slope

In 2006, the United States Geological Survey (USGS) completed a detailed analysis and interpretation of available 2-D and 3-D seismic data, along with seismic modeling and correlation with specially processed downhole well log data for identifying potential gas hydrate accumulations on the North Slope of Alaska. A methodology was developed for identifying sub-permafrost gas hydrate prospects within the gas hydrate stability zone in the Milne Point area. The study revealed a total of 14 gas hydrate prospects in this area.In order to validate the gas hydrate prospecting protocol of the USGS and to acquire critical reservoir data needed to develop a longer-term production testing program, a stratigraphic test well was drilled at the Mount Elbert prospect in the Milne Point area in early 2007. The drilling confirmed the presence of two prominent gas-hydrate-bearing units in the Mount Elbert prospect, and high quality well logs and core data were acquired. The post-drill results indicate pre-drill predictions of the reservoir thickness and the gas-hydrate saturations based on seismic and existing well data were 90% accurate for the upper unit (hydrate unit D) and 70% accurate for the lower unit (hydrate unit C), confirming the validity of the USGS approach to gas hydrate prospecting. The Mount Elbert prospect is the first gas hydrate accumulation on the North Slope of Alaska identified primarily on the basis of seismic attribute analysis and specially processed downhole log data. Post-drill well log data enabled a better constraint of the elastic model and the development of an improved approach to the gas hydrate prospecting using seismic attributes.     

Pore- and fracture-filling gas hydrate reservoirs in the Gulf of Mexico Gas Hydrate Joint Industry Project Leg II Green Canyon 955 H well

High-quality logging-while-drilling (LWD) downhole logs were acquired in seven wells drilled during the Gulf of Mexico Gas Hydrate Joint Industry Project Leg II in the spring of 2009. Well logs obtained in one of the wells, the Green Canyon Block 955 H well (GC955-H), indicate that a 27.4-m thick zone at the depth of 428 m below sea floor (mbsf; 1404 feet below sea floor (fbsf)) contains gas hydrate within sand with average gas hydrate saturations estimated at 60% from the compressional-wave (P-wave) velocity and 65% (locally more than 80%) from resistivity logs if the gas hydrate is assumed to be uniformly distributed in this mostly sand-rich section. Similar analysis, however, of log data from a shallow clay-rich interval between 183 and 366 mbsf (600 and 1200 fbsf) yielded average gas hydrate saturations of about 20% from the resistivity log (locally 50−60%) and negligible amounts of gas hydrate from the P-wave velocity logs. Differences in saturations estimated between resistivity and P-wave velocities within the upper clay-rich interval are caused by the nature of the gas hydrate occurrences. In the case of the shallow clay-rich interval, gas hydrate fills vertical (or high angle) fractures in rather than filling pore space in sands. In this study, isotropic and anisotropic resistivity and velocity models are used to analyze the occurrence of gas hydrate within both the clay-rich and sand dominated gas-hydrate-bearing reservoirs in the GC955-H well.     

Gas hydrate saturation from acoustic impedance and resistivity logs in the Shenhu area, South China Sea

During the China’s first gas hydrate drilling expedition -1 (GMGS-1), gas hydrate was discovered in layers ranging from 10 to 25 m above the base of gas hydrate stability zone in the Shenhu area, South China Sea. Water chemistry, electrical resistivity logs, and acoustic impedance were used to estimate gas hydrate saturations. Gas hydrate saturations estimated from the chloride concentrations range from 0 to 43% of the pore space. The higher gas hydrate saturations were present in the depth from 152 to 177 m at site SH7 and from 190 to 225 m at site SH2, respectively. Gas hydrate saturations estimated from the resistivity using Archie equation have similar trends to those from chloride concentrations. To examine the variability of gas hydrate saturations away from the wells, acoustic impedances calculated from the 3 D seismic data using constrained sparse inversion method were used. Well logs acquired at site SH7 were incorporated into the inversion by establishing a relation between the water-filled porosity, calculated using gas hydrate saturations estimated from the resistivity logs, and the acoustic impedance, calculated from density and velocity logs. Gas hydrate saturations estimated from acoustic impedance of seismic data are ∼10-23% of the pore space and are comparable to those estimated from the well logs. The uncertainties in estimated gas hydrate saturations from seismic acoustic impedances were mainly from uncertainties associated with inverted acoustic impedance, the empirical relation between the water-filled porosities and acoustic impedances, and assumed background resistivity.