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 with the Mohr-Coulomb perfectly plastic yield criterion is implemented. The solution is then tested on a geological scenario with simplified geometry but realistic rock parameters. Results are compared with the solution of a first-order FEM obtained from a commercial solver, showing a good agreement in terms of displacement and stress maps. In the current implementation, a stabilization term of the VEM projectors is tailored to the elastoplastic law, paving the way for the generalization to polyhedral grids and the introduction of domain discontinuities such as faults
A VEM application to geomechanical simulations for underground fluid storage / Raeli, Alice; Serazio, Cristina; Marcon, Francesca; Borio, Andrea; Benlalam, Nacer; Berrone, Stefano; Verga, Francesca. - In: GEAM. GEOINGEGNERIA AMBIENTALE E MINERARIA. - ISSN 1121-9041. - (2023), pp. 35-46. [10.19199/2023.170.1121-9041.035]
A VEM application to geomechanical simulations for underground fluid storage
Alice Raeli;Cristina Serazio;Francesca Marcon;Andrea Borio;Nacer Benlalam;Stefano Berrone;Francesca Verga
2023
Abstract
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 with the Mohr-Coulomb perfectly plastic yield criterion is implemented. The solution is then tested on a geological scenario with simplified geometry but realistic rock parameters. Results are compared with the solution of a first-order FEM obtained from a commercial solver, showing a good agreement in terms of displacement and stress maps. In the current implementation, a stabilization term of the VEM projectors is tailored to the elastoplastic law, paving the way for the generalization to polyhedral grids and the introduction of domain discontinuities such as faultsPubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
https://hdl.handle.net/11583/2991622
Attenzione
Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo