National Gas Hydrate Program Research Articles

Characterising the effect of gas hydrate in pore pressure and stress using NGHP-02 downhole logs in Krishna Godavari offshore Basin, India

The present study demonstrates the impact of gas hydrate concentration and its morphology on pore pressure and sediment stress. We analyse the well-logging data of the 2nd Indian National Gas Hydrate Program (NGHP-02) in Area B of the Krishna Godavari Basin. Overpressure zones are delineated by observing porosity with depth and log inferred compaction trends. Eaton, Bower, Miller and porosity-based models are used for predicting pore pressure. The predicted pore pressures are ranging between ∼3 and 4.5 MPa, and ∼2.8 and 4.4 MPa respectively at Hole 02–22A and 02–23A, where the estimated fracture pressures and horizontal stresses are nearly equal to vertical stress at both holes. Within the gas hydrate stability zone (GHSZ), overpressure is not found; rather, our findings show that pore pressure falls with increasing gas hydrates saturation at both holes. This could be due to a decrease in water concentration during hydrate formation. Consequently, the effective stress of sediments increases, causing gas hydrate-filled vertical fractures to heave apart in the minimum horizontal stress directions. We observe overpressure below the gas hydrate layers at two holes, which is due to the load of the overlying sediments, low permeability created by the high concentration of gas hydrate at the bottom of GHSZ (about 40% and 80% at Hole 02–22A and 02–23A, respectively), and may be the microbial gas being produced at a faster rate than they are being expelled due to the low permeability. It is noticed that the strike directions of fractures vary with depth at both holes, which shows the anisotropic effects of stress conditions on gas hydrate-filled fractures during the time of formation.

Machine-learning application to assess occurrence and saturations of methane hydrate in marine deposits offshore India

Artificial neural networks (ANN) were used to assess methane hydrate occurrence and saturation in marine sediments offshore India. The ANN analysis classifies the gas hydrate occurrence into three types: methane hydrate in pore space, methane hydrate in fractures, or no methane hydrate. Further, predicted saturation characterizes the volume of gas hydrate with respect to the available void volume. Log data collected at six wells, which were drilled during the India National Gas Hydrate Program Expedition 02 (NGHP-02), provided a combination of well-log measurements that were used as input for machine-learning (ML) models. Well-log measurements included density, porosity, electrical resistivity, natural gamma radiation, and acoustic wave velocity. Combinations of well logs used in the ML models provide good overall balanced accuracy (0.79 to 0.86) for the prediction of the gas hydrate occurrence and good accuracy (0.68 to 0.92) for methane hydrate saturation prediction in the marine accumulations against reference data. The accuracy scores indicate that the ML models can successfully predict reservoir characteristics for marine methane hydrate deposits. The results indicate that the ML models can either augment physics-driven methods for assessing the occurrence and saturation of methane hydrate deposits or serve as an independent predictive tool for those characteristics.

Semi-automated Rock Layer Recognition from Borehole Log Data Using Combined Wavelet and Fourier Transform: A Case Study in the KG basin, India

Abstract The accurate identification of rock layer interfaces from wireline logs is one of the crucial issue of geological interpretation as we face difficulties in traditional interpretations due to the presence of various kind of high and low frequencies in the log data. Determination of rock boundaries in shallow marine sediments is very important to understand several issues related to the exploration of conventional and non-conventional hydrocarbons. Presently a combined wavelet and Fourier transform-based approach was demonstrated to discriminate the rock layers from wireline log data in Krishna-Godavari basin, which were acquired during the India National Gas Hydrate Program Expedition 02 (NGHP-Exp.-02). Initially, we selected the suitable mother wavelets and optimum level of decomposition to perform the wavelet transform by computing the error between the original log signal(s) and approximation coefficients (cA’s) of the decomposed signal(s) at several decomposition levels. Subsequently, discrete wavelet transform (DWT) of a particular decomposition level was applied to the log signals to obtain the detail coefficients (cD’s), which contain the high-frequency components of the log signals. Further, these high-frequency components of the log signals were utilized in frequency spectrum-based filtering approaches. This approach was adapted to yield information about the abrupt changes occurring in the log signals across the studied depth interval. The combined application of Fourier and wavelet transforms to shallow marine log data demonstrates that this approach is an efficient tool for delineating rock layer boundaries.

