Damage Effects and Fractal Characteristics of Coal Pore Structure during Liquid CO 2 Injection into a Coal Bed for E-CBM

Pore structure has a significant influence on coal-bed methane (CBM) enhancement. Injecting liquid CO2 into coal seams is an effective way to increase CBM recovery. However, there has been insufficient research regarding the damage effects and fractal characteristics of pore structure at low temperature induced by injecting liquid CO2 into coal samples. Therefore, the methods of low-pressure nitrogen adsorption-desorption (LP-N2-Ad) and mercury intrusion porosimetry (MIP) were used to investigate the damage effects and fractal characteristics of pore structure with full aperture as the specimens were frozen by liquid CO2. The adsorption isotherms revealed that the tested coal samples belonged to type B, indicating that they contained many bottle and narrow-slit shaped pores. The average pore diameter (APD; average growth rate of 18.20%), specific surface area (SSA; average growth rate of 7.38%), and total pore volume (TPV; average growth rate of 18.26%) increased after the specimens were infiltrated by liquid CO2, which indicated the generation of new pores and the transformation of original pores. Fractal dimensions D1 (average of 2.58) and D2 (average of 2.90) of treated coal samples were both larger the raw coal (D1, average of 2.55 and D2, average of 2.87), which indicated that the treated specimens had more rough pore surfaces and complex internal pore structures than the raw coal samples. The seepage capacity was increased because D4 (average of 2.91) of the treated specimens was also higher than the raw specimens (D4, average of 2.86). The grey relational coefficient between the fractal dimension and pore structure parameters demonstrated that the SSA, APD, and porosity positively influenced the fractal features of the coal samples, whereas the TPV and permeability exerted negative influences.


Introduction
Coal is China's main energy source and accounted for 67.4% of the total energy consumption in the latest survey [1,2].In recent decades, the relationship between resource utilization and environmental protection has become increasingly tense [3,4].Coal resource development has generated massive amounts of CO 2 , CH 4 , and other pollutant gases, which threaten the safety of workers and have exacerbated the greenhouse effect [5,6].Coal-bed methane (CBM), a potential natural gas resource, is abundant in Chinese coal seams.Pre-extracting gas before coal mining can control gas emissions and has a positive influence on economic development and the environment [7].
Coal, an organic, fractured solid material with complex pores, supplies numerous gas storage sites and affects coal-bed methane extraction [8].However, owing to the low permeability, high gas content, deep mine pressure behaviors, and other issues, gas drainage productivity is lower in most Chinese mines.Therefore, improving the permeability of coal seams is essential for enhancing the efficiency of gas extraction [9].Injecting liquid CO 2 (LqCO 2 ) into the coal seam to enhance methane extraction is a common technology that has been used internationally [10,11].
Generally, coal pore structure is classified into five types according to IUPAC and the BB Hodot theory [12,13]: micropores (2-10 nm in diameter), transition pores (10-100 nm in diameter), mesopores (100-1000 nm in diameter), macropores (1000-10,000 nm in diameter), and microfractures (>10,000 nm in diameter).This classification method was adopted in this research.Several studies have demonstrated that micropores and transition pores are the main adsorption and storage sites of CBM, whereas mesopores and macropores are the main channels for CBM seepage and diffusion [14,15].Coal pore structure parameters include shape, size distribution, specific surface area, and pore volume.Currently, the typical methods of studying pore structure are centered on scanning electron microscopy (SEM), low-pressure nitrogen adsorption (LP-N 2 -Ad), mercury intrusion porosimetry (MIP), nuclear magnetic resonance (NMR), and small-angle X-ray scattering (SAXS) [16][17][18][19].Therefore, in this work, a combination of LP-N 2 -Ad (the measuring aperture rang is 0-300 nm) and MIP (the measuring aperture rang is 0.3-360 µm) were used to study the damage effect and fractal characteristics of pore structures during low-temperature freezing-thawing with LqCO 2 .
The relationship between pore structure and CBM properties has been studied by several researchers in recent decades [20].Xu [21] and Vishal [22] et al. adopted several methods to research the variation and permeability of coal pore structures in multiple LqCO 2 freezing-thawing cycles.Their results indicated that the multiple freezing-thawing cycles positively affect the porosity and permeability of coal pores.Zhang et al. [23,24] has been researched on the supercritical CO 2 injected into coal and concluded that the swelling stress in the coal matrix can fracture the inorganic mineral and discovered the supercritical CO 2 could closed the cleats and made some new fractures generated in the un-swelling phase by the swelling effect.Wen et al. [25] conducted an in situ fracturing test in the Yanzhou Coal Field by injecting LqCO 2 and reported the mechanism of damage during liquid CO 2 injection.Ma et al. [26] performed a liquid CO 2 injection test in the Hancheng Mining Area and observed an improved gas extraction efficiency.However, little research has focused on the structural characteristics and fractal geometry features of pores under low temperatures.It is difficult to accurately describe the heterogeneity of coal reservoirs using traditional geometric methods as coal pore structure is a complex, three-dimensional geometric model [27].A theory based on fractal geometry was presented by Mandelbrot in 1975 [28], and recent studies have stated that fractal theory has provided a useful method to quantitatively define the complexity of the structure, surface roughness, and heterogeneity of coal pores.
Therefore, this study will examine the variation in pore structure parameters before and after coal samples are treated by LqCO 2 , and the fractal dimensions D 1 and D 2 will be calculated by LP-N 2 -Ad, and D 4 by MIP.Furthermore, the grey relational theory will be used to describe the relationship between fractal dimension and pore structure parameters.This study could provide a theoretical basis for E-CBM by liquid CO 2 injection.

