Abstract
Coalbed methane (CBM) resources cannot be efficiently explored and exploited without a robust understanding of the permeability of fracture-size heterogeneities in coal. In this study, two sister coal samples were imparted with pre-developed cleat and connected fractures, and the permeability of the coal samples was measured under different conditions of controlled confining and gas pressures. Furthermore, the implications of the results for CBM exploration and exploitation were discussed. The permeability of coal with cleat development ranged from 0.001–0.01 mD, indicating ultra-low permeability coal. The gas migration in this coal changed from a linear flow to a non-linear flow, with the increase in gas pressure (>1 MPa). Thus, the permeability of the coal initially increased and then decreased. However, the Klinkenberg effect does not exist in this ultralow-permeability coal. For the coal sample with connected fracture, permeability ranged from 0.1–10 mD, which is larger by hundred orders of magnitude than that of the sample with cleat. For this coal, with a decrease in gas pressure (<1 MPa), the Klinkenberg effect significantly increased the permeability of the coal. With an increase in the applied confining pressure, both the Klinkenberg coefficient and permeability of the coal presented a decreasing trend. It is suggested that field fracture investigation is a prerequisite and indispensable step for successful CBM production. The coal beds that cleat network is well conductive to the connected fracture can be an improved target area for CBM production. During CBM production, a variety of flow regimes are available owing to the decrease in CBM reservoir pressure. In particular, under the low CBM reservoir pressure and low in situ geo-stress conditions, the gas migration in the CBM reservoir with connected facture development exhibits remarkable free-molecular flow. Thus, the reservoir permeability and predicted CBM production will be enhanced.
Introduction
Coalbed methane (CBM) is a rapidly growing source of clean energy worldwide. In addition to the USA, Canada, and Russia, the Asia-Pacific countries are presently the fastest-growing CBM markets, with a compound average growth of 14.9% between 2014 and 2020 (Bandyopadhayay et al., 2020). The efficient exploration and exploitation of this huge resource necessitate accurate reserve estimation and production forecasting. This cannot be achieved without a robust understanding of the permeability of coal (Kumar et al., 2018). This is also a prerequisite for N2 or CO2 enhanced CBM recovery (Anggara et al., 2016; Flores et al., 2019; Zhang and Ranjith, 2019).
Currently, significant research has been conducted on the permeability of coal, worldwide (Connell et al., 2016; Gensterblum et al., 2014; Li et al., 2019; Niu et al., 2020; Ramandi et al., 2015; Wang et al., 2017; Wang et al., 2018; Xie and Zhao, 2017). Gensterblum et al. (2014) studied the relationship between coal deformation and permeability under the action of cyclic loading and inferred that the volumetric strain of the coal gradually increases with the increase in temperature under applied loading, with a corresponding decrease in its permeability. Li et al. (2019) revealed that the permeability stress sensitivity of coal in the loading process was significantly higher than that in the unloading process. Ramandi et al. (2015) established that a direct relationship existed between the coal lithotype and permeability, that is, bright coal is more permeable than dull coal. Connell et al. (2016) inferred that with an increase in gas pressure, the permeability of coal decreases significantly owing to the adsorption-induced swelling strain. Wang et al. (2017) discussed the relationship between coal permeability and porosity under high-temperature and high-pressure conditions and concluded that an exponential relationship exists between the permeability and porosity of coal. With the increase in temperature, the porosity of low-rank coal increased. However, the porosity decreased for medium-rank coal. Xie and Zhao (2017) revealed that the change in pore size distribution in coal determines the permeability evolution behavior with temperature. Wang et al. (2018) investigated the relationship between coal permeability and fracture fractal dimension and concluded that the permeability of coal increases exponentially according to the fracture fractal dimension. Niu et al. (2020) discussed the permeability of coal from the perspective of vertical bedding and parallel bedding and found that the permeability of coal is anisotropic, with a higher permeability in the parallel bedding direction. Generally, research on the permeability of coal primarily has focused on the influences of stress, temperature, pore pressure, coal lithotype, bedding, and adsorption effect.
In addition, recent studies have exhibited that the permeability of coal is intimately related to the Klinkenberg effect, which is essential for enhancing the CBM productivity at lower reservoir pressures. Li et al. (2014) found that the Klinkenberg effect positively contributes towards improving the permeability of coal for all ranks, which becomes more influential for low-permeability coal. Gensterblum et al. (2014) indicated that the Klinkenberg effect in coal could, at least partly, offset the permeability decrease resulting from the increased effective stress during the CBM production process. This leads to a significant increase in the CBM productivity during the low reservoir pressure phase. Zhou et al. (2016) demonstrated that the Klinkenberg coefficient increases with increasing pore pressure at constant effective stress. For a given confining stress on the coal, the Klinkenberg coefficient reaches a maximum value for the competition relationship between the sorption-induced strain and the stress-strain of the coal.
However, coal is a fracture-size heterogeneous material (Du et al., 2020; Firouzi et al., 2014; Wang et al., 2018). In particular, for high-rank coal, cleats are massively developed as natural fractures, which are usually present in two sets that are perpendicular to the bedding and also mutually perpendicular. In terms of their size, the cleats are meso-fractures. In several situations, cleats are confined to the coal of particular maceral types, such as vitrain and clarain. In a coal bed, the connectivity of the cleat network is commonly limited owing to its termination at the interfaces of the coal types. However, the connection of the cleat network can be accomplished by the fracture of a larger size than the cleat network. This type of fracture crosses through the interbedded coal types, which are macro-fractures in size and are termed as connected fractures. These multiscale fractures act as channels for gas permeation, migration, and flow during CBM production. It not only determines the permeability of the CBM reservoir but also controls the productivity of the CBM (Connell et al., 2016; Laubach et al., 1998; Tan et al., 2018). Currently, limited fundamental knowledge exists on the permeability of cleats and connected fractures for fracture-size heterogeneities in CBM reservoirs. Therefore, this study aimed to provide the steps required for a better understanding of the phenomena involved in CBM production. It investigates the permeability of CBM reservoirs using two typical sister samples of coal with cleat and connected fracture development, implemented from the Permian coal measures, Southern Qinshui Basin, China, along with discussions regarding its implications for CBM exploration and exploitation.
Materials and methods
Sample collection and preparation
The Qinshui Basin of China is presently a fast-growing target block for CBM exploration and exploitation of high-rank coal in the world (Du et al., 2019; Ye et al., 2007; Zhou et al., 2016). After experiencing four-stage tectonic movements, the Qinshui Basin deformed into a synclinorium structure (Figure 1). For this coal-bearing basin, in the Permain and Pennsylvania coal-bearing measures, the No. 3 coal in the Shanxi Formation and No. 15 coal in the Taiyuan Formation are the primary minable beds. The coal samples for the experiment were collected from the Shanxi Formation of the Permian strata near the edge of the Qinshui Basin. The sampling method was utilized, as in previous studies (Cai et al., 2014). Coal blocks with a volume of 30 × 30 × 30 cm3 were collected from the mining faces of underground coal mines and were transported to the laboratory for further investigation.

