📄 chemistry

Effective diffusion in sharp-front reactive transport experiments

Cette étude présente une solution analytique de forme fermée pour le transport réactif à front net en milieux poreux, démontrant que le coefficient de diffusion effectif déduit expérimentalement diffère fondamentalement du coefficient de diffusion moléculaire intrinsèque en raison d'effets couplés d'advection-diffusion-thermodynamique, un cadre validé par des expériences de carbonatation de ciment Portland sous diverses conditions d'écoulement.

Auteurs originaux : Viktoriya M. Yarushina, Reinier van Noort, Yury Y. Podladchikov, Yannik F. Schneider

Publié 2026-09-01

📖 6 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Viktoriya M. Yarushina, Reinier van Noort, Yury Y. Podladchikov, Yannik F. Schneider

Article original sous licence CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Profondément sous terre, là où la roche et le fluide se rencontrent, une bataille chimique silencieuse est constamment engagée. Lorsque de l'eau ou du gaz riche en dioxyde de carbone circule à travers des matériaux poreux comme le ciment ou la roche, cela déclenche une réaction qui transforme la structure solide. Ce processus ne se produit pas uniformément dans l'ensemble du matériau ; au contraire, il crée une frontière distincte, un front net, qui progresse vers l'intérieur, séparant la zone altérée et réagie du matériau d'origine intact. Ce phénomène est crucial pour la sécurité du stockage géologique du carbone, où des scellés de ciment sont utilisés pour piéger les gaz à effet de serre profondément sous terre. Si ces scellés se dégradent trop rapidement, le gaz stocké pourrait s'échapper. Pendant des décennies, les scientifiques ont tenté de mesurer la vitesse à laquelle ce front se déplace pour prédire la durée de vie de ces scellés, en supposant souvent que la vitesse est dictée par la facilité avec laquelle les molécules dérivent à travers les minuscules pores du matériau, un processus connu sous le nom de diffusion.

Cependant, une nouvelle étude menée par des chercheurs de l'Institut de technologie de l'énergie et de l'Université de Lausanne remet en question la façon dont nous interprétons ces mesures. Ils ont découvert que la vitesse à laquelle ce front de réaction voyage n'est pas simplement le reflet de la vitesse à laquelle les molécules se déplacent dans le fluide. Au lieu de cela, le mouvement est régi par une interaction complexe entre l'écoulement du fluide, la diffusion des molécules et un effet de « tampon » chimique où le matériau solide absorbe la substance réagissante. Les chercheurs ont développé un nouveau cadre mathématique pour décrire ce mouvement, montrant que la vitesse « effective » mesurée lors des expériences est en réalité un résultat composé du mouvement du fluide et de la quantité massive de stockage chimique que le matériau solide peut contenir. Cette distinction est vitale car elle explique pourquoi différentes expériences, même utilisant des méthodes très différentes, produisent souvent des chiffres déroutants et divergents sur la vitesse de dégradation du ciment.

Pour tester cette théorie, l'équipe s'est tournée vers le ciment Portland, le matériau utilisé pour sceller les puits destinés au stockage du carbone. Ils ont soumis des échantillons de ciment identiques à deux environnements très différents. Dans un premier ensemble d'expériences, le ciment reposait calmement dans un récipient rempli d'eau saturée en dioxyde de carbone ou de CO2 supercritique humide. Il s'agit d'une expérience en mode « batch », où il n'y a pas de flux forcé, et le fluide repose simplement autour de l'échantillon. Dans le second ensemble, ils ont forcé ces mêmes fluides à circuler à travers les échantillons de ciment sous pression, simulant un scénario où le fluide circule activement à travers un puits. La sagesse conventionnelle suggérerait que le flux forcé pousserait le front de réaction beaucoup plus vite que les conditions calmes et statiques. Pourtant, lorsque les chercheurs ont mesuré la profondeur de la carbonatation après des semaines et des mois, les résultats étaient étonnamment similaires. Dans les deux scénarios, calme et à flux forcé, le front de réaction progressait à un rythme qui suivait le même schéma : la distance parcourue augmentait avec la racine carrée du temps. Ce schéma est la marque distinctive d'un processus contrôlé par la diffusion, et non par la vitesse du fluide en mouvement.

Les chercheurs ont utilisé cette similitude surprenante pour expliquer un puzzle de longue date dans le domaine. Lorsque les scientifiques mesurent la vitesse à laquelle un front de réaction se déplace lors d'une expérience, ils calculent souvent un « coefficient de diffusion effectif » et le traitent comme une propriété fixe du matériau, telle que la densité d'une roche. La nouvelle étude montre que c'est une erreur d'interprétation. La vitesse du front n'est pas déterminée uniquement par la vitesse à laquelle les molécules de dioxyde de carbone peuvent se faufiler à travers l'eau dans les pores. Elle est également déterminée par la capacité du ciment solide à absorber et stocker ce dioxyde de carbone. Le ciment agit comme une éponge, absorbant le dioxyde de carbone et le convertissant en un nouveau minéral. Parce que le solide contient beaucoup plus de carbone que le fluide, le front avance lentement, non pas parce que le fluide est léthargique, mais parce que la réaction doit « remplir » l'éponge solide avant de pouvoir progresser. Les chercheurs ont introduit une nouvelle façon de voir cela, en utilisant des nombres adimensionnels pour décrire l'équilibre entre le flux du fluide et la capacité chimique du solide. Ils ont découvert que même lorsqu'un fluide est forcé à travers le ciment, la capacité de stockage chimique est si grande que le flux ne fait presque aucune différence sur la vitesse du front.

