Coal Reservoir Characterization in a Tectonic Setting and the Effects of Tectonism on the Coalbed Methane ( CBM ) Content

Research on the relationships among tectonics, micropores, microfractures, and coalbed methane (CBM) content is important for the optimal selection of CBM production areas. In this study, micropore-microfracture structural parameters of coal samples from the Guojiahe coalfield are determined through the use of X-ray photography, an image recognition algorithm, and a liquid nitrogen adsorption method. ,e relationship between the micropore-microfracture characteristics of the reservoir and the gas content is quantitatively assessed using the Grassberger and Procaccia (GP) algorithm to calculate the correlation dimension of the parameters. Micropore-microfracture development varies in different tectonic zones. Additionally, the CBM content varies according to the characteristic parameters of hysteresis loops and the pore diameter. ,e correlation dimension is an effective indicator of the nonlinear relationship between reservoir micropore-microfracture characteristics and the gas content.


Introduction
Various geological factors, including tectonics, coal reservoir properties, and coal seam thickness, can strongly impact enhanced coalbed methane (CBM) extraction [1].In particular, coal reservoirs that have experienced multistage tectonic events show obvious property differences in terms of pore-fracture systems [2][3][4].
us, studies of the coal reservoir microfeatures associated with tectonics-induced effects on the CBM content could be conducted to optimize CBM extraction projects.Currently published results mainly focus on the structural characteristics of micropores and microfractures for use in the evaluation of the CBM exploitation potential [5].Furthermore, previous publications have largely focused on the pore characteristics of reservoirs and the reservoir structure without considering the connection between variations in pore characteristics and tectonics [6,7].Other studies have focused on CBM accumulation only from a tectonic perspective.Notably, previous studies have found that CBM is enriched above the neutral planes of synclines and below those of anticlines, due to compression forces [8].Additionally, extensional faults provide channels for gas emissions, and tectonism is the main factor that causes variations in the pore-fracture structures of coal reservoirs and results in nonuniform distributions of CBM.Specifically, the pore-fracture features vary based on structural forms and affect gas migration, accumulation, and preservation.us, it is worth considering pore-fracture characteristics in different types of tectonic zones and the associated effects on the gas content.
Coal reservoirs comprise pores and fractures [9], which can be subdivided into macrofractures, microfractures, macropores, and micropores according to pore-fracture size [10][11][12].e hybrid methods used to study reservoirs include liquid nitrogen adsorption or carbon dioxide adsorption [13,14], mercury injection experiments [15], scanning electron microscopy (SEM) [16], low-field nuclear magnetic resonance (NMR) [17][18][19].Among these, experiments involving liquid nitrogen adsorption and carbon dioxide adsorption can obtain estimates of different pore sizes and percentages.Adsorption pores strongly affect the gas storage capacity of coals through the pore structure system and gas diffusion through the fracture structure system.erefore, these results provide important guidance for gas extraction from coal seams and for gas-related hazard prevention during coal mining.SEM images can show pores or the fractures on the surface of coal samples.e internal fracture characterization of a coal mass can produce an estimate of the 3D structure [20,21], which can be quantitatively analyzed by means of X-ray penetration using digital radiography (DR) systems [20].Yao et al. [17] presented a new approach of combining NMR and X-ray computed tomography to quantitatively describe the types of size distribution and spatial arrangement of pores and fractures.
In this paper, liquid nitrogen adsorption and a DR system are combined with an X-ray image recognition algorithm to identify the micropore and microfracture structural characteristics of coal reservoir samples collected from geological structures.Based on the Grassberger and Procaccia (GP) algorithm [22], various micropore and microfracture system parameters and CBM contents are calculated to determine the micropore and microfracture structure of a coal reservoir, as well as the associated effects of this structure on the gas content of the reservoir.

