Caractérisation expérimentale et modélisation de la dispersion non-Fickiéenne dans les aquifères

Experimental characterization and modeling non-Fickian dispersion in aquifers

• Auteur / Autrice : Filip Gjetvaj
• Direction : Philippe Gouze
• Type : Thèse de doctorat
• Discipline(s) : Sciences de la Terre et de l'Univers
• Date : Soutenance le 12/11/2015
• Etablissement(s) : Montpellier
• Ecole(s) doctorale(s) : Systèmes Intégrés en Biologie, Agronomie, Géosciences, Hydrosciences, Environnement (Montpellier ; École Doctorale ; 2009-2015)
• Partenaire(s) de recherche : Laboratoire : Géosciences (Montpellier)
• Jury : Examinateurs / Examinatrices : Philippe Gouze, Branko Bijeljic, Michel Quintard, Anna Russian, Alberto Guadagnini, Séverin Pistre
• Jury : Rapporteur / Rapporteuse : Branko Bijeljic, Michel Quintard

Mots clés

Mots clés contrôlés

Aquifères

Analyse numérique

Matériaux poreux

Mots clés libres

Milieux poreux

Transport hydrodispersif

Expériences de laboratoire

Changement d’échelles

Multirate mass transfer

Méthodes numériques de marche aléatoire

Porous media

Hydrodynamic transport

Laboratory experiments.

Upscaling

Multirate mass transfer

Random walk numerical methods

Résumé

Ces travaux ont pour objectif de modéliser les mécanismes de dispersion dans les aquifères. L’hétérogénéité du champ de vitesse et le transfert de masse entre zones immobiles et mobiles sont deux origines possibles du comportement non-Fickéen, jusqu’alors étudiées de façon séparée. Notre hypothèse de départ est que ces deux mécanismes coexistent. Nos travaux comprennent : 1) des expériences de traçage sur colonnes de billes de verre et carottes de grès de Berea, en mode flow-through et push-pull, et 2) des simulations numériques réalisées à partir d’images en microtomographie RX segmentées en trois phases : solide, vide et microporosité. L’analyse du champ de vitesse (Stokes) montre l’importance de la discrétisation spatiale et de la prise en compte de la microporosité. Les résultats des simulations de transport (en utilisant la méthode time domain random walk) permettent de quantifier l’effet combiné de l’hétérogénéité du champ de vitesse et des transferts diffusifs dans la fraction micro-poreuse de la roche sur la dispersion non-Fickéenne, caractérisée à partir des courbes de restitution (BTC). Ces résultats sont cohérents avec les observations expérimentales. Nous concluons que ces deux effets doivent être pris en compte même si leur identification à partir de la forme des BTCs issues des traçages des milieux naturels (souvent caractérisés par de faible valeurs du nombre de Peclet ) reste difficile. Enfin, un modèle moyen macroscopique 1D est proposé dans le cadre d’une approche de type continuous time random walk dans laquelle des distributions spécifiques du temps de transfert des particules sont construites pour chacun des deux mécanismes de transport.

His work aims at modeling hydrodynamic dispersion mechanisms in aquifers. So far both flow field heterogeneity and mobile-immobile mass transfer have been studied separately for explaining the ubiquitously observed non-Fickian behaviors, but we postulate that both mechanisms contribute simultaneously. Our investigations combine laboratory experiments and pore scale numerical modeling. The experimental rig was designed to enable push-pull and flow through tracer tests on glass bead columns and Berea sandstone cores. Modeling consists in solving Stokes flow and solute transport on 3D X-ray microtomography images segmented into three phases: solid, void and microporosity. Transport is modeled using time domain random walk. Statistical analysis of the flow field emphasizes the importance of the mesh resolution and the inclusion of the microporosity. Results from the simulations show that both the flow field heterogeneity and the diffusive transport in the microporous fraction of the rock contribute to the overall non-Fickian transport behavior observed, for instance, on the breakthrough curves (BTC). These results are supported by our experiments. We conclude that, in general, this dual control must be taken into account, even if these different influences can hardly be distinguished from a qualitative appraisal of the BTC shape, specifically for the low values of the Peclet number that occurs in natural conditions. Finally, a 1D up-scaled model is developed in the framework of the continuous time random walk, where the influences of the flow field heterogeneity and mobile-immobile mass transfer are both taken into account using distinct transition time distributions.