Coal Samples Re-Preparation
To discuss the damage effects and fractal characteristics of coal pore structure during LqCO 2 injection, three coal samples were taken directly from previous industrial test mines in China: meager coals from the Hancheng Mining Area, Shaanxi Province; and 1/3 cocking coal (denoted as 1/3CC) from Huainan City, Anhui Province.One of the meager coals has a greater hardness (denoted as MHC), while the other is soft (denoted as MSC).Three coal samples were collected following the Chinese Standard Method GB/T 1922-2003.One kilogram of each type of coal was collected and then sent to the experimental laboratory immediately after classification and packing.The basic properties of the coal samples are presented in Table 1.To test the initial and LqCO 2 -treated coal samples with full-aperture pore structure-evolution tendencies at room temperature, the low-temperature N 2 isotherm adsorption/desorption method was followed using ASAP 2020 (Micromeritics Instrument, Shanghai, China) and the international standard, ISO 15901-3:2007.Auto pore IV 9500 equipment (Shanghai, China) was used for the MIP measurement in accordance with the international standard ISO-15901-1:2005.The coal sample was smashed to a particle size of 80-120 mesh (0.178-0.125 mm) for the LP-N 2 -Ad test and was crushed to approximately 3.0 mm for the MIP test.Prior to the LP-N 2 -Ad and MIP tests, the test coal samples were infiltrated by LqCO 2 , and the infiltration process of three coal samples were as follows.
(1) Three specimens were infiltrated in LqCO 2 at −50 • C for 4 h with self-developed equipment.
A previous study indicated that the cooling radius of LqCO 2 injected into coal seam was approximately 10 m, meanwhile the coal around the injection hole was fully infiltrated after 4 h of LqCO 2 injection [25,26].(2) After the infiltration process, raw and treated coal samples were added into a vacuum oven at 65 • C for 12 h until the mass changes in the initial and treated coal samples did not exceed 0.2% [5].(3) All the coal samples were desorbed in a vacuum chamber for 48 h to remove the mixed gases adsorbed on the surface of the coal samples.

Fractal Dimension from LP-N 2 -Ad Isotherms
Several studies have demonstrated that fractal Brunauer-Emmett-Teller (BET) and Frenkel-Halsey-Hill (FHH) models can be the most practical methods for describing the fractal features of the adsorbed pores of coal reservoirs [29,30].The FHH model is the most accurate and simplest method of calculating fractal dimensions, as follows [31].
where P is the gas equilibrium pressure, MPa; P 0 is the gas saturation pressure, MPa; V is the volume of the gas adsorbed at the equilibrium pressure, cm 3 /g; D is the fractal dimension, which can be calculated by Equation (2), whereas D = 3A + 3 is usually <2, and this calculation method is beyond the definition of fractal dimension; and A is the slope, which can be calculated by plotting the gas adsorption isotherm data in terms of lnV vs. ln[ln(P 0 /P)].

Fractal Dimension from MIP
The fractal dimension of seepage pores (>100 nm in diameter) can be calculated from the data provided by the MIP test.Several mathematical models for the fractal analysis of seepage pores have been proposed and discussed by scholars in detail [32].The seepage pore fractal model was initially presented by Mikula and Friesen [33], and was selected to calculate the fractal dimension of seepage pores in this study (Equation (3)): where P represents the absolute injection pressure, MPa; V is the cumulative mercury injection volume at the correlative pressure P, cm 3 /g; and D is the fractal dimension of the seepage pores that was obtained from the MIP data.

