Hydrogen–Rock Interactions in Carbonate and Siliceous Reservoirs: A Petrophysical Perspective

1. Introduction

In light of the rapid expansion of hydrogen production and its adoption as a clean energy carrier, there is a growing need for subsurface geological storage solutions capable of balancing fluctuations in supply driven by variable renewable electricity input and seasonal variations in demand, thereby enhancing long-term energy security [1]. While salt caverns are considered a promising option due to their operational flexibility [2,3], porous reservoirs such as depleted gas fields [4] and deep saline aquifers [5] despite their broader availability still require comprehensive investigation to assess their geochemical behavior and effectiveness as safe and efficient media for hydrogen storage[6,7].

The study of geochemical processes within underground hydrogen storage systems is of critical importance, particularly those that may contribute to hydrogen loss during storage operations [8]. Although the past decade has seen numerous investigations ranging from geochemical modeling and to experimental studies [9,10,11,12] across various rock types the current understanding of the physicochemical interactions between hydrogen and geological formations remains limited and requires further systematic development [13,14].

Hydrogen loss during high-pressure underground hydrogen storage (UHS) is primarily attributed to geochemical and microbial interactions that typically occur under reservoir conditions characterized by temperatures below 127 °C and pressures below 350 bar. Understanding these interactions is essential for evaluating hydrogen mobility and reactivity within geological media, which is key to enhancing storage monitoring and ensuring long-term operational stability [15].

Hydrogen-rock interactions are influenced by various factors, including mineral composition, storage pressure and temperature, pore structure, and the rock’s specific properties such as heterogeneity, mechanical stability, organic content, and wettability behavior [16]. These physical, chemical, and biological interactions can induce changes in permeability, porosity, and wettability [15,17,18], as well as trigger microfracture activation directly impacting the long-term safety and efficiency of geological hydrogen storage systems [19].

In this context, Al-Yaseri et al. [17] investigated carbonate cores exposed to hydrogen for 75 days (48 bar, 75°C), revealing that calcite expansion reduced effective porosity by 47%, significantly impairing limestone storage capacity. In contrast, dolomite showed a slight (~6%) improvement due to a balance between dissolution and grain expansion. In another study by Al-Yaseri et al. [9] exposed calcite-rich shale samples from the EFS and WCS formations to hydrogen at 9.65 MPa and 348 K for 108 days. Calcite precipitation occurred due to organic matter interactions, with no sulfate or pyrite reduction. Minimal pore changes were observed, indicating preserved rock integrity and suitability for geological hydrogen storage. Bo et al.[20] suggest that calcite-free reservoirs and caprocks are geochemically more suitable for underground hydrogen storage, as carbonate presence particularly calcite can induce hydrogen dissociation during dissolution, leading to losses of up to 9.5%. The study by Cheng et al.[21] demonstrated that hydrogen interaction with pyrite-free Buntsandstein sandstone is limited, even after exposure to 150 bar and 100°C for four weeks. No significant mineralogical transformations were observed, and no noticeable changes were detected in key petrophysical properties such as porosity and permeability, confirming the rock’s suitability for hydrogen storage. Flesch et al. [22]investigated how sandstone formations respond petrographically and petrophysically to hydrogen exposure. The study revealed a moderate decline in median porosity from 13.3% to 11.9% accompanied by a notable increase in permeability, rising from 48.78 mD to 61.53 mD. Interestingly, in some cases, an inverse behavior was observed, with porosity increasing while permeability declined. The observed porosity increase was primarily attributed to the dissolution of pore-filling phases such as anhydrite and carbonate minerals, whereas permeability remained relatively unaffected. In a recent study, Al-Yaseri et al. [23] examined the geochemical interaction of hydrogen with sandstone and limestone formations under 75°C and 96 bar for six months. Only minor morphological changes were observed, with slight reductions in porosity from 14.4% to 13.33% for sandstone and from 20.2% to 18.44% for limestone. These results indicate limited H₂ reactivity. Overall, the minimal hydrogen loss confirms the potential of sandstone and limestone as suitable candidates for underground hydrogen storage.

