Supercritical methane adsorption and storage in pores in shales and isolated kerogens

Shale gas is an important hydrocarbon resource in a global context. It has had a significant impact on energy resources in the US, but the worldwide development of this methane resource requires further research to increase the understanding of the relationship of shale structural characteristics to methane storage capacity. In this study a range of gas adsorption, microscopic, mercury injection capillary pressure porosimetry and pycnometry techniques were used to characterize the full range of porosity in a series of shales of different thermal maturity. Supercritical methane adsorption methods for shale under conditions which simulate geological conditions (up to 473 K and 15 MPa) were developed. These methods were used to measure the methane adsorption isotherms of Posidonia shales where the kerogen maturity ranged from immature, through oil window, to gas window. Subcritical methane and carbon dioxide adsorption studies were used for determining pore structure characteristics of the shales. Mercury injection capillary pressure porosimetry was used to characterize the meso and macro porosity of shales. The sum of the CO 2 sorption pore volume at 195 K and mercury injection capillary pressure pore volumes (1093–5.6 nm) were equal to the corresponding total pore volume (< 1093 nm) thereby giving an equation accounting for virtually all the available shale porosity. These measurements allowed quantification of all the available porosity in shales and were used for estimating the contributions of methane stored as ‘free’ compressed gas and as adsorbed gas to overall methane storage capacity of shales. Both the mineral and kerogen components of shale were studied by comparing shale and the corresponding isolated kerogens so that the relative contributions of these components could be assessed. The results show that the methane adsorption characteristics were much higher for the kerogens and represented 35–60% of the total adsorption capacity for the shales used in this study, which had total organic contents in range 5.8–10.9 wt%. Microscopy studies revealed that the pore systems in clay-rich, organic-rich and microfossil-rich parts of shale are very different, and also the importance of the inter-granular organic-mineral interface.


Introduction
The development of shale gas as a methane resource started in the US and it has become an increasingly important component of natural gas supply in the past twenty years [1]. There are large potential shale gas reserves in the rest of the world, which could be exploited in the future. Compared to coal, shale only contains relatively small amounts of methane, but is economic because of very large deposits. The economics of shale gas reservoirs depends on the Gas-in-Place (GIP) and the methane extraction rate that can be achieved [2]. This depends on the properties of the shale and the fracturing of the shale matrix induced by the hydraulic fracking and horizontal drilling processes.
Shales are complex heterogeneous materials comprising of inorganic mineral and organic kerogen components, which leads to wide variations in pore structures and surface chemistry. Methane is stored as either compressed gas or physisorbed gas in porous structures in organic-rich shales and the quantities of these components are very difficult to predict. The gas filled fraction of the porosity follows known relations with pressure and temperature, but it is difficult to predict the quantity. In comparison, the adsorbed gas contribution is poorly understood, but may amount to 50-60% of the gas in some shales [3]. Shale porous structure characteristics that are particularly relevant for assessing shale reservoir potential are pore volumes, pore size distribution and adsorption characteristics under simulated geological conditions. The economic potential of the resources needs to be established unequivocally, in order to justify extraction of shale gas.
In conventional reservoirs, gas is stored in large pores as 'free' , high pressure gas. There is very little gas adsorbed on the walls of macropores compared to free gas in the pores. In contrast, shale is a complex heterogeneous material typically consisting of > 90% minerals with < 10% organic kerogen material, in which almost all of the pores are less than a micrometre in diameter. Methane is generated from the decomposition of the organic kerogen phase in shale during geological processes over the temperature range 370-550 K. A gas rich phase is generated at temperatures > 430 K [4]. Methane is stored in pores, which range from the micropores (< 2 nm), through mesopores (2-50 nm) to larger macropores, with the adsorbed phase mainly in the micropores in equilibrium with compressed gas in larger meso and macropores through an interconnected pore structure. The adsorbed gas is slow to be released and could be an issue for shale gas exploitation. There have been only a limited number of studies of methane adsorption on shale because of the difficulty in measuring isotherms under supercritical conditions [5][6][7][8][9][10][11][12][13][14]. At the outset, measurements were limited to 338 K, but recently these have been extended to higher temperatures [15,16]. The amounts of methane in the adsorbed and compressed gas phases need to be quantified in order to understand the delivery mechanism of shale gas to the production well.
The challenge was to quantify the shale porosity, the methane stored in shale as 'free' compressed gas and the adsorbed gas phase, and the relative proportions of methane stored in both the organic and inorganic phases. Since shales contain pore structures ranging from micro to macro pores, a number of techniques must be used to quantify the different pore size ranges. These methods include high pressure supercritical gas adsorption to determine the methane adsorbed phase under simulated geological conditions, mercury injection capillary pressure (MICP) porosimetry for characterizing meso and macro pore characteristics, subcritical gas adsorption for characterizing micro and meso pore structure, helium and grain density measurements to measure the density of pore wall material and microscopy to provide ideas on the origin of kerogen and mineral particle interfaces and pore connectivity. Small-angle neutron scattering techniques may also be used to provide information on shale structure and porosity [17][18][19].
The porous structures of shales are affected by kerogen maturity [20], mineralogy [21], grain size [22], diagenesis and compaction [23][24][25]. There have been a number of studies of the relationship between total organic content (TOC) and adsorption characteristics for dry organic rich shales [5,7,8,14,26]. Gasparik et al. studied the effect of TOC on methane adsorption capacity for two suites of dry shales with a range of sources and maturity: 1) 6 shales [9] with TOC 1.0-10.5% did not show a correlation and 2) 28 shales [14] (vitrinite reflectance (V R ) = 0.7-2.4%) with TOC (0.4-14.1%) gave a positive correlation, but with some samples deviating from the trend. Wang et al. found [27] a positive correlation between methane sorption capacity and TOC for 10 Longmaxi shales (V R = 2.7-3.68 and TOC = 1.13-3.59%). Liu et al. [28] also investigated 28 Longmaxi shales (TOC 0.25-4.45%) and found a positive correlation for methane sorption at 303 and 363 K, with TOC. The complexity of the problem of the relation between methane sorption capacity and TOC over a wide range of shale sources with varying maturities and mineral contents is lack of knowledge of the methane fraction stored in the mineral and organic structures and also, the kerogen structure changes with maturity. An additional factor that needs to be considered is the presence of moisture in shale, which has a detrimental effect on methane sorption [29,30]. Clay mineralogy in the shale is an indicator of potential for water adsorption [29,31]. The decline in methane sorption capacity with increasing water content has been attributed to competitive adsorption, water cluster formation obstructing access to and occupying pore space. Therefore, methane adsorption capacities for dry shales can be regarded as maximum values whereas moisture equilibrated (relative humidity 97%) methane sorption capacities can be considered as minimum values [29]. Estimation of organic-rich shale potential for utilization as a gas resource presents major problems and challenges related to the quantification of the amounts of compressed and adsorbed gas and how this is related to shale porous structural characteristics. The aims of the study were: 1. to establish methods for the accurate measurements of supercritical methane adsorption capacity of dry shales at high temperatures and pressures, which simulate geological conditions (up to 15 MPa and 473 K); 2. to account for all the porosity (micro-, meso-and macro-) in order to estimate the amounts of sorbed gas and 'free' compressed gas stored in shale pores; 3. to determine the relative contributions of the organic kerogen and inorganic mineral phases to overall methane capacity; 4. to establish the role of kerogen content and maturity on methane capacity and 5. to integrate the porous structure characterization data with microscopic data to gain an insight into mineral and kerogen structural interfaces.

