Underground energy storage (UES) in porous and cavity reservoirs can be used to balance the mismatch between the production and demand of renewable energy, as well as for securing gas and oil supply during shortage or high demand periods. Understanding the geomechanical behavior of these reservoirs under different storage conditions, i.e., storage frequency and fluid pressure, is key in defining their capacity and effective lifetime. This thesis work presents a rigorous analysis performed on sandstones to unravel their geomechanical response under cyclic loading. This study includes, importantly, both experimental and numerical investigations under several conditions which are relevant to UES. The rock response was studied considering cyclic stress states above and below the onset of dilatant cracking, under different frequencies and amplitudes. Within the number of cycles studied, measurements of axial strains and acoustic emissions indicated that inelastic strains accumulated cycle after cycle following an exponentially decreasing rate. Five types of deformations were interpreted: elastic, plastic, viscoelastic, cyclic-plastic and brittle creep. Based on these novel experimental results and observations, Nishihara's constitutive model was used for simulating viscoelastic and brittle creep deformations, while dilatant plastic strains were modeled using a Hardening-Softening model. Finally, an extension of the Modified Cam-Clay model was proposed to account for cyclic-plastic compaction. This approach can be extended and improved to study cyclic sandstone deformation's implications on subsidence, fault reactivation and cap rock flexure, among other physical phenomena impacting a reservoir's storage capacity.

A novel 2D finite element method (FEM) on unstructured grid for nonlinear time-dependent deformation of materials is developed. The objective is to model the complex deformation behavior of the rock salt, inside which caverns are mined to store green fuels (such as hydrogen). The analyses and the developments of the present work allow for quantification of the state of the stress of the cavern, and assess the safety and reliability of the storage structure over time. The novelty of the approach is in using minimization of the potential energy principle in conjunction with nonlinear creep deformation physics. While the available FEM-based simulators offer tools only for solving linear and non-linear elastic models, the developed simulator takes into account cyclic loading and material’s damage evolution in time with possibility to predict the material’s failure. Apart from that, the stored product (hydrogen) density is taken into account, which affects cavern’s pressure variation with depth. Impurities,causing heterogeneous rock salt properties, are also considered in the developed model. Multivariate Gaussian distribution is utilized to generate distribution of the heterogeneous mechanical properties. Eulerian strains are introduced in the model to take into account deformation of the mesh. As such,the computational grid changes its geometry according to the deformation. Several numerical test cases are studied. Firstly, a consistency (verification) study is performed, to validate the linear elastic model. Remark that the linear elastic deformation model casts the basis of the non-linear model with creep physics. Then, several studies have been performed to analyze, quantify and approximate the deformation of the salt cavern under the gas pressure change, including the time-dependent creep physics.