Gas hydrate saturation from NGHP 02 LWD data in the Mahanadi Basin

During the Indian National Gas Hydrate Program (NGHP) Expedition 02, Logging-while-drilling (LWD) logs were acquired at three sites (NGHP-02-11, NGHP-02-12, and NGHP-02-13) across the Mahanadi Basin in area A. We applied rock physics theory to available sonic velocity logs to know the distribution of gas hydrate at site NGHP-02-11 and NGHP-02-13. Rock physics modeling using sonic velocity at well location shows that gas hydrate is distributed mainly within the depth intervals of 150–265 m and 100–215 mbsf at site NGHP-02-11 and NGHP-02-13, respectively, with an average saturation of about 4% of the pore space and the maximum concentration of about 40% of the pore space at 250 m depth at site NGHP-02-11, and at site NGHP-02-13 an average saturation of about 2% of the pore space and the maximum concentration of about 20% of the pore space at 246 m depth, as gas hydrate is distributed mainly within 100–246 mbsf at this site. Saturation of gas hydrate estimated from the electrical resistivity method using density derived porosity and electrical resistivity logs from Archie's empirical formula shows high saturation compared to that from the sonic log. However, estimates of hydrate saturation based on sonic P-wave velocity may differ significantly from that based on resistivity, because gas and hydrate have higher resistivity than conductive pore fluid and sonic P-wave velocity shows strong effect on gas hydrate as a small amount of gas reduces the velocity significantly while increasing velocity due to the presence of hydrate. At site NGHP-02-11, gas hydrate saturation is in the range of 15%–30%, in two zones between 150-180 and 245–265 mbsf. Site NGHP-02-012 shows a gas hydrate saturation of 20%–30% in the zone between 100 and 207 mbsf. Site NGHP-02-13 shows a gas hydrate saturation up to 30% in the zone between 215 and 246 mbsf. Combined observations from rock physics modeling and Archie’s approximation show the gas hydrate concentrations are relatively low (<4% of the pore space) at the sites of the Mahanadi Basin in the turbidite channel system.

Characterization of a complex sand-rich gas hydrate reservoir system in the Indian marine continental margin with downhole log and seismic data

Logging-while-drilling (LWD) and coring data were acquired in Areas A, B, C and E during the National Gas Hydrate Program Expedition 02 (NGHP-02). At Sites NGHP-02-16, −17, −20, −23 and −24 of Area B in the Indian Krishna-Godavari Basin, the gas hydrate-bearing sand-rich reservoirs directly in contact with an underlying water-saturated sediment were revealed near the seismic-inferred bottom simulating reflector (BSR) along a large anticlinal structure. Lithofacies analysis of cores show that the lower gas hydrate-bearing reservoir section at these sites is an interbedded unit with thin fine-sand layers alternating with clay-rich layers. We reanalyzed the LWD data including resistivity, compressional-wave and shear-wave acoustic velocities to assess and compare the spatial variation of sand-hosted gas hydrate and possible related free-gas accumulations. The presence of free gas-bearing sediments in direct contact with gas hydrate-bearing sediments was confirmed at Sites NGHP-02-20 and −24 along the northeastern flank of the anticline. The sand-rich reservoirs associated with the Area B anticlinal feature reveal different pore-filling constituents and variable boundary conditions, including gas hydrate/water contacts and gas hydrate/free gas contacts. We infer that the sedimentary depositional history at Sites NGHP-02-20 and −24, may have resulted in an upward shift of the base of gas hydrate stability zone, that could have led to in situ gas hydrate dissociation and the accumulation of free gas and coexistence interval below the log-inferred gas hydrate-bearing reservoir sections. This study demonstrates that reservoir heterogeneity and fluid migration were major factors affecting the distribution of gas hydrate and free gas in the identified sand-rich reservoirs.

Surface Drilling Parameters and Drilling Optimization Techniques: Are They Useful Tools in Gas Hydrate Detection?

