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.     

Scale-dependent gas hydrate saturation estimates in sand reservoirs in the Ulleung Basin,East Sea of Korea

Through the use of 2-D and 3-D seismic data, several gas hydrate prospects were identified in the Ulleung Basin, East Sea of Korea and thirteen drill sites were established and logging-while-drilling (LWD) data were acquired from each site in 2010. Sites UBGH2–6 and UBGH2–10 were selected to test a series of high amplitude seismic reflections, possibly from sand reservoirs. LWD logs from the UBGH2–6 well indicate that there are three significant sand reservoirs with varying thickness. Two upper sand reservoirs are water saturated and the lower thinly bedded sand reservoir contains gas hydrate with an average saturation of 13%, as estimated from the P-wave velocity. The well logs at the UBGH2–6 well clearly demonstrated the effect of scale-dependency on gas hydrate saturation estimates. Gas hydrate saturations estimated from the high resolution LWD acquired ring resistivity (vertical resolution of about 5–8 cm) reaches about 90% with an average saturation of 28%, whereas gas hydrate saturations estimated from the low resolution A40L resistivity (vertical resolution of about 120 cm) reaches about 25% with an average saturation of 11%. However, in the UBGH2–10 well, gas hydrate occupies a 5-m thick sand reservoir near 135 mbsf with a maximum saturation of about 60%. In the UBGH2–10 well, the average and a maximum saturation estimated from various well logging tools are comparable, because the bed thickness is larger than the vertical resolution of the various logging tools. High resolution wireline log data further document the role of scale-dependency on gas hydrate calculations.     

Gas hydrate saturation in the Krishna–Godavari basin from P-wave velocity and electrical resistivity logs

During the Indian National Gas Hydrate Program (NGHP) Expedition 01, a series of well logs were acquired at several sites across the Krishna–Godavari (KG) Basin. Electrical resistivity logs were used for gas hydrate saturation estimates using Archie’s method. The measured in situ pore-water salinity, seafloor temperature and geothermal gradients were used to determine the baseline pore-water resistivity. In the absence of core data, Arp’s law was used to estimate in situ pore-water resistivity. Uncertainties in the Archie’s approach are related to the calibration of Archie coefficient (a), cementation factor (m) and saturation exponent (n) values. We also have estimated gas hydrate saturation from sonic P-wave velocity logs considering the gas hydrate in-frame effective medium rock-physics model. Uncertainties in the effective medium modeling stem from the choice of mineral assemblage used in the model. In both methods we assume that gas hydrate forms in sediment pore space. Combined observations from these analyses show that gas hydrate saturations are relatively low (<5% of the pore space) at the sites of the KG Basin. However, several intervals of increased saturations were observed e.g. at Site NGHP-01-03 (S_h = 15–20%, in two zones between 168 and 198 mbsf), Site NGHP-01-05 (S_h = 35–38% in two discrete zone between 70 and 90 mbsf), and Site NGHP-01-07 shows the gas hydrate saturation more than 25% in two zones between 75 and 155 mbsf. A total of 10 drill sites and associated log data, regional occurrences of bottom-simulating reflectors from 2D and 3D seismic data, and thermal modeling of the gas hydrate stability zone, were used to estimate the total amount of gas hydrate within the KG Basin. Average gas hydrate saturations for the entire gas hydrate stability zone (seafloor to base of gas hydrate stability), sediment porosities, and statistically derived extreme values for these parameters were defined from the logs. The total area considered based on the BSR seismic data covers ∼720 km². Using the statistical ranges in all parameters involved in the calculation, the total amount of gas from gas hydrate in the KG Basin study area varies from a minimum of ∼5.7 trillion-cubic feet (TCF) to ∼32.1 TCF.     

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.     

印度克里希纳-戈达瓦里盆地天然气水合物饱和度特征

