Simulación numérica de terremotos en medios porosos: influencia del acoplamiento poroelástico en la ruptura dinámica y caracterización de la incertidumbre asociada

Las rupturas sísmicas en medios poroelásticos implican un conjunto de fenómenos complejos derivados de inestabilidades de fricción stick-slip y de los acoplamientos termohidromecánicos. En esta tesis, propongo un modelo de elementos finitos totalmente implícito, con un paso de tiempo adaptativo y monolíticamente acoplado para simular secuencias dinámicas de terremotos en medios poroviscoelásticos. Considero un material viscoelástico de Kelvin-Voigt y caracterizo el impacto de los efectos inerciales en los terremotos inducidos por inyección. En esta tesis, presento, por primera vez, simulaciones dinámicas de rupturas en fallas rate-and-state en medios poroelásticos. Mis simulaciones resuelven el ciclo sísmico completo incluyendo las fases intersísmica de nucleación espontánea del terremoto y de ruptura dinámica. Comparo las simulaciones dinámicas con las cuasi-dinámicas, en las que se desprecian los efectos inerciales y la singularidad del deslizamiento se resuelve mediante una aproximación de amortiguacin por radiación. La disipación viscosa modela el proceso físico de atenuación de las ondas sísmicas: a medida que aumenta el amortiguamiento viscoso, el tamaño del parche y el deslizamiento máximo de la falla se reducen, con lo que disminuye la magnitud esperada del terremoto. Desde un punto de vista computacional, la viscoelasticidad ayuda a evitar oscilaciones espurias de alta frecuencia durante la propagación de las ondas. Al incluir los efectos inerciales, el modelo dinámico tiene en cuenta las fluctuaciones transitorias en las presiones y en las tensiones del sólido durante la ruptura, que se desprecian en el enfoque cuasi-dinámico. La comprensión de estas perturbaciones transitorias puede arrojar luz sobre el papel de la presión de poro en el mecanismo del dynamic triggering. El enfoque dinámico poroviscoelástico es un buen compromiso entre el modelo sin viscosidad, totalmente dinámico, y el cuasi-dinámico. Una pequeña cantidad de amortiguación viscosa permite realizar cálculos más eficientes, preservando al mismo tiempo las características más relevantes de las rupturas dinámicas, en particular las velocidades de deslizamiento, el deslizamiento acumulado y el momento sísmico liberado. Además, la intensidad y daño potencial de los terremotos están relacionados con la velocidad de propagación de la ruptura a lo largo de las fallas de la corteza terrestre. La mayoría de los terremotos son sub-Rayleigh, con rupturas más lentas que las ondas de Rayleigh superficiales. En los terremotos supershear, la ruptura es más rápida que las ondas de corte, lo que da lugar a fuertes concentraciones de tensiones y a mayores intensidades en comparación con los más comunes sub- Rayleigh. A pesar de los importantes avances teóricos y experimentales de las dos últimas décadas, el control geológico y geomecánico de las transiciones de velocidad de ruptura sigue siendo poco conocido. En esta tesis, propongo que los fluidos de los poros desempeñan un papel importante en la explicación de la velocidad de ruptura de los terremotos: la presión de los poros puede aumentar bruscamente en el frente de compresión durante la propagación de la ruptura, promoviendo el fallo por cortante por delante del frente de ruptura y acelerando su propagación hacia el rango de las rupturas supershear. Caracterizo la transición de la ruptura sub-Rayleigh a la supershear en rocas saturadas, y muestro que el mecanismo de debilitamiento poroelástico propuesto puede ser un factor de control para las rupturas sísmicas intersónicas. Además, comprender el riesgo asociado a los terremotos antropogénicos es esencial para el desarrollo y la gestión de los procesos de ingeniería y las infraestructuras hidráulicas que pueden alterar las presiones y tensiones de los poros en profundidad. La posibilidad de que se produzcan terremotos provocados por el embalse, las mareas oceánicas y los fenómenos hidrológicos en la superficie de la Tierra (hidrosismicidad) ha sido objeto de un amplio debate. La relación entre la sismicidad inducida y los fenómenos hidrológicos se basa actualmente en correlaciones estadísticas mas que en mecanismos físicos. En este manuscrito, exploro las condiciones geomecánicas que podrían permitir que pequeños cambios de presión de poros debidos a la gestión de los yacimientos y a los cambios del nivel del mar se propaguen a profundidades compatibles con el desencadenamiento de terremotos en fallas sometidas a tensiones críticas (varios kilómetros). Considero una zona de falla dañada que está incrustada en una matriz de roca poroelástica, y realizo simulaciones hidromecánicas totalmente acopladas de la difusión de la presión y la deformación de la roca. Caracterizo las propiedades hidráulicas y geomecánicas de las zonas de falla que podrían permitir que pequeños cambios de presión y carga en la superficie del suelo (del orden de decenas o cientos de kPa) se propaguen con una atenuación relativamente pequeña hasta profundidades sismogénicas (hasta 10 km). Encuentro que la difusión de la presión a tales profundidades sólo es posible para zonas de falla altamente permeables y/o un fuerte acoplamiento poroelástico. ABSTRACT Earthquake ruptures in poroelastic media involve a suite of complex phenomena arising from stick-slip frictional instabilities and thermo-hydromechanical couplings. In this thesis, I propose a fully implicit, time-adaptive, and monolithically coupled finite element model to simulate dynamic earthquake sequences in poroviscoelastic media. I consider a Kelvin-Voigt viscoelastic material and characterize the impact of inertial effects on injection-induced earthquakes. In this thesis, I present, for the first time, dynamic simulations of ruptures in rate-and-state faults in poroelastic media. My simulations resolve the full earthquake cycle, including the interseismic, spontaneous earthquake nucleation, and dynamic rupture phases. I compare dynamic simulations with quasi-dynamic ones, in which inertial effects are neglected and the slip singularity is resolved through a radiation damping approximation. Viscous dissipation models the physical process of seismic wave attenuation: as viscous damping increases, the patch size and the maximum fault slip become smaller, hence decreasing the expected earthquake magnitude. From a computational perspective, viscoelasticity helps avoid spurious high-frequency oscillations during wave propagation. By including inertial effects, the dynamic model accounts for transient fluctuations of pressures and solid stresses during rupture, which are neglected in the quasi-dynamic approach. Understanding these transient perturbations may shed light on the role of pore pressure in the mechanism of dynamic earthquake triggering. The poroviscoelastic dynamic approach is a good compromise between the inviscid, fully dynamic model, and the quasi-dynamic one. A small amount of viscous damping allows us more efficient calculations, while preserving the most relevant features of dynamic ruptures, in particular slip velocities, accumulated slip, and seismic moment released. Moreover, the intensity and damage potential of earthquakes are linked to the speed at which rupture propagates along sliding crustal faults. Most earthquakes are sub-Rayleigh, with ruptures that are slower than the surface Rayleigh waves. In supershear earthquakes, rupture is faster than the shear waves, leading to sharp stress concentrations and larger intensities compared with the more common sub- Rayleigh ones. Despite significant theoretical and experimental advances over the past two decades, the geological and geomechanical control on rupture speed transitions remain poorly understood. In this thesis, I propose that pore fluids play an important role in explaining earthquake rupture speed: the pore pressure may increase sharply at the compressional front during rupture propagation, promoting shear failure ahead of the rupture front and accelerating its propagation into the supershear range. I characterize the transition from sub-Rayleigh to supershear rupture in fluidsaturated rock, and show that the proposed poroelastic weakening mechanism may be a controlling factor for intersonic earthquake ruptures. In addition, understanding the risk associated with anthropogenic earthquakes is essential in the development and management of engineering processes and hydraulic infrastructure that may alter pore pressures and stresses at depth. The possibility of earthquakes triggered by reservoir impoundment, ocean tides, and hydrological events at the Earth surface (hydro-seismicity) has been extensively debated. The link between induced seismicity and hydrological events is currently based on statistical correlations rather than on physical mechanisms. In this manuscript, I explore the geomechanical conditions that could allow for small pore pressure changes due to reservoir management and sea level changes to propagate to depths that are compatible with earthquake triggering at criticallystressed faults (several kilometers). I consider a damaged fault zone that is embedded in a poroelastic rock matrix, and conduct fully coupled hydromechanical simulations of pressure diffusion and rock deformation. I characterize the hydraulic and geomechanical properties of fault zones that could allow for small pressure and loading changes at the ground surface (in the order of tens or hundreds of kPa) to propagate with relatively small attenuation to seismogenic depths (up to 10 km). I find that pressure diffusion to such depths is only possible for highly permeable fault zones and/or strong poroelastic coupling.