#### Abstract

Hydraulic fracturing and premining gas drainage are important to safe mining and coalbed methane extraction. These technical processes cause the redistribution of in-situ stress and the regional variation of moisture contents within the affected zone. Therefore, we investigated the coupled effect of variable stresses (from 9 MPa to 27 MPa) and moisture contents (from 0.22% to 4.00%) on the permeability evolution of gas-saturated raw coal. The results show that (1) the relationship between the mean effective stress and the permeability can be described by a power function according to the permeability evolution model of the porous matrix established in this study. Besides, the influence mechanisms of moisture on fitting coefficients in the function were analyzed. (2) The permeability decreases with the increase of in-situ stress (e.g., confining pressure or volumetric stress) in a negative exponential manner. (3) The curves of permeability variations with moisture content are not always linear, and the permeability is more sensitive to the moisture content than the volumetric stress in the test range. Moreover, the sensitivity of permeability varies in different regions. These results would be beneficial for permeability prediction and surface well parameters design.

#### 1. Introduction

Coalbed methane (CBM) is an abundant valuable resource in underground coal mines. It is estimated that China bears approximately 36.8 trillion m^{3} of CBM in its reservoirs shallower than 2 km, ranking the third country in the world [1]. Qinshui Basin whose CBM reserve accounts for 1.08% of the total CBM reserve in China is one of the largest CBM reservoirs [2]. However, the reservoir is generally characterized by low permeability, which seriously inhibits efficient CBM extraction [3]. In-situ surface well extraction is a technique for enhancing the permeability and relieving the pressure of coal seams. Permeability is a key parameter for surface well design, and its temporal and spatial variations remarkably affect the occurrence state, migration, and extraction of CBM [4]. Scholars have conducted extensive researches on the factors influencing permeability, including in-situ stress field [5], gas pressure [6], Klinkenberg effect [7], geothermal field [8], geoelectric field [9], acoustic field [10], and physical properties [11]. These findings have allowed the mechanics of permeability evolution to be deliberated and also provided an adequate data volume and dimension for machine learning which has been proven to be a powerful tool for permeability prediction [12, 13]. Specifically, the process of fracturing, drainage, and gas extraction in surface wells destroys the original stress distribution balance and results in the existence of stress concentration or relief region [14]. Besides, the large amount of fracturing fluids leads to regional differences in water-bearing conditions. Therefore, special attention needs to be paid to in-situ stress and moisture content of coal seam for the sake of permeability evolution investigation and prediction, gas flow simulation around the wells, and parameter design of CBM surface wells.

The development of uniaxial and triaxial loading devices has actively promoted the knowledge of permeability evolution characteristics under different stress conditions [15, 16]. The permeability-strain curve of coal corresponds well to its full stress-strain curve. To be specific, the permeability drops first and then rises as a coal sample gets loaded, deformed, and destructed [5]. The peak value of coal permeability lags behind that of stress and strain, suggesting that the characteristics of CBM flow are closely related to the evolution of coal damage generated in the loading process. In addition, scholars have carried out abundant researches on seepage tests under complex loading/unloading paths and loading/unloading rates and on the establishment of permeability models in the case of multifield coupling [17], but these researches are mostly focused on dry coal. The permeability evolution of water-bearing coal in the case of moisture-stress coupling is rarely reported.

The influence of moisture on the permeability of CBM reservoirs is primarily reflected in coal deformation, gas desorption, and migration [18]. Pan et al. [19] believed that moisture in the coal matrix would cause coal swelling/shrinkage and mechanical property alteration that would impact on coal permeability under reservoir conditions. Zhao et al. [20] and Gai et al. [21] investigated the desorption law of coal subjected to high-pressure water injection and made a comparison with the natural desorption state. Guo and Su [22] conducted laboratory tests on the starting pressure gradient and permeability under different water saturation conditions. The test results showed that the permeability gradually declined and the starting pressure gradient gradually rose with the increase of water saturation degree. Nie et al. [23] explored the microscopic mechanism of gas adsorption on coal samples with different moisture contents in light of molecular thermodynamics and surface physicochemical theories. Wang et al. [24] probed into the relationship among moisture content, porosity, and permeability of fractured coal. Yin et al. [25], Liu et al. [26], Wei et al. [27], and Yuan and Jiang [28] analysed the seepage characteristics of gas-bearing coal with different moisture contents and found that the moisture content and permeability shared a linear negative relationship or negative exponential variation law. Hao et al. [29] performed axial and radial gas seepage experiments and revealed that the axial and radial permeabilities of coal first increased and then decreased with the rise of moisture content. Nevertheless, these researches failed to achieve a unified understanding, and the varied research results imply the difficulty in clearly explaining the effect of coupling terms on gas adsorption, diffusion, and seepage.

