Simulation of the Inelastic Deformation of Porous Reservoirs Under Cyclic Loading Relevant for Underground Hydrogen Storage

Subsurface geological formations can be utilized to safely store large-scale (TWh) renewable energy in the form of green gases such as hydrogen. Successful implementation of this technology involves estimating feasible storage sites, including rigorous mechanical safety analyses. Geological formations are often highly heterogeneous and entail complex nonlinear inelastic rock deformation physics when utilized for cyclic energy storage. In this work, we present a novel scalable computational framework to analyse the impact of nonlinear deformation of porous reservoirs under cyclic loading. The proposed methodology includes three diferent time-dependent nonlinear constitutive models to appropriately describe the behavior of sandstone, shale rock and salt rock. These constitutive models are studied and benchmarked against both numerical and experimental results in the literature. An implicit time-integration scheme is developed to preserve the stability of the simulation. In order to ensure its scalability, the numerical strategy adopts a multiscale fnite element formulation, in which coarse scale systems with locally-computed basis functions are constructed and solved. Further, the efect of heterogeneity on the results and estimation of deformation is analyzed. Lastly, the Bergermeer test case—an active Dutch natural gas storage feld—is studied to investigate the infuence of inelastic deformation on the uplift caused by cyclic injection and production of gas. The present study shows acceptable subsidence predictions in this feld-scale test, once the properties of the fnite element representative elementary volumes are tuned with the experimental data.