This paper examines the application of surface drilling parameters and drilling optimization techniques, such as mechanical specific energy (MSE) and equivalent strength (EST), in detecting hydrate gas-bearing sediments during drilling operations. Gas hydrates are usually detected from 3D seismic imaging and later confirmed with Measurement-While-Drilling (MWD)/Logging-While-Drilling (LWD) data and collected core samples. Here, we describe an analysis of the time-based surface drilling parameters recorded from two wells drilled during the International Ocean Discovery Program (IODP) Expedition 372A offshore of New Zealand and the Indian National Gas Hydrate Program Expedition 02 (NGHP-02) offshore of India. The investigation revealed that drilling parameters, as well as MSE/EST methods, can and should be used to monitor and optimize the drilling process and to detect lithological/tectonic (fractures, fault zones, rock hardness, etc.) changes in the drilled substrata and signs of the dynamic changes in the downhole environment (tool vibration, washouts, pack-offs, etc.). However, surface drilling parameters with MSE models cannot explicitly determine the hydrate gas-bearing sediments. This qualitative analysis of whether the gas-bearing sediments consist of hydrates can only be accomplished with the use of the MWD/LWD suite, preferably located as close as possible to the bit.

Optimal design of the field hydrate production test in the offshore India: Insights from the vertically heterogeneous hydrate reservoir model

Recently, the hydrate-bearing sand reservoirs were founded with high saturation in the offshore India by the National Gas Hydrate Program Expedition 02 (NGHP-02). Based on the detailed geological data at the NGHP-02-16 site, a more realistic reservoir model is constructed to analyze the unique multiphase flow and hydrate phase change behaviors from the geologically descriptive hydrate reservoirs. The main goals of this work are to optimize the adjustable production parameters and to predict the long-term gas production performance for guiding the upcoming field hydrate production test in India. The results indicate that gas production rates increase, but water production rates decrease slowly with time from the heterogeneous hydrate-reservoir. This is completely different from the previous prediction of multiphase flow behavior based on homogeneous hydrate-reservoir models. A non-uniform hydrate dissociation front is formed in the hydrate-reservoir due to the reservoir heterogeneity promote the horizontal dissociation of solid-hydrate within the high permeability sand layers. The maximum distance of free gas zone expands to outside more than 500 m from the production well. This indicates that the field fault system should be carefully concerned during the long-term production. Sensitivity analysis results suggest that the complete perforation in the hydrate-bearing layer and the underlying water saturated layer is beneficial for improving the hydrate production performance. Over the 5 years depressurization, a total of 2.08 × 107 ST m3 methane gas can be produced using the recommended well placement. The average gas-to-water ratio and energy ratio reach 10.22 ST m3 CH4/m3 H2O and 30.14, respectively, showing a high gas production potential from this hydrate site.

Determining the habitat depth of the planktic foraminifera Dentoglobigerina altispira in the eastern Arabian Sea during the middle Miocene

Tests of planktic foraminifera are important tracers for reconstructing past oceanic environments. It is essential to have accurate information of the habitat depth of the planktic foraminiferal species whose isotopic and elemental signatures are being used to infer past climatic and oceanographic conditions. The planktic foraminifera species Dentoglobigerina altispira is a very common component of Miocene and Pliocene tropical to warm subtropical planktic foraminiferal assemblages but different findings exist about the habitat depth of this species from shallow to deep water. For that reason, comparative stable carbon and oxygen isotope analyses of four coexisting planktic foraminifera species from National Gas Hydrate Program Expedition Site 01 (Site NGHP-01-01A) were performed in order to reconstruct the habitat depth of D. altispira in the eastern Arabian Sea for the middle Miocene (~16–12.5 Ma). The D. altispira oxygen isotope values are on average higher than those of the shallow water reference species Trilobatus sacculifer and Globigerinoides obliquus, but still significantly lower than those for the deeper dwelling species Dentoglobigerina venezuelana. This suggests that D. altispira inhabited a deeper habitat depth than the estimated habitat depth for modern T. sacculifer in the eastern Indian Ocean (~50 m). Paired D. altispira and T. sacculifer Mg/Ca data reveal lower Mg/Ca and thus colder temperatures for D. altispira than for T. sacculifer, also suggesting that D. altispira calcified at a greater water depth than T. sacculifer. The data of this study suggest that D. altispira inhabited the lower mixed layer and/or upper thermocline in the eastern Arabian Sea during the middle Miocene.