Shale samples
The sample of shale Alum 1 was obtained from the Skelbro-2 well in Bornholm, Denmark [32]. The Lower Toarcian Posidonia shales were obtained from the Wickensen (WIC), Harderode (HAR), and Haddessen (HAD) boreholes in Germany (see Fig. 1) [33]. These samples were core samples and the surface showed no significant evidence of contamination. The samples used were taken from the centre of the cores to avoid any possibility of surface contamination. The samples were dried at 105°, but no further treatment of the samples was considered necessary. These sampling and preparation methods were used to ensure that the samples were as far as possible representative of those in the geological formations. Two samples from each borehole were investigated in detail. The maturity of the shale samples increased progressively from the early oil-window Wickensen (V R = 0.53%), through the mid oil window Harderode (V R = 0.89%), to the gas-window Haddessen (V R = 1.45%) samples. Characterization data for the samples are given in Tables 1 and  2. The mineralogical compositions of shales in this sample series were similar over the range of maturity with the following ranges: calcite (31-55 wt%), phyllosilicates (23-37 wt%), quartz (8-16 wt%), pyrites (4-9 wt%), dolomite (0.3-6.4%), feldspars (1-5 wt%). TOC values were in the range 5.8-10.92 wt% [34]. The shale and kerogen sample set allowed the impact of kerogen maturity on shale and kerogen pore structure on methane capacity to be studied independently of TOC and mineral composition variables.

Kerogen samples
The kerogen samples were separated from the shales. Initially the shales were treated with HCl to dissolve carbonates and then treated with acidified CrCl 2 to remove pyrite. Silicates from this sample were removed by treatment with HF [34]. Kerogen purity was checked by powder X-ray diffraction [34].

Grain density
Samples (particle size range: 500-1180 μm) were pre-dried overnight at 105 °C in air. The sample was weighed in a pre-weighed pycnometer (50 mL). 10 mL of Teepol ® soap solution (concentration: 5%) was added to the pycnometer, which was then filled with degassed water and left overnight to equilibrate. The pycnometer + sample + water was weighed at 25 °C. The pycnometer was filled with Research Article SN Applied Sciences (2020) 2:780 | https://doi.org/10.1007/s42452-020-2517-6 de-gassed water and weighed at 25 °C to determine the volume of the pycnometer. The particle grain density was calculated from the above weights and the known density of water at 25 °C.

