Academic literature on the topic 'Transport en milieux poreux fracturés'

RésuméUne méthode d’éléments finis mixte hybride est appliquée pour l’approximation de l’écoulement associé au transport en milieu poreux non saturé. Le développement de ce modèle s’effectue dans le cadre du projet européen ARWET, lequel a pour objectif l’étude de nouvelles potentialités de dissipation des pesticides dans les zones humides. La formulation du modèle bidimensionnel est fondée sur les propriétés de l’espace de Raviart-Thomas. L’écueil numérique que posent les problèmes à convection dominante est surmonté par l’introduction d’un limiteur de flux alors qu’un limiteur de pente est généralement utilisé dans la littérature. Le limiteur suggéré est inconditionnellement stable et permet de conserver la précision des résultats à nombre de Peclet élevé.

AbstractIn natural environments, fluid density and viscosity can be affected by spatial and temporal variations of solute concentration, for example, due to saltwater intrusion in coastal aquifers, leachate infiltration from waste disposal sites, and upconing of saline water from deep aquifers. Potentially unstable situations may arise in which a dense fluid overlies a less dense fluid. This situation can produce instabilities manifested by dense plume fingers moving vertically downwards counterbalanced by vertical upward flow of the less dense fluid. The resulting free convection increases solute transport rates over large distances and times relative to constant-density flow. Unstable brine flow is further complicated if the porous medium is variably saturated. The results from a laboratory experiment of variably saturated variable-density flow and solute transport from Simmons et al. (2002) are used as the physical basis to define a new mathematical benchmark. This benchmark aims at realistically reproducing the experimental fingering patterns. Random hydraulic conductivity fields were used in the simulations as a numerical perturbation method to realistically mimic the observed dense plume fingering. The HydroGeoSphere code coupled with PEST are used to calibrate the parameter set that defines the benchmark. A grid convergence analysis is performed to obtain the adequate spatial and temporal discretizations. The new mathematical benchmark is useful for model comparison and testing of variably saturated variable-density flow in porous media. Simmons CT, Pierini ML, Hutson JL (2002) Laboratory investigation of variable-density flow and solute transport in unsaturated–saturated porous media. Transp Porous Media. 47(2): 215–244, 10.1023/A:1015568724369.

