Mixing and spreading of gravity currents in heterogeneous porous media

Cette étude utilise des simulations numériques haute fidélité pour démontrer que, bien que l'hétérogénéité de perméabilité accélère généralement la migration des courants de gravité et améliore la dissolution par des effets dispersifs, l'interaction entre la taille de l'instabilité, la longueur de corrélation et la stratification de densité détermine in fine l'efficacité de la dissolution, les milieux homogènes surpassant souvent les milieux hétérogènes, sauf à faibles nombres de Rayleigh.

L'essentiel

Imaginez que vous versez un sirop épais et lourd dans un verre d'eau. Dans un verre parfait et uniforme, le sirop glisserait le long du bord en une feuille lisse et prévisible. Mais que se passerait-il si le verre n'était pas vide ? Et s'il était rempli d'une éponge comportant des trous de la taille d'une épingle et d'autres de la taille d'une bille ?

Tel est le monde des courants de gravité dans les milieux poreux (comme la roche souterraine ou le sol) qu'Albert Jiménez-Ramos et Juan J. Hidalgo ont exploré dans leur article. Ils ont utilisé de puissantes simulations informatiques pour observer comment les fluides se mélangent et se déplacent à travers ces couches souterraines « spongieuses », en examinant spécifiquement comment l'hétérogénéité de la roche change la donne.

Voici l'histoire de leurs découvertes, décomposée en concepts du quotidien :

La Mise en Place : Deux Types de Fluides

Les chercheurs ont étudié deux scénarios, comme deux types de fêtes différents :

1. La Fête Stable : Imaginez verser de l'huile dans l'eau. L'huile est plus légère et reste au-dessus, ou si c'est de l'eau salée, elle coule doucement. Les fluides ne se battent pas ; ils se mélangent simplement lentement à la frontière.
2. La Fête Instable : C'est comme verser un fluide qui est lourd au centre mais léger sur les bords. C'est une situation chaotique où les parties lourdes veulent couler et les parties légères veulent flotter, créant des « doigts » de fluide qui plongent ou montent, mélangeant tout violemment.

L'Effet « Éponge » : Hétérogénéité

Dans le monde réel, la roche souterraine n'est pas uniforme. C'est un mélange hétérogène de zones dures et serrées (faible perméabilité) et de zones lâches et ouvertes (forte perméabilité). Les chercheurs ont traité cela comme une éponge avec des trous aléatoires.

Ce qu'ils ont découvert :

• L'Effet Barrière : Lorsque le fluide rencontre un endroit serré et dur dans la roche, il reste coincé. C'est comme essayer de courir à travers une foule ; s'il y a un mur (une zone de faible perméabilité), vous devez faire le tour. Cette « barrière » ralentit généralement le processus de mélange car le fluide ne peut pas passer facilement.
• Le Piège : Parfois, le fluide léger reste piégé dans une poche serrée entourée de fluide lourd. C'est comme une bulle prise dans un filet. Finalement, cette bulle piégée se dissout rapidement, créant une petite explosion de mélange.

La Grande Surprise : Chaos contre Ordre

La découverte la plus intéressante a été la façon dont l'« éponge » (la roche) interagissait avec les « doigts » (le mélange instable).

• Dans le Cas Stable : La roche irrégulière agissait comme un dispersant. Elle étalait le fluide, élargissant la zone de mélange et la ralentissant. C'était comme courir à travers une forêt ; vous vous dispersez et vous n'allez pas très loin très vite.
• Dans le Cas Instable : Vous pourriez penser que la roche ralentirait les doigts chaotiques, mais ce n'était pas le cas. Le « doigtage » chaotique était si fort qu'il dominait la tendance de la roche à disperser les choses. Les doigts perçaient les barrières de la roche.
  • Le Résultat : Le mélange devenait plus efficace dans le cas instable que dans le cas stable. Les doigts rendaient l'interface entre les fluides plus étroite et plus nette, leur permettant de se dissoudre l'un dans l'autre plus rapidement que si la roche était parfaitement lisse.

Le Compromis « Vitesse contre Mélange »

L'article met en lumière une lutte entre la vitesse de déplacement du fluide et la qualité de son mélange :

