Injection-rate effects on failure in a fluid-saturated granular fault gouge

Cet article combine la théorie analytique et les simulations numériques pour démontrer que le taux d'injection de fluide régit la rupture de la gouge de faille en créant une hétérogénéité de pression, où une injection lente provoque un affaiblissement uniforme tandis qu'une injection rapide préserve la résistance dans les régions distales, offrant ainsi un cadre raffiné pour prédire la sismicité lors des opérations géotechniques.

L'essentiel

Imaginez que la croûte terrestre soit comme un immense trottoir fissuré. À l'intérieur de ces fissures, il n'y a pas seulement du vide ; elles sont remplies de roche broyée, de sable et de terre. Les géologues appellent cela la « gouge de faille ». Maintenant, imaginez que quelqu'un commence à pomper de l'eau dans cette fissure. C'est une pratique courante pour des choses comme l'énergie géothermique ou l'élimination des déchets, mais cela comporte un risque : la pression de l'eau peut pousser la fissure à s'ouvrir, provoquant le glissement du sol et déclenchant un séisme.

Cet article pose une question simple mais complexe : Est-ce que la vitesse à laquelle vous pompez l'eau est importante ?

Les chercheurs ont découvert que, oui, cela compte énormément. Si vous pompez lentement, l'eau se répartit uniformément et la fissure glisse facilement. Mais si vous pompez très vite, vous avez en réalité besoin de plus de pression pour faire glisser la fissure.

Voici comment ils ont découvert cela, en utilisant un mélange de mathématiques et de simulations informatiques.

L'analogie de la « pièce bondée »

Considérez la gouge de la faille (la roche broyée) comme une pièce bondée remplie de gens (les grains).

• L'objectif : Vous voulez que tout le monde se déplace vers le côté (glissement).
• L'eau : La pression de l'eau est comme une « poussée » qui tente de faire écarter les gens.
• L'expansion : À mesure que les gens tentent de bouger, ils s'espacent naturellement et occupent plus d'espace (c'est ce qu'on appelle la « dilatation »).

Les deux scénarios

1. Le versement lent (Injection lente)
Imaginez verser de l'eau dans la pièce très lentement. L'eau a tout le temps de s'infiltrer à travers la foule et d'atteindre chaque personne. La pression devient uniforme. Tout le monde ressent la poussée en même temps, toute la foule se relâche ensemble, et la pièce glisse facilement. C'est ce que prédisaient les anciennes théories : plus de pression d'eau = glissement plus facile.

2. Le jet d'eau haute pression (Injection rapide)
Maintenant, imaginez projeter de l'eau dans la pièce d'un côté à grande vitesse, comme avec une lance à incendie.

• Le gradient : Les gens juste à côté de la lance sont immédiatement trempés et poussés. Mais les personnes à l'autre bout de la pièce ? Elles sont toujours au sec et tiennent bon.
• Le goulot d'étranglement : Même si les gens près de la lance sont prêts à bouger, les personnes à l'autre extrémité maintiennent la ligne. Toute la pièce ne peut pas glisser tant que les personnes les plus éloignées n'ont pas été poussées avec force.
• Le résultat : Pour faire glisser toute la pièce, vous devez augmenter la pression au niveau de la lance à un niveau bien plus élevé que si vous aviez fait un versement lent. L'injection rapide crée un « gradient de pression » où la poussée est forte à un endroit et faible à un autre.

L'effet du « sable meuble »

Il y a un second rebondissement. À mesure que les grains de roche tentent de glisser, ils ne font pas que glisser ; ils s'entremêlent et se réorganisent, rendant la couche légèrement plus épaisse (dilatation).

• Dans les simulations informatiques, les chercheurs ont augmenté la pression par paliers, laissant les grains se stabiliser après chaque étape.
• Ils ont découvert qu'à mesure que les grains se réorganisent et que la couche devient plus « meuble », le matériau devient en réalité plus faible.
• Cependant, parce que l'injection rapide crée ces zones de pression inégales (fortes près de la lance, faibles loin de celle-ci), la « faiblesse » du sable meuble n'aide l'ensemble de la faille à glisser que lorsque la pression est assez élevée pour surmonter les zones sèches et solides à l'extrémité opposée.

La « recette » mathématique