This study is focused on investigating the permeability evolution of gas-saturated raw coal samples with different moisture contents under varying stresses. Firstly, the experimental conditions were simplified by setting constant temperature and gas pressure, unidirectional loading path, and resaturated adsorption after the break of balance, so as to avoid irrelevant factors. Then, the relationship between effective stress and permeability was described by a new function that has a clear physical meaning based on the established permeability evolution model of the porous matrix. Furthermore, the mutual effect of moisture content and volumetric stress on permeability were obtained, and conclusions different from the abovementioned ones were drawn. Finally, a sensitivity analysis was performed in the test region.

#### 2. Materials and Methods

##### 2.1. Engineering Background

Yuwu Coal Mine is located in the south of Qinshui Basin, China. The buried depth of S2107 working face of the mine lies in the range of 480-543 m. No. 3 coal seam that is being mined belongs to a high-gas and low-permeability coal seam with a gas content of 7.71, a gas emission initial volume from a 100-meter-long borehole of 1.42, and a permeability coefficient of 0.28-0.42. The proximate analysis results and adsorption constants of coal samples from the No. 3 coal seam are listed in Table 1.

Surface wells #1-52 and #1-54 were constructed for conducting permeability-enhancing fracturing on S2107 working face before gas extraction. No. 1-No. 12 gas extraction drilling sites were arranged along the tailentry, and the single-hole gas flow rate was recorded within 50 d. During this period, the drilling sites were not affected by mining activities. The analysis results show that the influence radii of surface wells #1-52 and #1-54 are 85-100 m and 39-42 m, respectively. Coal samples YW1-YW4 were taken from the affected zone in Figure 1.

##### 2.2. Coal Sample Preparation

The steps of coal sample preparation are shown in Figure 2. (1) Standard raw coal samples with the size of *ϕ*50(100 ± 2) mm were prepared with the aid of a wire cutting machine (Figure 2(a)), and relative homogeneity of the samples was ensured through density measurement and wave velocity measurement (PDS-SW ultrasonic detector). Samples YW1-YW3 almost bear no cracks (wave velocity range 1921.8-2050.8), while Sample YW4 owns obvious cracks on its surface (wave velocity range 1826.8-2033.8 ). (2) Preparation of coal samples with different moisture contents : YW1 (), the samples that had experienced hydraulic measures, did not receive any treatment; YW2 (), the samples with an ultra-low moisture content (lower than the original moisture content 0.49%), were prepared by controlling the drying time of drying oven (Figure 2(b)); YW3 (), the samples with a high moisture content, were prepared by controlling the pressure holding time of pressure container (Figure 2(c)); YW4 () were the crack-containing sample with a high moisture content.

**(a)**

**(b)**

**(c)**

##### 2.3. Testing Apparatus

The testing apparatus is mainly composed of a loading frame, a servo hydraulic station, an air path system, a triaxial chamber, a constant-temperature oil bath, and a data acquisition system. The schematic diagram of the apparatus is displayed in Figure 3. The main technical parameters are as follows: axial stress range 10-800 KN, confining pressure range 0-15 MPa, gas pressure range 0-15 MPa, and temperature range from room temperature to 260. The chamber equipped with a temperature sensor (PT100, ±0.01) and a circumferential extensometer (Epsilon 3544, made in USA) was placed in the constant-temperature oil bath for maintaining a constant temperature.

##### 2.4. Experimental Contents and Procedures

