Modélisation 3D du transport particulaire asynchrone en simple et double continuum matrice-fractures : application au stockage de déchets nucléaires

• Main Author: Lam, Minh-Phuong
• Format: Others
• Published: 2008
• Online Access: http://oatao.univ-toulouse.fr/7738/1/lam.pdf

Ce mémoire présente les travaux réalisés dans le cadre d'une thèse CIFRE en collaboration avec le Laboratoire National d'Hydraulique et Environnement (LNHE) d'EDF R et D et l'Institut de Mécanique des Fluides de Toulouse (IMFT). La thèse a pour objectif in fine de fournir au LNHE un outil numérique précis lui permettant d'étudier la fiabilité d'un stockage de déchets nucléaires en milieu géologique profond. Les travaux réalisés portent plus particulièrement sur le développement et l'implémentation d'une méthode de suivi de particules par marche aléatoire pour modéliser le phénomène de transport 3D dans des milieux poreux hétérogènes (fracturés). Dans une première partie, le rapport présente la méthode particulaire utilisée pour modéliser le transport en milieu poreux, avec une approche discrète haute résolution. Celle-ci s'appuie sur la représentation de l'évolution d'un panache de concentration par le déplacement d'un nuage de particules. Les particules de traceur se déplacent selon un schéma lagrangien asynchrone, implémenté sur un modèle hydraulique 3D, sur un maillage non structuré (éléments finis tétraédriques). L'utilisation de pas de temps asynchrones permet d'éviter, en maillage non structuré, le phénomène de ‘dépassement' inhérent aux autres méthodes particulaires. En pré-traitement, afin de lisser les discontinuités des paramètres de transport, une méthode de filtrage est mise en place. En post-traitement, les concentrations de traceurs sont obtenues à partir des positions de particules en utilisant plusieurs méthodes de filtrage et d'échantillonnage de particules. Les méthodes développées sont illustrées et validées via des problèmes d'advection-diffusion de traceur en milieux homogènes et hétérogènes. Le modèle de suivi de particules à marche aléatoire asynchrone proposé s'avère efficace et satisfaisant. Les tests de validations permettent d'optimiser les paramètre de filtrage et donnent des résultats très proches des solutions de référence (analytique ou à partir de benchmark) dans les cas optimisés. Lorsque les milieux abordés possèdent des hétérogénéités spatialement denses, comme dans le cas des milieux fracturés, le modèle présenté dans la première partie est complété par une nouvelle approche permettant de représenter de façon adéquate le milieu étudié. La seconde partie du travail consiste donc à étendre le modèle particulaire au cas d'un milieu fracturé représenté à l'aide d'une approche double-continuum, dans laquelle la matrice poreuse et les fractures sont représentées par un continuum Matrice et un continuum Fracture. Les propriétés hydrogéologiques qui interviennent dans les équations d'écoulement et de transport des continua sont définies par des paramètres équivalents homogénéisés. Les deux continua sont couplés par des échanges de fluides et de traceurs. Les échanges inter-continuum de traceurs sont représentés par des échanges de particules et sont traités sous la forme de probabilités de transition de particules d'un continuum à l'autre. Des tests sur le modèle de suivi de particules avec l'approche double-continuum sont réalisés sur un réseau régulier de fractures. Les tests effectués pour des continua de contrastes modérés montrent l'adéquation de la méthode double continuum pour cette classe de modèles.