Machine learning in the classification of lithology using downhole NMR data of the NGHP-02 expedition in the Krishna-Godavari offshore Basin, India

We apply machine learning techniques to classify the lithology using downhole nuclear magnetic resonance (NMR) data of gas hydrate reservoirs in the Krishna-Godavari (KG) Basin, eastern Indian offshore. We choose two sites of the second expedition of the Indian National Gas Hydrate Program (NGHP-02) for our study, one with clay dominated site NGHP-02-01 from Area E and another with silt/sand dominated site NGHP-02-05 from Area C. We obtain five distinct classes in NMR signals at two sites using unsupervised techniques such as Davies-Bouldin index, silhouette, self-organizing map, k-means clustering and emergent self-organizing map. Obtained classes are interpreted using the NMR T2 distribution pattern with time in terms of clay, silt, sand and pebble/gravels with depth and reconciled with gamma-ray, density, porosity, resistivity, velocity and photoelectric factor. Our results demonstrate the lithology as silty-clay with minor sand at Hole NGHP-02-01A, and silt, sand with less clay and minor pebble/gravel at Hole NGHP-02-05A. The presence of gas hydrate is identified by observing an overall low amplitude T2 distribution curve, which is due to the reduction of pore space by solid gas hydrate. The maximum concentration of gas hydrate of about 43% of pore space estimated using NMR- and density-porosity are distributed mainly in silty-clay at Hole 01A and of about 75% in clayey-silt and silty-sand at Hole 05A. Low amplitude NMR T2 signals observed in clay-dominated sediments at higher relaxation time may be due to the fracture-filled gas hydrates at both holes. Interestingly, high porosity (90–95%), high permeability (4 mD), low bulk-density (1.2 gm/cm3) and low photoelectric factor (0.8 b/electron) are found in sandy layers within the gas hydrate-stability zone at Hole 05A, which possibly comprising buried channels. Comparison of the lithology obtained from other log data with that from NMR data and core samples illustrates that the NMR data are most effective for characterization of reservoirs, except for a high concentration of gas hydrate in a sandy/coarse-grained reservoir, which needs further studies.

Assessment of gas hydrate using prestack seismic inversion in the Mahanadi Basin, offshore eastern India

The Indian National Gas Hydrate Program Expedition-01 in 2006 has discovered gas hydrate in the Mahanadi offshore basin along the eastern Indian margin. However, well-log analysis, pressure core measurements, and infrared anomalies reveal that gas hydrates are distributed as disseminated within the fine-grained sediment, unlike massive gas hydrate deposits in the Krishna-Godavari Basin. The 2D multichannel seismic section, which crosses holes NGHP-01-9A and 19B located at approximately 24 km apart, indicates a continuous bottom-simulating reflector (BSR) along it. We aim to investigate the prospect of gas hydrate accumulation in this area by integrating well-log analysis and seismic methods with rock-physics modeling. First, we estimate gas hydrate saturation at these two holes from the observed impedance using the three-phase Biot-type equation. Then, we establish a linear relationship between the gas hydrate saturation and the impedance contrast with respect to the water-saturated sediment. Using this established relation and impedance obtained from prestack inversion of seismic data, we produce a 2D gas hydrate-distribution image over the entire seismic section. The gas hydrate saturation estimated from resistivity and sonic data at well locations varies within 0%–15%, which agrees well with the available pressure core measurements at hole 19. However, the 2D map of gas hydrate distribution obtained from our method indicates that the maximum gas hydrate saturation is approximately 40% just above the BSR between the common-depth points of 1450 and 2850. The presence of gas-charged sediments below the BSR is one of the reasons for the strong BSR observed in the seismic section, which is depicted as low impedance in the inverted impedance section. Closed sedimentary structures above the BSR are probably obstructing the movements of free-gas upslope, for which we do not see the presence of gas hydrate throughout the seismic section above the BSR.

Amount of gas hydrate estimated from rock physics analysis based on morphology and intrinsic anisotropy in area B, Krishna Godavari offshore basin, expedition NGHP-02

