THM numerical modeling of multiphase flow in fractured rock masses using zero-thickness interface elements

(English) To address the impacts of climate change, several measures have been suggested among which utilizing renewable energy sources such as geothermal energy and implementing technologies such as hydrogen storage or CCUS. Some of these solutions imply the injection of fluids into fractured rock...

Descripción completa

Detalles Bibliográficos
Autor: Barandiarán Villegas, Lucía Belén
Tipo de recurso: tesis doctoral
Estado:Versión publicada
Fecha de publicación:2024
País:España
Institución:CBUC, CESCA
Repositorio:TDR. Tesis Doctorales en Red
OAI Identifier:oai:www.tdx.cat:10803/692713
Acceso en línea:http://hdl.handle.net/10803/692713
https://dx.doi.org/10.5821/dissertation-2117-420576
Access Level:acceso embargado
Palabra clave:Àrees temàtiques de la UPC::Enginyeria civil
624
Descripción
Sumario:(English) To address the impacts of climate change, several measures have been suggested among which utilizing renewable energy sources such as geothermal energy and implementing technologies such as hydrogen storage or CCUS. Some of these solutions imply the injection of fluids into fractured rock masses (e.g. depleted O&G reservoirs) with the subsequent change in local stress conditions. This alteration may lead to the initiation of micro-seismic events by reactivating existing faults or the generation of new fractures, with the risk of ultimately establishing pathways for fluid leakage. Additionally, depleted reservoirs contain multiphase and multicomponent fluids that could react chemically with the new injection fluids. Therefore, the accurate prediction of operation conditions demands advanced geomechanical software capable of evaluating the viability and safety of potential extraction or storage locations. In this context, the main objective of this thesis is to develop a thermo-hydro-mechanical multiphase model in fractured porous medium with the Finite Element Method using zero-thickness interface elements to model discrete discontinuities. The main assumptions of the model are: the small strain theory, the system composition with two fluid phases (liquid and gas), two species (water and dry gas) and one solid phase. The gas phase follows the ideal gas law and the liquid phase evaporation is considered. The main variables are: displacements, phase pressures and temperature. The formulation takes into account: the equilibrium equation, the mass balance equation of the water species, the mass balance equation of the gas species and the balance of energy. The formulation of the discontinuities considers similar assumptions as the porous medium and takes into account a reduction factor of the water retention curve parameters based on the mechanical opening of the fracture. The model has been implemented and validated with academic verification examples in a FE code. However, some limitations were observed related to the time integration scheme, the treatment of storage terms of the governing equations and the mass flow. This issues have motivated the following new developments: new mass-conservative approach for diffusion equations using FE formulations, a new physically-based exponential model which considers the pore size or aperture distribution in case of the porous medium or discontinuities, respectively, and a new damage-frictional law which properly models discontinuity closure.