# Evaluation Method of the Gas Hydrate and Free Gas System and Its Application in the Shenhu Area, South China Sea

Authors: Pibo Su; Tingwei Li; Shurong Liang; Jinqiang Liang; Xiaoxue Wang; Xiaoming Wan; Feifei Wang; Fang Liu
Journal: Geofluids
 (2021)

---

## Abstract

As a new alternative energy source, gas hydrate has attracted wide attention all over the world. Since gas hydrate is always associated with free gas, the evaluation of the gas hydrate and free gas system is an important aspect of hydrate reservoir exploration and development. In this study, based on identifying gas hydrate and free gas by well logging, the seismic reflection characteristics of gas hydrate and free gas are determined by an accurate well-to-seismic calibration method. On account of seismic reflection characteristics, AVO attributes are used to identify gas hydrate and free gas qualitatively. Using prestack and poststack inversion to get the ratio ofP-wave impedance and P-wave-to-S-wave velocities, we determine the three-dimensional space distribution of gas hydrate and free gas, predict their effective porosity and saturation, and eventually achieve the meticulous depiction of gas hydrate and free gas in the body, which is necessary in subsequent estimation of gas hydrate and free gas resources. Results show that according to logging interpretation, gas hydrate of the B-well is located in the depth range of 1460–1510 mbsl and free gas is in 1510–1542 mbsl. Moreover, gas hydrate of the A-well is located in the depth range of 1425–1512 mbsl, and no obvious free gas is identified. Gas hydrate is located above free gas and distributed continuously. In plane form, gas hydrate and free gas both present subelliptical distribution in the NW-SE direction. Gas hydrate has an effective porosity of 0.30–0.40, an average saturation of 0.33–0.40, and an effective thickness of 3.0–10.5 m, whereas free gas possesses an effective porosity of 0.35–0.40, a saturation of 0.24–0.32, and an effective thickness of 2.0–5.0 m.


---

## Body

## 1. Introduction

Natural gas hydrate is a kind of crystalline compounds composed of natural gas and water under high pressure and low temperature, which is distributed in deep-water sediments and the permafrost [1–3]. Gas hydrate is considered one of the most potential new clean energies to replace coal, oil, and gas in the 21st century, and it is also a new energy with abundant reserves that have not been fully exploited [4–6]. Since natural gas hydrate is a kind of clean energy with great potential and has a positive impact on greenhouse effect and marine ecosystem, the United States, Canada, Japan, South Korea, and India are committed to the exploration and development of natural gas hydrate [4, 7–13].There is a great prospect of gas hydrate resources in China’s vast jurisdictional sea, exclusive economic zone, and permafrost region [14–16]. The South China Sea has favorable geological, geomorphic, and geochemical conditions for natural gas hydrate accumulation [17–24]. Regional marine geological survey and oil and gas exploration show that the South China Sea is rich in natural gas hydrate resources [17, 20, 23, 25–29].Recent studies have shown that gas hydrate and free gas coexist in the Hydrate Ridge in offshore Oregon [30], the Green Canyon in the Gulf of Mexico [31], the Blake Ridge in North America [32], and the Nankai Trough in Japan [33]. Researches on gas hydrate reservoirs in the Shenhu area, South China Sea, have systematically discussed the geological and geochemical characteristics of gas hydrate and associated free gas [34, 35].Although the above studies infer the phenomenon and formation mechanism of gas hydrate and associated free gas, there is no detailed description of natural gas hydrate layer and associated free gas layer. Therefore, in this study, we use an accurate well-to-seismic calibration method to fine depict gas hydrate and free gas reservoirs and compare the differences between them, which lays a foundation for subsequent estimation of gas hydrate and free gas resources.


## 2. Geological Setting

