Abstract

The pore-fracture structure of deep coal reservoirs is highly important for evaluating, exploring, and developing coalbed methane (CBM) resources. This study considers three coal samples from the DJ57 well in the Benxi Formation in the Daning–Jixian block on the eastern margin of the Ordos Basin as the research object. Based on the coal quality parameters of the coal samples, field emission scanning electron microscopy (FE-SEM), gas adsorption experiments, high-pressure mercury intrusion porosimetry (MIP), and microcomputed tomography (micro-CT) scanning were used to quantitatively characterize the nanoscale pores and microscale fractures in deep coal reservoirs and to evaluate the pore-fracture structure at different scales. The results reveal that the pore types in the Benxi Formation coal samples are diverse and include mainly organic matter (OM) pores, inorganic pores (intraparticle and interparticle pores), and microfractures. The organic pores are diverse in shape and mainly exhibit round, oval, and wedge shapes, while the microfractures exhibit slender stripes or serrated curves. The multiscale quantitative characterization of deep coal reservoir pores and fractures is based on a variety of pore characterization methods, and the pore and fracture volume distributions are mainly U-shaped, revealing the coexistence of micropores and microfractures. The volumes of micropores (0.3–2 nm), mesopores (2–50 nm), macropores (50 nm to 10 μm), and microfractures (>10 μm) account for 78.00%, 6.78%, 2.08%, and 13.14%, respectively, of the total pore volume (PV). Based on a full-scale pore-fracture splicing calculation, the total permeability of the Benxi Formation coal samples ranges from 5.77 to 28.22 mD. The observation results indicate that the microfractures are connected to each other, forming a network structure with strong connectivity. The microfractures are mainly associated with pore μm, accounting for approximately 95% of the total permeability. Moreover, micropores in deep coal reservoirs provide a large space for CBM adsorption, and microfractures enhance the seepage capacity of CBM.

1. Introduction

In recent years, China’s CBM exploration and development efforts have focused mainly on shallow and medium layers ( m). As exploration and development continue to progress, increasingly fewer blocks can be developed in a cost-effective manner, and deep CBM will become an important successor field for China’s next step in the exploration and development of CBM [1, 2]. China’s deep CBM resources are abundant. According to the latest round of CBM resource evaluation by the Ministry of Land and Resources, the CBM resources 1000–2000 m deep (medium-depth layer) are  m3, and it is estimated that the national CBM resources deeper than 2000 m (deep layer) are approximately  m3, of which the resources at depths of 2000–3000 m are  m3 [3, 4]. With the strong investment of state and major oil and gas companies, the exploration and development of deep CBM in recent years have been effective, and the emergence of local high-yield wells on the eastern edge of the Ordos Basin, the southern part of the Qingshui Basin, the eastern part of the Junggar Basin, etc., indicates that China’s deep CBM has great potential for development [1, 3].

The pore size distribution (PSD) of coal varies greatly, and the pore characteristics of different sizes and their contributions to CBM adsorption, desorption, diffusion, and migration also differ [57]. Investigators mainly classify pore size according to gas occurrence state and divide the pores in coal into micropores (0.3–2 nm), mesopores (2–50 nm), macropores (50 nm to 10 μm), and microfractures (>10 μm) [811]. CBM mainly exists in the adsorbed state on pore surfaces and in multiscale (nano-, micro-, and millimetre scales) and multiform connected pore and fissure systems, which act as gas migration and transmission channels. The gases are continuously produced through three stages—desorption, diffusion, and seepage—and diffusion and seepage in nanopores are the initial stages of gas desorption and migration [12, 13]. In CBM development, more attention has been given to the impact of the microfracture scale and connectivity on reservoir permeability, but the development characteristics and connectivity of nanoscale pores may constrain and potentially impact gas desorption and migration [13].

Previous studies on coal pore structure have been carried out in many ways and can be divided into three categories according to experimental principles: fluid injection methods, image analysis methods, and nonfluid injection methods [14, 15]. Fluid intrusion methods mainly include low-pressure CO2 adsorption (LP-CO2GA), low-temperature N2 adsorption (LT-N2GA), and high-pressure MIP methods [1620]. Image analysis methods mainly include scanning electron microscopy (SEM), atomic force microscopy (AFM), and transmission electron microscopy (TEM) [19, 21, 22]. Nonfluid injection methods mainly include computed tomography (CT), nuclear magnetic resonance (NMR), small-angle X-ray scattering (SAXS), and small-angle neutron scattering (SANS) [18, 2326]. Currently, the study of the pore-fracture structure of coal and shale is mainly based on the combination of high-pressure MIP, LT-N2GA, and LP-CO2GA methods [16, 18, 19]. A study based on this approach revealed that the pore diameter of marine shales in the Sichuan Basin is generally 10-50 nm and that the pore diameter of transitional shales in the eastern margin of the Ordos Basin is generally 10-20 nm, indicating that the PV of marine and transitional shales is mainly provided by mesopores [2729]. The pore diameters of coal samples from the Shanxi Formation and Taiyuan Formation are mainly 0.3-1.5 nm, indicating that the PV of these coal samples is mainly provided by micropores [13, 15]. However, the influence of macropores and microfractures in coal reservoirs has not been sufficiently emphasized, resulting in uncertainty in our understanding of the contribution of macropores and microfractures to the reservoir capacity of coal reservoirs. In addition, in high-pressure MIP experiments, due to the compressibility of coal samples, an increase in mercury intrusion pressure leads to compressive deformation of the coal matrix and pore damage, reducing the accuracy of the test results.

