Nanoindentation study of macerals in coals from the Ukrainian Donets Basin

This article discusses the mechanical properties of macerals in Carboniferous coals from the Ukrainian Donets Basin, covering a maturity range from 0.62 to 1.47 %Rr (vitrinite reflectance). The inherent inhomogeneity of geological samples requires characterization at the micro-/nanoscale, and hence material parameters, such as hardness H and reduced elastic modulus Er, were obtained from twelve coal specimens via nanoindentation tests. Different material properties and maturity trends were acquired for the individual maceral groups (vitrinite, inertinite, liptinite). The results indicate that apart from maturity, multiple factors like microstructural features and depositional environment control the mechanical properties of macerals in coals. H and Er of vitrinites show a complex relationship with maturity due to additional influence of microstructural features such as nanoporosity or mineral inclusions and complex changes in the organic matter structure with increasing rank. The mechanical behaviour of inertinites is mainly controlled by the prevailing conditions (temperature, exposure) during paleowildfires. Er and H of liptinites are strongly influenced by transformational processes related to devolatilization at low to medium rank, as well as later hydrocarbon generation. Variations in the depositional environment (e.g. increased organic S content due to marine influence) might play an additional role.


Introduction
Coal is still of particular importance as a natural resource today and remains one of the most important energy sources for the production of electricity or in various metallurgical processes.Moreover, it is considered a substantial source for methane which not only represents a mining hazard but can also be utilized by modern unconventional production techniques (e.g.coal bed methane; CBM).However, the producibility of coal bed methane depends on mechanical properties and fracture behaviour of the hosting coals.Coals are constituted by distinct organic components referred to as macerals (vitrinite, liptinite, inertinite), and obtaining the mechanical properties of individual maceral particles is challenging due to their usually small size (often microscaled).Nevertheless, knowledge of maceral-specific properties is critical for the understanding of controlling factors on the overall mechanical properties of a certain coal sample.Borodich et al. (2015) and Epshtein et al. (2015) proved the suitability of nanoindentation for the determination of mechanical properties of individual coal macerals.The great advantage of nanoindentation compared to other indentation testing approaches is the small required testing volume which allows analyses of material parameters like hardness (H ) and reduced elastic modulus (E r ) in areas that are considered highly inhomogeneous at the micro-scale (Oliver and Pharr, 1992).
In an attempt to determine material parameters of individual macerals and to study maturity-related changes, nanoindentation tests were carried out on twelve Carboniferous coal samples from the Ukrainian Donets Basin (Fig. 1).The Donets Basin hosts one of the most important coal fields worldwide and extends over an area of 60 000 km 2 (Sachsenhofer et al., 2012).The selected coals cover a broad maturity interval of 0.62-1.47%Rr (vitrinite reflectance) and were intensively investigated regarding their depositional environment, paleoecology, geochemistry and petrography dur- ing previous studies (e.g.Stovba and Stephenson, 1999;Stephenson et al., 2001;Sachsenhofer et al., 2003;Izart et al., 2006;Misch et al., 2016).Hence, they offer a perfect opportunity to evaluate multiple influencing factors on coal stability and to develop a nanoindentation testing routine for inhomogeneous, organic matter-rich samples.

Geological Setting
The Donets Basin (Fig. 1) is a Late Devonian rift structure and forms the south-eastern part of the Pripyat-Dniepr-Donets-Karpinsky Basin located within the eastern European craton (Stovba and Stephenson, 1999;Stephenson et al., 2001).A detailed chrono-and lithostratigraphic overview of the Donets Basin is given in Sachsenhofer et al. (2003) and Misch et al. (2016).Devonian pre-and syn-rift rocks with a thickness of 750 m form the margin of the basin, which reaches depths up to 5 km at its centre.The basement consists of Precambrian crystalline rocks of the eastern European craton, which are overlain by middle Devonian to Carbonif-erous rocks (Chekunov et al., 1993;Maystrenko et al., 2003).The Carboniferous succession is up to 14 km thick and subdivided into lithostratigraphic units named suites A, B, C, to P (Lutugin and Stepanov, 1913;Fig. 2).The coal-bearing sequence is formed by fluvial sandstone, marine limestone or mudstone, and deltaic mudstone and siltstone (Sachsenhofer et al., 2003).Marine deposits are more common in the eastern than in the western part of the basin.In general, coal seams and intercalations of coal exist throughout the Carboniferous succession.The lower Serpukhovian (suite C) and Moscovian successions (suites K to M) are particularly rich in coal (Levenshtein et al., 1991; see Fig. 2).Serpukhovian coal seams are concentrated in a narrow NW-SE trending zone along the southwestern basin margin (Shulga, 1981;Fig. 1) where accumulation occurred during a regressivetransgressive cycle in a roughly 10 km wide shore-zone dissected by several NE trending rivers discharging into a nearby shallow sea in the central Donets Basin (Sachsenhofer et al., 2012).The high amounts of inertinite and liptinite, as well as usually very low ash yields and low to moder-  (Stovba and Stephenson, 1999).For the Moscovian sequence a detailed stratigraphic column is given (after Sachsenhofer et al., 2003Sachsenhofer et al., , 2012)).The studied seams are marked with star symbols.ate sulphur contents are characteristic for Serpukhovian coals (Sachsenhofer et al., 2003).In contrast, Moscovian seams are generally rich in vitrinite macerals and show high ash yields and variable sulphur contents, related to changing depositional systems with varying marine influence (Sachsenhofer et al., 2003).
The Carboniferous succession in the Donets Basin is of great economic importance; it hosts about 130 economic coal seams, each with a thickness between 0.45 and 2.5 m, summing up to a total net thickness of about 60 m (Ritenberg, 1972).Most coal fields of the Donets Basin contain high amounts of methane (from 8-10 to 25-37 m 3 t −1 ; Sachsenhofer et al., 2003), which, apart from mine safety issues, constitute a high inherent potential for coal bed methane production (Privalov et al., 2007

Coal Samples
The investigated coal samples were collected from wellstudied seam profiles in the (south-)western part of the Donets Basin (Table 1; sampling locations are shown in Fig. 1, for stratigraphic assignment see Fig. 2).Twelve coal samples from Moscovian (11) and Serpukhovian (1) seams were selected based on their vitrinite reflectance (VR) obtained by Sachsenhofer et al. (2003).The selected coals cover a maturity range from 0.62 to 1.47 %Rr (Table 1).The coal specimen were embedded in epoxy resin and prepared as polished sections (final polishing steps with 9, 3, 1 µm diamond polishing suspension and 0.05 µm micropolish powder).Individual macerals were identified in an optical microscope under white reflected light (no oil immersion), and colormarked to ensure the correct identification during the nanoindentation procedure.No oil immersion was used to enable color-marking.

Nanoindentation
The nanoindentation tests were performed using a Hysitron (now Bruker) Triboscope system with the transducer being mounted on the scanner head of a Digital Instruments (now Bruker) D3100 atomic force microscope, allowing for surface scans.A Cube Corner diamond tip with a tip radius of approximately 200 nm was used as indenter.Surface scans generated with the Cube Corner tip provided sufficient information to identify suitably flat areas without scratches or other disturbing topography features that might influence the obtained material parameters.The area function of the Cube Corner was calibrated using a fused silica standard following the Oliver-Pharr method (1992).All indentation tests were performed in load control mode with a trapezoidal loading profile to a maximum load of 1000 µN held for 10 s and loading/unloading rates of 500 µN s −1 .To reduce the influence of thermal instrumental drift, fast loading and unloading sequences were selected, while the holding time of 10 s at maximum load was chosen to minimize the influence of creep during unloading (Fig. 3).To exclude surface effects or other measurement bias (e.g.due to unintentional indentation of mineral grains embedded in organic matter), all results were quality-checked based on the obtained load-displacement curves that can be regarded as material-characteristic finger prints.Hardness (H ) and reduced elastic modulus (E r ) were determined using the method proposed by Oliver and Pharr (1992).H is defined as the resistance of a material against external force.It is calculated following Eq.( 1): where P max is the maximum load applied on the projected contact area (A) between indenter and material.A is calculated from the tip area function and the indentation depth (h) into the sample surface at maximal load.E r is obtained by evaluation of the load-displacement curves according to Eq. ( 2).Therein, the contact stiffness S describes the slope of a tangent that is fitted to the elastic unloading curve segment.
During the indentation procedure the diamond tip is also elastically deformed, as it is not entirely rigid.This effect is considered trough Eq. ( 3): where E is the Young's modulus of the sample and ν denominates its Poisson's ratio.E i and ν i are the known material properties of the indenter tip.However, for soft materials E and E r are almost equal, as the deformation of the indenter tip can be neglected in this case.Therefore, and because the Poisson's ratios of macerals are not known, this study utilizes E r as material parameter.

Inertinite reflectance measurements
Three coal samples with different ranks (Table 1) were selected for inertinite reflectance measurements.The measurement of random inertinite reflectance in oil (IR) was performed using a magnification of 100× in non-polarized light at a wavelength of 546 nm (Taylor et al., 1998).Inertinite reflectance measurements were compared to nanoindentation results to investigate the relationship between IR, H and E r for various inertinite sub-macerals.4 Results

Vitrinite
The indentation results indicate that the obtained average H and E r values of vitrinites do not follow a clear maturity trend (Fig. 4; Table 2).The lowest H and E r averages are observed in the low rank coals (0.62 %Rr), with values of 0.43 ± 0.03 and 4.91 ± 0.37 GPa, respectively.Both parameters vary considerably at 0.62 %Rr.Medium rank coals (0.69-0.76 %Rr) show the highest H and E r averages.The obtained average H and E r values at 1.47 %Rr are 0.48 ± 0.01 and 5.06 ± 0.2 GPa, respectively, and hence comparatively lower compared to samples of medium rank.Generally, vitrinites show elasto-plastic behavior at all maturity levels (Fig. 5).

Inertinite
H and E r values obtained from inertinites do not follow a maturity trend (Table 3).However, a linear correlation between the obtained material parameters and measured inertinite reflectance (1.67-3.76%Rr) is observed (Fig. 6).
IR is mainly dependent on the investigated sub-maceral type.Semifusinite shows the lowest reflectance values, while macrinite shows intermediate reflectance.Inertodetrinite and pyrofusinite exhibit the highest reflectance values.Variable shapes of the load-displacement curves illustrate that the deformation behaviour changes from elasto-plastic to mainly elastic with increasing IR (Fig. 7).

Liptinite
Sporinite is the most abundant liptinite maceral in the selected coals and hence indentation tests were mainly conducted on sporinites.The lowest value for H was determined for the least mature sample (0.62 %Rr), whereas the highest values for both H and E r occur at 0.71 %Rr.After the peak, a decrease in the values of E r and (to a lesser degree) H is observed.E r and H averages at 0.83 and 0.92 %Rr are comparably low (Table 4; Fig. 8).At 0.93 %Rr, the obtained E r and (to a lesser degree) H values show a multimodal frequency distribution (Fig. 9).The material behavior of liptinites can be described as elasto-plastic with dominance of the plastic deformation part at all maturity-levels (Fig. 10).Loaddisplacement curves of liptinites exhibit the largest amount of energy dissipation and highest indentation depths of all maceral groups.

Vitrinite
In general, the coalification process of vitrinites at the bituminous rank is characterized by various changing parameters, including loss of hydrogen-and oxygen-bearing functional groups and volatile matter, increasing aromaticity, decreasing primary porosity and increasing density, which leads to generally increasing organic carbon contents and increased homogenization (e.g.Taylor et al., 1998;Tissot and Welte, 1984).However, the density of vitrinites does not change linearly with ongoing thermal maturation.At the boundary between sub-bituminous to bituminous coal rank, the true density decreases with rank as a result of the devolatilization of mainly heavy oxygen-rich compounds.The density increases again at advanced maturity because hydrogen-rich components are transformed during hydrocarbon generation and the increase of aromatic units leads to a tighter packing of the polymer structure (Taylor et al., 1998;Van Krevelen, 1961).Nevertheless, it might decrease again due to the generation of secondary (nano-)porosity.ties of vitrinite macerals investigated during this study reflect the complex effects described above, rather than a linear maturity trend (Fig. 4a, b).Samples at low to medium rank show largely varying mechanical properties, whereas less scatter-ing results at advanced maturity (> 1.0 %Rr) point to increasing homogeneity of the organic matter structure at this stage, reflected also by a more uniform geochemical composition (e.g.Misch et al., 2016).Apart from maturation-related processes such as devolatilization (decreases primary porosity) and hydrocarbon generation (increases secondary porosity), changing mineral matter contents, best approximated by the ash yield, might affect porosity and relative density of vitrinites especially at an early mature stage.The obtained E r and (to a lesser degree) H values correlate with ash yields (Fig. 11a, b; ash yield data from Sachsenhofer et al., 2003) at low to medium rank.Hence, the depositional setting can be considered an equally important influencing factor for the mechanical properties of vitrinites at this maturity stage.

Inertinite
The highest average values for H and E r were obtained from inertinites.Apparently, the mechanical properties of inertinites are not primarily controlled by the coal rank but correlate well with the measured inertinite reflectance.Re- flectance studies of recent wildfire charcoals and experimental pyrolysis studies on charcoalified wood prove that inertinite reflectance increases with burning temperature, regardless of its morphology (Scott and Glasspool, 2007).A wildfire might produce a spectrum of different temperatures up to 1000 • C depending on the grade of exposition (Pyne et al., 1996).Zickler et al. (2006) studied the influence of artificial pyrolysis treatments on the mechanical properties of wood cells via nanoindendation.Between 500 and 900 • C they described a continuous rise of E r .A continuous increase could also be observed in H at the pyrolysis temperature between 220 and 700 • C. Furthermore, a transition from elasto-plastic to mainly elastic material behaviour with rising temperature according to the load-displacement curves was described by the authors.These findings are reflected by the clear correlation of mechanical properties and deformation behaviour with the measured inertinite reflectance of the investigated samples.Zickler et al. (2006) interpreted this trend as an increase of covalent carbon-bonding and a glass-like structure characterized by randomly cross-linked aromatic carbon sheets.Hence, the main controlling factor for the material behaviour of inertinite is the temperature exposure during paleo-wildfire events.

Liptinite
The indentation results of liptinites indicate a non-linear relationship with maturity.Both H and E r of sporinites increase with maturity from low to medium rank (up to ∼ 0.7 %Rr; Fig. 8a, b).The sample at 0.71 %Rr exhibits a maximum of both parameters, likely related to progressing organization, dehydration and decarboxylation processes of the biopolymers, predominantly occurring at low to medium rank.A gradual decrease of both E r and H occurs from 0.76 to 0.83 %Rr, and the lowest E r value occurs at 0.93 %Rr (Fig. 8a).The sharp decrease in measured reduced elas-tic moduli might indicate the onset of thermal cracking related to hydrocarbon generation, associated with the loss of aliphatic chains and the formation of isolated aromatic rings (Yule et al., 2000).The multimodal distribution of measured E r and H (Fig. 9a, b) values points to an anisotropic behaviour that is caused by structural reorganisation and polymerisation of the aromatic units that leads to the formation of sheets (Yule et al., 2000).The outlier of E r at 0.69 %Rr is likely caused by changes in the depositional environment, as Misch et al. (2016) reported elevated organic S contents at the expense of C, related to an increased marine/brackish influence on the Novogrodovskaya seam during peat formation that agrees with findings from Sachsenhofer et al. (2003).
Various authors (e.g.Orr, 1986;Tomić et al., 1995) postulated earlier hydrocarbon generation for coals high in organic sulphur.Therefore, the low E r and H averages determined for the Novogrodovskaya sample can be interpreted as earlier onset of transformation reactions related to hydrocarbon generation.

Conclusion
The present results indicate that apart from maturity, multiple parameters control the mechanical properties of macerals.
In case of inertinites, the prevailing conditions (temperature, exposure) during paleo-wildfires represent the main controlling factor for the mechanical properties.E r and H of liptinites are strongly influenced by transformational processes related to devolatilization at low to medium rank, and later hydrocarbon generation.Furthermore, variations in the depositional environment (e.g.increased organic sulphur content due to brackish conditions) might play an additional role.The mechanical properties of vitrinites show no clear correlation with maturity, as a result of multiple influencing factors such as mineral matter inclusions, depositional conditions, and transformation of the organic matter structure and possibly www.adv-geosci.net/45/73/2018/Adv.Geosci., 45, 73-83, 2018 the generation of nanoporosity at advanced maturity.In order to achieve a deeper understanding of the interdependency between mechanical properties and microstructural features (e.g.nanoporosity, mineral inclusions), advanced characterization techniques such as transmission electron microscopy (TEM) are required, which will be a focus of future research.

Figure 1 .
Figure 1.(a) Geological map of the Donets Basin (after Sachsenhofer et al., 2003).(b) Detail map showing an overview of the sampled coal mines.Locations of samples investigated during this study are marked in bold.(c) Regional overview of the Donets Basin in Europe.

Figure 3 .
Figure 3. Example of a load-displacement curve, its trapezoidal loading profile and a 3-D-view of the resulting indent.

Figure 5 .
Figure 5. Averaged load-displacement curves for vitrinites at different coal rank.Slightly differing deformation behaviour is indicated by the individual curve shapes.However, a general elasto-plastic behaviour dominates at all maturity ranges.The insets show micrographs of typical indentation spots.

Figure 7 .
Figure 7. Averaged load-displacement curves for inertinites at different coal ranks.The deformation behaviour changes rapidly with increasing IR, as indicated by the dashed arrow.Insets show different investigated inertinite sub-macerals (SF = semifusinite; ID = inertodetrinite).

Figure 9 .
Figure 9.The frequency distribution of the obtained E r (a) and H (b) values at 0.93 %Rr.

Figure 11 .
Figure 11.Correlation between ash yield (data from Sachsenhofer et al., 2003), E r (a) and H (b) of vitrinite macerals.Note that samples of equal maturity (and hence chemical composition) are displayed for comparison.(CI = Confidence interval). ).

Table 2 .
Obtained H and E r averages for vitrinite, including standard deviations (SD) calculated from the respective number of measurements (n).

Table 3 .
Sachsenhofer et al. (2003)ned for various inertinite sub-macerals, including standard deviations (SD) calculated from the respective number of measurements (n).Vitrinite reflectance (VR), indicative for the thermal maturity of the whole coal sample, was adopted fromSachsenhofer et al. (2003).Inertinite reflectance values (IR) of investigated inertinite macerals were measured during this study and reflect the conditions during pre-depositional wildfires, rather than thermal maturity related to burial.Note that inertinite reflectance refers to the actually indented macerals.

Table 4 .
H and E r averages obtained from sporinites, including standard deviations (SD) calculated from the respective number of measurements (n).