The Shenhu area is geographically located in the northern slope area of the South China Sea and structurally situated in the Baiyun sag, Zhuer depression, and Peal River Mouth Basin (Figure1(a)). Its seafloor topography is relatively flat, with an average slope of 3°. There are many types of landforms in the Shenhu area, including sea knolls, sea valleys, erosion slots, and gullies [34, 36]. The water depth is between 1000 and 1700 mbsl, and the Shenhu area is divided into two parts, the north and south, with the water depth of 1350 mbsl as the boundary (Figure 1(b)). In the north part, the terrain is relatively steep, and there are 17 nearly north-south trending submarine trenches from west to east, which are arranged alternately with seamounts. In the south part, the terrain is flat and gradually becomes abyssal plain to the south. The A-B ore body is taken as the interested object in this study, which is located in the submarine canyon systems in the northern slope of the Baiyun sag (Figure 1(c)). The A-well and B-well are situated at the top of the eastern ridges, with water depths of 1309.75 mbsl and 1259.00 mbsl, total drilling depths of ~222 mbsf and ~230 mbsf, and geothermal gradients of 5.46°C/100 m and 4.43°C/100 m, respectively [3].Figure 1
Regional geological background of the Shenhu area in the north of South China Sea.During the Cenozoic, the Pearl River Mouth Basin experienced several regional tectonic activities, including Shenhu movement, Zhuqiong movement, Nanhai movement, Baiyun movement, and Dongsha movement (Figure2). The Cenozoic structural evolution of the Baiyun sag is characterized by two periods, i.e., rifting stage and depression stage. From the Paleocene to the early Oligocene, a half-graben or graben depression was formed, with the continental strata of the Shenhu formation in the Paleocene, Wenchang formation in the Eocene, and Enping formation in the lower Oligocene deposited. In the late Oligocene, it became a depression stage and deposited neritic facies strata of the Zhuhai formation. Subsequently, in the Miocene and Pliocene, it deposited neritic facies strata. There is an obvious unconformity between the lower Oligocene and the upper Oligocene. From the bottom to top, the Shenhu area develops successively continental facies, transitional facies, and marine facies strata, showing a trend of transgression [37].Figure 2
Comprehensive stratigraphic histogram of the Pearl River Mouth Basin.


## 3. Methods

### 3.1. Logging Identification Method of Natural Gas Hydrate and Free Gas

Logging response characteristics are the basis of identifying natural gas hydrate and free gas. Besides analyzing the response pattern of each logging curve, we adopt the curve overlap method to qualitatively identify natural gas hydrate and free gas. According to the overlap method, when gas was saturated in the reservoir, interval transit time increases, whereas density porosity and neutron porosity decrease. Through reasonable standardization, interval transit time and density porosity can overlap with neutron porosity in the water layer but separate in the gas layer [3].


### 3.2. Sensitivity Analysis of Gas Hydrate and Free Gas Elasticity Parameters

Elastic parameters are the link between well logging and seismic data, the sensitivity analysis of which is the basis of the seismic prediction method. Only when the elastic parameters of gas hydrate and free gas have a certain range, elastic inversion constrained by well logging is feasible [38].Through the intersection analysis of different elastic parameters, it can be found that gas hydrate is characterized by highP-wave impedance (P‐imp>3000g/cm3·m/s), while free gas has an obvious characteristic of low P-wave impedance (P‐imp<3000g/cm3·m/s) (Figure 3).Figure 3
Intersection plate of different elastic parameters.
(a)(b)(c)Therefore,P-wave impedance is one of the most sensitive elastic parameters for gas hydrate and free gas. Through P-wave impedance data, gas hydrate and free gas could be identified accurately.There is a strong difference inP-wave impedance values between gas hydrate and free gas. Using the constrained sparse spike inversion (CSSI) method in the software Jason, under the control of the gas hydrate-bearing zone and bottom simulating reflector (BSR), the entire logging data for the hydrate ore body in the study area was used to complete the poststacking CSSI processing and obtain the P-wave impedance, which revealed the distinct distribution of the GHBZ and underlying free gas. The top of gas hydrate is determined by strong amplitude reflection consistent with the polarity of the seafloor, and the bottom is identified by BSR.


### 3.3. Correlation Analysis of Physical Properties and Elasticity Parameters for Gas Hydrate and Free Gas

Through analyzing the correlation between physical properties and elasticity parameters of gas hydrate and free gas, it could be found that there is a good correlation between effective porosity andP-wave impedance, whereas saturation and P-S wave velocity ratio have a good relationship (Figure 4). The fitting relationship is as follows:
(1)Effectiveporosity=−0.43×lnP‐waveimpedance+3.85R2=0.82,Freegassaturation=0.84×P‐Swavevelocityratio2−3.84×P‐Swavevelocityratio+4.52R2=0.77,Gashydratesaturation=−0.29×P‐Swavevelocityratio2+1.47×P‐Swavevelocityratio−1.42R2=0.77.Figure 4
Correlation plate of physical properties and elastic parameters for gas hydrate and free gas.
(a)(b)(c)Using the above formulae, the spatial distribution of effective porosity, free gas saturation, and gas hydrate saturation of the A-B ore body could be well predicted based on theP-wave impedance and P-S wave velocity ratio.


## 3.1. Logging Identification Method of Natural Gas Hydrate and Free Gas