Grey System Theory
The grey relational space theory was proposed in Chinese by professor Deng JL and has been developed rapidly over the last few decades [34,35].Grey relational analysis is a systematic method that can quantitatively compare and describe the subject investigated [35].The core target of grey relational analysis is determining whether the influencing factors and the research objects are closely linked according to the degree of similarity between the shapes of the sequence curve.The degree of similarity between the research object and impact factors is also described using grey relational analysis [36,37].Using this method, the similarity between the coal pore fractal dimensions and structural parameters (including porosity, permeability, APD, SSA, and TPV) before and after LqCO 2 treatment can be analyzed, and the grey relational degree can be used to quantitatively describe the damage to the coal pore structure and changes in its fractal characteristics caused by LqCO 2 freezing-thawing.
The fractal dimensions of coal pore structure for the initial and LqCO 2 -treated coal samples can be taken as the system behavior sequences X 0 and X 0 , respectively.Additionally, the structural coal pore parameters were selected as the systematically internal correlation sequences X i and X i (i = 1, 2, 3, 4, 5).Therefore, the correlation coefficients ζ i (k) and ζ i (k) can be defined as follows [34]: The grey relational degree r i , r i (ranging from 0.0 to 1.0) among X 0 , X 0 and X i , X i can then be calculated by the following equation [34-37]: where k is a constant with a value of 1, 2, or 3; ρ is a correction factor with a value of 0.5; and n represents the three coal samples with a value of 3. The LP-N 2 -Ad isotherms for the studied coal samples are shown in Figure 1.All isotherms for the adsorption curves show that the hysteresis loops for the tested coal samples are of type B according to the International Union of Pure and Applied Chemistry.When the relative pressure is 0.5, the desorption curves have an inflection point that divides the curves into two regions (regions 1 and 2).At a relative pressure of 0 < P/P 0 < 0.5 (region 1), the adsorption curves increased more slowly, and the coal matrix surface is a monolayer of adsorbed nitrogen molecules [38].Furthermore, multilayer adsorption occurred at a relative pressure of 0.5 < P/P 0 < 0.95 (region 2), when the adsorption curves gradually increased.The adsorption curves then increased sharply as the relative pressure exceeded 0.95, indicating that capillary condensation occurred in accordance with the Kelvin equation, increasing the quantity of adsorbed nitrogen [39].

