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

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

Physical and Acoustic Properties of Gas-bearing Sediments in Jinhae Bay,the South Sea of Korea

High-resolution seismic survey and sediment core sampling were conducted to investigate acoustic characteristics of gas-bearing sediments in Jinhae Bay, the southeast of Korea. The sediment in Jinhae Bay is mostly homogenous mud deposited after the Holocene transgression. Along with the 410 km of chirp seismic profiling, five piston core samples were collected on the track lines.

Gassy sediments are common and occur widely in the bay. Core samples were analyzed for sediment texture, physical properties (porosity, water content, bulk density, and grain density), acoustic properties (compressional wave velocity and attenuation), and electrical resistivity. X-radiograph image analysis was also performed to observe the shape of degassing cracks. There is no significant downcore variation on physical and sediment textures regardless of existence of gas bubbles. However, compressional wave velocity dramatically decreases from average 1480 to 1380～739 m/s for the cores that penetrate the gas-bearing zones. This is probably due to degassying cracks that developed by escaping gases and free gas bubbles that are still trapped in the cores. Electrical resistivity is the only geotechnical property that increases in the gas-bearing zone where compressional wave velocity abruptly decreases. This indicates the possibility of using both electrical resistivity as an index variable as well as to compressional wave velocity to identify gassy sediment microstructure because there are little changes in texture and composition of sediment.     

Gas hydrate and free gas distribution from inversion of seismic data on the South Shetland margin (Antarctica)

Travel-time inversion is applied to seismic data to produce acoustic velocity images of the upper 800 m of the South Shetland margin (Antarctic Peninsula) in three different geological domains: (i) the continental shelf; (ii) the accretionary prism; (iii) the trench. The velocity in the continental shelf sediments is remarkably higher, up to 1000 m/s at 600–700 m below seafloor, than that of the other two geological domains, due to the sediment overcompaction and erosion induced by the wax and waning of a grounded ice sheet. Pre-stack depth migration was applied to the data in order to improve the seismic image and to test the quality of the velocity fields. Where the Bottom Simulating Reflector (BSR) is present, positive and negative velocity anomalies were found with respect to a reference empirical velocity profile. The 2D-velocity section was translated in gas hydrate and free gas distribution by using a theoretical approach. The analysis revealed that the BSR is mainly related to the presence of free gas below it. The free gas is distributed in the area with variable concentration and thickness, while the gas hydrate is quite uniformly distributed across the margin.     

Tomographic approach for gas hydrate investigation in Kerala-Konkan region,India

Seismic tomography is an effective means of estimating velocity and structure from multichannel seismic (MCS) reflection data. In this study we have followed a 2D approach to arrive at the probable velocity field configuration from multichannel seismic data and infer the presence of gas hydrates/free-gas in the offshore Kerala-Konkan region, along the eastern part of a seismic line on which a bottom simulating reflector (BSR) has previously been identified. Tomographic modeling consists of the identification of reflection phases and picking of respective travel times for various source-receiver positions. These picks were then utilized to arrive at a 2D velocity field following a forward and inversion approach using a ray tracing technique. The modeling for the first time brought out the finer scale velocity structure under the region of investigation. Modeling through the 2D approach shows lateral variation in velocity field along the studied segment of the seismic line. The results indicate a thin (∼50–60 m) sedimentary cover with velocity ranging from 1,770 to 1,850 m/s. A sedimentary layer with high P-wave velocity 1,980–2,100 m/s below the sea floor was interpreted as the hydrate layer. The thickness of this layer varies between 110 and 140 m. The hydrate layer is underlain by a low-velocity layer having velocities in the range 1,660–1,720 m/s. This low velocity may represent a free gas layer, whose thickness varies between 50 and 100 m located below the hydrated layer. The investigation suggests the occurrence of gas hydrate underlain by free gas in some parts of the Kerala-Konkan offshore region.     

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

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

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.     

Methodology for wide-angle seismic modeling of marine gas hydrate deposits