Building upon the existing knowledge gaps concerning the mineralogical and petrophysical responses of geological formations to hydrogen exposure, this study adopted a comprehensive experimental methodology to evaluate three representative lithologies: dolomite, calcite-rich limestone, and quartz-dominated siliceous rock. The investigations focused on quantifying hydrogen-induced changes in crystallographic structure and in characteristics related to: porosity, permeability, and elemental distribution. By integrating XRD, BET/BJH, SEM-EDS and poro-perm analyses, the study allowed a rigorous comparison of reactivity and structural integrity between the analyzed rock types, under reservoir-relevant conditions. The results discussed in the following section provide valuable insights into the physicochemical mechanisms governing hydrogen-rock interactions, with implications for assessing.

2. Materials and Methods

The research program included laboratory experiments that aimed to evaluate the effect of hydrogen gas (Linde, purity 5.0) on three distinct rock types: dolomite (CaMg(CO₃)₂), vaterite-rich limestone (CaCO₃), and a siliceous sample composed predominantly of quartz. These samples were labeled A, B, and C, respectively, as shown in Figure 1. The experiments were designed to simulate subsurface underground hydrogen storage conditions.

To identify geochemical alterations and variations in petrophysical properties induced by hydrogen exposure, a combination of analytical techniques was employed. The petrophysical properties of the rock cores were characterized both before and after hydrogen treatment to assess the extent of any changes.

Mineralogical analysis was performed using X-ray diffraction (XRD) with a D8 Advance diffractometer (Bruker-AXS, Karlsruhe, Germany) operating with Cu-Kα radiation (λ = 1.54 Å), equipped with a nickel filter, θ-θ geometry, and a Bragg–Brentano configuration. XRD measurements were carried out under the following conditions: 40 kV, 40 mA, a step size of 0.1°, a scan speed of 0.1° per 5 seconds, and a 2θ scan range from 5° to 60°.

Permeability and porosity measurements of the rock cores were conducted using a GASPERM permeameter (Vinci Technologies, France) and a helium porosimeter (He-Porosimeter, Vinci Technologies, France), respectively.

To analyze textural changes such as BET surface area, pore size distribution, and mean pore diameter, nitrogen adsorption–desorption isotherms were recorded using a Quantachrome NOVA 2200e gas adsorption analyzer (Boynton Beach, FL, USA). Isotherms were collected at 77.35 K over a relative pressure range (p/p₀) of 0.005 to 1.0, and the data were processed using NovaWin software, version 11.03.

Additionally, morphological changes and elemental composition of the rock samples before and after hydrogen treatment were examined using scanning electron microscopy (SEM). A Hitachi S-3400N SEM and an FEI SEM system (Hillsboro, OR, USA), both equipped with energy-dispersive X-ray spectroscopy (EDS), were utilized. SEM images were acquired at an accelerating voltage of 2 kV, while EDS analyses were conducted at 30 kV.

Initially, cylindrical rock samples (core plugs) were placed in a stainless-steel test cell. The cell had a diameter of 3.5 cm, a wall thickness of 1 cm, and a length of 50 cm, and was designed to withstand pressures up to 250 bar. It was equipped with a gas inlet fitted with a pressure gauge and a sealed outlet connected to a flowmeter for monitoring hydrogen flow.

Hydrogen gas was introduced into the cell until a pressure of 100 bar was reached, using a 50 L hydrogen cylinder pressurized to 200 bar and regulated with a pressure controller. The rock samples were exposed to the hydrogen atmosphere for 100 days under constant pressure (100 bar). The internal temperature of the cell was monitored using a thermocouple and maintained at 70°C throughout the experiment. After exposure, the samples were retrieved and stored under a hydrogen atmosphere for subsequent characterization analyses.

3. Results and Discussions

3.1. X-Ray Diffraction (XRD) Analysis

The X-ray diffraction (XRD) analysis conducted on the untreated dolomite sample (Core A) revealed distinct and intense reflections at 2θ value equal with 24.10°, 30.88°, 37.25°, 41.05°, 44.85° and 50.90°, corresponding to the characteristic crystallographic planes of dolomite [24,25,26], as presented in Figure 2. Additionally, minor peaks related to chlorite and calcite were identified at 6.2° and 29.5°, attributed to the (001) and (104) planes, respectively [26,27]. Quartz-related peaks were also observed at 20.8° and 26.6°, assigned to the (100) and (101) planes, respectively [28], alongside a weak reflection for illite at 8.9° (001 plane) [29].

