Les déversements des eaux des réseaux unitaires dans le milieu naturel ont pour conséquences en période d'étiage des perturbations du régime hydraulique et une forte dégradation de la qualité des eaux du milieu naturel. Cette thèse étudie les relations entre les débits et charges déversées et la qualité des eaux du milieu naturel et a pour but de définir le modèle le mieux adapté pour l'étude de ces phénomènes. La zone d'étude est limitée à la Seine entre les barrages de Port-à-l'Anglais à l'amont de Paris et de Suresnes à l'aval de Paris soit un tronçon de 25 km. Les déversoirs d'orage considérés sont ceux de la ville de Paris intra-muros.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal.

En France les filtres plantés à écoulement vertical sont utilisés pour le traitement des surverses de déversoir d’orage (CSO en anglais). Ils sont efficaces pour le traitement des polluants ainsi que pour assurer un rôle de tampon hydraulique permettant de protéger les milieux hydriques superficiels. Ils ont la particularité de recevoir des effluents bruts dont la fréquence, l’intensité, la durée et les concentrations sont stochastiques. En conséquence de quoi l’optimisation du dimensionnement est particulièrement délicate. Dans ce travail, le dimensionnement est abordé suivant deux approches. D’une part en définissant les processus clefs pour une optimisation du dimensionnement et, d’autre part, en construisant un modèle dynamique simplifié (Orage)permettant de prendre en compte, sur le long terme, le caractère stochastique des évènements et les contraintes locales. Le cœur du modèle a été calibré et validé sur une large gamme d’évènements mesurées sur un site en taille réelle (site de Marcy l’Etoile). La qualité de la calibration a été validé aussi bien visuellement que statistiquement et la robustesse du modèle étudiée par une analyse de sensibilité (méthode de Morris). L’optimisation automatique du dimensionnement a nécessité la réalisation d’une boucle itérative également paramétrée et testée dans le cadre du site de Marcy l’Etoile et pour des cas théoriques. Les dimensionnements proposés ont été estimés réalistes. La définition des processus clefs du système de traitement, et le développement du modèle simplifié, a également été possible grâce à la simulation d’expérimentations sur colonnes par un modèle mécaniste (HYDRUS / CW2D). La paramétrisation des processus du modèle simplifié a été réalisée sur la base d’expérimentations détaillées sur le site de Marcy l’Etoile. Le site a été suivi pendant trois ans aussi bien par des prélèvements automatiques, des traçages, que des mesures en lignes (hydraulique, azote, teneur en eau). Parmi les processus calibrés, nous retiendrons i) l’importance du niveau de saturation de la base du filtre au début d’un évènement vis-à-vis des court circuits hydrauliques, ii) la définition des capacités d’adsorption de l’azote ammoniacale de différents matériaux (pouzzolane, mélange de sable et de zéolite) et iii), la définition des cinétiques de nitrification suivant la température et quantité d’azote ammoniacal adsorbé. Si les performances de l’ouvrage sont hautes pour les paramètres majeurs (MES, DCO, N-NH4), l’impact du dimensionnement, de la gestion et des facteurs environnementaux sur les performances et l’accumulation de boue a été étudié. De plus, quelques suivis ont été réalisés sur les métaux et les HAP. Ces mesures ont permis non seulement de calibrer le logiciel Orage mais aussi de proposer des modifications pour l’amélioration du dimensionnement. / Constructed wetlands for combined sewer overflow treatment (CSO CWs) are vertical flow filters in France. They effectively remove pollutants and mitigate hydraulic peaks, which protects natural waters. They receive unsettled flows with stochastic return periods, volumes and concentrations, making design optimization difficult. In the presented work, design is targeted from two sides. First, the scientific basis is laid down for the design-support software Orage. The tool is model-based and as such, it considers the stochasticity of CSOs site-specifically. Its core model has been calibrated and then verified by a range of tests. Simulations fit to measured data from a full-scale CSO site (Marcy l’Etoile). The goodness of fit is evaluated visually and statistically. Model robustness has been confirmed by a sensitivity analysis (Morris method). After this, the iterative shell has been parameterized and tested (the software element doing automatic optimization). The design proposals of the tool have been found realistic. The other approach to target optimization has been a detailed field research at Marcy l’Etoile. The site at Marcy l’Etoile has been monitored for three years, using automatic samplers and online probes. The filling of the porous media at event start proves and explains shortcutting behaviour. Adsorption capacities have been quantified for pozzolana and a sand and zeolite mixture. A temperature-dependent equation is calibrated to calculate nitrification rate. System efficiency is high for target pollutants (TSS, COD, NH4-N). Design, operational and environmental factors have been analysed to seek potential effects on removal performances and sludge accumulation. Additionally, PAHs and metals are indicated for a few selected events. The field results were essential to calibrate Orage and to see options for design improvements. The understanding of CSO CWs and the development of Orage was promoted also by simulating lab-scale columns using the process-based model package HYDRUS / CW2D.

Mémoire numérisé par la Division de la gestion de documents et des archives de l'Université de Montréal

La modélisation à l’aide de l’approche RANS (Reynolds Averaged Navier Stokes) a été menée en trois dimensions, en considérant les régimes permanent et non permanent, dans le but de simuler les écoulements au niveau des jonctions, bifurcations et déversoirs d’orage. A travers ces trois exemples d’étude, plusieurs stratégies de modélisation (différentes combinaison lois de paroi/modèles de turbulence, différents algorithmes de couplage pression/vitesse, différents schémas de discrétisation spatiale, différentes conditions aux limites, différents types et tailles de maille, etc.) ont été testées et évaluées à l’aide de plusieurs indicateurs de performance (de type RMS –Root Mean Square) en s’appuyant sur les données de vitesses (vitesses moyennes in situ et champs de vitesses obtenus par PIV en laboratoire), hauteurs d’eau, débits (répartition de débits en bifurcation en laboratoire ou débits déversés in situ). Les résultats obtenus sont transposables aux autres cas de singularités et ouvrages spéciaux rencontrés en réseau d’assainissement et montrent que : i) les résultats des simulations 3D sont sensibles à la rugosité, aux conditions limites de hauteur et de vitesse ; ii) dans les trois cas d’étude, les schémas de discrétisation du second ordre et l’algorithme de couplage pression/vitesse PISO sont appropriés ; ii) la loi de paroi scalable couplée aux modèles de turbulence de type k-ε pour le cas des jonctions (avec un débit latéral inférieur ou égal au débit principal) ou des déversoirs semblent convenir, tandis que le modèle de turbulence RSM associé à la loi de paroi enhanced ou scalable permet de mieux représenter les écoulements à travers les bifurcations ou au niveau des jonctions lorsque le débit latéral est dominant. Sur la base de ces résultats, une méthodologie de modélisation plus générale définie en six étapes et fondée sur le guide proposée par Jakeman et al. (2006) a été mise au point. La méthodologie ainsi définie a été utilisée pour i) améliorer l’instrumentation du déversoir OTHU (Observatoire de Terrain en Hydrologie Urbaine) situé à Ecully, à partir de la simulation de sa courbe de fonctionnement et en fournissant les incertitudes sur les débits déversés obtenus ; ii) simuler l’implémentation des capteurs débitmétriques à l’aval d’une jonction. Elle a permis de concevoir et de dimensionner un nouveau dispositif de maîtrise des flux d’eau et de polluants déversés (technologie DSM – dispositif de surveillance et de maîtrise des flux déversés au milieu naturel). Ce dispositif a fait l’objet d’un dépôt de brevet international. Enfin, la mise en œuvre de cette méthodologie a été à l’origine de la découverte d’une nouvelle structure d’écoulement dans la branche latérale d’une bifurcation à 90°. L’analyse des résultats des simulations des écoulements mettant en évidence cette nouvelle structure a montré qu’il était possible de prédire l’apparition de cette dernière à partir du rapport d’aspect et du nombre de Froude. / The understanding of sewer flows behaviour is a key component for better urban drainage monitoring and management. However, these flows are conveyed across singularities (such as bends, drops, deviations, etc.) and special structures (combined sewer overflows –CSOs–, channels junction, dividing flow structures, etc.). These singularities and specific structures exhibit complex geometries, leading to open channel turbulent, three-dimensional and multiphase (pollutants and storm and sewer waters) flows. Using three-dimensional CFD (Computational Fluid Dynamics) platform allows a better understanding of mechanisms of contaminants transport through these structures and singularities, leading to a better sewer monitoring. In this thesis, 3D-RANS (Reynolds Averaged Navier Stokes) modelling approach under steady-state conditions is used in order to study flows within CSOs, junctions and bifurcations. Through these three structures, several modelling strategies (wall law/turbulence model combination, velocity/pressure coupling algorithm, spatial discretisation schemes, boundary conditions, computational mesh –shape and size of cells–, etc.) are tested and evaluated thanks to performance indicators (such as RMS –Root Mean Square– indicators) based on velocities (in situ mean velocities and PIV velocity fields obtained in laboratory), water depths and discharge (discharge repartition for bifurcation in laboratory or in situ overflow discharge, for CSOs). Results deriving from these tests are transposable to other singularities or special structures encountered in sewer network and suggest that: i) simulated CFD results are sensitive to the roughness coefficient; ii) for the three studied structures, second-order discretisation schemes and SIMPLE or PISO velocity/pressure coupling algorithm are appropriate; iii) scalable wall function associated to the group of k-ε turbulence model for junction flows (with a lateral inflow lower or equal to the main inflow) or for CSOs is appropriate, whereas RSM turbulence model associated to enhanced wall function allows a better representation of bifurcation flows or junctions flows when the lateral inflow is greater than the main inflow. Based on these results and on Jakeman et al. (2006) guideline, a six steps-methodology focused on the using of RANS approach modelling has been proposed. This six steps-methodology is used in order to i) enhance the monitoring of an OTHU (Observatoire de Terrain en Hydrologie Urbaine) CSO located at Ecully accounting for uncertainties on overflow discharge values ;ii) simulate the performance of flowmeters downstream of a junction, defining the best location for these sensors. This methodology is used to design the new overflow discharge measurement device. This device is an international patent. Finally, the application of the methodology led to point out a new flow structure, occurring in the downstream lateral branch of a 90° bifurcation.

