Immiscible two-phase flow in porous media: a statistical mechanics approach

Cet article présente une approche de mécanique statistique, fondée sur l'entropie informationnelle de Shannon, pour décrire l'écoulement diphasique immiscible dans les milieux poreux à l'échelle du continuum, offrant ainsi une formalisation thermodynamique gérable qui intègre la physique à l'échelle des pores et introduit de nouvelles variables émergentes comme la vitesse co-mobile.

L'essentiel

🌊 Le Grand Défi : Faire passer deux liquides qui ne se mélangent pas dans une éponge

Imaginez que vous essayez de faire couler de l'eau et de l'huile à travers une éponge très serrée. C'est un peu comme essayer de faire passer deux équipes de coureurs différents sur un sentier de montagne très étroit et rocailleux. L'eau veut rester collée aux rochers (c'est le fluide "mouillant"), tandis que l'huile essaie de passer au milieu (le fluide "non mouillant").

Le problème principal que les physiciens tentent de résoudre est le suivant : Comment décrire ce chaos microscopique à grande échelle ?

• À l'échelle de l'éponge (microscopique) : C'est un cauchemar de calculs. Les gouttes d'huile se cassent, se rejoignent, bloquent les trous, etc. C'est trop complexe pour faire des prédictions pratiques.
• À l'échelle du réservoir (macroscopique) : Les ingénieurs ont besoin d'une règle simple, comme une carte routière, pour dire : "Si je pousse avec telle force, l'eau sortira à telle vitesse".

Pendant 90 ans, les scientifiques ont utilisé une "théorie de la perméabilité relative". C'est un peu comme une vieille carte routière dessinée à la main : elle fonctionne souvent, mais elle est pleine de trous et de suppositions hasardeuses. Elle ne nous dit pas pourquoi ça marche, juste que ça marche.

🧠 La Révolution : Utiliser les maths de l'information (comme un détective)

C'est ici que l'article propose une idée géniale. Au lieu de regarder chaque goutte d'eau individuellement, les auteurs (Alex Hansen et Santanu Sinha) utilisent une méthode inventée dans les années 50 par un génie nommé Edwin Jaynes.

Imaginez que vous êtes un détective dans une foule. Vous ne pouvez pas connaître la position exacte de chaque personne (trop de détails !). Mais vous pouvez connaître des statistiques : "Combien de gens sont à gauche ?", "Combien sont à droite ?", "Quelle est la vitesse moyenne ?".

Jaynes a dit : "Si on ne connaît pas tout, on doit utiliser la théorie de l'information pour deviner la situation la plus probable." C'est comme si on utilisait la mécanique statistique (la physique des gaz et des atomes) pour étudier des fluides qui ne sont pas chauffés par la température, mais par le mouvement et la pression.

🎭 Le Nouveau Langage : Une "Thermodynamique" sans Chaleur

En appliquant cette méthode à l'éponge, les chercheurs ont découvert quelque chose de surprenant : le flux de ces deux liquides obéit à des règles qui ressemblent étrangement à la thermodynamique (la science de la chaleur et de l'énergie), mais sans la chaleur !

Ils ont créé de nouveaux mots pour décrire ce monde :

1. L'Agitation (Agiture) : Au lieu de la température (qui mesure l'agitation des atomes), ils utilisent l'agitation pour mesurer l'intensité du flux. Plus vous poussez fort (plus la pression est forte), plus le système est "agité". C'est comme si l'éponge devenait plus "chaude" d'un point de vue mouvement.
2. La Vitesse Co-mobile (Co-moving velocity) : C'est la grande découverte de l'article.

🚗 L'Analogie de la Route : La Vitesse Co-mobile

Pour comprendre la vitesse co-mobile, imaginez une autoroute très encombrée avec deux types de voitures : des camions (l'eau) et des motos (l'huile).

• Dans l'ancienne théorie, on disait : "Les camions roulent à 80 km/h et les motos à 120 km/h".
• Mais en réalité, les motos doivent souvent ralentir pour éviter les camions, et les camions accélèrent quand les motos passent. Ils s'influencent mutuellement.

La vitesse co-mobile est une sorte de "vitesse de compensation". C'est une vitesse cachée qui explique comment les deux fluides se poussent et se tirent l'un l'autre. C'est comme si les motos et les camions avaient un accord secret pour se déplacer ensemble, créant une dynamique globale que l'on ne voyait pas en regardant chaque véhicule séparément.

C'est une variable "émergente" : elle n'existe pas au niveau d'une seule goutte, mais elle apparaît quand on regarde l'ensemble du trafic.

🧊 Le Phénomène de "Verre" : Quand le trafic se fige

L'article montre aussi que selon la vitesse à laquelle on pousse les fluides, le système change de comportement, un peu comme l'eau qui gèle ou fond.