Helium density
Helium density measurements were carried out as part of the high pressure adsorption measurements since the information is necessary for calculating the surface excess values. Pre-dried kerogen and shale samples were loaded into a Hiden Isochema Intelligent Manometric Instrument (IMI). The samples were outgassed typically for ~ 4 h), at < 10 −4 Pa, at 110 °C in the IMI instrument before helium pycnometry. Helium dosing at a pressure of 2 MPa at 40 °C from a reference cell was used with successive evacuation and helium dosing cycles with the skeletal volume and sample density being determined each time. The mean of 3 repeats was used. A stainless steel cylinder was used as a volume calibration standard under matching conditions. The volume of gas injected into the sample cell was calculated using the NIST Standard Reference database 23 by using the REFPROP Version 9.0 software [37]. Helium density measurements were repeatable to ± 0.01 g cm −3 .

Low-pressure adsorption
An Intelligent Gravimetric Analyzer (IGA), supplied by Hiden Isochema Ltd., Warrington, UK; was used for low temperature nitrogen, carbon dioxide and methane adsorption studies of shales and kerogens. Samples of shales (500-1180 μm) and kerogens were out-gassed under ultra-high vacuum at 383 K until constant weight. Liquid nitrogen (77 K) and solid carbon dioxide/acetone cryogenic baths (195 K) were used for temperature control. Methane adsorption measurements at 112 K were  achieved using a Hiden Isochema cryostat system cooled by gaseous nitrogen generated from liquid nitrogen. A computer controlled recirculating bath of water/ethylene glycol was used for adsorption measurements at 273 K. Isotherm repeatabilities were typically ± 1.5% for CO 2 adsorption at both 196 and 273 K at 100 kPa and ± 1.1% for N 2 uptakes at 77 K and 99 kPa.

High-pressure supercritical methane sorption
A Hiden Isochema Intelligent Manometric Instrument (IMI) was used for high pressure CH 4 adsorption. Shales and kerogen samples were pre-dried at 383 K in a vacuum oven and loaded typically with ~ 10 g shale (particle size range: 500-1180 μm), 0.8-1.3 g kerogen (particle size range as produced from demineralization process) into the IMI sample cell. Displacers were employed to reduce the void volume, because only relatively small quantities of kerogen samples were available. The skeletal volumes of the samples were determined, which also allowed calculation of the corresponding helium densities. Methane isotherms were measured at 318, 338 and 358 K. Isotherm experimental repeatabilities were typically ± 5.0% for both shale and kerogens at 10 MPa. Helium isotherms were measured at all temperatures as blank determinations for no adsorption to monitor skeletal densities and leak-test the instrument.

Mercury injection capillary pressure (MICP) porosimetry
A Micromeritics Autopore IV Mercury Injection Porosimeter was used in the study [34]. MICP could only be used for the shales since kerogens only occur as very small particles in shale. Also, kerogens are plastic and deform under the very high pressures used in the technique. Initially shale samples were freeze-dried for 48 h. Samples with a volume of ~ 1 cm −3 were then outgassed under vacuum prior to MICP porosimetry studies. The mercury pressure was increased stepwise from 1.379 up to 268.9 MPa in order to characterize the meso/macroporosity. MICP bulk volumes/densities of shale were calculated from the known sample mass and the grain or helium density of the sample. The difference of the volume of mercury injected at 1.379 MPa and 268.9 MPa was used to calculate MICP pore volumes. Surface topography and micro fractures related to the de-stressing and drying of geological samples were excluded (see Table 1). The Washburn equation [41,42] is derived for cylindrical pores and calculations were carried out assuming a contact angle 141° between mercury and the particle surface, and a surface tension of 0.485 N m −1 [24,43,44]. The mercury pressures were used to predict pore throats (constrictions) with equivalent diameters between 1093 nm at 1.379 MPa and 5.6 nm at 268.9 MPa. These diameters should be regarded as equivalent pore diameters because of the variation in pore shape in the heterogeneous shale materials.

Microscopy studies
Scanning electron microscopy (SEM) was carried out using a Hitachi SU-70 High Resolution Analytical instrument equipped with an Oxford Instrument Energy Dispersive X-ray (EDX) microanalysis system (INCA Energy 700). Samples were prepared as carbon-coated polished thin sections. Back Scattered Electron (BSE) mode used a YAG detector. Some samples (3 mm diameter) were polished with an argon broad ion beam (BIB) in a GATAN 691 Precision Ion Polishing System (PIPS TM ) to reduce shale topography. Secondary Electron (SE) images were obtained using a through-the-lens detector (TLD) at magnifications of up to 10,000×. EDX analysis was used to generate elemental and phase distributions of selected areas.