印度国家天然气水合物计划(NGHP01)于2016年实施第1次钻探,证实了天然气水合物在印度大陆边缘的广泛分布。选择位于克里希纳-戈达瓦里盆地(KG盆地)的NGHP01-07D和NGHP01-15A钻孔,基于测井数据和岩心样品估算天然气水合物饱和度,分析天然气水合物赋存状态并探讨其形成机制。基于各向同性介质模型利用电阻率和声波测井计算NGHP01-15A钻孔的天然气水合物饱和度为0. 2%～33. 0%,平均值为9. 6%,在NGHP01-07D钻孔利用电阻率计算获得的天然气水合物饱和度高于岩心氯离子异常和气体释放获得的结果,但是基于各向同性岩石物理模型利用声波测井计算的天然气水合物饱和度与岩心结果一致,平均值为5. 0%。前人研究认为NGHP01-10D钻孔中天然气水合物以相对较高饱和度富集在高角度裂隙中。结合前人研究结果推断在克里希纳-戈达瓦里盆地存在3种不同的天然气水合物储层,即泥岩中砂质夹层各向同性储层、泥质/粉砂质高角度低连通性的低饱和度裂隙储层和泥岩中高角度高连通性的高饱和度裂隙储层,并提出对应的3种天然气水合物储层模型。     

Gas hydrate saturation in the Krishna–Godavari basin from P-wave velocity and electrical resistivity logs

During the Indian National Gas Hydrate Program (NGHP) Expedition 01, a series of well logs were acquired at several sites across the Krishna–Godavari (KG) Basin. Electrical resistivity logs were used for gas hydrate saturation estimates using Archie’s method. The measured in situ pore-water salinity, seafloor temperature and geothermal gradients were used to determine the baseline pore-water resistivity. In the absence of core data, Arp’s law was used to estimate in situ pore-water resistivity. Uncertainties in the Archie’s approach are related to the calibration of Archie coefficient (a), cementation factor (m) and saturation exponent (n) values. We also have estimated gas hydrate saturation from sonic P-wave velocity logs considering the gas hydrate in-frame effective medium rock-physics model. Uncertainties in the effective medium modeling stem from the choice of mineral assemblage used in the model. In both methods we assume that gas hydrate forms in sediment pore space. Combined observations from these analyses show that gas hydrate saturations are relatively low (<5% of the pore space) at the sites of the KG Basin. However, several intervals of increased saturations were observed e.g. at Site NGHP-01-03 (S_h = 15–20%, in two zones between 168 and 198 mbsf), Site NGHP-01-05 (S_h = 35–38% in two discrete zone between 70 and 90 mbsf), and Site NGHP-01-07 shows the gas hydrate saturation more than 25% in two zones between 75 and 155 mbsf. A total of 10 drill sites and associated log data, regional occurrences of bottom-simulating reflectors from 2D and 3D seismic data, and thermal modeling of the gas hydrate stability zone, were used to estimate the total amount of gas hydrate within the KG Basin. Average gas hydrate saturations for the entire gas hydrate stability zone (seafloor to base of gas hydrate stability), sediment porosities, and statistically derived extreme values for these parameters were defined from the logs. The total area considered based on the BSR seismic data covers ∼720 km². Using the statistical ranges in all parameters involved in the calculation, the total amount of gas from gas hydrate in the KG Basin study area varies from a minimum of ∼5.7 trillion-cubic feet (TCF) to ∼32.1 TCF.     

Formation pressure testing at the Mount Elbert Gas Hydrate Stratigraphic Test Well, Alaska North Slope: Operational summary, history matching, and interpretations

