Quantifying Injection-Driven Mass Transfer within Porous Media via Time-Elapsed X-ray micro-Computed Tomography

Cette étude évalue trois approches analytiques pour quantifier le transfert de masse par injection dans des milieux poreux à l'aide de microtomographie X, en introduisant une technique de filtrage pour réduire les biais et en démontrant que le choix de la méthode dépend d'un compromis entre la résolution physique souhaitée et les ressources de calcul disponibles.

L'essentiel

🫧 Le Grand Jeu de la Dissolution : Comment mesurer la disparition des bulles ?

Imaginez que vous avez un grand verre rempli de sable (c'est notre milieu poreux). Vous y avez piégé des bulles de gaz (comme de l'hydrogène) entre les grains de sable. Ensuite, vous commencez à verser de l'eau doucement à travers le sable. Petit à petit, l'eau "mange" les bulles de gaz : c'est ce qu'on appelle la dissolution.

Le but de cette étude est de répondre à une question simple : À quelle vitesse l'eau mange-t-elle les bulles ?

Pour le savoir, les chercheurs ont utilisé une machine incroyable : un scanner 3D aux rayons X (comme un scanner médical, mais pour des cailloux). Ils ont pris des photos de l'intérieur du sable à intervalles réguliers pour voir les bulles disparaître.

Mais il y a un problème : comment calculer la vitesse de disparition à partir de ces photos ? Il existe trois façons de faire, un peu comme trois façons différentes de compter les pièces d'un puzzle.

🛠️ Les Trois Méthodes (Les Trois Recettes)

Les chercheurs ont testé trois "recettes" mathématiques sur les mêmes données pour voir laquelle fonctionnait le mieux.

1. La Méthode "Moyenne du Quartier" (SAC)

• L'analogie : Imaginez que vous regardez une rivière. Au lieu de compter chaque goutte d'eau, vous prenez une moyenne de la vitesse de l'eau sur toute la largeur de la rivière à un endroit donné.
• Comment ça marche : Cette méthode regarde le sable par "tranches" (comme des tranches de saucisson). Elle calcule la moyenne de tout ce qui se passe dans une tranche.
• Le pour : C'est rapide et facile à calculer.
• Le contre : C'est un peu flou. On perd les détails. On ne sait pas exactement quelle bulle a disparu, juste que "globalement, il y a moins de gaz ici". C'est comme dire "il fait chaud dans la pièce" sans savoir si c'est près du radiateur ou près de la fenêtre.

2. La Méthode "Compteur de Bulles" (NPC)

• L'analogie : Ici, on compte chaque bulle individuellement, sans se soucier de son histoire. Si une bulle rétrécit, on note. Si elle grossit, on note aussi. On fait une énorme liste de toutes les bulles et on fait une moyenne statistique.
• Comment ça marche : On suit chaque bulle. On applique une formule à toutes les bulles qui changent de taille.
• Le pour : On a beaucoup de données.
• Le contre : C'est bruyant. Parfois, les bulles ne disparaissent pas parce que l'eau les mange, mais parce qu'elles bougent et fusionnent avec d'autres (comme deux gouttes de pluie qui se rejoignent). Si on compte ces mouvements comme de la dissolution, on se trompe. C'est comme si on comptait les gens qui changent de file dans un supermarché comme s'ils avaient acheté des produits !

3. La Méthode "Détective Intelligente" (CPC)

• L'analogie : C'est la méthode la plus pointue. Le détective ne regarde pas n'importe quelle bulle. Il ne regarde que les bulles qui ont complètement disparu. Il se dit : "Si une bulle a disparu, c'est qu'elle a été mangée par l'eau, pas juste déplacée."
• Comment ça marche : On classe les bulles. On ignore celles qui grossissent ou qui bougent bizarrement. On ne garde que celles qui ont été totalement dissoutes pour calculer la vitesse.
• Le pour : C'est très précis. On peut même voir la "frontière" de l'eau qui avance (le front de dissolution) et comprendre comment l'eau se comporte à l'échelle microscopique.
• Le contre : C'est très lent et demande beaucoup de puissance de calcul. C'est comme trier des millions de pièces de puzzle une par une.

🏆 Le Verdict des Chercheurs

Après avoir comparé les trois méthodes, voici ce qu'ils ont découvert :