Since 2019, deep CBM exploration in the Daning–Jixian block has advanced to a depth of more than 2000 m, and the Jishen 6-7-P-01 well has obtained a high production industrial gas flow of  m3 per day [2, 3]. In 2021, the known reserves of the Daning–Jixian block were  m3, making it the first large CBM field in China with high abundance below a depth of 2000 m and with reserves exceeding  m3, indicating that deep CBM in the Daning–Jixian block has promising development potential [1]. LP-CO2GA, LT-N2GA, high-pressure MIP, and micro-CT data tests were carried out on three coal samples from the DJ57 Jingbenxi Formation in the Daning–Jixian block on the eastern margin of the Ordos Basin to identify and quantitatively characterize the pores and fractures of coal reservoirs in the Daning–Jixian block. Additionally, physical experiments on the helium porosity and permeability of the samples were carried out, and the experimental results were compared with numerical results and analysed to determine the influence of microscopic pores and fissures on the occurrence and seepage of CBM in deep coal reservoirs to provide a reference for the exploration and development of deep CBM in the eastern margin of the Ordos Basin.

2. Geological Setting

The Ordos Basin is a typical cratonic marginal superimposed basin in the western part of the North China Plate. The basic tectonic framework was formed in the Yanshanian period and developed in the Himalayan period. Overall, the Ordos Basin features a large syncline with a gentle and long east wing, a steep and short west wing, and upturned north and south wings [30]. According to its underlying properties, structural morphology, and geological evolution history, the Ordos Basin is divided into six first-order tectonic units, namely, the Yimeng Uplift, West Margin thrust belt, Weibei Uplift, Tianhuan Depression, Jinxi fault belt, and Yishan slope (Figure 1(a)) [30]. The research area is located in the Daning–Jixian area on the eastern margin of the Ordos Basin. The Daning–Jixian block is located at the southern end of the Jinxi fault belt and the southeastern margin of the Yishan slope on the eastern margin of the Ordos Basin, and its southern neighbour is the Yanchuan block.

The coal-bearing strata in the study area are mainly in the upper Palaeozoic Carboniferous system and Permian system, and the corresponding sedimentary environment changes from marine to continental facies; these strata can be divided into the Shanxi Formation, Taiyuan Formation, and Benxi Formation. The Shanxi Formation is dominated by delta facies, while the Taiyuan Formation and Benxi Formation contain delta, lagoon, tidal flat, and shallow shelf deposits from north to south. In the vertical direction, mudstone is alternately distributed with coal and sandstone, and pure mudstone rarely develops [30, 31]. The main target layers are the Shanxi Formation coal seam, Taiyuan Formation coal seam, and Benxi Formation coal seam [31]. The thickness of the Shanxi Formation coal seam is only 1–3 m, the thickness of the Taiyuan Formation coal seam is approximately 1–6 m, and the thickness of the Benxi Formation coal seam is greater, ranging from 5 to 12 m; these seams are the main coal seams in the Ordos Basin. The no. 8 coal seam of the Benxi Formation has a wide and continuous distribution and large thickness; among these seams, the plane spreading area between 2000 and 3500 m is approximately  m2, which is the most favourable deep CBM exploration target of the Benxi Formation (Figures 1(b) and 1(c)) [31].

3. Samples and Methods

3.1. Sample Selection

Three fresh coal samples from newly drilled wells in the Daning–Jixian area in the eastern margin of the Ordos Basin were selected; all of these samples were obtained from the upper, middle, and lower parts of the Benxi Formation at the DJ57 well. To eliminate the influence of the heterogeneity on the experimental data, before the experiment, the core was cut by wire cutting equipment, and a cylindrical coal sample (50 mm high with a diameter of 25 mm) was cut for micro-CT scanning, helium porosity, and permeability measurements. Then, the cylindrical coal sample was cut and crushed. The resulting pieces were subjected to vitrinite reflectance () and maceral analysis, proximate analysis, FE-SEM, high-pressure MIP, LT-N2GA, and LP-CO2GA experiments (Figure 2). measurements (80 points) and maceral analysis (500 points) were performed randomly in oil-immersed reflected light according to the Chinese standards GB/T 6948–2008 and SY/T 6414–2014, respectively. The moisture content (), ash yield (), volatile yield (), and fixed carbon () content of the coal samples were analysed in accordance with the Chinese standard GB/T 212–2002. Helium porosity and permeability were measured in cylindrical coal samples according to the Chinese standard GB/T 34533–2017. The test results are shown in Table 1. The porosity and permeability tests of 3 coal samples from the Benxi Formation in the DJ57 well in the Daning–Jixian block reveal that the coal sample porosity ranges from 3.23% to 6.64% (mean 5.40%), and the permeability ranges from 1.42 to 11.01 mD (mean 6.93 mD).

3.2. Experimental Methods
3.2.1. FE-SEM

