L'hétérogénéité spatiale est un trait caractéristique des aquifères calcaires. Conditionnée par le karst et ses morphologies, cette hétérogénéité explique la fréquente sous-exploitation des réserves de ces hydrosystèmes mais également le degré lacunaire de connaissance de leurs fonctionnements. Nous avons voulu, dans cette recherche, approfondir cette connaissance en scrutant le fonctionnement des aquifères calcaires qui arment les bas plateaux du Nord-Est de la France. Trois de ces bas plateaux ont été retenus, ceux du Dogger et de l'Oxfordien (Jurassique) et celui du Muschelkalk (Trias). Afin de comprendre comment s'organisent les écoulements sur et dans ces trois systèmes, notre démarche s'articule autour de deux principaux axes ; le premier concerne l'identification des modalités de l'écoulement et le second est consacré à leur quantification. En pratique, la connaissance des modalités de circulation a été recherchée à travers l'étude des différents paramètres structurels (géologique, tectonique, karstique) et conjoncturels (précipitation, évapotranspiration) qui interagissent dans le cycle de l'eau. Nous montrons que la variété de fonctionnements hydrologiques observés, à l'échelle du domaine d'étude, est conditionnée principalement par les paramètres structurels ; l'hétérogénéité lithologique, notamment, apparaît comme une clef fonctionnelle déterminante, d'une part, par l'existence d'une multiplicité de nappes d'eau et d'autre part, par son rôle dans la karstogenèse. Les morphologies du karst s'observent principalement sur les secteurs de contact lithostratigraphique et nous avons pu établir, à travers une répartition spatiale des phénomènes du karst (en surface et en souterrain), que ces trois plateaux sont majoritairement déconnectés du karst et dominés par des processus de circulation lents, de l'infiltration à la vidange aquifère. Afin d'apporter des éléments de quantification sur les circulations, une analyse hydroclimatique et géochimique d'exutoires de petites dimensions a été ensuite menée. C'est ainsi qu'une petite dizaine de sources et de bassins versants, localisés dans des contextes physiographiques variés, ont fait l'objet de suivis sur plusieurs cycles hydrologiques. L'analyse nous a permis de préciser les caractères hydrodynamiques de deux grands types de comportement observés au sein des bas plateaux étudiés : le comportement fissural et le comportement karstique. Le comportement fissural se manifeste, à l'échelle annuelle, par un filtrage important des précipitations par les systèmes (forte inertie des valeurs de débits et de conductivités électriques) et également par de lents processus de mise en charge à l'échelle des évènements de crue. Le comportement karstique, quant à lui, est attesté, à l'échelle annuelle, par des pluies moyennement filtrées par les systèmes mais aussi par des variations hydrogéochimiques accusées lors des sollicitations pluviométriques majeures, en lien avec la transmissivité importante de ces milieux / The spatial heterogeneity is a feature of limestone aquifers. Conditioned by the karst and its morphologies, this heterogeneity explains the frequent under-exploitation of these hydrosystem reserves but it also explains the incomplete degree of knowledge of their functionings. In this research we wanted to deepen this knowledge by scrutinizing the functioning of the limestone aquifers which arm the low plateaus of the North-East of France. Three of these low plateaus held our attention, those of Dogger and Oxfordian (Jurassic) and the one of Muschelkalk (Triassic). In order to understand how the flows on and into these three systems are organized, our approach focuses on two main axes ; the first one concerns the identification of the modes of flow and the second one is devoted to quantification. In practice, the knowledge of the flow modes was investigated through the study of various structural parameters (geology, tectonics, karst) and conjunctural parameters (precipitation, evapotranspiration) interacting the water cycle. We show that the variety of hydrological functionings which are observed, on the scale of the study area, is mainly dependent on the structural parameters ; the lithological heterogeneity, in particular, appears as a determining functional key, on the one hand, by the existence of a multiplicity of groundwater and on the other hand, by its role in the karstogenesis. The morphologies of the karst are primarily obvious on the areas of lithostratigraphic contact and we were able to establish, through a spatial distribution of karst phenomena (on and underground), that these three plateaus are mainly disconnected from the karst and run by slow-flow processes, from infiltration to ground water discharge. In order to provide some quantifications of the circulations, a hydroclimatic and geochemical analysis of small outlets was then conducted. Thus a dozen springs and catchments, located in various physiographic settings have been surveyed for several water years. The analysis allowed us to specify the hydrodynamic characteristics of two types of behavior observed in the low plateaus we studied: the fissural behavior and the karstic behavior. The fissural behavior appears, on an annual scale, through a significant rainfall filtering of the systems (high inertia of the values of discharges and electrical conductivities) and also by slow charging processes when flood events occur. The karst behavior, meanwhile, is evidenced on an annual scale, by a moderately rainy filtering of the systems but also by hydrogeochemical variations accused during major rainfall solicitations in connection with the important transmissivity of these environments

