Revisiting the hydromechanical formulation of a micromechanics-based phase-field model for poro-elastoplastic media

Bien en dessous de la surface de la Terre, dans les roches poreuses qui renferment du pétrole, du gaz et de la chaleur géothermique, une bataille silencieuse se joue constamment entre la pierre solide et les fluides piégés dans ses minuscules pores. Lorsque les ingénieurs tentent d'extraire ces ressources, ils pompent souvent des fluides dans la roche pour l'écarter de force, un processus connu sous le nom de fracturation hydraulique. Pendant des décennies, les scientifiques ont modélisé ce processus en traitant la roche comme une simple éponge élastique qui s'étire et reprend sa forme. Cependant, les roches réelles sont plus complexes ; elles peuvent également se déformer de manière permanente, s'effriter et glisser le long de fissures internes, surtout lorsqu'elles sont comprimées par des pressions élevées. Comprendre exactement comment ces roches se brisent lorsque le squelette solide et le fluide à l'intérieur poussent et tirent simultanément est crucial pour forer des puits en toute sécurité, gérer le stockage des déchets et même prédire des événements géologiques naturels comme la fissuration des icebergs ou la poussée du magma à travers la croûte. Le défi consiste à créer une description mathématique capable de capturer la réalité désordonnée de la rupture de la roche sans se perdre dans des calculs impossibles.

Une équipe de chercheurs a revisité les équations fondamentales utilisées pour décrire ce processus de rupture, découvrant qu'un raccourci courant utilisé dans les modèles précédents conduit à une erreur significative. Dans le monde de la mécanique des roches, les scientifiques utilisent une approche par « champ de phase », qui traite une fissure non pas comme une ligne nette et invisible, mais comme une zone floue où le matériau s'affaiblit progressivement. Cela permet aux ordinateurs de simuler la croissance et le changement de forme d'une fissure. Pour ce faire, le modèle doit tenir compte de deux comportements distincts : lorsque la roche est étirée (tension) et lorsqu'elle est comprimée (compression). En tension, de minuscules fissures s'ouvrent ; en compression, elles se ferment et les surfaces rugueuses à l'intérieur glissent les unes contre les autres, créant de la friction. Les chercheurs ont découvert que de nombreux modèles existants, largement utilisés dans l'industrie et le milieu académique, confondent les variables utilisées pour décrire l'énergie du système. Plus précisément, ils ont constaté que traiter la pression du fluide comme une entrée directe dans l'équation d'énergie, plutôt que la quantité de fluide dans les pores, crée un « saut » ou une discontinuité mathématique. Ce saut signifie que le modèle prédit que la résistance de la roche change de manière abrupte et irréaliste lors du passage de la traction à la compression, conduisant à des prédictions inexactes sur la manière et l'endroit où une fracture va se former.

Pour corriger cela, les chercheurs ont développé une nouvelle formulation cohérente qui suit attentivement la relation entre la déformation de la roche et le contenu en fluide. Ils ont construit leur modèle sur une vue microscopique de la roche, l'imaginant comme une collection de minuscules fissures en forme de pièces de monnaie remplies de fluide. En analysant comment ces micro-fissures s'ouvrent et se ferment sous différentes contraintes, ils ont dérivé un ensemble de règles garantissant que la surface de résistance de la roche — la limite mathématique définissant quand elle se brisera — reste lisse et continue, quel que soit le changement de pression. Ils ont prouvé qu'en utilisant un type spécifique de règle de flux qui lie la déformation permanente de la roche directement aux forces qui agissent sur elle, et en incluant correctement un terme de couplage manquant dans l'équation d'énergie, le modèle peut décrire avec précision la transition entre la tension et la compression. Cette correction n'est pas un simple ajustement mineur ; elle modifie fondamentalement la façon dont le modèle calcule la force motrice qui pousse une fissure vers l'avant.

L'équipe a testé son nouveau modèle face à deux scénarios très différents pour voir s'il tenait la route. Premièrement, ils ont simulé un test de référence classique de fracturation hydraulique où un fluide est injecté dans une roche pour créer une fissure. Ils ont comparé leurs résultats à une solution analytique connue, qui est une réponse mathématique précise dérivée pour des conditions idéales. Le nouveau modèle correspondait beaucoup plus étroitement à cette solution que les anciennes formulations mixtes, prédisant correctement la pression nécessaire pour amorcer la fracture et la largeur de la fissure au fur et à mesure de sa croissance. Les anciens modèles avaient eu tendance à surestimer la pression de fluide requise, une faille qui pourrait conduire à des opérations de forage inefficaces ou dangereuses. Dans le second test, ils ont simulé un échantillon de roche étant compressé par deux côtés, une condition qui imite l'environnement souterrain profond où les roches sont soumises à des pressions immenses. Ici, la roche ne s'est pas contentée de se briser ; elle s'est déformée plastiquement, formant une bande de cisaillement où le matériau a glissé sur lui-même avant de se rompre. Le nouveau modèle a reproduit avec succès ce comportement complexe, capturant à la fois les fractures de cisaillement causées par la compression mécanique et les fractures de tension entraînées par l'injection de fluide.

Les résultats de ces simulations révèlent que la déformation plastique de la roche agit comme une barrière significative à la croissance de la fracture. Lorsque la roche est autorisée à se déformer de manière permanente, elle absorbe plus d'énergie, ce qui signifie que la pression du fluide à l'intérieur de la fissure doit être plus élevée pour maintenir la fracture ouverte et en croissance. Cette résistance supplémentaire change tout le profil de la fracture, la rendant plus étroite et nécessitant plus de force pour se propager qu'un modèle purement élastique ne le suggérerait. De plus, dans les tests de compression, le modèle a montré comment la pression interstitielle au sein de la roche change à mesure que la fracture se forme. À mesure que la roche se déforme et crée de l'espace pour le mouvement du fluide, la pression chute localement, attirant davantage de fluide dans la zone endommagée. Cette interaction crée une boucle de rétroaction où le flux de fluide et la défaillance mécanique de la roche sont étroitement couplés, une dynamique que le nouveau modèle capture avec une grande fidélité.

En corrigeant le fondement mathématique de ces modèles, les chercheurs ont fourni un outil plus fiable pour les ingénieurs et les scientifiques travaillant avec des milieux poreux. La nouvelle formulation garantit que la transition entre les différents types de rupture des roches est fluide et physiquement réaliste, éliminant les sauts artificiels qui ont entravé les tentatives précédentes. Cette précision est vitale pour des applications allant de l'optimisation de la production de pétrole et de gaz à la conception d'un stockage souterrain sûr pour la capture du carbone. L'étude démontre que même dans des systèmes complexes où le solide et le fluide interagissent, une approche rigoureuse et cohérente de la physique sous-jacente peut révéler la véritable nature de la rupture des matériaux. Le modèle est désormais prêt à aider à prédire le comportement de la croûte terrestre avec une plus grande confiance, comblant le fossé entre la mécanique théorique et la réalité plastique et désordonnée du sol sous nos pieds.