The qualitative characterization of coal samples was mainly realized through direct observation and description of pores via FE-SEM [14, 20]. For this, a small rectangular coal sample with dimensions of was prepared. Then, a HITACHI IM 4000 argon ion polishing instrument produced by Hitachi company was used to continuously bombard the cross-section of the vertical shale formation under a certain energy under vacuum conditions to obtain a super smooth surface. Argon ion polishing lasted approximately 3 hours. The sample preparation and polishing process can effectively avoid damage to the sample surface during the mechanical polishing process and can restore the true pore shape of the sample. The polished samples were imaged with a Helios NanoLab 660 double-beam scanning electron microscope produced by the FEI company to obtain high-resolution SEM images. Through a series of SEM images with different resolutions, the pore morphology, relationship between pores and particles, pore distribution, pore size, and pore connectivity can be directly observed. The experimental procedure was carried out in accordance with the Chinese Oil and Gas Industry Standard SY/T 5162–2014.

3.2.2. Gas Adsorption

The gas adsorption experiment was carried out by an AutoSorb-iQ-MP-C automatic physical chemical adsorption instrument. The particle size of the sample was 60-80 mesh. Before the analysis of the LT-N2GA experiment, the samples were degassed at 110°C for 12 hours to remove the water and volatile substances in the samples. Then, the degassed samples were moved to the analysis station, and the adsorption-desorption test was carried out at 77 K with high-purity N2 as the adsorbent. The structural parameters, including the PV, specific surface area (SSA), and PSD, were obtained by the Brunauer–Emmett–Teller (BET) and nonlocal density functional theory (NLDFT) models [3235]. The experimental procedure was carried out in strict accordance with GB/T 21650.2–2008. The sample pretreatment of the LP-CO2GA experiment was similar to that of the LT-N2GA experiment. A total of 1–2 g of 60–80 mesh powder sample was degassed for 16 hours, and the adsorption test was carried out at 273 K with high-purity CO2 as the adsorbent. The structural parameters, including the PV, SSA, and PSD, were obtained based on the NLDFT theoretical model [36, 37]. The experimental procedure was carried out in accordance with GB/T 21650.3–2011.

3.2.3. High-Pressure MIP

An AutoPore 9505 mercury injection instrument produced by Quantachrome Instruments, U.S., was used for the high-pressure MIP testing. During the test, the mercury injection pressure provided by this instrument reached 413 MPa, and the corresponding test pore diameter limit was approximately 3 nm. For this test, coal was processed into a cube of 1 cm3, and the surface of the coal sample was smoothed to avoid the “pitting effect.” The sample was also continuously dried for 8 hours at 105°C to ensure that the original pores of the coal sample were not destroyed by the high temperature, the impurities and gases inside the sample were effectively eliminated during the test, and the vacuum inside the instrument was maintained. According to the data automatically collected by the instrument during the process of mercury injection and mercury withdrawal, pore structure parameters such as the PV, SSA, and PSD were obtained by the Washburn equation [38]. The experimental process was carried out in accordance with GB/T 21650.1–2008.

3.2.4. Micro-CT Scanning

Micro-CT scanning is a commonly used test method that can quantitatively characterize the spatial distribution and connectivity of some macropores and microfractures in coal samples [22, 23]. A nanoVoxel3502E multiscale high-resolution X-ray three-dimensional microscopic imaging system produced by Sanying Precision Instruments Co., Ltd., was used to scan the coal samples, with a system resolution of more than 500 nm. Before the micro-CT scanning experiment, a cylindrical coal sample with a diameter of 25 mm and a length of approximately 5 cm was drilled along the vertical direction. After the cylindrical coal sample was fixed, it was placed vertically in the micro-CT scanning device. The X-ray source position was adjusted, and the sample was prepared for micro-CT scanning. The micro-CT scanning voltage was 60 kV, and the experimental temperature was 25°C. Due to the influence of the cylindrical coal sample thickness and exposure time, the minimum pore size characterized by the micro-CT approach was 10 μm.

4. Results

4.1. Pore-Fracture Morphology Analyses Based on FE-SEM Images

FE-SEM was used to observe and analyse the pore and fracture morphology of the Benxi Formation coal samples from the Daning–Jixian block. There are many nanometre-scale pores and microfractures in the Benxi Formation coal samples. At present, classification schemes for coal reservoir pores have not been unified [14, 20]. In this study, referring to the classification scheme proposed by Loucks et al. [7], the pores in the Benxi Formation coal samples are divided into OM pores, inorganic pores, and microfractures. The microfractures are shrinkage-induced fractures. OM pores are mainly distributed in the coal matrix (Figures 3(b)–3(d), 3(f), 3(h), and 3(i)), with various pore shapes, including regular round shapes and oval shapes (Figures 3(d), 3(f), and 3(i)) and irregular wedge shapes (Figures 3(h) and 3(c)), as well as dispersed and distributed spongy connected pores (Figure 3(b)). Moreover, there are isolated pores (Figures 3(d), 3(f), and 3(i)) with pore diameters ranging from 4.4 to 220 nm, which provide space for CBM and are important pore types for CBM reservoirs. Inorganic pores can be divided into intragranular pores and intergranular pores. The former is generally distributed between clay and mineral particles or between rigid particles (Figures 3(f) and 3(g)) and is mainly related to the preservation of primary pores, which are usually manifested as isolated pores with poor connectivity (Figures 3(c), 3(f), and 3(h)). Moreover, their pore diameter is larger, generally greater than 300 nm. The latter component usually develops in dissolved mineral particles and is poorly connected to the external environment (Figures 3(e) and 3(f)), and its pore diameter is generally small, between 20 and 50 nm. Microfractures are well developed in the Benxi Formation coal samples. Most of the microfractures are distributed inside the mineral particles or at the edges of the rigid particles (Figures 3(a), 3(d), and 3(e)). The microfractures are either elongated or zigzag-shaped and have a wide range of sizes. The microfracture lengths are approximately 3~50 μm, and the microfracture widths are approximately 20 nm to 1 μm.