Results and Discussion
The type of pore hysteresis loop indicated the presence of cylindrical-shaped pores with one closed side, slit-shaped pores with all sides open, and bottle-shaped pores in the studied coal samples [40].However, there are many open pores in the 1/3CC coal sample as the hysteresis loop was not closed because the relative pressure was below 0.5 [41].Thus, the pore structure of the studied coal sample is appropriate for CBM gathering, but detrimental for gas seepage [42].Moreover, Figure 1 indicates a large difference between the raw and LqCO 2 -frozen coal samples.The quantities of adsorbed nitrogen molecules in the untreated MHC, MSC, and 1/3CC coal samples were 2.46, 1.14, and 5.56 cm 3 /g (Figure 1a), respectively.However, when the coal samples were infiltrated by LqCO 2 , the quantities of adsorbed nitrogen molecules were 2.93, 1.54, 6.23 (Figure 1b), respectively.This indicates that the structure of coal pores was damaged by the frozen-heave force as the LqCO 2 infiltrated coal samples, causing changes in the relevant parameters.The coal pore structure was also damaged by the prestressing and shearing stress in the infiltration process, performed the coal matrix shrinkage and the surface was more rough and complex under low temperature by LqCO 2 infiltration [25].This result was different from the study of Zhang et al. [23,24], and Zhang's article mainly used Discrete Element Method (DEM) model to simulate the cleats' variation in the swelling process under high temperature and high pressure conditions [24].To accurately determine the effect of LqCO 2 infiltration, further analysis is based on the pore structure parameters of the coal samples [5].
Resources 2018, 7, x FOR PEER REVIEW 5 of 17 [40].However, there are many open pores in the 1/3CC coal sample as the hysteresis loop was not closed because the relative pressure was below 0.5 [41].Thus, the pore structure of the studied coal sample is appropriate for CBM gathering, but detrimental for gas seepage [42].Moreover, Figure 1 indicates a large difference between the raw and LqCO2-frozen coal samples.The quantities of adsorbed nitrogen molecules in the untreated MHC, MSC, and 1/3CC coal samples were 2.46, 1.14, and 5.56 cm 3 /g (Figure 1a), respectively.However, when the coal samples were infiltrated by LqCO2, the quantities of adsorbed nitrogen molecules were 2.93, 1.54, 6.23 (Figure 1b), respectively.This indicates that the structure of coal pores was damaged by the frozen-heave force as the LqCO2 infiltrated coal samples, causing changes in the relevant parameters.The coal pore structure was also damaged by the prestressing and shearing stress in the infiltration process, performed the coal matrix shrinkage and the surface was more rough and complex under low temperature by LqCO2 infiltration [25].This result was different from the study of Zhang et al. [23,24], and Zhang's article mainly used Discrete Element Method (DEM) model to simulate the cleats' variation in the swelling process under high temperature and high pressure conditions [24].To accurately determine the effect of LqCO2 infiltration, further analysis is based on the pore structure parameters of the coal samples [5].  Figure 2 shows the relationships of the APD, BET surface area, and Barrett-Joyner-Halenda (BJH) pore volume with the untreated and LqCO 2 -treated coal samples under relative pressure conditions.The APD, BET surface area, and BJH pore volume increased when the coal samples were frozen by LqCO 2 , which is in agreement with the results of Wen and Xu et al. [5,21].The APD of the adsorbed pores of the treated samples were 9.88, 6.34, and 6.96 nm larger than those of the initial samples, respectively.The SSA and TPV of the three LqCO 2 -treated coal samples were generally higher than those of the original samples.For instance, the SSA of the treated coal samples were 23.53 (MHC), 13.52 (MSC), and 4.08% (1/3CC) higher than that of the raw samples, and the TPV were 14.90 (MHC), 12.64 (MSC), and 11.87% (1/3CC) higher, respectively.It can be concluded that the pore structure of coal is damaged by the frost heaving force of liquid CO 2 at low temperature [41].The damage occurs due to the generation of countless micropores and transition pores, as the SSA and TPV of the adsorbed pores are higher.However, many of the adsorbed pores are converted into mesopores and macropores at low temperature [21].The SSA and TPV of the micropores and transition pores are presented in Table 2.These results indicated that the pore structures of the treated coal samples varied distinctly.
Resources 2018, 7, x FOR PEER REVIEW 6 of 17 Figure 2 shows the relationships of the APD, BET surface area, and Barrett-Joyner-Halenda (BJH) pore volume with the untreated and LqCO2-treated coal samples under relative pressure conditions.The APD, BET surface area, and BJH pore volume increased when the coal samples were frozen by LqCO2, which is in agreement with the results of Wen and Xu et al. [5,21].The APD of the adsorbed pores of the treated samples were 9.88, 6.34, and 6.96 nm larger than those of the initial samples, respectively.The SSA and TPV of the three LqCO2-treated coal samples were generally higher than those of the original samples.For instance, the SSA of the treated coal samples were 23.53 (MHC), 13.52 (MSC), and 4.08% (1/3CC) higher than that of the raw samples, and the TPV were 14.90 (MHC), 12.64 (MSC), and 11.87% (1/3CC) higher, respectively.It can be concluded that the pore structure of coal is damaged by the frost heaving force of liquid CO2 at low temperature [41].The damage occurs due to the generation of countless micropores and transition pores, as the SSA and TPV of the adsorbed pores are higher.However, many of the adsorbed pores are converted into mesopores and macropores at low temperature [21].The SSA and TPV of the micropores and transition pores are presented in Table 2.These results indicated that the pore structures of the treated coal samples varied distinctly.Notes: S BET , BET surface area, m 2 /g; S mi , micropore surface area, m 2 /g; S tra , transition pore surface area, m 2 /g; V BJH , BJH pore volume, cm 3 /g; V mi , micropore volume, cm 3 /g; V tra , transition pore volume, cm 3 /g.

Fractal Characteristics of the Adsorbed Pores from Raw and LqCO 2 -Treated Coal Samples
The FHH model was used to calculate the fractal dimensions of the original and LqCO 2 -treated coal samples, as shown in Figure 3.The ln(ln(P 0 /P)) vs. ln(V) fitting curves represented the LP-N 2 -Ad isotherms and can be divided into two segments with delivery points equal to approximately −1.5.
From the above analysis and isotherm descriptions, monolayer adsorption generally occurred in region 1, where nitrogen molecules are gradually adsorbed to the pores of the coal sample.Accordingly, the fractal dimension D 1 was calculated to be in region 1 and could be used to describe the fractal characteristics of the pores' surface areas [16].Multilayer adsorption and capillary condensation mainly occurred in region 2, and the fractal dimension D 2 was calculated in region 2 and can be used to describe fractal characteristics of internal pore structure [43].