Pour confirmer que le front de réaction était effectivement net et distinct, l'équipe a examiné les échantillons de ciment sous des microscopes puissants et a utilisé la lumière infrarouge pour cartographier les changements chimiques. Ils ont observé la limite entre le ciment réagi et non réagi avec une précision extrême. Les images ont révélé que la transition était incroyablement nette, se produisant sur une distance d'environ 50 à 100 micromètres, soit environ la largeur d'un cheveu humain. Dans cette zone étroite, la composition chimique changeait brusquement. L'hydrate de silicate de calcium, la colle qui maintient le ciment ensemble, disparaissait et était remplacé par du carbonate de calcium. Les chercheurs ont également observé que le calcium du ciment ne se dissolvait pas et ne s'échappait pas ; au contraire, il restait sur place, se réorganisant en une nouvelle structure minérale. Cette observation soutient l'idée que le front est une interface nette où la phase solide agit comme un tampon stationnaire, absorbant le fluide entrant et réagissant localement.

Les implications de ces découvertes sont significatives pour la conception et le test des matériaux destinés au stockage du carbone. L'étude démontre que le simple fait d'augmenter le flux de fluide à travers un puits ne va pas nécessairement accélérer la dégradation du scellé de ciment, car le processus est dominé par la capacité de stockage chimique du matériau plutôt que par la vitesse du fluide. Cela signifie que les expériences menées dans des conteneurs statiques et calmes peuvent être directement comparées à celles avec un flux forcé, à condition que les conditions chimiques soient similaires. La vitesse « effective » mesurée en laboratoire n'est pas une propriété fondamentale de la capacité du fluide à se déplacer, mais le résultat d'une danse chimique spécifique entre le fluide et le solide. En comprenant cela, les scientifiques pourront mieux prédire la durée de vie d'un scellé de ciment dans l'environnement hostile d'un puits profond, garantissant que le dioxyde de carbone reste en sécurité sous terre pendant des milliers d'années. Le travail clarifie que le mouvement de ces fronts chimiques est une histoire de stockage et d'équilibre, et non seulement de flux.

Résumé technique : Diffusion effective dans les expériences de transport réactif à front abrupt

Énoncé du problème
Le transport réactif dans les milieux poreux, tels que les systèmes de ciment exposés au CO₂, génère fréquemment des fronts d'altération abrupts séparant les régions réagies des régions non altérées. Bien que ces fronts soient essentiels pour comprendre les interactions fluide-roche dans le stockage géologique du CO₂ et l'intégrité des puits, leurs taux de propagation sont souvent interprétés à l'aide de lois de diffusion empiriques (par exemple, $x \propto \sqrt{t}$ ). Une ambiguïté significative subsiste dans la littérature concernant le « coefficient de diffusion effectif » ( $D_{eff}$ ) déduit de ces expériences. Les valeurs rapportées varient considérablement et sont souvent de plusieurs ordres de grandeur inférieures aux coefficients de diffusion moléculaire intrinsèques, ce qui crée une confusion sur la question de savoir si ces valeurs représentent de véritables propriétés matérielles ou des artefacts de processus de transport-réaction couplés. De plus, il reste difficile de comparer de manière cohérente des expériences menées sous différentes conditions aux limites (par exemple, en mode batch ou en flux forcé) ou de comprendre comment l'advection influence la propagation du front lorsque le tamponnage chimique est fort.

Méthodologie
Les auteurs développent une solution analytique en forme fermée pour le transport réactif advection-diffusion présentant un front de réaction mobile. Le modèle suppose :

Conservation de la masse : Le transport d'une espèce réactive est régi par l'advection et la diffusion dans la phase fluide, couplés à un échange avec la phase solide via un équilibre de partition.
Interface abrupte : Une frontière mobile (de type Stefan) sépare les zones réagies et non réagies, à travers laquelle la porosité, la densité et les concentrations présentent des sauts finis.
Condition de Rankine-Hugoniot : La vitesse du front est dérivée de la conservation de la masse à l'interface, tenant compte des discontinuités de la concentration totale.

La solution introduit deux nombres adimensionnels clés :

Nombre de Péclet ( $Pe$ ) : Contrôle l'importance relative de l'advection par rapport à la diffusion.
Nombre de Stefan Chimique ($Ste$) : Contrôle le taux de propagation du front en quantifiant le rapport entre le flux de concentration mobile moteur et l'accumulation totale de l'espèce requise à travers le front (intégrant le stockage de la phase solide et les sauts de concentration).