Figure 3 shows the morphology of the pores and microfractures, and it can be concluded that (1) the OM pores are mainly circular, elliptical, and wedge-shaped pores; (2) the inorganic pores are mainly closed pores, and some are open pores (such as parallel plate pores and conical pores); and (3) the microfractures are mainly elongated or zigzag-shaped. The pore types of the middle-shallow coal seam mainly include initially developed pores, cavity pores, intercrystalline pores, and dissolution pores. The microfractures of the middle-shallow coal seam are relatively developed, and no filling material is found in the microfractures; horizontal bedding fractures are common there, so the permeability of the middle-shallow coal seam is greater than that of the deep coal reservoir.

4.2. Analysis of LP-CO2GA Data

Figure 4 shows the LP-CO2GA curves of the coal samples from the study area, which exhibit similar characteristics. The amount of CO2 adsorbed by these samples steadily increases with increasing pressure. According to the IUPAC classification, these curves are type I isotherms [39]. The adsorption of CO2 by sample DJ57-3 is the highest, reaching 20.07 cm3/g, indicating that this sample has more micropores than the other samples. The amount of CO2 adsorbed by DJ57-2 is the lowest, at only 13.82 cm3/g, indicating that there are relatively few micropores in the sample. Based on the calculation of the LP-CO2GA data with the NLDFT model, the PSD curves and pore structure parameters of the coal samples in the pore size range of 0.3–1.5 nm were obtained (Table 2). The SSA of the micropores ranges from 141.64 to 207.37 m2/g (mean 180.46 m2/g), which is much greater than the result of the LT-N2GA analysis. This result shows that the proportion of micropores in coal is high. The volume of the micropores ranges from 0.042 to 0.061 cm3/g (mean 0.053 cm3/g). The mean pore diameter of all the samples is 0.501 nm. The distribution characteristics of micropore size based on LP-CO2GA data show that the micropore development of coal samples generally presents a multimodal peak distribution of 0.3–0.4 nm, 0.4–0.7 nm, and 0.8–0.9 nm (Figure 5). The increase in PV and SSA in the >1.0 nm pore diameter segment slows significantly, and the corresponding PV and SSA are very low. Because the micropores in the coal samples are relatively developed and have a large SSA, they can provide many adsorption sites for gas occurrence [13, 40].

4.3. Analysis of LT-N2GA Data

Figure 6 shows the N2 adsorption-desorption curves of the coal samples from the study area. The morphology of the adsorption curves of all the coal samples slightly varies, but the overall shape is a reverse “S.” According to the six types of physical adsorption isotherms proposed by the IUPAC, the N2 adsorption curve is close to the type IV adsorption isotherm [39]. In the low-pressure stage (), the gas adsorption capacity increases rapidly with increasing relative pressure. In this stage, micropore filling occurs, and single-molecule adsorption is the main adsorption mechanism, indicating that the gas adsorption capacity at this stage is mainly related to micropores and some mesopores [13]. When the relative pressure further increases (), the adsorption amount of N2 increases slowly, indicating that the single-molecular layer adsorption process has ended and that the multimolecular layer adsorption process has begun. When the relative pressure continues to increase to nearly 1.0, the adsorption capacity increases rapidly again, but there is no adsorption when the relative pressure approaches the saturated vapour pressure (), indicating that there is a certain amount of open macropores in the coal samples.

Previous studies have shown that the pore morphology of coal can be further evaluated by the type of N2 adsorption-desorption curve and the shape of the hysteresis loop [41]. For all the samples, the adsorption and desorption curves do not coincide in the region with higher relative pressure (), and a hysteresis loop forms. The LT-N2GA branch of the coal samples from the study area is steep at saturated vapour pressure (), the desorption curve is similar to the adsorption curve but steep at saturated vapour pressure (), and the adsorption and desorption curves tend to coincide, forming a narrow hysteresis loop between them. This type of curve is similar to the H3 type proposed by the IUPAC and has the characteristics of the H4 type, indicating that the coal samples from the study area contained many nanoscale pores and that the pore types are mainly conical pores, ink bottle pores, and slit pores.

Based on the NLDFT model, LT-N2GA data were calculated to obtain the PSD curve and pore structure parameters of the coal samples in the pore size range of 1.06–78 nm (Table 2). The DFT PV of the coal samples ranges from 0.004 to 0.007 cm3/g (mean 0.005 cm3/g). The DFT SSA ranges from 2.46 to 3.96 m2/g (mean 3.33 m2/g). The average pore diameter is 3.135 nm corresponding to a mesopore. The PSDs of the coal samples from the study area are shown in Figure 7. With increasing pore size, the PV of the stage decreases significantly. When the pore size is larger than 10 nm, the change in PV is small, indicating that the proportion of 10 nm pores is low.

4.4. Analysis of High-Pressure MIP Data

