Evidence of gas hydrate accumulation and its resource estimation in Andaman deep water basin from seismic and well log data

2D and 3D seismic reflection and well log data from Andaman deep water basin are analyzed to investigate geophysical evidence related to gas hydrate accumulation and saturation. Analysis of seismic data reveals the presence of a bottom simulating reflector (BSR) in the area showing all the characteristics of a classical BSR associated with gas hydrate accumulation. Double BSRs are also observed on some seismic sections of area (Area B) that suggest substantial changes in pressure–temperature (P–T) conditions in the past. The manifestation of changes in P–T conditions can also be marked by the varying gas hydrate stability zone thickness (200–650 m) in the area. The 3D seismic data of Area B located in the ponded fill, west of Alcock Rise has been pre-stack depth migrated. A significant velocity inversion across the BSR (1,950–1,650 m/s) has been observed on the velocity model obtained from pre-stack depth migration. The areas with low velocity of the order of 1,450 m/s below the BSR and high amplitudes indicate presence of dissociated or free gas beneath the hydrate layer. The amplitude variation with offset analysis of BSR depicts increase in amplitude with offset, a similar trend as observed for the BSR associated with the gas hydrate accumulations. The presence of gas hydrate shown by logging results from a drilled well for hydrocarbon exploration in Area B, where gas hydrate deposit was predicted from seismic evidence, validate our findings. The base of the hydrate layer derived from the resistivity and acoustic transit-time logs is in agreement with the depth of hydrate layer interpreted from the pre-stack depth migrated seismic section. The resistivity and acoustic transit-time logs indicate 30-m-thick hydrate layer at the depth interval of 1,865–1,895 m with 30 % hydrate saturation. The total hydrate bound gas in Area B is estimated to be 1.8 × 10¹⁰ m³, which is comparable (by volume) to the reserves in major conventional gas fields.     

Geothermal modeling of the gas hydrate stability zone along the Krishna Godavari Basin

A wide-spread bottom simulating reflector (BSR), interpreted to mark the thermally controlled base of the gas hydrate stability zone, is observed over a close grid of multichannel seismic profiles in the Krishna Godavari Basin of the eastern continental margin of India. The seismic data reveal that gas hydrate occurs in the Krishna Godavari Basin at places where water depths exceed 850 m. The thickness of the gas hydrate stability zone inferred from the BSR ranges up to 250 m. A conductive model was used to determine geothermal gradients and heat flow. Ground truth for the assessment and constraints on the model were provided by downhole measurements obtained during the National Gas Hydrate Program Expedition 01 of India at various sites in the Krishna Godavari Basin. Measured downhole temperature gradients and seafloor-temperatures, sediment thermal conductivities, and seismic velocity are utilized to generate regression functions for these parameters as function of overall water depth. In the first approach the base of gas hydrate stability is predicted from seafloor bathymetry using these regression functions and heat flow and geothermal gradient are calculated. In a second approach the observed BSR depth from the seismic profiles (measured in two-way travel time) is converted into heat flow and geothermal gradient using the same ground-truth data. The geothermal gradient estimated from the BSR varies from 27 to 67°C/km. Corresponding heat flow values range from 24 to 60 mW/m². The geothermal modeling shows a close match of the predicted base of the gas hydrate stability zone with the observed BSR depths.     

AVO analysis of BSR to assess free gas within fine-grained sediments in the Shenhu area,South China Sea

Gas hydrates have been identified from two-dimensional (2D) seismic data and logging data above bottom simulating reflector (BSR) during China’s first gas hydrate drilling expedition in 2007. The multichannel reflection seismic data were processed to be preserved amplitudes for quantitatively analyzing amplitude variation with offset (AVO) at BSRs. Low P-wave velocity anomaly below BSR, coinciding with high amplitude reflections in 2D seismic data, indicates the presence of free gas. The absolute values of reflection coefficient versus incidence angles for BSR range from 0 to 0.12 at different CMPs near Site SH2. According to logging data and gas hydrate saturations estimated from resistivity of Site SH2, P-wave velocities calculated from effective media theory (EMT) fit the measured sonic velocities well and we choose EMT to calculate elastic velocities for AVO. The rock-physics modeling and AVO analysis were combined to quantitatively assess free gas saturations and distribution by the reflection coefficients variation of the BSRs in Shenhu area, South China Sea. AVO estimation indicates that free gas saturations immediately beneath BSRs may be about 0.2 % (uniform distribution) and up to about 10 % (patchy distribution) at Site SH2.     

Multiple seismic attribute analyses for determination of bottom simulating reflector of gas hydrate seismic data in the Ulleung Basin of Korea