Logging response characteristics are the basis of identifying natural gas hydrate and free gas. Besides analyzing the response pattern of each logging curve, we adopt the curve overlap method to qualitatively identify natural gas hydrate and free gas. According to the overlap method, when gas was saturated in the reservoir, interval transit time increases, whereas density porosity and neutron porosity decrease. Through reasonable standardization, interval transit time and density porosity can overlap with neutron porosity in the water layer but separate in the gas layer [3].


## 3.2. Sensitivity Analysis of Gas Hydrate and Free Gas Elasticity Parameters

Elastic parameters are the link between well logging and seismic data, the sensitivity analysis of which is the basis of the seismic prediction method. Only when the elastic parameters of gas hydrate and free gas have a certain range, elastic inversion constrained by well logging is feasible [38].Through the intersection analysis of different elastic parameters, it can be found that gas hydrate is characterized by highP-wave impedance (P‐imp>3000g/cm3·m/s), while free gas has an obvious characteristic of low P-wave impedance (P‐imp<3000g/cm3·m/s) (Figure 3).Figure 3
Intersection plate of different elastic parameters.
(a)(b)(c)Therefore,P-wave impedance is one of the most sensitive elastic parameters for gas hydrate and free gas. Through P-wave impedance data, gas hydrate and free gas could be identified accurately.There is a strong difference inP-wave impedance values between gas hydrate and free gas. Using the constrained sparse spike inversion (CSSI) method in the software Jason, under the control of the gas hydrate-bearing zone and bottom simulating reflector (BSR), the entire logging data for the hydrate ore body in the study area was used to complete the poststacking CSSI processing and obtain the P-wave impedance, which revealed the distinct distribution of the GHBZ and underlying free gas. The top of gas hydrate is determined by strong amplitude reflection consistent with the polarity of the seafloor, and the bottom is identified by BSR.


## 3.3. Correlation Analysis of Physical Properties and Elasticity Parameters for Gas Hydrate and Free Gas

Through analyzing the correlation between physical properties and elasticity parameters of gas hydrate and free gas, it could be found that there is a good correlation between effective porosity andP-wave impedance, whereas saturation and P-S wave velocity ratio have a good relationship (Figure 4). The fitting relationship is as follows:
(1)Effectiveporosity=−0.43×lnP‐waveimpedance+3.85R2=0.82,Freegassaturation=0.84×P‐Swavevelocityratio2−3.84×P‐Swavevelocityratio+4.52R2=0.77,Gashydratesaturation=−0.29×P‐Swavevelocityratio2+1.47×P‐Swavevelocityratio−1.42R2=0.77.Figure 4
Correlation plate of physical properties and elastic parameters for gas hydrate and free gas.
(a)(b)(c)Using the above formulae, the spatial distribution of effective porosity, free gas saturation, and gas hydrate saturation of the A-B ore body could be well predicted based on theP-wave impedance and P-S wave velocity ratio.


## 4. Results

### 4.1. Logging Interpretation of Gas Hydrate and Free Gas

Applying the overlap method to the A-well and B-well, it can be found that the resistivity curve of the B-well is characterized by high value in the depth of 1460–1525 mbsl but medium value in the depth of 1525–1542 mbsl (Figure5(a)). Its interval transit time, density porosity, and neutron porosity curves overlap each other in the depth of 1460–1510 mbsl but separate obviously in the depth of 1510–1542 mbsl. Thus, the depth range of 1460–1510 mbsl in the B-well presents the response characteristics of gas hydrate, i.e., high resistivity and the overlap of interval transit time, density porosity, and neutron porosity. In contrast, in the 1510–1542 mbsl of the B-well, the response characteristics of free gas are characterized by medium-high resistivity, the separation of neutron porosity, and interval transit time with density porosity.Figure 5
Comprehensive log interpretation graphs of the A-well and B-well.
(a)(b)In addition, the depth range of 1425–1510 mbsl in the A-well shows the response characteristics of gas hydrate, i.e., high resistivity and the overlap of interval transit time, density porosity, and neutron porosity (Figure5(b)). However, no obvious free gas is identified on the well logs of the A-well, perhaps due to limited logging depth.


### 4.2. Seismic Reflection Characteristics of Gas Hydrate and Free Gas