The experimental conditions were simplified by setting constant temperature and gas pressure. The temperature was maintained at 20. In the process of gas pressure determination, the Klinkenberg effect needs to be avoided, because a small gas pressure difference can cause gas slippage. More specifically, the Klinkenberg effect refers to the slip phenomenon that the gas flow velocity at the wall does not equal zero when the average molecular free path of gas approximates the pore size of a porous medium. It is manifested in varying ways under different effective stresses [5, 6, 30, 31], moisture contents [27, 28, 32], and temperatures [33, 34]. In view of this fact, the permeabilities of coal samples from Qinshui coalfield were tested at 20 under variable gas pressures. As presented in Figure 4, the Klinkenberg effect disappears when the gas pressure exceeds 1.0-1.5 MPa, and its influence on the permeability becomes weaker with the rise of moisture content or volumetric stress.

**(a)**

**(b)**

Therefore, considering the in-situ stress value of No. 3 coal seam and the test results of Klinkenberg effect inflection point, the volumetric stress range and gas pressure were set as 9-27 MPa and 1.5 MPa, respectively, in this experiment. The prepared samples with different moisture contents were loaded in accordance with the path in Figure 5, which can ensure that the axial pressure is higher than the confining pressure and the applied unidirectional external stress is lower than the peak strength. Gas pressure at the outlet was 0.1 MPa. Saturated adsorption treatment was conducted on the samples when the stress condition was altered, in order to prevent the adsorption/desorption process from entangling the result analysis.

Experimental procedures: (1) Sample installation: First, the side of the sample was coated with 704 silicone rubber. Then, the sample was wrapped in a heat-shrinkable tube and placed in the chamber. Afterwards, the chamber was put in the 20°C oil bath. (2) Initial adsorption and seepage test: The in-situ stress was raised to the initial value, and 1.5 MPa of CH_{4} was injected. Next, the outlet valve was closed to allow the sample to adsorb CH_{4} until the equilibrium state was reached, i.e., until gas pressure in the chamber ceased changing. Finally, the outlet value was opened to record the steady flowrate. (3) The sample was allowed to readsorb CH_{4} for 2 h under an altered stress, after which the steady flowrate was recorded. Step (3) was repeated along the loading path. (4) Steps (1)-(3) were repeated for another sample.

#### 3. Permeability of Coal Samples: Theoretical Backgrounds

##### 3.1. Calculation Method for Permeability

According to the theory of geotechnical mechanics, the mean effective stress (MES) can be calculated by Eq. (1) [35, 36]: where is the MES, ; is the axial stress,; is the confining pressure, ; is the gas pressure at the inlet, ; is the gas pressure at the outlet, , andit takes the value 0.1 .

Regardless of the starting pressure gradient, the permeability of the coal sample is calculated by Eq. (2) according to Darcy’s law [37, 38]:
where is the permeability of coal sample, 10^{-3} ; is the gas flow rate at the outlet, ; is the atmospheric pressure, 0.1 MPa; is the dynamic viscosity coefficient of gas, ; is the length of deformed coal sample, ; is the area of deformed coal sample, ; is the gas pressure at the inlet, .
where are the length and area of deformed coal sample; and are the original height and diameter of coal sample, ; is the axial deformation of coal sample after the force loading, ; is the radial deformation of coal sample after the force loading, .

##### 3.2. Permeability Evolution Model of Porous Matrix

The Warren-Root model of coal is presented in Figure 6 [39]. The porous matrix, which comprises pore clusters supported by skeletons, is the main gas storage space of coal (Figure 6(a)), while cracks, which completely separate the matrix, are the primary gas migration channel [40]. Figure 6(b) shows the seepage model of pores, in which the matrix is equivalent to parallel capillary tubes with equal diameter and skeletons.

Space of cracks or pores will be significantly reduced if deformation or stagnant water exists. Effective porosity is not only a key indicator for measuring the size of the space, but also an important factor that determines the adsorption/desorption and permeability of coal. The effective porosity of coal is related to its structural deformation and bulk deformation. Structural deformation refers to compression deformation of the skeleton caused by external stress, or the relative dislocation between tubes which make the tubes more closely arranged. It is usually unrecoverable. Bulk deformation refers to tube expansion caused by thermal stress or adsorption swelling, or compression deformation of tubes under the action of gas pressure [41]. Such elastic deformation caused by internal stress can be recovered.

