Phd Position F/m Numerical Simulation of Coupled Seismo-hydro-mechanical Processes in Seismicity Induced By Subsurface Fluid Injection

Le descriptif de l’offre ci-dessous est en Anglais_ **Type de contrat**: CDD **Niveau de diplôme exigé**: Bac + 5 ou équivalent **Fonction**: Doctorant **A propos du centre ou de la direction fonctionnelle**: The Inria centre at Université Côte d'Azur includes 42 research teams and 9 support services. The centre's staff (about 500 people) is made up of scientists of different nationalities, engineers, technicians and administrative staff. The teams are mainly located on the university campuses of Sophia Antipolis and Nice as well as Montpellier, in close collaboration with research and higher education laboratories and establishments (Université Côte d'Azur, CNRS, INRAE, INSERM...), but also with the regiona economic players. With a presence in the fields of computational neuroscience and biology, data science and modeling, software engineering and certification, as well as collaborative robotics, the Inria Centre at Université Côte d'Azur is a major player in terms of scientific excellence through its results and collaborations at both European and international levels. **Contexte et atouts du poste**: The PhD is co-supervised by Roland Masson (LJAD, Inria) and Jean-Paul Ampuero (Géoazur, IRD, Inria), and will be carried out in collaboration with other members of the Inria project team Galets as well as with IFPEN. The PhD will take place at the J.A. Dieudonné Mathematics Laboratory (LJAD) at the University Côte d'Azur as part of the joint Galets project team, which involves Inria, LJAD, and Géoazur. **Mission confiée**: Subsurface resources such as deep geothermal energy, underground hydrogen storage, and geological CO2 sequestration have become crucial pillars of the energy transition and sustainable development. However, injecting or extracting fluids into or from the subsurface modify the pressure and stress state of the surrounding rocks, often extending far beyond the injection point. These changes can potentially trigger fault reactivation. Such phenomena pose seismic risks that must be better understood and mitigated to ensure operational safety, reduce risks to nearby populations, and enhance the social acceptance of these projects If it remains sufficiently low, this induced seismicity (then referred to as microseismicity) raises no issue and can even contribute to site monitoring. However, the ability to predict and mitigate the risks of induced seismicity is key for sustainable exploitation of the subsurface. In recent years, anthropogenic activities in storage reservoirs have triggered earthquakes with magnitudes up to ~6, high enough to cause damage, in regions where natural seismic activity was otherwise low. A recent example is a magnitude 3.9 earthquake in Vendenheim, Alsace, which led to the closure of this deep geothermal site in 2021. In this context, numerical simulation stands out as a key tool for better understanding, predicting, and managing these phenomena. It allows to account for coupled multiphysical processes underlying induced seismicity, evaluating the potential impacts of human activities on the subsurface, and designing strategies to minimize associated risks. At sufficiently large spatial scales, the models are based on a represention of the faults as interfaces equipped with mechanical properties (normal stiffness, friction laws, etc.) and hydraulic properties (aperture and hydraulic conductivity). The physical model couples Darcian fluid flow in the porous matrix and the fault network, the poro-elastic deformation of the rock in the matrix domain as well as the highly nonlinear frictional slip laws along the fault network. An essential feature is the dynamic nature of the friction law, commonly described by a Rate and State law [Pipping2016,Ampuero2011,Uphoff22] and for which the friction coefficient depends on the slip velocity (Rate) as well as on a State parameter accounting for the average contact age or the maturity of the fault asperities and allowing to model seismic cycles. - The goal of the PhD thesis is to develop numerical methods to efficiently simulate these strongly coupled Seismo-Hydro-Mechanical processes, taking into account 3D fault networks and dynamic rate and state friction laws. - We will first focus on quasi-dynamic models [Ampuero24, Ampuero2011, Uphoff22], which are based on a quasi-static contact mechanics (without acceleration term) combined with a "radiation damping" stabilization term along the fault network. We will then consider the extension to elastodynamic models capable of modeling the propagation of seismic waves in the rock. The spatial discretization will be polytopal to cope with the complexity of geological meshes. It will build upon our previous work in contact mechanics [Droniou23] which combines a Virtual Element Method (VEM) nodal discretization of degree 1 for the displacement field, a facewise constant approximation of the traction vector stabilized by a face bubble enrichment of the displacement space. This method