• À très faible vitesse : Les fluides sont bloqués. Les gouttes d'huile sont coincées dans les trous de l'éponge par la tension de surface. C'est comme un embouteillage total où personne ne bouge. Les physiciens appellent cela un état "vitreux" (comme du verre).
• À vitesse moyenne : Le trafic commence à bouger, mais de manière chaotique et imprévisible. C'est là que la nouvelle théorie brille : elle décrit parfaitement ce chaos.
• À très haute vitesse : Tout se met en ordre, et les fluides coulent de manière fluide et prévisible, comme un liquide normal.

🌟 Pourquoi est-ce important ?

Avant, les ingénieurs devaient deviner des paramètres pour faire fonctionner leurs modèles de forage pétrolier ou de filtration d'eau. C'était comme conduire avec les yeux bandés en espérant ne pas sortir de la route.

Grâce à cette approche :

1. On a une carte précise reliant ce qui se passe dans les minuscules pores de l'éponge à ce qui se passe à grande échelle.
2. On découvre que la physique des éponges est plus belle et plus ordonnée qu'on ne le pensait : elle suit des lois mathématiques élégantes, similaires à celles de la chaleur, mais appliquées au mouvement.
3. On comprend enfin ce mystérieux "mouvement ensemble" (vitesse co-mobile) qui permet de prédire exactement comment les fluides se comporteront.

En résumé

Cet article nous dit : "Arrêtons de voir l'écoulement des fluides dans les éponges comme un chaos incompréhensible. En utilisant les outils de l'information et de la statistique, nous pouvons voir que c'est en fait un système ordonné, régi par des lois de 'thermodynamique du mouvement', avec ses propres règles de température (agitation) et de trafic (vitesse co-mobile)."

C'est une nouvelle façon de voir le monde, où le désordre apparent cache en réalité une structure mathématique profonde.

Résumé technique

Voici un résumé technique détaillé de l'article « Immiscible two-phase flow in porous media: a statistical mechanics approach » (Écoulement diphasique immiscible en milieu poreux : une approche de mécanique statistique) par Alex Hansen et Santanu Sinha.

1. Le Problème : L'Échelle de Mise à l'Échelle (Scale-up)

Le problème central abordé dans cet article est la description physique de l'écoulement diphasique immiscible (deux fluides non miscibles, par exemple l'eau et l'huile) dans les milieux poreux à une échelle macroscopique (échelle de Darcy), où le milieu peut être considéré comme un continu.

• Limites de l'approche actuelle : La théorie dominante, la théorie de la perméabilité relative (développée depuis les années 1930-1940), repose sur des équations phénoménologiques. Elle suppose que la pression capillaire ($p_c$) et les perméabilités relatives ($k_{rw}, k_{rn}$) ne dépendent que de la saturation ($S_w$). Bien que pratique, cette approche présente des faiblesses théoriques majeures : elle ne relie pas explicitement la physique à l'échelle des pores (microscopique) à l'échelle de Darcy (macroscopique) et échoue à capturer des comportements complexes comme l'hystérésis et les transitions de phase dynamiques.
• Objectif : Développer une théorie rigoureuse à l'échelle de Darcy qui émerge directement de la physique à l'échelle des pores, sans recourir à des hypothèses phénoménologiques ad hoc, tout en restant gérable en termes de complexité.

2. Méthodologie : Mécanique Statistique Non-Thermique