In February 2007, the U.S. Department of Energy, BP Exploration (Alaska), and the U.S. Geological Survey, collected open-hole pressure-response data, as well as gas and water sample collection, in a gas hydrate reservoir (the BPXA-DOE-USGS Mount Elbert Gas Hydrate Stratigraphic Test Well) using Schlumberger's Modular Dynamics Formation Tester (MDT) wireline tool. Four such MDT tests, ranging from six to twelve hours duration, and including a series of flow, sampling, and shut-in periods of various durations, were conducted. Locations for the testing were selected based on NMR and other log data to assure sufficient isolation from reservoir boundaries and zones of excess free water. Test stages in which pressure was reduced sufficiently to mobilize free water in the formation (yet not cause gas hydrate dissociation) produced readily interpretable pressure build-up profiles. Build-ups following larger drawdowns consistently showed gas-hydrate dissociation and gas release (as confirmed by optical fluid analyzer data), as well as progressive dampening of reservoir pressure build-up during sequential tests at a given MDT test station.History matches of one multi-stage, 12-h test (the C2 test) were accomplished using five different reservoir simulators: CMG-STARS, HydrateResSim, MH21-HYDRES, STOMP-HYD, and TOUGH + HYDRATE. Simulations utilized detailed information collected across the reservoir either obtained or determined from geophysical well logs, including thickness (11.3 m, 37 ft.), porosity (35%), hydrate saturation (65%), both mobile and immobile water saturations, intrinsic permeability (1000 mD), pore water salinity (5 ppt), and formation temperature (3.3-3.9 °C). This paper will present the approach and preliminary results of the history-matching efforts, including estimates of initial formation permeability and analyses of the various unique features exhibited by the MDT results.     

Anisotropic models to account for large borehole washouts to estimate gas hydrate saturations in the Gulf of Mexico Gas Hydrate Joint Industry Project Leg II Alaminos Canyon 21 B well

Through the use of 3-D seismic amplitude mapping, several gas hydrate prospects were identified in the Alaminos Canyon (AC) area of the Gulf of Mexico. Two locations were drilled as part of the Gulf of Mexico Gas Hydrate Joint Industry Project Leg II (JIP Leg II) in May of 2009 and a comprehensive set of logging-while-drilling (LWD) logs were acquired at each well site. LWD logs indicated that resistivity in the range of ∼2 ohm-m and P-wave velocity in the range of ∼1.9 km/s were measured in the target sand interval between 515 and 645 feet below sea floor. These values were slightly elevated relative to those measured in the sediment above and below the target sand. However, the initial well log analysis was inconclusive regarding the presence of gas hydrate in the logged sand interval, mainly because large washouts caused by drilling in the target interval degraded confidence in the well log measurements. To assess gas hydrate saturations in the sedimentary section drilled in the Alaminos Canyon 21 B (AC21-B) well, a method of compensating for the effect of washouts on the resistivity and acoustic velocities was developed. The proposed method models the washed-out portion of the borehole as a vertical layer filled with sea water (drilling fluid) and the apparent anisotropic resistivity and velocities caused by a vertical layer are used to correct the measured log values. By incorporating the conventional marine seismic data into the well log analysis, the average gas hydrate saturation in the target sand section in the AC21-B well can be constrained to the range of 8–28%, with 20% being our best estimate.     

Rock physics-based seismic trace analysis of unconsolidated sediments containing gas hydrate and free gas in Green Canyon 955, Northern Gulf of Mexico

The gas hydrate petroleum system at the 2009 Gulf of Mexico Gas Hydrate Joint Industry Project Leg II (JIP Leg II) Green Canyon 955 (GC955) site shows a complex seismic amplitude and waveform response of highly negative and positive amplitudes with continuous and discontinuous character within inferred gas-hydrate- and gas-bearing sand reservoirs. Logging-while-drilling (LWD) data obtained during JIP Leg II and conventional 3-D seismic data allowed for the identification of thick highly concentrated hydrate layers by integrating rock physics modeling, amplitude and thin layer analysis, and spectral decomposition. Rock physics modeling with constraints from three JIP LWD holes allowed for the analysis of variations in acoustic amplitude characteristics as a product of hydrate saturation, gas saturation, and reservoir thickness. Using the well log-derived acoustic models, thick highly concentrated gas hydrate with and without underlying free gas accumulations have been identified. These results suggest that thick highly concentrated gas-hydrate-bearing sand units (with thicknesses greater than half of the seismic tuning thickness and gas hydrate saturations greater than 50%) underlain by gas can be differentiated from sands containing only gas, but thin gas-hydrate-bearing sand units with low gas hydrate concentrations (with thicknesses less than half of the seismic tuning thickness and gas hydrate saturations less than 50%) are difficult to identify from post-stack seismic amplitude data alone. Within GC955, we have identified six zones with seismic amplitude anomalies interpreted as being caused by gas hydrate deposits with variable lateral extent, thickness and saturation, and in some cases overlying free-gas-bearing intervals. Synthetic seismic images produced from well-log- and model-derived velocity and density distributions mimic similar reflection characteristics in the corresponding field seismic data.     