Regional Geotectonic Characteristics
In this paper, the Guojiahe coalfield in the Yonglong mining area, China, is selected as the case study.e coalfield is located in the western Miaobin Depression and north of the Weibei Flexure Belt tectonic unit (Figure 1).e regional strata consist of the Tongchuan, Fuxian, Yanan, Zhiluo, Anding, Yijun, Luohe, and Huachi Formations, as well as Neogene, Middle-Late Pleistocene, and Holocene units, in chronological order.e Middle Jurassic Yanan Formation includes coal-bearing strata, and the thickness of the main No. 3 coal seam ranges from 0.55 m to 26.83 m, with an average thickness of 11.88 m in the area.
e Yanan Coal Formation was formed by the Indosinian Movement and the first phase of the Yanshan Movement [23].e folded strata dip gently to the NW-NNW and strike NE-NEE.ese characteristics are controlled by the N-E (35 °) paleostructure created by the Indosinian Movement [24,25].e Liangting anticline extends inward to the coalfield.e anticline axis not only resulted in the nondeposition of the Yanan Formation but also in thin coal deposits, with a thickness of 2 to 5 m.However, the Yanan Formation exhibits good coal thickness of 15.2 to 27.9 m along the axis of the Caizigou-Zhangba syncline (Figure 1).us, the distribution of the coal seam thickness is controlled by the paleostructure and gradually thins toward the anticline and thickens toward the syncline.

Fracture Identification Method with the DR System.
DR is an effective method of identifying the microfeatures of coal samples [20]. is method can quantify fractures by projecting a three-dimensional fracture network onto a plane.In this paper, the microfeatures of twenty-five coal samples from eleven coal sampling locations are studied using an Angell-1500R DR system (Figure 2).Two to three samples were taken at each location, and twenty-five samples were taken from the different tectonic zones of the Guojiahe coalfield in China, as shown in Figure 3. e system parameters of the DR system are presented in Table 1.

Fracture Characterization of Coal Samples.
Quantitative analysis is the basic research method for determining the gas distribution in fractured coal.e fractal dimension [26], joint roughness coefficient (JRC) [27], and rock mass system (RMS) [28] methods have been commonly used to measure the macrofracture and microfracture characteristics of a coal-rock mass.Image processing software and the Kirsch operator image analysis method are used in this study to identify fracture characteristics and express them as parameter C i .e equation for C i is as follows: where h max (mm) is the maximum gap between tectonic fractures in a plane projection, l (mm) is the length of tectonic fractures in a plane projection, and i is the number of fractures.Image progressing software (Image-Pro) and the Kirsch operator image analysis method are used to detect the edges and identify the parameters of microfractures obtained by X-ray photography (Figure 4).en, the microfracture parameter C i is calculated.
e structure of microfractures obtained by the DR system appears as veins or clusters (Figure 5).Among the twenty-five coal samples, five samples from the points GJH3-X2, GJH3-01, and GJH3-02 exhibited no structural fractures.Additionally, three samples from points GJH3-03 and GJH3-04 had fracture lengths shorter than 20 mm and a maximum gap of less than 0.5 mm and are not listed because fractures with a reasonable size provide seepage channels rather than emission channels.
en, based on the proportion of fractures longer than 20 mm and wider than 0.5 mm, six coal samples out of the remaining seventeen were selected as representative samples for each point to calculate the fracture parameter.e maximum gap ranged from 0.54 to 6.94 mm, the length ranged from 23.16 to 130.1 mm, and the linear density ranged from 8.28 to 23.7 fractures per meter.Because tectonic fractures of reasonable size are the main channels for CBM migration [29], the maximum projected area of fractures was used to calculate C i .
A detailed analysis was performed based on the Kirsch method, and the pixel dimension of images of representative samples is shown in Figure 5.
e relationships between microfracture parameters and tectonic zones indicate that microfracture characteristics vary in different zones.Specifically, small C i values reflect a longer projected fracture length and a low level of microfracture development.
e C i values of GJH3-F1 to GJH3-F4 vary from 0.034 to 0.059, and the fracture density ranges from 16.93 to 23.7 in the fault area, as shown in Figure 6.ese values indicate that faults have a considerable influence on the density and 2 Advances in Materials Science and Engineering development of microfractures.e C i values di er in different tectonic zones, and the B1 anticline has a notable in uence on coal microfracture development.e fracture density is relatively small near the X1 syncline and B1 anticline, which indicates that tectonic movement may have little e ect on the microfeatures of coal samples in syncline and anticline zones.In conclusion, the development degree of microfractures varies based on the tectonic zone.

