Dans un champ captant, comme celui qui alimente l'agglomération lyonnaise, le fonctionnement de chaque " objet " (bassins d'infiltration, puits, rivières, nappe, zone non saturée) et leurs interactions sont complexes et mal connus. Dans un premier temps, une série d'essais d'infiltration à différentes échelles dans un bassin artificiel de réalimentation a permis de mieux appréhender le fonctionnement de cet ouvrage et de chiffrer les paramètres hydrodynamiques de la zone non saturée. Les résultats des interprétations par méthodes numériques fondées sur l'équation de Richards ont montré que les flux infiltrés dans les bassins dépendent principalement de la conductivité hydraulique à saturation d'une couche située directement sous le sable calibré couvrant le fond du bassin, identifiée comme étant moins perméable que la nappe. Cette couche conditionne l'existence d'une zone non saturée. La réalisation d'essais de pompage dans l'aquifère sur les forages d'exploitation et sur un dispositif spécialement mis en place durant ce travail a permis de déterminer les paramètres hydrodynamiques de la nappe. Une analyse des observations et une modélisation conceptuelle en 2D, puis en 3D ont permis d'identifier les mécanismes prépondérants (stratifications, apports et prélèvements) et de simuler correctement à la fois les flux infiltrés dans un des bassins d'infiltration et la remontée de la nappe. A l'échelle d'un bassin, les flux infiltrés sont variables dans le temps, ils dépendent de l'état de colmatage de la surface d'infiltration mais également de la température de l'eau infiltrée et de l'état hydrique initial du sol sous le bassin. Les analyses de sensibilité réalisées avec les modèles mis en place indiquent que la conductivité hydraulique à saturation de l'aquifère, mais également la proximité des conditions aux limites imposées dans la nappe (les rivières et les puits de pompage) influencent de manière prépondérante la remontée de la nappe. Une modélisation 3D d'un autre secteur du champ captant comprenant deux bassins d'infiltration, deux bras de rivière ainsi que des puits de pompage a été réalisée. La condition à la limite imposée sur les rivières est du troisième type en accord avec l'observation d'un sous-écoulement en nappe. Les échanges nappe/rivières sont calés sur des observations à partir d'une chronique de propagation d'une onde de crue dans la nappe. Des piézomètres en flûte de pan, spécifiquement implantés à proximité d'un bassin, ont permis d'observer des différences de charge hydraulique fortes à différents niveaux de l'aquifère lorsque le bassin d'infiltration est en eau. La modélisation 3D est conforme à ces observations. Elle a confirmé l'importance du rôle d'une hétérogénéité de type argilo-sableuse (de conductivité hydraulique à saturation inférieure aux autres lithologies présentes dans l'aquifère) dans les écoulements (direction et flux). Le modèle développé représente correctement les flux infiltrés via les bassins ainsi que les fluctuations de la piézométrie de la nappe. Il permet de vérifier l'inversion des écoulements par rapport aux infiltrations de la rivière, d'identifier les puits alimentés par les bassins d'infiltration et également de mettre en évidence les flux de nappe passant sous la rivière.

Dans un champ captant, comme celui qui alimente l'agglomération lyonnaise, le fonctionnement de chaque « objet » (bassins d'infiltration, puits, rivières, nappe, zone non saturée) et leurs interactions sont complexes et mal connus. Dans un premier temps, une série d'essais d'infiltration à différentes échelles dans un bassin artificiel de réalimentation a permis de mieux appréhender le fonctionnement de cet ouvrage et de chiffrer les paramètres hydrodynamiques de la zone non saturée. Les résultats des interprétations par méthodes numériques fondées sur l'équation de Richards ont montré que les flux infiltrés dans les bassins dépendent principalement de la conductivité hydraulique à saturation d'une couche située directement sous le sable calibré couvrant le fond du bassin, identifiée comme étant moins perméable que la nappe. Cette couche conditionne l'existence d'une zone non saturée. La réalisation d'essais de pompage dans l'aquifère sur les forages d'exploitation et sur un dispositif spécialement mis en place durant ce travail a permis de déterminer les paramètres hydrodynamiques de la nappe. Une analyse des observations et une modélisation conceptuelle en 2D, puis en 3D ont permis d'identifier les mécanismes prépondérants (stratifications, apports et prélèvements) et de simuler correctement à la fois les flux infiltrés dans un des bassins d'infiltration et la remontée de la nappe. A l'échelle d'un bassin, les flux infiltrés sont variables dans le temps, ils dépendent de l'état de colmatage de la surface d'infiltration mais également de la température de l'eau infiltrée et de l'état hydrique initial du sol sous le bassin. Les analyses de sensibilité réalisées avec les modèles mis en place indiquent que la conductivité hydraulique à saturation de l'aquifère, mais également la proximité des conditions aux limites imposées dans la nappe (les rivières et les puits de pompage) influencent de manière prépondérante la remontée de la nappe. Une modélisation 3D d'un autre secteur du champ captant comprenant deux bassins d'infiltration, deux bras de rivière ainsi que des puits de pompage a été réalisée. La condition à la limite imposée sur les rivières est du troisième type en accord avec l'observation d'un sous-écoulement en nappe. Les échanges nappe/rivières sont calés sur des observations à partir d'une chronique de propagation d'une onde de crue dans la nappe. Des piézomètres en flûte de pan, spécifiquement implantés à proximité d'un bassin, ont permis d'observer des différences de charge hydraulique fortes à différents niveaux de l'aquifère lorsque le bassin d'infiltration est en eau. La modélisation 3D est conforme à ces observations. Elle a confirmé l'importance du rôle d'une hétérogénéité de type argilo-sableuse (de conductivité hydraulique à saturation inférieure aux autres lithologies présentes dans l'aquifère) dans les écoulements (direction et flux). Le modèle développé représente correctement les flux infiltrés via les bassins ainsi que les fluctuations de la piézométrie de la nappe. Il permet de vérifier l'inversion des écoulements par rapport aux infiltrations de la rivière, d'identifier les puits alimentés par les bassins d'infiltration et également de mettre en évidence les flux de nappe passant sous la rivière. / In a well field of the Lyon metropolitan area, designed for drinking water supply, behaviour of each object (infiltration basins, wells, rivers, aquifer, and unsaturated zone) and their interactions are complex and not well-known. As a first step, infiltration tests at different spatial scales in one artificial basin were performed to better understand the basin operation and to estimate the hydrodynamic parameters of the unsaturated zone. Results of interpretation, using numerical methods based on Richards equation, reveal that infiltrated basin fluxes mainly depend on saturated hydraulic conductivity of a layer located just below the calibrated sand layer that cover the basin bottom. Indeed this layer has been estimated to be less permeable than the aquifer, which allows the existence of the unsaturated zone below. Pumping tests in the groundwater have been performed using production wells and a well specially implemented during this thesis work in order to estimate aquifer hydrodynamic parameters. Observations analysis and a conceptual modelling, in 2D and then in 3D, lead to a better understanding of the controlling mechanisms (stratification, input and output) and to simulate both basin infiltration rates and water table rise. Considering the whole basin scale, input fluxes are transient, related to the clogging statement of the infiltration area but also to the temperature of inflow water and the initial statement of the soil just below the basin. Sensibility analyses using the models highlight that the amount of the water table rise is mainly influenced by the aquifer saturated hydraulic conductivity and also by the location of imposed boundaries in the aquifer (rivers and pumping wells). The model properly accounts basin inflow fluxes and water table fluctuations. The model is able to verify if flows are reversed in relation to river exchanges, if wells are fed by infiltration basins and it highlights aquifer flows below the river. A 3D modelling has been realised in another area of the well field, comprising two infiltration basins, two river arms and pumping wells. In agreement with underflow in the aquifer, rivers are imposed in the model as third kind boundary conditions. Aquifer and river exchanges are calibrated with observed data of one aquifer flood-wave propagation. Significant differences of hydraulic heads have been observed at different depths of the aquifer using panpipes piezometers, specifically implemented, close to one infiltration basin. Theses differences are closely related to basin operation. These observations are properly calculated by the 3D model. Using the model, the effect of one sandy-clay heterogeneous layer (whose saturated hydraulic conductivity is lower than the ones of other aquifer lithologies) on aquifer flows (direction and flux) is notable. The model properly accounts basin inflow fluxes and water table fluctuations. The model is able to verify if flows are reversed in relation to river exchanges, if wells are fed by infiltration basins and it highlights aquifer flow below the river.