Afin de caractériser la nature et le rôle des particules fines et colloïdales (PFC) dans les milieux de subsurface, le suivi dans le temps des PFC contenues à la fois dans les retombées atmosphériques et dans les eaux d'infiltration et la détermination des facteurs contrôlant la distribution et l'évolution des PFC dans l'environnement a été réalisée. La démarche utilisée est basée essentiellement sur des observations in situ et consiste à prélever mensuellement, sur le terrain, les PFC contenues dans les retombées atmosphériques et dans les eaux d'infiltration d'un sol. L'étude en MET/EDX des PFC individuelles à permis de réaliser une caractérisation typologique et physico-chimiques et d'étudier leur abondance et leur réactivité face aux contaminants. Par ailleurs, la connaissance des paramètres environnementaux de chacun des sites et l'étude statistique détaillée de l'évolution de PFC a permis de proposer un modèle conceptuel de la l'évolution dynamique (distribution, réactivité) des PFC dans les environnements considérés. Les résultats obtenus en microscopie montrent que la très grande majorité des particules ont un diamètre inférieur à 0,45 µm. L'étude systématique au MET et à l'EDX des échantillons a permis de caractériser les particules présentes et de différencier plusieurs types récurrents dans des contextes environnementaux distincts, ainsi les PFC sont, (i) minérales et constituées par des argiles plus ou moins altérées, des oxy-hydroxydes et des sels (sulfures, chlorures) avec la présence ponctuelle de carbonates et de tectosilicates, ou (ii) organiques et constituées par des bactéries ou de la matière organique non-vivante (exsudats de bactérie et matière organique dégradée), ou encore (iii) des associations mixtes de particules organiques et minérales formant des micro-agrégats organo-minéraux. Au-delà de l'étude typologique, la caractérisation des propriétés physico-chimiques des différents types de PFC a permis d'en identifier les sources, l'abondance et la réactivité. L'étude de la distribution au cours du temps des PFC des retombées atmosphériques et la connaissances des facteurs environnementaux sur les deux sites étudiés a permis de mettre en évidence les facteurs clefs contrôlant la distribution des PFC. Ainsi, il apparaît, dans un premier temps, que l'intensité et la hauteur de pluie sont les facteurs principaux mis en jeu. Cependant, la localisation et le mode d'occupation des sols influent également en modifiant, en amont, les sources de PFC et en aval l'action de la pluie. En effet, si l'étude de la dynamique des PFC des retombées atmosphériques sur les deux sites montre une influence très nette de la pluie, les deux sites montrent une réponse différente. L'acteur majeur contrôlant les dépôts de PFC est la pluie, modifié par la présence de couvert végétal qui induit un enrichissement en éléments et l'agrégation des PFC sur les feuilles. L'acteur majeur contrôlant les PFC dans les eaux d'infiltration du sol étudié est la dynamique hydrique du sol et donc les fluctuations de son état de saturation en fonction de la profondeur considérée et de la pluviosité. Enfin, ce travail permet de souligner trois résultats importants que sont : la mise au point d'une méthodologie efficace, la mise en évidence d'un vecteur de contaminants essentiel et souvent négligé que sont les bactéries, via les observations in situ et des expériences en laboratoire et l'importance des mécanismes d'agrégation dans le transport des PFC.

