Modèle mécanique non linéaire d’infiltration et analyse de bifurcation pour le couplage fluide-sol dans les massifs rocheux fracturés

Pourquoi l’eau cachée dans les roches fissurées compte

Bien en dessous de nos pieds, l’eau circule à travers des couches rocheuses brisées dans les mines, les tunnels et les réservoirs d’hydrocarbures. La façon dont cette eau agit sur la roche — et dont la roche réagit en retour — peut préparer discrètement le terrain à des inondations soudaines, des effondrements ou des séismes. Cet article élabore une nouvelle manière de décrire et de prédire cette interaction bidirectionnelle entre l’écoulement de l’eau et la déformation de la roche dans les massifs fracturés, révélant pourquoi ces systèmes peuvent paraître calmes pendant longtemps puis céder abruptement.

L’eau qui serpente à travers la pierre brisée

Dans de nombreux ouvrages souterrains, le massif rocheux n’est pas solide comme une brique mais ressemble plutôt à une éponge fissurée : d’innombrables fractures et vides forment des voies pour l’eau. La pression externe des couches sus-jacentes et la pression interne de l’eau dans les pores remodèlent en permanence ces chemins. Lorsque la roche est comprimée, les pores se réduisent ou se referment ; quand la pression de l’eau augmente, ils peuvent se rouvrir ou s’élargir. Cette reconfiguration continue modifie la perméabilité, qui rétroagit ensuite sur la roche. Les auteurs soutiennent que pour comprendre des catastrophes telles que les arrivées d’eau en mine ou les fuites de réservoirs, il faut traiter le phénomène comme un système dynamique et couplé, et non comme une image statique.

Construire une image couplée de la roche et de l’eau

L’étude commence par étendre un concept classique en mécanique des sols appelé « contrainte efficace », qui décrit quelle part de la pression totale est réellement portée par l’ossature solide. Les auteurs réécrivent cette idée pour inclure explicitement la porosité — la fraction du volume rocheux occupée par des vides — de sorte que les variations de l’espace poreux soient directement liées à la répartition de la contrainte entre la roche et l’eau. Ils combinent ensuite cela avec des équations décrivant la réponse élastique d’une roche légèrement déformable et avec une description plus réaliste et non linéaire de l’écoulement de l’eau dans les fractures, dépassant la simple loi de Darcy utilisée dans beaucoup de modèles d’ingénierie.

Du flux lisse au changement soudain

Avec ce cadre, les auteurs se concentrent sur un cas unidimensionnel d’eau s’infiltrant verticalement à travers une couche de roche brisée. Ils dérivent une paire d’équations non linéaires qui suivent l’évolution dans le temps et l’espace de la pression de l’eau et du débit, tandis que la porosité s’ajuste au fur et à mesure que la roche se compacte. La résolution de ces équations montre que, dans certaines conditions, le système n’a pas un comportement unique et stationnaire : il présente ce que les mathématiciens appellent une bifurcation de type nœud-col (saddle‑node). En termes simples, lorsqu’un paramètre clé du flux varie, un état auparavant stable peut se scinder en une branche sûre et une branche dangereuse, ou disparaître complètement, provoquant un saut brutal du régime d’infiltration doux vers un écoulement incontrôlé.

Compression lente et stabilité retardée

Les auteurs examinent ensuite l’effet de variations lentes de la contrainte aux limites — par exemple dues au chargement progressif lié à l’exploitation minière en surface. Les simulations numériques montrent que lorsque cette contrainte externe évolue lentement, le système couplé roche–eau met lui aussi beaucoup plus de temps à atteindre un état d’équilibre. La pression de l’eau, le débit et la déformation volumique de la roche migrent lentement vers la stabilité au lieu de se stabiliser rapidement. Ce retard s’explique par le fait que l’ossature rocheuse doit continuellement réajuster sa structure poreuse pendant que l’énergie est injectée progressivement par la charge variable, allongeant le chemin vers l’équilibre.

Signaux d’alerte avant une inondation

Pour relier la théorie à la réalité, l’étude compare ses prédictions à un cas réel d’émission d’eau depuis une faille dans une mine de charbon. À l’approche de la faille par l’exploitation, un paramètre traduisant l’écart par rapport au comportement de Darcy simple est entré dans une plage critique où deux états d’écoulement pouvaient coexister : l’un stable, l’autre instable. Les mesures de terrain ont montré la vitesse de l’eau commençant à osciller entre deux niveaux distincts avant de basculer ensuite dans une montée brutale et catastrophique, exactement comme le suggérerait le diagramme de bifurcation du modèle. Ces fluctuations, estiment les auteurs, constituent un signe d’alerte plus net et plus précoce d’une arrivée d’eau que les indicateurs de sécurité classiques qui traitent le système comme linéaire et stationnaire.

Implications pour la sécurité souterraine

En résumé, l’article montre qu’un massif rocheux fracturé et saturé d’eau se comporte davantage comme un système complexe et non linéaire que comme un simple conduit. De faibles variations de contrainte ou de conditions d’écoulement peuvent le pousser au-delà de seuils critiques où son comportement change qualitativement, pas seulement en amplitude. En reliant explicitement la déformation rocheuse, la structure poreuse et l’écoulement non linéaire, le nouveau modèle peut rendre compte de plusieurs états stationnaires possibles, de transitions soudaines entre eux et d’une forte sensibilité aux conditions initiales. Pour les ingénieurs concevant mines, tunnels et réservoirs, cela signifie que surveiller l’évolution conjointe des flux et de la déformation — et guetter les fluctuations caractéristiques à double stabilité — pourrait fournir des avertissements plus précoces et plus fiables d’instabilités cachées avant qu’elles ne dégénèrent en catastrophes.»

Citation: Zhengzheng, C., Mengqi, X., Tao, R. et al. Nonlinear seepage mechanical model and bifurcation analysis for fluid-solid coupling in fractured rock mass. Sci Rep 16, 9578 (2026). https://doi.org/10.1038/s41598-025-25823-6

Mots-clés: roche fracturée, infiltration des eaux souterraines, couplage fluide–sol, dynamique non linéaire, arrivée d’eau en mine