Les auteurs adoptent une approche fondée sur la généralisation de la mécanique statistique proposée par Edwin T. Jaynes dans les années 1950. Au lieu de la mécanique statistique thermique classique (basée sur l'énergie et la température), ils appliquent le principe du maximum d'entropie informationnelle de Shannon à un système hors équilibre thermodynamique.

• Principe de Jaynes : L'entropie est définie comme une mesure quantitative de l'ignorance sur un système. En maximisant l'entropie informationnelle sous contraintes (valeurs moyennes connues des variables macroscopiques), on détermine la distribution de probabilité la plus probable des configurations microscopiques.
• Application au milieu poreux :
  • Le système est modélisé comme un écoulement stationnaire dans un échantillon poreux.
  • On définit une Zone Élémentaire Représentative (REA) traversée par les fluides.
  • Les variables microscopiques incluent les configurations des fluides (champs de vitesse, distribution des phases).
  • Les contraintes macroscopiques sont les débits volumiques ($Q_w, Q_n$) et les aires occupées par les fluides ($A_w, A_n$) au sein de la REA.
• Construction du Formalisme : En maximisant l'entropie configurationnelle sous ces contraintes, les auteurs dérivent une distribution de probabilité de type Boltzmann, introduisant des multiplicateurs de Lagrange qui jouent le rôle de variables intensives nouvelles.

3. Contributions Clés et Résultats

A. Une Thermodynamique « Non-Thermique » à l'Échelle de Darcy

Le formalisme dérivé ressemble structurellement à la thermodynamique classique, mais avec des variables émergentes spécifiques à l'écoulement :

• L'Agiture ($\theta$) : Analogue à la température, elle est proportionnelle au gradient de pression ($|\nabla p|$). Elle mesure l'intensité de l'agitation des interfaces fluides.
• Le Dérivé d'Écoulement ($\mu$) : Analogue au potentiel chimique, conjugué à l'aire mouillante ($A_w$) ou à la saturation.
• La Pression d'Écoulement ($\pi$) : Une nouvelle variable intensive sans équivalent direct en thermodynamique classique, conjuguée à l'aire des pores ($A_p$).
• Relation Fondamentale : Une équation d'état reliant le débit total, l'entropie et les variables intensives, permettant de décrire l'écoulement comme un fluide non newtonien dont les propriétés rhéologiques émergent de la statistique des interfaces.

B. La Vitesse Co-mobile ($v_m$)

Une découverte majeure de cette approche est l'introduction de la vitesse co-mobile ($v_m$).

• Définition : Les vitesses thermodynamiques ($\hat{v}_w, \hat{v}_n$) dérivées du formalisme ne sont pas les vitesses moyennes mesurées expérimentalement ($v_w, v_n$). La relation entre elles fait intervenir $v_m$ :
  $$v_w = \hat{v}_w - S_n v_m$$
  $$v_n = \hat{v}_n + S_w v_m$$
• Signification : $v_m$ représente le mouvement relatif des phases dû au réarrangement des interfaces et des clusters de fluides. Elle est essentielle pour fermer le système d'équations à l'échelle de Darcy.
• Loi Empirique : L'analyse des données montre que $v_m$ dépend linéairement du dérivé d'écoulement $\mu$ (qui est lié à la dérivée de la vitesse moyenne par rapport à la saturation) : $v_m = a + b\mu$.

C. Diagramme de Phase et Transition Vitreuse

L'application de ce formalisme permet de construire un diagramme de phase pour l'écoulement diphasique stationnaire en fonction du nombre capillaire ($Ca$) et de la saturation.

• Trois Régimes :
  1. Phase Ia (Gelée/Arrêt cinétique) : Les interfaces sont bloquées par les forces capillaires (débit nul ou très faible).
  2. Phase Ib (État Vitreux) : Une transition de phase critique apparaît. Les interfaces bougent de manière intermittente, créant un état « vitreux » caractérisé par une forte hystérésis et des fluctuations importantes. Le paramètre d'ordre d'Edwards-Anderson devient non nul.
  3. Phase II/III (Non-vitreux) : À haut nombre capillaire, le système revient à une relation linéaire entre vitesse et gradient de pression (comportement newtonien effectif).
• Résultat : La transition entre le régime de loi de puissance (non-linéaire) et le régime linéaire correspond à une transition de phase vitreuse, expliquant les comportements complexes observés expérimentalement qui échappaient à la théorie de la perméabilité relative.

D. Analogie avec la Thermodynamique Ordinaire

Les auteurs établissent un lien surprenant avec la thermodynamique des mélanges binaires. Ils montrent que la vitesse co-mobile dans les milieux poreux est l'analogue exact du volume molaire co ($\nu_m$) dans les mélanges de fluides (comme l'eau et l'éthanol), où le volume total n'est pas la somme des volumes partiels. Cette analogie valide la robustesse du concept de « co-variable » au-delà des systèmes hors équilibre.

4. Signification et Perspectives

• Révolution Théorique : Cette approche démontre que l'écoulement diphasique en milieu poreux n'est pas un problème « désordonné » réservé aux ingénieurs, mais un système physique riche possédant une structure mathématique thermodynamique bien définie, même hors équilibre.
• Résolution de l'Hystérésis : La présence d'un état vitreux (Phase Ib) fournit une explication physique fondamentale à l'hystérésis observée dans les courbes de perméabilité relative et de pression capillaire, attribuant ce phénomène à la dynamique des interfaces plutôt qu'à des variables manquantes.
• Suppression de la Pression Capillaire : Le formalisme suggère que la notion de pression capillaire ($p_c = p_n - p_w$) n'est pas nécessaire à l'échelle de Darcy. Les forces motrices peuvent être décrites uniquement par les gradients d'agiture et de dérivé d'écoulement, simplifiant potentiellement la modélisation.
• Applications Futures : Ce cadre ouvre la voie à l'utilisation de la thermodynamique hors équilibre pour modéliser des écoulements transitoires et à prédire le comportement des réservoirs avec une précision accrue, en intégrant la physique des pores directement dans les équations macroscopiques.

En résumé, cet article propose un cadre unificateur puissant qui remplace les lois phénoménologiques par une dérivation statistique rigoureuse, révélant l'existence de nouvelles variables d'état (agiture, vitesse co-mobile) et de transitions de phase vitreuses dans les écoulements multiphasiques.