According to logging curves and synthetic seismic records, it could be found that gas hydrate shows high P-wave impedance, and there is a positive reflection coefficient interface between the top of gas hydrate and the overlying strata with lowP-wave impedance. In the seismic section, the corresponding seismic wave at the top of gas hydrate is characterized by strong wave peak, which indicates that the seismic section belongs to a normal polarity section. If the bottom of gas hydrate is in direct contact with the underlying free gas with low P-wave impedance, a negative reflection coefficient interface will appear, and the corresponding seismic wave in the seismic section will show a strong trough. On the contrary, if the bottom of gas hydrate does not contact with free gas directly, the reflection intensity of seismic wave is weak.The top of gas hydrate in the B-well is characterized by strong wave peak, whereas the bottom shows strong wave trough, indicating that gas hydrate is in direct contact with free gas (Figure6(a)). In addition, the top of gas hydrate in the A-well is characterized by strong wave peak, but the seismic wave amplitude at the bottom of gas hydrate is relatively weak (Figure 6(b)). It is shown that gas hydrate does not contact with free gas directly in the A-well, so a transition zone may exist between gas hydrate and free gas. Gas hydrates, free gas, and water coexist in the transition zone. The gas hydrates in this area are in a critically stable state and thus will dissociate once the temperature slightly increases or the pressure slightly decreases [39].Figure 6
Logging curves and synthetic seismic records of the A-well and B-well.
(a)(b)


### 4.3. AVO Attribute Analysis of Gas Hydrate and Free Gas

According to well-to-seismic calibration, there are obvious anomaly areas in both intercept and gradient cross sections, showing a feature of bright points (Figure7). In view of this feature, root-mean-square attribute values of intercept and gradient are extracted. The high-value areas represent the favorable distribution zones of gas hydrate and free gas, which are nearly elliptical in the NW-SE direction (Figure 8).Figure 7
Well-tie AVO intercept and gradient cross section of the A-well and B-well: (a) AVO intercept; (b) AVO gradient.
(a)(b)Figure 8
Root-mean-square attribute distribution of AVO intercept and gradient in the A-B ore body: (a) root-mean-square attribute of AVO intercept for gas hydrate; (b) root-mean-square attribute of AVO intercept for free gas; (c) root-mean-square attribute of AVO gradient for gas hydrate; (d) root-mean-square attribute of AVO gradient for free gas.
(a)(b)(c)(d)


### 4.4. Seismic Inversion of Gas Hydrate and Free Gas

#### 4.4.1.P-wave Velocity

The well-tieP-wave velocity cross section of the A-well and B-well shows that gas hydrate possesses a high velocity of 1800–2000 m/s (Figure 9). On the contrary, the underlying free gas has abnormal low velocity of 1400–1600 m/s.Figure 9
Well-tieP-wave velocity cross section of the A-well and B-well (according to [3]).


#### 4.4.2.P-wave Impedance

According to seismic inversion, it could be found that gas hydrate is located above free gas and distributed continuously (Figure10). The underlying free gas is especially enriched at the boundary of gas hydrate. The inversion results are in good agreement with the actual logging data, which just proves the previous interpretation of logging data.Figure 10
Well-tieP-wave impedance cross section of the A-well and B-well.


#### 4.4.3. Thickness and Spatial Distribution

Based on the distribution characteristics ofP-wave impedance, gas hydrate and free gas are three-dimensionally carved to obtain the spatial distribution of the geologic body and calculate its thickness and volume. In plane form, gas hydrate shows nearly elliptical distribution in the NW-SE direction, and its thickness gradually thinned from the middle to the edge (Figure 11). Free gas has similar distribution characteristics with gas hydrate; that is, the distribution range, morphological feature, and thickness change are basically consistent with gas hydrate (Figure 12).Figure 11
Thickness and spatial distribution of gas hydrate in the A-B ore body: (a) thickness distribution; (b) top surface structure; (c) bottom surface structure.
(a)(b)(c)Figure 12
Thickness and spatial distribution of free gas in the A-B ore body: (a) thickness distribution; (b) top surface structure; (c) bottom surface structure.
(a)(b)(c)In space, three-dimensional carving technology is used to find out separate free gas reservoirs and quantitatively describe the surface structure, thickness, and volume of each free gas reservoir, so as to determine the most favorable free gas reservoir for exploitation. It is found that among all the free gas reservoirs, the free gas reservoir located directly below gas hydrate possesses the largest development scale, which is more than 40 times larger than other free gas reservoirs (Figure13).Figure 13
Three-dimensional spatial distribution of the main free gas reservoirs in the A-B ore body.
(a)(b)(c)(d)


#### 4.4.4. Effective Porosity