A wide angle seismic modeling methodology was developed for the purpose of studying marine gas hydrate deposits. The software for seismic modeling was selected on the basis of comparative analysis and testing of different computer programs. Five averaged prognostic two-dimensional models were developed. These models include the gas hydrate zone, the free gas zone, and the basement. The prognostic models suggested represent the structure of gas hydrate deposits for various regions of the World Ocean. Wave field calculations were made for various positions of ocean bottom seismometers with respect to the gas hydrate zone using the seismic tomography technique. Numerical experiments showed significant anomalies of the kinematical and dynamical characteristics of the refracted and wide-angle reflected waves. These anomalies are related to the gas hydrate and the free gas zones and to a possible channel of hydrocarbon supply.     

南海北部揭阳凹陷天然气水合物的地震异常特征分析

大量研究表明南海北部珠江口盆地是天然气水合物发育区,但是该盆地东部揭阳凹陷水合物研究较少。本文利用揭阳凹陷新采集三维地震资料,对该三维地震资料进行成像道集优化和叠前时间偏移处理,得到针对水合物的新处理地震数据体,并通过高精度网格层析反演得到层速度数据体。利用该数据开展叠后约束稀疏脉冲反演,获得含天然气水合物地层波阻抗异常,综合分析反演与地震属性识别水合物。从新处理地震资料看,该区域似海底反射(bottom simulation reflection,BSR)反射呈连续、不连续与地层斜交等特征,BSR发育在一个继承性小型水道上,且下部断裂和气烟囱发育。通过分析BSR特征及BSR上下地层的速度、波阻抗、振幅、频率、相干等属性异常,结合水合物成藏条件,发现了南海北部新的天然气水合物有利富集区,为该区域水合物勘探提供基础。     

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.     

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.     

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.     

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.     

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

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

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.     

A seismic survey to detect natural gas hydrate in the East Sea of Korea

Recently, several countries have conducted projects to explore and develop natural gas hydrate, which is one of the new alternative energy resources for the future. In Korea, a five-year national research project was initiated in 2000. As part of this project, a seismic survey was performed in the East Sea of Korea to quantify the potential magnitude and distribution of natural gas hydrates. Multi-channel seismic data and core samples have been acquired and recovered in the survey area. Analysis of seismic data show clear bottom simulating reflectors (BSRs), seismic blank zones (or wipe-out zones) with velocity pull-up structure, and pock-marks. In this study, we present the results of seismic surveys which indicate the existence of natural gas hydrates in Korean offshore areas. These results will be applied to select areas for coring (or drilling) and detailed exploration such as 2D seismic survey with long offset or 3D seismic in the future.     

High-resolution P-Cable 3D seismic imaging of gas chimney structures in gas hydrated sediments of an Arctic sediment drift

The newly developed P-Cable 3D seismic system allows for high-resolution seismic imaging to characterize upper geosphere geological features focusing on geofluid expressions (gas chimneys), shallow gas and gas hydrate reservoirs. Seismic imaging of a geofluid system of an Arctic sediment drift at the Vestnesa Ridge, offshore western Svalbard, provides significantly improved details of internal chimney structures from the seafloor to ∼500 m bsf (below seafloor). The chimneys connect to pockmarks at the seafloor and indicate focused fluid flow through gas hydrated sediments. The pockmarks are not buried and align at the ridge-crest pointing to recent, topography-controlled fluid discharge. Chimneys are fuelled by sources beneath the base of gas hydrate stability zone (GHSZ) that is evident at ∼160–170 m bsf as indicated by a bottom-simulating reflector (BSR). Conduit centres that are not vertically straight but shift laterally by up to 200 m as well as discontinuous internal chimney reflections indicate heterogeneous hydraulic fracturing of the sediments. Episodically active, pressure-driven focused fluid flow could explain the hydro-fracturing processes that control the plumbing system and lead to extensive pockmark formation at crest of the Vestnesa Ridge. High-amplitude anomalies in the upper 50 m of the chimney structures suggest formations of near-surface gas hydrates and/or authigenic carbonate precipitation. Acoustic anomalies, expressed as high amplitudes and amplitude blanking, are irregularly distributed throughout the deeper parts of the chimneys and provide evidence for the variability of hydrate and/or carbonate formation in space and time.     

南海北部荔湾3区块天然气水合物分布特征及目标识别

