Lattice-Boltzmann pour milieux poreux : plus de 100 millions d’heures GPU

Pourquoi les espaces minuscules dans la roche comptent

Lorsque nous stockons de l’énergie sous terre, capturons du dioxyde de carbone ou faisons fonctionner des piles à combustible, la manière dont deux fluides distincts s’entrelacent dans les micro‑espaces de la roche peut déterminer le succès de la technologie. Or, observer ces fluides dans des roches réelles est extrêmement lent et coûteux. Cet article présente un vaste jeu de données ouvert issu de simulations informatiques qui font l’étape suivante : elles recréent la façon dont l’eau et des fluides proches du pétrole se déplacent l’un par rapport à l’autre dans des échantillons de roche réalistes, en utilisant plus de 100 millions d’heures de calcul sur processeurs graphiques. Le résultat est une ressource partagée que n’importe quel chercheur peut utiliser pour tester des idées sur les écoulements souterrains sans disposer d’un accélérateur de particules ou d’un superordinateur personnel.

Explorer le mouvement des fluides dans des roches numériques

Les auteurs se concentrent sur l’« écoulement diphasique », où deux fluides non miscibles, comme l’eau et le pétrole, traversent ensemble un labyrinthe de pores dans la roche. Une grandeur clé pour les ingénieurs est la perméabilité relative, qui indique la facilité de déplacement de chaque fluide en présence de l’autre. Normalement, des mesures complètes exigent des semaines de travail de laboratoire pour chaque échantillon et chaque condition. À la place, l’équipe a utilisé un logiciel de simulation spécialisé, LBPM, pour calculer les écoulements directement sur des images 3D de roches réelles. Ces roches numériques proviennent de tomographies par rayons X à micro‑échelle de grès et de verre fritté, capturant des formes et des tailles de pores réalistes jusqu’à quelques micromètres.

Une expérience virtuelle massive

Réaliser des simulations réalistes à cette échelle fine demeure très coûteux. L’équipe a mobilisé des ordinateurs haute performance dotés de milliers de GPU pour balayer des conditions impraticables en laboratoire. Ils ont fait varier la préférence de la surface rocheuse pour un fluide ou l’autre (son comportement de mouillage) et la force avec laquelle les fluides sont poussés (capturée par un paramètre appelé nombre capillaire). Pour quatre matériaux poreux différents, ils ont exécuté des protocoles d’écoulement « en régime permanent » et « hors régime » qui imitent les essais standards de pompage de carottes utilisés dans l’industrie pétrolière et gazière. Au total, ils ont produit 50 courbes complètes de perméabilité relative et plus de 25 000 configurations fluides distinctes.

Voir les formes, pas seulement les moyennes

Au‑delà des débits moyens, les simulations suivent les formes détaillées et les connexions des fluides au fil du temps — des informations presque impossibles à obtenir à grande échelle par expériences. LBPM sépare les voies d’écoulement connectées des gouttelettes isolées, parfois appelées ganglions, et mesure leur volume, surface, courbure, connectivité, pression et mouvement. Ces grandeurs sont consignées à chaque pas de simulation dans des tables textuelles simples, de sorte que les utilisateurs peuvent reconstruire la formation de poches piégées, leur rupture, leur reconnexion ou leur drainage lent. En comparant différents types de roche et schémas de mouillage, le jeu de données révèle comment des variations subtiles de préférence de surface peuvent modifier l’emplacement des zones piégées et la facilité de déplacement des fluides, aidant à expliquer les tendances observées en laboratoire.

Enseignements sur la roche mouillante et les fluides piégés

À l’aide du jeu de données, les auteurs valident que leurs simulations se comportent de manière cohérente en vérifiant des schémas connus. Par exemple, lorsque la roche devient plus attirée par le pétrole, la quantité de pétrole résiduelle après une purge à l’eau a tendance à diminuer, en accord avec des expériences antérieures. Dans des cas plus hydrophiles avec des goulots d’étranglement mal résolus, les simulations montrent une fragmentation précoce du fluide non mouillant et des courbes d’écoulement anormalement aplaties, illustrant comment la résolution d’image peut induire en erreur. Dans d’autres roches aux pores mieux résolus, le déplacement des courbes d’écoulement avec le changement de mouillage s’aligne sur des études de laboratoire antérieures. Des mesures de connectivité confirment que, dans certaines conditions intermédiaires, les deux fluides restent fortement interconnectés, un état associé à des interfaces complexes et à des amas piégés de longue durée.

Une base partagée pour les travaux futurs

En termes simples, cet article fournit une carte détaillée de la manière dont deux fluides peuvent partager les mêmes micro‑espaces dans la roche, sous une grande variété de conditions, le tout encodé dans des données réutilisables. Les simulations sont comparées à des expériences d’imagerie par rayons X de pointe et organisées pour que d’autres puissent facilement calculer de nouvelles courbes récapitulatives, entraîner des modèles d’apprentissage automatique ou tester de nouvelles théories sur la façon dont les détails à l’échelle des pores se traduisent en écoulement à grande échelle. Pour quiconque travaille sur le stockage géologique du carbone, le stockage d’hydrogène, la dépollution des nappes phréatiques ou les composants de piles à combustible, ce jeu de données ouvert offre un raccourci puissant : au lieu de repartir de zéro, il permet de s’appuyer directement sur plus de 100 millions d’heures GPU d’expériences numériques soigneusement organisées.

Citation: Armstrong, R.T., Tavakkoli, O., Da Wang, Y. et al. Lattice-Boltzmann for Porous Media: 100M⁺ GPU Hours. Sci Data 13, 697 (2026). https://doi.org/10.1038/s41597-026-06823-1

Mots-clés: milieux poreux, écoulement diphasique, roche numérique, lattice boltzmann, perméabilité relative