• Haute Vitesse (Nombre de Rayleigh élevé) : Lorsque le fluide est très dense et se déplace rapidement, il a tendance à rester dans un courant étroit. Dans une roche uniforme, il se mélange bien. Mais dans une roche accidentée, l'« effet barrière » l'emporte. Le fluide est bloqué, se déplace plus vite le long des chemins faciles, mais se mélange moins globalement.
• Basse Vitesse (Nombre de Rayleigh faible) : Lorsque le fluide se déplace lentement, la diffusion (la tendance naturelle à s'étaler) fait le travail. Ici, la roche irrégulière aide en fait. Le chaos précoce causé par les bosses de la roche fait que le fluide se mélange mieux que dans une roche lisse et uniforme.

Le Facteur « Anisotropie » : La Direction Compte

Les chercheurs ont également examiné la direction des trous de la roche.

• Couches Horizontales (Comme un Gâteau à Étages) : Si la roche possède des couches horizontales de zones dures et molles, elle agit comme une série d'étagères. Les doigts qui coulent frappent une étagère et s'arrêtent. Cela arrête le mélange rapidement.
• Couches Verticales (Comme une Pile de Papiers) : Si les couches sont verticales, les doigts peuvent glisser le long d'elles facilement, mais tout le courant se déplace plus lentement car il doit naviguer à travers les murs verticaux.

La Conclusion

L'article conclut que l'efficacité du mélange dépend d'un équilibre délicat :

1. La Taille du Chaos : Quelle est la taille des « doigts » de mélange ?
2. La Taille des Bosses de la Roche : Quelle est la taille des trous dans l'éponge ?

Si les doigts sont petits et les bosses de la roche grandes (vitesse élevée, variance élevée), la roche agit comme une barrière, ralentissant le mélange et permettant au fluide de voyager plus loin.
Si les doigts sont grands et les bosses de la roche petites (vitesse faible), les bosses de la roche aident en fait à lancer le mélange, le rendant plus efficace qu'une roche lisse ne le serait.

En bref : Le sous-sol « spongieux » de la nature ne fait pas que ralentir les choses ; il change les règles du jeu. Parfois, il bloque l'écoulement, et parfois, si le fluide est suffisamment chaotique, il aide les fluides à se mélanger plus rapidement qu'ils ne le feraient dans un monde parfait et lisse.

Résumé technique

1. Énoncé du problème

L'étude examine l'interaction complexe entre les courants de gravité (fluides entraînés par des différences de densité) et l'hétérogénéité des milieux poreux. Alors que des recherches antérieures ont examiné les courants de gravité dans des milieux homogènes ou les effets de l'hétérogénéité sur un mélange simple, ce travail cible spécifiquement l'impact combiné de :

• Hétérogénéité de perméabilité : Représentée par des champs stochastiques log-normaux avec des variances et des longueurs de corrélation variables.
• Stratification de densité : Comparaison de scénarios stables (relation linéaire densité-concentration, par exemple intrusion d'eau salée) et de scénarios instables (relation densité non monotone, par exemple injection de CO₂ dans la saumure où se produisent des instabilités de doigts).
• Dissolution et mélange : Compréhension de la manière dont l'hétérogénéité modifie le flux de dissolution, la vitesse de migration et l'étalement du courant, ce qui est crucial pour des applications telles que la réhabilitation des aquifères et le stockage géologique du carbone.

2. Méthodologie

Les auteurs ont employé des simulations numériques haute fidélité pour résoudre les équations régissant l'écoulement miscible bidimensionnel dans les milieux poreux.

• Équations régissant : Le système est modélisé à l'aide de l'équation de Darcy pour l'écoulement et de l'équation d'advection-diffusion pour le transport de concentration. Les équations sont adimensionnalisées à l'aide du nombre de Rayleigh ($Ra$).
• Modélisation de la perméabilité : L'hétérogénéité est modélisée comme un processus stochastique log-normal. Les paramètres clés variés incluent :
  • Variance ($\sigma^2_{\log k}$) : 1, 3 et 5.
  • Longueurs de corrélation ($\lambda_x, \lambda_z$) : Utilisées pour créer des champs isotropes et anisotropes (stratification horizontale vs verticale).
  • Rapport d'anisotropie ($\gamma = \lambda_x / \lambda_z$) : Valeurs testées de 0,5, 1, 2 et 10.
• Lois de densité :
  • Cas stable : Relation linéaire ($\rho = -c$).
  • Cas instable : Relation non monotone mimant les systèmes CO₂-saumure, présentant un maximum de densité intermédiaire qui déclenche des instabilités de doigts.
• Configuration numérique :
  • Résolveur : SECUReFOAM (package OpenFOAM) utilisant la méthode des volumes finis.
  • Maillage : Grille structurée haute résolution ($16\,000 \times 800$ cellules) pour résoudre les doigts minces et les gradients de concentration.
  • Métriques : L'étude a quantifié la masse flottante, le taux de dissipation scalaire ($\chi$), le flux de dissolution, la longueur de l'interface et la largeur de l'interface.

3. Contributions clés

L'article fournit une analyse complète de la manière dont l'hétérogénéité modifie la physique des courants de gravité, remettant spécifiquement en question l'hypothèse selon laquelle l'hétérogénéité améliore toujours le mélange. Les contributions clés incluent :

• Mécanismes antagonistes : Démonstration que dans les régimes instables (formation de doigts), l'effet dispersif de l'hétérogénéité est contrebalancé par l'instabilité de doigts, conduisant à des résultats différents de ceux des régimes stables.
• Effets de l'anisotropie : Identification d'un « effet barrière » dans les milieux stratifiés horizontalement (fort $\gamma$) qui supprime la convection verticale et réduit la dissolution, contrastant avec le comportement dans les milieux stratifiés verticalement.
• Seuils de variance : Établissement que la variance du champ de perméabilité n'améliore significativement l'efficacité de dissolution par rapport aux cas homogènes que lorsque le nombre de Rayleigh est faible.
• Piégeage de « blobs » : Description détaillée du mécanisme par lequel le fluide flottant reste piégé dans des zones de faible perméabilité, formant des « blobs » qui se dissolvent rapidement et créent une formation de doigts ascendante localisée, même dans une stratification stable.

4. Résultats clés

A. Milieux homogènes vs hétérogènes

• Cas stable : L'hétérogénéité agit principalement comme un agent dispersif. Elle élargit l'interface de mélange, réduit le gradient de concentration et par conséquent diminue la dissolution. Cependant, elle augmente la vitesse de migration de la pointe du courant.
• Cas instable : L'hétérogénéité accélère le début de la convection. Bien que cela augmente initialement la dissolution, l'effet dispersif de l'hétérogénéité finit par rétrécir l'interface par rapport au cas stable, mais crée une interaction complexe où le cas homogène surpasse souvent le cas hétérogène en masse totale dissoute pour des $Ra$ élevés.

B. Impact de l'anisotropie ($\gamma$)

• Stratification horizontale ($\gamma > 1$) : Crée un effet barrière. Les couches horizontales de faible perméabilité piègent le fluide flottant, supprimant la formation de doigts et arrêtant localement la convection. Cela conduit à une réduction de la dissolution totale par rapport aux milieux homogènes, malgré un pic initial du flux de dissolution.
• Stratification verticale ($\gamma < 1$) : Agit comme une barrière à la migration horizontale mais favorise la descente des doigts verticaux. Cela résulte en un comportement de dissolution intermittent lorsque le courant traverse les barrières de perméabilité, avec un début plus précoce du régime d'« arrêt » (où l'accumulation de fluide dense au fond arrête le mélange).

C. Impact de la variance ($\sigma^2_{\log k}$)

• Variance élevée : Augmente l'effet dispersif, élargit l'interface et réduit le gradient de concentration entraînant la dissolution.
• Compromis d'efficacité : Pour des nombres de Rayleigh élevés (convection forte), les cas hétérogènes dissolvent généralement moins de masse que les cas homogènes car les effets barrière et dispersif dominent.
• Nombres de Rayleigh faibles : L'hétérogénéité peut améliorer la dissolution par rapport au cas homogène en raison du début précoce de la convection et de l'allongement de l'interface, car la diffusion joue un rôle plus dominant par rapport à la taille de l'instabilité.

D. Dynamique de l'interface

• Dans les cas instables, l'interface est généralement plus étroite que dans les cas stables en raison de la dominance des doigts sur la diffusion.
• La présence de « blobs » piégés dans les régions de faible perméabilité provoque un comportement oscillatoire de la longueur de l'interface et du flux de dissolution, en particulier dans les milieux anisotropes.

5. Importance et implications

Cette recherche est cruciale pour optimiser le stockage géologique du carbone (GCS) et la réhabilitation des aquifères :

• Stockage du carbone : Les résultats suggèrent que dans les scénarios de haute perméabilité et de forte convection (fort $Ra$), l'hétérogénéité naturelle pourrait en réalité entraver le piégeage par solubilité du CO₂ par rapport aux modèles homogènes idéalisés. Les ingénieurs doivent tenir compte de l'« effet barrière » de la stratification horizontale qui peut limiter le mélange vertical.
• Réhabilitation : Pour la réhabilitation des DNAPL, comprendre que l'hétérogénéité peut piéger des fluides dans des « blobs » qui se dissolvent rapidement aide à prédire la longévité des efforts de réhabilitation.
• Lignes directrices de modélisation : L'étude souligne qu'une simple mise à l'échelle des milieux hétérogènes est insuffisante. L'interaction entre la taille des instabilités convectives et la longueur de corrélation du champ de perméabilité est le facteur déterminant pour l'efficacité de la dissolution. Une dissolution optimale nécessite un équilibre spécifique entre le nombre de Rayleigh et l'intensité de l'hétérogénéité.

En conclusion, l'article établit que l'hétérogénéité n'améliore pas universellement le mélange ; son impact est conditionné par la stabilité de la stratification de densité et les échelles relatives des instabilités d'écoulement par rapport à l'hétérogénéité géologique.