Le travail présenté dans ce mémoire s'appuie sur quatre études de cas de récentes catastrophes françaises : le Grand-Bornand (1987), Nîmes (1988), Paris (1990) et Vaison-La-Romaine(1992). Après avoir fait le point sur les connaissances actuelles sur les ruissellements en milieu agricole, forestier et urbain (chapitre 1) et avoir analysé les quatre catastrophes citées et des catastrophes prévisibles (chapitre 2), nous abordons les facteurs qui nous semblent aggraver ou limiter le risque pour en arriver à une série de propositions techniques de lutte contre le ruissellement en milieu urbain (chapitre 3). Enfin, un guide méthodologique et des exemples,pour approcher le problème de la gestion de ces crues, sont proposés. Cette thèse, destinée aux aménageurs du territoire, Urbanistes, Architectes, Paysagistes et Ingénieurs de voirie, est une recherche multidisciplinaire : elle sert d'approche complémentaire aux efforts réalisés par les responsables des réseaux d'assainissement pour lutter contre les crues centennales. Le rapport pose les bases techniques élémentaires pour enclencher les moyens de lutte contre le ruissellement, sachant que cette lutte doit s'effectuer aussi bien contre la formation du ruissellement que contre ses effets, et qu'à chaque phase différente du processus de ruissellement (rétention, conduite d'écoulement, protection) correspondent différentes solutions techniques.