India's National Gas Hydrate Program (NGHP) Expedition-02 was conducted in 2015 to investigate the prospective locations of sand-rich gas hydrate reservoirs along the eastern offshore, India. We select Area-B among four areas A, B, C and E of expedition NGHP-02 as it was completely unexplored during NGHP-01. Logging while Drilling (LWD) data were collected at twelve sites in area B in Krishna Godavari (KG) offshore basin, where a huge amount of gas hydrate (up to 90% of the pore volume) is present as pore-fill, fracture-fill and both pore- and fracture-filled (mixed) morphology. The resistivity log generally overestimates gas hydrate saturation because of the increase of resistivity of fresh water near gas hydrate accumulation and also in the presence of fractured-filled gas hydrate. Though gas hydrate saturation obtained from pressure cores are assumed to be ground truth but they are discrete. Hence, the application of an advanced rock physics model is required to obtain continuous quantification of gas hydrate throughout the depth considering the gas hydrate morphology and microstructure of sediment. We apply a combination of self-consistent approximation (SCA) and differential effective medium (DEM) theories to model biconnected marine gas hydrate reservoir that can deal with inherent anisotropy due to the presence of clay platelets and pore-filling, fracture-filling or mixture of both morphologies of gas hydrate. First, we determine the distribution pattern of clay platelets from the observed background velocity (water-saturated sediment) and then estimate the gas hydrate saturation at four holes of area B, where pressure core measurements are available. We predict maximum gas hydrate saturation as 50–90% of the pore space above the BSR in this area. Comparison with other rock physics model shows that prediction from the combined SCA-DEM theory matches better with the pressure core measurements.

Timing and cause of the disappearance of some elongated taxa in the Indian Ocean: Study from NGHP Hole 17A

Disappearance of many elongated benthic foraminifera across the Mid-Pleistocene Transition (MPT) is popularly known as ‘Stilostomella extinction’. This event is reported from several sites across the world and it is believed that most of the elongated benthic families disappeared within 0.76–0.5 Ma. It is assumed that the presence of cold bottom water triggered their disappearance. The present study is pursued on sediment cores collected from National Gas Hydrate Program Hole 17A, Andaman Sea, Indian Ocean to examine the precise timing and cause of the disappearance of elongated taxa as well as their ecological preference. Our study reveals that elongated species prefer warm and low to intermediate food condition. This study documents the disappearance of 5 major species and 13 minor species from 0.6 to 0.2 Ma which is younger than the earlier estimates. We assume that change in 41–100 kyr climatic cycle along the MPT was responsible for a major decline in their number. However, amplified glacial–interglacial cycles along Mid-Brunhes (0.5–0.2 Ma) finally played a major role in their disappearance. Both of these consecutive events extended intensified cooling and dropped the deep-ocean temperature by which the shallow infaunal elongated species disappeared and deep-infaunal species survived.

Gas hydrate petroleum systems: What constitutes the “seal”?

The gas hydrate petroleum system (GHPS) approach, which has been used to characterize gas hydrates in nature, uses three distinct components: a methane source, a methane migration pathway, and a reservoir that not only contains gas hydrate, but also acts as a seal to prevent methane loss. Unlike GHPS, a traditional petroleum system (PS) approach further distinguishes between the reservoir, a unit with generally coarser sediment grains, and a separate overlying seal unit with generally finer sediment grains. Adopting this traditional PS distinction in the GHPS approach facilitates assessments of reservoir growth and production potential. The significance of the seal for the formation of a gas hydrate reservoir as well as for efficiency in methane extraction from the reservoir as an energy resource is evident in findings from recent offshore field expeditions, such as India’s second National Gas Hydrate Program expedition (NGHP-02). In regard to gas hydrate-bearing reservoir formations, the NGHP-02 gas chemistry data indicate a primarily microbial methane source. Fine-grained seal sediment in contact with coarser grained reservoir sediment can facilitate that microbial methane production. Logging-while-drilling and sediment core data also indicate that the overlying fine-grained seal sediment is less permeable than the underlying, highly gas hydrate-saturated reservoir sediment. The overlying seal’s capacity to act as a low-permeability boundary is important not only for preventing methane migration out of the reservoir over time, but also for preventing water invasion into the reservoir during methane extraction from the reservoir. Ultimately, the presence of an overlying, fine-grained, low-permeability “seal” influences how gas hydrate initially forms in a coarse-grained reservoir and dictates how efficiently methane can be extracted as an energy resource from the gas hydrate reservoir via depressurization.