The effective porosity distribution of the A-B ore body changes with P-wave impedance. In the section, the effective porosity of gas hydrate is mainly between 0.30 and 0.40, whereas the effective porosity of free gas is higher than that of gas hydrate, ranging from 0.35 to 0.40 and with small change and weak heterogeneity (Figure14). These results are consistent with those predicted by logging interpretation (Figure 5).Figure 14
Well-tie effective porosity cross section of the A-well and B-well.In the plane, the effective porosity of gas hydrate is between 0.31 and 0.35 and increases from the middle to the surrounding, so it has obvious high value at the boundary of the A-B ore body (Figure15(a)). In addition, the effective porosity of the B-well is higher than that of the A-well. On the whole, the effective porosity of free gas has little change in plane distribution, and its value is between 0.34 and 0.37 (Figure 15(b)). Different from gas hydrate, the effective porosity of free gas in the northern A-well area is slightly better than that in the southern B-well area.Figure 15
Plane distribution of effective porosity for gas hydrate and free gas in the A-B ore body.
(a)(b)


#### 4.4.5. Saturation

The saturation of the A-B ore body varies with theP-wave-to-S-wave velocity ratio. In the longitudinal section, the saturation of gas hydrate shows a trend of low value in the upper part and high value in the lower part, mainly in the range of 0.30–0.50, while the free gas saturation is generally low, as low as 0.20 (Figure 16).Figure 16
Well-tie saturation cross section of the A-well and B-well.In the plane, the average saturation of gas hydrate is distributed in the range of 0.33–0.40, with a small variation range (Figure17(a)). The value of gas hydrate saturation in the southern B-well area is slightly larger than that in the northern A-well area. On the whole, the saturation of free gas is smaller than that of gas hydrate, which primarily ranges from 0.24 to 0.32 with little change (Figure 17(b)). There is a NW-trending low saturation zone through the B-well, which divides the A-B ore body into two distinct free gas area.Figure 17
Plane distribution of average saturation for gas hydrate and free gas in the A-B ore body.
(a)(b)


#### 4.4.6. Average Effective Thickness

Based on the above study, the average effective thicknesses of gas hydrate and free gas could be calculated according to the thickness, effective porosity, and saturation. It can be seen in the plane that the average effective thicknesses are consistent with the trend of gas hydrate reservoir thicknesses (Figure18). The average effective thicknesses in the north are better than those in the south, and the low values are primarily distributed at the boundary of the A-B ore body. The average effective thickness of gas hydrate primarily varies from 3.0 m to 10.5 m (Figure 18(a)), and that of free gas is mainly between 2.0 m and 5.0 m (Figure 18(b)).Figure 18
Plane distribution of average effective thickness for gas hydrate and free gas in the A-B ore body.
(a)(b)


## 4.1. Logging Interpretation of Gas Hydrate and Free Gas

Applying the overlap method to the A-well and B-well, it can be found that the resistivity curve of the B-well is characterized by high value in the depth of 1460–1525 mbsl but medium value in the depth of 1525–1542 mbsl (Figure5(a)). Its interval transit time, density porosity, and neutron porosity curves overlap each other in the depth of 1460–1510 mbsl but separate obviously in the depth of 1510–1542 mbsl. Thus, the depth range of 1460–1510 mbsl in the B-well presents the response characteristics of gas hydrate, i.e., high resistivity and the overlap of interval transit time, density porosity, and neutron porosity. In contrast, in the 1510–1542 mbsl of the B-well, the response characteristics of free gas are characterized by medium-high resistivity, the separation of neutron porosity, and interval transit time with density porosity.Figure 5
Comprehensive log interpretation graphs of the A-well and B-well.
(a)(b)In addition, the depth range of 1425–1510 mbsl in the A-well shows the response characteristics of gas hydrate, i.e., high resistivity and the overlap of interval transit time, density porosity, and neutron porosity (Figure5(b)). However, no obvious free gas is identified on the well logs of the A-well, perhaps due to limited logging depth.


## 4.2. Seismic Reflection Characteristics of Gas Hydrate and Free Gas