An approach to estimate gas hydrate saturation from 3-D heterogeneous resistivity model: A study from Krishna-Godavari basin,Eastern Indian offshore

A method is proposed to estimate gas hydrate saturation from three dimensional (3-D) heterogeneous model of resistivity simulated using resistivity log. Pure gas hydrates are highly resistive compared to the host sediments, and their presence in the pore space of sediments increase the resistivity of the formation. The anomalous increase of resistivity is used as a proxy for the delineation of gas hydrates using the resistivity log. A 3-D heterogeneous resistivity model has been constructed from one dimensional resistivity log in Krishna Godavari (KG) basin, eastern Indian offshore. The simulated model contains all small scale variation in resistivity of the reservoir and maintains all properties associated with covariance, like root mean square fluctuation, characteristic scales and fractal dimension of the observed log. We have estimated volumetric hydrate saturation using the three dimensional simulated model. The porosity used for estimating hydrate saturation is calculated from the simulated density field generated using the observed density log. Estimated average gas hydrate saturation is about 9.84% of the pore volume over a 1000 m × 1000 m x 131 m cubic meters area.     

Large-scale depositional characteristics of the Ulleung Basin and its impact on electrical resistivity and Archie-parameters for gas hydrate saturation estimates

Gas hydrate saturation estimates were obtained from an Archie-analysis of the Logging-While-Drilling (LWD) electrical resistivity logs under consideration of the regional geological framework of sediment deposition in the Ulleung Basin, East Sea, of Korea. Porosity was determined from the LWD bulk density log and core-derived values of grain density. In situ measurements of pore-fluid salinity as well as formation temperature define a background trend for pore-fluid resistivity at each drill site. The LWD data were used to define sets of empirical Archie-constants for different depth-intervals of the logged borehole at all sites drilled during the second Ulleung Basin Gas Hydrate Drilling Expedition (UBGH2). A clustering of data with distinctly different trend-lines is evident in the cross-plot of porosity and formation factor for all sites drilled during UBGH2. The reason for the clustering is related to the difference between hemipelagic sediments (mostly covering the top ∼100 mbsf) and mass-transport deposits (MTD) and/or the occurrence of biogenic opal. For sites located in the north-eastern portion of the Ulleung Basin a set of individual Archie-parameters for a shallow depth interval (hemipelagic) and a deeper MTD zone was achieved. The deeper zone shows typically higher resistivities for the same range of porosities seen in the upper zone, reflecting a shift in sediment properties. The presence of large amounts of biogenic opal (up to and often over 50% as defined by XRD data) was especially observed at Sites UBGH2-2_1 and UBGH2-2_2 (as well as UBGH1-9 from a previous drilling expedition in 2007). The boundary between these two zones can also easily be identified in gamma-ray logs, which also show unusually low readings in the opal-rich interval. Only by incorporating different Archie-parameters for the different zones a reasonable estimate of gas hydrate saturation was achieved that also matches results from other techniques such as pore-fluid freshening, velocity-based calculations, and pressure-core degassing experiments. Seismically, individual boundaries between zones were determined using a grid of regional 2D seismic data. Zoning from the Archie-analysis for sites in the south-western portion of the Ulleung Basin was also observed, but at these sites it is linked to individually stacked MTDs only and does not reflect a mineralogical occurrence of biogenic opal or hemipelagic sedimentation. The individual MTD events represent differently compacted material often associated with a strong decrease in porosity (and increase in density), warranting a separate set of empirical Archie-parameters.     

