Next Article in Journal
Evaluation of Temperature on the Methane Hydrates Formation Process Using Sodium Surfactin and Rhamnolipids
Previous Article in Journal
Potential of Offshore Wind Energy in Malaysia: An Investigation into Wind and Bathymetry Conditions and Site Selection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characteristics of Permeability Evolution and Pore Structure of Coal with High Gas

School of Mechanics & Civil Engineering, China University of Mining & Technology, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Submission received: 22 October 2023 / Revised: 1 December 2023 / Accepted: 12 December 2023 / Published: 21 December 2023
(This article belongs to the Section H: Geo-Energy)

Abstract

:
To study the influence of gas pressure on coal permeability evolution, we conducted experiments on coal samples from the No. 9 coal seam in Tangshan Coal Mine, Hebei Province, China. Different gas pressures (helium and nitrogen) were applied, and nitrogen-induced deformations were measured. We also analyzed the coal samples’ pore structure using mercury injection porosimetry, obtaining pore surface fractal dimensions. The increase in nitrogen pressure from 0.3 MPa to 3 MPa resulted in an elevation of adsorption strain from 0.168 × 10−3 to 1.076 × 10−3, with a gradual decrease observed in the extent of this increase. However, the permeability of coal samples initially decreased from 16.05 × 10−18 m2 to 4.91 × 10−18 m2 and subsequently rose to 5.69 × 10−18 m2. Helium showed similar trends to nitrogen, with average permeability 1.42–1.88 times higher under the same pressure. The lowest permeability occurred at 1.5 MPa for helium and 2.5 MPa for nitrogen. Gas absorptivity plays a crucial role in coal permeability evolution. Additionally, we observed coal’s compressibility to be 7.2 × 10−11 m2/N and corrected porosity to be 53.8%, considering matrix compression. Seepage pores larger than 100 nm accounted for 80.4% of the total pore volume, facilitating gas seepage. Surface fractal dimension Ds1 correlated positively with micropore volume, while Ds2 and Ds3 correlated negatively with pore volume and gas permeability.

1. Introduction

Natural gas is currently regarded as an efficient and clean energy source. As a kind of unconventional natural gas, coalbed methane (CBM) has received more and more attention, and it has gradually become an important emerging energy source.
Previous research and practice have shown that during CBM mining, permeability is the main factor affecting CBM production [1,2]. It is of great significance to study the evolution of permeability to improve CBM production. The main factors affecting permeability include the physical properties of coal, pore structure, reservoir depth, in situ stress, reservoir pressure, temperature, gas adsorption and desorption capacity, among others. During the production of coalbed methane, permeability varies over time and conditions during the mining process. In the early stages of coalbed methane production, the effective stress increased due to the decrease in fluid pressure, which resulted in a decrease in gas reservoir pressure [3]. When the reservoir pressure dropped below the critical desorption pressure, the desorption of methane in the adsorbed state caused the shrinkage of the coal matrix, the pore fissure opening increased, and the permeability increased [4,5]. When the fluid pressure dropped below 1 MPa, gas slippage increased permeability. Therefore, the main factors for permeability evolution include effective stress, desorption-induced shrinkage of the coal matrix, and gas slippage. In addition to the fact that reservoir pressure affects changes in permeability, the physical properties of coal will also affect permeability during coalbed methane production. Coal permeability is influenced by in situ factors like pore structure (cleavage distance, pore size, heterogeneity) and stress conditions [6]. Experimental studies have shown that as gas pressure decreases, permeability also decreases due to the increase in effective stress, and the effective stress and permeability have an exponential relationship [7,8]. McKee et al. [9] found that as the depth of a coal seam increases, porosity is negatively correlated with the effective stress, and permeability decreases following a negative exponential curve. Seidle et al. [10] established a model that explained the exponential relationship between permeability and effective stress. Previous studies indicate that not only effective stress but also adsorption and desorption affect permeability [11]. In the experimental study, comparing permeability in coal with adsorbed and non-adsorbed gas revealed how gas properties affect coal adsorption/desorption strain [12,13]. Shi and Durucan [3] proposed a series of permeability models to describe the dynamic changes in permeability. Connell [14] established a linear relationship between adsorption amount and strain in his model.
Accurate identification of coal structure is of great significance to CBM exploration and development [15,16,17]. As a porous material, coal rock contains abundant pores and fractures. Pore and fracture structure affect the permeability and deformation of coal. The methods used to measure the pore structure of coal mainly include the mercury injection method, low-temperature liquid nitrogen adsorption method, low-temperature carbon dioxide adsorption method, small angle X-ray scattering method, atomic force microscope, scanning electron microscope, and other physical and photoelectromagnetic methods. The pore structure is studied by the high-pressure mercury injection method in this paper [18,19,20,21,22,23,24,25,26,27]. The pores measured, by pore diameter d, are divided into the micropore (d < 10 nm), transition pore (10 nm < d < 100 nm), mesopore (100 nm < d < 1 μm), macropore (1 μm < d < 10 μm) and micro fracture (d > 10 μm) [28]. The fractal dimension is used to quantificationally measure the features of pore surface and structure space, revealing the characteristics of space or structures. The volume fractal dimension is used to characterize the spatial heterogeneity of pore structure, and the surface fractal dimension is used to characterize the roughness of pore surface [20,29,30,31,32,33,34,35,36].
The relationship between pore properties and permeability has been extensively investigated by numerous researchers. Zhang et al. [37] combined the mercury injection experiment, computed tomography (CT) experiment, gas adsorption and desorption experiment, and the seepage experiment, and concluded that the degree of pore accumulation is the dominant factor affecting permeability change. Cheng [38] and Qin [39] obtained pore throat parameters at the scale of porous media with the mercury injection porosity (MIP) digital core reconstruction model and micro-computed tomography (micro-CT) images, respectively. The permeability contribution was evaluated according to pore connectivity, connectivity parameters, fractal dimension, and pore throat parameters. In addition, Hu et al. [40] established a multiple regression model of the influence of pore structure on the porosity and permeability of low-permeability sandstone reservoirs using multiple stepwise regression analysis; the multiple regression equation has high reliability and practicability.
This paper used the developed coal rock seepage–adsorption–deformation synchronous test device to test coal permeability evolution. Using nitrogen and helium as experimental gases, the coal permeability at different pressures was measured. The relationship between adsorption strain and permeability was analyzed, and the mechanism of adsorption on permeability was described in detail. In addition, the pore structure characteristics of coal samples were tested, and the roughness of the pore surface was measured by surface fractal dimension.