Les rejets urbains par temps de pluie constituent un problème majeur pour les gestionnaires de l'eau en région parisienne. Une meilleure caractérisation des matières organiques rejetées par temps de pluie permet de mieux évaluer la désoxygénation du milieu récepteur, la Seine. Un bassin versant expérimental a été choisi en proche banlieue parisienne pour étudier les eaux usées par temps de pluie. Les résultats montrent que les flux sont une fonction du cycle journalier des eaux usées et de l'intensité des pluies. En journée par temps sec, les teneurs en polluants et les débits sont plus élevés que pendant la nuit. Les flux des polluants déversés pendant la journée sont ainsi supérieurs aux flux rejetés la nuit à pluie égale. Dans les réseaux de petite taille et de faibles pentes, un dépôt de temps sec peut être remis en suspension par des pluies intenses et augmenter temporairement la charge polluante. Les matières organiques du réseau en temps sec sont mieux et plus rapidement dégradables qu'en temps de pluie. La biodégradabilité des matières en suspension dans le réseau diminue avec l'intensité des pluies en raison d'un apport de carbone faiblement dégradable par les eaux de ruissellement et les dépôts en réseau. L'effet toxique des rejets urbains de temps de pluie sur le milieu récepteur a été étudié au laboratoire sur une population mixte du phytoplancton de la Seine. Nos études montrent que les rejets peuvent diminuer la production primaire de la Seine de l'ordre de 60% après une dilution au cinquième, cette inhibition est principalement due aux micropolluants organiques et partiellement aux métaux lourds, dont le zinc. Les effets varient fortement avec la composition des rejets, la toxicité des métaux en particulier dépendant de leur état de complexation. Le comportement de la matière organique et des bactéries hétérotrophes rejetées par temps de pluie dans le milieu récepteur a été étudié en Seine après d'importants rejets du principal déversoir de la région parisienne. Lorsque la dilution du rejet est inférieure à 10%, les effets sur le milieu sont peu sensibles. Les masses d'eau usée ont pu être suivies in situ à l'aide de la conductivité et des concentrations d'oxygène. Une relation entre la quantité de carbone organique dissous biodégradable rejetée et le déficit d'02 a pu être établie. La matière organique fixée aux matières en suspension sédimente très rapidement et ne participe pas directement à la désoxygénation. La désoxygénation in situ est principalement due aux bactéries allochtones, apportées par le rejet, de taille supérieure à 1µm. Leur croissance est plus élevée que celle des bactéries fluviales, mais leur vitesse de disparition l'est aussi. Les polluants toxiques du rejet sont susceptibles de diminuer la production du phytoplancton et d'aggraver ainsi les déficits de près d'un milligramme d'oxygène voire plus en période de forte production. Ces données permettent des simulations détaillées de l'évolution de l'oxygène dissous en Seine suite aux principaux rejets de l'agglomération parisienne.

Cette thèse comporte deux thèmes principaux: l'un concerne l'accélération non adiabatique des ions dans la couche de plasma proche de la Terre, l'autre les pulsations géomagnétiques dues à la décroissance de la pression dynamique du vent solaire. L'accélération nonadiabatique des ions de la couche de plasma est importante pour comprendre la formation du courant annulaire et les injections énergétiques des sous-orages. Dans la première partie de cette thèse, nous présentons des études de cas d'accélération non adiabatique des ions de la couche de plasma observés par Cluster et TC-1 Double Star dans la queue magnétique proche de la Terre (par exemple à (X,Y, Z)=(-7.7, 4.6, 3.0) RT lors de l'évènement du 30 octobre 2006), beaucoup plus près de la Terre que ceux qui avaient été précédemment observés. Nous trouvons que les variations du flux d'énergie des ions, qui sont caractérisées par une décroissance entre 10 eV et 20 keV et une augmentation entre 28 keV et 70 keV, sont causées par l'accélération non adiabatique des ions associée étroitement avec les fortes fluctuations du champ électromagnétique autour de la gyrofréquence des ions H+. Nous trouvons aussi que les ions après l'accélération non adiabatique sont groupés en phase de gyration et c'est la première fois que ceci est trouvé dans la couche de plasma alors que cet effet a été observé dans le vent solaire dans les années 1980. Nous interprétons les variations de flux d'énergie des ions et les groupements en phase de gyration en utilisant un modèle non adiabatique. Les résultats analytiques et les spectres simulés sont en bon accord avec les observations. Cette analyse suggère que cette acceleration nonadiabatique associée aves les fluctuations du cham magnétique est un mécanisme efficace pour l'accélération des ions dans la couche de plasma proche de la Terre. Les structures de flux d'énergie présentées peuvent être uilisées comme un proxy pour identifier ce processus dynamique. Dans le deuxième thème, nous étudions la réponse du champ géomagnétique à une impulsion de la pression dynamique (Psw) du vent solaire qui a atteint la magnétosphère le 24 août 2005. En utilisant les données du champ géomagnétique à haute résolution fournies par 15 stations au sol et les données des satellites Geotail, TC-1 et TC-2, nous avons étudié les pulsations géomagnétiques aux latitudes aurorales dues à la décroissance brusque de Psw dans la limite arrière de l'impulsion. Les résultats montrent que la décroissance brutale de Psw peut exciter une pulsation globale dans la gamme de fréquence 4.3-11.6 mHz. L'inversion des polarisations entre deux stations des latitudes aurorales, la densité spectrale de puissance (PSD) plus grande près de la latitude de résonance et la fréquence augmentant avec la diminution de la latitude indiquent que les pulsations sont associées avec la résonance des lignes de champ (FLR). La fréquence résonante fondamentale (la fréquence du pic de PSD entre 4.3-5.8 mHz) dépend du temps magnétique local et est maximale autour du midi local magnétique. Cette caractéristique est due au fait que la dimension de la cavité magnétosphérique dépend du temps et est minimale à midi. Une onde harmonique est aussi observée à environ 10 mHz, qui est maximale dans le secteur jour, et est fortement atténuée lorsque on va vers le coté nuit. La comparaison entre les PSDs des pulsations provoquées pas l'augmentation brutale et la décroissance brutale de Psw montre que la fréquence des pulsations est inversement proportionnelle à la dimension de la magnétopause. Puisque la résonance des lignes de champ (FLR) est excitée par des ondes cavité compressée/guide d'onde, ces résultats indiquent que la fréquence de résonance dans la cavité magnétosphérique/guide d'onde est due non seulement aux paramètres du vent solaire mais est aussi influencée par le temps magnétique local du point d'observation.