Pore Characterization of Coal Samples.
e majority of CBM accumulation occurs in coal pores, and the pore structure parameters are key indicators that can be used to evaluate coal seam adsorption.Speci cally, liquid nitrogen experiments can be performed to investigate coal pore characteristics [14].
e method is suitable for pore diameters ranging from 1 to 20 nm, and this method has a wider measuring range than the method based on carbon dioxide.ASAP JW-BK122 W is used to obtain the relative pressure and adsorption curve under isothermal conditions and classify the adsorption-desorption isotherms.e slope K i and length L i of a straight line passing through an initial point (P′/P 0 , V′) and end point (P * /P 0 , V * ) of an isotherm can be calculated as follows: where K i represents the adsorption/desorption value (Δcm 3 /g) associated with a relative pressure change during the capillary and pore condensation process.L i represents the amount of condensation from capillary pores relative to that from all pores.e adsorption volume associated with a relative pressure depends on the minimum pore size at the   Liquid nitrogen experiments were conducted on seven samples from the syncline, anticline, and fault areas, and 7 isothermal curves with clear hysteresis loops related to the pore shape and size were obtained.us, the pore shape, diameter, and distribution were determined by analyzing the hysteresis loop features and calculating the proportion of di erent pore sizes.
Figure 7 shows the hysteresis loops in di erent tectonic zones, and the feature parameters K i and L i are calculated based on the above equation.e second type includes small closed loops with initial points at a relative pressure of 0.49 and small K i values.is type of loop is associated with anticline samples.e third type of loop exhibits start points at a relative pressure of 0.17, and K i varies from 3.23 to 14.57.e third type of loop is observed for fault samples. is classi cation has also been used in previous studies of hysteresis loops for di erent coal body structures [30][31][32].In some cases, isothermal curves may not form closed loops because liquid nitrogen remained in open-ended micropores with a decrease in relative pressure.us, it is scienti cally reasonable to de ne a point on the adsorption curve where the adsorption volume is similar to the desorption value at low relative pressure as the initial point of K i.
According to the Wpept, C. C. classi cation, the pores are mainly micropores with diameters of 3.55 nm to 8.99 nm (average of 5.09 nm).e Brunauer-Emmett-Teller (BET) speci c surface area ranges from 2.26-15.69m 2 •g −1 , with an average of 10.16 m 2 •g −1 .e Barrett-Joyner-Halenda (BJH) pore capacity ranges from 0.00475 to 0.02122 cm 3 •g −1 , with an average of 0.0101 cm 3 •g −1 .e pore diameters of coal samples in syncline and anticline areas exhibited a small di erence, but those in the fault area were considerably di erent.Speci cally, the pore diameters of coal samples from the fault zone were larger than those from syncline or anticline areas.e proportions of di erent pore sizes in di erent tectonic samples are shown in Table 3. e pore capacity and diameter class of greater than 2 nm is abundant in the fault zone area but not in the syncline and anticline areas.e pore size distribution may be related to the tectonic forces in the area because such forces may change the pore and fracture structures by destroying primary pores and developing apertures, mesopores, and macropores.In conclusion, di erent tectonic zones in uence the pore diameter and distribution in di erent ways.

Correlation Dimension Calculation with the GP Algorithm.
e GP algorithm is one e ective method of calculating the correlation dimension of chaotic time series [33].Correlation dimension analysis is an important branch of fractal theory, and this dimension is used as an indicator to evaluate complex trends in data.e microfracture parameters, pore characteristic parameters, and gas content were investigated based on a series of i collected samples: Selecting a subset of values from the above equation yields an m-dimensional phase space N 1 : e number of values in the m-dimensional phase space is n − m + 1. us, the correlation function is as follows: In equation ( 5), H(x) represents the Heaviside function: where r has a given value of r kr 0 .In this relation, k is a segment of radius r, and r is equal to 8 in this paper.Additionally, the minimum embedding dimension is 1, the maximum embedding dimension is 10, and the delay time is 2. e correlation dimension is the slope of the straight line t to ln (r) and ln C(r) plotted in log-log coordinates:

CBM Content Analysis in Di erent Types of Tectonic
Zones.Structural geologic factors, such as coal thickness, burial depth, and fault distribution, are important factors in uencing CBM accumulation [34][35][36].Coal microfeature evolution is in uenced by multiphase tectonic activity, which can lead to gas migration, di usion, and concentration variations.In this paper, panel I of Guojiahe coal eld is used as an engineering example, and the microporemicrofracture parameters and gas content are analyzed.
Based on the GP algorithm, the correlation dimension is calculated to quantitatively analyze the relationship between tectonics and the gas content.e results of this analysis can be used to optimize CBM production in di erent zones.
Coal microfeatures are related to the gas content.Based on previous experimental data collected from a borehole during the exploitation period, the contours of the CBM content are illustrated in Figure 3. Microfeature data obtained using the DR system and via liquid nitrogen experiments and the gas contents of coal samples are listed in Table 4. Gaussian tting is an e ective method of analyzing the relationships among these data, and the tting results are shown in Figure 8.
e microfracture parameter C i re ects the development of coal fractures, which in uence the migration and emission of CBM.Speci cally, gas migrates and spreads along coal fractures, and this process often decreases the gas content.
e relationship between the microfracture parameter C i and the gas content was analyzed based on Gaussian tting, and the results are shown in Figure 8(a).C i reaches a peak at 0.154, when the gas content is 2.365 m 3 /t, and the average value of C i is 0.055.When the gas content is greater than 2.52 m 3 /t, C i remains small.is result suggests that the tectonic fractures in the studied plane projection are relatively long, and long fractures are generally associated with low gas contents.Figure 8(a) shows that the gas content is high in more than half of the samples, but the C i values are less than 0.06.e Gaussian curve exhibits a small width and high peak, which indicates that data are centrally distributed and that C i is highly correlated with the gas content of samples from di erent tectonic zones.Advances in Materials Science and Engineering when the gas content is 2.24 m 3 /t.L i reaches a peak at 9.37 when the gas content is 2.19 m 3 /t.e hysteresis loop parameters K i and L i must be sufficiently small.erefore, the distance between the initial point and the ending point must be decreased so that the gas content remains relatively high.In conclusion, the development of microfractures inhibits gas accumulation.Specifically, the more expansive the fracture development is, the lower the gas content is.Conversely, the development of micropores is good for gas accumulation, yields small values for the hysteresis loop parameters K i and L i , and results in high gas contents.

Correlation Dimension for Multiparameter Calculations.
e correlation dimension is a key factor used to quantify changes in chaotic time series, and correlation dimension analysis is an important branch of fractal theory.e GP algorithm uses a correlation integral to calculate the correlation values of variables and determine the regularity and degree of change.us, this approach reconstructs the phase space.resholds must be defined for different parameters.In this approach, two parameters are considered related when the difference in dimension is smaller than the threshold, and the opposite is true as well.Correlated parameters are identified and grouped in phase space.e greater the number of parameters in a group, the higher the associated correlation is.e microfracture parameter C i , hysteresis loop parameters K i and L i , and gas content q i are selected for analysis.In addition, 5 of the 103 characteristic parameters plotted in the Gaussian curve are added to the database.MATLAB 2014R software is used to analyze the relationship between the time-series correlation integral C(r) and distance r in the reconstructed phase space.en, the correlation dimensions are calculated after analyzing the time series using the delay embedding method.e correlation dimensions are used as indicators to quantitatively evaluate the effects of the microfeatures of the coal reservoir on the gas content.
e GP algorithm and delay embedding method combined with the saturation correlation dimension are used to calculate the minimum embedding dimension.Specifically, the minimum embedding dimension is expressed as m � m 0 + 1, when the estimated correlation dimension D m hardly varies with an increase in the embedding dimension.Based on the relationship between correlation integral C(r) and distance r in the reconstructed phase space, the slope of the nonscaling interval on the fitting line is the correlation dimension based on the least square method.
Figure 9(a) shows that the estimated correlation dimension linearly increasing when m 0 is 1, 2, 3, or 4. e value of the estimated correlation dimension becomes stable when m 0 is larger than 7. us, the embedding dimension is 8. Based on dividing the lnC(r) -ln(r) curve into sections, the slope of the A-B interval defined as the correlation dimension D is 2.17, as shown in Figure 9(b).
e relationship between micropore and microfracture parameters and CBM content can be quantitatively defined by the correlation dimension.Figure 10(a) shows that the correlation dimension is 0.425 in the reconstructed phase space consisting of the microfracture parameter C i and the CBM content q i .e correlation dimension is 3.57 in the reconstructed phase space consisting of the microfracture parameters K i and L i and the CBM content q i as shown in Figure 10(b).
e results show that the hysteresis loop parameters K i and L i have a dramatic impact on the CBM content q i .e calculation results provide a quantitative assessment of the relationships among tectonic structures, the microfeatures of the coal reservoir, and the CBM content.