The bottom simulating reflector (BSR), the boundary between the gas hydrate and the free gas zone, is considered to be the most common evidence in seismic data analysis for gas hydrate exploration. Multiple seismic attribute analyses of reflectivity and acoustic impedance from the post-stack deconvolution and complex analysis of instantaneous attribute properties including the amplitude envelope, instantaneous frequency, phase, and first derivative of the amplitude of seismic data have been used to effectively confirm the existence of a BSR as the base of gas hydrate stability zone. In this paper, we consider individual seismic attribute analysis and integrate the results of those attributes to locate the position of the BSR. The outputs from conventional seismic data processing of the gas hydrate data set in the Ulleung Basin were used as inputs for multiple analyses. Applying multiple attribute analyses to the individual seismic traces showed that the identical anomalies found in two-way travel time (TWT) between 3.1 and 3.2 s from the results of complex analyses and l ₁ norm deconvolution indicated the location of the BSR.     

Canyon-infilling and gas hydrate occurrences in the frontal fold of the offshore accretionary wedge off southern Taiwan

We utilized reflection seismic and bathymetric data to infer the canyon-infilling, fold uplift, and gas hydrate occurrences beneath the frontal fold at the toe of the accretionary wedge, offshore SW Taiwan. The lateral migrating paleo-Penghu canyons has cut across the frontal fold with six distinct canyon/channel incisions marked by channel infills. The longitudinal bathymetric profile along the modern canyon course shows a knickpoint of ~300 m relief at this frontal fold, indicating that the rate of fold uplift is greater than that of canyon incision. The age for the initial thrusting of this fontal fold is around 240 kyr ago, as estimated by using the maximum thickness of growth strata of this fold divided by the sedimentation rate obtained from a nearby giant piston core. Bottom simulating reflector (BSR) on seismic sections indicates the base of gas hydrate stability zone. Beneath the frontal fold, there is a widespread occurrence of BSRs, suggesting the highly probable existence of substantial quantities of gas hydrates. A seismic flat spot and a few push-down reflectors below BSR are found lying beneath the anticlinal axis with bathymetric four-way dip closure. The flat spot, cutting across a series of dipping reflections beneath BSR, may indicate the contact between free gas and its underlying formation water. The push-down reflectors beneath BSRs are interpreted to result from abundant free gas hosted beneath the gas hydrate stability zone. The multiple paleo-canyon infills seen along and beneath the frontal fold and above BSRs may provide thick porous sands to host gas hydrates in the frontal fold.     

Multi-frequency seismic study of gas hydrate-bearing sediments in Lake Baikal,Siberia

In this paper we present and discuss the frequency-dependent behaviour of the acoustic characteristics of methane hydrate-bearing sediments in Lake Baikal, Siberia. Five different types of seismic sources (airgun-array, two types of single airguns, watergun and sparker) are used, encompassing a frequency bandwidth from 10 up to 1000 Hz. On low-frequency airgun-array data, the base of the hydrate stability zone (HSZ) is observed as a high-amplitude bottom-simulating reflection (BSR) with reversed polarity. The amplitude and continuity of the BSR decrease or even disappear on medium- to high-frequency data, a feature explained in terms of vertical and horizontal resolution. The increasing reflection amplitude of the BSR with increasing offset, the calculated reflection coefficient of the BSR and the occurrence of enhanced reflections below the BSR suggest the presence of free gas below the HSZ. The observation of some enhanced reflections extending above the BSR may be interpreted as an indication for free gas co-existing with hydrates within the HSZ. Amplitude blanking above the BSR is highly variable while the BSR itself appears to act as a low-pass frequency filter for medium- to high-frequency data.New single-channel airgun profiles provide the first seismic information across the Baikal Drilling Project (BDP-97) deep drilling site, at which hydrate-bearing sediments were retrieved at about 200 m above the base of the local HSZ. At the drilling site there are no seismic characteristics indicative of the presence of hydrates. Combination of the drilling and seismic information has allowed us to make a rough estimation of the volume of hydrates and carbon stored in the sediments of Lake Baikal, which lead us to conclude that the Lake Baikal gas hydrate reservoirs do not form a prospective energy resource.     

Seismic evidence and formation mechanism of gas hydrates in the Zhongjiannan Basin,Western margin of the South China Sea