La gestion en temps réel (GTR) est considérée comme une solution économiquement efficace pour réduire les déversements par temps de pluie car elle optimise la capacité disponible des réseaux d'assainissement. La GTR permet d'éviter la construction de volumes de rétention supplémentaires, d'augmenter l'adaptabilité du réseau aux changements de politiques de gestion de l'eau et surtout d'atténuer l'impact environnemental des déversoirs d'orage. À la suite de l'intérêt croissant pour la GTR fondée sur la qualité de l'eau (QBR), cette thèse démontre une stratégie simple et efficace pour les charges polluantes déversées par temps de pluie. La performance de la stratégie QBR, basée sur la prédiction des courbes masse-volume (MV), est évaluée par comparaison avec une stratégie typique de GTR à base hydraulique (HBR). Une étude de validation de principe est d'abord réalisée sur un petit bassin versant de 205 ha pour tester le nouveau concept de QBR en utilisant 31 événements pluvieux sur une période de deux ans. Par rapport à HBR, QBR offre une réduction des charges déversées pour plus d'un tiers des événements, avec des réductions de 3 à 43 %. La stratégie QBR est ensuite mise en oeuvre sur le bassin versant de Louis Fargue (7700 ha) à Bordeaux, France et comparée à nouveau à la stratégie HBR. En implémentant QBR sur 19 événements pluvieux sur 15 mois, ses performances sont constantes et apportent des avantages précieux par rapport à HBR, 17 des 19 événements ayant une réduction de charge variant entre 6 et 28.8 %. La thèse évalue en outre l'impact de l'incertitude de prédiction de la courbe MV (due à l'incertitude de prédiction du modèle) sur la performance de la stratégie QBR, en utilisant un événement pluvieux représentatif. La marge d'incertitude qui en résulte est faible. En outre, l'étude de sensibilité montre que le choix de la stratégie QBR ou HBR doit tenir compte des dimensions réelles des bassins et de leur emplacement sur le bassin versant. / Real time control (RTC) is considered as a cost-efficient solution for combined sewer overflow (CSO) reduction as it optimises the available capacity of sewer networks. RTC helps to prevent the need for construction of additional retention volumes, increases the network adaptability to changes in water management policies, and above all alleviates the environmental impact of CSOs. Following increasing interest in water quality-based RTC (QBR), this thesis demonstrates a simple and nothing-to-lose QBR strategy to reduce the amount of CSO loads during storm events. The performance of the QBR strategy, based on Mass-Volume (MV) curves prediction, is evaluated by comparison to a typical hydraulics-based RTC (HBR) strategy. A proof-of-concept study is first performed on a small catchment of 205 ha to test the new QBR concept using 31 storm events during a two-year period. Compared to HBR, QBR delivers CSO load reduction for more than one third of the events, with reduction values from 3 to 43 %. The QBR strategy is then implemented on the Louis Fargue catchment (7700 ha) in Bordeaux, France and similarly compared with the HBR strategy. By implementing QBR on 19 storm events over 15 months, its performance is consistent, bringing valuable benefits over HBR, with 17 out of 19 events having load reduction varying between 6 and 28.8 %. The thesis further evaluates the impact of MV curve prediction uncertainty (due to model prediction uncertainty) on the performance of the QBR strategy, using a representative storm event. The resulting range of uncertainty is limited. Besides, results of the sensitivity study show that the choice of the QBR or HBR strategy should take into account the current tank volumes and their locations within the catchment.