Supercritical methane adsorption on shale
Initially it was necessary to provide information on supercritical methane adsorption characteristics under simulated geological conditions in order to gain an insight into the variation of sorption capacity with temperature and pressure.

Isotherms
At the start of the study it was apparent that the accuracy of supercritical gas adsorption measurements at high pressures was not established. Subsequently, an investigation of the accuracy of high pressure supercritical CO 2 adsorption measurements on highly porous activated carbons was reported [45]. Similar studies have been reported recently for CO 2 adsorption on zeolite ZSM-5 [46]. Initial validation of isotherm reproducibility successfully used a small sample (~ 500 mg) of F400 activated carbon. The development and validation of high pressure adsorption techniques for measuring supercritical adsorption isotherms for shales is much more challenging than for activated carbon since the adsorption is two orders of Research Article SN Applied Sciences (2020) 2:780 | https://doi.org/10.1007/s42452-020-2517-6 magnitude lower. This led to the use of larger sample amounts (typically 10 g of shale and ~ 1 g of kerogen) being used to obtain the required accuracy. These techniques were used to investigate the variation of isotherms under a range of simulated geological pressure and temperature conditions and validate the consistency of the data. Figure 2 shows methane adsorption isotherms for shale Alum 1 over the temperature range 299.6-473.2 K and pressures up to ~ 140 bar [35]. These conditions cover the temperature and pressure range that simulate geological conditions. The reproducibility of the measurements were confirmed in a 'round robin' study of an immature Posidonia shale and an over-mature Upper Chokier formation shale [47]. The Newcastle results were in the middle of the set of results confirming reasonable agreement. The surface excess isotherm can be converted to the isotherm for the absolute amount using the equation below where n abs is the absolute amount adsorbed, n ex is the surface excess determined experimentally, ρ b is the bulk gas phase density and V ad is the adsorption pore volume. This equation has been used for crystalline materials, where the structure has been determined from crystallographic studies [48][49][50][51]. There is good agreement between crystallographic and adsorption pore volumes. However, shales are heterogeneous materials with a wide pore size distribution and not all larger pores have significant gas adsorption. The issue for shales is to identify the pore size ranges and volumes associated with both adsorption and high pressure compressed gas, which are the phases in equilibrium for methane storage in shale.
(1) n abs = n ex + b * V ad 3.1.1.1 Isosteric enthalpy of adsorption The isotherm data were used to calculate the isosteric enthalpy of adsorption using the Clausius-Clapeyron equation and to validate the consistency of the data [35]. The derivation of this equation assumes that the gas is ideal and the adsorbed phase volume is neglected. These two assumptions may lead to errors at high pressures [52], which can be examined using graphs ln(Fugacity) versus 1/T at various surface coverages for the set of data shown in Fig. 3. The linearity of the graphs demonstrate the consistency of the data over simulated geological conditions.
The determination of isosteric enthalpy of adsorption (Q st ) graph is limited by the availability of data at constant surface coverage. The Clausius-Clapeyron graph for the lowest 5 temperatures and 2.85-7.62 MPa and ~ 70% of maximum surface excess (0.125 mmol g −1 ) gives Q st of 20.5 ± 1.3 kJ mol −1 . This compares with a Q st of 19.0 ± 0.2 kJ mol −1 for a surface excess of 0.01 mmol g −1 .
Graphs of ln(Pressure) versus 1/T gave very similar values for Q st . These values are typical for adsorption of CH 4 on shales [8], activated carbon [53][54][55], coal [56] and metal organic framework materials [57]. There is no evidence for any significant variation in isosteric enthalpy of adsorption even though shale is heterogeneous with a variety of inorganic mineral and organic kerogen components.

Isobars
The variation of CH 4 surface excess with temperature at constant pressure is important in assessing shale adsorption characteristics under simulated geological conditions. The adsorption isotherm data in Fig. 2   adsorption on temperature can be determined. Figure 4 shows a comparison of surface excess (SE) isobars for a phenomenological model [SE = k/T where k is a constant and T is temperature (K)] and the equation proposed by Illin (SE = k 1 exp(− k 2 T 0.5 )) where k 1 and k 2 are constants and T is the temperature (K) [58,59]. Comparison of the model over the pressure range 0.5-13.5 MPa shows that both equations fit the data, but the simple phenomenological model provides the slightly better fit to the data. Previous studies of supercritical methane adsorption on a metal organic framework have shown that these isobar equations provide good descriptions of the data [57].