An analysis of 3D seismic data from the Zhongjiannan Basin in the western margin of the South China Sea (SCS) reveals seismic evidence of gas hydrates and associated gases, including pockmarks, a bottom simulating reflector (BSR), enhanced reflection (ER), reverse polarity reflection (RPR), and a dim amplitude zone (DAZ). The BSR mainly surrounds Zhongjian Island, covering an area of 350 km² in this 3D survey area. The BSR area and pockmark area do not match each other; where there is a pockmark developed, there is no BSR. The gas hydrate layer builds upward from the base of the stability zone with a thickness of less than 100 m. A mature pockmark usually consists of an outside trough, a middle ridge, and one or more central pits, with a diameter of several kilometers and a depth of several hundreds of meters. The process of pockmark creation entails methane consumption. Dense faults in the study area efficiently transport fluid from large depths to the shallow layer, supporting the formation of gas hydrate and ultimately the pockmark.     

Seismic attribute study for gas hydrates in the Andaman Offshore India

Seismic data from the Andaman offshore region has been examined to investigate for the presence of gas hydrates. The seismic data displays reflection characteristics such as blanking, enhanced reflection patterns, shadows in instantaneous frequency, and increase in amplitude with the offset, which are indicative of gas hydrates and underlying free gas. A prominent bottom-simulating reflection, BSR, coupled with reverse polarity is observed around 650–700 ms. Seismic attributes such as the reflection strength and instantaneous frequency are computed along this reflector in order to probe for the presence of gas hydrates or free gas in this region. The reflection plot shows a strong reflector paralleling the seafloor. In addition, attenuation of the high frequency signal is noticed, indicating the presence of free gas below the BSR.     

Reservoir characterization of gas hydrate in the Northwestern part of the Ulleung Basin,East Sea

An analysis of 3D seismic data from the northwestern part of the Ulleung Basin, East Sea, revealed that the gas hydrate stability zone (GHSZ) consists of five seismic units separated by regional reflectors. An anticline is present that documents activity of many faults. The seismic indicators of gas hydrate occurrence included bottom simulating reflector (BSR) and acoustic blanking in the gas hydrate occurrence zone (GHOZ). By the analysis of the seismic characteristics and the gradient of the sedimentary strata, the GHOZ was divided into four classes: (1) dipping strata upon strong BSR, (2) dipping strata below strong BSR, (3) parallel strata with acoustic blanking, and (4) parallel strata below weak BSR. Seismic attributes such as reflection strength and instantaneous frequency were computed along the GHOZ. Low reflection strength and high instantaneous frequency were identified above the BSR, indicating the occurrence of gas hydrate. A remarkably high reflection strength and low instantaneous frequency indicated the presence of free gas below the BSR. Considering the distribution of the gas hydrate and free gas, two gas migration processes are suggested: (1) stratigraphic migration through the dipping, permeable strata and (2) structural migration from below the GHSZ along faults.     

The enigma of double BSRs: indicators for changes in the hydrate stability field?

 A classical bottom simulating reflector (BSR) and a presently unknown double BSR pattern are detectable in reflection seismic profiles from the Storegga Slide area west of Norway. Pressure and temperature modeling schemes lead to the assumption that the strong BSR marks the base of a hydrate stability zone with a typical methane gas composition of 99%. The upper double BSR may mark the top of gas hydrates and the lower double BSR may represent a relict of former changes of the hydrate stability field from glacial to interglacial times or the base of gas hydrates with a gas composition including heavier hydrocarbons.     

Seismic Insights into Bottom Simulating Reflection (BSR) in the Krishna-Godavari Basin,Eastern Margin of India

The multichannel seismic data along one long-offset survey line from Krishna-Godavari (K-G) basin in the eastern margin of India were analyzed to define the seismic character of the gas hydrate/free gas bearing sediments. The discontinuous nature of bottom simulating reflection (BSR) was carefully examined. The presence of active faults and possible upward fluid circulation explain the discontinuous nature and low amplitude of the BSR. The study reveals free gas below gas hydrates, which is also indicated by enhancement of seismic amplitudes with offsets from BSR. These findings were characterized by computing seismic attributes such as the reflection strength and instantaneous frequency along the line. Geothermal gradients were computed for 18°C and 20°C temperature at the depth of BSR to understand the geothermal anomaly that can explain the dispersed nature of BSR. The estimated geothermal gradient shows an increase from 32°C/km in the slope region to 41°C/km in the deeper part, where free gas is present. The ray-based travel time inversion of identifiable reflected phases was also carried out along the line. The result of velocity tomography delineates the high-velocity (1.85–2.0 km/s) gas hydrate bearing sediments and low-velocity (1.45–1.5 km/s) free gas bearing sediments across the BSR.     

Seismic character of gas hydrates on the Southeastern U.S. continental margin

