1 Introduction
China has abundant shale gas resources. Exploration and development of shale gas has become a hot spot due to its wide distribution, large resources, and long stable production period (
Hao et al., 2013;
Liang et al., 2014;
Zou et al., 2016). In recent years, China’s shale gas production has grown rapidly, and remarkable accomplishments have been achieved in several marine shale gas blocks in the southern Sichuan Basin, including the blocks of Fuling, Weiyuan, Changning, Zhaotong and Fushun-Yongchuan, etc. (
Guo et al., 2017;
Dong et al., 2018;
Jin et al., 2018;
Jiao, 2019;
Yi et al., 2019;
Ma et al., 2020). Shale gas reservoirs, one of the unconventional natural gas reservoirs, commonly contain connate water (i.e., pore water), and has a coexisting characteristic of gas and water under
in situ condition (
Ahmad and Haghighi, 2013;
Li et al., 2016;
Li et al., 2019a). A large amount of connate water is bound in shale pores, with a water saturation of 10%−90%, which affects the accumulation and seepage of adsorbed gas and free gas in shale reservoirs (
Cheng et al., 2017;
Li et al., 2019b). Field practice shows that shale reservoirs with a higher connate water saturation have a lower gas content and lower single-well production (
Fang et al., 2014;
Gasparik et al., 2014;
Xu et al., 2022).
As the material basis of shale gas production,
in situ contents of adsorbed and free gases are significantly affected by the connate water in shale reservoirs (
Dehghanpour et al., 2013;
Gasparik et al., 2014;
Xu et al., 2022). The influence of water content on the adsorption capacity of shale gas is mainly manifested in two aspects: a) water molecules occupy some of the effective gas adsorption sites on the pore surface, reducing the amount of gas adsorption (
Chalmers and Bustin, 2007); b) water molecules obstruct some small pore-throats, reducing the gas storage space (
Tang et al., 2017). The impact of water on free gas content is mainly to reduce the gas storage space. Previous scholars have reported that the adsorption capacity of methane in argillaceous shale under the condition of equilibrium water is reduced by 65%−90%, and the adsorption capacity in organic-rich shale is reduced by 20%−60% compared with the dry condition (Ross and Marc Bustin, 2009;
Gasparik et al., 2014). Moisture has a great negative impact on the contents of adsorbed gas and gases in shales.
The migration process of shale gas is also under the restrictions of pore water (
Pan et al., 2010;
Yuan et al., 2014). The transport mechanism of adsorbed gas is surface diffusion, which is affected by pressure gradient and gas adsorption concentration (
Yuan et al., 2014;
Yang et al., 2016). The transport mechanism of free gas is a continuous flow, gas slippage flow, Knudsen diffusion, and is also affected by factors such as effective stress, matrix shrinkage, and gas desorption during the development process (
Wu et al., 2016). The surface diffusivity of adsorbed gas will be affected by the existence of adsorbed water in shale, and the extent of the influence is determined by the amount of adsorbed water (
Yuan et al., 2014). In the process of free gas migration, the free path of gas molecules is close to the characteristic scale of pores, and nanoscale effects (such as gas slippage, Knudsen diffusion, etc.) are gradually significant (
Loucks et al., 2009;
Jiang et al., 2020;
Sheng et al., 2020;
Yang et al., 2020). For water-bearing shale, when the pore water occupies part of the space, the effective migration channel for gas flow becomes narrower. The pore-throat radius and the thickness of the adsorbed water layer, particularly in the situation of low flow rate (
Jiang et al., 2020), are in the same order of magnitude, and the adsorbed water becomes the main factor affecting the free gas seepage (
Li et al., 2017a;
Sang et al., 2018,
2020). The diffusion coefficient of methane in shale water-containing conditions is significantly lower than that in dry conditions, indicating that water molecules hinder the migration of methane in shale (
Yuan et al., 2014). Wu et al. (2014) analyzed the gas-water two-phase flow in nano-pores (100 nm) and showed that when the pore water saturation is 20%, the gas-phase flow capacity is reduced by about 10% compared with the dry condition; In the case of 40%, the gas-phase flow capacity is reduced by about 20%. The permeability of gas generally decreases with increasing water saturation (
Zhang et al., 2017;
Li et al., 2019b;
Wang et al., 2019), therefore the pore water of shale has a significant effect on the migration of shale gas.
Water-bearing characteristics of shales are mainly determined by macro parameters such as initial water saturation, irreducible water saturation, and movable water saturation (
Yuan et al., 2018;
Yang et al., 2021). Considering the different impacts of pore water (adsorbed water and free water) on shale gas in different states, it is significant to quantitatively estimate the amounts of adsorbed and free water in shale matrix pores. Additionally, in recent years there has been a growing amount of interest in deep shale gas with a burial depth of 3500–4500 m (
Jiang et al., 2017;
Long et al., 2018;
Zhang et al., 2022). Pore water in deep shales was quantitatively evaluated in this study. The amounts of adsorbed and free water in saturated shales were obtained based on the saturation-centrifugation experiments. Further, the impact of shale material composition and microstructure on the pore water was analyzed. The research provides a new idea for the quantitative characterization of shale matrix pore water, which is significant to the exploration and extraction of shale gas.
2 Experiments and methodologies
2.1 Experimental section
2.1.1 Saturation and centrifugation
Six samples, with a range of burial depths of 3978.13−4047.04 m, were obtained from the deep marine shale reservoirs of the Wufeng-Longmaxi Formation at Well LU206. The well is located in the Luzhou block from the southern Sichuan Basin, China (
Zhu et al., 2021). First, the core was prepared into a core column with a 2.5 cm diameter. Secondly, each core column was put into a DZF vacuum drying oven and was dried under a vacuum and 110°C environment for at least 24 h to ensure the residual water in the core is completely removed. After taking it out, the mass of the dry core is tested as
m1. Thirdly, the dry core was then put into the ZYB-II saturation instrument and then the core was saturated with distilled water. During the saturation process, 10 MPa fluid pressure and −0.1 MPa vacuum pressure were applied to the experimental core for 12 days. Then, the core was taken out and weighed as
m2. Next, the saturated water core was put into the GL-21M refrigerated centrifuge to perform a centrifugation test by setting different centrifugal rotation speeds (3000−10000 r/min, interval 1000 r/min) under the experimental condition of 30°C. Centrifugation has a duration of at least 4 h at each centrifugal rotation speed. Then, the mass was recorded as
mx to determine the movable water amount.
The total water amount can be obtained by
where Qt is the total water amount in shale, mg/g; m1 is the mass of the dry core, g; m2 is the mass of the saturated core, mg/g.
And movable water amount (Qm, mg/g) under each centrifugation condition can be obtained by
where Qm is the movable water amount under a certain centrifugal force, mg/g; mx is the mass of core after centrifugation balance, g.
2.1.2 other experiments
Petrology, geochemistry, microstructure, and other analysis tests were also carried out to analyze the influencing factors of adsorbed and free water amounts, including total organic carbon (TOC) content test, X-ray diffraction (XRD) whole-rock mineral analysis, low-temperature nitrogen adsorption, and scanning electron microscope analysis.
2.2 Evaluation method
The main states of water in the pores of the shale matrix are adsorbed and free, with free water having the potential to flow and adsorbed water keeping a nearly immobile situation (
Li et al., 2019a). However, due to the influence of the pore-throat microstructure of shale, not all free water can flow effectively under certain conditions. Under external conditions (such as centrifugal force), some free water can flow, which is movable water; the part of free water that cannot flow is capillary-bound water. With the change in external conditions, capillary-bound water can be converted into movable water. When all capillary-bound water is converted into movable water, the amount of movable water reaches the maximum value and is equal to the amount of free water (
Li et al., 2019a).
Theoretically, as long as the centrifugal force is continuously increased, the capillary-bound water will be continuously converted into movable water, to obtain the maximum movable water amount. However, it is unrealistic to convert all capillary-bound water into movable water by infinitely increasing the centrifugal force of the centrifuge under the current experimental condition. Since the amount of free water cannot be directly obtained through experiments, it is extremely important to find the functional relationship between centrifugal force and movable water amount. In this way, the maximum movable water amount, that is free water amount, can be theoretically calculated through a series of centrifugal that can be realized in the laboratory.
The relationship between movable water amount and centrifugal force is as follows (
Li et al., 2019a):
where Qf is the maximum movable water amount when the centrifugal force reaches infinity, that is, the free water amount, mg/g; ΔP is the centrifugal pressure difference (that is, the centrifugal force), MPa; ΔPL is the median pressure difference, that is the centrifugal force when the movable water volume reaches half of the free water amount, MPa.
To facilitate the analysis of parameters in Eq. (3), the equation can be converted into a linear function of 1/Qm and 1/ΔP:
By fitting the experimental data and establishing a linear relationship between 1/Qm and 1/ΔP, the free water amount can be obtained.
Further, the amount of adsorbed water is equal to the difference between saturated water amount and free water amount and can be written as
According to the amounts of adsorbed water and free water, the adsorption ratio can also be calculated as follows:
where ra is the ratio of adsorbed water to total water, %.
3 Results and discussion
3.1 Adsorbed and free amounts of pore water
The results of the saturation-centrifugation experiment are shown in Fig.1(a) and Tab.1, showing the movable water amount under different centrifugal force conditions. It can be seen that the movable water volume gradually increases with the increase in centrifugal force, but the growth rate can be seen to slow down significantly. Thus, it can be inferred that when the centrifugal force increases to infinity, the movable water amount tends to a constant value, this value is the amount of free water. According to the relationship between the reciprocal of movable water amount and the reciprocal of centrifugal force (Fig.1(b)) it can be noted that there is a good linear correlation (R2 > 0.99). The linear relationship can be used to estimate the maximum movable water amount, that is, the free water amount. For the studied shales under the experimental conditions (30°C, distilled water), the amount of free water varies from 9.5511 to 19.802 mg/g (mean 13.9541 mg/g). Further, the amount of adsorbed water can be calculated by subtracting the amount of free water from the total water amount and varies from 1.7967 to 9.8218 mg/g (mean 6.4501 mg/g). Comparatively, the amount of adsorbed water is larger than free water and accounts for 15.83%–42.61% (mean 30.45%) of total water (Tab.2). This demonstrates that part of the water stored in the matrix pores can potentially flow, and will have a certain contribution to the flowback fluid of shale gas well.
Statistics show that the total water amount in shales has a good positive correlation with both adsorbed water and free water amounts (Fig.2), which implies the coexistence of adsorbed water and free water in the pores of the shale matrix. Vertically, the free water and total water amounts in Well LU206 show an overall trend of decrease, and the change in the amount of adsorbed water is small with the increase in burial depth (Fig.3). The Wufeng Formation shale has the lowest water-bearing amount and adsorption ratio.
3.2 Influencing factors of pore water
3.2.1 Coupling relationship between shale compositions and pore microstructure
The occurrence of pore water has an internal relationship with the material composition and pore-throat system (reservoir space) of shale (
Hu et al., 2015;
Wang et al., 2019;
Tian et al., 2020). Various organic/inorganic minerals constitute the rock skeleton, in which a pore-throat system develops, and pore water occurs in the pore-throat. Therefore, in order to reveal the microscopic occurrence mechanism of pore water, it is necessary to clarify the coupling relationship between shale material composition and pore-throat microstructure.
The inorganic mineral components of the shale in the study area are mainly quartz and clay minerals, accounting for 43.7%−79.3% (mean 65.5%) (Fig.4(a)). The quartz contents of the shales range from 31.5% to 61.7% (mean 43.9%); the contents of clay minerals range from 11.6% to 36.1%, with an average of 21.6%. Shales develop ferrodolomite (mean 16.5%) and calcite (mean 8.9%) minerals, but the ferrodolomite content in the L206-1 shale sample is abnormally high with a content of 47.6%. In addition, shales contain a small amount of plagioclase and pyrite. Pyrite is an important geochemical indicator mineral, which can indicate the formation environment of organic matter. The pyrite content in the study area is 1.1%−5.6%, with an average of 3.9%. With the increase in burial depth, the content of clay minerals gradually decreased, and the quartz minerals showed a general trend of increase.
Through scanning electron microscope observation, more mineral components can be seen. In addition to the main minerals such as clay, quartz, calcite, ferrodolomite, and pyrite, a small amount of albite, fluorapatite, and titanium dioxide can also be seen (Fig.5). Crystals of gypsum were even found in localized shale samples. Under the microscope, it can be observed that the intragranular pores of clay particles are generally developed, and some calcites develop intercrystalline pores. In addition, the intercrystalline pores in the gypsum crystals are very developed. However, pores in clay minerals provide major pore space (Fig.5(a)).
The TOC content of the studied shales is 2.21%−4.34%, with an average of 3.21% (Fig.4(b)). Through the identification under the microscope, the organic matter of shale is mainly composed of solid bitumen, graptolite, and vitrinite-like, of which solid bitumen is the main organic micro-component, accounting for 48%−80% (mean 64.8%) (Tab.3). Graptolite accounts for 15%−46%, with an average of 31.6%. Vertically, the bottom of the Longmaxi Formation of Well LU206 has higher TOC content and solid bitumen ratio, indicating favorable exploration and development intervals. Under field emission scanning electron microscope (Fig.6), the solid bitumen has an irregular shape with clear edges and corners and a large number of nano-pores developed in it. These pores may be the storage space for pore water.
The Luzhou block is located in the center of deep-water continental shelf deposition, and the bottom of the Longmaxi Formation is rich in biosiliceous minerals. The statistical results show that the shale TOC content has a good positive correlation with the quartz mineral content (Fig.7(a)), showing the characteristics of biogenic silicon, but the pores related to the quartz mineral are not developed (Fig.5(b)) with a weak contribution to water storage. In addition, there is a certain positive correlation between TOC content and pyrite content (Fig.7(b)). On the one hand, pyrite indicates a reducing environment, which is conducive to the enrichment of organic matter; on the other hand, pyrite is closely related to hydrocarbon generation (release of sulfur) from organic matter (
Berner, 1984;
Chen et al., 2021). Thus, organic matter (especially solid bitumen) is usually in a symbiotic state with pyrite (Fig.5(f)).
Based on low-temperature nitrogen adsorption/desorption experiments (
Li et al., 2017b,
2018), the microstructure of shale was determined, including pore volume, specific surface area, pore size, etc. As shown in Fig.8, shale pores are dominated by mesopores (2−50 nm), accounting for 73%−78%, with an average of 75%. The statistical results show that the pore volume of shale is positively correlated with the specific surface area (Fig.9(a)). It may be that the more pores, the higher the pore volume and the specific surface area when the pore size is close. In addition, the higher the proportion of mesopores, the higher the pore volume and the larger the specific surface area (Fig.9(b)−9(c)).
Finally, the relationship between shale material composition and microstructure is summarized. The positive correlations of TOC and pyrite contents with specific surface area reflect that shale organic pores (mainly developed in solid bitumen) contribute to shale-specific surface area (Fig.10(a) and 10(b)). Clay and pyrite contents are positively correlated with pore volume (Fig.10(c) and 10(d)). The correlation reflects that the shale clay mineral intragranular pores and solid bitumen intragranular pores comprehensively contribute to the shale pore volume.
3.2.2 Controls of shale compositions and pore microstructure on pore water occurrence
The adsorbed water amount in the pores of the shale matrix is mainly controlled by the specific surface area of shale, and the free water amount is mainly controlled by shale pore volume (Fig.11). It can also be found that although the amount of adsorbed water is positively correlated with the specific surface area, and the amount of free water is positively correlated with the pore volume, the regression line has an intercept, indicating that some of the pores determined by low-temperature nitrogen have no adsorbed water or free water. A portion of voids do not contain water. Therefore, it is speculated that there may be a lower limit for the occurrence of pore water in the shale matrix. The adsorbed water in the study area mainly occurs on the pore surface, and the lower limit of the specific surface area is approximately 9 m2/g; the free water mainly occurs inside the pores, and the lower limit of the pore volume is approximately 5 × 10−3 cm3/g. After the lower limits of specific surface area and pore volume are determined, the corresponding lower limit values can be subtracted to calculate the specific surface area and pore volume of water-bearing pores in the shale matrix, which lays the foundation for the subsequent evaluation of adsorption parameters. In addition, the lower limit of shale in different regions may be different, which is closely related to the regional depositional environment, shale pore development characteristics, and pore surface properties.
Shale compositions (organic matter and inorganic matter) are complex and the hydrophilic capacity of each component is different (
Xie et al., 2015), resulting in significant differences in the water-bearing characteristics of shale organic matter pores and inorganic mineral pores. It is generally believed that the surface of inorganic minerals is usually strongly hydrophilic, especially clay minerals, which have charges on their surfaces, and water molecules and clay particles can be tightly bound by hydrogen bonds, electrostatic forces, and intermolecular forces (Wang and Li, 2006). However, whether there is water in the pores of organic matter is still controversial, and the water may be controlled by factors such as kerogen type, maturity, and functional groups. The organic pores formed during hydrocarbon generation are generally considered to be oil-wet (
Odusina et al., 2011), and the pores hardly contain water. However, some studies have shown that the functional group structure in the pores of kerogen makes it exhibit a certain hydrophilic ability, and the water molecules in the pores can be adsorbed on the surface of the functional group (
Larsen and Aida, 2004), and distributed in the shale organic matter within pores (
Cheng et al., 2017;
Li et al., 2019a).
Therefore, the current key geological problem is to clarify in which pore spaces water is stored in shale with different water saturation and reveal the control mechanism of shale material composition and pore-throat system on the occurrence of pore water. The study found that the contents of TOC and pyrite in the well LU206 shales are positively correlated with the amount of adsorbed water; the contents of clay and pyrite are positively correlated with the amount of free water (Fig.12). This shows that the organic matter (solid bitumen) pores of shale contribute to the amount of adsorbed water and free water, and the pores of clay minerals contribute to the amount of free water.
Finally, the microscopic occurrence mechanism of shale matrix pore water is summarized. The specific surface area of shale controls the amount of adsorbed water, and the pore volume controls the amount of free water; organic pores in shale solid asphalt contribute specific surface area and pore volume, and inorganic pores in clay mineral contribute pore volume. Therefore, the pores of shale solid asphalt contribute to the amount of adsorbed water and free water, and the pores of clay minerals mainly contribute to free water.
4 Conclusions
1) This study has shown a quantitative evaluation of pore water in deep shales. As a case study, shale samples were obtained from the Wufeng-Longmax Formations in well LU206, located in the Luzhou area, southern Sichuan Basin, China. The amounts of adsorbed and free water in saturated shales were obtained based on the saturation-centrifugation experiments combined with a new pore water classification scheme proposed by
Li et al. (2019a). Under the experimental conditions (30°C, distilled water), the amount of free water varies from 9.5511 to 19.802 mg/g (mean 13.9541 mg/g); the amount of adsorbed water is calculated by subtracting the amount of free water from the total water amount and varies from 1.7967 to 9.8218 mg/g (mean 6.4501 mg/g). The ratio of adsorbed water to total water is 15.83%–42.61% (mean 30.45%).
2) A coupling relationship between shale compositions and pore microstructure was established. The positive correlations of TOC and pyrite contents with specific surface area reflect that shale organic pores (mainly developed in solid bitumen) contribute to shale-specific surface area. Clay and pyrite contents are positively correlated with pore volume reflecting that the shale clay mineral intragranular pores and solid bitumen intragranular pores comprehensively contribute to the shale pore volume.
3) The impact of shale material composition and microstructure on the pore water was analyzed. Based on this, the microscopic occurrence mechanism of shale matrix pore water is revealed. The specific surface area of shale controls the amount of adsorbed water, and the pore volume controls the amount of free water. The pores of shale solid asphalt contribute to the amount of adsorbed water and free water, and the pores of clay minerals mainly contribute to free water.