2. Experimental Tests

2.1. Coal Sample Preparation

Some coal blocks were collected from No. 9 coal seam of Tangshan Mine, situated in the northwestern region of the Kaiping coalfield in Hebei Province, China. The samples were taken at a depth of 900 m and exhibited an initial coalbed methane content of approximately 8 m3/t, classifying the mine as high gas. The measured average maximum vitrinite reflectance (R0,max%) of coal is 0.97%; considering this, No. 9 coal in Tangshan Mine belongs to the middle rank of coal. Coal samples were processed as cylinders with a diameter of approximately 2.5 cm and a height of approximately 5.0 cm, oriented along the parallel bedding direction. Following a 30 h dehumidification and drying process, the samples were stored in a dry environment. Cylindrical coal samples were used for simultaneous adsorption–permeation–deformation testing. Additionally, coal fragments were gathered for pore structure analysis using MIP tests.

2.2. Adsorption-Permeability Experiment

2.2.1. Experimental Equipment

The adsorption–permeability experiment was conducted using a “Permeation and Adsorption Simultaneous Test System”, which primarily consists of four integral components: a gas flowrate quantification module, a coal deformation measurement system, a gas pressure regulation mechanism, and an autonomous stress control unit, as depicted in Figure 1. This system is capable of simultaneously measuring the permeability, coal deformation, and gas adsorption capacity of samples. The permeability assessment was conducted in compliance with the Chinese National Standard SY/T 5336-2006 [41], pertaining to Core Analysis Methods. Given the low permeability of the Tangshan coal samples, the gas flow at the device’s outlet was minimal during the permeability test. Therefore, a gas collection device, indicated as element C in Figure 1, was employed to measure the outlet gas flow using the drainage method. The reference cylinder was connected to the stress sensor to monitor the inlet pressure and the gas pressure in the reference cylinder in real time.

2.2.2. Experimental Method

With the injection of gas into the coal samples, the gas adsorbs or fills the pores and cracks of the sample, thus causing deformation of the sample. This deformation is attributable to pore pressure in cases where the coal does not adsorb gas. Conversely, when the coal adsorbs gas, the deformation arises from both pore pressure and adsorption phenomena. To isolate the impact of gas adsorption on coal deformation, it is imperative to negate the deformation instigated solely by pore pressure. To this end, helium (He), a non-adsorptive gas [42], was initially injected into the coal sample, inducing strain solely through pore pressure influence. Subsequently, the sample was subjected to vacuum conditions before the introduction of an adsorbable gas. Given nitrogen’s (N2) adsorptive properties and its weaker adsorption capacity relative to methane (CH4), it serves as an effective proxy for CH4 in adsorption–permeability experiments on coal samples, in accordance with the principle of inclusivity. Moreover, the inert nature of N2 enhances the safety parameters of such experimental procedures. Therefore, N2 was injected into the coal as the adsorbable gas [43]. The experimental procedures are as follows:
(1)
The dried coal samples were set in the core holder along with two sets of axial–radial resistance strain gauges. With a flex flat cable (FFC), the strain gauges were joined to the axial–radial strain transducer.
(2)
To simulate the real stress of coal seam, the axial pressure and confining pressure were kept at 7 MPa, and a vacuum pump was used to vacuum the coal sample, and then the air tightness of the pipeline was checked.
(3)
Helium gas (He) was injected into the reference cylinder at a predetermined pressure from the gas source. Helium injection was performed with incremental pressure levels, including 0.3 MPa, 0.5 MPa, 1 MPa, 1.5 MPa, 2 MPa, 2.5 MPa, and 3 MPa. At each injection pressure, the pressure was kept constant while the outlet was sealed before injection. The valve connecting the gas source to the reference cylinder was then closed, while the outlet of the reference cylinder was opened to inject He into the core gripper. This process was maintained for 72 h to calibrate the pore volume of the coal sample.
(4)
The outlet valve of the core gripper was kept open to the atmosphere until all specified injection pressures were completed, while measuring the radial strain of the sample under different injection pressures. The flow rate of helium (He) was measured using a gas collection device to calculate the permeability of the coal sample. The gas flow volume was measured by a gas collecting device for a certain period of time, thereby calculating the flow rate, and calculating the coal sample permeability based on Darcy’s law, as demonstrated as Equation (1).
k = 2 Q p 0 μ L A p 1 2 p 2 2
where k is coal permeability (10−3 μm2); Q is the flow rate (mL/s) through coal sample; p0 is atmospheric pressure (take 0.1 MPa); μ is the hydrodynamic viscosity (Pa·s); L is the length of coal sample (cm); A is the cross-sectional area of coal sample (cm2) of the cylinder; and p1 and p2 are the gas pressures (MPa) at both ends of the inlet and outlet of coal sample.
(5)
We removed the gas pressure and closed the outlet, and then pumped the vacuum for more than 12 h.
(6)
Following the same procedure as the helium injection, nitrogen (N2) was introduced for a duration of 72 h to obtain the volume of gas adsorbed by the coal sample. The flow rate was then measured to determine the sample’s permeability.
(7)
For each experimental set, the gas pressure was varied, starting from 0.3 MPa and incrementally increasing to 3 MPa, with the entire experimental procedure repeated seven times. The effective stress of the coal samples varied with the increase in pore pressure, calculated by Equation (2).
σ e = ( 2 σ r + σ z ) / 3 p 1 + p 2 2
Here, σe represents effective stress, MPa; σr represents radial stress, MPa; σz represents axial stress, MPa. The inlet gas pressure and effective stress changes are listed in Table 1.