Les auteurs ont créé une nouvelle formule mathématique pour prédire exactement quand la faille va glisser. Leur formule indique que la pression nécessaire pour provoquer un glissement dépend de trois facteurs :

1. La vitesse de pompage : Un pompage plus rapide = une pression plus élevée nécessaire.
2. La longueur de la faille : Des failles plus longues = une pression beaucoup plus élevée nécessaire (car la pression doit voyager plus loin pour atteindre les zones faibles).
3. La vitesse à laquelle l'eau circule à travers la roche : Si la roche est très poreuse (l'eau circule vite), la pression s'égalise rapidement et vous avez besoin de moins de pression.

Pourquoi cela importe (selon l'article)

L'article conclut que nous ne pouvons pas simplement utiliser les anciennes règles simples qui supposent que la pression de l'eau est la même partout.

• Le compromis : Si vous injectez du fluide très rapidement, vous pourriez avoir besoin d'une pression plus élevée pour déclencher un glissement, ce qui semble plus sûr. Cependant, comme vous pompez rapidement, vous pourriez atteindre ce seuil de pression dangereux plus tôt dans le temps.
• Le facteur de taille : La longueur de la faille est un facteur énorme. Une failte courte se comporte comme le « versement lent » (pression uniforme), mais une faille longue se comporte comme le « jet d'eau » (pression inégale), ce qui rend la prédiction de son glissement beaucoup plus difficile.

En résumé, l'article montre que la vitesse change les règles. Pomper rapidement crée des zones de pression inégales qui agissent comme un « mode de maintien », nécessitant une force nettement plus grande pour briser la faille que de pomper lentement. Cela aide les ingénieurs à comprendre que la taille de la faille et la vitesse d'injection sont des facteurs critiques pour prévenir ou prédire les tremblements de terre induits.

Résumé technique

Résumé technique : Effets du taux d'injection sur la rupture dans une gouge de faille saturée en fluide

Énoncé du problème
L'injection de fluides dans le sous-sol pour l'extraction d'énergie, l'élimination des déchets et le développement des ressources est un moteur connu de la sismicité induite, réactivant souvent les gouges de failles remplies. Bien que la réduction de la contrainte normale effective via la diffusion de la pression interstitielle soit le mécanisme standard de l'affaiblissement des failles (loi de la contrainte effective de Terzaghi), le rôle du taux d'injection dans le contrôle de la rupture reste mal compris. Les modèles classiques supposent souvent une pression interstitielle uniforme, mais des études récentes indiquent que des distributions de pression hétérogènes peuvent modifier considérablement les seuils de rupture. Plus précisément, une injection rapide peut créer de forts gradients de pression, laissant les régions distales de la faille fortes et retardant la rupture, tandis qu'une injection lente permet un affaiblissement uniforme. De plus, contrairement aux failles à surfaces nues où la connectivité hydraulique prédomine, les gouges de failles impliquent des réponses granulaires complexes, notamment la dilatation (expansion volumique) lors de la phase de pré-rupture. Cette dilatation a des effets divergents : elle peut provoquer des chutes de pression interstitielle transitoires (durcissement dilatant) ou, dans le contexte d'un réarrangement progressif des grains, conduire à une réduction de la densité de compactage et de la résistance au cisaillement (affaiblissement). L'interaction entre le taux d'injection, la diffusion de la pression et la dilatation granulaire dans une gouge saturée en fluide nécessite un cadre quantitatif.

Méthodologie
Les auteurs emploient une double approche combinant théorie analytique et simulations numériques couplées :

1. Simulations numériques : Un code solide-fluide couplé est utilisé, intégrant une méthode des éléments discrets (DEM) pour les grains solides (modélisés comme des disques 2D avec des contacts frictionnels élastiques linéaires) et un solveur de pore-fluide continu (différences finies sur une grille eulérienne). Le dispositif simule une couche granulaire pré-contrainte et saturée en fluide représentant une gouge de faille, confinée entre des parois imperméables. Le fluide est injecté aux limites latérales à des taux contrôlés (20–200 kPa/min) par paliers de pression discrets pour permettre l'équilibrage entre pressurisation et dilatation, évitant ainsi les chutes de pression transitoires. Les simulations sont exécutées selon divers niveaux de pré-contrainte ($\tau/\sigma_n$) et de longueurs de couche normalisées ($L/d$).
2. Théorie analytique : Un modèle analytique est dérivé d'une équation de diffusion de la pression interstitielle qui inclut un terme de « puits dilatant ». Ce terme rend compte de l'expansion volumique du squelette granulaire à l'approche de la rupture. Le modèle suppose une relation linéaire entre le taux de dilatation et le taux de pressurisation. La solution de cette équation de diffusion avec des conditions aux limites spatialement uniformes et un terme de puits prédit l'évolution des profils de pression à travers la faille.

Résultats clés

• Dépendance de la pression de rupture vis-à-vis du taux : Les simulations numériques démontrent que la pression de bordure requise pour déclencher la rupture de la faille ($P^b_{fail}$) augmente systématiquement avec le taux d'injection. Cela confirme qu'une injection rapide retarde la rupture par rapport à une injection lente.
• Distribution de pression hétérogène : La solution analytique révèle que l'injection rapide génère de forts gradients de pression spatiaux. La pression est la plus élevée près des limites d'injection et la plus faible au centre de la faille ($x=L/2$). Par conséquent, la rupture globale est contrôlée non pas par la pression moyenne ou la pression de bordure, mais par la région la moins pressurisée (la plus forte) de la couche.
• Lois d'échelle : Le modèle analytique prédit que la pression de rupture suit une mise à l'échelle linéaire avec le taux d'injection ($\dot{P}_b$) et une mise à l'échelle quadratique avec la longueur de la faille ($L^2$), tout en variant inversement avec la diffusivité hydraulique ($\alpha$) et en dépendant du produit du module de compressibilité effectif et de la capacité de stockage du fluide ($K\beta$).
• Rôle de la dilatation : Les simulations montrent que la dilatation augmente avec la pression, suivant une mise à l'échelle en loi de puissance ($\Delta h/h_0 \sim (\Delta P^b/\sigma_n)^{3/2}$). Le module de compressibilité effectif ($K$) diminue à mesure que la couche se dilate, reflétant l'affaiblissement progressif du compactage granulaire. Le modèle analytique limite les résultats des simulations en traitant $K$ comme une constante effective, bien que les simulations présentent un comportement sous-linéaire attribué à la réduction progressive du coefficient de friction à mesure que la dilatation augmente.
• Validation : Les prédictions analytiques encadrent avec succès les résultats numériques sur trois ordres de grandeur de rigidité du compactage des grains ($K\beta$). Le modèle reproduit les observations expérimentales que la théorie classique de la contrainte effective à pression uniforme ne parvient pas à capturer.

Signification et affirmations
L'article affirme fournir un cadre quantitatif reliant la physique à l'échelle du grain (dilatation granulaire et réarrangement) à la rupture de la faille à l'échelle de la faille sous injection de fluide. En dérivant un critère de rupture dépendant du taux, l'étude explique pourquoi des taux d'injection plus élevés nécessitent des pressions plus élevées pour déclencher la réactivation dans les failles remplies de gouge.

Les auteurs soulignent que cette dépendance au taux provient de l'hétérogénéité de la pression interstitielle contrôlée par la diffusion. Dans les scénarios d'injection lente, la pression s'équilibre, conduisant à un affaiblissement uniforme. Dans les scénarios d'injection rapide, de forts gradients persistent, et la faille reste stable jusqu'à ce que la région centrale, la moins pressurisée, atteigne le seuil de rupture.

Concernant les implications pratiques, l'article suggère que les opérateurs pourraient potentiellement gérer les risques de réactivation en contrôlant les protocoles d'injection. Bien qu'une injection plus rapide augmente le seuil de pression de rupture, elle réduit également le temps avant la rupture. Les auteurs notent que ce compromis entre taux et pression pourrait être avantageux dans certains contextes opérationnels (par exemple, l'injection cyclique), mais ils précisent explicitement que les applications sur le terrain doivent tenir compte de l'hétérogénéité naturelle des failles, des géométries complexes et des effets thermiques, qui sont simplifiés ou négligés dans leur modèle 2D. L'étude conclut que les critères de rupture dérivés des échelles de laboratoire ne peuvent pas être directement extrapolés aux échelles de terrain sans tenir compte des effets d'échelle contrôlés par la diffusion, particulièrement en ce qui concerne la longueur de la faille.