1. Le résultat global est le même : Que vous utilisiez la méthode rapide (moyenne), la méthode brute (compteur) ou la méthode intelligente (détective), vous obtenez une estimation de la vitesse de dissolution qui est du même ordre de grandeur. C'est-à-dire que si l'une dit "10 bulles par seconde", les autres diront "8" ou "12". C'est bien !
2. Le détail change tout :
   • La méthode rapide (SAC) est excellente pour avoir une idée générale, mais elle lisse trop les détails.
   • La méthode brute (NPC) est pleine de "bruit" à cause des bulles qui bougent sans se dissoudre.
   • La méthode intelligente (CPC) révèle des phénomènes cachés, comme le fait que l'eau avance de manière irrégulière (comme un doigt qui traverse le sable) et permet de voir exactement où se trouve la zone de dissolution maximale.

💡 Pourquoi est-ce important ?

Ces recherches ne servent pas juste à compter des bulles. Elles sont cruciales pour des applications réelles comme :

• Le stockage d'énergie : On veut injecter de l'hydrogène ou du CO2 dans le sous-sol pour le stocker. Il faut savoir comment ces gaz vont se dissoudre dans l'eau souterraine pour ne pas qu'ils remontent à la surface.
• La dépollution : Pour nettoyer les nappes phréatiques polluées.

En résumé :
Si vous voulez une réponse rapide et approximative, utilisez la méthode simple (SAC). Mais si vous voulez comprendre exactement ce qui se passe à l'intérieur de la roche, avec tous ses détails et ses pièges, vous devez utiliser la méthode intelligente (CPC), même si cela demande beaucoup plus de travail et d'ordinateurs puissants.

C'est un peu comme choisir entre regarder une photo floue d'un match de foot (vous voyez qui gagne) et regarder le match en 4K avec des ralentis (vous voyez chaque passe et chaque faute). Les deux vous disent qui a marqué, mais l'un vous donne l'histoire complète !

Résumé technique

1. Problématique et Contexte

La compréhension de la dynamique du transfert de masse interphasique dans les milieux poreux est cruciale pour de nombreuses applications, allant de la dépollution des eaux souterraines au stockage géologique de l'énergie (hydrogène, carbone). Bien que la tomographie micro-ordinateur X (µCT) permette d'observer in situ ces processus, l'analyse quantitative du transfert de masse se heurte à des limitations spatiales et temporelles inhérentes à l'imagerie.

Le défi principal réside dans le choix de la méthode d'analyse pour résoudre l'équation de transfert de masse (Eq. 3) à partir des données µCT. La littérature propose actuellement trois cadres analytiques distincts :

1. Approche à concentration moyenne par tranche (SAC - Slice-Averaged Concentration) : Utilise un modèle d'advection-diffusion 1D.
2. Approche par amas non classés (NPC - Non-Classified per-Cluster) : Traite tous les changements de volume des amas de gaz comme des événements de transfert de masse.
3. Approche par amas classés (CPC - Classified per-Cluster) : Filtre les événements pour ne considérer que ceux correspondant à une dissolution complète au front de dissolution.