According to logging curves and synthetic seismic records, it could be found that gas hydrate shows high P-wave impedance, and there is a positive reflection coefficient interface between the top of gas hydrate and the overlying strata with lowP-wave impedance. In the seismic section, the corresponding seismic wave at the top of gas hydrate is characterized by strong wave peak, which indicates that the seismic section belongs to a normal polarity section. If the bottom of gas hydrate is in direct contact with the underlying free gas with low P-wave impedance, a negative reflection coefficient interface will appear, and the corresponding seismic wave in the seismic section will show a strong trough. On the contrary, if the bottom of gas hydrate does not contact with free gas directly, the reflection intensity of seismic wave is weak.The top of gas hydrate in the B-well is characterized by strong wave peak, whereas the bottom shows strong wave trough, indicating that gas hydrate is in direct contact with free gas (Figure6(a)). In addition, the top of gas hydrate in the A-well is characterized by strong wave peak, but the seismic wave amplitude at the bottom of gas hydrate is relatively weak (Figure 6(b)). It is shown that gas hydrate does not contact with free gas directly in the A-well, so a transition zone may exist between gas hydrate and free gas. Gas hydrates, free gas, and water coexist in the transition zone. The gas hydrates in this area are in a critically stable state and thus will dissociate once the temperature slightly increases or the pressure slightly decreases [39].Figure 6
Logging curves and synthetic seismic records of the A-well and B-well.
(a)(b)


## 4.3. AVO Attribute Analysis of Gas Hydrate and Free Gas

According to well-to-seismic calibration, there are obvious anomaly areas in both intercept and gradient cross sections, showing a feature of bright points (Figure7). In view of this feature, root-mean-square attribute values of intercept and gradient are extracted. The high-value areas represent the favorable distribution zones of gas hydrate and free gas, which are nearly elliptical in the NW-SE direction (Figure 8).Figure 7
Well-tie AVO intercept and gradient cross section of the A-well and B-well: (a) AVO intercept; (b) AVO gradient.
(a)(b)Figure 8
Root-mean-square attribute distribution of AVO intercept and gradient in the A-B ore body: (a) root-mean-square attribute of AVO intercept for gas hydrate; (b) root-mean-square attribute of AVO intercept for free gas; (c) root-mean-square attribute of AVO gradient for gas hydrate; (d) root-mean-square attribute of AVO gradient for free gas.
(a)(b)(c)(d)


## 4.4. Seismic Inversion of Gas Hydrate and Free Gas

### 4.4.1.P-wave Velocity

The well-tieP-wave velocity cross section of the A-well and B-well shows that gas hydrate possesses a high velocity of 1800–2000 m/s (Figure 9). On the contrary, the underlying free gas has abnormal low velocity of 1400–1600 m/s.Figure 9
Well-tieP-wave velocity cross section of the A-well and B-well (according to [3]).


### 4.4.2.P-wave Impedance

According to seismic inversion, it could be found that gas hydrate is located above free gas and distributed continuously (Figure10). The underlying free gas is especially enriched at the boundary of gas hydrate. The inversion results are in good agreement with the actual logging data, which just proves the previous interpretation of logging data.Figure 10
Well-tieP-wave impedance cross section of the A-well and B-well.


### 4.4.3. Thickness and Spatial Distribution

Based on the distribution characteristics ofP-wave impedance, gas hydrate and free gas are three-dimensionally carved to obtain the spatial distribution of the geologic body and calculate its thickness and volume. In plane form, gas hydrate shows nearly elliptical distribution in the NW-SE direction, and its thickness gradually thinned from the middle to the edge (Figure 11). Free gas has similar distribution characteristics with gas hydrate; that is, the distribution range, morphological feature, and thickness change are basically consistent with gas hydrate (Figure 12).Figure 11
Thickness and spatial distribution of gas hydrate in the A-B ore body: (a) thickness distribution; (b) top surface structure; (c) bottom surface structure.
(a)(b)(c)Figure 12
Thickness and spatial distribution of free gas in the A-B ore body: (a) thickness distribution; (b) top surface structure; (c) bottom surface structure.
(a)(b)(c)In space, three-dimensional carving technology is used to find out separate free gas reservoirs and quantitatively describe the surface structure, thickness, and volume of each free gas reservoir, so as to determine the most favorable free gas reservoir for exploitation. It is found that among all the free gas reservoirs, the free gas reservoir located directly below gas hydrate possesses the largest development scale, which is more than 40 times larger than other free gas reservoirs (Figure13).Figure 13
Three-dimensional spatial distribution of the main free gas reservoirs in the A-B ore body.
(a)(b)(c)(d)


### 4.4.4. Effective Porosity