Following hydrogen exposure at 100 bar and 70 °C, the XRD patterns of the treated dolomite sample (Core AH) (Figure 2) revealed significant mineralogical and structural transformations. A marked reduction in the intensity of the characteristic reflections at 24.1° and 41.05°, corresponding to the (012) and (113) planes, was observed. Additionally, reflections associated with calcite, chlorite, and illite diminished or disappeared entirely. These changes are attributed to the ability of hydrogen to induce partial dissolution and destabilization of the mineral’s crystalline structure, possibly due to phase transformation or local rearrangement of the crystal lattice [30] . The process is illustrated by the dissolution mechanism of dolomite, which leads to the release of constituent ions and rearrangement within the mineral matrix, the following reaction describes the hydrogen-induced breakdown of dolomite in the presence of moisture [17,30,31,32]:

$$CaMg(C{O}_3{)}_{2\left(s\right)}+4{H_2}_{\left(g\right)}\rightleftharpoons CaC{O}_{3\left(s\right)}+Mg(OH{{)}_2}_{\left(s\right)}+C{H_4}_{\left(g\right)}+H_2{O}_{\left(l\right)}$$

(1)

$$CaMg(C{O}_3{)}_{2\left(s\right)}+10{H_2}_{\left(g\right)}\rightleftharpoons C{a}^{+2}+M{g}^{+2}+2C{H_4}_{\left(g\right)}+6H_2{O}_{\left(l\right)}+4e^{-}$$

(2)

These pathways suggest that dolomite undergoes hydrogen-driven destabilization, releasing Ca²⁺ and Mg²⁺ ions, hydroxides, and methane, thus altering its crystalline framework.

In the case of the limestone sample (Core B), the XRD pattern depicted in Figure 3 confirmed the dominance of carbonates (calcite and dolomite) as the primary crystalline phases. Calcite (CaCO3) shows sharp peaks at 29.4°, 36.10°, 39.38°, 43.2°, 47.4°, 48.50° matching the (104), (110), (113), (202), (018), and (116) planes. On the other hand, dolomite show significant peaks at 30.88°, 56.2°, and 57.3°, corresponding to (104), (018), and (116) planes [24,26,27]. Furthermore, vaterite (CaCO3 biogenic polymorph) was detected through additional peaks at 24.8°, 27.1°, 32.7°, 43.8°, and 55.7°, corresponding to the (110), (112), (114), (300), and (224) planes [27].

Following hydrogen exposure, the treated sample Core BH (Figure 3) exhibited the complete disappearance of several calcite-associated peaks, including those at 31.3°, 43.2°, and 57.3° ((006), (202), and (122) planes), along with the 20.9° peak corresponding to the (004) plane of vaterite. These changes imply strong reactivity between hydrogen and carbonate phases, particularly calcite, resulting in progressive dissolution [17]. The reduction of carbonate by H₂ likely increases the acidity of the pore fluid[9,15], thereby accelerating the dissolution process and possibly leading to methane formation, as outlined in reaction (3) [32].

$$CaC{O}_3+5H_2\leftrightarrow C{a}^{+2}+C{H}_4+3H_{2}O+2e^{-}$$

(3)

The siliceous rock sample (Core C), prior to hydrogen exposure, exhibited a mineralogical composition dominated by quartz, accompanied by minor amounts of chlorite and albite, as shown in Figure 4. The diffractogram revealed prominent quartz peaks at 20.8°, 26.7°, 36.45°, 39.5°, 42.5°, 45.9°, 50.2°, 54.9 and 60.2°, corresponding to the (100), (101), (110), (102), (200), (201), (112), (202), and (211) planes [28,33]. In addition, chlorite reflections were noted at 6.2° and 12.5° ((001) and (002) planes), with a distinct illite signal at 8.9° (100) [29], and peaks at 27.92° (albite (112)) [33] and 29.4° (calcite (104)) [27], indicating minor accessory phases.

Following hydrogen exposure, the XRD analysis of the siliceous rock sample Core CH (Figure 4) showed that the diffraction peaks associated with quartz remained unchanged after hydrogen exposure compared to other minerals such as feldspars [30]. This observation confirms the exceptional chemical stability of quartz in reactive hydrogen-rich environments, thereby reinforcing its suitability as a long-term host matrix for underground hydrogen storage.

These findings establish the baseline mineralogical composition of the three core samples, serving as a reference for comparative analysis after hydrogen exposure.

Table 1, Table 2 and Table 3 provide quantitative insights into the compositional changes in core samples A, B, and C before and after hydrogen exposure, highlighting the extent of mineral dissolution and secondary transformations as reflected in the variation of mineral content (%). These variations, more than peak position shifts, directly illustrate the chemical reactivity and structural resilience of different lithologies under high-pressure hydrogen conditions.

In Table 1, the dolomitic core (A) exhibits a clear reduction in dolomite content after hydrogen treatment (from 12.20% to 10.53% for the dominant (104) plane, and complete disappearance of other dolomite components), confirming selective dissolution of CaMg(CO₃)₂. This process is facilitated by hydrogen in the presence of moisture, which triggers partial decarbonation and release of Ca²⁺ and Mg²⁺ ions. Concurrently, the complete loss of calcite (from 1.32% to 0%) and albite (2.03% to 0%) suggests a destabilization of accessory phases, while the increase in the (116) dolomite plane (from 2.68% to 2.84%) may indicate local recrystallization or lattice reorganization, as mentioned above. These compositional shifts align with the dissolution–reprecipitation mechanism and support the formation of secondary Fe-rich carbonates (ankerite), as also confirmed by EDS.

In Table 2, the limestone core (B) demonstrates a sharp decrease in vaterite content (e.g., from 6.05% to 4.82% at 27°, and complete loss at 20.9°), as well as the disappearance of several calcite planes (e.g., from 2.35% to 0% at 31.2° and from 1.11% to 0% at 43.1°). These losses indicate intensive carbonate dissolution, likely promoted by pore fluid acidification due to hydrogen-carbonate interactions. The slight increase in some vaterite and calcite reflections (e.g., 43.8° vaterite from 2.98% to 3.50%, 50° calcite from 2.68% to 2.94%) suggests secondary carbonate precipitation possibly clogging pore throats and contributing to the observed drop in porosity and permeability. This mineralogical behavior reinforces the idea that calcitic formations are structurally vulnerable under reductive hydrogen environments.

Table 3 shows minimal compositional changes in the siliceous core (C), particularly for quartz, whose dominant reflections remain almost unchanged (e.g., 14.64% to 13.87% for 26.7°, and 3.57% to 3.78% for 20.9°). Minor losses in illite (2.62% to 0%) and chlorite variations (6.87% to 5.94% for 6.3°, but 2.19% to 3.22% for 12.6°) are within the experimental margin and may reflect superficial rearrangements rather than true dissolution. This overall mineralogical stability strongly supports the chemical inertness of quartz-rich formations under hydrogen storage conditions.

3.2. Porosity, Permeability and Textural Characteristics

The nitrogen adsorption–desorption isotherms obtained at 77 K for core samples A and AH (Figure 5) exhibit the characteristics of mesoporous materials, showing a type IV profile in accordance with the IUPAC classification. The presence of hysteresis loops confirms capillary condensation in pores with intricate geometries. Specifically, core A displays a hysteresis loop of the H2 type, associated with “Ink-Bottle” shaped pores, whereas core AH exhibits an H3-type loop, indicative of Groove pores [34,35,36].

Hydrogen exposure at 100 bar for 100 days led to pronounced morphological alterations in core AH, resulting in a 61.57% increase in the adsorbed nitrogen volume at high relative pressures from 14.376 cm³/g (A) to 23.35 cm³/g (AH) as shown in Figure 5. Although the total porosity variation was negligible (23.95% for A vs. 23.87% for AH), the pore network in AH appears significantly more extended and interconnected.

BET and BJH analyses corroborate this transformation: the specific surface area increased from 3.270 m²/g (A) to 4.023 m²/g (AH), and the total pore volume rose from 0.0214 cm³/g to 0.245 cm³/g. Additionally, the mean pore diameter (Dv(d)) increased from 3.3 nm to 5.09 nm (Figure 6), reflecting a more open and fluid-accessible mesoporous architecture.

The pore size distribution further reveals a higher density of small pores in AH, contributing to enhanced adsorption capacity and internal network expansion. These structural modifications positively influence permeability, which increased from 574.84 mD (A) to 592.44 mD (AH), suggesting improved pore connectivity and possibly reduced tortuosity.