Le cadre est validé et appliqué à des données expérimentales provenant d'échantillons de ciment Portland exposés à des fluides riches en CO₂, tant en mode batch qu'en flux forcé. Les expériences ont utilisé du CO₂ supercritique humide et de l'eau saturée en CO₂. La caractérisation a été réalisée par microscopie électronique à balayage avec électrons rétrodiffusés (SEM-BSE), spectroscopie de dispersion d'énergie (EDS) et microscopie infrarouge à transformée de Fourier (FTIR) afin de résoudre la morphologie des fronts et les gradients chimiques.

Contributions clés

Cadre analytique : L'article fournit une solution analytique de type Stefan qui prend explicitement en compte les discontinuités de porosité, de densité et de concentration à une interface mobile. Cette solution unifie les régimes dominés par la diffusion et l'advection grâce à un paramètre composite $\Lambda = Pe / (Ste \cdot R)$ , où $R$ est le facteur de retard.
Redéfinition de la diffusion effective : L'étude démontre que le « coefficient de diffusion effectif » déduit de la propagation du front n'est pas une propriété intrinsèque du fluide poreux. Il s'agit d'un coefficient de propagation composé dépendant du coefficient de diffusion moléculaire ( $D_f$ ) mis à l'échelle par le nombre de Stefan chimique ( $D_{eff} = D_f \cdot Ste$ ).
Résolution des divergences : Le cadre explique pourquoi les coefficients de diffusion déduits expérimentalement sont souvent de plusieurs ordres de grandeur inférieurs aux valeurs moléculaires : un faible nombre de Stefan implique qu'une petite fraction du flux diffusif contribue à l'avancement du front, la majorité étant consommée par le stockage de la phase solide.
Validation expérimentale : Les auteurs montrent que, malgré des méthodologies d'exposition fondamentalement différentes (batch vs flux forcé), les taux de propagation du front peuvent être similaires si les paramètres adimensionnels directeurs sont comparables.

Résultats

Propagation contrôlée par la diffusion : Dans les deux types d'expériences (batch et flux forcé) impliquant du ciment Portland, la propagation du front de carbonatation a suivi une loi de racine carrée du temps ( $x \propto \sqrt{t}$ ). Les nombres de Péclet calculés étaient faibles ( $Pe \approx 0,2$ pour le CO₂ supercritique humide et encore plus bas pour l'eau saturée en CO₂), indiquant que même sous un flux imposé, le transport restait dominé par la diffusion.
Coefficients effectifs vs intrinsèques : Les coefficients de diffusion effectifs déduits étaient de $1,4 \times 10^{-11}$ m²/s (eau saturée en CO₂) et $2,5 \times 10^{-11}$ m²/s (CO₂ supercritique humide). Ils sont de plusieurs ordres de grandeur inférieurs à la diffusion moléculaire intrinsèque du CO₂ dans l'eau ( $\sim 10^{-8}$ m²/s). Les auteurs attribuent cela à de faibles nombres de Stefan ( $Ste \sim 10^{-3}$ ), reflétant la haute capacité de carbonatation de la phase solide par rapport au contraste de concentration du fluide.
Morphologie du front : Les observations par SEM-BSE, EDS et FTIR à haute résolution ont confirmé que les fronts de réaction sont extrêmement abrupts, avec des zones de transition localisées à des échelles micrométriques (environ 50–110 µm).
- L'analyse FTIR a révélé que l'apparition des carbonates et la disparition du gel C-S-H se produisent sur des zones étroites et distinctes.
- Les profils EDS ont montré des transitions nettes dans les rapports Ca/Si, soutenant l'hypothèse d'une interface abrupte et d'un équilibre local.
- Dans les expériences en flux forcé, une zone d'enrichissement étroite en Ca/Si a été observée immédiatement derrière le front, mais le front global est resté abrupt, validant l'approximation de l'immobilité du calcium dans le régime à faible $Pe$ .

Signification et affirmations
L'article affirme fournir une base générale et physiquement transparente pour interpréter la propagation des fronts d'altération dans les milieux poreux réactifs. En distinguant la diffusion moléculaire intrinsèque du coefficient effectif de propagation du front, l'étude résout la contradiction apparente entre les très variables valeurs de diffusion expérimentales et les propriétés de diffusion moléculaire constantes. Les auteurs soutiennent que leur cadre permet une comparaison cohérente d'expériences réalisées sous différentes conditions limites (batch vs flux traversant) et compositions de fluides. Les résultats soutiennent l'applicabilité des formulations d'interface abrupte et d'équilibre local pour modéliser le transport réactif dans des systèmes allant des barrières de ciment technique aux zones d'interaction fluide-roche naturelles, à condition que les nombres adimensionnels directeurs ( $Pe$ , $Ste$ et $\Lambda$ ) soient correctement caractérisés.

Articles similaires