The effective porosity distribution of the A-B ore body changes with P-wave impedance. In the section, the effective porosity of gas hydrate is mainly between 0.30 and 0.40, whereas the effective porosity of free gas is higher than that of gas hydrate, ranging from 0.35 to 0.40 and with small change and weak heterogeneity (Figure14). These results are consistent with those predicted by logging interpretation (Figure 5).Figure 14
Well-tie effective porosity cross section of the A-well and B-well.In the plane, the effective porosity of gas hydrate is between 0.31 and 0.35 and increases from the middle to the surrounding, so it has obvious high value at the boundary of the A-B ore body (Figure15(a)). In addition, the effective porosity of the B-well is higher than that of the A-well. On the whole, the effective porosity of free gas has little change in plane distribution, and its value is between 0.34 and 0.37 (Figure 15(b)). Different from gas hydrate, the effective porosity of free gas in the northern A-well area is slightly better than that in the southern B-well area.Figure 15
Plane distribution of effective porosity for gas hydrate and free gas in the A-B ore body.
(a)(b)


### 4.4.5. Saturation

The saturation of the A-B ore body varies with theP-wave-to-S-wave velocity ratio. In the longitudinal section, the saturation of gas hydrate shows a trend of low value in the upper part and high value in the lower part, mainly in the range of 0.30–0.50, while the free gas saturation is generally low, as low as 0.20 (Figure 16).Figure 16
Well-tie saturation cross section of the A-well and B-well.In the plane, the average saturation of gas hydrate is distributed in the range of 0.33–0.40, with a small variation range (Figure17(a)). The value of gas hydrate saturation in the southern B-well area is slightly larger than that in the northern A-well area. On the whole, the saturation of free gas is smaller than that of gas hydrate, which primarily ranges from 0.24 to 0.32 with little change (Figure 17(b)). There is a NW-trending low saturation zone through the B-well, which divides the A-B ore body into two distinct free gas area.Figure 17
Plane distribution of average saturation for gas hydrate and free gas in the A-B ore body.
(a)(b)


### 4.4.6. Average Effective Thickness

Based on the above study, the average effective thicknesses of gas hydrate and free gas could be calculated according to the thickness, effective porosity, and saturation. It can be seen in the plane that the average effective thicknesses are consistent with the trend of gas hydrate reservoir thicknesses (Figure18). The average effective thicknesses in the north are better than those in the south, and the low values are primarily distributed at the boundary of the A-B ore body. The average effective thickness of gas hydrate primarily varies from 3.0 m to 10.5 m (Figure 18(a)), and that of free gas is mainly between 2.0 m and 5.0 m (Figure 18(b)).Figure 18
Plane distribution of average effective thickness for gas hydrate and free gas in the A-B ore body.
(a)(b)


## 4.4.1.P-wave Velocity

The well-tieP-wave velocity cross section of the A-well and B-well shows that gas hydrate possesses a high velocity of 1800–2000 m/s (Figure 9). On the contrary, the underlying free gas has abnormal low velocity of 1400–1600 m/s.Figure 9
Well-tieP-wave velocity cross section of the A-well and B-well (according to [3]).


## 4.4.2.P-wave Impedance

According to seismic inversion, it could be found that gas hydrate is located above free gas and distributed continuously (Figure10). The underlying free gas is especially enriched at the boundary of gas hydrate. The inversion results are in good agreement with the actual logging data, which just proves the previous interpretation of logging data.Figure 10
Well-tieP-wave impedance cross section of the A-well and B-well.


## 4.4.3. Thickness and Spatial Distribution

Based on the distribution characteristics ofP-wave impedance, gas hydrate and free gas are three-dimensionally carved to obtain the spatial distribution of the geologic body and calculate its thickness and volume. In plane form, gas hydrate shows nearly elliptical distribution in the NW-SE direction, and its thickness gradually thinned from the middle to the edge (Figure 11). Free gas has similar distribution characteristics with gas hydrate; that is, the distribution range, morphological feature, and thickness change are basically consistent with gas hydrate (Figure 12).Figure 11
Thickness and spatial distribution of gas hydrate in the A-B ore body: (a) thickness distribution; (b) top surface structure; (c) bottom surface structure.
(a)(b)(c)Figure 12
Thickness and spatial distribution of free gas in the A-B ore body: (a) thickness distribution; (b) top surface structure; (c) bottom surface structure.
(a)(b)(c)In space, three-dimensional carving technology is used to find out separate free gas reservoirs and quantitatively describe the surface structure, thickness, and volume of each free gas reservoir, so as to determine the most favorable free gas reservoir for exploitation. It is found that among all the free gas reservoirs, the free gas reservoir located directly below gas hydrate possesses the largest development scale, which is more than 40 times larger than other free gas reservoirs (Figure13).Figure 13
Three-dimensional spatial distribution of the main free gas reservoirs in the A-B ore body.
(a)(b)(c)(d)