In this study, the permeability evolution model of porous matrix considering moisture is established because Samples YW1-YW3 contain few cracks. Since the matrix will undergo structural deformation (external stress compression) and bulk deformation (adsorption swelling and gas pressure extrusion) under the condition of constant temperature and gas pressure, the dynamic evolution model of porosity of the dry matrix is [14, 41–44]: where is the porosity of dry matrix, %; is the initial porosity of matrix, %; is the volumetric strain, of which the value is negative under compression; is the adsorption expansion strain per unit volume of coal; is the volumetric compression coefficient, .

The unit in the matrix is taken and regarded to be equivalent to the capillary tubes (Figure 6(b)). Given the fact that water distributed in macropores () and mesopores () occupies the space, the porosity of the water-containing matrix is [14]:
where is the density of matrix, ; is the moisture content of matrix pores; *n* is the number of capillary tubes per unit area , tubes/m^{2}; is the density of water in capillary tubes, ; and is the radius of a tube, .

By integrating Eq. (5) with Eq. (6), the effective porosity evolution model of water-containing matrix can be obtained:

The relationship between permeability and porosity is given in the Kozeny-Carman equation established on the basis of the capillary tubes model [44]. Ignoring the change in the total surface area of coal particles per unit volume, a permeability evolution model that takes into account the moisture content in pores can be obtained:
where is the initial permeability of coal sample, 10^{-3}.

#### 4. Results and Discussion

##### 4.1. Effect of Mean Effective Stress on Permeability

Changes in flow rates and permeabilities of coal samples with moisture contents of 0.22%, 1.98%, and 3.24% under mean effective stress (MES) values of 4.2-8.3 MPa are disclosed in Figures 7(a)–7(c). The flow rates of coal samples decrease nonlinearly at a reduced rate with the increase of MES, and the relationship between flow rate and MES can be expressed by a quadratic function where , , and are fitting coefficients.

**(a)**

**(b)**

**(c)**

**(d)**

The permeabilities and flow rates of samples with different moisture contents vary in similar trends (Figure 7(d)). The variation trend of permeability with MES can fitted by the power function which agrees with the form of Eq. (8). In the function, , , and are the fitting coefficients with clear physical meanings (see Table 2 for their values), and the influence mechanisms of moisture on the values of , , and differ.

The fitting coefficient is related to the initial permeability of the matrix with specified moisture content. The value of a drops sharply with the rise of moisture content, and the effect of water on the value of is mainly reflected in the following two nondeformation aspects of matrix: (1) Free water, adhered water, and film water in macropores and mesopores occupy the gas flow channel, thus lowering the effective porosity of coal. In addition, a water film formed on the pore surface generates a certain vapor pressure, thus raising the viscous resistance of gas migration [45, 46]. (2) Within the main space for CH_{4} adsorption and desorption, namely, micropores (0.01-0.1 um) and molecules structural pores (<0.01 um), capillary resistance formed by water hinders gas desorption and exhibits a water-locking effect.

The fitting coefficient is related to the bulk deformation, structural deformation, and initial porosity of the coal matrix. By conducting wave velocity measurement, the initial permeabilities of Samples YW1-YW3 are considered to be consistent, so the variation of value is mainly affected by structural deformation (external stress) and bulk deformation (adsorption swelling and gas pressure extrusion) of the matrix in this test. As the moisture content increases, the volumetric compressive strain goes up (Figure 8); the volumetric compression coefficient grows [47]; and the adsorption expansion strain does down [26, 42, 43, 48], because water molecules possess stronger adsorption competitiveness on the coal matrix surface compared with CH_{4} molecules [49, 50]. Thereby, with the rise of moisture content, the effective adsorption sites of gas on the pore surface decrease, so does the adsorption expansion strain . As revealed by the analysis on the values of , , , and with the change of moisture content, with the rise of moisture content, the effective porosity of matrix increases in terms of bulk deformation, but meanwhile the existence of external stress makes the structural deformation of the matrix with a high moisture content more obvious. As a result, the effective porosity is reduced because of the larger structural deformation.