Fractal Characteristics of the Adsorbed Pores from Raw and LqCO2-Treated Coal Samples
The FHH model was used to calculate the fractal dimensions of the original and LqCO2-treated coal samples, as shown in Figure 3.The ln(ln(P0/P)) vs. ln(V) fitting curves represented the LP-N2-Ad isotherms and can be divided into two segments with delivery points equal to approximately −1.5.
From the above analysis and isotherm descriptions, monolayer adsorption generally occurred in region 1, where nitrogen molecules are gradually adsorbed to the pores of the coal sample.Accordingly, the fractal dimension D1 was calculated to be in region 1 and could be used to describe the fractal characteristics of the pores' surface areas [16].Multilayer adsorption and capillary condensation mainly occurred in region 2, and the fractal dimension D2 was calculated in region 2 and can be used to describe fractal characteristics of internal pore structure [43].The fractal dimensions D 1 and D 2 for the initial and treated coal samples are presented in Table 3.The D 1 and D 2 of the three treated coal samples were typically higher than those of the untreated samples.D 1 ranged from 2.46 to 2.63 (avg.2.55), while D 2 ranged from 2.83 to 2.93 (avg.2.87) for the raw coal samples.After the coal samples were frozen treated by LqCO 2 , D 1 ranged from 2.47 to 2.65 (avg.2.58) and D 2 ranged from 2.85 to 2.94 (avg.2.90).It can be known that the coal matrix's surface and internal pore structure of the treated samples were more rough and complex due to the damaged effect by the frost-heave force as the coal samples infiltrated in LqCO 2 , which caused the increase in the fractal dimension [22,26,42].There is a positive linear correlation between D 1 and D 2 (Figure 3), which is consistent with the results of previous research [5,16,43].The fractal characteristics of the adsorbed pores were clearer when the coal samples were frozen by LqCO 2 [32].Research has demonstrated that the mesopores and macropores influence coal pore seepage properties [44], exceeding the test scope of LP-N 2 -Ad.Therefore, other investigations should be conducted.

MIP Curves
The MIP method is widely used for analyzing the structural characteristics of the seepage pores of porous materials [15,31].Thus, this method was applied in this study to test the structural properties of the seepage pores and their fractal characteristics.The MIP curves for the initial and LqCO 2 -treated samples are shown in Figure 4. MIP curves have some clear traits that have been researched previously [5,19], and can be divided into three sections, as shown in Figure 4.The cumulative mercury volume increased rapidly in section 1, before the intrusion pressure increased from 0.0013 to 0.035 MPa, the intruded mercury volume for the raw coal samples were 0.015 (MHC), 0.034 (MSC), and 0.039 mL/g (1/3CC), respectively, while those for the three treated coal samples were 0.027 (MHC), 0.041 (MSC), and 0.044 mL/g (1/3CC), respectively.At this point, the mercury molecules come into contact with the coal's surface and enter the pores (pore size: 10-20 µm).When the pressure exceeded 0.0013 MPa, the curves changed as the intrusion pressure continued to increase, the cumulative mercury volume in Section 2 increased slowly until the pressure reached 50 MPa.During this period, the mercury molecules entered the internal pore spaces (pore size: 0.03-0.15µm), and the mercury intrusion volumes of the untreated coal samples were 0.025, 0.042, and 0.051 mL/g, respectively.However, those of the treated coal samples were 0.051, 0.052, and 0.058 mL/g, respectively.When the pressure reached 50 MPa, the yield point was reached and coal matrix compression began, particularly for the seepage pores (i.e., mesopores and macropores).The capillary condensation phenomenon would occur in Section 3 as the injected mercury pressure increased from 50 to 100 MPa.Furthermore, the coal matrix would be compressed by the largest injection pressure.Therefore, the intrusion volumes in this section were 0.038, 0.063, 0.0701 mL/g for the initial samples, respectively, and 0.051, 0.072, 0.077 mL/g for the infiltrated coal samples, respectively.From the MIP curves, it can be concluded that pore volume is significantly enhanced by treatment with LqCO2.The total pore volume was increased, based on the results of the MIP test.Some new adsorbed and seepage pores may have been generated, and the coal pore structures may have converted between each other [19].Therefore, more detailed investigation of the next step is required.

Pore Structure Parameters Analysis of Seepage Pores
The variation in pore structure parameters obtained from the MIP for the coal samples is listed in Table 4. Before the infiltration of LqCO2, the porosity of the initial samples was 5.47, 7.63, and 8.77%, respectively.However, those of the treated coal samples increased to 7.03, 8.57, and 9.44%, respectively.Thus, the porosity increased by an average of 12.82% after the coal samples were treated by LqCO2.Similarly, the permeability of the treated samples was increased at an average rate of 7.64%.The growth rates of the APD, SSA, and TPV were 18.50, 16.16, and 20.11% (avg.18.2%), 6.34, 4.42, and 11.38% (avg.7.38%), and 8.60, 13.95, and 25.20% (avg.18.26%) higher than those of the raw coal samples, respectively.It is revealed that new pores were generated, including more macropores and some micro-fractures from the transformation of prime pores during LqCO2 treatment [5,27].These results exhibit similar variations in micropores and transition pores to those obtained using LP-N2-Ad.
The pore size distribution and incremental pore volume of the untreated and treated coal samples are given in Figure 5.The incremental seepage pore volume for the treated coal specimens is larger than that of the original coal samples, and the adsorbed pore volume of the treated specimens was slightly higher because the structure of coal pore was damaged by LqCO2 infiltration under the low temperature (i.e., −50 °C), indicating that the methane seepage and diffusion ability of the treated  From the MIP curves, it can be concluded that pore volume is significantly enhanced by treatment with LqCO 2 .The total pore volume was increased, based on the results of the MIP test.Some new adsorbed and seepage pores may have been generated, and the coal pore structures may have converted between each other [19].Therefore, more detailed investigation of the next step is required.