2.3. Mercury Injection Porosimetry Experiment

2.3.1. Experimental Method

The MIP experiment was conducted at Tsinghua University in China with the Auto Pore IV 9520 device, according to the standard GB/T 21650.1-2008 [44] of China. The coal specimens were cube-shaped with 1 cm long sides. The surface tension was 0.485 N/cm and the contact angle between the mercury and coal was 130°, respectively. Washburn’s equation [45] can describe the relationship between mercury intrusion pressure and the pore size of the coal sample, as shown as Equation (3).
d = 4 σ cos θ P
where d represents pore diameter injected, nm; σ represents the surface tension of mercury, MPa; θ represents the angle of coal surface and mercury, °; P represents mercury injection pressure, MPa.
The pore size decreased as the injection pressure increased. The highest intrusion pressure reached 413 MPa, and the corresponding pore diameter was 3 nm. The largest pore diameter injected was about 370 μm at a low mercury injection pressure. Furthermore, we can obtain the pore volume at different pore sizes according to the mercury quantity injected.

2.3.2. Surface Fractal Dimension Calculation Method

It is universally acknowledged that the surface of coal is highly uneven [46]. Surface fractal dimension Ds can quantitatively estimate the roughness and complex characteristics of coal surface. The value of Ds, generally speaking, is between 2 and 3. Coal surface tends to be smoothy when Ds is close to 2. Coal surface is more complex and rougher when Ds approaches 3. Based on mercury intrusion volume modified data, Zhang’s model was adopted in this work [35,36], calculated by Equation (4).
D s = d ln W n r n 2 d ln V n 1 / 3 r n
where Wn represents the accumulative surface energy, J; rn represents the smallest pore radius injected by mercury pressure Pi, nm; Vn represents cumulative mercury injecting volume, mL/g.
Where
W n = i = 1 i = n P i V i
The surface fractal dimension Ds is equal to the slope of the ln V n 1 / 3 r n vs. ln W n r n 2 as the ordinate.

3. Experimental Results and Analysis

3.1. Nitrogen Adsorption Strain of Coal Samples

Nitrogen is an adsorbent gas. The deformation of three coal samples at various gas pressure points is influenced not only by effective stress changes, but also by gas adsorption. The adsorption of nitrogen in coal causes the coal matrix to expand, resulting in a negative volumetric strain on the nitrogen-containing coal sample under identical experimental circumstances, indicating that each sample creates a considerable expansion strain. Helium has no absorbability, and its strain value reflects other stress conditions except adsorption effect. Therefore, the adsorption strain is equal to the difference between the strain values of nitrogen and helium, which is obtained from Equation (6). The adsorption strain results are listed in Table 2.
ε s = ε N 2 ε He
The adsorption strain of three coal samples rises as the gas pressure increases, as seen in Figure 2. The results demonstrate that when the gas pressure increases, the change rate of the coal sample strain slows, indicating that the adsorption strain is suppressed by effective stress [47,48]. To evaluate the law of coal adsorption gas, a Langmuir-type strain model was employed to fit the relationship between sorption-induced deformation and injection pressure. The strain induced by adsorption in coal is intrinsically linked to its gas adsorption capacity. Consequently, an equation analogous to the Langmuir equation was utilized in this study to characterize the adsorption deformation across varying coal ranks.
ε s = ε L p p + p L
where εL denotes the Langmuir strain, defined as the theoretical maximum axial, radial, or volumetric strain of coal, also in percentages; pL is the gas pressure at which the strain attains half of its theoretical maximum value, analogous to the Langmuir pressure, measured in MPa; p signifies the gas pressure, in MPa.
To fit the experimental data, Equation (7) is employed, and the fitting parameters are provided in Table 3. The solid points in Figure 2 represent experimental data, while the solid lines are fitting curves. The fitting correlative parameters are more than 0.99.

3.2. Coal Permeability under Different Gas Pressure

Based on the outlet flow rate measured at each specified gas pressure point, the permeability curves of three coal samples were calculated by Equation (1), as shown in Figure 3. The maximum permeability of the coal samples was determined to be 70 × 10−18 m2. Tangshan Coal Mine is clearly associated with a poor-permeability coal seam.
The confining pressure and axial pressure remain constant, and the effective stress decreases with the increase in gas pressure. Coal permeability changes with the increase in gas pressure, as described in Figure 3. As the gas pressure increases, the permeability decreases first and then increases. The turning point corresponding to lowest permeability is called the permeability rebound point [49]. The permeability rebound point of Tangshan coal samples is approximately equal to 1.0 MPa. The lowest permeabilities of the three samples are found to be 0.98 × 10−18 m2, 6.76 × 10−18 m2, and 15.0 × 10−18 m2, suggesting that, despite the fact that the coal samples are from the same mining location, the internal pore structure is different, resulting in considerable variances in permeability.
Previous studies held that effective stress was the primary factor influencing coal permeability variations. The relationship between permeability k and effective stress increments ∆σe is commonly expressed using Equation (8) [3].
k = k 0 exp - 3 C f σ e
where k0 is the initial permeability, (10−3 μm2); Cf is the pore compression factor.
We can observe from Equation (8) that permeability decreases as the effective stress increases. However, the permeability–effective stress curve recorded in our experiment is V-shaped, which contradicts Equation (8), suggesting that the effective stress is not the only factor affecting the permeability change. The permeability variation of helium is V-shaped too, indicating that the Klinkenberg effect exists.

3.3. Pore Structure Characteristics Detected by MIP

Mercury intrusion and extrusion curves of Tangshan coal samples are deployed in Figure 4. Mercury intrusion curves are in a reverse “S” shape. Mercury filled the interspaces of pores rapidly when the injecting pressure is less than 0.01 MPa, and the corresponding pore diameter is 12 μm. The mercury injection volume increases slowly as the mercury increases. The slopes of the cumulative mercury intrusion curves rise faster when the mercury injection pressure is greater than 100 MPa. The mercury injection curves also reveal that the pores with a diameter greater than 10 μm, and less than 100 nm, are more developed. There is a narrow hysteresis loop formed between the mercury intrusion and extrusion curves of each sample, indicating that seepage pores mainly consist of semi-closed pores. Meanwhile, a certain amount of open pores also exist [50]. Overall, Tangshan coal samples mainly consist of micropores, transport pores, and micro fractures. Some pore feature parameters of the coal samples are listed in Table 4. What is noteworthy is that the coal matrix produces elastic compression when the mercury injection pressure is greater than 20 MPa, resulting in the pore volume measured the pore in this range being larger than the actual volume [24,51]. An intrusion pressure of 20 MPa corresponds to a pore diameter of 60 nm in this work. So, the pore volume measured should be modified when the size less than 60 nm.

4. Discussion

4.1. The Relationship between Permeability and Adsorption Strain

For convenience, a dimensionless permeability ratio k/k0 was used to show the trend of permeability. As shown in Figure 5, the three coal samples all show a trend of k/k0 decreasing sharply at first and then rising slowly with the increase in gas pressure, and the volumetric strain increases gradually with the increase in gas pressure, and the permeability is basically negatively correlated with the volumetric strain. The permeability of the coal samples is affected by a variety of factors. The external confining pressure and axial pressure compress the coal sample, and the internal pores of the coal body are squeezed inward, resulting in a decrease in coal permeability [7]. After gas fills the coal, the gas pressure makes the pores expand, resulting in an increase in coal penetration and permeability. At the same time, the adsorption effect also makes the coal matrix swell [28,29,30]. There are competition mechanisms with several functions, which makes gives permeability change complex rules. When the external confining pressure and axial pressure remain unchanged, gas pressure keeps increasing, and the pores under the action of gas pressure will expand, and the permeability will increase with the increase in coal permeability. However, the experimental data show that the permeability basically shows a downward trend with the increase in gas pressure, indicating that the adsorption effect has a negative effect on permeability.
The volumetric strain of the coal sample increasing with the increase in adsorption capacity is shown in Figure 6. The more gas is adsorbed by the coal matrix, the greater the expansion deformation. The volumetric strain has a linear relationship with the adsorption capacity when the gas pressure is in low-pressure conditions (less than 2 MPa). The increase in volumetric strain slows down when the gas pressure is greater than 2 MPa.

4.2. Pore Size Distributions