Ce travail de thèse expose le développement de deux modèles numériques destinés à résoudre le transport de masse en milieux poreux et dans une fracture unique lisse ou rugueuse. Le modèle pour milieu poreux décrit, à l’échelle microscopique, les processus de transport d'une particule soluble que sont la convection, la dispersion et la diffusion moléculaire. La convection et la dispersion sont décrites au travers d'un vecteur (dans les trois composantes de l'espace) dont le module et la direction sont liés à des lois probabilistes destinées à rendre compte de la nature aléatoire d'un milieu poreux. La diffusion moléculaire est simulée par un bruit aléatoire ajouté aux composantes de ce vecteur disperso-convectif. Une telle description des processus de transport permet de simuler de nombreux cas hydrogéologiques de façon souple : interaction solute/matrice ou transport de colloides. Le modèle en fracture autorise le transport de solutes et de colloides dans plusieurs cas de figure : en 2d ou 3d et pour une fracture lisse ou rugueuse, la rugosité de la fracture étant générée par une méthode fractale. Le transport des particules se fait par l’intermédiaire d'une marche aléatoire, la distribution des vitesses de la fracture suivant un profil de poiseuille. La diffusion moléculaire, calculée en fonction de la taille des particules par la loi de Stokes-Einstein, permet aux particules d’échantillonner l'ensemble des vitesses du profil. Plusieurs cas de transport sont traités : interaction particules/parois ou déplacement simultané de particules de différentes tailles. Outre des vérifications numériques par comparaison à des solutions analytiques, ces deux modèles sont confrontés à diverses expériences de la littérature et du projet européen Caress (1996-1999), cadre technique de cette thèse.

Ce travail traîte de l'écoulement et du transport réactif de soluté dans des milieux poreux fracturés à l'échelle mésoscopique par simulation numérique directe dans un modèle de fractures discrètes (DFM). Le modèle est tridimensionnel. Le schéma numérique est basé sur une formulation aux volumes finis, avec mise en œuvre d'un limiteur de flux de type Superbee afin de minimiser la dispersion numérique. Une attention particulière est portée à la formulation du triple volume de contrôle et des échanges entre la matrice poreuse et les fractures qui permet une description précise de l'écoulement et des processus de transport, même lorsque le contraste entre les transmissivités des fractures et de la matrice poreuse est grand. Diverses situations physiques et applications sont étudiées, de l'analyse de cas simples comme le transport passif ou réactif dans la matrice non fracturée, le transfert entre une fracture - en forme de disque – et de la matrice poreuse à des cas plus complexes, comme l'induction entre fractures dans le milieu poreux. Enfin, des simulations complètes de transport de soluté passif et réactive dans de longs échantillons de la matrice poreuse contenant des réseaux de fractures percolants ou non percolants sont réalisées et analysées
This work addresses flow and reactive solute transport in fractured porous media on a mesoscopic scale by direct numerical simulations in a discrete fracture model (DFM). The model is three-dimensional. The numerical scheme is based on a finite volume formulation, with implementation of a flux limiting scheme of the SuperBee type in order to minimize numerical dispersion. Particular attention is paid to the triple control volume formulation and exchanges between the porous matrix and the fractures which allow for a precise description of the flow and the transport processes even when the contrast between the fractures and porous matrix transmissivities is high. Various physical situations and applications are studied, ranging from the analysis of the simple cases of passive and reactive transport in the unfractured matrix, transfer between disk shaped fracture to the embedding porous matrix, to more complex ones such as induction between fractures in the porous medium. Finally, full simulations of passive and reactive solute transport in long rods of porous matrix containing percolating or non percolating fracture networks are conducted and analyzed

Comprendre l'altération du verre nucléaire dans un réseau de fracture au sein d'un bloc de verre vitrifié est important pour la sûreté du conditionnement des déchets nucléaires (quantification des risques associés au relâchement des radionucléides). L’évaluation de la performance du stockage géologique des déchets nucléaires passe obligatoirement par la modélisation de l’altération aqueuse d’un bloc de verre nucléaire fracturé, l’échelle de temps envisagée (plusieurs milliers d’années) dépassant toute possibilité d’expérience directe. Cette thèse vise donc à combler le fossé entre les simulations d'écoulement et de transport à l'échelle du réservoir et la modélisation à l'échelle micrométrique des processus interfaciaux verre-eau, en apportant l'évaluation quantitative de la dégradation aqueuse du verre à l'échelle d’un bloc.Pour aborder ce problème, les objectifs principaux de cette thèse ont été fixés comme suit : (i) la reproduction des résultats expérimentaux obtenus précédemment (pour quelques fractures modélisées de manière discrète en mode diffusif), (ii) l’analyse de l'impact des géométries de fractures sur la quantité de verre altéré pour quelques fissures modélisées de manière discrète, (iii) l’étude de la possibilité d'adaptation du modèle géochimique à la modélisation dans le cadre de l’approche milieu équivalent, (iv) la mise au point d'une méthodologie de caractérisation, (v) la modélisation géostatistique et géométrique de réseau de fractures à l’échelle d’un conteneur de verre, (vi) le calcul des paramètres équivalents diffusifs, hydrauliques et les paramètres qui contrôlent la cinétique de dissolution de verre, et au final, (vii) la modélisation de transport réactif à l’échelle d’un conteneur.À titre illustratif, la méthodologie de la caractérisation de réseau fracturé proposée, basée sur le traitement des images, a été appliquée aux images bidimensionnelles (2D) de haute résolution de deux blocs de verre. Cette application a permis de mettre en œuvre à la fois les données directes obtenues par mesures des paramètres d’un réseau fracturé de verre vitrifié et les données indirectes explicatives issues des simulations thermomécaniques. L’application a abouti à la création de multiples réalisations de tessellation de réseaux fracturés équivalents qui ont ensuite été utilisées comme représentations physiques pour les calculs de la perméabilité équivalente, de la diffusion équivalente et des paramètres contrôlant la cinétique de dissolution de verre borosilicaté. L'évolution de la quantité de verre altéré obtenue en effectuant la modélisation de transport réactif appliquée à plusieurs réalisations de la tessellation de réseau fracturé équivalent a été comparée aux données expérimentales d’un essai d'altération aqueuse d'un conteneur non radioactif de verre nucléaire. Les résultats montrent que la méthodologie conçue offre une opportunité pour mieux comprendre l'impact de la fracturation sur l'altération aqueuse du verre vitrifié et constitue un outil fiable permettant de prendre en compte différents scénarios d'évolution du stockage
Understanding the alteration of nuclear glass in a fracture network of a vitrified glass block is important for the safe conditioning of nuclear waste (quantification of the risks associated with radionuclide release). Performance assessment of geological nuclear waste repositories entails modelling of the long-term evolution of the fractured nuclear glass block aqueous alteration, because the considered time scale, of several thousands of years, is beyond the range of any direct experimental perspectives. This dissertation aims then to bridge the gap between the reservoir-scale flow and transport simulations and the micron-scale modeling of the glass-water interfacial processes, by bringing the quantitative evaluation of the glass aqueous degradation at the block scale.To tackle this issue, the main objectives of this thesis were fixed as follows: (i) reproduction of the experimental results previously obtained (for some fractures modeled in a discrete way in the diffusive mode),(ii) analysis of the impact of fractures geometries on the quantity of altered glass at the scale of some fractures modeled in a discrete way, (iii) investigation of the possibilities of the geochemical model adaptation for the equivalent homogenous modeling, (iv) establishment of a methodology for glass block fracture network characterization, (v) geostatistical and geometric modeling, (vi) calculation of the equivalent diffusive, hydraulic and glass dissolution kinetics controlling properties and (vii) upcoming reactive transport modeling at the scale of one canister.As an illustrative example, the proposed image processing-based fracture network characterization methodology was applied to two-dimensional (2D) high-resolution images of two blocks of vitrified glass. This application brought into service both hard data obtained by direct measurement of the fracture network and soft physics-based explanatory data and resulted in the creation of multiple realizations of fracture network equivalent tessellation that were further used as physical representation for the calculation of the equivalent hydraulic, diffusive, and alteration kinetics - controlling properties. The evolution of the quantity of altered glass obtained by conducting reactive transport modeling applied to several realizations of the equivalent fracture network tessellation was compared with the experimental data of the aqueous alteration test of a non-radioactive full-scale nuclear glass canister. The results show that implementation of the devised procedure presents an opportunity for better understanding the impact of fracturing on aqueous alteration of borosilicate glass and provides a reliable tool enabling different scenarios of repository evolution to be accounted for

Comprendre l'altération du verre nucléaire dans un réseau de fracture au sein d'un bloc de verre vitrifié est important pour la sûreté du conditionnement des déchets nucléaires (quantification des risques associés au relâchement des radionucléides). L’évaluation de la performance du stockage géologique des déchets nucléaires passe obligatoirement par la modélisation de l’altération aqueuse d’un bloc de verre nucléaire fracturé, l’échelle de temps envisagée (plusieurs milliers d’années) dépassant toute possibilité d’expérience directe. Cette thèse vise donc à combler le fossé entre les simulations d'écoulement et de transport à l'échelle du réservoir et la modélisation à l'échelle micrométrique des processus interfaciaux verre-eau, en apportant l'évaluation quantitative de la dégradation aqueuse du verre à l'échelle d’un bloc.Pour aborder ce problème, les objectifs principaux de cette thèse ont été fixés comme suit : (i) la reproduction des résultats expérimentaux obtenus précédemment (pour quelques fractures modélisées de manière discrète en mode diffusif), (ii) l’analyse de l'impact des géométries de fractures sur la quantité de verre altéré pour quelques fissures modélisées de manière discrète, (iii) l’étude de la possibilité d'adaptation du modèle géochimique à la modélisation dans le cadre de l’approche milieu équivalent, (iv) la mise au point d'une méthodologie de caractérisation, (v) la modélisation géostatistique et géométrique de réseau de fractures à l’échelle d’un conteneur de verre, (vi) le calcul des paramètres équivalents diffusifs, hydrauliques et les paramètres qui contrôlent la cinétique de dissolution de verre, et au final, (vii) la modélisation de transport réactif à l’échelle d’un conteneur.À titre illustratif, la méthodologie de la caractérisation de réseau fracturé proposée, basée sur le traitement des images, a été appliquée aux images bidimensionnelles (2D) de haute résolution de deux blocs de verre. Cette application a permis de mettre en œuvre à la fois les données directes obtenues par mesures des paramètres d’un réseau fracturé de verre vitrifié et les données indirectes explicatives issues des simulations thermomécaniques. L’application a abouti à la création de multiples réalisations de tessellation de réseaux fracturés équivalents qui ont ensuite été utilisées comme représentations physiques pour les calculs de la perméabilité équivalente, de la diffusion équivalente et des paramètres contrôlant la cinétique de dissolution de verre borosilicaté. L'évolution de la quantité de verre altéré obtenue en effectuant la modélisation de transport réactif appliquée à plusieurs réalisations de la tessellation de réseau fracturé équivalent a été comparée aux données expérimentales d’un essai d'altération aqueuse d'un conteneur non radioactif de verre nucléaire. Les résultats montrent que la méthodologie conçue offre une opportunité pour mieux comprendre l'impact de la fracturation sur l'altération aqueuse du verre vitrifié et constitue un outil fiable permettant de prendre en compte différents scénarios d'évolution du stockage
Understanding the alteration of nuclear glass in a fracture network of a vitrified glass block is important for the safe conditioning of nuclear waste (quantification of the risks associated with radionuclide release). Performance assessment of geological nuclear waste repositories entails modelling of the long-term evolution of the fractured nuclear glass block aqueous alteration, because the considered time scale, of several thousands of years, is beyond the range of any direct experimental perspectives. This dissertation aims then to bridge the gap between the reservoir-scale flow and transport simulations and the micron-scale modeling of the glass-water interfacial processes, by bringing the quantitative evaluation of the glass aqueous degradation at the block scale.To tackle this issue, the main objectives of this thesis were fixed as follows: (i) reproduction of the experimental results previously obtained (for some fractures modeled in a discrete way in the diffusive mode),(ii) analysis of the impact of fractures geometries on the quantity of altered glass at the scale of some fractures modeled in a discrete way, (iii) investigation of the possibilities of the geochemical model adaptation for the equivalent homogenous modeling, (iv) establishment of a methodology for glass block fracture network characterization, (v) geostatistical and geometric modeling, (vi) calculation of the equivalent diffusive, hydraulic and glass dissolution kinetics controlling properties and (vii) upcoming reactive transport modeling at the scale of one canister.As an illustrative example, the proposed image processing-based fracture network characterization methodology was applied to two-dimensional (2D) high-resolution images of two blocks of vitrified glass. This application brought into service both hard data obtained by direct measurement of the fracture network and soft physics-based explanatory data and resulted in the creation of multiple realizations of fracture network equivalent tessellation that were further used as physical representation for the calculation of the equivalent hydraulic, diffusive, and alteration kinetics - controlling properties. The evolution of the quantity of altered glass obtained by conducting reactive transport modeling applied to several realizations of the equivalent fracture network tessellation was compared with the experimental data of the aqueous alteration test of a non-radioactive full-scale nuclear glass canister. The results show that implementation of the devised procedure presents an opportunity for better understanding the impact of fracturing on aqueous alteration of borosilicate glass and provides a reliable tool enabling different scenarios of repository evolution to be accounted for

Cette thèse porte sur la simulation numérique de modèles couplés pour l'écoulement et le transport dans les milieux poreux. Nous présentons une nouvelle méthode de couplage entre les réactions chimiques et le transport en utilisant une méthode de Newton-Krylov, et nous étudions également un modèle d'écoulement en milieu fracturé qui traite l'intersection des fractures par une méthode de décomposition de domaine.
Ce travail est divisé en trois parties : la première partie contient une analyse de différents schémas numériques pour la discrétisation des problèmes d'advection-diffusion, notamment par une technique de séparation d'opérateurs, ainsi que leur mise en oeuvre informatique, dans un code industriel.
La deuxième partie, qui est la contribution majeure de cette thèse, est consacrée à la modélisation et à l'implémentation d'une méthode de couplage globale pour le transport réactif. Le système couplé transport-chimie est décrit, après discrétisation en temps, par un système d'équations non linéaires. La taille du système sous-jacent, à savoir le nombre de points de grille multiplié par le nombre d'espèces chimiques, interdit la résolution du système linéaire par une méthode directe. Pour remédier à cette difficulté, nous utilisons une méthode de Newton-Krylov qui évite de former et de factoriser la matrice Jacobienne.
Dans la dernière partie, nous présentons un modèle d'écoulement dans un milieu fracturé tridimensionnel, basé sur une méthode de décomposition de domaine, et qui traite l'intersection des fractures. Nous démontrons l'existence et l'unicité de la solution, et nous validons le modèle par des tests numériques.

Les milieux poreux fracturés sont des milieux composés d'une roche présentant des zones de fracturation. Géologiquement, ces milieux sont caractérisés par la présence de fractures sur plusieurs échelles avec des propriétés hétérogènes créant des zones fortement perméables à comparer de la roche environnante. Hydrauliquement, ces milieux sont caractérisés par des temps de réponses courts correspondant aux sollicitations des structures fortement perméables et des temps de réponses longs correspondant aux sollicitations des structures faiblement perméables. Ce type de milieux est impliqué dans de nombreux enjeux sociétaux tels que l'étude de sites contaminés, le stockage d'éléments et l'exploitation de ressources. Le principal défi de la modélisation des milieux poreux fracturés réside en la représentation des hétérogénéités géométriques et physiques caractérisant ces milieux. Une représentation exacte du milieu naturel étant impossible, il s'agit de déterminer quelles sont les propriétés caractéristiques de ces milieux, c'est-à-dire les propriétés principales "responsables" de leur comportement. Mes travaux de thèse ont consisté à étudier l'impact des propriétés géométriques et physiques des fractures et de la matrice de l'échelle locale de la fracture et du bloc matriciel à l'échelle globale du réseau de fractures. Une première partie du travail correspond à la mise au point de méthodes d'évaluation des effets des structures et de leurs propriétés sur les échanges et une seconde partie en l'exploitation de ces méthodes. Au final, je propose un modèle double-porosité discret prenant en compte les propriétés du milieu identifiées comme caractéristiques de son comportement.

Le point de départ de cette thèse est le constat suivant : les modèles prédictifs de transfert dans les milieux souterrains, souvent basés sur une représentation de type milieu poreux homogène, ne sont pas adaptés aux milieux naturels, qui sont caractérisés au contraire par une forte hétérogénéité des vitesses écoulement à toutes les échelles. La grande incertitude qui existe actuellement sur la prédiction des transferts et des temps de résidence dans ces systèmes ne provient pas seulement de l'incertitude sur la paramétrisation des modèles mais également de l'incertitude sur leur formulation même, liée aux hypothèses effectuées sur la représentation des distributions de vitesses d'écoulement. Pour aborder ce problème, les objectifs principaux de cette thèse étaient, d'une part de développer des méthodes de mesure in situ adaptées à la caractérisation de l'hétérogénéité des écoulements et de leur organisation, et d'autre part de mener une étude théorique sur la représentation de l'organisation des écoulements dans les modèles prédictifs. Cette thèse a été effectuée dans le cadre des projets européens SALTRANS (2000-2004) et ALIANCE (2001-2005) et du réseau national de sites hydrogéologiques de l'Observatoire de Recherche en Environnement H+.
Trois résultats principaux ont été obtenus. (1) Nous avons tout d'abord apporté une démonstration du caractère fractal des écoulements dans les aquifères fracturés à partir de tests hydrauliques. Cette approche permet de définir un modèle équivalent permettant de rendre compte du comportement hydraulique moyen d'un site naturel sur une large gamme d'échelles de temps et d'espace. (2) Du point de vue de la caractérisation expérimentale de la variabilité des vitesses d'écoulement, nous avons mis en oeuvre des méthodes de mesure de vitesse d'écoulement en forage, en collaboration avec Frederick Paillet (University du Maine, USA). Nous avons proposé et testé une méthode d'inversion de ces mesures permettant d'imager les chemins d'écoulement et de quantifier leurs propriétés hydrauliques à différentes échelles sur les sites naturels. (3) Enfin, à partir d'une étude numérique destinée à simuler les écoulements et le transport de soluté dans des champs de perméabilité hétérogènes, nous avons proposé une nouvelle méthode permettant de quantifier l'organisation des écoulements. Les résultats obtenus concernant l'analyse détaillée de la structure de corrélation du champ de vitesse permettent de mettre en question les hypothèses utilisées dans les modèles de transport de solutés, en vue de progresser vers une représentation plus réaliste des distributions de vitesse dans ces modèles.

Dans le cadre de la problématique de l'enfouissement des déchets nucléaires, la barrière géologique constitue la dernière zone de transfert des radio-éléments. Cette barrière pouvant être fracturée des recherches sont menées pour faire progresser la modélisation des transferts en milieu fracturé. La complexité du milieu, les contrastes existant entre les différentes zones et les incertitudes liées aux propriétés physiques rendent cette tâche complexe. De plus, après fermeture du site, les écoulements lents dans le milieu favorisent des phénomènes diffusifs dans la roche qui contribuent à augmenter le temps de transit des radio-éléments. Dans ce contexte, une approche Smeared Fractures a été développée pour un schéma en Eléments Finis Mixtes Hybrides et implémentée dans le code Cast3M. Cette approche permet une représentation explicite des fractures principales alors que la fracturation de plus petite échelle est homogénéisée. L'utilisation d'un maillage régulier permet, en outre, d'éviter un maillage explicite coûteux. La présence des fractures est prise en compte par un champ hétérogène de propriétés. Ces propriétés sont affectées de manière à respecter les critères de conservation des flux (hydraulique et massique) à l'échelle de la fracture. Pour l'écoulement, l'approche Smeared Fractures présente des performances comparables à celles obtenues avec des approches discrètes tandis que le caractère 3D de la géométrie des blocs matriciels est respecté lors de la résolution du transport. Le choix des discrétisation spatiale et temporelle doit respecter des critères qui ont été établis. Néanmoins, à l'intérieur de ces limites, et suivant la précision désirée, l'approche permet de réduire les temps de calcul. Les résultats de validation et de qualification de l'approche appliquée à des géométries 2D et 3D, synthétiques et réalistes, sont présentés pour différents jeux de paramètres physiques. Des applications de l'approche au site d'Äspö (Suède) clôturent ce travail.