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

Michel Quintard

Description

Niveau: Supérieur, Doctorat, Bac+8
N° d'ordre :……………… THESE présentée pour obtenir LE TITRE DE DOCTEUR DE L'INSTITUT NATIONAL POLYTECHNIQUE DE TOULOUSE École doctorale : Science de l'Univers, de l'Environnement et de l'Espace Spécialité : Modélisation en Hydrogéologie et Sûreté Nucléaire Par Minh-Phuong LAM Titre de la thèse Modélisation 3D du transport particulaire asynchrone en simple et double continuum matrice-fractures : Application au stockage de déchets nucléaires Soutenue le 05 Juin 2008 devant le jury composé de : M. Michel QUINTARD Président M. Rachid ABABOU Directeur de thèse M. Amvrossios BAGTZOGLOU Rapporteur M. Frederic DELAY Rapporteur M. Francisco Javier ELORZA TENREIRO Rapporteur Mme Regina NEBAUER Membre M. Jean-Michel HERVOUET Membre M. Iraj MORTAZAVI Membre

• moments de folies douces…


• application au stockage de déchets nucléaires


• science de l'univers, de l'environnement et de l'espace spécialité

Sujets

• Institut national polytechnique de Toulouse
• Karlsruhe Institute of Technology
• Vincent Hervouët
• Université polytechnique de Madrid
• Renaud

Informations

• Publié par: profil-zyak-2012
• Publié le: 01 juin 2008
• Nombre de lectures: 68
• Langue: Français
• Poids de l'ouvrage: 6 Mo

Extrait

N° d’ordre :………………





THESE


présentée
pour obtenir

LE TITRE DE DOCTEUR DE L’INSTITUT NATIONAL POLYTECHNIQUE DE
TOULOUSE


École doctorale : Science de l’Univers, de l’Environnement et de l’Espace
Spécialité : Modélisation en Hydrogéologie et Sûreté Nucléaire


Par Minh-Phuong LAM


Titre de la thèse

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



Soutenue le 05 Juin 2008 devant le jury composé de :


M. Michel QUINTARD Président
M. Rachid ABABOU Directeur de thèse
M. Amvrossios BAGTZOGLOU Rapporteur
M. Frederic DELAY
M. Francisco Javier ELORZA TENREIRO Rapporteur
Mme Regina NEBAUER Membre
M. Jean-Michel HERVOUET Membre
M. Iraj MORTAZAVI Membre ii iii


























A ma famille,



iv REMERCIEMENTS



Cette thèse a été réalisée dans les locaux du Laboratoire National d’Hydraulique et
Environnement (LNHE) à Chatou, au sein de la division Recherche et Développement d’EDF
et en collaboration avec l’Institut de Mécanique des Fluides de Toulouse (IMFT) et l’Institut
National Polytechnique de Toulouse (INPT).

Mes premiers remerciements vont à mon directeur de thèse Rachid Ababou pour avoir accepté
d’encadrer mes travaux et pour ses conseils avisés. Un grand merci également à Regina
Nebauer pour son encadrement efficace, son soutien dans les moments difficiles et sa grande
patience tout au long de ces années chez EDF.

Je tiens aussi à remercier Israël Cañamon de l’Université Polytechnique de Madrid,
JeanPhilippe Renaud de l’Université de Bristol, Matthias Beyer de l’Université de Karlsruhe et
Gérald Debenest de l’IMFT: nos échanges et débats d’idées m’ont beaucoup aidé.

Merci également à Messieurs Amvrossios Bagtzoglou, Frederic Delay, Fco. Javier Elorza
Tenreiro, Michel Quintard, Jean-Michel Hervouet et Iraj Mortazavi pour leur participation au
jury de cette thèse.

Ces trois années passées au LNHE m’ont parues bien courtes, tant l’ambiance chaleureuse et
conviviale a fait son effet. Activités sportives, repas « gastronomiques » et Comité des Fêtes
sont aussi des mots que j’associerai à mes souvenirs de thèse… A tous les jeunes et moins
jeunes, stagiaires et thésards, agents et prestataires: merci pour tous ces fous rires et ces
moments de folies douces… capturés sur pellicules!

Plus personnellement, je tiens à remercier ma petite « bande bordelaise » pour les soirées et
délires en tout genre et pour l’amitié sans faille partagée depuis notre enfance.

Et bien sûr une tendre pensée pour ma famille pour son soutien indéfectible.

v RESUME



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&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
intercontinuum 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.

Particules/Marche aléatoire/ Asynchrones/ Maillage non-structuré/ Filtrage,
Fractures/ Double-continuum/ Paramètres équivalents homogénéisés/ Echanges

vi ABSTRACT



This PhD research was conducted as a collaboration between Laboratoire National
d’Hydraulique et Environnement (LNHE) from EDF R&D and the Institut de Mécanique des
Fluides de Toulouse (IMFT) in the frame of a CIFRE contract. This PhD thesis aims at
providing LNHE a reliable numerical model to study the feasibility of a nuclear waste storage
in deep geological structures. The main focus of the thesis is put on developing and
implementing a Random Walk Particle Method (RWPM) to model contaminant transport in
3D heterogeneous and fractured porous media.

In its first part, the report presents the lagrangian particle tracking method used to model
transport in heterogeneous media with a direct high resolution approach. The solute plume is
discretized into concentration packets: particles. The model tracks each particle based on a
time-explicit displacement algorithm according to an advective component and a diffusive
random component. The method is implemented on a hydraulic model discretized on a 3D
unstructured tetrahedral finite element mesh. We focus on techniques to overcome problems
due to the discontinuous transport parameters and the unstructured mesh. First, we introduce
an asynchronous time-stepping approach to deal with the numerical and overshoot errors that
occur with conventional RWPM. Then, a filtering method is applied to smooth discontinuous
transport parameters (pre-proce