针对天然气水合物钻探与取样难以解决的水合物矿体空间展布等问题,利用白云-荔湾凹陷高密度分析重新处理的三维地震资料,首先基于模糊数学的多属性融合技术对水合物分布进行刻画;再通过高分辨率速度场对浅层开展高分辨率宽频无井反演技术,提高了水合物层分辨率;最后,利用岩石物理方法及多种模型对水合物饱和度进行定量预测,实现了对5～6m厚水合物层的有效辨别,进而形成了一套适合于孔隙充填型的水合物矿藏目标识别评方法。结果表明:应用该技术可有效对荔湾3水合物富集区第四条带水合物空间刻画,揭示出该区水合物饱和度最高可超40%,同时薄层与厚层水合物具有明显互层分布特征,在水合物矿体刻画及饱和度预测基础上,进一步对该区实施了井位优选,该方法预测的水合物层与实际钻探H1和H2站位吻合较好。这些结果说明常规三维油气地震数据在经过宽频处理后可应用于高分辨率水合物勘探,节约经济成本,同时提高了常规地震在水合物勘探中精度与实用性。     

Evidence for gas hydrate occurrences in the Canadian Arctic Beaufort Sea within permafrost-associated shelf and deep-water marine environments

The presence of a wedge of offshore permafrost on the shelf of the Canadian Beaufort Sea has been previously recognized and the consequence of a prolonged occurrence of such permafrost is the possibility of an underlying gas hydrate regime. We present the first evidence for wide-spread occurrences of gas hydrates across the shelf in water depths of 60–100 m using 3D and 2D multichannel seismic (MCS) data. A reflection with a polarity opposite to the seafloor was identified ∼1000 m below the seafloor that mimics some of the bottom-simulating reflections (BSRs) in marine gas hydrate regimes. However, the reflection is not truly bottom-simulating, as its depth is controlled by offshore permafrost. The depth of the reflection decreases with increasing water depth, as predicted from thermal modeling of the late Wisconsin transgression. The reflection crosscuts strata and defines a zone of enhanced reflectivity beneath it, which originates from free gas accumulated at the phase boundary over time as permafrost and associated gas hydrate stability zones thin in response to the transgression. The wide-spread gas hydrate occurrence beneath permafrost has implications on the region including drilling hazards associated with the presence of free gas, possible overpressure, lateral migration of fluids and expulsion at the seafloor. In contrast to the permafrost-associated gas hydrates, a deep-water marine BSR was also identified on MCS profiles. The MCS data show a polarity-reversed seismic reflection associated with a low-velocity zone beneath it. The seismic data coverage in the southern Beaufort Sea shows that the deep-water marine BSR is not uniformly present across the entire region. The regional discrepancy of the BSR occurrence between the US Alaska portion and the Mackenzie Delta region may be a result of high sedimentation rates expected for the central Mackenzie delta and high abundance of mass-transport deposits that prohibit gas to accumulate within and beneath the gas hydrate stability zone.     

Characteristics of Gas Hydrate and Free Gas Offshore Southwestern Taiwan from a Combined MCS/OBS Data Analysis

In this study, we present the results of the combined analyses of ocean bottom seismometer and multi-channel seismic reflection data collection offshore southwestern Taiwan, with respect to the presence of gas hydrates and free gas within the accretionary wedge sediments. Estimates of the compressional velocities along EW9509-33 seismic reflection profile are obtained by a series of pre-stack depth migrations in a layer stripping streamlined Deregowski loop. Strong BSR is imaged over most of the reflection profile while low velocity zones are imaged below BSR at several locations. Amplitude versus angle analysis that are performed within the pre-stack depth migration processes reveal strong negative P-impedance near the bottom of the hydrate stability zone, commonly underlain by sharp positive P impedance layers associated with negative pseudo-Poisson attribute areas, indicating the presence of free gas below the BSR. Ray tracing of the acoustic arrivals with a model derived from the migration velocities generally fits the vertical and hydrophone records of the four ocean-bottom seismographs (OBS). In order to estimate the Poisson’s ratios in the shallow sediments at the vicinity of the OBSs, we analyze the mode-converted arrivals in the wide-angle horizontal component. P-S mode converted reflections are dominant, while upward P-S transmissions are observed at large offsets. We observe significant compressional velocity and Poisson’s ratio pull-down in the sediment below the BSR likely to bear free gas. When compared to Poisson’s ratio predicted by mechanical models, the values proposed for the OBSs yield rough estimates of gas hydrate saturation in the range of 0–10% in the layers above the BSR and of free gas saturation in the range of 0–2% just below the BSR.