Gas hydrates are stable at relatively low temperature and high pressure conditions; thus large amounts of hydrates can exist in sediments within the upper several hundred meters below the sea floor. The existence of gas hydrates has been recognized and mapped mostly on the basis of high amplitude Bottom Simulating Reflections (BSRs) which indicate only that an acoustic contrast exists at the lower boundary of the region of gas hydrate stability. Other factors such as amplitude blanking and change in reflection characteristics in sediments where a BSR would be expected, which have not been investigated in detail, are also associated with hydrated sediments and potentially disclose more information about the nature of hydratecemented sediments and the amount of hydrate present.Our research effort has focused on a detailed analysis of multichannel seismic profiles in terms of reflection character, inferred distribution of free gas underneath the BSR, estimation of elastic parameters, and spatial variation of blanking. This study indicates that continuous-looking BSRs in seismic profiles are highly segmented in detail and that the free gas underneath the hydrated sediment probably occurs as patches of gas-filled sediment having variable thickness. We also present an elastic model for various types of sediments based on seismic inversion results. The BSR from sediments of high ratio of shear to compressional velocity, estimated as about 0.52, encased in sediments whose ratios are less than 0.35 is consistent with the interpretation of gasfilled sediments underneath hydrated sediments. This model contrasts with recent results in which the BSR is explained by increased concentrations of hydrate near the base of the hydrate stability field and no underlying free gas is required.     

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.     

Variability of the bottom-simulating reflector (BSR) and its association with tectonic structures in the Chilean margin between Arauco Gulf (37°S) and Valdivia (40°S)

Multichannel seismic reflection data recorded between Arauco Gulf (37°S) and Valdivia (40°S), on the Chilean continental margin, were processed and modeled to obtain seismic images and sub-surface models, in order to characterize the variability of the bottom-simulating reflector (BSR), which is a geophysical marker for the presence of gas hydrates. The BSR is discontinuous and interrupted by submarine valleys, canyons, as well as by faults or fractures. The BSR occurrence is more common south of Mocha Island due to moderate slopes and greater organic matter contribution by rivers in that area. Tectonic uplift and structural instability change the stability gas hydrate zone and consequently the BSR position, creating in some cases missing or double BSRs. Our modeling supports the presence of gas hydrate above the BSR and free gas below it. Higher BSR amplitudes support higher hydrate or free gas concentrations. In the study area, gas hydrate concentration is low (an average of 3.5%) suggesting disseminated gas hydrate distribution within the sediments. Also higher BSR amplitudes are associated with thrust faults in the accretionary prism, which serve as conduits for gas flow from deeper levels. This extra gas supply produces a wider thickness of gas hydrates or free gas.     

天然气水合物BSR的识别与地震勘探频率

地震勘探是调查天然气水合物广泛使用而有效的方法,而BSR是水合物赋存的主要标志.通过对实际调查资料的分析对比,结合国外的调查研究成果,探讨了地震勘探频率在BSR识别中的影响和作用,提出了在我国海洋天然气水合物的地震调查中有利于BSR识别的合适的频率范围.     

Amplitude-versus-offset modeling of the bottom simulating reflection associated with submarine gas hydrates

Velocity analysis of multi-channel seismic (MCS) data and amplitude-versus-offset (AVO) modeling provides an efficient way of identifying gas hydrate and free gas, and therefore the nature of the bottom-simulating reflector (BSR). Additionally, AVO modeling also yields estimates of the hydrate concentration and free gas saturation across the BSR in terms of velocity distribution. In the present study, we apply directivity correction in order to accentuate the AVO behavior. Modeling for AVO pattern of the observed BSR over the Kerala–Konkan Offshore Basin may provide the probable velocity distribution across the BSR and thereby infer whether hydrate or hydrate/free gas model governs the AVO observations. Initial results indicate the possible presence of free gas underlying the gas hydrates-saturated sediments in this region.     

Application of AVO attribute inversion technology to gas hydrate identification in the Shenhu Area,South China Sea

AVO (Amplitude Versus Offset) is a seismic exploration technology applied to recognize lithology and detect oil and gas through analyzing the feature of amplitude variation versus offset. Gas hydrate and free gas can cause obvious AVO anomaly. To find geophysical evidence of gas hydrate and free gas in Shenhu Area, South China Sea, AVO attribute inversion method is applied. By using the method, the multiple seismic attribute profiles and AVO intercept versus gradient (I-G) cross plot are obtained. Bottom-simulating reflector (BSR) is observed beneath the seafloor, and the AVO abnormal responses reveal various seismic indicators of gas hydrate and free gas. The final AVO analysis results indicate the existence of gas hydrate and free gas in the upper and lower layers of BSR in the study area.     

Qualitative assessment of gas hydrate and gas concentrations from the AVO characteristics of the BSR in the Ulleung Basin,East Sea (Japan Sea)

