La protection et la remédiation des ressources en eau sont un enjeu sociétal majeur, ainsi il est nécessaire de comprendre l’évolution de solutés, tels les polluants, au sein de la zone saturée et non saturée. Dans ce but, de nombreux travaux ont été consacrés à la modélisation du transport réactif en milieux poreux. Son déroulement à l’échelle de Darcy dépend des hétérogénéités microscopiques du milieu. Les modèles de réseau de pores qui simplifient la géométrie en un ensemble de pores reliés par des liens de sections constantes, permettent de se placer à une échelle mésoscopique, faisant le lien entre l’échelle porale et l’échelle de Darcy. Sur de telles formes géométriques, l’écoulement admet un traitement analytique. En ce qui concerne le transport réactif des solutés, nous proposons une solution analytique dans les liens qui permet de calculer le débit de masse entre pores. Le modèle de transport se formule alors comme un système d’équations de Volterra de secondes espèces dont les noyaux de convolution sont des séries d’exponentielles décroissantes (hormis le premier terme qui est constant). Leurs temps de relaxation sont pilotés essentiellement par le temps de dispersion td. Dans la limite où td tend vers 0 à Péclet constant, les termes transitoires des noyaux se réduisent à un Dirac, débouchant sur un premier modèle simplifié à réponse instantanée c'est-à-dire un modèle de transport quasi-statique. Dans le cas où les volumes des pores sont suffisamment grands, les noyaux se réduisent à leur premier terme. Ces formulations du transport généralisent celles de la littérature. En particulier pour des Péclet petit ou grand on retrouve respectivement les modèles usuels en régime dispersif et convectif. Numériquement, la décroissance exponentielle des noyaux permet d’optimiser le calcul des convolutions avec une précision arbitrairement fixée, réduisant drastiquement le temps de résolution. / Protection and remediation of ground water resources are a major societal challenge. It implies to understand the evolution of solutes as pollutant in the saturated and non-saturated zones. For that purpose numerous studies have been conducted for modeling the reactive transport in a porous media. At Darcy scale, the behavior of solutes depends on microscopic heterogeneity for the media. The Pore Network Models (PNM) simplifies drastically its geometry and considers pores linked by straight throats the section of which is constant. They give a description which is in between the macroscopic and the pore descriptions. With such geometry it is possible to use a Poiseuille flow modeling the flux. With respect to the reactiontransport equation, we seek the analytical solution of the CDE in throats, which in turn allows computing the mass flux in pores. The transport solution consists of a Volterra equation system. Its convolution kernels result in a summation of time function which is decreasing exponentially with time (except the first term which still constant). The time constant is driven by the diffusion time td. As td goes to zero, keeping the Peclet number fixed, each term of the summation reduces to a Dirac. The response of the system is then instantaneous. When the volume of the pore is large enough it is possible to neglect all the term of the kernel except the constant one. In the limit where the Peclet number goes to zero, usual models are recovered. Numerically, the exponential time decreasing of the kernel allow to optimize their computational time up to an arbitrary fixed precision.

Cette thèse propose d'étudier en détail une méthode générale de construction d'un réseaude pores interconnectés à partir d'une image 3D d'un matériau poreux afin de calculer laperméabilité absolue. Un squelette de l'espace poral est tout d'abord utilisé pour définir laposition des pores. Ce squelette est transformé en graphe, puis plusieurs étapes de fusiondes nœuds sont réalisées pour obtenir un réseau exploitable. Une ligne de partage des eauxprenant comme marqueurs les nœuds du graphe permet d'attribuer une géométrie auxpores. Le système d'équation linéaire à résoudre est construit en intégrant les équations deStokes sur des éléments de l'espace poral. Une méthode originale de calcul du coefficient deproportionnalité existant entre les pressions au centre de deux pores et le débit traversantla surface les séparant est proposée. Une application de l'intégralité de l'approche estégalement présentée sur un matériau réel.

Cette thèse propose d'étudier en détail une méthode générale de construction d'un réseaude pores interconnectés à partir d'une image 3D d'un matériau poreux afin de calculer laperméabilité absolue. Un squelette de l'espace poral est tout d'abord utilisé pour définir laposition des pores. Ce squelette est transformé en graphe, puis plusieurs étapes de fusiondes nœuds sont réalisées pour obtenir un réseau exploitable. Une ligne de partage des eauxprenant comme marqueurs les nœuds du graphe permet d'attribuer une géométrie auxpores. Le système d'équation linéaire à résoudre est construit en intégrant les équations deStokes sur des éléments de l'espace poral. Une méthode originale de calcul du coefficient deproportionnalité existant entre les pressions au centre de deux pores et le débit traversantla surface les séparant est proposée. Une application de l'intégralité de l'approche estégalement présentée sur un matériau réel. / A general method to build a pore network model from a 3D image of a porous material ispresented in this work in order to compute its permeability. A skeleton of the void spaceis first used to define pores position. This skeleton is converted to a first graph. Mergingsteps are necessary to obtain a relevant network. A volume is allocated to the pores using awatershed algorithm, starting with markers defined from the nodes of the graph. The linearsystem of equations to solve is deduced by integrating the Stokes equations on pore spacepartitioning elements. The proportionality coefficient linking the pressure drop betweentwo pore centers and the flow going through the surface separating them is calculated withan original method. Application of the complete process to a real material is presented.