Le réchauffement climatique est lié aux augmentations des concentrations de gaz à effet de serre dans l'atmosphère terrestre et en particulier aux émissions anthropiques de CO2. La séquestration géologique a la capacité et la longévité potentielles pour diminuer de façon significative les émissions anthropiques de CO2. Cette séquestration à grande profondeur induit des risques de fuite des réservoirs géologiques. Parmi les scénarios de fuite envisagés, celui d'une fuite diffuse est le plus inquiétant puisque sans surveillance, cette fuite pourrait perdurer et entrainer des séquelles sur l'environnement ainsi que des risques pour les populations. Des outils et protocoles de surveillance doivent donc être mis au point pour la surveillance en proche surface. Ce travail de thèse s'inscrit dans le cadre de cette problématique. Il a pour objectif la caractérisation, la quantification et la modélisation du transport et des interactions du CO2 dans une zone non saturée carbonatée. Ce travail a suivi une approche expérimentale sur un site pilote naturel à Saint-Emilion (Gironde, France), avec la réalisation de fuites diffuses en ZNS carbonatée. Cette étude aborde plusieurs thématiques: la description et l'instrumentation du site pilote naturel ; la caractérisation physico-chimique de l'hétérogénéité du réservoir carbonaté ; l'étude du fonctionnement naturel de la ZNS carbonatée et en particulier la mise en place d'une ligne de base des concentrations en CO2 ; la caractérisation de l'extension des panaches de gaz suite à des expériences de fuite diffuse dans la ZNS carbonatée et l'étude par simulation numérique des interactions gaz-eau-roche lors d'une fuite diffuse de CO2 dans une ZNS carbonatée. Les résultats de ces travaux montrent l'importance de la caractérisation de l'hétérogénéité du réservoir carbonaté ainsi que des techniques d'échantillonnage et d'analyse des différentes phases en présence. L'établissement de la ligne de base a une importance particulière pour permettre de distinguer les variations naturelles de celles induites par une fuite diffuse de CO2 dans la ZNS carbonatée. Les modes de transport du CO2 vont évoluer en fonction des paramètres physico-chimiques. Ce transport se fait par advection et/ou par diffusion. L'utilisation de gaz inertes au niveau du site de séquestration géologique est très importante puisque la détection de ces traceurs permettrait de prédire les arrivées de panaches de CO2 en proche surface. Par ailleurs, les interactions chimiques doivent être prises en compte dans les modèles de transport afin de pouvoir définir les facteurs de retard et l'impact d'une fuite diffuse de CO2 sur une ZNS carbonatée.

