Abstract
1. Introduction
Coalbed methane is a serious safety hazard in coal mines [1, 2], while it is also an important economical natural gas supply [3, 4]. The mechanism of gas occurrence, emission, and outburst is not only controlled by macroscopic geological factors [5, 6] but also closely correlated to the microscopic structure in coal reservoirs [7, 8]. The effects of geological factors on gas occurrence and outburst have been studied extensively [9–11]. However, it is still necessary to further study geological significance from the perspective of porosity and adsorptivity in coal reservoirs. The solid-gas coupling between coal and methane is essentially the interaction between gas and inner pores [12], and pore structure directly influences the occurrence of adsorbed gas and free gas in coal seams [13–19]. Therefore, the investigation into pore structure and methane adsorption capacity of coal is the key to analyzing the microscopic mechanism of coalbed gas occurrence and outburst, which is critical to safety mining and coalbed methane resource evaluation.
Coal is a type of porous medium, in which the well-developed pore network provides the storage place and seepage channel for gas [20–24]. The decimal classification is widely adopted in coal pore research [25], namely, micropores (<10 nm), transition pores (10–100 nm), mesopores (100–1000 nm), and macropores (>1000 nm). Generally, pore structure in coal reservoirs can be characterized by image analysis [18, 26], fluid intrusion [27, 28], gas adsorption [13, 15, 29], small-angle scattering of X-rays (SAXS), and neutrons (SANS) [30, 31]. Among them, high-pressure mercury intrusion porosimetry (HP-MIP) has been commonly used in pore characterization for obtaining quantitative information in a wide pore size range (3 nm–100
Due to the complexity of the pore network in coal, the heterogeneity of pore structure is difficult to be described in light of the traditional Euclidean geometry theory, while with the application of modern fractal theory, fractal dimension can quantitatively characterize the complexity of pore structure [16, 34]. Adopting Frenkel-Halsey-Hill (FHH) and Menger models to analyze the data from LP-N2GA and HP-MIP, Yao et al. [15, 28] described the fractal characteristics of nanoscale pores (<100 nm) and micron-scale pores (>100 nm) for different-rank coal, respectively. The fractal heterogeneity of coal pores leads to differences in SSA and gas adsorption capacity [35]; thus, it is necessary to perform a comparative study on fractal pore characteristics between outburst and nonoutburst coals.
Methane is mainly stored in coal as patterns of adsorbed gas, free gas, and dissolved gas [36]. Adsorbed gas accounts for 80%–90% of the total gas content in coal reservoirs [14, 37]. The methane adsorptivity of coal is vital to gas content and gas pressure. Usually, adsorption isotherms are obtained from methane adsorption experiments in the laboratory [20, 38], and adsorption parameters are calculated by fitting an appropriate adsorption model to quantify the adsorption capacities of coal samples. The gas adsorption capacities for coal reservoirs are not only affected by external conditions such as pressure, temperature, and moisture content [39, 40] but also closely related to coal properties including coal rank, maceral composition, pore parameters, and fractal dimension [16, 41, 42]. Tectonically deformed coal (TDC) is prevalent in China, which is usually associated with coal and gas outburst [12, 43]. The coal strength, pore structure, and fractal heterogeneity all changed significantly from the native structure (primary coal) to TDC [12, 43, 44]. However, the pore structure and fractal heterogeneity caused by tectonic deformation and their influence on gas adsorption and occurrence still need to be further studied.
This study focuses on the fractal pore and gas adsorption capacity of the outburst coal and their geological significance to coalbed gas occurrence and outburst. HP-MIP and LP-N2GA were performed to jointly characterize the full-scale pore structure in coal samples. According to the branches of mercury intrusion and nitrogen adsorption, the fractal dimensions were calculated to study the fractal heterogeneity of microscopic pore structure. Isothermal adsorption experiments were performed to determine methane adsorption capacities of coal samples and discover the effects of pore structure and fractal dimension on gas adsorptivity. Moreover, we compared the difference of fractal pores and adsorbability between the outburst coal and nonoutburst coal and revealed their geological significance to gas occurrence and outburst. This study helps to understand coalbed gas geological occurrence and prevent coal and gas outburst.
2. Samples and Methodology
2.1. Geological Background and Sampling
The study area is located in the Qianjiaying coal mine, Kaiping coalfield, China (Figure 1(a)). Main coal seams are distributed in the Upper Carboniferous to Lower Permian Damiaozhuang Formation, and the study area is structurally located in the southeastern Kaiping syncline (Figure 1(b)). The burial depth of coal seam increases from southeast to northwest. The content and pressure of coalbed gas are generally consistent with the variation trend of burial depth, as well as gas occurrence closely controlled by folds and faults. There are two high-risk areas of coal and gas outburst delineated based on gas content, gas pressure, structural complexity, and variation coefficient of coal thickness, comprehensively. The study area had experienced a gas anomaly and outburst, with a volume of 7380 m3 gas emission and 489 t of pulverized coal and rock block ejection at the main crosscut on the -850 m level (Figure 1(c)), which provides a good comparative object for the present study.

Location of the study area (a, b), structural outline map of the coal mine (c), and macroscopic and microscopic characteristics of the coal samples (d).
Eight coal samples were collected from the gas anomaly and outburst area and nonoutburst area (Table 1). The coal samples from gas anomaly and outburst areas mainly belong to scaly coal and wrinkle coal because of the tectonic deformation and coal seam rheology, whereas the samples from nonoutburst areas are mainly primary coal and cataclastic coal (Figure 1(d)). Compared with the primary coal, TDC in the outburst area is characterized by low strength, low firmness, and poor stability through the observation on sample specimen. According to microscopic observation, there are rare microfractures and a few crossed microfractures in primary coal and cataclastic coal, respectively. However, abundant and various microfractures formed a complex network connecting pores in the scaly coal and wrinkle coal from the gas outburst area (Figure 1(d)). The coal samples belong to coking coal with measured vitrinite reflectance (
Basic information and testing data of the coal samples.
Notes: n/a: no available data;
2.2. Experimental Methods
2.2.1. High-Pressure Mercury Intrusion Porosimetry
HP-MIP experiments were conducted using an AutoPore 9510 mercury porosimeter (Micrometrics Instrument, USA) at the Key Laboratory of Coalbed Methane Resources and Reservoir Formation Process, China University of Mining and Technology. According to the China National Standard GB/T 21650.1-2008, coal samples were broken down to about
2.2.2. Low-Pressure N2 Gas Adsorption
LP-N2GA measurements were carried out at the Key Laboratory of Coal Preparation & Purification, China University of Mining and Technology, using a specific surface and pore size analyzer (Autosorb IQ, Quantachrome, USA), which can detect the pore size ranging within 0.35–200 nm. The resolution of the microporous test can reach 0.02 nm, and the minimum detectable SSA and PV are 0.0005 m2/g and 0.0001 cm3/g, respectively. Based on the China National Standard GB/T19587-2004, the coal samples were crushed into 0.18–0.25 mm (60–80 mesh). Then, 5–8 g coal powder samples were dried at 110°C for 12 h in an oven to remove impurity fluids including water vapor. After sample drying, all coal powders were evacuated for complete outgassing in a high vacuum system within 5 h. Using N2 gas with a purity of 99.99%, adsorption and desorption isotherms of coal powders were measured successively with the relative pressures (
2.2.3. Methane Adsorption Experiment
Isothermal adsorption experiments were carried out adopting an IS-300 volumetric adsorption and desorption analyzer (TerraTek, USA) at the Key Laboratory of Coalbed Methane Resources and Reservoir Formation Process, China University of Mining and Technology. Firstly, the coal samples were crushed into 0.18–0.25 mm (60–80 mesh). Secondly, 120–150 g pulverized coal powder was prepared as the equilibrium water sample at the relative humidity of 98%. Lastly, using the pure methane gas (purity: 99.99%), adsorption measurements were performed with pressures ranging within 0–8 MPa at a constant temperature (30°C). The pressure data of both the sample cell and reference cell was recorded to calculate the adsorbed methane content of the coal samples. The Langmuir model (Equation (1)) [51] was chosen to fit the adsorption data and further determine the relevant parameters and isotherms for the coal samples:
2.3. Fractal Dimension Analysis
Based on the fractal theory, the heterogeneity and complexity of pore structure can be quantitatively evaluated with fractal dimension. In this study, the fractal dimensions with different pore size ranges were determined based on the HP-MIP and LP-N2GA data, denoted as
As long as the fractal characteristics exist in the pore structure,
The Frenkel-Halsey-Hill (FHH) model was widely adopted for fractal dimension determination, which is shown in the following equation:
If the pore structure is characterized by fractal heterogeneity,
3. Results
3.1. Pore Characterization Using HP-MIP
The mercury intrusion and extrusion curves are illustrated in Figure 2 for all the coal samples. For the nonoutburst coal samples, the amount of mercury intrusion increased slowly at the low pressures but increased rapidly while pressure is over 5000 psi, indicating that micropores are relatively rich compared to larger-size pores. Both mercury intrusion branch and extrusion branch are concave-up curves, and hysteresis loops are narrow for nonoutburst coal samples. It is speculated that semiclosed pores are dominant in the nonoutburst coal samples, resulting in poor openness and connectivity. For the outburst coal samples, the amount of mercury intrusion increases consistently with pressure from the initial stage to the high-pressure stage, indicating that the PSDs are relatively uniform in the samples. With the decrease of pressure, the mercury extrusion curves initially decreased rapidly, then tended to be flat (sample TC-3 and TC-4), and finally presented a sharp decline with the upper convex shape (sample TC-1 and TC-2). The hysteresis loops are relatively wide for outburst coal samples as a result of the considerable difference in the volume between mercury intrusion and extrusion. Hence, the outburst coal mainly contains open pores, but the inkbottle pores with narrow throats occupy a substantial proportion, which leads to poor pore connectivity.

Mercury intrusion and extrusion curves of the coal samples.
The parameters of pore structure obtained from HP-MIP are shown in Table 2. The measured porosity of coal samples varies between 2.24% and 7.29%. The PV of macropores and mesopores range within 0.0012–0.0174 cm3/g and 0.0006–0.0170 cm3/g, respectively. And the SSAs of macropores and mesopores are 0.001–0.015 m2/g and 0.012–0.249 m2/g, respectively. The relationship between incremental intrusion volume and pore size is plotted in Figure 3. The PSDs are unimodal for the nonoutburst coal samples (Figure 3(a)), and incremental intrusion volume possesses a peak at the range of micro- and transition pores. However, the PSDs are relatively uniform for the outburst coal samples (Figure 3(b)), showing a multimodal distribution, and two major peaks appear in micropores and mesopores. It should be noted that HP-MIP has advantages in macro- and mesopore characterization, though the accuracy is questionable in characterizing pores with a diameter below 50 nm. So the PSDs of pores (<50 nm) should be further analyzed in combination with LP-N2GA.
Pore structure parameters of the coal samples measured by HP-MIP and LP-N2GA.

Relationship between incremental intrusion volume and pore size for nonoutburst coal samples (a) and outburst coal samples (b).
3.2. Pore Characterization Using LP-N2GA
LP-N2GA results show that N2 adsorption curves are in reverse “S” shape for all the coal samples (Figure 4). The shape of adsorption curves depends on the interaction between nitrogen and pore surface in different pressure ranges [54]. When

N2 adsorption and desorption isotherms of the coal samples.
According to the results measured by LP-N2GA (Table 2), the PV of micro- and transition pores in the coal samples are 0.0001–0.0022 cm3/g and 0.0005–0.0036 cm3/g, respectively. The SSAs of micro- and transition pores are 0.052–2.706 m2/g and 0.074–0.602 m2/g, respectively. LP-N2GA has an advantage in the characterization of PSDs of the micro- to transition pores. According to the PSDs calculated using the DFT model presented in Figure 5, the PSDs of the nonoutburst coal samples show multimodal with the peaks in the pore size ranges of 18–20 nm, 30–50 nm, and 70–80 nm (Figure 5(a)). Although the multimodal feature of the PSDs is apparent for the outburst coal samples (Figure 5(b)), the peaks are distributed over a relatively wide pore size range (3–80 nm).

Relationship between incremental adsorbed volume and pore size for nonoutburst coal samples (a) and outburst coal samples (b).
3.3. Pore Fractal Heterogeneity of the Coal Samples
Based on the data from mercury intrusion branches of the coal samples, the correlations between

Plots of
Calculated results of fractal dimensions based on HP-MIP and LP-N2GA.
Notes:
To overcome the limit of fractal dimension calculation in pores ranging within 3–50 nm from HP-MIP data, the fractal dimension

Plots of
3.4. Adsorption Isotherms of the Coal Samples
As illustrated in Figure 8 and Table 4, the methane adsorption capacities of all the coal samples were determined through methane adsorption experiments on the equilibrium water samples at 30°C. The adsorption isotherms were well fitted adopting the Langmuir model (Equation (1)), with

Methane adsorption isotherms of the coal samples.
Fitting results of adsorption parameters for the coal samples.
4. Discussion
4.1. Effect of Pore Structure on Methane Adsorption Capacity
The adsorption capacity of coal is an important factor affecting the gas content in the coal mine. As the place where gas interacts with coal matrix, the pore properties are the most direct factor affecting the methane adsorption in coal reservoir. Also, coal rank, maceral, and other factors can indirectly affect the adsorption performance of coal by controlling the pore structure [13, 41]. According to the scatter diagrams of total PV and SSA versus the maximum methane adsorption capacities of the coal samples (Figure 9), there is a strong positive correlation between

Correlations between the maximum methane adsorption capacity and total pore volume (a) and total specific surface area (b).

Correlations between the maximum methane adsorption capacity and pore volume of micro- and transition pores (a) and meso- and macropores (b).
The pores with different sizes differ in the contribution to the PV and SSA. As shown in Figure 11(a), the total PV of the coal samples is mainly provided by macropores and mesopores. Macropores occupy the largest percentage of the total PV, i.e., 41.4%–59.7% (average: 48.3%), followed by mesopores with the percentage of 19.9%–47.2% (average: 31.3%). However, the total SSA is mainly contributed by micropores and transition pores (Figure 11(b)). The micropores dominated the total SSA with the percentage of 37.3%–82.0% (average: 67.2%), followed by transition pores with the percentage of 13.5%–53.3% (average: 26.6%). On the one hand, given the same total PV, the micro- and transition pores can provide a larger surface area than meso- and macropores for methane adsorption. On the other hand, the distance is relatively short between the pore walls in micro- and transition pores. It may cause the van der Waals force field to overlap between the pore walls, resulting in the larger adsorption potential in nanopores [16, 55]. Therefore, micropores and transition pores are the governing factors on gas adsorption capacity in coal seams. The well-developed micro- and transition pores in the outburst coal are accompanied by a rich internal surface, which enhances the gas adsorption capacity and expansion potentiality.

Pore volume (a) and specific surface area (b) distributions for different pore sizes in the coal samples.
4.2. Effect of Fractal Dimension on Methane Adsorption Capacity
The complexity of the pore structure is another key factor affecting the gas adsorption capacity [13, 35, 42, 56]. Even if coal samples have the same PV or SSA, the pore complexity may be completely different, which is expressed by fractal dimensions

Relationship between the maximum methane adsorption capacity and fractal dimensions.
We further analyze the effect of the fractal dimension on the methane adsorption capacities by plotting the

Relationship between total specific surface areas and fractal dimension
4.3. Fractal Pore Structure of Outburst Coal and Its Geological Significance
4.3.1. Comparison of Fractal Pore Characteristics between Outburst and Nonoutburst Coals
The pore structure and fractal heterogeneity of the coal reservoir are closely related to the gas occurrence since they affect the methane adsorption capacities of the coal samples. The differences of pore structure (including porosity, PV, SSA, PSD, and pore shape) and fractal characteristics between the outburst coal and nonoutburst coal were compared, and their geological significance was analyzed. The porosity, PV, and SSA of the outburst coal samples vary from 4.69% to 7.29%, 0.0188–0.0360 cm3/g, and 2.388–3.294 m2/g, respectively (Table 2). The porosity, PV, and SSA of the nonoutburst coal samples are 2.24%–3.98%, 0.0024–0.0078 cm3/g, and 0.139–0.640 m2/g, respectively. Compared to the nonoutburst coal samples, the porosity, PV, and SSA of outburst coal samples are significantly higher (Figure 11). Quantitatively, the porosity, PV, and SSA of the outburst coal are 1.82 times, 4.56 times, and 5.77 times those of the nonoutburst coal on average, respectively. The well-developed porosity of the outburst coal has two aspects of geological significance to coalbed gas occurrence and outburst. On the one hand, the larger PV leads to the thinner coal skeleton in unit volume, and thus, the coal body is characterized by low strength and fragile structure. On the other hand, the rich surface area inside the outburst coal provides sufficient adsorption sites for methane molecules, and the large pore volume provides sufficient space for gas storage. Therefore, gas is easy to accumulate in the outburst coal seam, which increases gas pressure and outburst tendency.
The pore shape of coal can be analyzed based on the hysteresis loops determined by the mercury intrusion/extrusion and the N2 adsorption/desorption [35, 54]. Yao and Liu [57] concluded that the open pore and inkbottle pore can retain mercury and lead to a hysteresis loop, while the semiclosed pore cannot. Taking samples ZC-S8 and TC-4 as an example, the difference of pore shape between the outburst coal and nonoutburst coal can be analyzed. As shown in Figure 14(a), the hysteresis loops of the nonoutburst coal are narrow or inconspicuous. Accordingly, the nonoutburst coal is dominated by semiclosed pores, including narrow slit-like pores, wedge-shaped pores, and conical/tubular pores with a closed end. This type of pore system is of poor openness and connectivity, which is not conducive to gas occurrence and enrichment. However, the hysteresis loops of the outburst coal are wide, and the N2 desorption branch presents a sharp decline (Figure 14(b)). Therefore, the pores in the outburst coal are mainly open pores, including parallel plate-like pores open on all sides and cylindrical pores open at both ends. In addition, there are a certain number of inkbottle pores with narrow throats and wide bodies, which limit the pore connectivity. This type of pore structure is beneficial to gas enrichment and storage, and abundant gas can be stored in the wide bottle body of the inkbottle pores. More importantly, the narrow pore throat restricts the gas seepage. Thus, the gas does not migrate but accumulates locally to form a high-pressure gas zone.

Schematic diagram of pore shape for the nonoutburst coal (a) and outburst coal (b) based on the HP-MIP and LP-N2GA data.
PSD is directly related to the occurrence and migration of coalbed gas. Micro- and transition pores provide the primary space for adsorbed gas, while meso- and macropores mainly constitute the diffusion and seepage space for free gas [15, 28]. LP-N2GA and HP-MIP were combined to characterize PSD in pore sizes ranging from 1 to 50 nm and >50 nm, respectively, according to their advantages. As the full-scale PSDs are shown in Figure 15, the nonoutburst coal is dominated by micropores. However, the pore size is widely distributed in the outburst coal, in which not only micropores are dominant, but also, transition pores and mesopores are developed to a certain extent. The difference of PSDs between outburst and nonoutburst coals is mainly affected by tectonic deformation. As the degree of tectonic deformation increases (sequence from primary coal, cataclastic coal, scaly coal, to wrinkle coal), the micropores in TDC increase significantly because of the squeezing action [12], which corresponds to the narrow throat of the inkbottle pores in the outburst coal. In addition, new macropores and microfractures were formed in TDC under the tectonic deformation, which provides a favorable channel for rapid and instantaneous gas release in gas outburst.

The pore size distribution of the coal samples from HP-MIP and LP-N2GA data.
According to the comparison of the fractal dimensions between the nonoutburst coal and outburst coal (Figure 16),

Fractal dimensions of the coal samples with different tectonically deformed types.
4.3.2. Geological Implication of Fractal Pore Structure on Gas Occurrence and Outburst
The coal and gas outburst occurred in the zone of TDC including scaly coal and wrinkle coal in the study area. Meanwhile, there are significant differences in pore structure between TDC and primary coal as mentioned in Section 4.3.1. To reveal the relationship between tectonic deformation and gas outburst, the combined effect of tectonic deformation, pore structure, and adsorption capacity of coal reservoir on gas occurrence and outburst was summarized. Under the tectonic deformation, the primary coal transits to TDC accompanied by the direct effect on the physical properties and pore structure (Figure 17). Firstly, the coal strength is reduced yielding low hardness and brittle structure, thus forming the zone of soft and broken coal. Secondly, TDC is characterized by abundant PV and SSA and developed microporosity and strong fractal heterogeneity, resulting in a large number of internal surfaces for methane adsorption. Accordingly, the gas adsorption capacities of TDCs are remarkably enhanced. Thirdly, inkbottle pores with narrow throats emerge under the tectonic deformation, resulting in poor pore connectivity, which is unfavorable for gas migration and extraction. In addition, the TDC area is also a geostress concentration zone with high stress and high energy. Cheng and Lei [58] also suggested that in situ coal reservoirs in the TDC area possess low permeability as the fracture system was severely squeezed. Therefore, TDC is of strong adsorption capacity but poor permeability, and coalbed gas is enriched to form a high-pressure gas sealing zone.

The combined effect of tectonic deformation, pore structure, and gas adsorptivity on coalbed gas occurrence and outburst.
In fact, coal and gas outburst can be regarded as a combined action of gas, the porous medium of coal/rock, and underground stress according to the integrated factor outburst hypothesis [11]. Li [59] proposed that, for coal and gas outburst, high-pressure gas is an essential factor, while coal structure is an obstacle, and underground stress acts as a dynamic factor. The place, where TDC is well-developed, is in an unstable state due to the combination of soft and deformed coal zone, high-pressure gas sealing zone, and underground stress concentration zone (Figure 17). Once the triggering factors (excavating disturbance, rock burst, active geological structure, temperature anomaly, etc.) break the original equilibrium state, the adsorbed gas is desorbed to free gas rapidly, resulting in a gas expansion in the coal reservoir. Ultimately, high-pressure free gas and crushed coal powder/block are released instantaneously driven by stress energy, which leads to abnormal gas emission, even continuous coal fragmentation, and gas outburst. In summary, weak coal structure, abundant gas enrichment, and concentrating underground stress are the prerequisites for coal and gas outburst. The risk of gas outburst in the zone of TDC is higher than that in the zone of primary coal.
5. Conclusions
In this work, qualitative description, quantitative experiments, and fractal dimension calculation have been carried out to compare the pore structure and methane adsorption capacity between outburst and nonoutburst coals. Further, the effects of fractal pore on gas adsorption capacity and its geological significance to coalbed gas occurrence and outburst were discussed. The following conclusions can be drawn from the present study.
On average, the PV and SSA of outburst coal samples are 4.56 and 5.77 times those of nonoutburst coal samples, respectively. The pore size is widely distributed in the outburst coal compared with the nonoutburst coal. Developed porosity in the outburst coal is beneficial to gas enrichment but leads to weak coal strength. Meanwhile, the considerable inkbottle pores in the outburst coal restrict gas migration and extraction, which is an important microscopic factor leading to gas outburst Pore structure in coal reservoir exhibits fractal characteristics. The fractal dimensions The methane adsorption capacities of the coal samples are governed by micropores and transition pores and affected by the fractal dimension The outburst coal is mostly scaly coal and wrinkle coal with tectonic deformation in the study area. The TDC is characterized by weak coal strength, complex pore structure, and strong adsorption capacity but limited permeability. Therefore, the area of TDC is at high risk for gas outburst, where there is a combination of soft coal zone, high-pressure gas sealing zone, and underground stress concentration zone