These textural and transport enhancements are primarily attributed to selective dolomite (CaMg(CO₃)₂ dissolution under reductive hydrogen exposure [17,30,31,32], as described by reactions (1) and (2), along with the concurrent precipitation of ankerite (Fe-carbonate). The elevated presence of ankerite likely originates from the reduction of siderite, while the cationic substitution of Mg²⁺ by Fe²⁺ in the dolomite matrix led to the formation of ferroan carbonates (reaction 4) [31].

$$CaMg(C{O}_3{)}_2+F{e}^{+2}+H_2\rightleftharpoons Ca(Fe,Mg\left)\right(C{O}_3{)}_2+H_{2}O$$

(4)

Accordingly, dolomitic formations emerge as promising candidates for subsurface hydrogen storage. These findings extend current literature [17,37], highlighting that not only the geochemical reactivity but also the physical structure of the host rock is critical in governing mineralogical transformation kinetics in reductive environments.

The nitrogen adsorption–desorption isotherms at 77 K for all studied core samples (B/BH, and C/CH) reveal mesoporous characteristics with type IV profiles and H3-type hysteresis loops and H4-type hysteresis respectively according to IUPAC classification, denoting capillary condensation within groove pores and silt pores, indicates the presence of a developed but poorly connected pore network typical of carbonate rocks with irregular medium porosity [34,35,36]. However, the impact of hydrogen exposure (100 bar for 100 days) differed significantly depending on mineralogical composition

Upon hydrogen exposure (100 bar for 100 days), core BH exhibited a considerable reduction in adsorption capacity throughout the relative pressure range. At the maximum relative pressure, the nitrogen adsorption volume dropped by approximately 19%, from 26.11 cm³/g (B) to 21.17 cm³/g (BH) As observed in Figure 7 suggesting significant structural alterations in the pore network caused by the hydrogen treatment.

Textural parameter analysis confirmed this trend: the BET specific surface area decreased from 6.93 m²/g (B) to 5.80 m²/g (BH), while the total pore volume (BJH) declined from 0.0175 cm³/g to 0.0143 cm³/g. Additionally, the mean pore diameter (Dv(d)) experienced a slight reduction from 3.023 nm to 2.853 nm. These variations align with the decrease in total porosity (from 19.89% to 16.33%) and permeability (from 132.16 mD to 107.43 mD), indicating diminished pore connectivity and a potential partial collapse of the porous framework.

The pore size distribution (Figure 8) further supports this observation, revealing a reduced density of active pores especially within the 3–6 nm range in the hydrogen-treated sample. Such transformations are likely attributable to mineral reprecipitation phenomena or localized pore blockage, possibly induced by the formation of secondary carbonate phases (e.g., redeposited calcite), triggered by solid gas reactions under high hydrogen pressure [32,38].

These findings underscore the negative impact of prolonged hydrogen exposure on the structural integrity of calcitic rocks. Unlike dolomite, which exhibited pore network enhancement, calcite underwent compaction and occlusion, leading to reduced storage potential. Consequently, calcitic formations appear less favorable for long-term hydrogen containment due to their susceptibility to morphological degradation under reductive conditions.

In contrast, core samples C and CH, composed primarily of quartz, exhibited remarkable structural resilience following exposure to hydrogen at 100 bar for 100 days (Figure 9). However, the structural response of quartz was notably more stable.

The adsorbed nitrogen volume showed only a marginal decrease of ~1.5%, from 27.81 cm³/g (C) to 27.40 cm³/g (CH), significantly lower than the changes observed in dolomite (+61.57%) and calcite (−19%). This indicates that hydrogen has a minimal effect on the microstructural integrity of quartz.

This conclusion is further supported by the textural parameters: the BET surface area slightly decreased from 19.048 m²/g to 18.87 m²/g, the total pore volume (BJH) declined marginally from 0.0291 cm³/g to 0.0284 cm³/g, and the mean pore diameter (Dv(d)) increased moderately from 3.107 nm to 3.305 nm. Total porosity and permeability also exhibited only slight reductions, from 24.5% to 24.13%, and from 581.96 mD to 567.43 mD, respectively.

The pore size distributions (Figure 10) for samples C and CH are nearly superimposed, indicating negligible structural changes. The stable pore population within the 3–10 nm range confirms the absence of collapse, obstruction, or secondary phase precipitation, in contrast to the behavior observed in dolomite and calcite.

Quartz, as a silica-dominated mineral phase, demonstrates outstanding stability under high-pressure hydrogen conditions, maintaining both its structural and petrophysical integrity. Unlike dolomite, which undergoes pronounced pore enlargement, or calcite, which shows susceptibility to partial network collapse, quartz exhibits negligible microstructural degradation, thereby reinforcing its candidacy as a robust matrix component in subsurface hydrogen storage environments [15,20,37,39].

This stability is quantitatively reflected in the minimal variations recorded across key textural descriptors: the BET surface area decreased only slightly (from 19.048 to 18.87 m²/g), total pore volume showed a marginal reduction (0.0291 to 0.0284 cm³/g), and the average pore diameter increased modestly (3.107 to 3.305 nm). Additionally, total porosity declined from 24.5% to 24.13%, and permeability from 581.96 to 567.43 mD. These minor shifts, detailed in Table 4, underscore the resilience of the quartz framework and further validate its functional neutrality in the face of hydrogen-induced chemical and physical stress.

3.3. SEM-EDS Analysis

The SEM micrographs acquired at magnifications of 100× and 200× reveal substantial morphological differences between untreated and hydrogen-exposed carbonate core samples. As illustrated in Figure 11, core A (dolomite, untreated) displays a relatively compact granular structure, composed of well-defined dolomite crystals with sharp boundaries and tight grain contacts features typical of a preserved stoichiometric carbonate phase with a uniform crystalline framework [40].

Conversely, the hydrogen-exposed counterpart (core AH) exhibits a more developed porous texture, with partially eroded grains, irregular margins, and expanded intergranular voids. Fine particles are observed to accumulate between larger crystals, indicating localized dolomite dissolution and possible reprecipitation of secondary phases [17,30,31,32]. These morphological observations corroborate the BET/BJH adsorption data and confirm partial reorganization of the pore network following hydrogen exposure.

A similar trend is observed for calcitic samples. As shown in Figure 12, SEM images of core B reveal a relatively dense microstructure formed by large, well-faceted calcite crystals with regular faces and compact intergranular contacts. The pore network is poorly developed, with minimal connectivity. After hydrogen treatment (core BH), the microstructure becomes noticeably irregular, characterized by fragmented crystal surfaces, fine fissures, and increased grain boundary roughness. While intercrystalline voids are more apparent, they do not translate into enhanced effective porosity, likely due to partial pore clogging by redeposited carbonate phases (Figure 12).

Quartz-based samples exhibit a markedly different behavior. As seen in Figure 12, core C shows a compact granular architecture with large, smooth-faced quartz crystals and limited intergranular space. Following hydrogen exposure, core CH maintains this structure with only minor surface roughening and the appearance of microcracks. The pore network remains largely unaffected, with no signs of structural collapse or significant mineralogical alteration, confirming the high chemical and structural stability of quartz under reductive conditions [15,20,37,39].

Figure 13. Scanning Electron Microscopy (SEM) Images of Core Samples C and CH.

Figure 13. Scanning Electron Microscopy (SEM) Images of Core Samples C and CH.

EDS analysis provides further insights into the chemical evolution associated with hydrogen treatment. For dolomitic cores, EDS spectra and elemental maps Figure 14 show that core A is dominated by calcium and magnesium with uniform distribution, consistent with stoichiometric CaMg(CO₃)₂. Iron is present only in trace amounts. After hydrogen exposure, core AH exhibits a notable increase in iron concentration and a relative decrease in magnesium. The heterogeneous distribution of Fe suggests cation exchange, where Fe²⁺ likely sourced from siderite reduction partially replaces Mg²⁺ in the dolomite lattice, leading to ankerite-type phase formation [31].

Other elements (Si, Al, K) show no significant compositional variation between A and AH, indicating that these elements are not actively involved in hydrogen-induced reactions (Figure 14). This SEM–EDS evidence confirms that dolomite is susceptible to both structural and chemical alteration under hydrogen exposure, primarily via selective dissolution, Fe for Mg substitution, and secondary carbonate reprecipitation. These changes reshape the pore framework and influence the rock’s petrophysical behavior, which is critical when assessing its suitability for hydrogen storage applications.

EDS analysis of the calcitic cores (Figure 14) reinforces these trends. Core B exhibits a composition dominated by Ca and O with a uniform distribution of minor elements. In core BH, a slight increase in Fe is detected, with a more heterogeneous spatial pattern, suggesting secondary Fe-rich phase formation driven by hydrogen-induced redox reactions. Redistribution of Mg and minor elements is also observed, indicating localized elemental migration within the calcite matrix.

Figure 15. EDS Spectrum and Elemental Distribution for Core Samples B and BH.

Figure 15. EDS Spectrum and Elemental Distribution for Core Samples B and BH.

In contrast, EDS results for the quartz samples (C and CH, Figure 16) reveal consistent compositional stability. Both spectra show Si and O as dominant elements with uniform distributions of Al, K, Mg, and Fe. No notable differences are observed in elemental composition or spatial distribution after hydrogen treatment, further affirming the chemical inertness of quartz. Meanwhile, major elemental redistributions were observed in AH and BH samples, linked to ankerite formation (in dolomite) and secondary calcite precipitation (in calcite).

These microstructural and compositional findings are consistent with textural parameter evolution. In calcite, hydrogen exposure reduced BET surface area from 6.93 to 5.80 m²/g, total pore volume from 0.0175 to 0.0143 cm³/g, mean pore diameter from 3.023 to 2.853 nm, porosity from 19.89% to 16.33%, and permeability from 132.16 to 107.43 mD highlighting partial pore collapse and occlusion due to carbonate reprecipitation.

In contrast, the dolomitic system demonstrated a distinct structural response: after hydrogen exposure, total pore volume increased from 0.0214 to 0.245 cm³/g, average pore diameter expanded from 3.3 to 5.09 nm, and permeability improved alongside surface area growth. SEM images of core AH confirmed severe grain disaggregation and development of an interconnected pore system, enriched with Fe-bearing phases attributed to ankerite formation.

These contrasting behaviors clearly illustrate that calcite and dolomite respond differently to hydrogen treatment. While dolomite undergoes beneficial structural modifications enhancing porosity and permeability, calcite exhibits degradation of the pore network through mineral reprecipitation and occlusion. This distinction is crucial for evaluating the petrophysical suitability of carbonate rocks in geological hydrogen storage, with dolomite emerging as a more promising candidate than calcite.

4. Conclusions

This section is not mandatory but can be added to the manuscript if the discussion is unusually long or complex.

This experimental study – which was based on the simulation of geological hydrogen storage conditions (100 bar, 70°C, for 100 days) - revealed significant differences in the mineralogical and petrophysical responses of three representative rock types dolomite, calcite-rich limestone, and quartz-rich siliceous rock. The results demonstrated that dolomite exhibited moderate but constructive geochemical reactivity, leading to a remarkable enhancement in its pore architecture. Permeability increased from 574.84 to 592.44 mD (+3.1%), BET surface area rose from 3.27 to 4.02 m²/g (+22.9%), and the average pore diameter expanded from 3.3 to 5.09 nm (+54.2%). Most notably, total pore volume increased from 0.0214 to 0.245 cm³/g, marking a dramatic +1045% improvement. Importantly, these changes occurred without significant variation in total porosity (23.95% to 23.87%), indicating improved internal pore connectivity rather than simple void expansion. These enhancements were driven by selective dolomite dissolution under moist hydrogen conditions and the formation of ankerite-type ferroan carbonates via Fe²⁺ ↔ Mg²⁺ substitution mechanisms.

In contrast with the dolomite, the calcite-rich limestone samples exhibited a pronounced deterioration in storage-relevant properties. Porosity declined from 19.89% to 16.33% (−17.9%), permeability dropped from 132.16 to 107.43 mD (−18.7%), surface area decreased from 6.93 to 5.80 m²/g (−16.3%), and total pore volume was reduced from 0.0175 to 0.0143 cm³/g (−18.3%). These reductions are attributed to intense carbonate dissolution followed by secondary carbonate precipitation, leading to partial pore clogging and loss of effective pore space connectivity.

Research conducted on quartz-dominant siliceous rock has shown that it maintained near-complete structural and chemical stability. The porosity of these samples remained nearly unchanged (24.50% to 24.13%), and permeability dropped only slightly (581.96 to 567.43 mD, −2.5%). Changes in BET surface area and pore volume were minimal (within ±1%), confirming the mineral’s inert nature and structural resilience under reductive hydrogen conditions.

In summary, dolomite emerges as the most promising and practically relevant lithology among those tested. It balances controlled chemical reactivity with favorable enhancements in pore structure and permeability without compromising integrity making it a strong candidate for underground hydrogen storage applications. In contrast, calcite-rich formations present structural instability risks, while quartz may serve as a reliable inert matrix within mixed lithological systems. These findings contribute valuable scientific insight for the strategic selection of geological formations suited for long-term green hydrogen storage.