The MIP test results for mesopores and macropores should have relatively high reliability. However, coal’s skeleton structure will be damaged under very high injection pressure. So, the results of mercury injection experiment are not suitable to be directly used as the results of nano-aperture, so it is necessary to revise the results of mercury injection experiment. In this paper, only mercury injection data were used for correction, which is the same as Li [26]. The matrix compressibility kc can be expressed as Equation (9) [52].
k c = d V m V m d P
Here, Vm represents the volume of the matrix skeleton, m3.
For a more practical approach to MIP data correction, it is assumed that the same effective compressibility can be applied as a component-independent parameter. In other words, the organic matrix, inorganic minerals, pores, and fissures without mercury filling have the same effective compressibility, which represents a bulk behavior. Once the pores and fractures are filled with mercury, they are not compressed. The compressed coal matrix volume exhibits the following relationship with the observed mercury volume (see Equation (10) [26].
V m c P i = V b u l k V o b s P i 1
Here, V m c P i represents the coal matrix compressive volume at pressure Pi, Vbulk represents the bulk volume of the coal sample, and Vobs (Pi−1) represents the observed mercury volume at pressure Pi−1.
According to coal matrix compressibility, the pore volume modified at Pi can be expressed by Equation (11).
V p o r e P i = V p o r e P i 1 + V o b s P i V o b s P i 1 V b u l k V o b s P i 1 × k c × P i P i 1
Here, V p o r e P i and V p o r e P i 1 represent the pore volume at Pi and Pi−1, respectively.
The compressibility of the coal matrix and other related parameters are listed in Table 5. The average compressibility of coal samples is 7.199 × 10−11 m2/N, and the mean volume correction rate is 53.85%. The larger the micropores’ volume, the smaller the matrix compressibility km.
The pore size distributions of the coal samples before being modified and after being modified are displayed by Figure 7. The dominant pores of Tangshan coal samples are micro fractures (pore diameter > 10,000 nm), accounting for 65.93%, which were convenient for coalbed gas seepage. The content of the other types of pores is low, which is not conducive to the storage and diffusion migration of coalbed methane. As shown in Figure 7, pore content fluctuates significantly at different pore size stages, indicating that pore size and pore volume are independent of each other.

4.3. Fractal Analysis on Pore Structure

As shown in Figure 8, the surface fractal dimension curves consist of three obvious regions. Namely, Ds1, micropore, and transport pore region (r < 50 nm); Ds2, mesopore, and macropore region (50 nm < r < 5000 nm); Ds3 and micro fracture region (r > 5000 nm). Compared to cement-based porous materials, Tangshan coal samples have distinct fractal characteristics in the mesopore and macropore regions [30,31]. The pore fractal features of the Tangshan coal samples in all regions are obvious (R2 > 0.99). The fractal dimensions of different regions are listed in Table 6. For the surface fractal dimensions of sample No. 1, Ds1 > Ds2 > Ds3, indicating that the pore surface with a smaller pore size is coarser and more complex. At the same time, the pore surface roughness of each pore size segment of sample No. 1 is greater than that of sample No. 2 and sample No. 3.
In Table 6, the adsorption pores make contributions to the corresponding pore surface fractal dimension Ds1, and the larger the adsorption pore volume is, the rougher the pore surface is and the larger the fractal dimension Ds1 is. Therefore, the surface fractal dimension Ds1 is positively correlated with the adsorption pore content. The more pore seepage content, the smoother the porous channel, the smaller pore surface dimensions Ds2 and Ds3, and the more conducive the sample is to gas seepage. Numerically, pore surface fractal dimensions Ds2 and Ds3 are negatively correlated to seepage pore content and permeability.

5. Conclusions

(1)
Under the condition that the axial and confining stresses are both 7 MPa, with the increase in gas pressure, the nitrogen adsorption strain aggrandizes and follows the laws of the Langmuir trend. The trend of permeability for helium is similar to a “V” shape, revealing that the permeability of helium is affected by effective stress and the slippage effect. The evolution of nitrogen permeability is similar to that of nitrogen. The permeability difference between them shows the influence of nitrogen adsorption effect on gas permeability, and adsorption effect reduces the permeability of coal samples.
(2)
The pore structure characteristics of coal samples with pore sizes ranging from 3 nm to 370 μm was investigated with a high-pressure mercury injection experiment. The coal matrix produces compressibility when the mercury injection pressure is greater than 20 MPa, and the mercury injection data are used merely for compressibility correction in this work. The dominant pores of Tangshan coal samples are micro-cracks, accounting for 66% of the total pore volume. The pores are mainly composed of semi-open pores with one end closed, which is convenient for gas seepage.
(3)
Surface fractal dimension Ds can reveal the roughness and complexion of the pore surface, and its value is between two and three. The larger Ds is, the rougher and more complex the pore surface is. The surface fractal dimension Ds1 is larger than Ds2 and Ds3, indicating that the surface roughness of adsorbed pores is higher than that of seepage pores. Ds1 is positively correlated with micropore volume content, while Ds2 and Ds3 are negatively correlated with pore volume content and gas permeability.

Author Contributions

Conceptualization, J.Z.; Data curation, Z.C., Y.Z., Q.W. and L.L.; formal analysis, T.S.; funding acquisition, J.Z.; investigation, T.L., Z.C., Y.Z. and Q.W.; methodology, T.S.; resources, J.Z.; software, T.L.; validation, T.S., T.L. and L.L.; visualization, Y.Z.; writing—original draft, T.S., Q.W. and L.L.; writing—review and editing, J.Z. and T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (52074297) and the Fundamental Research Funds for the Central Universities (2022YJSLJ02).

Data Availability Statement

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

Conflicts of Interest

We declare that we have no financial and personal relationships with other people or organizations that can inappropriately influence our work. In addition, there is no professional or other personal interest of any nature or kind in any product, service, and/or company that could be construed as influencing the position presented in, or the review of, the entitled manuscript.

Abbreviations

Across-sectional area of coal sample, cm2
Cfpore compression factor
dpore diameter injected, nm
Dssurface dimension
kcoal permeability, 10−3 μm2
k0initial permeability, 10−3 μm
kcmatrix compressibility coefficient, m2/N
Llength of coal sample, cm
psignifies the gas pressure, in MPa.
p0atmospheric pressure, MPa
p1gas pressures at both ends of the inlet of coal sample, MPa
p2gas pressures at both ends of the outlet of coal sample, MPa
pLLangmuir pressure, MPa
Pmercury injection pressure, MPa
Pimercury injected pressure at step i, MPa
Qflow rate through coal sample, mL/s
rnthe smallest pore radius injected by mercury pressure Pi, nm
Viincrement of mercury injected volume at pressure Pi, mL/g
Vmvolume of matrix skeleton, m3
V m c P i coal matrix compressive volume at pressure Pi, m3
Vbulkbulk volume of coal sample, m3
Vobs (Pi−1)represents the observed mercury volume at pressure Pi−1, m3
V p o r e P i pore volume at Pi, m3
V p o r e P i 1 pore volume at Pi−1, m3.
Wnaccumulated surface energy, J
εstrain
εsadsorption strain
εLLangmuir strain
εvvolumetric strain
σsurface tension of mercury, MPa
σeeffective stress, represents, MPa
σrradial stress, MPa
σzaxial stress, MPa
μhydrodynamic viscosity, Pa·s
θangle of coal surface and mercury, °

References

  1. Lv, Y.; Tang, D.; Xu, H.; Luo, H. Production characteristics and the key factors in high-rank coalbed methane fields: A case study on the Fanzhuang Block, Southern Qinshui Basin, China. Int. J. Coal Geol. 2012, 96–97, 93–108. [Google Scholar] [CrossRef]
  2. Tao, S.; Wang, Y.; Tang, D.; Xu, H.; Lv, Y.; He, W.; Li, Y. Dynamic variation effects of coal permeability during the coalbed methane development process in the Qinshui Basin, China. Int. J. Coal Geol. 2012, 93, 16–22. [Google Scholar] [CrossRef]
  3. Shi, J.Q.; Durucan, S. Drawdown Induced Changes in Permeability of Coalbeds: A New Interpretation of the Reservoir Response to Primary Recovery. Transp. Porous Media 2004, 56, 1–16. [Google Scholar] [CrossRef]
  4. Gray, I. Reservoir Engineering in Coal Seams: Part 1—The Physical Process of Gas Storage and Movement in Coal Seams. SPE Reserv. Eng. 1987, 2, 28–34. [Google Scholar] [CrossRef]
  5. Meng, Y.; Li, Z.; Lai, F. Experimental study on porosity and permeability of anthracite coal under different stresses. J. Pet. Sci. Eng. 2015, 133, 810–817. [Google Scholar] [CrossRef]
  6. Gentzis, T.; Deisman, N.; Chalaturnyk, R.J. Geomechanical properties and permeability of coals from the Foothills and Mountain regions of western Canada. Int. J. Coal Geol. 2007, 69, 153–164. [Google Scholar] [CrossRef]
  7. Lin, B.; Zhou, S. Experimental Investigation on the Permeability of the Coal Samples Containing Methane. J. China Univ. Min. Technol. 1987, 16, 21–28. [Google Scholar]
  8. Pan, Z.; Connell, L.D. Modelling permeability for coal reservoirs: A review of analytical models and testing data. Int. J. Coal Geol. 2012, 92, 1–44. [Google Scholar] [CrossRef]
  9. McKee, C.R.; Bumb, A.C.; Koenig, R.A. Stress-Dependent Permeability and Porosity of Coal and Other Geologic Formations. SPE Form. Eval. 1988, 3, 81–91. [Google Scholar] [CrossRef]
  10. Seidle, J.P.; Jeansonne, M.W.; Erickson, D.J. Application of Matchstick Geometry to Stress Dependent Permeability in Coals. In Proceedings of the SPE Rocky Mountain Regional Meeting, Casper, WY, USA, 18–21 May 1992; pp. 433–444. [Google Scholar]
  11. Connell, L.D.; Mazumder, S.; Sander, R.; Camilleri, M.; Pan, Z.; Heryanto, D. Laboratory characterisation of coal matrix shrinkage, cleat compressibility and the geomechanical properties determining reservoir permeability. Fuel 2016, 165, 499–512. [Google Scholar] [CrossRef]
  12. Schraufnagel, R.A.; Sawyer, W.K.; Paul, G.W. Development and Application of a 3D Coalbed Simulator Petroleum Society of Cim/Society of Petroleum Engineers Paper Development and Application of a 3D Coalbed Simulator By. In Proceedings of the Annual Technical Meeting, Calgary, AB, Canada, 10–13 June 1990. [Google Scholar]
  13. Harpalani, S.; Chen, G. Influence of gas production induced volumetric strain on permeability of coal. Geotech. Geol. Eng. 1997, 15, 303–325. [Google Scholar] [CrossRef]
  14. Connell, L. Coupled flow and geomechanical processes during gas production from coal seams. Int. J. Coal Geol. 2009, 79, 18–28. [Google Scholar] [CrossRef]
  15. Zhengguang, Z.; Yong, Q.; Geoff, W.; Hansen, S.; Zhenjiang, Y.; Jun, J.; Zhaobiao, Y. Evaluation of Coal Body Structures and Their Distributions by Geophysical Logging Methods: Case Study in the Laochang Block, Eastern Yunnan, China. Nat. Resour. Res. 2021, 30, 2225–2239. [Google Scholar] [CrossRef]
  16. Zhang, Z.; Qin, Y.; Yi, T.; You, Z.; Yang, Z. Pore Structure Characteristics of Coal and Their Geological Controlling Factors in Eastern Yunnan and Western Guizhou, China. ACS Omega 2020, 5, 19565–19578. [Google Scholar] [CrossRef]
  17. Chen, H.; Tian, W.; Chen, Z.; Zhang, Q.; Tao, S. Genesis of Coalbed Methane and Its Storage and Seepage Space in Baode Block, Eastern Ordos Basin. Energies 2022, 15, 81. [Google Scholar] [CrossRef]
  18. Li, X.; Chen, D.; Kang, Y.; Meng, X. Characterization of Pores and Fractures of Coal Based on CT Scan. Coal Geol. Explor. 2016, 44, 58–62+70. [Google Scholar] [CrossRef]
  19. Shadrin, A.V. Geophysical criterion of pre-outburst crack propagation in coal beds. J. Min. Sci. 2016, 52, 670–682. [Google Scholar] [CrossRef]
  20. Bai, P. Analysis of Pore Comprehensive Fractal Feature of Medium and Low Rank Coal in Xinjiang Region. Shaanxi Coal 2019, 38, 11–27. [Google Scholar]
  21. Groen, J.C.; Peffer, L.A.; Pérez-Ramirez, J. Pore size determination in modified micro- and mesoporous materials. Pitfalls and limitations in gas adsorption data analysis. Microporous Mesoporous Mater. 2003, 60, 1–17. [Google Scholar] [CrossRef]
  22. Yu, S.; Bo, J.; Ming, L.; Chenliang, H.; Shaochun, X. A review on pore-fractures in tectonically deformed coals. Fuel 2020, 278, 118248. [Google Scholar] [CrossRef]
  23. Wang, C.; Zhai, P.; Chen, Z.; Liu, J.; Wang, L.; Xie, J. Experimental study of coal matrix-cleat interaction under constant volume boundary condition. Int. J. Coal Geol. 2017, 181, 124–132. [Google Scholar] [CrossRef]
  24. Suuberg, E.M.; Deevi, S.C.; Yun, Y. Elastic behaviour of coals studied by mercury porosimetry. Fuel 1995, 74, 1522–1530. [Google Scholar] [CrossRef]
  25. Shao, P.; Wang, X.; Song, Y.; Li, Y. Study on the characteristics of matrix compressibility and its influence factors for different rank coals. J. Nat. Gas Sci. Eng. 2018, 56, 93–106. [Google Scholar] [CrossRef]
  26. Li, Y.; Song, D.; Liu, S.; Ji, X.; Hao, H. Evaluation of pore properties in coal through compressibility correction based on mercury intrusion porosimetry: A practical approach. Fuel 2021, 291, 120130. [Google Scholar] [CrossRef]
  27. Sun, B.; Yang, Q.; Zhu, J.; Shao, T.; Yang, Y.; Hou, C.; Li, G. Pore size distributions and pore multifractal characteristics of medium and low-rank coals. Sci. Rep. 2020, 10, 22353. [Google Scholar] [CrossRef]
  28. Sing, K.S.W.; Everett, D.H.; Haul, R.A.W.; Moscou, L.; Pierotti, R.A.; Rouquérol, J.; Siemieniewaka, T. Reporting physisorption data for gas/solid systems with special reference to the determination of surface area and porosity. Pure Appl. Chem. 1985, 57, 603–619. [Google Scholar] [CrossRef]
  29. Ma, Y.-K.; Nie, B.-S.; He, X.-Q.; Li, X.-C.; Meng, J.-Q.; Song, D.-Z. Mechanism investigation on coal and gas outburst: An overview. Int. J. Miner. Met. Mater. 2020, 27, 872–887. [Google Scholar] [CrossRef]
  30. Zhu, J.; Zhang, R.; Zhang, Y.; He, F. The fractal characteristics of pore size distribution in cement-based materials and its effect on gas permeability. Sci. Rep. 2019, 9, 17191. [Google Scholar] [CrossRef] [PubMed]
  31. Zeng, Q.; Luo, M.; Pang, X.; Li, L.; Li, K. Surface fractal dimension: An indicator to characterize the microstructure of cement-based porous materials. Appl. Surf. Sci. 2013, 282, 302–307. [Google Scholar] [CrossRef]
  32. Wang, J. Multi-Fractal Characteristics of Pore Structure of Coal Rock in Low and Medium Rank CBM Reservoirs. Shiyou Huagong Gaodeng Xuexiao Xuebao/J. Petrochem. Univ. 2019, 32, 26–32. [Google Scholar] [CrossRef]
  33. Mahamud, M.M. Textural characterization of active carbons using fractal analysis. Fuel Process. Technol. 2006, 87, 907–917. [Google Scholar] [CrossRef]
  34. Kravchenko, A.N.; Boast, C.W.; Bullock, D.G. Multifractal Analysis of Soil Spatial Variability. Agron. J. 1999, 91, 1033–1041. [Google Scholar] [CrossRef]
  35. Wang, F.; Li, S. Determination of the Surface Fractal Dimension for Porous Media by Capillary Condensation. Ind. Eng. Chem. Res. 1995, 34, 1383–1386. [Google Scholar] [CrossRef]
  36. Zhang, B.; Liu, W.; Liu, X. Scale-dependent nature of the surface fractal dimension for bi- and multi-disperse porous solids by mercury porosimetry. Appl. Surf. Sci. 2006, 253, 1349–1355. [Google Scholar] [CrossRef]
  37. Zhang, B.; Zhu, J.; Jiang, Y. Multiscale Experiments and Multifractal Analysis of Bituminous Coals: Dominant Pore Properties for Gas Permeability. J. Test. Eval. 2019, 47, 2847–2866. [Google Scholar] [CrossRef]
  38. Cheng, Y.; Luo, X.; Hou, B.; Zhang, S.; Tan, C.; Xiao, H. Pore Structure and Permeability Characterization of Tight Sandstone Reservoirs: From a Multiscale Perspective. Energy Fuels 2023, 37, 9185–9196. [Google Scholar] [CrossRef]
  39. Qin, X.; Cai, J.; Wang, G. Pore-scale modeling of pore structure properties and wettability effect on permeability of low-rank coal. Int. J. Min. Sci. Technol. 2023, 33, 573–584. [Google Scholar] [CrossRef]
  40. Zuowei, H.; Sijing, H.; Donghuan, W.; Yongkun, M.; Xiaoning, L. Application of Multiple Stepwise Regression to Influential Evaluation of Pore-Throat Size on Low-Permeability Sandstone Reservoirs. J. Guilin Univ. Technol. 2013, 33, 21–25. [Google Scholar]
  41. SYT 5336-2006; Practices for Core Analysis. Petroleum Industry Press: Beijing, China, 2006.
  42. Wang, S.; Elsworth, D.; Liu, J. Permeability evolution in fractured coal: The roles of fracture geometry and water-content. Int. J. Coal Geol. 2011, 87, 13–25. [Google Scholar] [CrossRef]
  43. Battistutta, E.; van Hemert, P.; Lutynski, M.; Bruining, H.; Wolf, K.-H. Swelling and sorption experiments on methane, nitrogen and carbon dioxide on dry Selar Cornish coal. Int. J. Coal Geol. 2010, 84, 39–48. [Google Scholar] [CrossRef]
  44. GB/T 21650.1-2008; Pore Size Distribution and Porosity of Solid Materials by Mercury Porosimetry and Gas Adsorption-Part 1: Mercury Porosimetry. China Standard Press: Beijing, China, 2008.
  45. Washburn, E.W. The Dynamics of Capillary Flow. Phys. Rev. B 1921, 17, 273–283. [Google Scholar] [CrossRef]
  46. Farmer, I.; Pooley, F. A hypothesis to explain the occurrence of outbursts in coal, based on a study of west wales outburst coal. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1967, 4, 189–193. [Google Scholar] [CrossRef]
  47. Pan, Z.; Connell, L.D. A theoretical model for gas adsorption-induced coal swelling. Int. J. Coal Geol. 2007, 69, 243–252. [Google Scholar] [CrossRef]
  48. Liu, S.; Harpalani, S. A new theoretical approach to model sorption-induced coal shrinkage or swelling. AAPG Bull. 2013, 97, 1033–1049. [Google Scholar] [CrossRef]
  49. Meng, Y.; Liu, S.; Li, Z. Experimental study on sorption induced strain and permeability evolutions and their implications in the anthracite coalbed methane production. J. Pet. Sci. Eng. 2018, 164, 515–522. [Google Scholar] [CrossRef]
  50. Yu, S.; Bo, J.; Jie-Gang, L. Nanopore Structural Characteristics and Their Impact on Methane Adsorption and Diffusion in Low to Medium Tectonically Deformed Coals: Case Study in the Huaibei Coal Field. Energy Fuels 2017, 31, 6711–6723. [Google Scholar] [CrossRef]
  51. Debelak, K.A.; Schrodt, J. Comparison of pore structure in Kentucky coals by mercury penetration and carbon dioxide adsorption. Fuel 1979, 58, 732–736. [Google Scholar] [CrossRef]
  52. Li, Y.; Lu, G.; Rudolph, V. Compressibility and Fractal Dimension of Fine Coal Particles in Relation to Pore Structure Characterization Using Mercury Porosimetry. Part. Part. Syst. Charact. 1999, 16, 25–31. [Google Scholar] [CrossRef]
Figure 1. Permeation and adsorption simultaneous test system. A: monitor and acquisition data system; B: axial–radial strain transducer; C: gas collection device; D: core gripper; E: axial and confining pressure controllers; F: reference cylinder; G: radial stress transducer; H: axial stress transducer; I: thermostatic controller; J: vacuum pump; K: gas cylinder; 1: double-way valve; 2: double-way valve; 3: double-way valve; 4: double-way valve; 5: pressure relief valve.
Figure 1. Permeation and adsorption simultaneous test system. A: monitor and acquisition data system; B: axial–radial strain transducer; C: gas collection device; D: core gripper; E: axial and confining pressure controllers; F: reference cylinder; G: radial stress transducer; H: axial stress transducer; I: thermostatic controller; J: vacuum pump; K: gas cylinder; 1: double-way valve; 2: double-way valve; 3: double-way valve; 4: double-way valve; 5: pressure relief valve.
Energies 17 00066 g001
Figure 2. Adsorption strain as a function of gas pressure.
Figure 2. Adsorption strain as a function of gas pressure.
Energies 17 00066 g002
Figure 3. Coal permeability curves with helium and nitrogen gases.
Figure 3. Coal permeability curves with helium and nitrogen gases.
Energies 17 00066 g003
Figure 4. The pore volume and pore size distributions by MIP.
Figure 4. The pore volume and pore size distributions by MIP.
Energies 17 00066 g004aEnergies 17 00066 g004b
Figure 5. The relationship between k/k0 and adsorption volumetric strain εv of nitrogen.
Figure 5. The relationship between k/k0 and adsorption volumetric strain εv of nitrogen.
Energies 17 00066 g005
Figure 6. The relationship between adsorption strain and adsorption capacity.
Figure 6. The relationship between adsorption strain and adsorption capacity.
Energies 17 00066 g006
Figure 7. The pore size distributions of coal samples.
Figure 7. The pore size distributions of coal samples.
Energies 17 00066 g007
Figure 8. Surface fractal dimensions of coal samples.
Figure 8. Surface fractal dimensions of coal samples.
Energies 17 00066 g008aEnergies 17 00066 g008b
Table 1. Permeability experiments for different stress conditions.
Table 1. Permeability experiments for different stress conditions.
Stress ConditionValue
Gas inlet pressure (MPa)0.30.511.522.53
Effective stress (MPa)6.86.76.456.25.955.75.45
Table 2. Testing adsorption strain data at different pressure points.
Table 2. Testing adsorption strain data at different pressure points.
Pressure/MPaNo. 1 (10−3)No. 2 (10−3)No. 3 (10−3)
0.30.1420.1120.251
0.50.2310.1730.391
1.00.3920.3080.727
1.50.5430.4240.954
2.00.6310.5991.143
2.50.7040.7981.258
3.00.7570.9891.483
Table 3. Langmuir fitting parameters.
Table 3. Langmuir fitting parameters.
Sample NumberNo.1No.2No.3
εL (10−3)2.5162.5244.513
pL (MPa)6.06.36.0
Table 4. Pore parameters from MIP test.
Table 4. Pore parameters from MIP test.
Sample No.Porosity/%Median Pore
Diameter in Volume/nm
Average Pore
Size/nm
Total Pore
Volume/cm3/g
The Efficiency of Mercury Withdrawal/%
No. 13.96311.49.70.0374.7
No. 25.586127.513.80.033663.1
No. 35.45152.612.90.04254.8
Table 5. The parameters of compressibility of coal samples.
Table 5. The parameters of compressibility of coal samples.
Sample No.Volume of Coal Matrix Skeleton (cm3/g)Bulk Volume of Coal Sample (cm3/g)km
(×10−11 m2/N)
Pore Volume before Modified
(cm3/g)
Pore Volume after Modified
(cm3/g)
Volume Correction Rate
(%)
No.10.72690.75696.96160.03000.013654.557
No.20.71350.75577.19950.03360.014158.124
No.30.72830.77027.43720.04200.021548.868
Table 6. The results of pore fractal dimensions and other parameters.
Table 6. The results of pore fractal dimensions and other parameters.
Sample No.Ds1Ds2Ds3Va
(mL/g)
Vs
(mL/g)
Permeability N2 * (10−18 m2)
No. 12.992.762.560.00500.00870.32
No. 22.892.532.550.00200.01215.73
No. 32.892.532.550.00180.019713.76
Note: Va represents adsorption pore volume of r < 50 nm; Vs represents seepage pore volume of r > 50 nm; * the permeability measured in 3.0 MPa gas pressure.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhu, J.; Shao, T.; Lan, T.; Cheng, Z.; Zhang, Y.; Wang, Q.; Lin, L. Characteristics of Permeability Evolution and Pore Structure of Coal with High Gas. Energies 2024, 17, 66. https://0-doi-org.brum.beds.ac.uk/10.3390/en17010066

AMA Style

Zhu J, Shao T, Lan T, Cheng Z, Zhang Y, Wang Q, Lin L. Characteristics of Permeability Evolution and Pore Structure of Coal with High Gas. Energies. 2024; 17(1):66. https://0-doi-org.brum.beds.ac.uk/10.3390/en17010066

Chicago/Turabian Style

Zhu, Jie, Tangsha Shao, Tianxiang Lan, Zhiyuan Cheng, Yubo Zhang, Quanqi Wang, and Li Lin. 2024. "Characteristics of Permeability Evolution and Pore Structure of Coal with High Gas" Energies 17, no. 1: 66. https://0-doi-org.brum.beds.ac.uk/10.3390/en17010066

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop