The combined effect of fractures and mineral content on coal hydromechanical response

The hydromechanical behaviour of fractured coal is a complex function of interaction between coal bulk and fracture deformation driven by fluid pressure and external stress. Despite the research studies conducted to date, the combined effect of mineral content and fracture structure on hydromechanical behaviour of sorptive fractured coal remains unexplored. To study this combined effect, we performed a series of X-ray computed tomography (XRCT) imaging on a range of coal specimens with non-sorbing (helium) and sorbing (CO2) gases at different effective stress paths using a newly developed X-ray transparent triaxial system. The compressibility of system components was obtained from processed 3D XRCT images which were used to interpret the results. The results of this study show that coal matrix/solid compressibility has a positive nonlinear relation with mineral content irrespective of mineral type. Effective stress coefficient is also a strong function of both mineral content and fracture porosity. Furthermore, the increase in mineral content leads to less fracture opening by an increase in helium pressure. Interestingly, the effect of mineral content on the bulk strength of coal is more significant than the effect of fracture porosity. Finally, coal with more open fractures shows less bulk swelling by gas adsorption under external stress due to damping effect of fracture volume on developed internal volumetric swelling strain.


Introduction
The enhanced coalbed methane recovery and CO 2 storage in coal seams has gained a broad interest in recent decades (Detheridge et al. 2019;Karacan et al. 2011;Shi and Durucan 2005). In the coalbed methane recovery, the pressure-driven fluid flow often combined with hydraulic fracturing and N 2 -CO 2 injection is used to produce the methane. The geological storage of carbon in coalbed seams can be also achieved through this technique (Bachu 2003;Gou et al. 2021). The CO 2 injection for enhanced coalbed methane leads to swelling of the coal matrix which can in turn reduce the coal permeability significantly (Bertrand et al. 2019;Li and Fang 2014). The interaction between bulk of solid with fractures and fluid constitutes is thus of critical importance during coalbed methane production. Such interaction is complex, especially when sorbing gas is involved. Amongst different aspects of these interactions, the hydromechanical response of fractured coal in the presence of non-sorbing and sorbing gases needs special attention. The study of hydromechanical response of geomaterials goes back to the work of Biot (1941) and Terzaghi (1943). Their theories have been extensively used in rock and soil mechanics for analysis of the rock hydromechanical behaviour. Following Biot's work, poroelastic behaviour of geomaterials was further explored by many researchers (Detournay and Cheng 1993;Gassmann 1951;Rice and Cleary 1976;Zimmerman et al. 1986).
The Biot linear poroelasticity was defined for porous media consisting of grains and pores (Lv et al. 2019). The coal structure, however, contains fractures with micrometers aperture (cleats) and matrix that contains a range of pores with majority sitting in nanometer scale. In a recent study, Lv et al. (2021a) showed that the Biot poroelasticity can be used for fractured coal considering matrix as the solid and fractures as void constitute. Recent studies (Lv et al. 2021b(Lv et al. , 2022 also indicated that coal fractures are highly sensitive to external stress and pore pressure even for stresses lower than 1 MPa, while matrix micropores require much higher level of stress to deform (above 10 MPa). When sorbing gas is considered, the gas stays at bulk condition in the fractures, while it is mostly adsorbed in matrix micropores; thus, hydromechanical response (no sorption involved) is dictated by the fractures' behaviour which depends on stress and pore pressure of the system.
Hydromechanical responses of coal especially from porosity and bulk volume perspective have been studied in the literature (Cheng et al. 2021;Tan et al. 2019;Wang et al. 2020;Zhang et al. 2019) concluding that the higher porosity leads to higher compressibility of the coal bulk. For instance, Shao et al. (2018) showed that the matrix compressibility has a strong correlation with porosity and mineral content in coal. Golsanami et al. (2019) differentiated the fracture porosity from matrix pores and studied its relation to bulk and Young's modulus calculated from acoustic velocities. The impact of mineral matter on hydromechanical responses of coal can be attributed to mineral content, mineralogical composition and morphology (Ahamed et al. 2019). The effect of mineral content and composition on mechanical properties of coal has been investigated in the literature showing that the high-density minerals can stiffen the coal matrix and significantly decrease its bulk compressibility. In an early study, Barbour and Ko (1979) investigated the effect of carbonate and clay contents on the Young's modulus of coal and found that an increase in carbonate content of coal results in an increase in Young's modulus while the clay content has an opposite effect. Guo et al. (2014) tested 21 coal samples of different ranks and concluded that porosity and organic matter content were positively correlated to the matrix compressibility, except for high-volatile group with low minerals content. This important finding from experimental observation was also confirmed by Wetzel et al. (2018) through simulation on synthetic digital rocks, where a positive correlation between calcite content and the bulk modulus was found. Cai et al. (2018) concluded that clay mineral content shows negative exponential correlation to matrix compressibility from the obtained experimental results of coal samples with mineral contents less than 20%. Following these experimental results, Song and Liu (2020) used carbonate solution to improve the coal strength. In addition, Miao et al. (2020) performed hydromechanical experiments on different coal samples leading to similar conclusion that the coal matrix compressibility decreases with an increase in the minerals content. The possible effect of mineral structure on coal bulk compressibility, however, has not been fully understood due to limitations of traditional methods in tracking mineral structures, especially when the mineral content is low. To investigate the coal mechanical properties through morphological characteristics, nondestructive methods such as X-ray Computed Tomography (XRCT), Scanning Electron Microscopy (SEM), and ultrasonic transmission, are used.
In particular, the rapid development of laboratory XRCT over the past decades has enabled studying different processes in coal (Mathews et al. 2017). The XRCT is occasionally combined with other methods such as SEM (Gupta 2007) to attain more detailed information. This information is generally used to model permeability Ramandi et al. 2021;Wang et al. 2019b;Zhao et al. 2018;Zhu et al. 2021), sorption effect (Karacan and Okandan 2001;Pirzada et al. 2018;Shi et al. 2020) and mechanical deformation of coal (Dong et al. 2021;Lv et al. 2019;Stappen et al. 2019;Wang et al. 2021;Zhao et al. 2019). To investigate the mechanical deformation and fluid transport in porous media, XRCT images are used in simulation exercises (Blaheta et al. 2015;Chen et al. 2019;Roslin et al. 2019;Wang et al. 2019a;Wei et al. 2021;Zhang et al. 2017;Zhao et al. 2019). An example of recent hydromechanical simulation on XRCT images is the development of a novel method to calculate Biot's coefficient in a granular media (2021). Using image segmentation, fractures, coal organic matter and mineral contents can be differentiated to investigate their interactions with each other. For example, Zhang et al. (2018) reported that Young's modulus of coal is highly influenced by the mineral content through XRCT and acoustic measurements. Li et al. (2020) observed a noticeable difference between mechanical properties of coal matrix and mineral through 2D and 3D XRCT image analysis.
Sorption effect is another distinct natural behaviour of coal; the sorption capacity is closely related to coal composition, especially carbon content and micro-structure influencing gas accessibility (Cheng et al. 2017;Crosdale et al. 1998;Hou et al. 2017). Due to microscopic structure of carbon molecules, organic matter content has a positive effect on the sorption capacity (Chalmers and Bustin 2007;Mastalerz et al. 2004;Meng and Liu 2018). The mineral matters, on the other hand, have a negative impact on the sorption capacity (Laxminarayana and Crosdale 1999) due to their low adsorption capacity compared to organic matters (Clarkson and Bustin 1996) although minerals like clay can still contribute to sorption capacity (Deng et al. 2015). Pone et al. (2009) suggested that mineral distribution can influence the connectivity of the microstructures and continuity of the pore networks thus affecting gas sorption process. It is clear from above discussion that investigating the combined effect of mineral content and fractures on hydromechanical response of coal requires hydromechanical testing coupled with methods such as XRCT.
In this study, we therefore employed micro-XRCT (µCT) technique to investigate the hydromechanical response of different coal samples though analysing bulk and fracture volumes with non-sorbing and sorbing gases using a newly designed X-ray transparent triaxial system (Roshan et al. 2019a). To investigate the contribution of different mineral contents on hydromechanical behaviour of coals, images were segmented into fracture, organic matter, and mineral domains and the volumetric strain of the bulk and fractures were extracted to calculate the compressibilities of fracture and matrix as well as the effective stress coefficient. The imaging results were then compared with a measured effective stress coefficient of one of the coal specimens to assess its validity. The obtained compressibilities of fracture and matrix as well as the effective stress coefficient were then linked to mineral content and fracture porosity. The XRCT imaging was also performed with an adsorption test to investigate the effect of adsorption on the hydromechanical process.

Characterization of hydromechanical properties
In order to systematically analyse the effect of mineral content and structure on hydromechanical response of fractured coals, we follow the definition of compressibility of the porous rock constitutes introduced by Zimmerman et al. (1986). This methodology allows us to link the mineral content to different compressibilities of the fractured coal. Knowing that the fractures are much more compressible than micropores in the coal matrix, especially under a low level of external stress, and by analogy to a pore-solid system, the volume of the bulk of fractured coal ( V b ) is defined as sum of the fracture volume, V f and matrix volume, V m including solids and its micropores.
Fracture porosity, f is then defined as the ratio of fracture volume to bulk volume: Having two phases (matrix and fracture), the system compressibility functions for the fractured coal are defined according to Zimmerman et al. (1986) for porous rock (compressive stress and strain are assumed positive) including: (a) bulk compressibility at different boundary conditions (BCs): where and p are the total stress and pore pressure, respectively. Also C bc and C bp are bulk compressibilities measured at constant pore pressure and constant external stress, respectively. (b) fracture compressibility at different BCs: where C fc and C fp are the fracture compressibilities measured at constant pore pressure and constant external stress, respectively. In addition to above compressibilities of system components, the solid (matrix) compressibility, C s is also defined as: In a laboratory measurement, C s is measured from an unjacketed test where the differential pressure P d is kept constant. P d is defined as the difference between the total stress and the pore pressure.
Biot's concept of effective stress law (Biot 1941) defines the effective stress, ′ of a porous medium by relating the total stress, to pore pressure, p through the effective stress coefficient, (or ESC): While the Biot poroelasticity was defined around solid constitute and pores, recent study by Lv et al. (2021a) (1) showed that it can be used for fractured coal, where coal matrix is considered as the solid constitute and fractures as voids. The can be measured using jacketed-unjacketed test (Biot and Willis 1957;Geertsma 1957;Skempton 1961). In Eq. (3), C bp is measured from a hydrostatic compression on a jacketed specimen under Δp = 0 . C s is measured from a hydrostatic compression under ΔP d = 0 as in Eq. (5). The then describes the effect of void volume on the bulk deformation which ranges from 0 to 1.
Another way of measuring the is through a hydrostatic compressive test under Δp = 0 . From Eq. (8), is determined by the change in fracture volume, ΔV f and bulk volume, ΔV b leading to (Franquet and Abass 1999):

Sample preparation and characteristics
All coal samples used in this study were collected from different boreholes from a coal seam in Sydney Basin, New South Wales, Australia (Zhao 2012). Before coring, a layer of resin was applied on the surface of coal samples allowing at least one day for plasticising. Then, a half-inch (12.7 mm) coring bit was used to obtain specimens with a diameter of about 12.5 mm. The coal specimens were then cut into a length of approximately 2 times of the diameter (between 19 to 28 mm depending on the brittleness of the coal sample). Finally, specimens were evacuated in the oven at 105° Celsius for at least 2 days to remove any moisture.
Six coal specimens labelled as SY1 (one specimen), SY2 (one specimen) and SY3 (four specimens labelled as SY3-1 to SY3-4) with varying fracture porosity and mineral content were used for the analysis. Fracture porosity refers to the ratio of fracture volume per bulk volume of the specimen and mineral content refers to minerals volume per bulk volume of the specimen. The SY2 specimen had the lowest mineral content of 2.49% volume ratio, followed by SY1 with 8.76%. The SY3 specimens had a high mineral content of average above 30% showing a distinct minerals content difference with SY1 and SY2. In addition, the unconfined fracture porosity was around 0.8% and 0.7% for SY1 and SY2, respectively, whereas SY3 specimens had different fracture porosity ranging from approximately 0.3% to 0.8%. Specimens with different fracture porosity and mineral content were used for different tests that will be further explained in the experimental methodology. The leftover powder from the samples was used for X-ray diffraction to identify the mineral composition of each category (SY1 to SY3).

X-ray diffraction
X-ray diffraction (XRD) was carried out on three coal samples from each category. Samples were measured in PANalytical MPD Powder diffractometer, with Cu X-ray source, at 45 kV and 40 mA X-Ray generation condition. XRD results of three samples are summarised in Table 1 showing four groups of minerals in coal samples, including clay, carbonate, quartz and phosphates. These minerals have different properties in coal seams. As shown in Table 1, in SY1, the dominant mineral was clay (kaolinite) with minor quartz and carbonates (dolomite and aklimaite). Similarly, in SY2, the major mineral was kaolinite with minor quartz, and traces of siderite, aklimaite and hydroxylapatite. The mineral content of SY3 was different mainly composed of carbonates, including the dominant ankerite, traces of siderite, and minimal content of quartz. In general, kaolinite is a commonly occurring mineral in coal as epigenetic cleat infillings (Spears and Caswell 1986). Quartz is the most common oxide mineral occurring in coal and largely appears as discrete grains of detrital origin. Phosphates often coexist with clay minerals, mainly present in the form of cell and pore infillings. In addition, a wide range of carbonate minerals have been reported in coal with most found as siderite, calcite, dolomite and ankerite. Ankerite appears in cleat infillings and broken cell cavities within organic matters. It is noted that the volume of whole minerals in each experimental specimen was measured from µCT images.

Triaxial μCT imaging
The schematic view of the X-ray imaging triaxial system (Roshan et al. 2019a) is shown in Fig. 1. The system has independent controls for axial and confining stresses with fluid flow access to the specimen through the inner hole at both end platens. Similar to a traditional triaxial cell, the coal specimen was hosted in a rubber sleeve to isolate it from confining oil. The axial load was applied through a micro-ram which placed on the top platen of the triaxial cell. In the axial direction, the system has both displacement and load control, but only load control was used in this study. Three ISCO pumps were used to control axial, radial and fluid pressures and the pressure transducers were connected to the pressure lines. The fluid flow section allowed both drained and undrained testing conditions by opening or closing the downstream outlet at the bottom platen. The X-ray scanning used in the study is in helical mode, meaning that the triaxial cell rotates along its axis while the X-ray source and the detector maintain stationary. An X-ray source with a voltage of 100 kV and a current of 120 μA were used for scanning. This setup yields an image resolution of 12.4 μm with about 2-h scanning time on the specimen (Chen et al. 2020;Pirzada et al. 2021). To study the relation between mineral content and hydromechanical properties, images should be processed and segmented into three domains: mineral, organic matter and fracture. The image processing will be discussed later. To measure the bulk compressibility, C bc for example, we scanned the specimen twice under different confining pressures, and then measured the volume difference of the two obtained images. The fracture porosity reported in the plots is those at the second stage of testing (specimen compressed). Table 2 summarizes the stress-pressure conditions of XRCT imaging to obtain the hydromechanical properties (compressibility of system components including bulk, fracture and solid constitutes). We conducted three groups of tests to determine the system compressibilities and find their relationships with samples' mineral content and fractures including Unjacketed BC, Drained BC1 and Drained BC2. In unjacketed test, C s was determined by maintaining a constant P d as in Eq. (5). The first drained test was conducted at a constant p and the second drained test was performed at a constant σ. The bulk and fracture compressibilities were studied in two drained tests with different boundary conditions (Eqs. 3 and 4). An extra set of tests were carried out with both sorbing and non-sorbing gases to investigate the swelling effect on hydromechanical properties.

Unjacketed BC
The original unjacketed test (Biot and Willis 1957) involves immersing the sample in pressurised fluid, so the fluid can penetrate through pores of the sample. This method was later modified to an unjacketed test using a triaxial cell under p and Gas type undrained conditions. As shown in Fig. 2, a constant differential pressure P d is maintained throughout the test, with both and p increasing by the same amount and stay close to each other. To use this method in XRCT imaging experiment, two scans were performed (see Fig. 2). By measuring the applied volumetric strain from both images, the solid compressibility, C s can be calculated using Eq. (5) which is the gradient of volumetric strain to stress. Three specimens (SY1, SY2 and SY3-1) were used for these experiments under identical BCs, as shown in Table 3. The three specimens had a range of different mineral contents (~ 2%-26%) and type of minerals (clay-rich versus calcite-rich). Different mineral types were chosen for the analysis to capture their effect on overall solid compressibility of the system. It is also noted that the fracture porosity of the specimens is not central in analysis of solid compressibility.

Drained BC-1
This group of tests aims to acquire the bulk and fracture compressibilities, C bp and C fp using the boundary condition specified in Table 2. In these experiments, p was increased, while was maintained constant, so the volume change of bulk and fractures due to pore pressure increase could be characterized. Figure 3 shows the test procedure including the scanning sequence. It is noted the hydrostatic confinement was kept constant but different in each test while the pore pressure was increased. Three specimens were tested in this group as shown in Table 4. All specimens were from category SY3 (SY3-1, SY3-3, and SY3-4) having different mineral contents (~ 25%-35%) and fracture porosity (~ 0.3%-0.8%) but the same mineral constituents (calcite-rich). As the fracture compressibility is of particular importance in this test, a range of specimens with high fracture porosity variation was chosen.

Drained BC-2
In the Drained BC-2, the gas was allowed to flow out of the sample during external loading ( ) while keeping the pore pressure, p constant. The change of the bulk and fracture volumes were then measured. As the fluid pressure Fig. 2 Test design in Unjacketed BC using two scans to measure the solid compressibility at constant differential pressure  Fig. 3 Test design in Drained BC-1 using two scans to measure the bulk and fracture compressibility at constant hydrostatics pressure throughout the loading is identical, Δ V f represents the change of fracture volume resulted from the incremental loading. From this group of tests, the drained bulk compressibility ( C bc ), the fracture compressibility ( C fc ) and ESC ( ) can be determined. Figure 4 shows the loading path and the sequence of two scans at this BC. The outlet line was opened to atmospheric pressure thus the pore pressure was maintained constant. Hydrostatic confinement was increased from the first scan to the second, to compress the specimen. Comparing the volume difference in the bulk and fractures from the two scans, C bc and C fc can be calculated. Two specimens (SY3-1 and SY3-4) were used in the experimental design of this group, as shown in Table 5. These two specimens have similar mineral-type (calcite-rich), fracture porosity (0.66% and 0.62%, respectively) and mineral content (25.4% and 29.3%, respectively) but different mineral structures (isolated versus spread) chosen with the view of investigating mineral structure.

Adsorption BC (sorption-hydromechanical effect)
An additional group of tests was carried out on three specimens (SY2, SY3-2 and SY3-3) with sorbing gas to investigate the effect of swelling on hydromechanical response. It is noted that these three specimens had different fracture porosity (~ 0.3%-0.8%), different mineral contents (~ 2.5%-32%) and different minerals type (clay-rich and calcite-rich) thus a combination of different factors on swelling behaviour was considered. It is well known that carbon has an adsorptive nature and coal shows swelling when exposed to sorbing gas. The specimens were first saturated by helium and then by CO2 and they were scanned in turns at the BCs shown in Fig. 5. As the gas adsorption amount is dependent on the pore pressure, a consistent pore pressure of 2.5 MPa was applied to eliminate the pressure effects. Time (1-2 days) was given to CO2 to saturate the specimen and reached an adsorption equilibrium. Three specimens used in this group and their experimental BCs are presented in Table 6. 3.5 SY3-4 1 5 MPa 0 2 4.5

Fig. 4
Test condition in Drained BC-2 using two scans to measure the bulk and fracture compressibility at the constant pore pressure

Image processing
After performing experiments, the obtained 3D XRCT images were processed and analysed. In this study, the commercial software AVIZO (FEI Group 2013) was used for volume editing, registration, denoising, segmentation, etc. As a first step, the images were denoised using multiple filter, e.g. the non-local means filter was used to replace the grey value of a pixel with an average of the grey values of nearby pixels (Buades et al. 2005;Matsuoka et al. 2017). The images from each set of experiment were then aligned geometrically for registration purpose which can additionally enable investigating variation in the volume of each component. Images were next segmented into three domains of mineral, organic matter and fracture. Segmentation was performed by combined methods of grayscale thresholding, watershed algorithm and magic wand tool (Thermo-Fisher-Scientific 2018). Grayscale thresholding is a simple approach to segment different domains using grayscale values. Watershed is an interactive thresholding which provides a statistical algorithm to assign each pixel to a region or a watershed (Beucher and Meyer 2018;Preim and Botha 2013). Magic wand tool is a quick method to segment irregularly shaped regions based on edges. Suitable application of combined processing methods yields an ideal segmentation. Figure 6 shows typical segmented domains of a coal specimen from a 3D image.

Bulk deformation calculation
To measure the specimen deformation at different BCs, the diameter and length of specimen at the corresponding reconstructed 3D images were measured, as shown schematically in Fig. 7. A pair of tomography images from orthogonal planes were extracted from the 3D image, while the specimen diameter and height were measured from the horizontal (slice A in Fig. 7) and vertical (slice B in Fig. 7) cross sections, respectively. Multiple measurements were performed, and the average value was considered in the analysis. The accuracy of the measurement sits within the obtained image resolution of 12.4 μm.

Validation study
To assess the validity of measured hydromechanical properties from XRCT images, the bulk and solid compressibilities, as well as ESC of the specimens' category SY3, were measured using conventional triaxial experiments. It is noted that this validation exercise with one specimen is to ensure that the measurements and calculation are of high quality as the accuracy of image-extracted properties including deformation has been already proven (Pone et al. 2010;Salemi et al. 2021;Sufian and Russell 2013;Viggiani and Hall 2008). To do so, we modified the previously in-house designed triaxial cell (Roshan et al. 2019b) to adapt the specimen with a diameter of 0.5 inch (~ 12.5 mm) as shown in Fig. 8. The axial load on the sample was applied through a servo-controlled loading frame with capacity of 500 kN. Axial deformation was measured through the LVDT along the axial direction. Radial deformation was measured using the internal LVDTs. These deformations were logged using a LabView interface and recorded in the computer. A manual pump with hydraulic oil was used for confining pressure and gas cylinder with a regulator for the fluid pressure. To obtain the ESC using the triaxial system, we measured bulk modulus in a drained experiment and solid bulk modulus in an unjacketed experiment, as summarised in Table 7 and detailed in Fig. 9. It is noted that the bulk modulus is inverse of bulk compressibility. In the drained experiment, the confining pressure was increased by Δ , while the pore pressure p was maintained constant. The confining pressure was increased by applying a loading rate of 0.18 mm/min, and the volumetric strain was continuously measured. The derivative of with respect to V b is the bulk modulus, K b . The solid modulus, K s was measured in a jacketed experiment, with a boundary condition of constant differential pressure, P d as previously discussed. As shown in Fig. 9, after each pressure increase ( Δp ), the difference between the confining pressure and pore pressure was maintained unchanged (close to each other) which allowed the volumetric strain reaching equilibrium, enabling the measurement of ΔV b . The pressure was plotted versus volumetric strain obtained at each step, and then, the effective stress coefficient was determined following Eqs. ((3), (5)- (7)).

Comparison of measurements using triaxial cell and XRCT imaging
As noted above, the bulk modulus was measured through a uniform loading where the K is the gradient of stress to volumetric strain. The solid bulk modulus was measured in a stepwise test, where a constant differential pressure was applied to calculate Ks. These two gradients measured from both measurements are presented in Fig. 10. The mechanical tests yielded a bulk modulus of 1.019 GPa and a solid bulk modulus of 8.8 GPa. Therefore, the effective stress coefficient was calculated to be 0.88.
The same set of above parameters were calculated from XRCT images. As shown in Fig. 11, the first image was taken without stress and the second under 5 MPa uniform stress. Comparing 2D tomograph from two images, as highlighted in red circle on images, the fracture porosity and fracture size in the first image are higher than that of the second image. From the raw image slices shown, the fracture aperture reduces from around 50 microns to less than 13 microns, and the number of fractures reduces from more than 10 to only 2 small discrete fractures. After segmentation, the fracture domain was extracted from two images showing an evident closure of fractures induced by stress.
Having the fracture and bulk volumetric strains, the effective stress coefficient is calculated as 0.89 using Eq. 8. Similarly, bulk modulus can be calculated using Eq. 3 through the uniform compression and solid bulk modulus can be inferred from Eq. 7. Results from both approaches are presented in Table 8 showing a high level of confidence in the methodology adopted.

Hydromechanical response
Unjacketed BC (constant differential pressure) In the unjacketed testing group, two scans were taken under the same differential pressure for specimens SY1, SY2 and SY3-1 having different mineralogy and minerals content. Figure 12 illustrates the stress and helium pore pressure conditions of the tests plus their tomography images for each specimen, i.e. the first scan was taken at 0.5 MPa stress with no pore pressure, and the second scan was taken at a total stress of 3 MPa with helium pressure of 2.5 MPa for all specimens. To precisely calculate the fracture volume and investigate the relation between the sample compressibility and mineral content, images were segmented into three domains: fracture, mineral and organic matter (Fig. 13). As specimens were cored from different samples, they showed distinct mineral compositions and structure. Figure 12 shows that SY1 and SY2 have relatively low amount of mineral (8.7% and 2.5% volume ratio, respectively, from XRCT analysis) with kaolinite clay being the dominant mineral from XRD analysis (Table 1). Minerals in SY1 have a lamellar structure, while SY2 has scattered and dotlike minerals Fig. 11. 2D tomograph images and segmented fractures of specimen SY3-4 (atmospheric pore pressure) distribution. From the tomography images (see Fig. 12), kaolinite appears to be well ordered and in bands and laminae due to its in situ leaching and reprecipitation within the peat swamp (Ward 1989). SY3-1, however, contains large amounts of carbonates with authigenic structure (25.2% volume ratio from XRCT analysis (see Fig. 13)). They were formed during deposition and coalification process (Taylor et al. 1998;Ward 1992) and seem to be infillings of fractures as well as veining in the matrix structure. This has a high potential to stiffen the bulk of coal and reduces its compressibility. From these experiments, the solid compressibility, C s or its inverse, the solid bulk modulus, K s of different categories were measured and presented in Fig. 14 as a function of mineral content. Interestingly, the K s shows a positive nonlinear trend with an increase in the mineral content irrespective of the mineral type. It emphasizes that organic matter compressibility is higher than other minerals whether it is composed of soft (clay) or hard (carbonate) minerals. Different from those of SY1 and SY2, minerals in SY3-1 are majority carbonate with high stiffness and low compressibility compared to clay minerals in other samples. Despite the higher stiffness of calcite compared to clay minerals and that SY3-1 contains much higher minerals than SY1 and SY2 (25.2%, 2.5% and 8.7% for SY3-1, SY2 and SY1, respectively), the rate of increase in compressibility with mineral content reduces at higher mineral volumes. It has been inferred that ions in the carbonates originated from organic matters during coal rank advancement mainly fill in the coalbed fractures through fluid migration (Spears and Caswell 1986).

Drained BC-1 (constant hydrostatic stress)
Tomograph slices along with 3D segmented images from Drain BC-1 for three specimens of the same category SY3 (SY3-1, SY3-3 and SY3-4) are shown in Fig. 15, where different pore pressures of helium at a constant hydrostatic pressure were applied on each specimen. It is noted that the pore pressure and stresses were different for each specimen but the condition of constant confining (hydrostatic) stress at varying pore pressures was maintained. All the specimens were prepared from the same sample to present a similar type of mineral (calcite-rich). Figure 15 shows that SY3-3 has the highest mineral content (34.21% volume ratio) compared to others (25.2% and 29.3% volume ratio for SY3-1 and SY3-4, respectively). An increase in the pore pressure leads to opening of fractures in all specimens but to a different extent. Interestingly, the fracture opening by fluid pressure follows the mineral content; higher the mineral content, the less fracture opening by increase in pore pressure. The fracture opening in SY3-1 was the highest, followed by SY3-4 and SY3-3. A gradual decrease was observed in both C fp and C bp with an increase in the mineral content, while an increasing trend was observed with fracture porosity (Fig. 16). Similar positive correlation between C fp and fracture porosity has also been reported by Zhu et al. (2018) and Bagherzadeh et al. (2021). The results, however, suggest that fracture compressibility is much more sensitive to increase in pore pressure than the bulk compressibility. This behaviour also emphasises that the effect of minerals on the hydromechanical behaviour of fracture and bulk is different. In fact, the findings of this study evidently show that an increase in mineral content stiffens the coal skeleton more efficiently than the absence of fracture.

Drained BC-2 (constant pore pressure)
Similar to the Drained BC-1, both fracture and matrix compressibilities, C fc and C bc were measured but at a constant pore pressure (atmospheric) with varying hydrostatic stress (see Figs. 17 and 18). The ESC was also calculated in this analysis. Both specimens used in this analysis, were highly cleated with high fracture porosity (0.66% and 0.62%) and a minor difference in mineral content (volume ratio of 25.4% and 29.3% for SY3-1 and SY3-4, respectively, Fig. 17). The increase in hydrostatic stress resulted in a decrease in the fracture porosity of both specimens, but at different rates. In SY3-1, the volume of fractures was uniformly reduced while in SY3-4, two main fractures were remained partially open after compression. This difference resulted in a greater value of bulk compressibility in SY3-1. The fracture compressibility of both specimens was, on the other hand, considerably different. The non-uniform fracture closure in SY3-4 was controlled by the structure of minerals. While fractures in SY3-1 were more evenly distributed in the bulk of the specimen, in SY3-4, fractures were preferably aligned and relatively isolated. Thus, the reduction of volume due to compression in SY3-4 was localised. The ESC was also found to be dependent to both the fracture porosity and mineral content. The higher fracture porosity and lower mineral content leads to a higher ESC (Fig. 18). This is intuitive as the increase in the mineral content and reduction of fracture porosity will lead to an increase in the load bearing capacity of the skeleton and consequently reduction of the bulk compressibility.

Adsorption (sorption-hydromechanical effect)
To investigate the combined effect of minerals and fractures on bulk and matrix deformation in presence of a sorbing gas, a set of experiments were conducted on specimens SY2 and SY3-2 and SY3-3 (Table 9). It is important to note that the mineral content of specimen SY2 (2.5% Fig. 14. a Ks as a function of mineral content (volume ratio) and b Cs as a function of mineral content (volume ratio) in unjacketed BC.

Fig. 15
Fracture domain in three specimens used in Drained BC-1 volume ratio) was considerably lower than the two specimens from category SY-3 (31.5% and 34.2% volume ratios for SY3-2 and SY3-3, respectively). The gas pore pressure of 2.5 MPa was first applied on each specimen, using a non-sorbing gas (helium), followed by injection of sorbed gas (CO 2 ) at constant hydrostatic pressure of 3 MPa. The 2D tomograph images and segmented fractures are illustrated in Fig. 19 showing that the SY3-3 specimen contains the highest number of fractures while the others contain a single main fracture through the bulk of the  specimens (the specimens have fracture porosity ranging from ~ 0.3%-0.8%). The sorbing gases cause the coal to swell, resulting in an increase in the matrix and bulk volumetric strain and a reduction in internal fracture volume. By calculating solid volume excluding fractures, we obtained the bulk and matrix deformation between the treatment of helium and CO 2 and plotted it versus mineral content in Fig. 20. It can be observed from this figure that SY2 has the largest swelling strain both in bulk and matrix as it contains a negligible amount of mineral (2.5% volume ratio) and high content of organic matter (96.8% volume ratio) with a single fracture (with high porosity due to its fracture aperture). Specimens from SY3 category showed lower matrix-bulk swelling. This is attributed to presence of high volume of carbonate minerals that reduces the micropore volume (Mastalerz Fig. 18 Results in Drained BC-2; a1) C bc as a function of mineral content, a2) C bc as a function of fracture porosity, b1) C fc as a function of mineral content, b2) C fc as a function of fracture porosity, c1) effective stress coefficient as a function of mineral content, and c2) effective stress coefficient as a function of fracture porosity This is in agreement with the findings of previous studies (Zhang et al. 2020).
It is interesting to note that the greatest difference in swelling strain between the bulk and matrix was found in SY3-2, while minimal differences were observed in other specimens. The aperture (opening) of fractures in SY3-2 was much greater than those of other specimens (highest fracture porosity). It implies that specimens with widely open fractures show less bulk swelling as fractures tend to absorb the coal matrix swelling internally thus leading to less bulk swelling when under external stress. Overall, the volumetric strain is highly influenced by both mineral content and fracture structure in the adsorption process.

Conclusions
In this study, the effect of mineral content and structure together with their combined effect of fracture porosity on the hydromechanical response of fractured coal specimens was experimentally evaluated. With the combined method of XRCT visualization and triaxial testing, we performed unjacketed and drained tests and calculated the effective stress and a range of system compressibilities through image processing. This image-based calculations were validated with a standard triaxial experiment. The results of this study show that the solid bulk modulus of coal K s has a positive nonlinear relationship with mineral content irrespective of the mineral type. In general, the coal fracture compressibility is negatively correlated to the mineral content, while the matrix compressibility is positively correlated to the fracture porosity. Specially, carbonates as infillings strengthen the coal matrix. Specimens in this study show that kaolinite are well ordered and appears in bands and laminae while carbonate appears in authigenic pattern. As a result, specimens with low content of minerals containing dominant kaolinite are more compressible than those containing high content of carbonate. It was also found that the fracture opening by internal fluid pressure follows the mineral content; the higher mineral content, the less fracture opening due to increase in pore pressure. Also, the results suggest that the change in fracture compressibility is much more sensitive to pore pressure increase than bulk compressibility. Interestingly, the results further indicate that the minerals can considerably stiffen the skeleton much more efficiently than the absence of fracture in the coal specimens. The effective stress coefficient is also found to be a complex function of not only the fracture porosity but also the mineral content. The higher fracture porosity and lower mineral content leads to a higher ESC. Finally, coal samples with more open fractures in the presence of adsorbing gas show less bulk swelling despite the higher swelling potential of their matrix. This shows that specimens with widely open fractures experience less bulk swelling as fractures tend to absorb the coal matrix swelling internally thus leading to less bulk swelling when under external stress. The findings of this study have significant implications in field studies such as gas production and injection in coal seams where the hydromechanical response of coal is important.
Funding Open Access funding enabled and organized by CAUL and its Member Institutions.

Fig. 20
Effect of mineral content and fracture porosity on swelling strain of three specimens; a) matrix and bulk swelling as a function of mineral content, and b) differential swelling strain between matrix and bulk as a function of fracture porosity

Declarations
Competing Interests Authors declare that there is no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.