[en] 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.

[es] Para mitigar las consecuencias del cambio climático se han propuesto diversas medidas, entre las cuales se encuentran el uso de energías renovables como la geotermia y el desarrollo de tecnologías como el almacenamiento de hidrógeno o CO2, que requieren la inyección de fluidos dentro de macizos rocosos fracturados (como los reservorios abandonados), lo cual supone un cambio en el estado tensional del mismo. Dicha alteración puede ocasionar el desencadenamiento de eventos micro-sísmicos mediante la reactivación de fallas geológicas o la creación de nuevas fracturas que en última instancia formarían caminos preferenciales para la filtración de los fluidos. Además, los reservorios abandonados contienen fluidos multifase con más de una especie que podrían reaccionar químicamente con los nuevos fluidos inyectados. Por lo tanto, la predicción precisa de las condiciones de operación requiere un análisis profundo de la problemática, así como disponer un software geomecánico avanzado capaz de evaluar la viabilidad y la seguridad de los posibles lugares de extracción o almacenamiento. En este contexto, el objetivo principal de esta tesis es el desarrollo de un modelo termo-hidro-mecánico (THM) acoplado multifase en un medio poroso fracturado mediante el método de los Elementos Finitos empleando elementos junta de espesor nulo con doble nodo para representar las discontinuidades. Las principales hipótesis del modelo son: la teoría de pequeñas deformaciones, un sistema compuesto por dos fases fluidas (líquida y gaseosa), dos especies (agua y gas seco) y una fase sólida, la fase gaseosa se comporta de acuerdo a la ley de los gases ideales y se considera la evaporación de la fase líquida. Las variables primarias son: desplazamientos, presiones de fase y temperatura. La formulación contempla: la ecuación de equilibrio, la ecuación de balance de masa de la especie agua, la ecuación de balance de masa de la especie gas y el balance de energía. La formulación de las discontinuidades considera las mismas hipótesis que el medio poroso y toma en cuenta un factor de reducción de los parámetros de la curva de retención basado en la apertura mecánica de la fractura. El modelo ha sido implementado y validado con ejemplos de verificación académicos en un código 3D de Elementos Finitos. Sin embargo, durante la fase de validación se han encontrado algunas limitaciones del al modelo, relacionadas con el esquema de integración temporal, en cuanto al tratamiento de los términos de almacenamiento de las ecuaciones de gobierno y los flujos másicos. Estas limitaciones han motivado los siguientes nuevos desarrollos: una nueva formulación de conservación de la masa aplicada a ecuaciones de difusión en el medio poroso y en las discontinuidades, una nueva curva de retención de tipo exponencial que toma como punto de partida la porosimetría en el caso del medio poroso y la distribución de aperturas en las discontinuidades, y, por último, una nueva ley constitutiva cohesivo-friccional para las juntas que permite la modelación del cierre en descarga.