A Comprehensive Coal Reservoir Classification Method Base on Permeability Dynamic Change and Its Application

Due to the unique adsorption and desorption characteristics of coal, coal reservoir permeability changes dynamically during coalbed methane (CBM) development. Coal reservoirs can be classified using a permeability dynamic characterization in different production stages. In the single-phase water flow stage, four demarcating pressures are defined based on the damage from the effective stress on reservoir permeability. Coal reservoirs are classified into vulnerable, alleviative, and invulnerable reservoirs. In the gas desorption stage, two demarcating pressures are used to quantitatively characterize the recovery properties of permeability based on the recovery effect of the matrix shrinkage on permeability, namely the rebound pressure (the pressure corresponding to the lowest permeability) and recovery pressure (the pressure when permeability returns to initial permeability). Coal reservoirs are further classified into recoverable and unrecoverable reservoirs. The physical properties and influencing factors of these demarcating pressures are analyzed. Twenty-six wells from the Shizhuangnan Block in the southern Qinshui Basin of China were examined as a case study, showing that there is a significant correspondence between coal reservoir types and CBM well gas production. This study is helpful for identifying geological conditions of coal reservoirs as well as the productivity potential of CBM wells.


Introduction
Coal reservoir permeability is a key property controlling fluid migration. Many studies and production experiences have shown that permeability is crucial for coalbed methane (CBM) production [1,2]. Coal reservoir permeability varies significantly with variation in effective stress, matrix shrinkage, and gas slippage during CBM development [3][4][5][6]. Specifically, coal reservoirs are sensitive to pressure, and reservoir permeability decreases exponentially with increasing effective stress during the complete depressurization process [7,8]. When the reservoir pressure drops below the critical desorption pressure, CBM desorbs from the coal matrix, and the coal pore structure deforms, resulting in increasing permeability [9][10][11]. Gas slippage also increases the apparent permeability when reservoir permeability is low, but the effect is small [12,13].
Multiple permeability models have been developed to describe coal permeability behavior. Coal permeability models considering the effect of effective stress were first proposed [14][15][16], followed by models considering the effect of sorption-induced strain on coal permeability evolution [17][18][19]. Some models have proven popular with practitioners, such as the Seidle model, P&M model, S&D model, and Pan & Connell model [12,14,[20][21][22][23]. These models either consider more complicated conditions or use different interpretations and have been widely used to match historical field production data. Two important motivations for using coal permeability behavior models include the development of models that can provide a practical means for explaining permeability behavior when analyzing reservoir behavior or during laboratory tests and improve our fundamental understanding permeability control mechanisms, which are more theoretically focused [24]. Most scholars obtained a parabolic relation between permeability and reservoir pressure. The permeability of coal reservoirs changes significantly during gas production, commonly initially decreasing but then increasing as the reservoir pressure and gas content are drawn down [20,24,25]. However, in some cases, reservoir permeability does not increase after gas desorption, but it decreases continuously with decreasing reservoir pressure, or, conversely, reservoir permeability increases monotonously during the complete production process. Therefore, there are still some crucial questions, for example: what are the specific characteristics of the permeability dynamic behavior; what parameters can be used to quantitatively characterize this dynamic behavior; what is the physical significance of these parameters; how do these parameters affect coalbed methane (CBM) development? Additionally, the roles of coal physical and chemical properties on permeability behavior have been investigated [25][26][27][28][29][30][31]. These contributions primarily analyze the influence of a single factor on the dynamic change in reservoir permeability based on laboratory testing. Although much work has been done on the influencing factors of permeability with remarkable results, there is still a lack of quantitative research on the comprehensive impact of various factors on reservoir permeability. In summary, the advanced characterization of dynamic change in coal reservoir permeability and its influence are still urgent problems to be solved.
Presently, studies on permeability dynamics are primarily used in CBM numerical simulation, including history matching, productivity prediction, and well pattern optimization. Few studies have applied permeability dynamic behavior to reservoir classification. Meanwhile, most previous studies on coal reservoir classification primarily focus on initial geological conditions such as geologic structure, hydrodynamic conditions, and initial physical properties of coal reservoirs [2,[31][32][33][34][35][36][37]. These works have made a significant contribution to the early exploration of CBM. However, coal reservoir classification needs to remain relevant and modern. On the one hand, the initial geological conditions may not be the key factor determining the productivity potential of wells with CBM development for decades, and reservoir classification based on static geological conditions is often inconsistent with gas production; on the other hand, the dynamic behavior of coal reservoir physical properties has become more and more important for CBM production. Therefore, it is a very meaningful work to classify coal reservoirs based on the dynamic behavior of reservoir permeability.
In this contribution, the effects of effective stress and sorption-induced strain, which have been fully considered relative to the permeability of a high-rank coal reservoir in China [12,38], were applied to examine the permeability dynamic characterization in the single-phase water flow stage and gas desorption stage. Some demarcating pressures and corresponding permeabilities were proposed to quantitatively characterize the dynamic permeability in various stages. Furthermore, the influence of geologic parameters on these pressures was analyzed in detail, and the type of coal reservoir was determined based on the permeability dynamic behavior. This method is applied to classify reservoirs in the Shizhuangnan Block, and the results are consistent with the actual production, which show the validity of this method and its important role in reservoir classification and CBM production.

Reservoir Properties in the Single-Phase Water Flow Stage
For unsaturated coal reservoirs, when reservoir pressure is greater than the critical desorption pressure, CBM cannot desorb from the coal matrix, and pores are filled with coal seam water, which defines CBM production in the single-phase water flow stage. Coal reservoirs are only affected by effective stress, and the mathematical model for the permeability change induced in the in situ coal condition was established (Equation (1)) [12,39]. The permeability damage rate is defined as the slope of a dynamic permeability curve, which can be obtained from the first derivative of the dynamic permeability formula. Its physical meaning refers to the change degree of reservoir permeability caused by the drop of reservoir pressure in the single-phase water flow stage. This parameter is used to reasonably characterize the effect of a drop in reservoir pressure on permeability (Equation (2)).
where k 1 is the dynamic permeability in the single-phase water flow stage, mD; k i is the initial permeability, mD; P i is the initial reservoir pressure, MPa; P is the reservoir pressure, MPa; C f is the cleat-volume compressibility, MPa −1 ; and v is Poisson's ratio, dimensionless. In order to analyze the trend of the curve more intuitively, the initial values of each parameter are set to Reservoir permeability and its damage rate decrease monotonously with decreasing pressure (Figure 1). The effect of effective stress on coal reservoir permeability is significant at high reservoir pressure. Then, the permeability damage rate decreases with decreasing reservoir pressure, indicating that the effect of effective stress on permeability gradually decreases. When the reservoir pressure drops to a low level, the permeability damage rate is almost zero. Since coal reservoirs with permeability less than 0.1 mD are impermeable [40], the invulnerable pressure (P n ) is defined as the pressure at which the permeability damage rate is equal to 0.1 mD/MPa. When reservoir pressure is less than the invulnerable pressure, the effective stress has a minimal effect and can be ignored (Equation (3)).
Energies 2020, 13, x FOR PEER REVIEW  3 of 17 For unsaturated coal reservoirs, when reservoir pressure is greater than the critical desorption pressure, CBM cannot desorb from the coal matrix, and pores are filled with coal seam water, which defines CBM production in the single-phase water flow stage. Coal reservoirs are only affected by effective stress, and the mathematical model for the permeability change induced in the in situ coal condition was established (Equation (1)) [12,39]. The permeability damage rate is defined as the slope of a dynamic permeability curve, which can be obtained from the first derivative of the dynamic permeability formula. Its physical meaning refers to the change degree of reservoir permeability caused by the drop of reservoir pressure in the single-phase water flow stage. This parameter is used to reasonably characterize the effect of a drop in reservoir pressure on permeability (Equation (2)).
where is the dynamic permeability in the single-phase water flow stage, mD; is the initial permeability, mD; is the initial reservoir pressure, MPa; is the reservoir pressure, MPa; is the cleat-volume compressibility, MPa -1 ; and is Poisson's ratio, dimensionless. In order to analyze the trend of the curve more intuitively, the initial values of each parameter are set to = 4, = 1, = 0.3, = 4.
Reservoir permeability and its damage rate decrease monotonously with decreasing pressure (Figure 1). The effect of effective stress on coal reservoir permeability is significant at high reservoir pressure. Then, the permeability damage rate decreases with decreasing reservoir pressure, indicating that the effect of effective stress on permeability gradually decreases. When the reservoir pressure drops to a low level, the permeability damage rate is almost zero. Since coal reservoirs with permeability less than 0.1 mD are impermeable [40], the invulnerable pressure ( ) is defined as the pressure at which the permeability damage rate is equal to 0.1 mD/MPa. When reservoir pressure is less than the invulnerable pressure, the effective stress has a minimal effect and can be ignored (Equation (3)). In order to quantitatively characterize the dynamic change in reservoir permeability in the single-phase water flow stage, the curvature of the dynamic permeability ( ) was calculated (Equation (4)). Generally, the curvature positively correlates with the permeability damage rate [41]. The shape of the curve is similar to a parabola; that is, the curvature of the dynamic permeability curve rises first and then decreases with the depletion of pressure, indicating that the permeability damage rate changes from rapid decrease to a slow decrease ( Figure 2). The turning pressure ( ) is In order to quantitatively characterize the dynamic change in reservoir permeability in the single-phase water flow stage, the curvature of the dynamic permeability (k q ) was calculated Energies 2020, 13, 644 4 of 17 (Equation (4)). Generally, the curvature positively correlates with the permeability damage rate [41]. The shape of the curve is similar to a parabola; that is, the curvature of the dynamic permeability curve rises first and then decreases with the depletion of pressure, indicating that the permeability damage rate changes from rapid decrease to a slow decrease ( Figure 2). The turning pressure (P t ) is defined as the reservoir pressure corresponding to the maximum curvature, and its value is the result of the first derivative of k q being zero (Equations (5) and (6)). (4) where k and k are the first and second derivatives of the dynamic change formula of permeability, respectively. k q is the curvature of the dynamic curve of permeability.
Energies 2020, 13, x FOR PEER REVIEW 4 of 17 defined as the reservoir pressure corresponding to the maximum curvature, and its value is the result of the first derivative of being zero (Equations (5) and (6)).
where and are the first and second derivatives of the dynamic change formula of permeability, respectively.
is the curvature of the dynamic curve of permeability. The curve also has two extrema, and the corresponding pressures are defined as the vulnerable pressure ( ) and alleviate pressure ( ), respectively. The extrema can be obtained by calculating = 0 (Equation (7)).
In summary, according to coal permeability dynamic characterization in the single-phase water flow stage, the coal permeability undergoes five stages during depressurization: the vulnerable stage, vulnerable transition stage, alleviative transition stage, alleviative stage, and invulnerable stage. The vulnerable pressure, turning pressure, alleviate pressure, and invulnerable pressure are the demarcating pressures corresponding to the five stages. Based on these stages, the coal reservoir can be further classified such that if > , the reservoir is a vulnerable reservoir; if < < , The k q curve also has two extrema, and the corresponding pressures are defined as the vulnerable pressure (P v ) and alleviate pressure (P a ), respectively. The extrema can be obtained by calculating k q = 0 (Equation (7)).
Energies 2020, 13, 644 5 of 17 In summary, according to coal permeability dynamic characterization in the single-phase water flow stage, the coal permeability undergoes five stages during depressurization: the vulnerable stage, vulnerable transition stage, alleviative transition stage, alleviative stage, and invulnerable stage. The vulnerable pressure, turning pressure, alleviate pressure, and invulnerable pressure are the demarcating pressures corresponding to the five stages. Based on these stages, the coal reservoir can be further classified such that if P i > P t , the reservoir is a vulnerable reservoir; if P n < P i < P t , the reservoir is an alleviative reservoir; if P i < P n , the reservoir is an invulnerable reservoir ( Figure 3).
Energies 2020, 13, x FOR PEER REVIEW 5 of 17 the reservoir is an alleviative reservoir; if < , the reservoir is an invulnerable reservoir ( Figure   3).

Reservoir Properties in the CBM Desorption Stage
For saturated coal reservoirs or unsaturated coal reservoirs when reservoir pressure drops below the critical desorption pressure, CBM desorbs from the coal matrix with decreasing pressure, and pores are filled with not only coal seam water but also desorbed gas. At this time, CBM development is in the gas desorption stage, and reservoir permeability is not only damaged by the effective stress effect, but also restored by the matrix shrinkage and gas slippage effects. The effect of gas slippage on permeability is very small (the permeability increase due to the gas slippage effect is 1/10 th of the matrix shrinkage effect), and the gas slippage effect is commonly significant at low pressure (<1 MPa) [12,13]; hence, the gas slippage effect on permeability is not considered in this study. The dynamic change formula for reservoir permeability in the CBM desorption stage is as follows [42,43]: where is the dynamic permeability in the CBM desorption stage, mD; is the critical desorption pressure, MPa; is the reservoir permeability corresponding to the critical desorption pressure, mD; is the specific surface area, m 2 /kg; is the matrix particle radius, m; is the initial porosity, dimensionless; ( ) is the equivalent matrix particle radius, m; is the Langmuir volume, m 3 /t; is the Langmuir pressure, MPa; and is the coal density, g/cm 3 . The relationship between permeability and reservoir pressure is parabolic ( Figure 4). Coal reservoir permeability changes significantly during gas production, commonly initially decreasing due to the dominate effective stress but then increasing due to the dominate matrix shrinkage effect as reservoir pressure and gas content are drawn down. This turning point is defined as the rebound pressure ( ) [21,22]. is given by the solution to the equation = 0, which can be solved using Matlab (Equation (13)).

Reservoir Properties in the CBM Desorption Stage
For saturated coal reservoirs or unsaturated coal reservoirs when reservoir pressure drops below the critical desorption pressure, CBM desorbs from the coal matrix with decreasing pressure, and pores are filled with not only coal seam water but also desorbed gas. At this time, CBM development is in the gas desorption stage, and reservoir permeability is not only damaged by the effective stress effect, but also restored by the matrix shrinkage and gas slippage effects. The effect of gas slippage on permeability is very small (the permeability increase due to the gas slippage effect is 1/10 th of the matrix shrinkage effect), and the gas slippage effect is commonly significant at low pressure (<1 MPa) [12,13]; hence, the gas slippage effect on permeability is not considered in this study. The dynamic change formula for reservoir permeability in the CBM desorption stage is as follows [42,43]: where k 2 is the dynamic permeability in the CBM desorption stage, mD; P cd is the critical desorption pressure, MPa; k cd is the reservoir permeability corresponding to the critical desorption pressure, mD; S v is the specific surface area, m 2 /kg; r i is the matrix particle radius, m; ϕ i is the initial porosity, dimensionless; R(P) is the equivalent matrix particle radius, m; V L is the Langmuir volume, m 3 /t; P L is the Langmuir pressure, MPa; and ρ is the coal density, g/cm 3 . The relationship between permeability and reservoir pressure is parabolic (Figure 4). Coal reservoir permeability changes significantly during gas production, commonly initially decreasing due to the dominate effective stress but then increasing due to the dominate matrix shrinkage effect as reservoir pressure and gas content are drawn down. This turning point is defined as the rebound Energies 2020, 13, 644 6 of 17 pressure (P rb ) [21,22]. P rb is given by the solution to the equation k 2 = 0, which can be solved using Matlab (Equation (13)).
Energies 2020, 13, x FOR PEER REVIEW 6 of 17 rebound in the CBM desorption stage. The effective stress effect on permeability is greater than the matrix shrinkage effect in the early stage, and the matrix shrinkage effect is greater than the effective stress effect in the later stage. However, if the rebound pressure is less than or equal to 0, the effective stress is greater than the matrix shrinkage effect throughout pressure depletion, and the permeability decreases continually. If the rebound pressure is greater than or equal to , the matrix shrinkage effect is greater than the effective stress effect throughout pressure depletion, and the permeability increases continually.

Reservoir Properties in the Whole Production Stage
For unsaturated coal reservoirs, dynamic changes in reservoir permeability during the complete drainage process include the dynamic change in the single-phase water flow stage and the CBM desorption stage. Therefore, the reservoir permeability dynamic change formula is given by Equation (14). In order to quantitatively analyze the permeability recovery, the recovery pressure ( ) is defined as the reservoir pressure, where reservoir permeability is recovered to initial permeability (Equation (15)). If the recovery pressure is greater than 0, the coal reservoir is a recoverable reservoir. On the contrary, if the recovery pressure is less than or equal to 0, the coal reservoir is an unrecoverable reservoir.
In summary, based on the dynamic change in the permeability during CBM well development, coal reservoirs can be classified as vulnerable reservoirs, alleviative reservoirs, and invulnerable reservoirs in the single-phase water flow stage, as well as recoverable reservoirs and unrecoverable reservoirs in the CBM desorption stage ( Figure 5). Reservoir classification is helpful for studying the geological conditions of CBM development and optimizing CBM well locations. The first derivative of the reservoir permeability equation in the CBM desorption stage increases monotonically with increasing reservoir pressure, indicating that the rebound pressure has at most one solution [13]. If the rebound pressure is between 0 − P cd , reservoir permeability will rebound in the CBM desorption stage. The effective stress effect on permeability is greater than the matrix shrinkage effect in the early stage, and the matrix shrinkage effect is greater than the effective stress effect in the later stage. However, if the rebound pressure is less than or equal to 0, the effective stress is greater than the matrix shrinkage effect throughout pressure depletion, and the permeability decreases continually. If the rebound pressure is greater than or equal to P cd , the matrix shrinkage effect is greater than the effective stress effect throughout pressure depletion, and the permeability increases continually.

Reservoir Properties in the Whole Production Stage
For unsaturated coal reservoirs, dynamic changes in reservoir permeability during the complete drainage process include the dynamic change in the single-phase water flow stage and the CBM desorption stage. Therefore, the reservoir permeability dynamic change formula is given by Equation (14). In order to quantitatively analyze the permeability recovery, the recovery pressure (P rc ) is defined as the reservoir pressure, where reservoir permeability is recovered to initial permeability (Equation (15)). If the recovery pressure is greater than 0, the coal reservoir is a recoverable reservoir. On the contrary, if the recovery pressure is less than or equal to 0, the coal reservoir is an unrecoverable reservoir. e In summary, based on the dynamic change in the permeability during CBM well development, coal reservoirs can be classified as vulnerable reservoirs, alleviative reservoirs, and invulnerable

Advanced Characterization of the Four Demarcating Pressures in the Single-Phase Water Stage
The above research shows that the coal permeability undergoes five stages in the single-phase water flow stage: vulnerable stage, vulnerable transition stage, alleviative transition stage, alleviative stage, and invulnerable stage. Vulnerable pressure, turning pressure, alleviate pressure, and invulnerable pressure are the demarcating pressures corresponding to the five stages. Therefore, the physical properties of these demarcating pressures are further examined.
Firstly, the permeability may not experience all five of these stages in the single-phase water stage depending on the coal reservoir properties. For example, if > , permeability can experience from the vulnerable stage. It is worth noting that if > is to be matched, In other words, if < < , the permeability may not experience the vulnerable stage, but it will experience the vulnerable transition stage, and √ < ( ) < √ is matched in this condition. Similarly, the conditions for the other stages are shown in Table 1.
The permeability and permeability damage rate at the demarcating pressures can be obtained by substituting the formula of the demarcating pressure into the corresponding dynamic change formulas. The reservoir permeability of the demarcating pressures is independent of the initial permeability but is related to the reservoir properties, specifically, cleat-volume compressibility and Poisson's ratio. The reservoir permeability damage rate corresponding to the demarcating pressures is constant. The lower the pressure, the lower the permeability damage rate.

Advanced Characterization of the Four Demarcating Pressures in the Single-Phase Water Stage
The above research shows that the coal permeability undergoes five stages in the single-phase water flow stage: vulnerable stage, vulnerable transition stage, alleviative transition stage, alleviative stage, and invulnerable stage. Vulnerable pressure, turning pressure, alleviate pressure, and invulnerable pressure are the demarcating pressures corresponding to the five stages. Therefore, the physical properties of these demarcating pressures are further examined.
Firstly, the permeability may not experience all five of these stages in the single-phase water stage depending on the coal reservoir properties. For example, if P i > P m , permeability can experience from the vulnerable stage. It is worth noting that if P i > P m is to be matched, . In other words, if P z < P i < P m , the permeability may not experience the vulnerable stage, but it will experience the vulnerable transition stage, and is matched in this condition. Similarly, the conditions for the other stages are shown in Table 1.
The permeability and permeability damage rate at the demarcating pressures can be obtained by substituting the formula of the demarcating pressure into the corresponding dynamic change formulas. The reservoir permeability of the demarcating pressures is independent of the initial permeability but is related to the reservoir properties, specifically, cleat-volume compressibility and Poisson's ratio. The reservoir permeability damage rate corresponding to the demarcating pressures is constant. The lower the pressure, the lower the permeability damage rate. Table 1. Characteristics of reservoir permeability in single-phase water flow stage.

Demarcating Pressure and Stages Conditions k (mD) k (mD/MPa)
Vulnerable stage Alleviative transition stage (P a < P i < P t ) Advanced Characterization of P rb and P rc P rb and P rc can be calculated based on the coal reservoir parameters using Equations (13) and (15). In order to study the influence of geological factors on P rb and P rc , an orthogonal test was designed, which is one of the most effective and time-saving methods for studies involving multiple variables to determine which factors (or variables) primarily influence the properties of the target product [44]. The basic feature of the orthogonal test involves replacing the comprehensive experiments with some characteristic experiments and examining the comprehensive experiments by analyzing some experimental results [45]. Through the analysis of a part of the experimental results, the comprehensive experiment is studied, and the optimal level combination is determined.
Target geological parameters include cleat-volume compressibility, Poisson's ratio, critical adsorption pressure, porosity, Langmuir volume, Langmuir pressure, and coal density ( Table 2). The Langmuir volume, Langmuir pressure, and critical desorption pressure can reflect the coal reservoir adsorption performance, while the cleat-volume compressibility, Poisson's ratio, porosity, and density can reflect the physical properties of coal reservoirs. The parameters are independent of each other, and the joint collocation between the parameters has little effect on the experimental results; hence, the interaction between the parameters is ignored. Finally, an orthogonal test is designed according to the standard orthogonal array L9 (3 7 ). It is worth noting that the initial reservoir pressure is equal to the critical desorption pressure in Equation (13) to calculate the rebound pressure. However, when the properties of P rc are studied, the dynamic change of permeability in the single-phase water flow stage cannot be ignored; so, P i ≥ P cd and P i = 4 MPa in this test (Table 3).   The influence degree of geological parameters on the results can be intuitively analyzed ( Table 4). The influence degree of geological parameters on the rebound pressure from large to small are critical desorption pressure, cleat-volume compressibility, porosity, and Langmuir volume. The influence degree of geological parameters on the ratio of rebound pressure to critical desorption pressure from large to small are cleat-volume compressibility, porosity, Langmuir volume, and critical desorption pressure. The influence degree of geological parameters on the recovery pressure from large to small are critical desorption pressure, density, cleat-volume compressibility, and Langmuir volume. In addition, cleat-volume compressibility, porosity, and Langmuir pressure are inversely proportional to rebound pressure and recovery pressure. In contrast, at large critical desorption pressures, Langmuir volume and coal density are conducive to the rebound and recovery permeability. The critical desorption pressure positively correlates with rebound pressure and recovery pressure but negatively correlates with the ratio of rebound pressure to critical desorption pressure.
The above experimental results can reasonably explain the problems faced in actual CBM production. (1) An increased cleat-volume compressibility indicates that a coal reservoir is more easily compacted, resulting in coal reservoirs being more susceptible to damage. (2) The larger the porosity, the greater the water content of the coal seam. Therefore, the reason for low gas production in high water-yield wells is closely related to the difficult recovery of permeability. (3) A high Langmuir volume indicates that the gas content in the coal seam is high, so large CMB desorption is beneficial to the easy recovery of permeability. (4) The higher the critical desorption pressure, the greater the rebound pressure and recovery pressure, but the smaller the ratio of rebound pressure to critical desorption pressure, indicating that there is a relatively large gap between the rebound pressure and critical desorption pressure for coal seams with high critical desorption pressure. Therefore, gas production should be controlled to prevent gas lock in the early stage of production. For coal seams with low critical desorption pressure, permeability will gradually increase after gas production.

Regional Geology
The Shizhuangnan Block is located in the southern region of the Qinshui Basin, Shanxi Province (Figure 6a,b). At present, the minefield primarily produces CBM in the No.3 coal reservoir, which primarily consists of anthracite with a vitrinite reflectance (R o max ) ranging from 2.92% to 3.02%. The No.3 coal seam has a simple structure, and it is generally deeper in the northern and the central regions, while it is shallower in the southern and the eastern regions, with depths ranging from 451 to 1030 m. There are more than 10 secondary high-angle normal faults in the northern study area near the Sitou Faults, which cut the coal strata and limit CBM production. The coal seam is stable, and the thickness ranges from 4.45 to 8.75 m, with an average of 6.35 m. Coal gas content ranges from 7 to 21 m 3 /t. The gas content in the west and north is greater than in the east and south [35,36]. There are 937 wells in the entire field, which are primarily distributed in the north, central, and southwest regions, and there are no production wells in the east and southeast regions. The CBM wells have been produced for more than four years and more than eight years in the central region. The distribution map of the maximum gas production of CBM wells is drawn using the production data from all wells. The western and central regions are the primary gas-producing areas, while most of the CBM wells in the north are low gas-yield wells (Figure 9a).
There are 937 wells in the entire field, which are primarily distributed in the north, central, and southwest regions, and there are no production wells in the east and southeast regions. The CBM wells have been produced for more than four years and more than eight years in the central region. The distribution map of the maximum gas production of CBM wells is drawn using the production data from all wells. The western and central regions are the primary gas-producing areas, while most of the CBM wells in the north are low gas-yield wells (Figure 9a).

Characteristic Analysis of the Target Wells
The reservoir dynamic characterization of two CBM wells located in different areas was analyzed. Using basic experiment, logging data analysis, and simulation historical matching, the actual geologic parameters of each well are obtained. v is equal to 0.3 (experimental data), S v is equal to 1295 m 2 /kg (cite from [12]), and ρ is equal to 1.4 g/cm 3 (logging data) with small differences, while other geological parameters vary significantly (Table 5). Based on the above methods, six demarcating pressures for each well were calculated, and the recovery factor of CBM wells (Rc) was obtained by combining the numerical simulation and volumetric method. Additionally, a reservoir pressure prediction model based on the material balance equation (MBE) for coal reservoirs was used, which considers the self-regulating effects of coal reservoirs and the dynamic change of the equivalent drainage area, such that the actual reservoir average pressure and permeability can be calculated based on the reservoir conditions and actual production data [38]. The initial reservoir pressure of well Z26 is greater than the invulnerable pressure, while the rebound and recovery pressure is 2.7 MPa and 2.3 MPa, respectively ( Table 5), indicating that permeability transitions from the alleviative stage and gradually increases in the later period. The reservoir at well Z26 is slightly damaged in the single-phase water flow stage and recoverable in the gas desorption stage. Although low initial permeability is not conductive to CBM production, the effective depressurization of a coal reservoir benefits from artificial fractures and reasonable drainage strategies (Figure 7a). Reservoir permeability can be increased to more than 1 mD by the matrix shrinkage effect (Figure 8a). This type of reservoir has excellent development potential, and artificial construction should receive more attention. The correlation between drainage strategies and reservoir physical properties in the early production cycle as well as the recovery factor and residual gas content in the late stage should also be thoroughly examined.

Coal Reservoir Classification
In order to study the types of coal reservoirs in the Shizhuangnan Block, 26 wells were selected to further analyze the dynamic change in coal reservoir permeability (Figure 6c). These wells are distributed throughout the region and include both high and low gas production wells. The specific parameters of these wells are shown in Table 6.

Coal Reservoir Classification
In order to study the types of coal reservoirs in the Shizhuangnan Block, 26 wells were selected to further analyze the dynamic change in coal reservoir permeability (Figure 6c). These wells are distributed throughout the region and include both high and low gas production wells. The specific parameters of these wells are shown in Table 6. The reservoir at well T57 is invulnerable in the single-phase water flow stage and unrecoverable in the gas desorption stage. Additionally, due to the low reservoir permeability, the reservoir pressure cannot be effectively dropped (Figure 7b), so desorption of the coal reservoir is primarily concentrated in the fracturing area near the well. More importantly, reservoir permeability cannot increase significantly throughout the field life (Figure 8b), so this type of reservoir has little development potential.

Coal Reservoir Classification
In order to study the types of coal reservoirs in the Shizhuangnan Block, 26 wells were selected to further analyze the dynamic change in coal reservoir permeability (Figure 6c). These wells are distributed throughout the region and include both high and low gas production wells. The specific parameters of these wells are shown in Table 6. Table 6. Geological parameters and calculation results of the target wells in the Shizhuangnan Block. The coal reservoirs in this block were classified according to the methods described above (Figure 9). For the damaged nature of permeability, due to the influence of faults and the deep burial depth, the coal seam in the northern region is mostly tectonic coal with low permeability, which has led to an invulnerable permeability; hence, the coal reservoirs in the northern region is classified as invulnerable. In contrast, the coal reservoirs in the central and southwestern region are shallow and less affected by tectonism. The permeability of coal seam is relatively high and slightly damaged (Figure 9b). For the recovery nature of permeability, the dynamic permeability in the central region can restore to the initial permeability during pressure depletion due to its characterization of high gas content and low water content, so the reservoir is recoverable. However, because of the low gas content and high water content in the northern coal reservoirs, the permeability cannot be restored to its initial stage, so the reservoir is unrecoverable (Figure 9c). less affected by tectonism. The permeability of coal seam is relatively high and slightly damaged (Figure 9b). For the recovery nature of permeability, the dynamic permeability in the central region can restore to the initial permeability during pressure depletion due to its characterization of high gas content and low water content, so the reservoir is recoverable. However, because of the low gas content and high water content in the northern coal reservoirs, the permeability cannot be restored to its initial stage, so the reservoir is unrecoverable (Figure 9c).  Combining the reservoir type and the gas production characterization of CBM wells, it is concluded that although coal reservoir permeability in the central region decreases gradually during the development process, reservoir permeability can return to the initial level or even greater after the CBM desorption. Therefore, the CBM wells in the central region have high gas production and excellent development potential. On the contrary, reservoir permeability in the northern region is invulnerable, and it cannot increase after CBM desorption. In addition, the reservoir is controlled by low initial permeability and natural faults, so the gas production is low (Figure 9a).
A high-rank coal reservoir in the Shizhuangnan Block is employed as an example to analyze in detail. However, this model is universally applicable, regardless of the conditions of any coal rank.
The key is to obtain the accurate geological parameters of coal reservoirs as shown in the paper, and researchers can use this model to analyze the types of coal reservoirs in different blocks with complete geological parameters.

Conclusions
The permeability dynamic characterization of coal reservoirs in the single-phase flow stage and the gas desorption stage are comprehensively analyzed, and coal reservoir geologic conditions in the Shizhuangnan Block are examined.
(1) The damage of effective stress on reservoir permeability in the single-water flow stage is analyzed.
The coal permeability damage rate can be divided into the vulnerable stage, vulnerable transition stage, alleviative transition stage, alleviative stage, and invulnerable stage based on the law of decreasing permeability. Vulnerable pressure, turning pressure, alleviate pressure, and invulnerable pressure are defined as the demarcating pressures. The permeability damage rate of these demarcating pressures is constant, and their permeability is only related to the reservoir physical properties, but not to initial reservoir permeability. (2) The influence of geologic factors on the rebound pressure and recovery pressure is quantitatively analyzed. The critical desorption pressure has the greatest influence on rebound pressure, followed by cleat-volume compressibility, porosity, and Langmuir volume. The influence order of recovery pressure with the geological parameters is as follows: the critical desorption pressure, density, formation compressibility, and Langmuir volume. Cleat-volume compressibility, porosity, and Langmuir pressure are inversely proportional to rebound pressure and recovery pressure. At large critical desorption pressure, Langmuir volume and coal density are conducive to permeability rebound and recovery. The higher the critical desorption pressure, the smaller the ratio of rebound pressure to critical desorption pressure, indicating that gas lock should be prevented in coal reservoirs with high gas content. (3) Coal reservoirs are classified according to the coal reservoir dynamic characteristics. When the initial reservoir pressure is greater than the critical desorption pressure, coal reservoirs can be classified as vulnerable, alleviative, or invulnerable based on the effect of effective stress. When the reservoir pressure is less than the critical desorption pressure, coal reservoirs can be classified as recoverable and unrecoverable reservoirs on the permeability recovery properties. In addition, this study on the types of coal reservoirs in the Shizhuangnan Block found that there is a significant correspondence between the types of coal reservoirs and CBM well gas production.