The presence of gas hydrate in the Ulleung Basin, East Sea (Japan Sea), inferred by various seismic indicators, including the widespread bottom-simulating reflector (BSR), has been confirmed by coring and drilling. We applied the standard AVO technique to the BSRs in turbidite/hemipelagic sediments crosscutting the dipping beds and those in debris-flow deposits to qualitatively assess the gas hydrate and gas concentrations. These BSRs are not likely to be affected by thin-bed tuning which can significantly alter the AVO response of the BSR. The BSRs crosscutting the dipping beds in turbidite/hemipelagic sediments are of low-seismic amplitude and characterized by a small positive gradient, indicating a decrease in Poisson’s ratio in the gas-hydrate stability zone (GHSZ), which, in turn, suggests the presence of gas hydrate. The BSRs in debris-flow deposits are characterized by a negative gradient, indicating decreased Poisson’s ratio below the GHSZ, which is likely due to a few percent or greater gas saturations. The increase in the steepness of the AVO gradient and the magnitude of the intercept of the BSRs in debris-flow deposits with increasing seismic amplitude of the BSRs is probably due to an increase in gas saturations, as predicted by AVO model studies based on rock physics. The reflection strength of the BSRs in debris-flow deposits, therefore, can be a qualitative measure of gas saturations below the GHSZ.     

Qualitative assessment of gas hydrate and gas concentrations from the AVO characteristics of the BSR in the Ulleung Basin, East Sea (Japan Sea)

The presence of gas hydrate in the Ulleung Basin, East Sea (Japan Sea), inferred by various seismic indicators, including the widespread bottom-simulating reflector (BSR), has been confirmed by coring and drilling. We applied the standard AVO technique to the BSRs in turbidite/hemipelagic sediments crosscutting the dipping beds and those in debris-flow deposits to qualitatively assess the gas hydrate and gas concentrations. These BSRs are not likely to be affected by thin-bed tuning which can significantly alter the AVO response of the BSR. The BSRs crosscutting the dipping beds in turbidite/hemipelagic sediments are of low-seismic amplitude and characterized by a small positive gradient, indicating a decrease in Poisson’s ratio in the gas-hydrate stability zone (GHSZ), which, in turn, suggests the presence of gas hydrate. The BSRs in debris-flow deposits are characterized by a negative gradient, indicating decreased Poisson’s ratio below the GHSZ, which is likely due to a few percent or greater gas saturations. The increase in the steepness of the AVO gradient and the magnitude of the intercept of the BSRs in debris-flow deposits with increasing seismic amplitude of the BSRs is probably due to an increase in gas saturations, as predicted by AVO model studies based on rock physics. The reflection strength of the BSRs in debris-flow deposits, therefore, can be a qualitative measure of gas saturations below the GHSZ.     

Source and accumulation of gas hydrate in the northern margin of the South China Sea

The northern South China Sea (NSCS) experienced continuous evolution from an active continental margin in the late Mesozoic to a stable passive continental margin in the Cenozoic. It is generally believed that the basins in the NSCS evolved as a result of Paleocene–Oligocene crustal extension and associated rifting processes. This type of sedimentary environment provides a highly favourable prerequisite for formation of large-scale oil- and gas–fields as well as gas hydrate accumulation. Based on numerous collected data, combined with the tectonic and sedimentary evolution, a preliminary summary is that primitive coal-derived gas and reworked deep gas provided an ample gas source for thermogenic gas hydrate, but the gas source in the superficial layers is derived from humic genesis. In recent years, the exploration and development of the NSCS oil, gas and gas hydrate region has provided a basis for further study. A number of 2D and 3D seismic profiles, the synthetic comparison among bottom simulating reflector (BSR) coverage characteristics, the oil-gas area, the gas maturity and the favourable hydrate-related active structural zones have provided opportunities to study more closely the accumulation and distribution of gas hydrate. The BSR has a high amplitude, with high amplitude reflections below it, which is associated with gas chimneys and pockmarks. The high amplitude reflections immediately beneath the BSR are interpreted to indicate the presence of free gas and gas hydrate. The geological and geochemical data reveal that the Cenozoic northern margin of the NSCS has developed coal-derived gas which forms an abundant supply of thermogenic gas hydrate. Deep-seated faults and active tectonic structures facilitate the gas migration and release. The thermogenic gas hydrate and biogenic gas are located at different depths, have a different gas source genesis and should be separately exploited. Based on the proven gas hydrate distribution zone, we have encircled and predicted the potential hydrate zones. Finally, we propose a simple model for the gas hydrate accumulation system in the NSCS Basin.