Figure 8 shows the intrusion-extrusion curves of coal samples from the study area. The shape of this curve can reflect the distribution characteristics of the pore throat sizes and the quality of the pore connectivity [15]. Overall, the curve shapes are similar, indicating that the coal samples in the study area have relatively similar structures. The high-pressure MIP curve has the characteristics of a gentle front end and a steep back end. When the mercury inlet pressure is less than 100 MPa, the inflow is slow, and the corresponding pore size mainly increases, indicating that the pores are relatively undeveloped. When the pressure reaches approximately 100 MPa, the mercury intake increases rapidly, reflecting the development of many nanoscale pores in the coal sample. The mercury removal efficiency of all the coal samples is high, which reflects that the pores developed in the coal samples are mostly open pores with good connectivity.

Based on the Washburn equation, the pore size distribution and pore structure parameters of the coal samples with pore diameters ranging from 3 nm to 10 μm were obtained by calculating high-pressure MIP experimental data (Table 2). The PV of the coal samples varies from 0.014 to 0.026 cm3/g (mean 0.019 cm3/g), and the SSA varies from 7.42 to 16.31 m2/g (mean 11.27 m2/g). The PSDs of the coal samples from the study area are shown in Figure 9. The highest mercury intake occurs within the pore size range of less than 30 nm, nearly 80% of the mercury saturation, indicating that the coal samples from the Benxi Formation have mainly nanoscale pores.

4.5. Micro-CT Analysis Data

A micro-CT was used to identify three different components in coal, namely, the pores and fractures, the minerals, and the coal matrix. High-density minerals are shown in white in the corresponding images, low-density pores and fractures are shown in black in the corresponding images, and the coal matrix is shown in grey in the corresponding images (Figure 10) [4244]. The superposition of the pores and fractures, coal matrix, and minerals constitutes the micro-CT three-dimensional pore-fracture structure model of the coal sample (Figure 11). Figures 10 and 11 show that microfractures are mainly developed in all the coal samples, and microfractures are mainly manifested as isolated points, interconnected slender strips, or zigzag curves. Large-scale microfractures constitute the coal reservoir fracture system, which is the main seepage channel controlling the production of CBM. In addition, there are micron-scale pore systems in all the coal samples, which are distributed in almost every area of the coal sample, mainly in irregular ellipsoid shapes and multilateral shapes, and some in relatively concentrated flake shapes, which are conducive to the diffusion of free gas in the coal reservoirs. The fracture system and micron-scale pore system of the coal reservoir together constitute the complex dual-porosity structure of coal [4547]. In general, in three-dimensional space, the microfractures are connected to each other, forming a network structure, and the overall connectivity is good (Figure 11). The total volume, microfracture volume, and SSA of the coal samples were determined via micro-CT scanning and three-dimensional reconstruction. The coal density was calculated with the micro-CT data to compare it with the measured density of the coal samples. Finally, the microfracture volume of the coal samples was calculated to be 0.0031–0.0116 cm3/g. The SSA of the microfractures ranges from 0.0002 to 0.0009 m2/g, among which DJ57-1 has the largest microfracture volume, and DJ57-2 has the smallest volume. As shown in Figure 11, the relationship between the equivalent pore diameter and volume of micron-scale pores and fractures identified via CT scanning is dominated by μm.

5. Discussion

5.1. Full-Scale Pore-Fracture Distribution Characteristics
5.1.1. Principle of Dominant Pore Size Segments

Due to the different pore size ranges measured by different methods, it is difficult to fully characterize pores and fractures by a single method [4852]. According to the test principle and calculation model, micro-CT scanning, high-pressure MIP, LT-N2GA, and LP-CO2GA all have their own advantageous pore size segments, in which the corresponding method can accurately characterize the pore-fracture structure [35, 48, 52]. As shown in Figure 12, rectangles of different widths represent pores and fractures of different scales. According to the analysis of the above test methods, CO2 is an ideal micropore test probe. At 273 K, the molecular kinetic diameter of CO2 is only 0.33 nm, and its saturated vapour pressure is relatively high (3480 KPa) [35, 48]. At a higher relative pressure, the micropores can be filled, the diffusion rate is fast, and the pores in the pore size range of 0.3–1.5 nm can be accurately characterized; this is the advantageous pore size segment of the LP-CO2GA experiment. However, because CO2 molecules do not undergo capillary condensation in meso- or macropores, the Kelvin equation is no longer applicable, and the traditional Barrett–Joyner–Halenda (BJH) model cannot be used to test meso- or macropores [4951]. For the LT-N2GA method, the larger the pore size filled with gas during the test process, the closer the relative pressure () required is to 1. Due to the accuracy of the sensor and control mechanism of the instrument, it is impossible to accurately test macropores with a pore  nm. For pores with pore  nm, N2 molecules must be filled under very low relative pressure (10-7–10-5). At this time, the N2 diffusion rate is extremely slow, and adsorption equilibrium can take a long time to reach or even may never be reached. As a result, it is difficult to test pores with pore  nm via this method. Therefore, the pore size range of 1.5–50 nm is the dominant pore size segment of the LT-N2GA experiment. The maximum theoretical pore size of the high-pressure MIP method ranges from 3.6 nm to 10 μm, covering the distribution of mesopores and macropores. However, due to the compressibility of coal, fractures easily expand under the action of high-pressure mercury; on the other hand, when the mercury inlet pressure is greater than 30 MPa (pore diameter of approximately 50 nm), high-pressure mercury destroys micropores [13, 35, 52]. Therefore, high-pressure MIP experiments are not suitable for characterizing mesoporous and large-scale fractures, but they can more accurately characterize macropores ranging in size from 50 nm to 10 μm, which is the dominant aperture segment of high-pressure MIP experiments. The resolution of a micro-CT scan is related to the size of the coal sample. For the cylindrical sample tested in this work (50 mm high with a diameter of 25 mm), the scanning resolution was 10 μm. Although micro-CT scans are limited by their resolution, they can accurately characterize fractures in aperture segments larger than 10 μm, thus compensating for the inability of high-pressure MIP experiments to accurately characterize fractures.