Predicting lithology using neural networks from downhole data of a gas hydrate reservoir in the Krishna–Godavari basin, eastern Indian offshore

SUMMARYWe use the unsupervised and supervised neural network methods together to predict lithology of a gas hydrate reservoir from downhole data in the Krishna–Godavari (KG) offshore basin, India. In this study, we successfully identify the host litho-units of gas hydrate and show its effects in the identification of lithology using neural network techniques, which is not reported earlier. We use well log data acquired at three holes (10A, 03A and 04A) in 2006 during the first expedition of the Indian National Gas Hydrate Program (NGHP-01). Five different logging while drilling data (e.g. density, neutron porosity, gamma ray, resistivity and sonic) are considered for the mapping of lithology and gas hydrate. In the presence of gas hydrate, the resistivity and sonic velocity of the host sediments increase significantly, whereas density, neutron porosity and gamma ray are negligibly affected. Therefore, we calculate resistivity and sonic velocity for water-saturated sediment (without gas hydrate) theoretically to remove the effects of gas hydrate. At first, we apply the seven unsupervised classification methods (i.e. elbow, dendrogram, K-means, 3-D clustering, principal component analysis, Devies–Bouldin index and self-organizing map) to the data with gas hydrate (e.g. observed) and without gas hydrate (i.e. water-saturated/theoretical) to assess the data dimensionality and the number of clusters/litho-units. Each of the unsupervised schemes has its own pros and cons, and may provide different number of cluster/litho-units; sometimes, it is difficult to interpret from only one method. However, all seven methods provide same number of clusters in our study. Then, we apply the supervised classification method (i.e. Bayesian neural networks optimized by hybrid Monte Carlo searching technique) to the training data to refine the defined litho-units and map them with depth. Our approach identifies four types of litho-units and illustrates that the lithology in this area is dominated by clay (∼64 per cent) with some amount of silty clay, silt and minor sand. Gas hydrate is found in clay, silty clay and silt and not in sand. Results also show that, if gas hydrate is not considered as a separate unit, it is distributed as lithology in its hosts (i.e. clay, silty clay and silt). The method is very stable up to ∼15 per cent of random noise added to the data and results are well matched with the analysis of recovered core data. Identified lithologies at three wells correlate very well with seismic section crossing the wells. Very low permeability (&amp;lt;0.1 mD) estimated at three wells also indicates the clay-dominated lithology in our study area.

Bed boundary identification from well log data using Walsh transform technique: A case study from NGHP Expedition-02 in the Krishna–Godavari basin, India

Identification of thin stratigraphic beds is often dealt with using the scale of core data analysis. The downhole wireline log data also provide accurate information on the characteristics of lithological boundaries. In the present study, we have used a Walsh low pass filter and bed boundary detection algorithm to develop an automated bed boundary identification approach. Initially, we ran a Walsh low pass filter against the wireline log responses. Afterwards, a bed boundary detection algorithm has been applied to the Walsh low pass filtered version of the wireline log data. The Walsh domain filter creates a stepped version of the well log data having a constant step width over the entire low pass version of the log signal. As the Walsh function basically represents the binary waveforms, the Walsh low pass filtering can appropriately identify the thin lithological units within a complex succession of sedimentary strata. Further, the proposed approach is an efficient way of resolving and understanding subsurface inhomogeneity from wireline log data. The feasibility of the technique is successfully tested on the downhole wireline log data acquired from the Krishna–Godavari basin during the Expedition 02 of Indian National Gas Hydrate Program (MGHP).

India National Gas Hydrate Program Expedition 02 Summary of Scientific Results: Gas hydrate systems along the eastern continental margin of India