Pore Structure Parameters Analysis of Seepage Pores
The variation in pore structure parameters obtained from the MIP for the coal samples is listed in Table 4. Before the infiltration of LqCO 2 , the porosity of the initial samples was 5.47, 7.63, and 8.77%, respectively.However, those of the treated coal samples increased to 7.03, 8.57, and 9.44%, respectively.Thus, the porosity increased by an average of 12.82% after the coal samples were treated by LqCO 2 .Similarly, the permeability of the treated samples was increased at an average rate of 7.64%.The growth rates of the APD, SSA, and TPV were 18.50, 16.16, and 20.11% (avg.18.2%), 6.34, 4.42, and 11.38% (avg.7.38%), and 8.60, 13.95, and 25.20% (avg.18.26%) higher than those of the raw coal samples, respectively.It is revealed that new pores were generated, including more macropores and some micro-fractures from the transformation of prime pores during LqCO 2 treatment [5,27].These results exhibit similar variations in micropores and transition pores to those obtained using LP-N 2 -Ad.
The pore size distribution and incremental pore volume of the untreated and treated coal samples are given in Figure 5.The incremental seepage pore volume for the treated coal specimens is larger than that of the original coal samples, and the adsorbed pore volume of the treated specimens was slightly higher because the structure of coal pore was damaged by LqCO 2 infiltration under the low temperature (i.e., −50 • C), indicating that the methane seepage and diffusion ability of the treated samples were increased, while the adsorbed capacity was decreased [42].Furthermore, the application site of LqCO 2 injection has revealed that the permeability and gas extraction efficiency of the coal seam was significantly improved [25,26].However, the fracture damage effect of coal with the LqCO 2 injection needs further research, and the result of Zhang's research will provide great enlightenment for future research.samples were increased, while the adsorbed capacity was decreased [42].Furthermore, the application site of LqCO2 injection has revealed that the permeability and gas extraction efficiency of the coal seam was significantly improved [25,26].However, the fracture damage effect of coal with the LqCO2 injection needs further research, and the result of Zhang's research will provide great enlightenment for future research.

Fractal Features of Seepage Pores between the Initial and LqCO2-Treated Coal Samples
The calculation processes for the fractal dimensions of the seepage pores in the tested coal samples are shown in Figure 6.The fractal dimensions D3, D4, and D5 were calculated from the three parts that were divided by the relationships between the ln( / ) dV dP and ln(P) for the three coal samples, based on Equation (3).Previous studies demonstrated that D3, D4, and D5 were related to coal particle shape, internal seepage pore structure, and coal matrix compressibility [15,45], respectively.The values of the three fractal dimensions for all coal samples are listed in Table 5.From the fractal theory model, the fractal dimension Di ranged from 2 to 3 and is a suitable parameter for discussing the fractal characteristics of the coal matrixes surface [16].The value for D3 was below 2 (avg.1.80) except for the untreated 1/3CC coal sample, while all values for D5 exceeded 3 (avg.3.79).D3 and D5 were outside the range of the denoted fractal theory, therefore these dimensions are meaningless in this section.However, fractal dimension D4, which corresponds to intra-particle filling and can describe the mesopore and macropore matrix surface features, was between 2 and 3 for all samples, and thus D4 can be used to determine the characteristics of the seepage pores [45].The value of D4 for the original coal sample ranged from 2.84 to 2.89 (avg.2.86), while that for the treated samples increased from 2.88 to 2.97 (avg.2.91) at an average rate of increase of 1.62%.