Geological map of the qinshui basin, including coal sampling location and coal bed profile (modified from Cai et al., 2015).
Figure 2 displays an illustration of the coal samples collected. This exhibits that cleats in the coal occur in two sets: face cleats and butt cleats, which are mutually perpendicular and perpendicular to the bedding (Figure 2(a)). For this high-rank coal, cleats are primarily confined to particular coal types, such as vitrain and clarain. An investigation from the field displays that the cleat network among the coal types can be connected by the fracture of a larger size. This type of fracture crosses through the interbedded coal types, and even the coal bed, which is important for enhancing the permeability of the CBM reservoir. To distinguish with the cleat, this type of fracture is termed as a connected fracture herein (Figure 2(b)).

Collected coal samples: (a) coal blocks, (b) cleat and connected fracture patterns in the coal blocks.
Knowledge regarding the permeability of cleats and connected fractures is essential for permeability estimation of CBM reservoirs and production forecasting of CBM, as they influence the control of gas migration. To prepare these sister samples with cleat and connected fractures, the coals were cored vertical to the bedding from the block samples, following standards of the International Society for Rock Mechanics (ISRM) (Figure 3).

Illustration of coal core preparation: (a) with cleat, (b) with connected fracture.
Figure 4 displays the sister coals with cleat and connected fractures, which are numbered C-1 and C-2, respectively. Their dimensions are 10 cm in length and 5 cm in diameter.