In summary, this study examined the test data for the optimal aperture range of each gas adsorption method and selected a suitable calculation model. The LP-CO2GA experimental data were calculated based on the NLDFT model, and the effective aperture measurement range is 0.3–1.5 nm. The LT-N2GA experimental data can be calculated based on the NLDFT model, and the effective aperture measurement range is 1.06–78 nm. The overlap between the two models is 1.06–1.50 nm. Due to the characteristics of the test technology and the calculation model itself, the overlap interval was characterized by the results of the LP-CO2GA experiment. Based on the Washburn equation, the effective aperture measurement range of 3 nm to 10 μm was calculated. The overlap between the results of this method and the results of the LT-N2GA experiment is 3.0–50 nm. Due to the matrix compression effect when the pressure in the high-pressure MIP experiment was greater than 30 MPa, the overlap interval was characterized by the results of the LT-N2GA experiment. The resolution of the micro-CT scans is approximately 10 μm, which is close to the upper limit of the experimental results for high-pressure MIP in this paper. To increase the accuracy and precision of the analysis results, the advantages of certain tests in the corresponding intervals were fully utilized to ensure the continuity of the pore-fracture results in the entire aperture range (Figure 12).

5.1.2. Comprehensive Characteristics of Full-Scale Pores and Fractures

According to the joint characterization results of the pore and fracture volume of the coal samples from the study area (Figure 13), the pore and fracture volume distributions are mainly bimodal, reflecting the coexistence of micropores and microfractures, mainly concentrated in the ranges of 0.3–1.5 nm and >100 μm, respectively. The pore and fracture volumes of the coal samples range from 0.057 to 0.075 cm3/g (mean 0.067 cm3/g). Among them, the PVs are mostly associated with the micropores, which provide volumes of 0.042–0.061 cm3/g (mean 0.053 cm3/g). The second-most common types of fractures and pores are microfractures and mesopores, with the volume of microfractures ranging from 0.003 to 0.012 cm3/g and that of mesopores ranging from 0.003 to 0.006 cm3/g. The macropore PV is the smallest, between 0.001 and 0.002 cm3/g. The proportions of pore and fracture volumes at different scales are different (Figure 13 and Table 3). Among them, the pore and fracture volumes in the coal samples are dominated by micropores, with an average proportion of 78.00%, followed by mesopores and microfractures, with average proportions of 6.78% and 13.14%, respectively; additionally, macropores are the least abundant, with an average proportion of 2.08%. The joint characterization results of pore size variation with SSA (Figure 14) reveal that the SSA distribution is mainly unimodal and that the peak SSA value is mainly concentrated in the pore size range of 0.3–1.5 nm. The SSAs of all the coal pores and fractures range from 145.0 to 209.3 m2/g (mean 182.5 m2/g). Among these components, the SSA of micropores makes the largest contribution, ranging from 143.7 to 207.4 m2/g (mean 180.5 m2/g). The mesopore contribution is the second highest, with an SSA ranging from 1.33 to 2.86 m2/g (mean 2.03 m2/g). The contributions of the macropores and microfractures are the smallest. The SSA of the macropores is between 0.01 and 0.03 m2/g. The SSA of the microfractures is between 0.0002 and 0.0009 m2/g. The SSAs of pores and fissures at different scales are different (Figure 14 and Table 3). Among these components, micropores are the dominant pore SSA in the coal samples, accounting for 98.89% of the total pore SSA on average, and the SSA provided by mesopores, macropores, and microfractures can be ignored. Since the coal reservoirs of the Benxi Formation in the Daning–Jixian blocks are buried deep (>2000 m) and the associated vitrinite reflectance indicates that the OM is in the overmature stage, with increasing coal rank, a large amount of gas is generated, resulting in abundant nanoscale pores in the OM and a high proportion of micropores overall. According to SEM observations of the coal samples from the study area, OM is widely distributed in the coal samples, and the presence of OM pores is one of the important reasons for the increase in porosity in the high-maturity stage of this coal. In addition, the contribution of microfractures to the total pore-fracture volume is approximately 13.1%. The main reasons include the following: (1) during the diagenetic evolution of the coal matrix or clay minerals in the Benxi Formation, many shrinkage-induced fractures formed due to geological processes such as mechanical compaction, dehydration, and degassing, which result in a large fracture volume; and (2) the Daning–Jixian block was tectonically located at the eastern margin of the Yishan Slope in the Ordos Basin and the southern end of the Jinxi Fault fold belt. This region has experienced intense periods of regional deformation and microstructural development. Affected by compressive stress, local high-point tensile fractures of folds are widely developed, and abundant secondary microfractures develop, resulting in a high volume proportion of microfractures [1, 2]. The PV and SSA of the coal samples are mainly attributed to micropore characteristics. Studies have shown that micropores provide many adsorption sites and occurrence spaces for CBM and are also important channels for the initial migration and diffusion of CBM after desorption; therefore, further study and analysis of the development and structural characteristics of CBM are necessary [13].

5.2. Contribution of the Pore-Fracture Structure to Permeability