Low-amplitude BSRs and gas hydrate concentration on the northern margin of the South China Sea

The passive northern continental margin of the South China Sea is rich in gas hydrates, as inferred from the occurrence of bottom-simulating reflectors (BSR) and from well logging data at Ocean Drilling Program (ODP) drill sites. Nonetheless, BSRs on new 2D multichannel seismic reflection data from the area around the Dongsha Islands (the Dongsha Rise) are not ubiquitous. They are confined to complex diapiric structures and active fault zones located between the Dongsha Rise and the surrounding depressions, implying that gas hydrate occurrence is likewise limited to these areas. Most of the BSRs have low amplitude and are therefore not clearly recognizable. Acoustic impedance provides information on rock properties and has been used to estimate gas hydrate concentration. Gas hydrate-bearing sediments have acoustic impedance that is higher than that of the surrounding sediments devoid of hydrates. Based on well logging data, the relationship between acoustic impedance and porosity can be obtained by a linear regression, and the degree of gas hydrate saturation can be determined using Archie’s equation. By applying these methods to multichannel seismic data and well logging data from the northern South China Sea, the gas hydrate concentration is found to be 3–25% of the pore space at ODP Site 1148 depending on sub-surface depth, and is estimated to be less than values of 5% estimated along seismic profile 0101. Our results suggest that saturation of gas hydrate in the northern South China Sea is higher than that estimated from well resistivity log data in the gas hydrate stability zone, but that free gas is scarce beneath this zone. It is probably the scarcity of free gas that is responsible for the low amplitudes of the BSRs.     

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.     

Evidence of gas hydrate from downhole logging data in the Ulleung Basin, East Sea

The Gas Hydrate Research and Development Organization (GHDO) of Korea successfully accomplished both coring (hydraulic piston and pressure coring) and logging (logging-while-drilling, LWD, and wireline logging) to investigate the presence of gas hydrate during the first deep drilling expedition in the Ulleung Basin, East Sea of Korea (referred to as UBGH1) in 2007. The LWD data from two sites (UBGH1-9, UBGH1-10) showed elevated electrical resistivity (>80 Ω-m) and P-wave velocity (>2000 m/s) values indicating the presence of gas hydrate. During the coring period, the richest gas hydrate accumulation was discovered at these intervals. Based on log data, the occurrence of gas hydrate is primarily controlled by the presence of fractures. The gas hydrate saturation calculated using Archie’s relation shows greater than 60% (as high as ∼90%) of the pore space, although Archie’s equation typically overestimates gas hydrate saturation in near-vertical fractures. The saturation of gas hydrate is also estimated using the modified Biot-Gassmann theory (BGTL) by Lee and Collett (2006). The saturation values estimated rom BGTL are much lower than those calculated from Archie’s equation. Based on log data, the hydrate-bearing sediment section is approximately 70 m (UBGH1-9) to 130 m (UBGH1-10) in thickness at these two sites. This was further directly confirmed by the recovery of gas hydrate samples and pore water freshening collected from deep drilling core during the expedition. LWD data also strongly support the interpretation of the seismic gas hydrate indicators (e.g., vent or chimney structures and bottom-simulating reflectors), which imply the probability of widespread gas hydrate presence in the Ulleung Basin.     

Estimating saturation of gas hydrates using conventional 3D seismic data,Gulf of Mexico Joint Industry Project Leg II