Two sister coals: (a) C-1 sample with cleat: The top view (left) and the side view (right), (b) C-2 sample with connected fracture: The top view (left) and the side view (right).
Device for experimentation
Permeability measurements were performed using the following device (Figure 5). This device consists of six subsystems, as shown in the illustration. Part 1 provides gases including N2, CO2, and He; Part 2 is equipped with a temperature control apparatus to achieve a given temperature, Part 3 is the holder for placing the coal core, Part 4 is the hydraulic pump to provide confining stress, Part 5 is utilized to measure the gas flow through the core, and Part 6 contributes towards vacuuming the device.

Illustration of the device for permeability measurement.
Method of experimentation
For the sister coals, the permeability experiment was conducted at a constant temperature (30 °C) with controlled confining pressures and gas pressures. For the permeability measurement experiment, two schemes were performed (Table 1): (a) permeability measurement under low gas pressure (
The experiment scheme for the permeability measurement.
Notably, for the C-1 coal, because of its poor permeability, the gas cannot penetrate the coal core as the gas pressure is less than 1 MPa. Therefore, the C-1 sample was used for permeability measurement under high gas pressure (> 1 MPa). For the C-2 coal, owing to its higher permeability, the gas penetrates the coal core effortlessly, and the gas flow exceeds the maximum measurement range of the flowmeter as the gas pressure exceeds 1 MPa. Therefore, the C-2 sample was utilized for permeability measurement under low gas pressure (< 1 MPa).
For the C-1 sample, inert gas (He) was utilized for permeability measurements. For the C-2 sample, in addition to the inert gas (He), the sorbing gases (N2 and CO2) were also utilized for permeability measurements. For each given confining pressure, the gas pressure increased step-by-step. As all the gas pressure steps were completed, the confining pressure was increased to the following step. The measured gas permeability of the coal was acquired implementing the modified Darcy's law as follows:
Results and analysis
Method of analysis
In previous studies, several methods have been proposed and adopted to evaluate the permeability of coal reservoirs (Cui and Bustin, 2005; Levine, 1996; Seidle and Huitt, 1995; Shi et al., 2014; Wu et al., 2010). Generally, these methods can be divided into two types: (1) porosity-permeability models, and (2) stress-permeability models. Herein, the following simplified stress-permeability model was utilized to analyze the permeability characteristics of the sister coals (Seidle and Huitt, 1995):
According to Eq. 2, under the condition of constant gas pressure with variable confining pressure, the relationship between the permeability and confining pressure is as follows:
Similarly, under the condition of constant confining pressure with variable gas pressure, the relationship between the permeability and gas pressure is as follows:
Permeability of coal with cleat
Figure 6 displays the measured permeability results for the C-1 sample. For the high-rank coal, the permeability of the coal with cleat ranges is observed to be in the order of 0.001–0.01 mD, which is an ultra-low permeability but agrees with the field result from the Qinshui Basin (Li et al., 2014). This indicates that, for the CBM reservoir of high-rank coal, if the reservoir consists of coal with only cleat development, the reservoir will be nearly impermeable. With the increase in gas pressure, the permeability of the coal exhibits a trend of initial increase and later decrease, under each of the applied confining pressures (Figure 6(a)). Incorporating the gas pressure corresponding to the peak permeability as the critical pressure, the permeability of the coal exhibits an opposite evolutionary trend before and after the gas pressure. Figure 6(b) demonstrates the relationship of the gas flow to the gas pressure for the C-1 sample, where the gas flow is measured by the high-precision gas flowmeter (0.01mL·min−1). Similar to the permeability, the gas flow also displayed an opposite evolutionary trend before and after the peak flow at a given gas pressure. Before the peak flowed through the coal core, the gas flow exhibited a nearly linear increase in the gas pressure. The larger the applied confining pressure on the coal, the larger the peak flow and the critical pressure. However, after the peak flow, the gas flow through the coal decayed significantly with an increase in gas pressure.

Experimental results for the coal with cleat (C-1 sample): The relationship of permeability to gas pressure (a), the relationship of gas flow to gas pressure (b).
In Eq. 1, parameters such as the dynamic viscosity coefficient of the gas, coal core cross-sectional area, and length are constant in terms of permeability. Compared to gas flow, permeability is only an indirect indicator that reflects the ability of gas to penetrate through the coal under a certain confining pressure and gas pressure. Previous studies have revealed the relationship between the fluid flow and fluid pressure gradient (Zhou and Wang, 2004):
Figure 7 exhibits the relationship between the permeability and the confining pressure for the coal with cleat. Generally, the permeability of coal displays a decreasing trend with an increase in the applied confining pressure.

Relationship between the permeability and the confining pressure for the coal with cleat.
Table 2 displays the fitting relationship between the permeability and confining pressure for the coal with cleat using Eq. 3. In Eq. 3, the parameter
Fitting relationship between the permeability and the confining pressure for the coal with cleat.
Permeability of coal with connected fracture
Figure 8 displays the relationship between the permeability and confining pressure for C-2 coal. It can be observed that, for the high-rank coal, the permeability of the coal with connected fracture is in the order of 0.1–10 mD. Compared to coal with cleats, coal with connected fractures has an ultra-high permeability. Generally, with an increase in the applied confining pressure on the coal, the permeability of the coal exhibits a downward trend. In particular, compared to the initial confining pressure, the permeability of the coal decreased sharply with an increase in the applied confining pressure.

Relationship between permeability and confining pressure for the coal with connected fracture: He(a), N2(b), CO2(c).
Interestingly, in terms of sorbing gas N2, as the gas pressure is 0.1, 0.2, and 0.3 MPa, the relationship between the permeability of the coal and the confining pressure presents a convex curve trend (Figure 8(b)). As the gas pressure is 0.1 MPa, the coal exhibits a maximum permeability of 1.95 mD for an applied confining pressure of 6 MPa. This is higher than the neighboring permeability of 0.56 mD for an applied confining pressure of 4 MPa. Subsequently, the permeability of the coal exhibited a decreasing trend with an increase in the applied confining pressure. This may be due to: (1) the experimental error during the operation; (2) less data (three data points) which do not reflect the true regularity; (3) absorption effect; and (4) the gas molecular slippage effect. In practice, during the process of permeability measurement, regardless of the non-absorbing gas (He) or the sorbing gases (N2 and CO2), all the operations conform to the identical procedure. Thus, the first three reasons were excluded. However, it can be observed that, at a given confining pressure of 6 MPa, with the increase of gas pressure from 0.1 MPa to 0.4 MPa, the permeability of the coal decreases. In this situation, although the effective stress acting on the coal decreases, owing to the gas molecular slippage effect the permeability of the coal is higher at lower gas pressures. Besides, with the increase of gas pressure from 0.1 MPa to 0.4 MPa, the permeability convex curve for N2 gradually converts to the concave curve, consistent with that of He and CO2. Therefore, it is concluded that the gas molecular slippage effect results in a complicated permeability relationship for ultra-high permeability coal with connected fractures under low gas pressure. According to Klinkenberg, the gas molecular slippage effect is closely related to the gas molecular mean free path and fracture aperture (Klinkenberg, 1941). For confined coal, the fracture aperture of the coal is closely related to the applied confining pressure. Thus, the gas molecular slippage effect was determined by the gas pressure and confining pressure. At lower pressure conditions, the gas molecular mean free path increases, and the slippage effect and permeability of the coal are enhanced.
For the C-2 coal, the relationship between the permeability and confining pressure under each gas pressure was fitted by Eq. 3 (Table 3). It can be observed that under different gas pressures, the relationship between the permeability of the coal and the applied confining pressure satisfies an exponential law. The measured permeability by the non-absorbing gas (He) is between 1.7 to 10.8 mD, the measured permeability by the sorbing gas (N2) is 5.1 mD, and by the sorbing gas (CO2) is between 2.1 to 48.6 mD. Among these, the CO2 measured permeability was the largest Several studies have demonstrated that there are obvious absorption-induced swelling effects for CO2 by coal, which will narrow the fracture aperture and decrease the permeability. However, this further implies that, for coal with connected fractures or ultra-high permeability, under low gas pressures (
The fitting relationship of the permeability and the confining pressure for the coal with connected fracture.
Note: There is no fitting relationship for N2 at gas pressures of 0.1, 0.2, and 0.3 MPa.
Figure 9 displays the permeability characteristic parameters of the C-2 coal. With the increase in gas pressure, the fracture compression coefficient of the coal exhibits a decreasing trend. Incorporating instance from the non-adsorbing gas (He), under each gas pressure, the fracture compression coefficient of the coal with connected fracture (0.0690–0.1429MPa−1) is larger than that of the coal with cleat (0.0511–0.2742 MPa−1). This suggests that compared with the coal with cleat, fracture closure via applied confining pressure is more extensive for coal with connected fractures. For this reason, the permeability stress sensitivity of the coal with cleat is inapparent. Furthermore, with an increase in gas pressure, the permeability stress sensitivity coefficient of the coal decreases. Compared with non-absorbing gas (He) and sorbing gas (N2), an approximately negative linear relationship exists between the gas pressure and the permeability stress sensitivity coefficient for the sorbing gas (CO2). The permeability stress sensitivity coefficient of the sorbing gas (CO2) is larger than that of the sorbing gas (N2). This can be attributed to the more significant sorption-induced fracture closure effect for the sorbing gas (CO2).

The relationship of permeability characteristic parameters to the gas pressure for the C-2 sample under low gas pressures (< 1 MPa).
Figure 10 displays the relationship between permeability and gas pressure for C-2 coal under each confining pressure. In addition to the condition of applied confining pressure at 4 MPa in Figures 10(a) and 10(b), the measured permeability of the coal by the non-absorbing gas (He) and the sorbing gases (N2 and CO2) presents a decreasing trend with increasing gas pressure.

The relationship of the permeability to the gas pressure for the coal with connected fracture: He(a), N2(b), CO2(c).
In fact, with the change in gas pressure and applied confining pressure, the gas flow in the fracture-size heterogeneities material can be variable. Under lower pressure conditions, the gas flow regime for the gas molecules can be characterized by studying the Knudsen number (
Knudsen numbers and flow types.
For example, as the mean free path at a given gas pressure is 10 nm and the average fracture aperture is 50 nm, the Knudsen number is 0.2, thereby implying a transition flow regime. Clearly, a variety of flow types can occur in the fracture-size heterogeneities of coal with variable confining pressures and gas pressures. With the decrease in gas pressure, as the mean free path is far more than the fracture aperture, the gas will exhibit a free-molecular flow during the flow regimes. When the gas velocity at the walls is non-zero, the gas transport is enhanced as the “gas slippage” reduces the viscous drag and increases the apparent permeability. For this situation, the gas flow exhibits the Klinkenberg effect, and the permeability of the coal can be given by (Klinkenberg, 1941):
From Table 5, it can be observed that, for the non-absorbing gas (He) or the sorbing gases (N2 and CO2), the relationship of the measured permeability by the gas to the gas pressure conforms to the expression of Eq. 8. Compared to the adsorbing gases (N2 and CO2), the non-adsorbing gas (He) has a larger Klinkenberg coefficient in general (Figure 11). For a given confining pressure, the value of the Klinkenberg coefficient is determined by the combined effect of the stress–strain and sorption-induced strain (Zhang et al., 2018). This indicates that the gas slippage phenomenon for the non-absorbing gas is more prominent than that for the adsorbing gas. In addition, for all the gases, with an increase in the applied confining pressure, the Klinkenberg coefficient of the coal exhibits a decreasing trend. The increased confining pressure not only weakened the Klinkenberg effect but also decreased the Klinkenberg permeability of the coal.

Permeability characteritic parameter for the coal with connected fracture under the variable confining pressure.
The fitting relationship of the permeability and the gas pressure for the coal with connected fracture by Eq. 8.
For the C-2 coal, the fitting relationship of the fracture compression coefficient to the gas pressure (
The fitting relationship of permeability and the gas pressure for the coal with connected fracture.
Implication for CBM exploration and exploitation
Coal is a fracture-size heterogeneous material, which includes macro-, meso-, and micro-fractures. Cleats act as meso-fractures, usually referred to as natural fractures, which are formed and developed by the coupling mechanism of in-situ geo-stress, coal lithotype, and coal bed thickness during the process of coalification, and are closely related to the coal rank. However, the connected fracture can be of different sizes, including meso -and macro-fractures. This type of fracture usually forms and develops owing to geologic movements, especially in the coal basin suffering from the uplift movement or in the tectonic regions due to the brittle geologic deformation action.
As mentioned previously, in high-rank coal, the cleat (0.001–0.01 mD) is far less permeable than connected fracture (0.1–10 mD). In addition, field investigations exhibit that cleats are usually confined to particular lithotypes of coal, such as vitrain or clarian, whose connectivity is limited by its termination at the interfaces of coal types. Thus, to determine the target blocks for CBM exploration and exploitation, the locations where the cleat network of the coal bed is well connected to the connected fracture may be the best target. Thus, field fracture investigation is a prerequisite and indispensable step for successful CBM production.
In addition, during the process of CBM production, with the downdraw of reservoir pressure, owing to the gas pressure change in the CBM reservoir, a variety of flow regimes are present. In particular, under the conditions of lower reservoir pressure and lower in situ geo-stress, the gas migration in the CBM reservoir exhibits remarkable free-molecular flow. The reservoir permeability and predicted CBM production were enhanced.
Conclusions
The permeability of high-rank coal with only cleat development is 0.001–0.01 mD, which represents ultra-low permeability coal. No Klinkenberg effect was observed in this coal. With an increase in the gas pressure (> 1 MPa), the gas migration in the coal changes from linear flow to nonlinear flow. This caused the permeability of the coal to initially increase and then decrease with the increase in gas pressure.
For high-rank coal with connected fracture development, its permeability is in 0.1 of 10 mD, which is much larger than that of the coal with only cleat development. Under low gas pressure conditions (<1 MPa), with the decrease in gas pressure, the Klinkenberg effect significantly increased the permeability of the coal. The Klinkenberg effect for the non-absorbing gas was more prominent than the adsorbing gas. In addition, with an increase in the applied confining pressure, the Klinkenberg coefficient and permeability of the coal presented a decreasing trend.
The stress sensitivity of permeability for coal with cleats was weaker than that of coal with connected fractures. However, because of the sorption-induced fracture closure effect, the permeability stress sensitivity for the sorbing gas (CO2) was larger than that of non-adsorbing (He) and sorbing gas (N2). With the increase in gas pressure, the stress sensitivity of the permeability for both coal dropped significantly.
Field fracture investigation is a prerequisite for successful CBM production. It may be the best block for CBM production where the cleat network is well conductive to the connected fracture. During the CBM production process, with the downdraw of reservoir pressure, a variety of flow regimes are present. In particular, under the lower reservoir pressure and lower in situ geo-stress conditions, the gas migration in the CBM reservoir exhibits remarkable free-molecular flow, and the reservoir permeability and predicted CBM production were enhanced.
Footnotes
Funding
This study has been sponsored financially by the Science and Technology Development Project of Luoyang (Project No. 2101025A), and the Projects of National Natural Science Foundation of China (Grant Nos. 51804100, 51774113, U1904126); Besides, this study has also been funded by the Key Science and Technology Research Projects of Henan Province (Grant No. 202102310612), and Science and Technology Innovation Leading Talents of Zhongyuan (Grant No. 214200510030).
Author contributions
Zhigang Du and Xiaodong Zhang have contributed to the experiments and the writing of this paper. Qiang Huang and Yawen Tao have contributed to the mapping of this paper. Wuxiu Ding has provided language help. Each author has contributed to this study.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (grant number 51774113, 51804100, U1904126, 202102310612, 214200510030, 2101025A).
