@article{fdi:010086975,
  title = {{F}luid flow migration, rock stress and deformation due to a crustal fault slip in a geothermal system : a poro-elasto-plastic perspective},
  author = {{S}aez-{L}eiva, {F}. and {H}urtado, {D}. {E}. and {G}erbault, {M}uriel and {R}uz-{G}inouves, {J}. and {I}turrieta, {P}. and {C}embrano, {J}.},
  editor = {},
  language = {{ENG}},
  abstract = {{G}eothermal systems are commonly genetically and spatially associated with volcanic complexes, which in turn, are located nearby crustal fault systems. {F}aults can alter fluid flow in their surroundings, potentially acting as barriers or conduits for fluids, depending on their architecture and slip-rate. {H}owever, this fundamental control on fluid migration is still poorly constrained. {M}ost previous modeling efforts on volcanic and hydrothermal processes consider either only fluid flow in their formulations, or only a mechanical approach, and seldom a full, monolithic coupling between both. {I}n this work, we present a poro-elasto-plastic {F}inite {E}lement {M}ethod ({FEM}) to address the first-order, time-dependent control that a strike-slip crustal fault exerts on a nearby geothermal reservoir. {F}or the model setting, we selected the {P}lanchon-{P}eteroa geothermal system in the {S}outhern {A}ndes {V}olcanic {Z}one ({SAVZ}), for which the geometry and kinematics of a potentially seismogenic fault and fluid reservoir is constrained from previous geological and geophysical studies. {W}e assess the emergence and diffusion of fluid pressure domains due to fault slip, as well as the development of tensile/dilational and compressive/contractional domains in the fault' surroundings. {M}ean stress and volumetric strain magnitudes in these domains range between +/- 1 [{MP}a] and +/- 10-4 [-], respectively. {O}ur results show the appearance of negative and positive fluid pressure domains in these dilational and contractional regions, respectively. {W}e also investigate the spatial and temporal evolution of such domains resulting from changes in fault permeability and shear modulus, fluid viscosity, and rock rheology. {T}hese variations in fluid pressure alter the trajectory of the reservoir fluids, increasing migration to the eastern half of the fault, reaching a maximum fluid flux of 8 to 70 times the stationary flux. {P}ressure-driven fluid diffusion over time causes fluid flow to return to the stationary state between weeks to months after fault slip. {T}hese results suggest that the mechanism that exerts a first-order control is similar to a suction pump, whose duration heavily depends on fault permeability and fluid viscosity. {W}e also show how a von {M}ises plasticity criterion locally enhances fluid flow. {T}he transient process analyzed in this work highlights the importance of addressing the solid-fluid coupling in numerical models for volcano-tectonic studies.},
  keywords = {finite element method ; poromechanics ; geothermal system ; fault zone ; suction pump ; {CHILI} ; {ANDES}},
  booktitle = {},
  journal = {{E}arth and {P}lanetary {S}cience {L}etters},
  volume = {604},
  numero = {},
  pages = {117994 [15 ]},
  ISSN = {0012-821{X}},
  year = {2023},
  DOI = {10.1016/j.epsl.2023.117994},
  URL = {https://www.documentation.ird.fr/hor/fdi:010086975},
}