Le réchauffement climatique est lié aux augmentations des concentrations de gaz à effet de serre dans l'atmosphère terrestre et en particulier aux émissions anthropiques de CO₂. La séquestration géologique a la capacité et la longévité potentielles pour diminuer de façon significative les émissions anthropiques de CO₂. Cette séquestration à grande profondeur induit des risques de fuite des réservoirs géologiques. Parmi les scénarios de fuite envisagés, celui d'une fuite diffuse est le plus inquiétant puisque sans surveillance, cette fuite pourrait perdurer et entrainer des séquelles sur l'environnement ainsi que des risques pour les populations. Des outils et protocoles de surveillance doivent donc être mis au point pour la surveillance en proche surface. Ce travail de thèse s'inscrit dans le cadre de cette problématique. Il a pour objectif la caractérisation, la quantification et la modélisation du transport et des interactions du CO₂ dans une zone non saturée carbonatée. Ce travail a suivi une approche expérimentale sur un site pilote naturel à Saint-Emilion (Gironde, France), avec la réalisation de fuites diffuses en ZNS carbonatée. Cette étude aborde plusieurs thématiques: la description et l'instrumentation du site pilote naturel ; la caractérisation physico-chimique de l'hétérogénéité du réservoir carbonaté ; l'étude du fonctionnement naturel de la ZNS carbonatée et en particulier la mise en place d'une ligne de base des concentrations en CO₂ ; la caractérisation de l'extension des panaches de gaz suite à des expériences de fuite diffuse dans la ZNS carbonatée et l'étude par simulation numérique des interactions gaz-eau-roche lors d'une fuite diffuse de CO₂ dans une ZNS carbonatée. Les résultats de ces travaux montrent l'importance de la caractérisation de l'hétérogénéité du réservoir carbonaté ainsi que des techniques d'échantillonnage et d'analyse des différentes phases en présence. L'établissement de la ligne de base a une importance particulière pour permettre de distinguer les variations naturelles de celles induites par une fuite diffuse de CO₂ dans la ZNS carbonatée. Les modes de transport du CO₂ vont évoluer en fonction des paramètres physico-chimiques. Ce transport se fait par advection et/ou par diffusion. L'utilisation de gaz inertes au niveau du site de séquestration géologique est très importante puisque la détection de ces traceurs permettrait de prédire les arrivées de panaches de CO₂ en proche surface. Par ailleurs, les interactions chimiques doivent être prises en compte dans les modèles de transport afin de pouvoir définir les facteurs de retard et l'impact d'une fuite diffuse de CO₂ sur une ZNS carbonatée. / Global warming is related to atmospheric greenhouse gas concentration increase and especially anthropogenic CO₂ emissions. Geologic sequestration has the potential capacity and the longevity to significantly diminish anthropogenic CO₂ emissions. This sequestration in deep geological formation induces leakage risks from the geological reservoir. Several leakage scenarios have been imagined. Since it could continue for a long period, inducing environmental issues and risks for human, the scenario of a diffusive leakage is the most worrying. Thus, monitoring tools and protocols are needed to set up a near-surface monitoring plan. The present thesis deals with this problematic. The aims are the characterisation, the quantification and the modelling of transport and interactions of CO₂ in a carbonate unsaturated zone. This was achieved following an experimental approach on a natural pilot site in Saint-Emilion (Gironde, France), where diffusive gas leakage experiments were set up in a carbonate unsaturated zone. Different aspects were investigated during the study: natural pilot site description and instrumentation; the physical and chemical characterisation of carbonate reservoir heterogeneity; the natural functioning of the carbonate unsaturated zone and especially the set-up of a CO₂ concentrations baseline; the characterisation of gas plume extension following induced diffusive leakage in the carbonate unsaturated zone and the study of gas-water-rock interactions during a CO₂ diffusive leakage in a carbonate unsaturated zone through numerical simulations. The results show the importance of the carbonate reservoir heterogeneity characterisation as well as the sampling and analysing methods for the different phases. The baseline set-up is of main interest since it allows discrimination between the induced and the natural CO₂ concentrations variations. The transfer of CO₂ in a carbonate unsaturated zone is varying in function of physical and chemical properties. This transfer is done by diffusion and/or advection. Because the detection of the noble gases allows the prediction of CO₂ plume arrival, the use of tracers in the sequestration site is of main importance. The chemical interactions have to be taken under account in transport models in order to predict delay factors and the impact of a CO₂ leakage in a carbonate unsaturated zone.

Le stockage géologique du carbone est l’une des techniques les plus prometteuses pour réduire le taux de CO₂ dans l’atmosphère. La séquestration géologique possède la capacité et la longévité potentielles pour diminuer les émissions de CO₂ vers l’atmosphère. Dans le cadre d’injections à l’échelle industrielle, les réservoirs carbonatés peuvent faire partie des sites aptes à stocker du CO₂. Toutefois, ces injections à grandes profondeurs sont sujettes à des risques de fuites du piège géologique lui-même ou des infrastructures liés à l’exploitation du site de stockage. Ainsi, il existe principalement deux types de fuite : brutale et diffuse. Dans les deux cas, elles sont susceptibles d’entrainer des risques pour l’environnement et de mettre en danger les populations. Il est ainsi nécessaire de développer des outils capables de prévenir une fuite de CO₂ quel que soit son type. Par ailleurs, il est particulièrement indispensable de comprendre les mécanismes de transport réactif qui rentrent en jeu lors de l’arrivée de cette fuite en contexte de proche surface (zone vadose) et ainsi d’essayer d’étudier comment cette fuite peut s’amortir. Ces travaux de thèse traitent donc de la caractérisation, de la quantification et de la modélisation des processus de transferts et des interactions CO₂-H₂O-CaCO₃ dans la zone vadose en contexte de fuite à partir d’un puits de forage. Cette problématique a été d’abord abordée par une approche expérimentale sur un site pilote à Saint-Emilion. Puis, les interactions CO₂-H₂O-CaCO₃ ont été étudiées au travers d’une approche expérimentale à l’échelle de la carotte en laboratoire. L’approche expérimentale a conduit à la réalisation de deux fuites dans la zone vadose du site pilote : une fuite diffuse et une fuite ultra diffuse. Elles furent réalisées dans la continuité des expériences qui avaient déjà eu lieu auparavant. Une comparaison de l’ensemble des fuites a montré la nécessité d’utiliser des gaz nobles comme précurseurs de l’arrivée en surface du CO₂. Selon le type de fuite, l’hélium peut servir de précurseur temporel du CO₂, tandis que le krypton prévient de l’étendue du panache de gaz durant la fuite. Plus la pression d’injection du CO₂ est importante et plus le gaz migre par advection. Par ailleurs, une pression d’injection importante favorise l’existence de passage préférentiel dans la zone vadose. L’utilisation d’isotopes tels que ceux de l’hélium et du carbone permet de mettre en évidence la présence locale de phases aqueuses dans le massif et de déterminer l’origine biologique ou anthropique du CO₂. Les expériences à l’échelle de la carotte permettent d’estimer le pouvoir tampon des calcaires oligocènes en fonction du faciès de la roche. La perméabilité et la porosité de celle-ci conditionnent la dissolution des calcaires. De même, la réactivité des carbonates en contexte de fuite dépend du pH de la phase aqueuse, du débit qui traverse le réseau poreux, de la saturation en eau et des caractéristiques pétro-physiques des carbonates. / Carbon storage is one of the most encouraging methods to decrease CO₂ concentration into the atmosphere. Carbon storage provides the longevity and the capacity needed to decrease CO₂ emissions toward the atmosphere. When dealing with storage on an industrial scale, carbonated reservoirs can be among the most suitable storage sites. However, these high depth injections are subject to leakage risks from the geologic trap itself or from the framework created by the establishment of the site. Two main types of leakage exist: brutal and diffusive leakage. In both cases, they are likely to endanger the environment and the population. Therefore, it is essential to develop tools that are able to anticipate any types of CO₂ leakage. Furthermore, it is also necessary to understand the reactive transport mechanism that take place when the leakage arrives in the shallow subsurface (vadose zone)and to see how the leakage can be buffered. This work deals with the characterization, the quantification and the modelling of transfer processes and CO₂-H₂O-CaCO₃ interactions into the vadose zone in a context of a leakage from a drilling well. This issue was first dealt through field experiment on the site of Saint Emilion. Then, the CO₂-H₂O-CaCO₃ interactions were studied through an experimental approach in laboratory. Two leakage experiments were performed on the site: a diffusive leakage and an ultra-diffusive leakage. They were performed as a sequel of former experiments carried on the pilot site. A comparison of all the leakage experiments revealed the necessity to use noble gases as precursor of the CO₂ arrival at the surface. Depending of the type of the leakage, helium can be a temporal precursor while krypton can anticipate the spread of the CO₂ gas plume. The higher the injection pressure, the more the gas migrates through advective flux. Moreover, a high injection pressure favors the existence of preferential paths in the vadose zone. The use of helium and carbon isotopes makes it possible to reveal the presence of a local aqueous phase within the porous media and to identify the origin of CO₂. The core scale experiments lead to the estimation of the buffering power of Oligocene limestone according to the rock facies. The permeability and the porosity influence the dissolution of the limestone. The reactivity of carbonates during a leakage depends on the pH of the aqueous phase, the flow rate that goes through the porous media, the water saturation and petrophysical characteristics of the carbonates.

Les bassins d’infiltration font partie intégrante des techniques alternatives de gestion des eaux pluviales en milieu urbain. Néanmoins, la potentialité de transfert de polluants vers la nappe est accrue en cas d’écoulements préférentiels dans les sols sous-jacents. Une bonne compréhension du couplage entre processus d’écoulements préférentiels en zone vadose et mécanismes géochimiques est requise. Cette thèse s’inscrit dans le cadre du suivi d’un bassin d’infiltration depuis plusieurs dizaines d’années de fonctionnement. Le site d’étude est situé sur le dépôt fluvioglaciaire hétérogène couvrant une grande partie de la région lyonnaise. Des auscultations sur une fosse sous le bassin (section 13.5m*2.5m) ont mis en évidence une régionalisation particulière de la pollution dans le sol. Cette étude s’appuie sur une étude numérique visant à identifier l’origine de la régionalisation des polluants et à la relier aux écoulements préférentiels résultant des hétérogénéités lithologiques. En amont de l’étude numérique, les lithofaciès sont complètement caractérisés aux regards de leurs propriétés hydrodynamiques, hydrodispersives et géochimiques. La modélisation numérique permet de souligner l’établissement de cheminements préférentiels en lien avec le contraste de propriétés hydrodynamiques, notamment lorsque de faibles débits sont appliqués en surface. Le rôle de chaque lithofaciès et de l’architecture du dépôt (stratification et inclusions) est clairement identifié. Les répercussions de tels écoulements sur les transferts non réactifs sont ensuite investiguées en combinant l’influence des écoulements préférentiels et le fractionnement de l’eau en fractions mobile et immobile résultant de l’hétérogénéité intrinsèque au sein de chaque lithofaciès. Enfin, ces processus physiques sont couplés à la réactivité géochimique pour le cas d’un polluant modèle (le cuivre) en prenant en compte la réactivité différentielle des lithofaciès. Ces résultats permettent de générer un modèle conceptuel d’écoulements préférentiels et de transfert de masse en milieu fortement hétérogène. / An infiltration basin is a stormwater best management practice (BMP) designed to infiltrate runoff volumes in urban areas. Nevertheless, preferential flow paths in the underlying soil may cause rapid migration of pollutants, thus contributing to groundwater contamination. Understanding the coupling between preferential flow processes in the vadose zone and geochemical mechanisms is then required. This thesis is a part of the follow-up of an infiltration basin for several decades of exploitation. The study site was settled over a highly heterogeneous glaciofluvial deposit covering much of the Lyon region. The investigation of an excavated section of the basin (13.5m long and 2.5m deep) pointed out a specific regionalization of pollution in the soil. This research is based on a numerical study to identify the origin of such a pollutant pattern and link this with preferential flow resulting from lithological heterogeneities. Different lithofacies were fully characterized regarding their hydraulic, hydrodispersive and geochemical properties. The numerical study proves that the high contrast in hydraulic properties triggers the establishment of preferential flow (capillary barriers and funneled flow). Preferential flow develops mainly for low initial water contents and low fluxes imposed at surface. The role of each lithofacies and architecture of deposit (stratification and inclusions) is clearly identified. The impact of such flows on non-reactive transfers is then investigated by combining the influence of preferential flow and pore water fractionation info into mobile and immobile fractions, resulting from the intrinsic heterogeneity within each lithofacies. Finally, these physical processes are coupled to the geochemical reactivity for a pollutant model (copper), taking into account the differential reactivity of lithofacies. These results generate a conceptual model of preferential flow and mass transfer in strongly heterogeneous media.