L'objectif de cette étude est d'évaluer et de comparer ces trois méthodes appliquées au même jeu de données expérimentales (dissolution d'hydrogène dans un milieu poreux) pour déterminer leurs forces, faiblesses et domaines d'application respectifs.

2. Méthodologie

Données Expérimentales :
L'étude utilise des données de tomographie µCT temporelles issues d'une étude précédente (Patmonoaji et al., 2023). Le système consiste en un empilement granulaire de plastique (porosité ~45 %) saturé en eau, contenant des amas d'hydrogène (H2) piégés. Des injections d'eau non saturée ont été réalisées à quatre débits différents (0,10 ; 0,25 ; 0,50 et 1,0 mL/min). La résolution voxel est de 16,472 µm.

Traitement des Données et Filtrage :
Un défi majeur identifié est la remobilisation des amas (fusion ou fragmentation) induite par la dissolution, qui fausse les mesures de transfert de masse en créant une perte de volume apparente non liée à la dissolution pure.

• Nouvelle contribution méthodologique : Les auteurs introduisent un filtre basé sur le ratio de volume. Ils calculent le ratio entre le volume total gagné par les amas (croissance/fusion) et le volume total perdu (dissolution) par intervalle de temps.
• Seuil : Les intervalles où ce ratio dépasse 0,2 sont exclus de l'analyse pour garantir que les données restantes sont dominées par la dissolution et non par la remobilisation. Ce filtrage est appliqué uniformément aux trois approches pour assurer la comparabilité.

Approches Analytiques Comparées :

1. SAC : Résout l'équation d'advection-diffusion sur des tranches transversales pour obtenir un profil de concentration moyen.
2. NPC : Applique la loi du transfert de masse à tous les changements de volume observés, en supposant un gradient de dissolution maximal partout.
3. CPC : Classe les amas selon leur évolution (dissolution complète, partielle, croissance, etc.). Seuls les amas ayant dissous complètement sont utilisés pour calculer le coefficient de transfert, supposant qu'ils se trouvent au front de dissolution où le gradient est maximal.

3. Résultats Clés

Coefficients de Transfert de Masse ($K$) :

• Les trois méthodes produisent des estimations du coefficient de transfert de masse moyen ($K_{seq}$) qui sont du même ordre de grandeur pour un débit d'injection donné.
• Relation avec le débit : Tous les coefficients augmentent avec le débit d'injection (nombre de Reynolds).
• Comparaison des valeurs :
  • À des débits plus faibles, l'ordre est généralement : $K_{CPC} > K_{NPC} > K_{SAC}$.
  • À des débits élevés (1,0 mL/min), l'approche NPC donne des valeurs légèrement supérieures à la CPC.
  • L'approche SAC est moins coûteuse en calcul mais tend à sous-estimer légèrement les valeurs par rapport aux approches par amas.

Profil de Concentration et Résolution Spatiale :

• SAC : Produit des profils de concentration très dilués (sous-estimés) en raison de l'hypothèse d'uniformité radiale qui dilue artificiellement la concentration sur toute la section transversale. Elle ne permet pas de résoudre les détails microscopiques.
• NPC : Génère des distributions de concentration très variables, incluant des valeurs physiquement impossibles (concentration > saturation ou < 0) car elle ne filtre pas les événements de croissance ou de remobilisation. Elle perd la résolution spatiale fine.
• CPC : Offre la meilleure résolution phénoménologique. Elle permet de visualiser le front de dissolution non uniforme et de localiser spatialement les zones de concentration négative (sous-estimation du coefficient) au bord d'attaque du front. Cependant, elle est plus sensible aux erreurs de mesure sur les petits amas et nécessite un filtrage rigoureux.

Impact du Filtrage de Remobilisation :
Le filtrage basé sur le ratio de volume a permis d'éliminer les biais significatifs, en particulier dans les premières phases des séquences d'injection où la remobilisation était prédominante. Sans ce filtrage, les estimations de transfert de masse auraient été surestimées.

4. Contributions Principales

1. Comparaison systématique : Première étude à appliquer et comparer directement les trois cadres analytiques (SAC, NPC, CPC) sur le même jeu de données temporelles µCT.
2. Filtrage de la remobilisation : Développement et validation d'une technique de filtrage basée sur le ratio volume gagné/volume perdu pour isoler les événements de dissolution pure des artefacts de remobilisation d'amas.
3. Évaluation des compromis (Trade-offs) : Démonstration claire que le choix de la méthode dépend du compromis entre la résolution physique souhaitée et les ressources computationnelles disponibles.
   • SAC : Faible coût, résolution macroscopique, perte des détails microscopiques.
   • NPC : Coût moyen, robustesse statistique mais perte de la résolution spatiale et des phénomènes complexes.
   • CPC : Coût élevé, haute résolution des phénomènes (fronts non uniformes, remobilisation), mais sensibilité accrue aux erreurs de mesure et aux valeurs aberrantes.

5. Signification et Conclusion

Cette étude fournit un cadre de référence pour les chercheurs travaillant sur le transfert de masse en milieux poreux via l'imagerie 3D. Elle démontre que, bien que les méthodes donnent des résultats globalement cohérents (même ordre de grandeur) pour les paramètres globaux comme le coefficient de transfert de masse, elles divergent considérablement lorsqu'il s'agit de décrire des phénomènes complexes comme les profils de concentration locaux ou la dynamique des fronts de dissolution.

La conclusion majeure est qu'il n'existe pas de méthode "meilleure" universelle. Le choix doit être guidé par l'objectif de l'étude :

• Pour des paramètres systémiques généraux, l'approche SAC est suffisante et efficace.
• Pour comprendre la dynamique fine, la morphologie des fronts de dissolution et les mécanismes de remobilisation, l'approche CPC est indispensable, malgré son coût computationnel et sa complexité de mise en œuvre.

Les auteurs soulignent la nécessité de futures études sur différents milieux poreux et gaz pour généraliser davantage ces résultats, car la variabilité du milieu poreux reste un facteur dominant dans la dispersion des coefficients de transfert mesurés.