1 Introduction

Pores in coal are the main spaces for coalbed methane (CBM) storage, while fractures are the main channel for gas and water flow (Liu et al. 2015, 2018a). The pore structure, e.g., porosity, scale, morphology and connectivity, plays a key role in CBM production (Liu et al. 2015, 2018a; Lai et al. 2016, 2018). The pores in coal are generally classified into four categories, i.e., primary pores, epigenetic pores, coalification-related pores and mineral-related pores (Chen et al. 2015; Lai et al. 2016, 2018; Liu et al. 2017). The genesis and morphology of the pores are closely related to the microlithotype, metamorphism and deformation of the coal (Alexeev et al. 2007; Zhang et al. 2010, 2018; Zhu et al. 2016; Li et al. 2018). The pores in coal directly impact gas generation, gas accumulation and coal permeability, and control the adsorption, diffusion and seepage of CH4 in coal (Alexeev et al. 2007; Zhang et al. 2010; Zhu et al. 2016; Du et al. 2019). Fractures are the main channel for fluid migration and output and are the key factor for coal permeability (Yao et al. 2010a; Liu et al. 2016). Fractures can be divided into macroscopic fractures and microscopic fractures (Chen et al. 2015; Lai et al. 2017; Zhang et al. 2017; Zhou et al. 2018). Microscopic fractures can be further divided into endogenous fractures and exogenous fractures (Chen et al. 2015; Zhang et al. 2017; Zhou et al. 2018). Fractures are generally believed to form from internal forces and external forces (Yao et al. 2010b; Li et al. 2014; Liu et al. 2016). Tectonic evolution, metamorphism and microlithotype play a key role in the formation of fractures (Yao et al. 2010a; Li et al. 2014; Liu et al. 2016; Ai et al. 2018).

Research on pores and fractures is focused on their morphology, genetic type, scale, syntagmatic relation, connectivity, and their relationships with coal porosity and permeability, fluid occurrence, migration and output (Liu et al. 2016, 2017; Li et al. 2017, 2019). Characteristics of the coal structure parallel and perpendicular to the bedding plane are rarely considered. For vertical CBM wells, the fluid flow in coal can be regarded as radial fluid flow; thus, fluid migration perpendicular to the bedding plane can be ignored, and the coal structure, connectivity and permeability in this orientation can also be disregarded (Clarkson and Bustin 1999; Wei et al. 2007; Aminian and Ameri 2009). However, during horizontal CBM well production, fluid flow is mainly sphere-shaped flow in the horizontal section; thus, fluid migration perpendicular to the bedding plane cannot be ignored, and the coal structure, connectivity and permeability in this orientation are important to CBM production (Keim et al. 2011; Wen et al. 2011; Fang et al. 2019). Therefore, the investigation of the growth characteristics and connectivity of pores and fractures in coal parallel and perpendicular to the bedding plane is significant to horizontal CBM well production and can provide a better understanding of fluid migration in coal.

Pores and fractures in coal are mainly studied by traditional methods, e.g., optical microscopy, electron microscopy, mercury intrusion porosimetry and nitrogen/CO2 adsorption tests (Bolin 2014; Okolo et al. 2015). These traditional methods have limitations, e.g., difficulty in studying pore-fracture connectivity and the destruction of pores and fractures during sample preparation and testing (Okolo et al. 2015; Liu et al. 2018a). Nuclear magnetic resonance (NMR) and X-ray computed tomography (CT) scanning are becoming important technologies to study pores and fractures in coal due to their measurement speed and nondestructive testing (Mathews et al. 2011; Yao and Liu 2012; Smolná et al. 2013; Zhao et al. 2017). In this paper, an integrated approach combining low-field nuclear magnetic resonance (LF-NMR), X-ray CT, and field emission scanning electron microscopy (FESEM) was employed to describe the coal porosity, pore size distribution, genetic types, and connectivity parallel and perpendicular to the bedding plane, and the reasons for the differences in connectivity between the two orientations were further explained. This study aims to provide a better understanding of fluid migration and output during horizontal CBM well production.

2 Materials and methodology

2.1 Sample

The two coal samples used in this study were systematically collected from coal seam #3 of the Shanxi Formation in the southern Qinshui Basin in China. Anthracite was collected from the Sihe Mine, and semi-anthracite coal was collected from the Yuwu Mine (Table 1). The coal samples were collected in large blocks (approximately 300 × 300 × 300 mm). The coal pillars used in the study were drilled parallel and perpendicular to the bedding plane for each coal sample (Fig. 1). Coal pillars with a diameter of approximately 25 mm and a length of 40–60 mm were drilled by a mechanical drill, and their ends were polished with a grinder (Fig. 1). The names, key properties and sizes of these coal pillars are shown in Table 1.

Table 1 Physical properties of the coal samples used
Fig. 1
figure 1

Modified from Liu et al. (2018b)

Sketch map of coal sample preparation.

2.2 LF-NMR

2.2.1 LF-NMR analyzer

The NMR experiment was performed with a MicroMR12-025 V LF-NMR analyzer made by Suzhou Niumag Analytical Instrument Corporation, China. The resonance frequency of the LF-NMR analyzer is 11.897 MHz, and the magnet temperature is 1.99–32.01 °C.

2.2.2 LF-NMR experimental process

  1. A.

    The coal pillars were first vacuum dried at 50 °C for 24 h and then immediately saturated with distilled water for more than 24 h.

  2. B.

    LF-NMR experiments on 100% water-saturated coal pillars were carried out to obtain T2 (the transverse relaxation time, ms) spectra under 100% water-saturated conditions (Sw).

  3. C.

    Then, the coal pillars were centrifuged at a high speed using a centrifugal pump at 200 Psi for 1.5 h to achieve irreducible water conditions (Sir).

  4. D.

    LF-NMR experiments on coal pillars at Sir were carried out to obtain T2 spectra at Sir.

2.2.3 T2 spectra analysis

A. Effective porosity


The NMR signals of coal pillars at Sw can be measured and scaled using a standard sample. Then, NMR signals can be converted to porosity using Eq. 1.

$$\varphi = \varPhi \cdot \frac{s}{S} \cdot \frac{{N_{s} }}{{n_{\text{s}} }} \cdot 10^{{\frac{1}{20}(RG1 - rg1)}} \cdot 2^{(RG2 - rg2)}$$
(1)

where Φ, RG1, RG2 and S are the porosity, simulation gain, signal strength and NMR signal value of the standard sample, respectively; φ, ns, rg1, rg2, and s are the porosity, cumulative collection time, simulation gain, signal strength and NMR signal value of the coal pillars, respectively.

The effective porosity (φpr, also called producible porosity, %) is obtained by the cumulative T2 spectra method proposed by Yao et al. (Yao et al. 2014; Zhao et al. 2017). The maximum peak amplitudes of cumulative T2 spectra at Sw and Sir are equal to cumulative porosity at Sw (φw, %) and Sir (φir, also called irreducible porosity, %), respectively. φpr is the difference between the φw and φir.

B. Pore size distribution

The “T2 cutoff value” (T2C is the threshold for φpr and φir, ms) is used to calculate the pore size distribution of the coal samples. Assuming that the pore radius (r) corresponding to T2C is constant for all coal samples, the pore radius (rci) corresponding to i relaxation time (T2i) can be expressed as Eq. 2 (Yao and Liu 2012; Yao et al. 2014).

$$r_{ci} = r\frac{{T_{2i} }}{{T_{2C} }}$$
(2)

In this study, the pore radius corresponding to the T2C was chosen to be 0.1 μm.

2.3 Three-dimensional (3D) coal structure modeling

2.3.1 X-ray CT scanning

X-ray CT scanning was performed by the X-ray CT scanner Xradia 520 Versa produced by the Carl Zeiss Foundation Group, Germany. The X-ray CT scanning samples were small coal pillars of approximately 2 mm in diameter and 2 mm in height (Fig. 2). The samples were drilled using a mechanical sampler.

Fig. 2
figure 2

Photographs of coal samples for X-ray CT

The X-ray CT scanning area was 2 mm in diameter and 1 mm in height. The total of 3600 scans were performed at a spatial resolution of 1 μm. X-ray CT scanning provides structural information for the development of a 3D pore-fracture network model of the coal samples at a cubic scale of approximately 400 μm on each side. The process of establishing the 3D pore-fracture network model using the X-ray CT scanner includes several steps, and threshold selection is the key to identifying pores and fractures during this process. The distribution intervals of the thresholds for minerals, organic matter and pore fractures were 3000 Hounsfield units (HU), 1000–1800 HU and < 600 HU, respectively (Yao et al. 2010a).

2.3.2 3D pore-fracture network modeling

A 3D pore-fracture network model is established via Avizo 9, which is professional software used to produce a 3D digital core using 3D X-ray CT scanning images. The main steps include skeletonization of the pore-fracture space, description of the pore-fracture system, segmentation of the pore-fracture space, establishment of a new 3D pore-fracture network model and analysis of the spatial structure parameters. The process of establishing the 3D pore-fracture network model is beyond the scope of this paper, and the details can be referenced elsewhere (Hughes and Blunt 2001; Knackstedt et al. 2005; Ma et al. 2014a, b).

2.4 Scanning electron microscopy analysis

Pores, fractures and minerals were investigated using an FESEM SUPRA 55 produced by the Carl Zeiss Foundation Group, Germany, with amplification from 103 to 105. Bulk coal samples were analyzed by FESEM. The coal samples were polished into small lumps approximately 10 mm across and 2–3 mm high using a polishing and burnishing machine. Then, the small lumps of coal were polished using a cross-section polisher. Coal is known as a nonconducting substance. Therefore, to achieve better experimental results, coal samples were usually spluttered with a thin gold coating. The pores were observed using the high-vacuum mode.

3 Results

3.1 Porosity of the coal samples

According to the NMR results, the development and connectivity of the pores in the Sihe coal samples are better than those in the Yuwu coal samples. The φN, φw, φir and φpr of the Sihe coal samples are higher than those of the Yuwu coal samples, where the φN of the Sihe coal samples is approximately 2 times that of the Yuwu coal samples (Table 2). The φN and φir of the coal samples parallel to the bedding plane are lower than those perpendicular to the bedding plane, whereas the φw and φpr of the coal samples parallel to the bedding plane are higher than those perpendicular to the bedding plane (Table 2). This finding indicates that the coal samples have better pore connectivity parallel to the bedding plane.

Table 2 Porosity, effective porosity and other parameters based on NMR and X-ray CT

The X-ray CT results also show that the pore development in the Sihe coal samples is better than that of the Yuwu coal samples due to the higher pore volumes and porosities (Table 2). The porosity of coal SH-B is lower than that of coal SH-V, and the Yuwu coal samples are the opposite.

3.2 Pore size distribution of the coal samples

This study used a pore structure classification system proposed by the International Union of Pure and Applied Chemistry (IUPAC), classifying the pores into micropores (< 2 nm in width), mesopores (2–50 nm in width) and macropores (> 50 nm in width) (Moore 2012; Liu et al. 2017).

3.2.1 Pore size distribution by LF-NMR

The Sihe and Yuwu coal samples have similar pore size distributions. Mesopores (48.21%–76.61%) are most abundant, followed by macropores (23.37%–49.38%) (Fig. 3). The pores widths of the coal samples are mainly in the range of 2–500 nm, and the pore volume percent decreases with increasing pore width (Fig. 3). The pore size distribution parallel to the bedding plane is also similar to that perpendicular to the bedding plane. Mesopores are the most abundant in both orientations, followed by macropores of 50–500 nm in width, with the smallest volume percent of macropores > 1 μm (Fig. 3). In particular, 76.61% of the pores in the coal YW-V are mesopores, and the volume percent of macropores decreases rapidly (Fig. 3). The volume percent of mesopores and macropores > 1 μm in the coal samples perpendicular to the bedding plane are slightly higher than those in the samples parallel to the bedding plane, indicating that mesopores and macropores > 1 μm are more abundant perpendicular to the bedding plane.

Fig. 3
figure 3

Pore size distribution determined by NMR

3.2.2 Pore size distribution by X-ray CT

The X-ray CT results show that the number and volume of macropores > 1 μm in the Sihe coal samples are higher than those in the Yuwu coal samples (Fig. 4). The number of macropores > 1 μm in the coal samples perpendicular to the bedding plane are larger than that in the coal samples parallel to the bedding plane. The volumes of macropores 1 μm–6 μm and > 50 μm in width in coal SH-V are significantly greater than those in coal SH-B, while the volumes of macropores 6 μm–30 μm in width in coal SH-V are slightly lower (Fig. 4). The volumes of macropores 1 μm–4 μm in width in coal YW-V are significantly greater than those in coal YW-B (Fig. 4). The volumes of macropores > 4 μm in coal YW-V are lower than those in coal YW-B, and with increasing pore width, the differences between the pore volumes of coal YW-B and coal YW-V become more significant (Fig. 4). Therefore, pores with smaller widths are more abundant in the coal samples perpendicular to the bedding plane, while microscopic fractures and pores with larger width are more abundant in the coal samples parallel to the bedding plane.

Fig. 4
figure 4

Pore size distribution determined by X-ray CT

3.3 Connectivity of the pores and fractures

3.3.1 Pores and fractures identified by T2 spectra

The mesopores, macropores and microscopic fractures can be identified based on the T2 spectra at Sw (Yao et al. 2010a; Yao and Liu 2012). Coal SH-B and coal SH-V have similar T2 spectra at Sw. Their T2 spectra at Sw have three peaks consisting of two large peaks and one small peak (Fig. 5). The two large peaks located at 0.01–0.2 ms and 0.2–5.0 ms are the peaks of the most abundant pores, which represent pores with width ranges of 2–15 nm and 15–380 nm in coal SH-B, respectively, and represent pores with width ranges of 2–12 nm and 12–230 nm in coal SH-V, respectively (Fig. 5). The small peak is located at 20–100 ms (Fig. 5), which is the peak of the macropores > 1 μm, indicating that macropores > 1 μm are not abundant. In addition, the T2 spectra at Sw of the Sihe coal samples lack peaks of microscopic fractures (Fig. 5), indicating that microscopic fractures are not abundant.

Fig. 5
figure 5

NMR T2 distribution

Similar to the Sihe coal samples, coal YW-B and coal YW-V also have similar T2 spectra at Sw. Their T2 spectra at Sw have one large peak and three small peaks (Fig. 5). The large peak is located at 0.02–1.0 ms and represents pores with width ranges of 10–300 nm in coal YW-B and 3.5–160 nm in coal YW-V (Fig. 5). The three small peaks are located at 2–25 ms, 25–500 ms and 500–1000 ms (Fig. 5), which are the peaks of macropores, microfissures and microfissures and cleats, respectively. The results show that microscopic fractures are abundant in the Yuwu coal samples, which are related to the tectonic deformation of coal in the Yuwu Mine.

3.3.2 Connectivity of the pores and fractures by LF-NMR

By comparing the T2 spectra at Sw and Sir, the connectivity of the coal samples can be further analyzed (Yao et al. 2010a; Yao and Liu 2012). The peak of 0.01–0.2 ms in the Sihe coal samples has good continuity with the peak of 0.2–5.0 ms, indicating that pores with width ranges of 50–500 nm have good connectivity (Fig. 5). However, the peak of 0.2–5.0 ms is almost discontinuous with the peak of 20–100 ms due to the low number of pores > 1 μm (Fig. 3). Coal SH-B has a different peak of 0.01–0.2 ms at Sir for coal SH-V. The T2 spectra at Sir for SH-B is lower than those at Sw because of the good connectivity of the mesopores 2–15 nm allowing the water to be centrifuged. The T2 spectra at Sir for coal SH-B are significantly lower than those at Sw in the range of 0.01–0.1 ms and higher in the range of 0.1–0.2 ms (Fig. 5), because pores with widths of 2–12 nm have good connectivity, while the connectivity of pores with widths of 2–230 nm is poor. Under the centrifugal effects, the saturated water migrates from pores of 2–12 nm to larger-scale pores and accumulates in pores with widths of approximately 12 nm. Almost no change occurs between the T2 spectra at Sir and those at Sw for coal SH-V (Fig. 5). The connectivity of pores with widths of 12–230 nm is extremely poor, resulting in the connectivity and effective porosity of coal SH-V being lower than those of coal SH-B. The T2 spectra at Sir for coal SH-B are slightly higher than those at Sw in the range of 0.2–5.0 ms because of the low connectivity between pores of 15–380 nm and larger-scale pores. The migration of saturated water to larger-scale pores is difficult under the centrifugal effects. In addition, the small peak of macropores > 1 μm disappears after centrifugation, which is most obvious in coal SH-B (Fig. 5), indicating that macropores > 1 μm have a certain connectivity even though they are not abundant.

The continuity between the peaks of 0.02–1.0 ms and 2.0–25 ms, 2.0–25 ms and 25–500 ms of coal YW-B are weak, while that of coal YW-V is almost discontinuous (Fig. 5). This result shows that the connectivity between the pores of the Yuwu coal samples is poor, which also causes the characteristics of the T2 spectra at Sw and Sir. The T2 spectra at Sir for coal YW-B are significantly lower than those at Sw in the range of 0.02–1.0 ms due to the good connectivity of pores with widths of 10–300 nm. The T2 spectra at Sir for coal YW-B in the range of 2.0–25 ms and 25–500 ms move to the smaller pores, indicating that saturated water is centrifuged from pores corresponding to the peak of 0.02–1.0 ms and accumulates in pores corresponding to 2.0–6.0 ms. Macropores > 1 μm and microfissures in coal YW-B have connectivity, and the saturated water in pores corresponding to the peaks of 2–25 ms and 25–500 ms is centrifuged. The T2 spectra at Sir for coal YW-B in the range of 25–500 ms almost disappeared, indicating that microfissures and cleats have good connectivity. In summary, pores and fractures in coal YW-B have the greatest connectivity among the coal samples. The T2 spectra at Sir for coal YW-V in the range of 25–500 ms almost disappeared, while the other 3 peaks disappeared compared to the T2 spectra at Sw. This result indicates that the connectivity of coal YW-V is extremely poor, and only the microfissures and cleats had connectivity.

3.3.3 Connectivity of the pores and fractures by X-ray CT

NMR has a low resolution of macropores > 1 μm. Therefore, the connectivity of pores and fractures at the micron-scale is further studied based on the X-ray CT results. Pore coordination numbers based on the X-ray CT results can quantitatively express pore connectivity (Lindquist et al. 2000; Sok et al. 2002). Figure 6 shows that the pore connectivity of coal SH-B is significantly better than that of coal SH-V, and the number of pores in each coordination number is larger than that in coal SH-V. The number of pores with a coordination number of 0 in coal YW-B and coal YW-V is basically the same, whereas the number of pores with a coordination number > 1 in coal YW-B is substantially larger than that in coal YW-V. Moreover, the coordination number of interconnected pores in coal SH-B is primarily 1–2, followed by 3–5, while the coordination number of pores in coal YW-B is slightly higher, which is mainly 5–6, followed by 3–4 (Fig. 6). However, there are almost no interconnected pores in coal SH-V and coal YW-V according to the coordination number (Fig. 6), indicating that the coal connectivity is extremely poor perpendicular to the bedding plane.

Fig. 6
figure 6

Coordination number of pores and interconnected pores of the coal samples based on X-ray CT. Notes: a, c Coordination number of all pores in Sihe and Yuwu coal samples, respectively; b, d coordination number of interconnected pores in Sihe and Yuwu coal samples, respectively

4 Discussion

4.1 Comparison of pore-fracture types

4.1.1 Morphology of pores and fractures identified by FESEM and X-ray CT images

Comparing the FESEM photographs (Fig. 7) and the X-ray CT images (Fig. 8), three main categories of pores and microscopic fractures are present in the coal samples, i.e., pores in the matrix (Figs. 7a, b; 8b, e), mineral-related pores (Figs. 7b, c; 8a) and microscopic fractures (Fig. 8c, d).

Fig. 7
figure 7

Pores and fractures in coal samples based on FESEM. Notes: a Secondary gas pores in organic matter, Yuwu Mine; b secondary gas pores in organic and mineral-related pores in minerals, Yuwu Mine; c mineral-related pores in minerals and shrinkage-induced pores around minerals, Sihe Mine; and d shrinkage-induced pores around minerals, Sihe Mine

Fig. 8
figure 8

3D digital model of coal samples based on X-ray CT. Notes: In the cubes, pores and fractures are carmine, minerals are blue, and organic matter is gray. In the X-ray CT images, pores and fractures are black, minerals are white and organic matter is light gray

A. Pores in the matrix

The pores in the matrix develop in the coal organic matter and merge to form clusters with poor connectivity. These pores generally form irregular ellipses, and the long axes of the elliptical pores are directionally oriented due to paleotectonic stress (Figs. 7a, b; 8b, e). The pore widths of these pores are usually < 5 μm and mainly 50–500 nm (Figs. 7a, b; 8b, e). These pores correspond to the peaks located at 0.2–5.0 ms (Sihe Mine) and 0.02–1.0 ms (Yuwu Mine) in the T2 spectra at Sw and are major contributors to pores 1–5 μm in width according to X-ray CT. Based on the morphology, the pores in the matrix are mainly secondary gas pores.

B. Mineral-related pores

Mineral-related pores are formed around or in minerals and merge to form clusters with good connectivity (Figs. 7b, c; 8a). Some of these pores connect to each other to form a larger “hole” (Fig. 7c). Mineral-related pores do not have a fixed shape but are usually in the form of irregular elliptic and linear shapes and are also directionally oriented by the long axis (Figs. 7b; 8a). The pore size distribution of mineral-related pores is wide from mesopores to macropores. The pore widths of these pores are usually < 20 μm and mainly 50–500 nm and 1–10 μm (Figs. 7b; 8a). Mineral-related pores are abundant in the coal samples, correspond to the peaks located at 0.01–0.2 ms, 0.2–5.0 ms (Sihe Mine) and 0.02–1.0 ms (Yuwu Mine) in the T2 spectra at Sw and are major contributors to pores 1–20 μm in width according to X-ray CT. Based on the morphology and distribution area, these pores are mainly shrinkage-induced pores around minerals and pores in minerals, with more abundant shrinkage-induced pores.

C. Microscopic fractures

The microscopic fractures in the coal samples are mainly macropore-scale, and some of the fractures are mesopore-scale. They are “pores” > 50 μm according to the T2 spectra and X-ray CT. Microscopic fractures with mesopores and macropores 50–500 nm in width are mainly linear shrinkage-induced pores (Fig. 7d). Microscopic fractures with widths > 50 μm can be both cleats (C in Fig. 8) and microfissures (Fig. 8d). Cleats are typical endogenous fractures, and they are formed by the combined compaction, dehydration and contraction of devolatilization of the gelification materials during coalification (Liu et al. 2015). There are two groups of predominant cleat directions that are nearly orthotropic, e.g., face cleats and butt cleats. The cleats have good connectivity, and their widths are between a dozen micrometers and dozens of micrometers. The cleats are mostly closed or filled by minerals, which reduces the connectivity (Fig. 8c). The microfissures are relatively narrow (50–10 μm in width) and generally short (less than 1 mm in length), curved and dense without directionality (Fig. 8d). Although the microfissures have good connectivity, they are also filled by minerals, which reduces the connectivity and is the main reason for the lack of microscopic fracture peaks in the T2 spectra at Sw in the Sihe coal samples (Fig. 8d).

4.1.2 Influence of pore-fracture types on connectivity

The 3D digital model and pore model show that coal SH-B is dominated by pores with few microscopic fractures, while coal SH-V has a large number of pores and microscopic fractures (Figs. 8; 9a, b). The abundant microscopic fractures result in a large volume of macropores > 50 μm, as shown in Fig. 4. Opposite to the Sihe coal samples, the pores and microscopic fractures in coal YW-B are more abundant than those in coal YW-V. The interconnected pore models show that coal SH-B and coal YW-B have connectivity, while the coal connectivity is extremely poor perpendicular to the bedding plane (Fig. 9c). That is identical with coordination number of interconnected pores shown in Fig. 6. Coal SH-B and coal YW-B have significantly different connectivity modes. Coal SH-B samples are mainly connected by pores, while microscopic fractures are not abundant. Therefore, coal SH-B has mainly “pore connectivity,” where the throats are mainly pores, with single connecting channels and low coordination numbers. The coal YW-B are connected by microfissures, while the pores are not abundant in coal YW-B, and the connected area in coal YW-B is actually a microfissure with pores around it. Therefore, coal YW-B has “microfissure connectivity,” and the microfissures are the main channels connecting the pores, with more complicated connecting path.

Fig. 9
figure 9

Pore model and interconnected pore model of coal samples based on X-ray CT. Notes: a Pore models; b skeleton models of pores and throats; and c interconnected pore models. In the skeleton models and interconnected pore models, the pores are red, and the throats are green

4.2 Comparing with the throats of coal samples

The connectivity of coal samples indicates that throats are the key of coal connectivity, and the main reasons for the difference in the coal sample connectivity parallel and perpendicular to the bedding plane.

Figures 6b and 10 show that the number, width and length of throats in coal SH-B are larger than those in coal SH-V, which is the main reason for the better connectivity in coal SH-B. The throat width of the Sihe coal samples is mainly 2–3 μm and is consistent with the pore width (Fig. 10), which further confirms that the throats of the Sihe coal samples are mainly pores. However, the 3D digital models (Fig. 8) and interconnected pore models (Fig. 9) show that the connectivity between cleats and pores, microscopic fractures and pores in coal SH-B is weak, and the throats in coal SH-B are mainly shrinkage-induced pores around minerals. Different from the Sihe coal samples, the number, width and length of throats in YW-B are smaller than those in coal YW-V (Figs. 6b, 10). The throats in coal YW-B are mainly microfissures, and the throat connectivity is determined by the connectivity between the microfissures and pores. However, the microfissures are not abundant in coal YW-V, and the throats in coal YW-V are mainly pores, resulting in the low coordination number of pores (0–1) (Fig. 6).

Fig. 10
figure 10

Throat widths and lengths of the coal samples based on X-ray CT

In conclusion, the pores are connected by throats, which determine the connectivity of high-rank coal. The contribution of microfissures to coal connectivity is superior to that of pores. Generally, microfissures are the main throats, which determine the coal connectivity. Then, shrinkage-induced pores become the decisive factor of the pore connectivity when microfissures are not abundant.

4.3 Differences in spatial distribution of pores and fractures

4.3.1 Spatial distribution of pores on Coal SH-B

The pores, microscopic fractures and minerals are directionally oriented due to the influence of paleotectonic stress, forming a strip-like distribution in the coal samples (Figs. 7, 8). The orientations of the pores and minerals in the coal samples perpendicular to the bedding plane are significantly stronger than those of the samples parallel to the bedding plane (Fig. 8), indicating that the influence of the maximum and minimum horizontal principal stresses on the pore and mineral distributions is more significant than the vertical principal stress, causing the long axis of the pore to extend along the minimum horizontal principal stress. Due to the strip-like distribution of pores in coal SH-B, the pore connectivity is continuous parallel to the bedding plane (Fig. 11). Therefore, the pore connectivity in the coal SH-B is better than the coal SH-V. It’s important to note that microscopic fractures in coal SH-V are filled by minerals, which is an important factor for the poor connectivity of coal SH-V. We can assume that without mineral filling, the connectivity of coal SH-V may be better than coal SH-B, because that the contribution of microscopic fractures to coal connectivity is superior to that of pores.

Fig. 11
figure 11

Schematic diagram of the fractures in coal seam #3 in the southern Qinshui Basin (modified from Liu et al. 2018b). Notes: a Fractures parallel to the bedding plane; b fractures perpendicular to the bedding plane; σmax, maximum horizontal principal stress; σmin, minimum horizontal principal stress; σv, vertical principal stress

4.3.2 Influence spatial distribution of fractures on connectivity

The cleats and microfissures of coal seam #3 in the southern Qinshui Basin are also directional and are usually perpendicular to or nearly perpendicular to the bedding and extend parallel to the bedding plane (Fig. 8). Therefore, the cleats and microfissures usually extend long distances parallel to the bedding plane, even throughout the coal pillar (Fig. 11a), which greatly improves the connectivity in the coal samples. These findings cause the effective porosity of coal YW-B to be comparable to that of coal SH-B or even better, although the NMR porosity and X-ray CT porosity of coal YW-B is only 1/2 of that of coal SH-B. Perpendicular to the bedding plane, the height of the cleats and microfissures in the coal samples are small and cannot intersect the coal pillar. In addition, the macropores are not abundant, and forming connections between the cleats and microfissures is difficult (Fig. 11b), resulting in the extremely low connectivity of coal YW-V. In conclusion, the connectivity of the Yuwu coal samples is mainly dependent on the directionality of the microfissures.

5 Conclusions

In this paper, an integrated approach consisting of LF-NMR, X-ray CT, and FESEM was utilized to study the characteristics and causes of high-rank coal structure and connectivity parallel and perpendicular to the bedding plane. The following conclusions can be drawn from this study:

  1. 1.

    High-rank coal predominantly contains mesopores both parallel and perpendicular to the bedding plane, and the pore volume decreases with increasing pore width. Mesopores in high-rank coal are mainly shrinkage-induced pores and secondary gas pores. Due to the significant influence of the maximum horizontal principal stress, the morphology and spatial distribution of pores, fractures and minerals are strongly oriented perpendicular to the bedding plane.

  2. 2.

    The difference in the connectivity of high-rank coal is mainly due to the throat development and the orientation of pores and fractures. The pores extend along the minimum horizontal principal stress direction, while the cleats and microfissures mainly extend parallel to the bedding plane, resulting the connectivity parallel to the bedding plane is significantly greater than that perpendicular to the bedding plane.