Coal is a complex porous medium. According to the theory of dual-porosity structures, a porous medium includes both pores and fractures [53]. The pore-fracture system includes not only the occurrence space of CBM but also the migration channels of CBM. The characteristics of coal reservoirs, such as adsorbability, gas-bearing capacity, and permeability, are related to the development degree and structural characteristics of pores and fractures [5355]. In this paper, a formula similar to that of Sinn et al. is used to calculate the deep coal reservoir permeability [56]. An ideal capillary bundle is assumed, and laminar flow through a cylindrical pipe is described by combining Darcy’s law and the Hagen–Poiseuille equation to express permeability. Based on LT-N2GA, the LP-CO2GA method, and the high-pressure MIP method, the pore characteristics (PV, SSA, and pore diameter) of the coal samples can be quantitatively characterized; thus, the porosity corresponding to each pore can be calculated. Combined with the pore radius (), the absolute permeability of a single pore can be further calculated () [56]: where is the porosity corresponding to the ith pore and is the radius of the th pore. is the absolute permeability of the th pore. is the calculated total pore permeability of the coal sample.

Based on micro-CT scanning and three-dimensional reconstruction technology, we quantitatively characterize the volume of microfractures and the total volume of coal samples after three-dimensional reconstruction, calculate the pore (fracture) clearance () corresponding to each microfracture, and then combine the microfracture width () to further calculate the absolute permeability () of each microfracture [56]: where is the microfracture degree of the th microfracture and is the tortuosity of microfractures. In this study, was used for calculation [56, 57], and is the absolute permeability of the th microfracture. is the calculated total fracture permeability of the coal sample.

With LP-CO2GA, LT-N2GA, high-pressure MIP, and micro-CT scanning data, the pore and fracture volume of the coal samples can be quantitatively characterized (Table 3), and the pore (fracture) clearance corresponding to each pore and each microfracture can be calculated in combination with the pore radius and microfracture width. The total permeability (sum of the pore permeability and fracture permeability) of the Benxi Formation ranges from 5.77 to 28.22 mD. There are three reasons for the slight difference between the calculated total permeability and the measured steady-state Farcke permeability (Figure 15). (1) When the steady-state Farcke permeability in a cylindrical sample (50 mm high with a diameter of 25 mm) is tested under a certain confining pressure (5 MPa), the microfractures may be partially or completely closed under the stress, which reduces the permeability of the coal samples. 2) Although the tortuosity () of microfractures is considered in the permeability calculation using the Hagen–Poiseuille equation (Equation (2)), the actual pore-fracture structure in coal samples is extremely complex, resulting in a relatively small measured permeability. (3) Through X-ray scanning and imaging, micro-CT technology can be used to identify not only connected microfractures but also isolated and closed microfractures [24, 47, 58]. Because these experimental procedures cannot be applied, the calculated fracture permeability is too high.

Figure 16 shows the contribution distribution of the pore (fracture) permeability to the calculated total permeability at the different pore segments. The permeability of the Benxi Formation coal samples is mainly controlled by microfractures with pore μm, accounting for approximately 95% of the total permeability. The contribution of pores with μm to the permeability of coal is relatively low, at only approximately 2% (Figure 16). The main reasons are as follows: the extent of μm is large, and more relatively independent pores and smaller microfractures can be connected spatially, thus forming a relatively developed microfracture network (Figures 10 and 11), which greatly improves the seepage capacity of coal samples [42, 43]. Although there are many large pores and small microfractures in the coal samples, these pores and microfractures are relatively spatially independent and have poor connectivity with each other (Figure 11); thus, they contribute little to the permeability of coal reservoirs.

5.3. Effect of the Pore-Fracture Structure on Gas Migration

The adsorption and desorption characteristics of CBM in nanoscale pores in microporous and fractured structures are studied. The microscale pores and connected pores in the structure are the main channels for CBM migration and play a decisive role in the flow migration of CBM and the permeability of coal reservoirs [55, 59]. In this paper, by combining the experimental methods of LP-CO2GA, LT-N2GA, high-pressure MIP, and micro-CT scanning, it was found that the pore-fracture structure of coal reservoirs is heterogeneous and spans multiple scale characteristics. The wide distribution of fractures from the nanoscale to the microscale results in the formation of a complex pore-fracture network (Figure 11). Therefore, based on a variety of pore characterization methods for the full-scale quantitative characterization of coal sample pores and fractures, there are obvious differences in the distributions of pores and fractures at different scales. The coal is mainly composed of 0.3–1.5 nm and >100 μm particles, in which the volume proportions of micropores (<2 nm), mesopores (2–50 nm), macropores (50 nm to 10 μm), and microfractures (>10 μm) are 78.00%, 6.78%, 2.08%, and 13.14%, respectively (Figure 13). This multiscale approach is highly important for the comprehensive evaluation of deep coal reservoir characteristics.