Fractal Features of Seepage Pores between the Initial and LqCO 2 -Treated Coal Samples
The calculation processes for the fractal dimensions of the seepage pores in the tested coal samples are shown in Figure 6.The fractal dimensions D 3 , D 4 , and D 5 were calculated from the three parts that were divided by the relationships between the ln(dV/dP) and ln(P) for the three coal samples, based on Equation (3).Previous studies demonstrated that D 3 , D 4 , and D 5 were related to coal particle shape, internal seepage pore structure, and coal matrix compressibility [15,45], respectively.The values of the three fractal dimensions for all coal samples are listed in Table 5.From the fractal theory model, the fractal dimension D i ranged from 2 to 3 and is a suitable parameter for discussing the fractal characteristics of the coal matrixes surface [16].The value for D 3 was below 2 (avg.1.80) except for the untreated 1/3CC coal sample, while all values for D 5 exceeded 3 (avg.3.79).D 3 and D 5 were outside the range of the denoted fractal theory, therefore these dimensions are meaningless in this section.However, fractal dimension D 4 , which corresponds to intra-particle filling and can describe the mesopore and macropore matrix surface features, was between 2 and 3 for all samples, and thus D 4 can be used to determine the characteristics of the seepage pores [45].The value of D 4 for the original coal sample ranged from 2.84 to 2.89 (avg.2.86), while that for the treated samples increased from 2.88 to 2.97 (avg.2.91) at an average rate of increase of 1.62%.

Fractal Features of Seepage Pores between the Initial and LqCO2-Treated Coal Samples
The calculation processes for the fractal dimensions of the seepage pores in the tested coal samples are shown in Figure 6.The fractal dimensions D3, D4, and D5 were calculated from the three parts that were divided by the relationships between the ln( / ) dV dP and ln(P) for the three coal samples, based on Equation (3).Previous studies demonstrated that D3, D4, and D5 were related to coal particle shape, internal seepage pore structure, and coal matrix compressibility [15,45], respectively.The values of the three fractal dimensions for all coal samples are listed in Table 5.From the fractal theory model, the fractal dimension Di ranged from 2 to 3 and is a suitable parameter for discussing the fractal characteristics of the coal matrixes surface [16].The value for D3 was below 2 (avg.1.80) except for the untreated 1/3CC coal sample, while all values for D5 exceeded 3 (avg.3.79).D3 and D5 were outside the range of the denoted fractal theory, therefore these dimensions are meaningless in this section.However, fractal dimension D4, which corresponds to intra-particle filling and can describe the mesopore and macropore matrix surface features, was between 2 and 3 for all samples, and thus D4 can be used to determine the characteristics of the seepage pores [45].The value of D4 for the original coal sample ranged from 2.84 to 2.89 (avg.2.86), while that for the treated samples increased from 2.88 to 2.97 (avg.2.91) at an average rate of increase of 1.62%.The treated coal sample had a more complex morphological pore structure under cryogenic freezing conditions, and the roughness and irregularity of the coal pore surfaces were the most prominent characteristics of the treated coal sample.Combined with the pore structure results from Table 4, the larger the APD, SSA, and TPV, the higher the value of D 4 .These features indicate that the coal pore structure was transformed, and the seepage pore structure was more complicated [5,15].

Calculated Degrees of Correlation
From Tables 4 and 5, the grey systematic behavior sequences X 0 , X 0 can be calculated by the equalization and dimensionless data transformation of the fractal dimensions D j (j = 1, 2, 4).
The internal effect factor sequences X i , X i (i = 1, 2, 3, 4, 5) can also be calculated from the coal pore structure parameters in Table 4 following the same method.The results are as follows [34].
The differences in the absolute values of X 0 , X 0 and the corresponding X i , X i are obtained (Table 6), and these values are defined as ∆ i (k) and ∆ i (k), respectively.Table 6 shows that the maximum and minimum values of the absolute difference between |X 0 − X i | and X 0 − X i are as follows [36]: The absolute difference between the grey systematic behavior sequence and its internal correlation factors sequence.

Sample
No. In grey relational theory, ρ is used as a constant and its value is usually 0.5, therefore, by bringing the data for the highest and lowest absolute difference values into Equation (4), we can calculate the correlation coefficients ζ i (k) and ζ i (k), respectively: From combining the data from Table 5 with the above formulae, the results are as follows [34]:   5), the systematic grey correlation degrees γ i and γ i can be calculated as follows, respectively [37]: The section above presents the process of calculating the grey correlation degrees γ i and γ i between the fractal dimension and pore structure parameters for the initial and LqCO 2 -treated coal samples.Next, we will describe the relationship between the fractal dimension and pore structure parameters based on the grey relational theory.γ i and γ i can be sorted in descending order, as follows: From comparing the sizes of γ i and γ i , it can be seen that the degree of correlation between the specific surface area of the raw and LCO 2 treated coal samples was largest, followed by the that of the full aperture, porosity and total pore volume, and permeability.Similarly, when the tested coal samples were treated by LqCO 2 , γ i was larger than the initial coal sample's grey correlation degree γ i .According to the grey relational theory, of correlation degree γ is higher, the internal behavioral factors will have a greater influence on the system [36].The above analysis indicates that γ i and γ i are positively correlated with D i and D i , respectively.Therefore, the coal matrix surface and internal pore structure of the treated coal samples became more rough and complex (i.e., fractal characteristics are more evident).
The size of γ for homologous pore structures shows that the pore specific surface area, average pore diameter, and porosity have the greatest effect on the roughness and irregularity of the coal matrix surface.Therefore, total pore volume and permeability have less of an effect on the coal, and this conclusion fully considers the rate of increase for each pore structure parameter (i.e., specific surface area, average pore diameter, porosity, total pore volume, and permeability).