## 4.4.4. Effective Porosity

The effective porosity distribution of the A-B ore body changes with P-wave impedance. In the section, the effective porosity of gas hydrate is mainly between 0.30 and 0.40, whereas the effective porosity of free gas is higher than that of gas hydrate, ranging from 0.35 to 0.40 and with small change and weak heterogeneity (Figure14). These results are consistent with those predicted by logging interpretation (Figure 5).Figure 14
Well-tie effective porosity cross section of the A-well and B-well.In the plane, the effective porosity of gas hydrate is between 0.31 and 0.35 and increases from the middle to the surrounding, so it has obvious high value at the boundary of the A-B ore body (Figure15(a)). In addition, the effective porosity of the B-well is higher than that of the A-well. On the whole, the effective porosity of free gas has little change in plane distribution, and its value is between 0.34 and 0.37 (Figure 15(b)). Different from gas hydrate, the effective porosity of free gas in the northern A-well area is slightly better than that in the southern B-well area.Figure 15
Plane distribution of effective porosity for gas hydrate and free gas in the A-B ore body.
(a)(b)


## 4.4.5. Saturation

The saturation of the A-B ore body varies with theP-wave-to-S-wave velocity ratio. In the longitudinal section, the saturation of gas hydrate shows a trend of low value in the upper part and high value in the lower part, mainly in the range of 0.30–0.50, while the free gas saturation is generally low, as low as 0.20 (Figure 16).Figure 16
Well-tie saturation cross section of the A-well and B-well.In the plane, the average saturation of gas hydrate is distributed in the range of 0.33–0.40, with a small variation range (Figure17(a)). The value of gas hydrate saturation in the southern B-well area is slightly larger than that in the northern A-well area. On the whole, the saturation of free gas is smaller than that of gas hydrate, which primarily ranges from 0.24 to 0.32 with little change (Figure 17(b)). There is a NW-trending low saturation zone through the B-well, which divides the A-B ore body into two distinct free gas area.Figure 17
Plane distribution of average saturation for gas hydrate and free gas in the A-B ore body.
(a)(b)


## 4.4.6. Average Effective Thickness

Based on the above study, the average effective thicknesses of gas hydrate and free gas could be calculated according to the thickness, effective porosity, and saturation. It can be seen in the plane that the average effective thicknesses are consistent with the trend of gas hydrate reservoir thicknesses (Figure18). The average effective thicknesses in the north are better than those in the south, and the low values are primarily distributed at the boundary of the A-B ore body. The average effective thickness of gas hydrate primarily varies from 3.0 m to 10.5 m (Figure 18(a)), and that of free gas is mainly between 2.0 m and 5.0 m (Figure 18(b)).Figure 18
Plane distribution of average effective thickness for gas hydrate and free gas in the A-B ore body.
(a)(b)


## 5. Conclusions

Through evaluating the gas hydrate and free gas system of the A-B ore body in the Shenhu area, the following conclusions are obtained.(1)
According to logging interpretation, gas hydrate of the B-well is located in the depth range of 1460–1510 mbsl and free gas is in 1510–1542 mbsl. Moreover, gas hydrate of the A-well is located in the depth range of 1425–1512 mbsl, and no obvious free gas is identified, due to limited logging depth(2)
Through the intersection analysis of elastic parameters, gas hydrate has the characteristics of highP-wave impedance, while free gas is characterized by low P-wave impedance. P-wave impedance is one of the most sensitive elastic parameters for gas hydrate and free gas(3)
Based on the correlation analysis of physical and elastic parameters, the effective porosity has a good correlation withP-wave impedance and so has saturation with the P-wave-to-S-wave velocity ratio. Therefore, according to the ratio of P-wave impedance and P-wave-to-S-wave velocity ratio, the effective porosity and saturation of hydrate and free gas can be calculated, finally predicting the spatial distribution of physical parameters(4)
According to logging curves and seismic synthetic record, gas hydrate in the B-well is in direct contact with free gas, while a transition zone may exist between gas hydrate and free gas in the A-well. According to well-to-seismic calibration, the anomaly areas of root-mean-square attribute values in intercept and gradient cross sections represent the favorable distribution zones of gas hydrate and free gas(5)
Gas hydrate has an effective porosity of 0.30–0.40, an average saturation of 0.33–0.40, and an effective thickness of 3.0–10.5 m, whereas free gas possesses an effective porosity of 0.35–0.40, a saturation of 0.24–0.32, and an effective thickness of 2.0–5.0 m



---

*Source: 5514263-2021-10-05.xml*