A Nonmatching-Grid Discontinuous Galerkin Method for Nonlinear Reactive Transport in Hybrid-Dimensional Fractured Media

Dans le monde caché sous nos pieds, l'eau et les produits chimiques dissous circulent à travers un labyrinthe complexe de roche. Cette roche est rarement uniforme ; il s'agit souvent d'un bloc solide, appelé matrice, sillonné par de fines fissures appelées fractures. Ces fractures sont les autoroutes du sous-sol. Elles peuvent agir comme des canaux ultra-rapides qui font défiler les fluides à travers la roche, ou comme des barrières qui bloquent totalement le mouvement, selon leur nature. Comme ces fissures sont souvent de largeur microscopique mais peuvent s'étendre sur des kilomètres de long, modéliser le déplacement des substances à travers elles représente un défi colossal pour les scientifiques. S'ils tentent de construire un modèle informatique qui capture chaque détail infime de l'ouverture de la fissure, le calcul devient si lourd que même les supercalculateurs les plus puissants ne peuvent le résoudre. Pour contourner cela, les chercheurs utilisent un raccourci ingénieux appelé modèle hybride-dimensionnel. Au lieu de traiter la fissure comme une tranche tridimensionnelle mince, ils la réduisent en une ligne bidimensionnelle ou une ligne unidimensionnelle située à l'intérieur de la roche. Cela permet de simuler l'écoulement sans avoir besoin de résoudre l'épaisseur minuscule de la fissure elle-même.

La difficulté surgit lorsque les scientifiques tentent de relier les mathématiques de la roche solide avec celles de ces fissures aplaties. Dans le monde réel, la roche et la fissure sont en constante conversation, échangeant des fluides et des produits chimiques à travers leur frontière commune. Pour simuler cela, les ordinateurs divisent généralement la roche et la fissure en petites pièces, comme une mosaïque, afin de calculer l'écoulement. Pendant des années, l'approche standard exigeait que les pièces de la roche et les pièces de la fissure s'alignent parfaitement le long de leur bord partagé, comme deux pièces de puzzle s'emboîtant l'une dans l'autre. Si la roche était divisée en grands carrés et la fissure en segments minuscules, l'ordinateur devait les forcer à correspondre, souvent en créant une couche de données complexe et artificielle pour combler l'écart. Cette exigence rendait les modèles rigides et difficiles à utiliser lorsqu'il fallait affiner différentes parties d'un système à des échelles différentes.

Une nouvelle étude d'Alvaro H. Salas et Carlos D. Acosta, de l'Universidad Nacional de Colombia, s'attaque à cette rigidité en posant une question simple mais profonde : et si la roche et la fissure n'avaient pas besoin de s'aligner du tout ? Ils ont développé une méthode où la roche et la fracture peuvent être divisées en pièces de tailles et de formes complètement différentes, générées indépendamment l'une de l'autre, et pourtant échanger des informations parfaitement. Imaginez que vous essayiez de verser de l'eau d'un grand seau grossier dans un tamis fin et détaillé. Habituellement, vous auriez besoin d'un entonnoir spécial pour que l'eau soit transférée sans déversement. Les chercheurs ont trouvé un moyen de verser l'eau directement, en calculant la quantité exacte qui passe par chaque point où le bord du seau touche le tamis, quels que soient le décalage ou l'inadéquation des trous.

L'équipe a créé un cadre mathématique qui traite l'échange de produits chimiques entre la roche et la fracture comme une interaction physique directe. Au lieu de forcer les deux grilles à correspondre ou d'utiliser un intermédiaire pour traduire les données, ils ont construit un système qui calcule le transfert sur les sections minuscules et chevauchantes où le bord de la roche et la surface de la fracture se rencontrent réellement. Cette approche est appelée une méthode de grille non appariée (nonmatching-grid). Dans leurs simulations, ils ont testé ce système avec des réactions chimiques complexes, où les substances changent au fur et à mesure qu'elles se déplacent, un scénario notoirement difficile à modéliser. Ils ont prouvé mathématiquement que cette connexion directe est stable et précise, même lorsque la roche est divisée en très gros blocs et la fracture en segments très fins, ou inversement. La stabilité de leur méthode ne dépend pas du rapport de taille de ces blocs, ce qui signifie que l'ordinateur ne plante pas ou ne produit pas d'erreurs aberrantes simplement parce qu'une grille est beaucoup plus fine que l'autre.