Subcritical CH 4 adsorption studies
Since physisorption increases with decreasing temperature, assessment of the maximum amount of methane that can be adsorbed on a shale was obtained using low temperature subcritical adsorption isotherm. A comparison of the subcritical isotherms for nitrogen, methane and carbon dioxide on a relative pressure basis for the Alum 1 shale are shown in Fig. 5. The total sorption pore volume obtained from subcritical CH 4 adsorption at 112 K represents the maximum pore volume available for methane adsorption since the low temperature allows measurements with a relative pressure close to 1, where condensation occurs. There is an experimental difficulty in measuring subcritical CH 4 adsorption because of the lack of a suitable cryogenic fluid. The measurement was carried out using a special cell with circulation of cooled gaseous nitrogen, but this is experimentally challenging. The Dubinin-Radushkevich (DR) ultramicropore volume obtained for the Alum 1 shale from CO 2 adsorption at 273 K was 12.9 mm 3 g −1 . The sorption pore volumes for subcritical CO 2 at 195 K (16.8 mm 3 g −1 ) and CH 4 at 112 K (18.0 mm 3 g −1 ) are similar as expected from the Gurvitch Rule [60]. This rule states that the adsorption uptakes at relative pressure ≈ 1 on a given adsorbent when expressed as a volume of liquid, using the liquid density, are similar for all adsorptives. Therefore, the CO 2 sorption pore volume was used as a substitute for the CH 4 total pore volume in subsequent studies of shale maturity, since it is much easier to use solid CO 2 /acetone as the cryogenic liquid [34,35]. Therefore, when studying the effect of shale and kerogen maturity on supercritical methane adsorption, correlations between subcritical CO 2 sorption pore volume and high pressure supercritical methane adsorption were examined to validate this approach.

Porosity characterization
Porosity characterization methods are well developed with the exception of the measurements of supercritical methane adsorption capacities for shales and kerogens  under a range of conditions, which simulate geological conditions of high pressure and temperature. The methane adsorbed phase is expected to be almost exclusively stored in the smallest pores of shales and kerogens. The porosity in which compressed gas is stored in equilibrium with the adsorbed phase mainly involves the larger pores.
The strategy for the characterization of shale porosity was to quantify all the porosity from macro-through mesoto micropores to complement the supercritical methane adsorption studies. A range of microscopy, gas adsorption, mercury injection and pycnometry techniques were used in order to explore the full pore size range in shales so that the relative amounts of adsorbed and compressed CH 4 could be established. Small-angle neutron scattering techniques may also be used to characterize shale structure and porosity [17][18][19] but these techniques were not available for this study. The measurements techniques available and size ranges covered are shown in Fig. 6. Mercury injection capillary pressure (MICP) porosimetry is an established technique for assessing porosity in minerals and it was important in providing measurements of macro and meso (> 5 nm) porosity in shales. MICP porosimetry measures pore throats (smaller passages or constrictions in the pore structure), which are connections between larger pore bodies. The pore throats are obvious restrictions in the pore structure and are typically a fraction of the pore size. The Washburn equation [42] was used to calculate the equivalent pore diameters for mercury entering pores at specific pressures [24,41,43,44]. At a pressure of 1.379 MPa mercury enters pores of > 1093 nm diameter and this allows the shale bulk density (ρ HgBulk ) to be determined. Therefore 1/ρ HgBulk is the volume of shale and pores < 1093 nm. Helium is not adsorbed to any significant extent under ambient conditions and enters all accessible pores. Therefore, the volume of shale is 1/ρ He , where ρ He is the helium density measured by pycnometry [34]. As expected grain density and helium density have a strong linear correlation and therefore ρ Grain can be used in place of ρ He . Therefore, the shale total pore volume (TPV) for equivalent pore diameters < 1093 nm and Total Porosity are given by the equations below The Washburn equation indicates that a capillary pressure of 1.379 MPa fills pore throats up to 1093 nm and 268.9 MPa fills up to 5.6 nm. Therefore, the MICP pore volume provides a measure of the macro + mesopore volume between equivalent pore diameters of 1093 and 5.6 nm from the amounts injected at pressure of 1.379 and 268.9 MPa.
Subcritical gas adsorption studies are routinely used for characterizing micro and mesoporous structures of porous materials using a variety of isotherm equations [61]. Adsorption occurs in micropores (< 2 nm) and some mesopores (> 2 nm). The Brunauer Emmett and Teller (BET) equation [62] accounts for multilayer adsorption and is usually applied to N 2 adsorption isotherm data at 77 K. N 2 adsorption isotherm data at 77 K has also been used to determine pore size distributions in order to gain insights into organic matter in coal and shales [63,64]. However, characterization parameters such as surface areas, obtained from N 2 adsorption at 77 K may have issues due to activated diffusion effects in ultramicroporous materials at very low temperatures [65][66][67] and may overestimate surface areas in microporous materials. In this case CO 2 adsorption is commonly used to characterize activated carbons [61,68] because of the higher adsorption temperature. CO 2 adsorption on shale is also of interest in relation to storage of CO 2 in depleted shales for sequestration applications [64]. The Dubinin-Radushkevich (DR) equation [69] is often applied to CO 2 isotherms at 273 K to determine the ultramicropore volumes (< 0.7 nm) [70][71][72]. Adsorption at 195 K was used to provide total CO 2 sorption pore volumes. Nitrogen adsorption at 77 K is used to provide BET surface areas and DR micropore volumes (0.3-2 nm) [72].
These techniques were used to characterize the porous structure in shales and kerogens (see Fig. 7) [34]. The CO 2 sorption pore volumes were obtained from adsorption isotherms at 195 K. Ultramicropore volumes were obtained using the Dubinin-Radushkevich equation from CO 2 adsorption at 273 K. The relationships of these parameters with supercritical methane adsorption is important in the context of accounting for all porosity present in shales by providing pore volumes in relation to pore size for shales. Since shales are heterogeneous materials, microscopic methods were used to visualize macro and meso porosity in order to understand the storage of methane in pores and flow properties of gases through pores. Microscopic methods showed that the suite of Posidonia shale samples did not have significant amounts of porosity above ~ 1 μm in diameter.

The role of shale structure and kerogen maturity
Supercritical CH 4 isotherms for Posidonia shales and the corresponding isolated kerogens at 338 K are shown in Fig. 8 [34]. The isotherms have a Type I shape in the IUPAC Classification scheme [73]. The adsorption isotherms for shales show that the surface excess goes through a minimum for the Harderode (oil window) shales. The adsorption isotherms are similar for Wickensen (immature) and Haddessen (gas window) shales. Similar trends were observed [34] for adsorption on shales at 318 and 358 K. The supercritical CH 4 isotherms for isolated kerogens show different trends (See Fig. 8b). The Haddessen (gas window) kerogens have higher surface excess uptakes than the Wickensen (immature) shales. The Harderode and Wickensen kerogens have similar surface excess isotherms. A similar trend was observed for the adsorption isotherms for isolated kerogens at 338 K. Comparison of the CO 2 adsorption isotherms at 195 K for the suite of shales and kerogens show that they have similar trends for the differences in the order of individual Wickensen, Harderode and Haddessen shales and kerogens as shown for supercritical methane adsorption described above. The CO 2 sorption pore volumes for kerogens (71.8-113 mm 3 g −1 ) are x 5-8 those of shales (8.3-16.4 mm 3 g −1 ). Thus the trends for shales and kerogens are confirmed. The high pressure shale and kerogen CH 4 isotherms can be parameterized and modelled using Langmuir and modified Langmuir equations [34,35].
Comparison of low temperature subcritical CH 4 and CO 2 adsorption studies at 112 K and 195 K at relative pressure = 1 showed the sorption pore volumes for both gases were very similar in agreement with Gurvitch's rule [35]. This represents the methane maximum sorption volume, which includes micropores (< 2 nm) and some mesopores. Research Article SN Applied Sciences (2020) 2:780 | https://doi.org/10.1007/s42452-020-2517-6 Figure 9 shows the variation of CH 4 surface excess at 11.5 MPa and 338 K with CO 2 sorption pore volume for both Posidonia shales and kerogens [34]. The correlation for shales alone had a regression coefficient (R 2 ) of 0.96 whereas there was more scatter for the kerogens with a R 2 value of 0.76 It is evident that there is a good overall linear correlation (R 2 = 0.956) for the combined shale and kerogen data showing the linear relationship passing through the origin between supercritical methane adsorption and CO 2 sorption pore volume [34]. Similar correlations were observed for CH 4 supercritical adsorption on shales and kerogens at 318 K. This observation supports the use of CO 2 sorption pore volume in the assessment of porosity available for high pressure CH 4 adsorption to account for the full range of shale porosity in order to assess the storage of compressed and adsorbed gas.
The interrelations between high-pressure CH 4 sorption capacity, maturity, and pore structure characteristics were investigated for shales and isolated kerogens. Mercury injection capillary pressure (MICP) porosimetry measurements could only be studied for larger particle sizes, which are only available for shales. Kerogens only occur as very small particles distributed throughout shales.
Microscopy studies showed that the shales only have sub-micrometre pores. The maximum pore throat size in the shale samples were also determined from MICP. The Washburn equation [42] provides the equivalent pore diameters as a function of injection pressures. The values for maximum pore throat size in the shales pass through a minimum (21 nm for HAR7038 and 24 nm for HAR7060) for the oil window shales. The largest maximum pore throat size was 547 nm in gas-window shale HAD7199. The maximum pore throat sizes for WIC7145, WIC7155 and HAD7090 were all very similar (~ 156 nm) [34]. The results show that the MICP range (< 1093 nm) studied covers the upper macropore range in this suite of Posidonia shale samples. Figure 10 shows the variation of Total Pore Volume (< 1093 nm) versus the sum of the CO 2 Sorption Pore Volume and MICP pore volume (5.6-1093 nm). It is apparent that there is a linear relationship (R 2 = 0.99) intersecting the origin, thereby accounting for all the porosity in the shales. This relationship allows predictions of the volumes available for adsorbed and compressed methane and hence stored amounts of adsorbed and compressed gas using an estimated density for the adsorbed phase and known equations of state for the gas phase. The variation of high pressure supercritical CH 4 surface excess adsorption at 11.5 MPa and 338 K for Posidonia shales and kerogens with subcritical CO 2 sorption pore volume [34] Vol.:(0123456789) SN Applied Sciences (2020) 2:780 | https://doi.org/10.1007/s42452-020-2517-6 Research Article The DR micropore, CO 2 sorption and MICP porosimetry total pore volumes for shales pass through a minimum in the oil window. Wickensen (immature) shales have CO 2 sorption pore volumes that are ~ 25% of the total pore volume. HAR (oil window) shales have CO 2 sorption pore volumes that are 46-66% of total pore volume and there is less porosity present as larger pores. This is due to compaction and bitumen filling. HAD (gas window) shales have CO 2 sorption pore volume that are 21-32% of total pore volume. The DR CO 2 micropore volumes as a percentage of the CO 2 sorption pore volume were 41-62% and 40-52% for the shales and kerogens, respectively. This suggests that there is ultramicroporosity in the mineral component of shale as well as the kerogen. In comparison, the DR micropore volumes of coals are in the range 0.014-0.052 cm 3 g −1 [79][80][81]. The trend observed for CO 2 DR micropore volume with coal maturity has a minimum for medium volatile bituminous coals. This is ascribed pore filling by low volatile hydrocarbon species followed by decomposition during coalification [82]. Figure 11 shows the variation of DR micropore volume, CO 2 sorption pore volume and Hg total pore volume with vitrinite reflectance for Posidonia shales and Fig. 12 shows the corresponding variation of DR micropore volume and CO 2 sorption pore volume for isolated kerogens. Comparison of Figs. 11 and 12 shows that the minimum in CO 2 sorption pore volume observed for the oil window shales is not observed for the corresponding kerogens indicating the significant contribution of porosity in the mineral phase of the shale. The clay content of the Posidonia shales is relatively small with illite/smectite and kaolinite in the range 18.5-26.2 wt% and 1.1-8.9 wt%, respectively. These clays adsorb CH 4 and are present in dry shales [76,83]. Since the kerogens adsorb 5-8 times more CH 4 than dry Posidonia shales, a correlation with TOC with some scatter associated with differences in maturity and mineral content would be expected (Fig. 12). Clays are present in the Posidonia shales and dry clays adsorb methane [76,83]. The isotherms observed for shale, clay content and kerogen and the composition of the shales do not always account for all the adsorbed CH 4 . This suggests that some porosity occurs in other mineral components and possibly occurs between organic and mineral matter interfaces [34]. This observation and the different trends for the series of shales and kerogen adsorption isotherms in Fig. 8, may explain the relatively large scatter in the points in the linear correlation between maximum shale capacity at 318 K and TOC shown in Fig. 13. Porosity distribution between the mineral and kerogen phases is important in the context of identifying correlations between methane storage and geochemical characterization data. Taking into consideration the TOC, 30-60% of the shale sorption pore volume is due to adsorption in the kerogen even though in these shales the TOC is 5.8-10.9% of the shale (see Fig. 14a) [34]. Since the CO 2 sorption pore volume has a linear correlation with supercritical CH 4 adsorption (Fig. 9) and the correlation in Eq. (4), the shale pore volume (5.6-1093 nm) determined using MICP porosimetry contains almost exclusively compressed gas which should be easily accessible. Figure 14b shows a typical comparison of supercritical shale adsorption of methane normalized to TOC with the kerogen methane isotherm [34]. These isotherms have a very similar shape suggesting that the kerogen in the shale has similar adsorption characteristics to the isolated kerogen indicating that the kerogen extraction process had little or no effect on the kerogen adsorption characteristics. Comparisons for other Posidonia shales and kerogens show the same similarity in isotherm shapes. This suggests that the kerogen adsorption properties are unchanged by the extraction/isolation process and implies that adsorption also takes place in the inorganic mineral phase. This is not surprising since clays are present. Some differences might also arise from accessibility within the rock matrix, but there is no direct evidence to support this proposition.
The research has been extended further to examine the microporous region by microscopic visualization and quantification of pores, which is significant in understanding the sorbed gas versus compressed gas capacity as well as the flow properties. Microscopy was used to study a wider range of 26 Posidonia shales [36]. Typical micrographs are shown in Fig. 15. In the case of immature shales, pores quantifiable by SEM (> ca. 50 nm diameter) only account for small part of total porosity. Visible pores are mainly associated with biogenic calcite. There was essentially no macroporosity in the clay matrix or kerogen components. Porosity loss in the oil window shales is due to a) compaction, b) carbonate cementation and possible swelling of kerogen by retained oil. In the gas window shales, porosity increases due to the formation of meso and macro pores in organic matter related to thermal cracking and gas generation. The macroporosity varies greatly and this reflects the different physicochemical characteristics of individual organic phase particles and the surrounding rigid mineral matrix preventing compaction of the more ductile organic phases. Helium and grain densities have a linear relationship and increase with increasing maturity for shales (see Fig. 16) [34]. The development of organic macroporosity cannot alone account for the porosity increase from oil to gas window. There is an accompanying change in the volume of micro and meso porosity. Intra-organic matter pores are a characteristic feature of gas window shales, but macroporosity is most commonly observed at the kerogenmineral interface (see Fig. 15). Broad Ion Beam (BIB) -SEM, MICP, and gas sorption methods individually only generate partial descriptions of shale pore systems. In combination, they generate a quantitative description of the size distribution and connectivity of pores ranging in size from nanometer to micrometer [36,84]. Pore size and volume data have been placed into a mineralogical clay-rich, microfossil-rich and organic matter and textural framework in order to be more predictive about the pore systems of shales generally. Clay-rich domains are both micro-and meso-porous (pore diameters < 40 nm) at all maturities. At temperatures below gas generation macropores are mainly 50-300 nm diameter. During oil generation carbonate-related macropores are traps for bitumen with porosity reduced due to bitumen retention and physicochemical compaction.
Gasification of organic matter increases porosity due to the loss of both solid phase and liquid phase organic matter. The porosity generated does not only occur within the organic matter but is physically associated with organic matter within carbonate porosity. Evolution of porosity with increasing maturity is thus linked to changes in both carbonate and organic matter domains. By quantifying the nature of pores in the various mineral and organic building blocks of shales we can consider the larger scale flow properties of shales and thus their potential either to transmit fluid (useful for a shale gas reservoir) or to retain fluid (useful for a CO 2 or nuclear waste storage site). These findings are consistent with a recent study [85]. Research Article SN Applied Sciences (2020) 2:780 | https://doi.org/10.1007/s42452-020-2517-6

Conclusions
The storage of methane in shale geological reservoirs is complex due to the heterogeneous nature of shale, which contains both mineral and amorphous organic kerogen components. Adsorption of methane under simulated geological conditions presented major challenges due to the low shale methane adsorption capacities and the wide range of temperatures (300-473 K) and pressures (< 150 bar). Therefore, detailed validation of the accuracy and reproducibility of the experimental methods was required as a prerequisite to studying other aspects of shale gas storage. High pressure methane sorption capacities and pore characteristics of bulk shales and isolated kerogens were investigated for immature, oil-window and gas-window samples from the lower Toarcian Posidonia formation. The equivalent maximum pore throat diameter for Posidonia shale samples were in the range 21-547 nm and optical microscopy did not reveal any visible porosity > 1 μm. However, only two samples at each maturity level have been studied. Total pore volumes obtained from MICP porosimetry and subcritical CO 2 sorption volumes (195 K) of organic-rich Posidonia shales decrease through the oil-window and then increase into the gas-window. The sum of porosities measured by CO 2 sorption at 195 K and mercury injection capillary pressure porosimetry [1.379-268.9 MPa (1093-5.6 nm)] were equal to the corresponding total pore volume (< 1093 nm) for shales thereby giving an equation that accounts for virtually all the available porosity in this series of shales. The mercury injection at 268.9 MPa occupies pores with constrictions larger than ca. 5.6 nm indicating that porosity measured by CO 2 adsorption at 195 K is largely within pores smaller than 5.6 nm. The CO 2 sorption pore volume represented 21-66% of the total pore volume in these shales, while the meso/macropore volume from MICP was 40-79% of the total pore volume. Around 40-50% of the CO 2 sorption pore volume was within the DR ultramicroporosity studies (< 0.7 nm) in the shales studied. Oil window shales have the smallest meso/macropore fraction and the largest CO 2 sorption pore volume fraction of the total porosity. It can be inferred that the adsorbed methane is predominantly in pores < 5.6 nm and the 'free' compressed gas is in the larger pores (5.6-1093 nm). The 5.6-1093 nm meso-macropore volume represents the major part in the Wickensen (78%) immature and Haddessen (63%) gas window shales, but only 41% of Harderode oil window shales. Hence the correlation allows an estimate of the contribution of the adsorbed and 'free' gas phases to the overall methane capacity.
Shale and kerogen adsorption and total organic content considerations indicate that approximately half of the CH 4 sorption on dry shales occurs within the organic matter and this indicates the significance of the mineral phase, including the role of clays. Microscopic studies of the nature of pores show the importance of the intra-granular organic-mineral interface and this should be useful in understanding larger scale flow characteristics. Free compressed gas will flow much faster than sorbed gas, where desorption rates decrease with decreasing pore size.
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/.