The fitting coefficient is related to the structural deformation of the coal matrix caused by external stress. Figure 8 displays the curves of volumetric strain variation measured in the test. It can be seen that as the moisture content rises, the absolute value of volumetric strain increases, so does the absolute value of .

##### 4.2. Effect of In-Situ Stress on Permeability

Figure 9(a) demonstrates the permeability variations of samples with moisture contents of 0.22%, 1.98%, and 3.24% under a constant axial pressure of 5.5 MPa and confining pressures of 3.9, 4.1, 4.3, 4.5, and 4.7 MPa. Figure 9(b) shows the permeability variations of samples with moisture contents of 0.22%, 1.98%, and 3.24% under volumetric stresses of 13.3-27.3 MPa.

**(a)**

**(b)**

Obviously, permeability goes downward nonlinearly at a reduced rate with the increase of in-situ stress (e.g., confining pressure or volumetric stress ). The above downward trend can be fitted with a negative exponential function where and are the fitting coefficients and is related to the initial permeabilities of coal samples with different moisture contents. This finding corresponds to previous research results [27, 35, 51]. The reason can be explained as follows. The initial permeabilities of coal samples with different moisture contents differ. Specifically, the higher the moisture content is, the lower the initial permeability is. First, the stagnant water and the seal off the effect of water reduce the effective porosity of coal [29]. Second, with the increase of external stress, the structural deformation of coal brings about further shrinkage of gas seepage channel.

Sample YW4 () contains evident microcracks (Figure 6). The test result shows that its permeability is up to 60% higher than the value in Figure 9 under the same in-situ stress condition. The result fully proves that although the moisture content seriously weakens coal permeability, cracks, as the dominant factor affecting coal permeability [52], will effectively weaken the negative impact of water on coal permeability.

##### 4.3. Mutual Effect of Moisture Content and Volumetric Stress on Permeability

###### 4.3.1. Relationship between Permeability and Moisture Content under Applied In-Situ Stress

Figure 10 exhibits the variation of permeability with the moisture content and the volumetric stress . The curves are not alway linear. In the whole region (=0.22-3.24%), under low volumetric stresses (<15 MPa), decreases nonlinearly with the increase of ; as the volumetric stress grows (e.g., =18 MPa), it becomes linearly correlated with, which agrees with the conclusions of Liu et al. [26], Yin et al. [25, 53], and Yuan and Jiang [28]; under high volumetric stresses (>18 MPa), it is no longer linearly correlated with, the slope of curve altering evidently in the vicinity of =1.98%. Moisture occupies pores and external stress causes structural deformation of the pore, both reducing the effective porosity of matrix and thus resulting in a decrease in permeability.

In different regions ( and ), the slope of curve decreases with the increase of , because influence of on weakens with the increase of . It is noteworthy that in Zone I-1 where is high and is low, the slope of curve remains basically constant with the increase of ; in Zone II-1 where is low and is high, the slope of curve remains basically unchanged with the increase of . The variation trend reflects the differences in permeability sensitivity to volumetric stress and moisture content in different regions.

###### 4.3.2. Sensitivity Analysis of the Test Regions

To determine the influence degrees of moisture content and volumetric stress on permeability , the sensitivity coefficients of permeability to moisture content and volumetric stress were defined as and , respectively. The meanings of and were defined in the same way as above.

First, the correlation between variables was analysed. The Pearson correlation coefficient (PPMCC) shown in Eq. (9), which is the quotient of the product of covariance and standard deviation of two variables [54], is widely used for measuring the correlation between two variables. Through calculation, the correlation between the moisture content and the permeability is -0.782 (highly significant, ), while that between the volumetric stress and the permeability is -0.470 (significant, ). In addition, and are two independent variables. where is the covariance; and is the mathematical expectation.