CBM mainly exists in the pore-fracture system of coal samples in adsorbed and free states. The adsorbed gas is mainly adsorbed on the surface of the coal matrix and clay mineral particles, and the free gas is mainly distributed in the large pores and fractures of the coal. Yang et al. and Li et al. found that the deep coal seam in the eastern margin of the Ordos Basin is in a state of supersaturation with respect to gas content overall and that the proportion of free gas is 17-34%, with an average of approximately 20%, which is characterized by “high gas content, high gas saturation, and rich free gas” [6063]. Tang et al. and Ma et al. studied the relationship between the pore structure at different scales and the gas content characteristics of deep coal reservoirs in the Daning–Jixian block. PV and SSA at different scales have different mechanisms of action on the CBM content in different phase states (Figure 17) [62, 64]. The adsorbed gas volume is closely related to the PV and SSA of the micropores, and the two variables have a strong linear relationship ( and , respectively; Figure 17). The results show that with the development of micropores, the adsorbed gas volume increases gradually, and micropores can provide many adsorption sites, which are conducive to gas adsorption and provide a place for the adsorption and occurrence of deep CBM. However, the free gas volume is closely related to the PV and SSA of microfractures, and these two variables are linearly positively correlated ( and , respectively; Figure 17), indicating that free gas mainly occurs in microfractures. The more developed the microfractures are, the more gas storage space is provided, and the greater the free gas volume is. The adsorbed gas in deep coal seams mainly occurs in micropores, which are mainly affected by the large adsorption SSA, while the free gas mainly occurs in microfractures, which are mainly affected by the gas volume.

According to the distribution characteristics of the Benxi Formation coal samples from the study area, the micropore PV of the coal samples ranges from 0.042 to 0.061 cm3/g (mean 0.053 cm3/g) (Table 3), which is slightly lower than that of Taiyuan Formation coal samples from the Qinshui Basin (0.046–0.062 cm3/g) [15] and much greater than that of the shale of the Longmaxi Formation in the Sichuan Basin (mean 0.0104 cm3/g) [64, 65]. The micropore SSA of the Benxi Formation coal samples ranges from 143.7 to 207.4 m2/g (mean 180.5 m2/g), which is similar to that of the coal samples from the Taiyuan Formation in the Qinshui Basin (139.6–209.129 cm3/g) [15]. The Benxi Formation micropore SSA is dozens of times greater than that of the Longmaxi Formation shales (7.57–30.59 m2/g) in the southeastern Sichuan Basin and the Shanxi Formation shales (0.655–11.42 m2/g) in the eastern Ordos Basin [64, 65]. Because micropores can provide a large SSA, they can provide adsorption sites for methane molecular adsorption and indirectly control the methane adsorption capacity in deep coal reservoirs. Moreover, micropores and microfractures, which are the main reservoir spaces available for free gas occurrence, have relatively low SSAs but play decisive roles in determining the permeability of coal reservoirs and CBM production [66, 67]. With the continuous desorption of gas in coal and the decrease in coal reservoir pressure, the amount of gas in pores and fractures in coal seams changes from increasing to decreasing, and multistage diffusion and penetration occur through interconnected nanoscale pores and microscale fractures [68]. Moreover, the form and scale of gas migration in pores and fractures at different scales are different, and interconnected pores and microfractures are effective channels for fluid transport (Figure 18). LP-CO2GA, LT-N2GA, high-pressure MIP, and micro-CT scanning are used to characterize the full-scale pore-fracture structure of coal reservoirs. The advantageous pore distribution intervals determined by the respective experiments are utilized to quantitatively characterize the micropores, mesopores, macropores, and microfractures. Among them, the development characteristics of micropores and microfractures control the adsorption capacity and development potential of deep CBM, respectively. Research on the development characteristics of nanoscale pores and microscale fractures plays a decisive role in the enrichment and exploitation potential of deep CBM, and further research on the migration mechanism, scale, and controlling factors of deep CBM is needed.

6. Conclusions

The results of this study on the characteristics of pore structure in a deep coal reservoir and its influence on CBM occurrence, with a case study of the Daning–Jixian block, support the following conclusions: (1)The coal samples from the Benxi Formation are mainly composed of OM pores, inorganic pores, and microfractures. The OM pores have various shapes, mainly round, elliptical, and wedge-shaped. The inorganic pores can be divided into intraparticle and interparticle pores. The diameters of the intraparticle pores are small, between 20 and 50 nm, while the diameters of the interparticle pores are larger, generally greater than 300 nm. The microfractures are mainly distributed within the margins of the coal matrix and inside the clay minerals and are elongated or zigzag-shaped(2)According to the full-scale quantitative characterization of the pore and fracture volume distributions with respect to pore size, the PSD is mainly U-shaped, exhibiting a bimodal state of micropores and microfractures, which are mainly concentrated in the ranges of 0.3–1.5 nm and >100 μm, respectively. The pore-fracture volume is dominated by micropores (average volume proportion of 78.0%), followed by mesopores and microfractures (average volume proportions of 6.8% and 13.1%, respectively). Macropores were the least common (mean 2.1%)(3)The microfractures are interconnected in three-dimensional space, forming a network structure with strong connectivity. The permeability is mainly contributed by microfractures with pore μm, accounting for approximately 95% of the total permeability. These findings indicate that microfractures are the main channels for deep CBM migration and play a decisive role in the flow of deep CBM(4)By combining experimental methods such as LP-CO2GA, LT-N2GA, high-pressure MIP, and micro-CT scanning, we determine that the pore structure of deep coal reservoirs is heterogeneous and spans multiple scales. The development characteristics of nanoscale pores and micrometre fractures control the enrichment capacity and development potential of deep CBM, respectively

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Funding

This research was jointly supported by Grant No. 2022KT1401 and Grant No. 2021DJ2302.

Acknowledgments

We thank PetroChina Coalbed Methane Company, Ltd., for collecting the samples.