L'une des découvertes les plus significatives de leur travail est que cette méthode directe est supérieure aux anciennes approches indirectes qui tentent de projeter les données d'une grille sur l'autre. Lorsque l'échange de produits chimiques est fort, les anciennes méthodes de projection perdent de l'information. Elles lissent essentiellement les détails de la fracture, écartant les variations à petite échelle qui sont cruciales pour une prédiction précise. Les chercheurs ont démontré que leur méthode directe conserve toutes les informations nécessaires. Lors de tests où le transfert de produits chimiques était intense, l'ancienne méthode produisait des erreurs vingt et un pour cent plus grandes que leur nouvelle méthode. De plus, lorsqu'ils affinaient la grille de la fracture pour la rendre encore plus détaillée tout en maintenant la grille de la roche fixe, l'ancienne méthode stagnait. Leur nouvelle méthode, cependant, continuait de devenir plus précise, prouvant qu'elle capture la véritable physique de l'interaction sans limitations artificielles.

Les chercheurs ont également montré que cette méthode fonctionne efficacement sur les ordinateurs. Ils ont écrit un programme informatique pour tester leur théorie, effectuant des simulations sur un ordinateur de bureau standard. Ils ont comparé le temps nécessaire pour préparer les calculs avec des grilles correspondantes par rapport aux grilles non appariées. Les résultats ont montré que le travail supplémentaire requis pour gérer les pièces dépareillées était négligeable. Le temps nécessaire pour assembler les équations était presque identique, ce qui signifie que les scientifiques peuvent désormais utiliser des grilles beaucoup plus flexibles et réalistes sans payer de pénalité en temps de calcul. Cette flexibilité est vitale pour les applications réelles, où différentes parties d'une formation géologique pourraient être modélisées par différentes équipes utilisant des logiciels différents, ou là où le réseau de fractures est si complexe qu'une grille unique et parfaitement alignée est impossible à créer.

L'étude a également exploré comment la méthode se comporte dans des conditions extrêmes. Ils ont testé des scénarios où les réactions chimiques étaient très fortes et le transfert de fluides était rapide. Dans ces cas difficiles, la méthode est restée robuste, convergeant vers une solution sans échouer. Ils ont même testé la géométrie des fractures, simulant des fissures qui se courbent ou s'intersectent entre elles. Bien que les tests principaux aient utilisé une fissure droite, la logique sous-jacente de leur méthode a tenu bon lorsqu'elle a été appliquée à ces formes plus complexes, suggérant que l'approche pourrait être étendue aux réseaux emmêlés que l'on trouve dans la nature. Les chercheurs ont pris soin de noter que leur travail actuel se concentre sur des types spécifiques de roches et de géométries de fractures, et que l'extension de la méthode à des réseaux tridimensionnels complets, courbes et intersectés est une tâche pour l'avenir.

Ce qui rend ce travail particulièrement important, c'est qu'il lève une barrière de longue date dans la simulation des milieux fracturés. Pendant des décennies, la nécessité de grilles correspondantes a limité la capacité des scientifiques à modéliser ces systèmes, les forçant souvent à faire des compromis sur le détail de la fracture ou de la roche pour que les mathématiques fonctionnent. En prouvant que les deux peuvent être couplés directement et de manière conservatrice sans correspondance parfaite, les chercheurs ont ouvert la voie à des simulations plus réalistes et plus flexibles. Cela signifie qu'à l'avenir, les modèles de contamination des eaux souterraines, de stockage de carbone ou d'extraction d'énergie géothermique pourront être construits avec plus de liberté, permettant à la roche et aux fissures d'être représentées exactement telles qu'elles sont, plutôt que telles qu'elles doivent être pour s'adapter à une contrainte mathématique rigide. La méthode fournit une base stable et nouvelle pour comprendre comment les fluides et les produits chimiques circulent à travers la terre fracturée, garantissant que la connexion entre la roche solide et les fissures est aussi fluide dans l'ordinateur que dans le sol.