*(1) Sensitivity analysis in the whole test region.* Based on the correlation analysis, a multiple regression model was established to calculate the sensitivity coefficients and . The data in Figure 10 were subjected to -Score standardization via Stata software and then substituted into the regression equation where , , and are the regression coefficients and is the error term.

The permeability is more sensitive to the moisture content in the whole test region. The overall sensitivity coefficients and were defined. Their values =-0.782 (highly significant, ) and =-0.470 (highly significant, ) can be found in Table 3. In the whole test region, , that is, permeability is more sensitive to moisture content on the whole.

*(2) Sensitivity analysis in different regions.* The interaction variable [55–57] was added into the regression equation for comparing the values of and in different regions. The regional sensitivity coefficients and were defined.

As given in Table 4, the interaction term coefficient is 1.909 (highly significant, ), suggesting that although moisture content and volumetric stress are independent variables, they exert a mutual effect on permeability. and demonstrate that the sensitivity coefficient decreases gradually with the increase of volumetric stress while decreases gradually with the increase of moisture content (Figure 11). Therefore, the sensitivity coefficient decreases with the increase of volumetric stress so that the slope of curve becomes smaller and smaller in the ranges of and (Figure 10), which indicates the weakening influence of moisture content on permeability.

Figure 12(a) is a three-dimensional diagram of the data in Figure 10. As can be observed in Figure 12(b), in Zone I (), that is, permeability is more sensitive to moisture content in this region; in Zone II (), that is, permeability is more sensitive to volumetric stress in the region. The constant slopes of curves in Zone I-1 and Zone II-1 in Figure 10 indicate the absolute influence of high volumetric stress or high moisture content on permeability.

**(a)**

**(b)**

The value of is larger than that of when the volumetric stress is lower than 15 MPa (Figure 11). At this time, coal permeability is not simply linearly correlated but shares a negative exponential relationship with moisture content (Figure 12(a)). This finding differs from the research results of Yin et al. [25, 53], probably because the briquette samples used by Yin et al. possess a homogeneous internal structure while the complex pore structure of raw coal samples used in this test decides the more complex influence mechanism of moisture on raw coal permeability. Under a low volumetric stress, the pores in coal are slightly affected by stress, and the coal is strongly hydrophilic. As a result, the increased water adsorbs on the pore surface easily and meanwhile occupies the effective gas seepage channel. Since the porosity of coal is an important parameter affecting gas permeability, the permeability of coal falls rapidly with the rise of moisture content under a low volumetric stress. With the increase of volumetric stress, the sensitivity of permeability to moisture content decreases, so that permeability becomes linearly correlated with moisture content. After volumetric stress continues to increase to 24 MPa, the coal sample gets more sensitive to volumetric stress (Zone II in Figure 12(b)), which is reflected by the phenomenon that the increase of moisture content fails to notably lower permeability under a low moisture content.

#### 5. Conclusions

In this study, a series of experiments were performed for grasping the permeability evolution of gas-saturated raw coal under the coupled influences of moisture contents and in-situ stress. The main conclusions reached are as follows. (1)A permeability evolution model of a porous matrix was established. On the basis of this mode, it is found that the power function can be used to describe the relationship between MES and permeability . In this function, the fitting coefficient is related to the initial permeability with specified moisture content; is related to the bulk deformation, structural deformation, and initial porosity of coal matrix; and is related to the structural deformation of coal matrix caused by external stress. The influence mechanisms of moisture on these coefficients were analysed according to the experimental results afterwards(2)Permeability decreases nonlinearly with the increase of in-situ stress (e.g., confining pressure or volumetric stress ), which can be fitted by a negative exponential function , where is related to the initial permeabilities of coal samples with specified moisture content (3)The mutual influences of moisture content and volumetric stress on permeability were analysed. In different regions, the slope of curve decreases with the increase of , because the increase of would weaken the influence of on . Overall, the curves are not always linear. The variation reflects the differences in permeability sensitivity to volumetric stress and moisture content in different regions

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare no conflicts of interest.

#### Acknowledgments

The authors are grateful to Professor Yaojiang Zhao from Taiyuan University of Technology for his support and contribution in this study. This work has been supported by the Fundamental Research Funds for the Central Universities (2017XKZD06).