Abstract The primary objectives of the India National Gas Hydrate Program Expedition 02 (NGHP-02) were to obtain new data on the occurrence of gas hydrate systems and to advance the understanding of the controls on the formation of gas hydrate accumulations in the offshore of India. In accordance with the ultimate overall goal of the NGHP effort to assess the energy resource potential of marine gas hydrates in India, particular focus was placed on the exploration and evaluation of gas hydrate occurrences at high saturations in sand-rich systems. NGHP-02 operations were conducted from 3-March-2015 to 28-July-2015 off the eastern coast of India and included logging while drilling (LWD) operations at 25 locations, and coring and wireline logging operations at 10 locations, in the Krishna-Godavari and Mahanadi Basins. The formation of highly concentrated gas hydrate accumulations, which are more suitable for energy extraction, requires the presence of relatively coarse-grained sediments with porosity needed to support the migration and accumulation of gas, and the nucleation of gas hydrate. The results of downhole logging, coring and formation pressure testing operations during NGHP-02 have confirmed the presence of extensive sand-rich depositional systems throughout the deepwater portions of the Krishna-Godavari and Mahanadi Basins. Two areas of the Krishna-Godavari Basin, referred to as Areas B and C, contain substantial gas hydrate accumulations in sand-rich systems and therefore represent ideal candidate sites for future gas hydrate production testing. This summary and technical report includes a comprehensive synthesis of the geologic, geophysical, geochemical, and physical property data acquired during NGHP-02 as it relates to the controls on gas hydrate occurrence, particularly with regards to sand-hosted accumulations. In the Mahanadi Basin, despite the confirmation of extensive reservoir capacity, gas supply at the NGHP-02 sites was insufficient to charge the reservoirs with gas hydrates. In the Krishna-Godavari Basin, extensive reservoir systems were confirmed with sediment grain-sizes ranging from coarse-silts to gravels. These reservoirs range from fully- to partially-filled with gas hydrate. The gas is determined to be from only microbial sources, and in part migrated into the reservoirs from deeper systems. The controls on gas hydrate occurrence are complex and varied; and include substantial reservoir heterogeneity and sufficient permeability throughout the reservoirs and seals that allowed pervasive fluid flow into and through the hydrate-bearing systems. These discoveries are the most significant confirmation of the exploration approach that focuses on direct detection of hydrate reservoirs supported by comprehensive petroleum systems analyses.

Estimation of gas hydrate saturation using isotropic and anisotropic modelling in the Mahanadi basin

A base of gas hydrate stability zone was established after coring and drilling under the National Gas Hydrate Program (NGHP) Expedition-01 in the Mahanadi basin. At two sites, logging-while-drilling log data, and, at one site, wireline log data, were acquired during the NGHP Expedition-01. Gas hydrate reservoirs modelling can be performed in two different ways. One way is isotropic (load bearing) and, on the other hand, anisotropic media (fracture filling with gas hydrate). Here, we have performed anisotropic modelling and estimated gas hydrate saturation using P-wave velocity, assuming an incidence angle of 75\(^{\circ }\) represents the vertical fracture. The estimated gas hydrate saturation at sites NGHP-01-08 and NGHP-01-09, assuming anisotropic media, reduces the estimate by half compared to the saturation estimation by assuming isotropic media. The saturation at site NGHP-01-19 estimated from the isotropic and anisotropic P-wave velocity models are more or less similar except in the zone (175–210 m) just above the bottom simulating reflector depth, and this zone shows similar reduction in saturation as estimated at sites NGHP-01-08 and NGHP-01-09. Observations show that average gas hydrate saturations are relatively low (up to 5% of the pore space). The saturation of a gas hydrate estimated from an isotropic P-wave model varies from 5% to 20%. However, the saturation estimated from the anisotropic P-wave model shows a variation up to 10% of the pore spaces at three sites.

Lithofacies characterization of gas hydrate prospects discovered during the National Gas Hydrate Program expedition 02, offshore Krishna-Godavari Basin, India

Abstract Lithofacies of host sediments associated with gas hydrate reservoirs may have a significant influence on the formation of hydrates. The type of sediments and the depositional processes responsible for the gas hydrate-bearing sediments along the eastern margin of the India have been studied by integrating core data and logging-while-drilling downhole logs. Borehole data correlated with seismic data were used for the interpretation of geological process responsible for gas hydrate formation. The acquired seismic profiles were correlated with the downhole logs by identifying the shape of the well logs curves that indicate finning up and blocky patterns. Turbidite sands and debris flows were also identified from the log and core sediment studies. This study identified two prominent types of depositional systems: (1) turbidity currents and (2) debris flows were mainly responsible for the deposition of gas hydrate-bearing reservoirs in the study area. Based on these studies two well-developed gas hydrate systems have been identified in the Krishna-Godavari Basin and hydrates were determined to be distributed mostly in sand dominated facies. Fine-grained sand in the identified turbidite deposits were shown to be a good hydrate reservoir. The hydrate-bearing sediments in the central part of the Krishna-Godavari Basin (Area C in D6 block) consist of mostly debris flows, which is consistent with scattered gravel- and sand/silt-sized grains in channel-levee systems.

