Caractérisation de la dynamique transitoire des écoulements d'un site de filtration sur berge localisé dans un système géologique complexe

La filtration sur berge (FSB) est un terme qui décrit des installations de puits qui exploitent un mélange d’eau souterraine et d’eau de surface en raison de leur position à côté d’un lac ou d’une rivière. Ces systèmes ont de nombreux avantages, notamment dans un contexte où les eaux de surface sont de plus en plus contaminées, puisqu’ils permettent d’accéder à des quantités d’eaux de surface filtrées naturellement et diluées avec l’eau souterraine de meilleure qualité. Les problématiques de gestion de la qualité sont spécifiques pour ces systèmes à cause du caractère hautement dynamique et transitoire des écoulements, hérités des séquences de pompage et de la variation naturelle dans l’eau de surface. Les éléments essentiels pour comprendre l'évolution de la qualité de l'eau sont les taux de mélange entre la rivière et les eaux souterraines, qui influent sur la dilution, ainsi que les temps de séjour, qui affectent la capacité d'atténuation. Classiquement, les taux de mélange sont calculés en considérant 2 pôles d’eau (souterrain et surface) par des modèles hydrodynamiques ou des approches géochimiques et isotopiques. Aucune étude ne considère l’impact de plusieurs types d’eau souterraine sur les taux de mélange et les effets de dilution. Pourtant, les lieux de prédilection de l’installation des sites de FSB sont les aquifères alluviaux, par définition hétérogènes et fréquemment exutoires d’aquifères régionaux de plus grande envergure. En outre, si certaines études s’intéressent à l’impact de la variation temporelle du signal de l’eau de surface sur les taux de mélange des eaux pompées, la majorité la considère comme fixe. D’autre part, l’infiltration induite de l’eau de surface a pour conséquence de perturber les équilibres physico-chimiques et biologiques dans la zone d’infiltration entre l'eau de surface et les puits de pompage. La perturbation des équilibres est à double tranchant, elles peuvent être à l’origine de l’amélioration de l’atténuation des contaminants comme à leur relargage. La forte dynamique transitoire des écoulements, liée à des schémas de pompage qui s’adaptent à la demande en eau, est susceptible de rompre la connexion hydraulique avec la rivière. Les bénéfices occasionnés par cette connexion sont alors perdus. L’objectif général de cette thèse est de proposer une approche à même de rendre compte de la dynamique transitoire des écoulements d’eau dans un site de FSB dont les sources d’eau sont multiples. Pour ce faire, une méthode d’investigation combinant des séries temporelles géochimiques et isotopiques à différentes échelles de temps et réparties dans l’espace a été appliquée sur un site de prélèvement multipuits à Lachute (Québec, Canada). Les objectifs spécifiques se déclinent en quatre : 1) Caractériser les structures séparant les aquifères régionaux de la nappe alluviale par une analyse de l’histoire géologique et par des méthodes géophysiques; 2) Qualifier et quantifier l’impact des écoulements souterrains régionaux sur la dynamique de recharge d’une nappe alluviale utilisée pour la FSB par l’élaboration d’un modèle conceptuel à partir d’outils géochimiques et isotopiques; 3) Déterminer l’impact du caractère transitoire et multisource des écoulements d’un site de FSB sur la quantification des taux de mélange des eaux pompées avec des équations de bilan de masse; et 4) Déterminer l’impact d’une stratégie de pompage sur la représentativité d’un échantillon dans la zone d’infiltration et dans le puits de pompage. Les deux premiers objectifs spécifiques sont nécessaires pour déterminer le modèle conceptuel des écoulements du site de FSB étudié. L’hypothèse avancée est que la nappe alluviale ne peut être représentée par un pôle d’eau unique dans un modèle de mélange en FSB en étant l’exutoire des aquifères adjacents. Il apparaît primordial de caractériser les masses d’eau susceptibles de s’écouler dans la nappe alluviale et d’évaluer leur qualité respective. L’approche scientifique est de définir les écoulements d’eau souterraine depuis les aquifères régionaux vers la nappe alluviale à partir des données dont les conditions se rapprochent le plus d’un fonctionnement en régime « naturel » pré-opératoire où il n’y avait pas encore d’infiltration induite de la rivière. Les travaux ont démontré qu’il existe 4 pôles d’eau souterraine affectant la variabilité spatiale de la nappe alluviale : 3 naturels et 1 anthropique, lié à la présence d’un ancien site de dépôt à neige. Ils se distinguent par la combinaison des faciès chimiques et de la signature isotopique de la molécule d’eau. Ces travaux ont également démontré que la nappe alluviale n’est pas à l’origine de la minéralisation et de la signature isotopique de l’eau souterraine qui y circule. Elles sont héritées des eaux qui rechargent cette nappe. Ce constat permet d’utiliser certains ions majeurs pour distinguer les zones de mélange dans la nappe alluviale et de calculer les taux de mélanges. Parmi les 4 sources d’eau souterraines affectant la nappe alluviale, le pôle anthropique est celui qui a une qualité médiocre avec des concentrations en Ca et Mg élevées. L’utilisation de la base de données de la nappe alluviale qui n’est pas sous l’influence de l’infiltration de la rivière a permis de localiser la source anthropique à l’extrémité ouest du champ captant. La partie ouest du champ captant est également sous l’influence d’un pôle naturel. Aucun taux de mélange n’a pu être calculé entre ces deux pôles dans la partie ouest du champ captant par manque de données. Cependant, un mélange binaire a permis de calculer, par des systèmes d’équations de bilan de masse, les taux de mélange dans la partie est du champ captant. Elle est constituée à 50 – 40 % des eaux souterraines provenant de l’amont et à 50 - 60% des eaux souterraines provenant de l’aquifère sous-jacent. Devant ce constat, considérer un seul pôle d’eau souterraine pour représenter la nappe alluviale dans les calculs de taux de mélange en FSB ne paraît pas approprié. Des modèles de mélange ont été définis pour discuter les taux de mélange dans les puits d’exploitation calculés par des systèmes d’équations de bilan de masse. Le pôle d’eau de surface (qui varie dans le temps) a dû être considéré en plus des pôles d’eau souterraine considérés précédemment. Dans la partie est du champ captant, la même tendance des taux de mélange provenant des eaux souterraines s’observe pour les 2 systèmes d’équations : % nappe sous-jacente > % nappe amont. La part d’eau de surface est celle qui est la plus importante entre les trois. Il y a une variation temporelle des taux de mélange pour la nappe sous-jacente et pour l’eau de surface infiltrée, la part de la nappe amont est plus stable dans le temps. Ainsi, les pompages ont tendance à solliciter prioritairement les eaux souterraines de la nappe sous-jacente, avant celle de la nappe amont. Ce qui est un atout pour la qualité, car la nappe sous-jacente est captive et plus protégée des contaminations s’infiltrant par le toit de la nappe. Dans la partie ouest du champ captant, les eaux provenant du pôle anthropique représentent 1 % à 17 % du mélange dans les 3 puits à l’extrémité ouest du champ captant. C’est le taux le plus bas devant le pôle d’eau de surface et le pôle de la nappe amont, pourtant son impact est important sur la qualité dans les puits. Les puits ayant le plus haut taux provenant du pôle anthropique sont aussi ceux qui ont la qualité la plus faible. Devant le défaut de qualité de la source anthropique et l’apparente hétérogénéité des écoulements dans la nappe alluviale, il apparaît important d’évaluer l’impact sur les taux de mélange si une seule source d’eau souterraine était considérée dans les calculs. L’étude de sensibilité appliquée au système d’équations binaire avec la conductivité électrique a montré que les paramètres temps de séjour et pôle d’eau souterraine unique ont un impact important sur les taux de mélange. Une erreur d’un pas de temps de mesure, ici le mois, peut entraîner des différences jusqu’à 40 %. Une mauvaise approximation de la valeur représentant le pôle d’eau souterraine a un impact variable en fonction de la valeur du traceur dans chacun des pôles. Dans le cas d’une contamination anthropique au sel, les erreurs peuvent être jusqu’à 50 %. Ces erreurs sont très problématiques pour estimer la dilution susceptible d’atténuer la contamination. Le pôle d’eau de surface a été estimé à chaque instant à partir de temps de séjour moyen défini en comparant graphiquement les pics de variations dans les chroniques de la rivière et des puits d’exploitation. Les temps de séjour moyens dans chaque puits sont compris entre 2 mois et 6 mois. Une variation temporelle des temps de séjour moyen est suspectée malgré les chroniques discontinues. Ce constat était attendu dans ce contexte hautement transitoire qu’est la FSB où les écoulements sont bouleversés au rythme des séquences de pompage. Une hypothèse de travail est que les variations de temps de séjour sont liées en partie à une rupture de la connexion hydraulique entre la rivière et les puits de pompage. Cette rupture en plus d’entraîner des conséquences sur les temps de séjour est susceptible d’impacter les sources d’eau sollicitées. Des conditions transitoires de pompage ont été recréées en interrompant durant une semaine l’activité de tous les puits de la partie ouest du champ captant, avant d’en redémarrer un. L’objectif étant de rompre la connexion hydraulique entre la rivière et le puits de pompage. Durant les deux semaines suivantes de pompage, l’évolution des charges hydrauliques, de la géochimie et de la signature isotopique des eaux pompées et autour a été suivie. Les eaux sont constituées à 54-67% d’eau souterraine et 46-54% d’eau de surface anciennement infiltrée. Ces tendances évoluent drastiquement la première journée de l’ordre 60 % de l’augmentation de la CE. Puis, elle évolue plus lentement. Le dernier jour de pompage, les eaux prélevées par le puits d’exploitation sont un mélange avec 97% d’eaux souterraines amont et 3% d’eaux de recharge. Cette évolution valide que la rupture de la connexion hydraulique entraîne un bouleversement des écoulements dans la zone d’infiltration et que le rétablissement de la connexion hydraulique n’est pas instantané. Les conséquences techniques directes sont 1) l’augmentation du temps de séjour dans la zone d’infiltration et 2) la définition de la représentativité d’un échantillon doit s’effectuer en tenant compte de la dynamique de pompage. La rupture de la connexion hydraulique entraîne une perte des bénéfices de qualité liée. Cependant, la connaissance des conditions de pompage nécessaires pour rompre la connexion hydraulique est un atout dans la gestion des schémas de pompage. Par exemple, il serait possible d’éviter qu’une contamination ponctuelle dans l’eau de surface s’infiltre dans la zone d’infiltration en rompant la connexion hydraulique. Pour conclure, la représentativité des échantillons est une notion fondamentale dans la compréhension de la dynamique transitoire d’un site de FSB. Il est primordial d’obtenir des données réparties spatialement pour caractériser les eaux souterraines et des chroniques continues des traceurs dans les eaux de surface. D’une part, il est recommandé d’évaluer toutes les sources d’eau souterraine alimentant la nappe alluviale et leur qualité respective. D'autre part, il est recommandé d’évaluer les conditions de pompage lors d’un échantillonnage pour évaluer les sources d’eau impliquées. Pour ce faire, il est nécessaire de caractériser les conditions de pompage qui coupe la connexion hydraulique entre la rivière et les puits de pompage, et le temps nécessaire pour que la connexion se rétablisse. La connaissance de ces éléments permettra d’optimiser les stratégies de pompage pour limiter les variations dans la qualité de l’eau pompée.