Conclusion
(1) e microfracture characteristics of coal samples can be quantitatively described based on the Kirsch recognition algorithm and X-ray images.Advances in Materials Science and Engineering Additionally, C i is presented based on microfracture features, and the results indicate that the microfracture degree of coal samples obviously differs in tectonic zones.(2) e N 2 adsorption/desorption isotherms of coal samples were defined as three types associated with syncline, anticline, and fault areas.Notably, the pore size and pore distribution of coal specimens from fault areas were observably different than those of coal samples from syncline and anticline areas.ese results suggest that the pore structure system evolves with tectonic stress.(3) e correlation dimension was calculated with the GP algorithm.e relationships among the experimental time series of the microfracture parameters of coal samples from different tectonic zones were analyzed.Additionally, the hysteresis loop parameters of adsorption/desorption isotherms and the CBM content were considered.e micropore and microfracture characteristics of coal samples around different tectonic zones exhibited a nonlinear relationship with the CBM content, and this relation could be quantitatively defined as a correlation dimension.

Data Availability
e data used to support the findings of this study are included within the article.Figure 10: Correlation dimension of micropore and microfracture parameters: (a) lnC(r) and lnr curve in log-log coordinates in the reconstructed phase space of C i and q i ; (b) lnC(r) and lnr curve in log-log coordinates in the reconstructed phase space of K i , L i , and q i .Advances in Materials Science and Engineering
e calculation results are shown in Table 2. e hysteresis loops are classi ed into three types based on their isotherm features.e rst type includes large closed loops with initial points at relative pressures of 0.4 to 0.42 and K i values from 4.35 to 6.94 (Figures 7(a) and 7(b)). is type of loop is observed for syncline samples.

Figure 4 :Figure 5 :
Figure 4: Microfracture parameter identi cation for coal specimens based on an X-ray image.

Figure 6 :
Figure 6: Microfracture parameters of coal specimens in di erent tectonic zones.

R 2 =Figure 8 :
Figure 8: Analysis of the relationships of the CBM content and parameters of micropores and microfractures: (a)relations of parameter C i and CBM content; (b) relations of adsorption, desorption parameters, and CBM content.

Figure 9 :
Figure 9: Correlation dimension in the reconstructed phase space: (a) relationship between the reconstructed phase space and embedding dimension; (b) lnC(r) and lnr curve in log-log coordinates.

Table 1 :
Sampling locations and voltage parameters of the DR system.
initial point.e length of the process depends on the pore sizes and distribution of pores.

Table 2 :
Parameter variations based on the hysteresis loops.

Table 3 :
Analysis results for the specific surface area and low-temperature nitrogen adsorption.: the letter X in the numbers of coal samples represents the syncline area, B represents the anticline area, and F represents the fault area.
Note* Calculated by pore volume distribution.

Table 4 :
Gas content and micropore-microfracture parameters in different tectonic zones.