In this paper, we discuss how to describe the geomechanical behavior of geological formations used for underground fluid storage through the application of the Virtual Element Method (VEM) on conforming polyhedral meshes for the solution of stress-strain equilibrium equations. Under the assumption of small deformations, the solution algorithm for an Isotropic Linear Elastic (ILE) constitutive law coupled...