Assessment of gas hydrate accumulations using velocities derived from vertical seismic profiles and acoustic log data in Krishna-Godavari Basin, India

The second Indian National Gas Hydrate Program Expedition (NGHP-02) was executed in 2015 in four areas termed as Areas A, B, C and E. During this expedition, twenty five research sites were drilled and/or cored in Krishna-Godavari (KG) and Mahanadi Basins, eastern Indian offshore. During NGHP-02, zero offset vertical seismic profile (VSP) data were acquired at three sites: Sites NGHP-02-17, −19 and −22 in Area B of the KG Basin. In this study, we focus on the three sites in Area B of the KG Basin, where, zero-offset VSP and downhole acoustic log data are used to assess and characterize the gas hydrate deposits. Zero-offset VSP data are correlated with wireline log, surface seismic and synthetic seismic data to characterize and delineate gas hydrate accumulations in the KG Basin. Low velocities ranging from 1500 to 1650 m/s are observed in the unconsolidated shallow sedimentary section above gas hydrate-bearing units, whereas, very high velocities are observed in the acquired acoustic log and VSP data. In the gas hydrate-bearing sedimentary sections, the VSP derived interval velocities range from 2000 to 3000 m/s in the depth interval between 267 and 287 m below sea floor (mbsf) at Site NGHP-02-17, between 1650 and 1750 m/s in the depth interval between 306 and 366 mbsf at Site NGHP-02-19 and between 1700 and 1800 m/s in the depth interval between 198 and 290 mbsf at Site NGHP-02-22. Gas hydrates are distributed both as pore-filling and fracture-filling at all three sites, however, high concentrations are only observed as pore-filling morphology. We estimate the amount of gas hydrate considering both isotropic (pore-filling) and anisotropic (fracture-filling) acoustic reservoir models. Our estimations of gas hydrate saturation match well with the available pressure core data and suggest that high gas hydrate concentrations (to nearly 85% of pore space) are distributed in a load bearing morphology.

Lithological properties of natural gas hydrate–bearing sediments in pressure-cores recovered from the Krishna–Godavari Basin

Abstract We investigated the sediment properties and gas hydrate saturations in pressure-cores recovered during the National Gas Hydrate Program Expedition 02 (NGHP-02) in Areas B, C, and E of the Krishna–Godavari Basin (K–G Basin) off the eastern margin of India. These investigations included the analyses of grain size, mineralogy, and grain density and employed scanning electron microscopy (SEM). Sediments with median grain sizes below 90 μm were generally observed within most of the pressure-cores, whereas sediments with median grain sizes exceeding 250 μm were only observed in a few pressure-cores. Higher concentrations of gas hydrate were found in sediment layers with median grain sizes of 10–90 μm; the occurrence of gas hydrates was lower in sediment layers with median grain sizes below 10 μm. Three types of gas hydrate-bearing sediments containing above 50 vol.% gas hydrate saturations were observed: (1) sediments with silt grain size ratios above 70 vol.%, (2) sediments with sand grain size ratios above 35 vol.%, and (3) coarser sandy sediments with gravel. Sandy sediments with high gas hydrate saturations (above 50 vol.%) were very minor constituents in the pressure-core sediments. The gas hydrate saturations were lower in sediments with a clay fraction above 30 vol.%. SEM images revealed microfossil debris, fragments of sedimentary rock, and silt/sand-sized-agglomerated grains throughout the pressure-cores, which indicate that the sediments were deposited at and/or transported to the foot of the continental slope through reworking of upper-slope deposits. The geological setting and SEM observations indicate that the pressure-core sediments are mainly a terrigenous source, with a small component from pelagic or hemipelagic sedimentation. The results of the pressure-core analyses confirm that natural gas hydrates can occur in silt-rich sediments. The occurrence of gas hydrates in the pressure-core sediments recovered during the NGHP-02 is controlled by the lithology of the sediments, especially their clay content.