Underground hydrogen storage in porous formations is a promising solution for large-scale energy storage. Understanding hydrogen flow and trapping at the pore-scale is crucial for assessing storage capacity and recovery efficiency. While pore-scale flow visualisation experiments provide realistic insights, they are resource intensive and technically challenging. Pore-network models offer a computationally efficient tool for simulating multiphase flow in porous media and can serve as a valuable complement to pore-scale experiments. However, their accuracy remains a key uncertainty and must be evaluated for future application. This study evaluates the performance of a quasi-static pore-network model by comparing its predictions against three-dimensional pore-scale hydrogen flow visualisation experiments in a homogeneous Bentheimer sandstone and a layered Clashach sandstone. The model was calibrated to match experimental end-state saturations, and its performance was evaluated through comparisons of spatial saturation profiles and pore occupancy. The novelty of this study lies in the direct comparison of hydrogen displacement between pore-scale experimental observations and pore-network model simulations, providing an assessment of model performance under varying degrees of rock heterogeneity relevant to underground hydrogen storage. The pore-network model shows good agreement with experimental observations for the homogeneous rock, particularly during drainage, and is subsequently used to analyse additional scenarios, including cyclic hydrogen injection and withdrawal, and wettability variations. These simulations provide insights into capillary pressure behaviour and residual saturation trends. In contrast, for the heterogeneous and layered Clashach sandstone, the model fails to capture the trapping and fluid redistribution observed experimentally during imbibition, revealing limitations in modelling fine-scale heterogeneity.

Understanding pore-scale hydrogen displacement and trapping is crucial for developing subsurface hydrogen storage facilities. While pore-scale flow visualization experiments provide critical insights, they are complex and re source-intensive. Quasi-static pore-network models (PNMs) offer a faster alternative for simulating multiphase flow. This study uses a widely employed PNM to simulate hydrogen flow in sandstones, comparing results with pore-scale flow visualization experiments at reservoir conditions.

Two sandstone samples were used: homogeneous Bentheimer and heterogeneous Clashach. Pore networks were extracted comprising pores and throats, and hydrogen-water flow was simulated, modelling drainage and imbibition processes. Results were analysed for fluid saturations and pore occupancies.

For the homogeneous rock, the PNM matches experimental results for both drainage and imbibition, enabling simulations of different wettability conditions and multiple injection and production cycles. For the heterogeneous rock, the PNM reasonably predicts the hydrogen flow path during drainage but fails to accurately predict imbibition. This discrepancy highlights the limitations of PNMs in predicting pore-scale flow in complex rocks.

In conclusion, while PNMs offer a computationally efficient means to simulate hydrogen flow, they cannot currently replace experimental observations for complex rocks. Further validation against experimental findings is necessary to refine these models and expand their applicability for underground hydrogen storage.

Subsurface porous rocks hold significant hydrogen (H₂) storage potential to support an H₂-based energy future. Understanding H₂ flow and trapping in subsurface rocks is crucial to reliably evaluate their storage efficiency. In this work, we perform cyclic H₂ flow visualization experiments on a layered rock sample with varying pore and throat sizes. During drainage, H₂ follows a path consisting of large pores and throats, through a low permeability rock layer, substantially reducing H₂ storage capacity. Moreover, due to the rock heterogeneity and depending on the experimental flow strategy, imbibition unexpectedly results in higher H₂ saturation compared to drainage. These results emphasize that small-scale rock heterogeneity, which is often unaccounted for in reservoir-scale models, plays a vital role in H₂ displacement and trapping in subsurface porous media, with implications for efficient storage strategies.

Hypothesis

Underground hydrogen (H2) storage is a potentially viable solution for large-scale cyclic H2 storage; however, the behavior of H2 at subsurface pressure and temperature conditions is poorly known. This work investigates if the pore-scale displacement processes in H2-brine systems in a porous sandstone can be sufficiently well defined to enable effective and economic storage operations. In particular, this study investigates trapping, dissolution, and wettability of H2-brine systems at the pore-scale, at conditions that are realistic for subsurface H2 storage.

Experiments

We have performed in situ X-ray imaging during a flow experiment to investigate pore-scale processes during H2 injection and displacement in a brine saturated Bentheimer sandstone sample at temperature and pressure conditions representative of underground reservoirs. Two injection schemes were followed for imbibition: displacement of H2 with H2-equilibrated brine and with non-H2-equilibrated brine. The results from the two cycles were compared with each other.

Findings

The sandstone was found to be wetting to the brine and non-wetting to H2 after both displacement cycles, with average contact angles of 54° and 53°, for H2-equilibrated and non-H2-equilibrated brine, respectively. We also found a higher recovery of H2 (43.1%) when displaced with non-H2-equilibrated brine compared to that of H2-equilibrated brine (31.6%), indicating potential dissolution of H2 in the unequilibrated, imbibing brine at reservoir conditions. Our results suggest that underground H2 storage may indeed be a suitable strategy for energy storage, but considerable further research is needed to fully comprehend the pore-scale interactions at reservoir conditions.