We present a methodology for generating pre-drill estimations of hydrate saturations using conventional 3D seismic data. These seismic-based estimates will be compared with well log derived saturations from the subsequently drilled wells of the Gulf of Mexico Gas Hydrate Joint Industry Project Leg II (JIP Leg II) expedition.Predicting saturation of gas hydrates (S_h-seismic) combines pre-stack seismic inversion, rock physics modeling and stratigraphic interpretation. Before the wells were drilled, no nearby sonic and density logs were available to define and calibrate the elastic property trends for the shallow target interval containing the gas hydrate stability zone. Therefore, rock property trends were established by applying principles of rock physics and shallow sediment compaction, constrained by known regional geological parameters. S_h-seismic volumes were generated by inverting pre-stack data to acoustic impedance (PI) and shear impedance (SI) volumes, and then analyzing deviations from modeled impedance trends. To enhance the quality of the inversion, the signal-to-noise ratio of the offset data was maximized by conditioning the seismic prior to inversion. Seismic stratigraphic interpretation plays an important role by identifying the more promising strata and structures for the presence of gas hydrates.The pre-drill S_h-seismic results are compared with saturations calculated from log data, S_h-log, of the wells drilled in the JIP Leg II campaign. Due to weaker seismic reflections, predictions may be less accurate for low concentrations, such as saturations less than 40%, and for thin intervals below the vertical resolution of the seismic data (about 15 m). However, the integrated geophysical workflow is very effective for identifying and quantifying significant hydrate concentrations, making the method a promising prospecting technique.     

南海神狐海域天然气水合物饱和度的数值模拟分析

珠江口盆地神狐海域是天然气水合物钻探和试验开采的重点区域,大量钻探取心、测井与地震等综合分析表明不同站位水合物的饱和度、厚度与气源条件存在差异。本文利用天然气水合物调查及深水油气勘探所采集的测井和地震资料建立地质模型,利用PetroMod软件模拟地层的温度场、有机质成熟度、烃源岩生烃量、流体运移路径以及不同烃源岩影响下的水合物饱和度,结果表明:生物成因气分布在海底以下1500 m范围内的有机质未成熟地层,而热成因气分布在深度超过2300 m的成熟、过成熟地层。水合物稳定带内生烃量难以形成水合物,形成水合物气源主要来自于稳定带下方向上运移的生物与热成因气。模拟结果与测井结果对比分析表明,稳定带下部生物成因气能形成的水合物饱和度约为10%,在峡谷脊部的局部区域饱和度较高;相对高饱和度(>40%)水合物形成与文昌组、恩平组的热成因气沿断裂、气烟囱等流体运移通道幕式释放密切相关,W19井形成较高饱和度水合物的甲烷气体中热成因气占比达80%,W17井热成因气占比为73%,而SH2井主要以生物成因为主,因此,不同站位甲烷气体来源占比不同。     

Estimation of seismic velocities and gas hydrate concentrations: A case study from the Shenhu area,northern South China Sea

In the Shenhu area of the northern South China Sea (SCS), canyon systems and focused fluid flow systems increase the complexity of the gas hydrate distribution in the region. It also induces difficulties in predicting the hydrate reservoir characteristics and quantitatively evaluating reservoir parameters. In this study, several inversion methods have been executed to estimate the velocities of strata and gas hydrate concentrations along a profile in the Shenhu area. The seismic data were inverted to obtain the reflection coefficient of each stratum via a spectral inversion method. Stratigraphic horizons were then delineated by tracking the inverted reflectivities. Based on the results of spectral inversion, a low-frequency velocity field of the strata was constructed for acoustic impedance inversion. Using a new iterative algorithm for acoustic impedance inversion, reflection coefficients were converted into velocities, and the velocity variations of the strata along a 2D seismic line were then obtained. Subsequently, gas hydrate saturations at well SH2 were estimated via the shale-corrected resistivity method, the chloride ion concentration method and three different rock physics models. The results were then compared to determine the optimal rock physics model, and the modified Wood equation (MWE) was found to be appropriate for this area. Finally, the inverted velocities and MWE were used to predict the distribution and concentrations of gas hydrates along the seismic line. The estimated spatial distribution of gas hydrates is consistent with that from sonic logging and resistivity data at well SH2, and with the drilling results. Therefore, this method is applicable in areas with no well data, or with few wells, and provides an effective tool for predicting and evaluating gas hydrates using seismic data.