La compréhension des écoulements préférentiels dans la zone non saturée en fond de bassins d’infiltration est une condition préalable à une gestion optimale de ces ouvrages. Le basin d’infiltration étudié repose sur un dépôt fluvioglaciaire hétérogène constitué d’un lithofaciès majoritaire contenant des lentilles de sable et de gravier. Une étude géophysique (radar) a permis de caractériser cette hétérogénéité sédimentologique le long de radar-grammes (sections de quelques mètres de profondeurs et de plusieurs dizaines de mètre de long). Les écoulements ont été modélisés pour toutes ces sections ainsi qu’une section fictive remplie uniquement du lithofaciès majoritaire (témoin) et pour divers évènements pluvieux type incluant une chronique météorologique complète. Les résultats montrent clairement que l'hétérogénéité lithologique a un effet significatif sur l’écoulement avec la mise en place d’écoulements hétérogènes du fait de phénomènes de barrière capillaire au niveau des lentilles de sable et de gravier en non-saturé. Cet effet se répercute sur les chemins de l’eau dans le dépôt et sur les flux en surface (infiltration et ruissellement). Les résultats montrent aussi que cet effet dépend fortement de la section considérée. La variabilité spatiale de l’hétérogénéité lithologique induit donc une variabilité spatiale des écoulements préférentiels. Certains paramètres architecturaux et géométriques relatifs à l’hétérogénéité lithologique ont été proposés pour expliquer l’effet de l’hétérogénéité lithologique sur les écoulements préférentiels, de sorte à proposer une première tentative de modélisation des écoulements à méso-échelle. / Understanding preferential flows in the vadose zone is a prerequisite for a proper management of infiltration basins. The studied infiltration basin lies on a strongly heterogeneous glaciofluvial deposit made of a predominant lithofacies that embeds lenses of sand and gravel. Geophysics (GPR) were used to characterize the lithological heterogeneity along radar-grammes (sections few meters deep and dozens of meters long). Flow was modelled for all these sections and compared to a synthetic section filled only with the predominant lithofacies (control). Numerical results show that, under unsaturated conditions, the lithological heterogeneity has a strong impact on flow with the establishment of preferential flows induced by capillary-barrier effects in the vicinity of lenses of sand and gravel. These effects have also consequences on water fluxes at surface (infiltration and runoff) and depend strongly on the section considered for numerical modelling. It can thus be concluded that the spatial variability of lithological heterogeneity induces a spatial variability of preferential flow. Geometrical and architectural parameters were defined to quantity lithological heterogeneity and link it to preferential flows, with the aim to propose a first step towards modelling flow in the basin at the meso-scale. / O entendimento do fluxo na zona não saturada em formações heterogêneas é um prérequisito para uma gestão adequada de poluentes em bacias de infiltração inseridas nessas formações. Este estudo aborda o efeito da heterogeneidade litológica nos processo de escoamento numa bacia de infiltração com geologia inserida num depósito fluvioglacial. A bacia de infiltração é a Bacia Django Reinhardt, existente em Lyon. Esta bacia foi objeto de vários estudos que demonstraram algumas avaliações do impacto da heterogeneidade do solo. Os objetivos desses trabalhos foram todos baseados apenas no estudo sedimentológico de uma trincheira. Nenhum estudo abordou a potencial variabilidade espacial dos atributos do solo. Neste trabalho, foi modelizado o fluxo na zona vadosa através de três casos: drenagem interna, infiltração de água durante um evento de precipitação pluviométrica constante e solicitação hidrológica a partir de um fluxo variável escoamento superficial. Estes cálculos foram realizados para várias seções, que foram previamente caracterizadas quanto a sua heterogeneidade estrutural usando a técnica não destrutiva do GPR com uma frequência de 200MHz e 100MHz. O comportamento hidráulico dessas seções foi comparado com a de uma seção homogênea. Foi comparada o efeito da resolução da imagem nas estimativas das variáveis de fluxo. Alguns indicadores arquiteturais, geométricos e texturais foram propostos para esclarecer e avaliar a variabilidade espacial do efeito da heterogeneidade do solo sobre o fluxo em condições não saturadas. Os resultados mostram claramente que a heterogeneidade gera impactos nas variáveis de fluxo em condições não saturadas. Os impactos dependem da seção considerada.Os resultados mostram também que: a) a heterogeneidade estrutural do subsolo na meso-escala apresenta um efeito sobre a geração de escoamento preferencial e isso independe do tipo de solicitação aplicada na seção; b) O efeito da heterogeneidade se mostrou independentende das condições de contorno aplicadas na superfície; c) Nos períodos de ocorrência de chuvas intensas quem controla o comportamento das variáveis de fluxo é a intensidade de precipitação. No entanto, nos períodos em que ocorre ausência de precipitação, as variáveis de fluxo são controladas pela heterogeneidade da zona não saturada; d) A utilização de diferentes frequências influencia na estimativa das variáveis de fluxo. Isso é decorrente do impacto na quantificação da heterogeneidade decorrente do efeito da resolução de cada antena; e) Os indicadores não apresentaram relação direta com as variáveis de fluxo estudadas. Isso demonstra que o efeito da heterogeneidade não pode ser explicado por um único fator, mas por um conjunto de fatores.