Conclusions
The structure of coal pores significantly influences the adsorption and seepage of CBM as liquid CO 2 was injected into coal seam.The LP-N 2 -Ad and MIP methods were combined to investigate the full aperture pore structure evolution for original and treated coal samples.Furthermore, based on the grey theory, the relationship between the fractal dimension and pore structure parameters were quantitively analyzed.Some principal conclusions were as follows.
(1) The adsorption isotherms of three coal specimens were of type B, which illustrates that the coal samples contained numerous cylindrical shaped-pores with one closed side and slit-shaped and bottle-shaped pores.From the LP-N 2 -Ad and MIP tests, the APD (average growth rate of 18.20%), SSA (average growth rate of 7.38%), and TPV (average growth rate of 18.26%) were higher after the coal specimens were infiltrated by liquid CO 2 .This is because of the large number of new pores generated and plenty of micropores and transition pores transferred into mesopores and macropores.Therefore, the adsorption ability was lower while the seepage capacity was higher, which is suitable for CBM recovery.(2) Fractal dimensions D 1 (average of 2.58), D 2 (average of 2.90), and D 4 (average of 2.91) exhibited the same tendencies for the treated coal samples, which were typically higher than those of the original specimens (D 1 , average of 2.55, D 2 , average of 2.87, and D 4 , average of 2.86), indicating that the coal surfaces are rougher and the internal pore structures are more complex after the coal samples were treated by LCO 2 .(3) The grey relational theory was applied to analyze the relationship between fractal dimension and coal pore structure parameters.The correlation degree γ i was higher than γ i , indicating that the fractal features of the treated coal specimens are more evident.The degree of correlation between the fractal dimension and pore structure parameters show that the SSA of the raw and treated coal samples was largest, followed by APD, porosity, TPV, and permeability.The SSA, APD, and porosity positively influence the fractal characteristics of coal samples, but TPV and permeability exert negative influences.

Figure 1 .
Figure 1.Nitrogen adsorption and desorption curves for (a) the raw coal sample, and (b) the liquid CO2-infiltrated coal samples.

Figure 1 .
Figure 1.Nitrogen adsorption and desorption curves for (a) the raw coal sample, and (b) the liquid CO 2 -infiltrated coal samples.

Figure 4 .
Figure 4. Mercury intrusion and extrusion curves for the untreated and treated coal samples, (a) MHC coal sample, (b) MSC coal sample, (c) 1/3CC coal sample.

Figure 4 .
Figure 4. Mercury intrusion and extrusion curves for the untreated and treated coal samples, (a) MHC coal sample, (b) MSC coal sample, (c) 1/3CC coal sample.

Figure 5 .
Figure 5. Incremental pore volume distributions between the untreated and LqCO2-treated coal samples obtained from the MIP test, (a) MHC coal sample, (b) MSC coal sample, (c) 1/3CC coal sample.

Figure 6 .Table 5 .
Figure 6.Fractal characterization using the MIP results for (a) raw and (b) treated coal samples.

2 .
Characterization of the Fractal Dimension and Coal Pore Structure Parameters by the Grey Correlation Degree

Table 1 .
General information for tested samples.

Table 2 .
Structural coal pore parameters obtained using LP-N 2 -Ad.

Table 3 .
Fractal dimensions obtained from the LP-N 2 -Ad test.

Table 4 .
Porosity and pore parameter (specific surface area, pore volume, average pore diameter) results obtained using MIP.porosity; K, permeability; S me+ma , total specific surface area of mesopores and macropores; S me , specific surface area; S ma , specific surface area; V me+ma , total pore volume of mesopores and macropores per unit mass; V me , pore volume of mesopores per unit mass; V ma , pore volume of macropores per unit mass.

Table 4 .
Porosity and pore parameter (specific surface area, pore volume, average pore diameter) results obtained using MIP.Sme+ma, total specific surface area of mesopores and macropores; Sme, specific surface area; Sma, specific surface area; Vme+ma, total pore volume of mesopores and macropores per unit mass; Vme, pore volume of mesopores per unit mass; Vma, pore volume of macropores per unit mass.
Then, by taking ζ i (k) and ζ i (k) into Equation (