TP 6
Modélisation du fonctionnement du réseau d'assainissement
pluvial: Application sur le logiciel SWMM (Storm Water
Management Model )
1
Modélisation du fonctionnement du réseau d'assainissement
pluvial avec SWMM
1. But:
Le but est de modéliser le fonctionnement du réseau d'assainissement pluvial d'une
zone résidentielle de 4,86 ha. Ce réseau schématisé sur la figure 1 ci- dessous comprend les
objets suivants :
• 3 sous bassins versants (Subcatchments) S1 à S3,
• 4 collecteurs (Conduits) C1 à C4
• 4 nœuds de jonction (Junctions) J1 à J4
• 1 exutoire (Outfalls) Out1
Figure 1: Schéma du réseau
Les étapes de la modélisation seront :
1) Création des objets sur le plan (Study Area Map) et entrée des données de chaque
objet (Properties),
2) Simulation de la réponse quantitative du système à un événement pluvieux ponctuel de
76,2 mm en 6 heures.
2. Création d'un nouveau projet:
La première tâche consiste à créer un nouveau projet SWMM, lui donner un titre et à
l'enregistrer dans un fichier de sauvegarde. Pour cela :
1- Exécuter EPA SWMM et sélectionner File>>New dans la barre de menu principal pour
créer un nouveau projet.
2- Sélectionner la rubrique "Title/Notes" dans l'explorateur de données à gauche de l'écran et
cliquer sur le bouton .
2
3- Dans la fenêtre qui apparaît alors, taper "Exercice d'initiation" comme titre de projet et
cliquer sur OK pour fermer la fenêtre.
4- Dans le menu, sélectionner File>>Saves as
5- Dans la boîte de dialogue, sélectionner un dossier et un nom de fichier sous lequel
enregistrer le projet, par exemple "Exo_initiation.inp" (l'extension .inp sera
automatiquement ajoutée au nom de fichier).
6- Cliquer sur "Enregistrer" pour sauvegarder le projet dans le fichier.
Au fur et à mesure de la saisie des données, il faudra penser à enregistrer régulièrement les
modifications en cliquant sur le bouton ou en sélectionnant File>>Save.
Ensuite, nous allons choisir et renseigner certaines options par défaut. En effet l'utilisation de
ces options par défaut nous simplifiera plus tard le travail.
1- Sélectionner Project>>Defaults pour ouvrir la boîte de dialogue des options par défaut
du projet.
2- Dans le tableau de l'onglet ID Label, renseigner la colonne ID Prefixes comme indiqué
sur la figure 2. Cela aura pour effet d'étiqueter chaque objet avec un Identifiant composé d'un
préfixe et d'un numéro automatique.
Figure 2: Étiquetage des objets
3- Dans l'onglet Subcatchments, entrer les valeurs suivantes :
3
Valeur à
Property Unité SI Signification rentrer
Area ha Surface du bassin versant 1.62
Paramètre du modèle de ruissellement
(modèle de l'onde cinématique) : largeur du
canal rectangulaire auquel est assimilée la
surface du BV pour y calculer la vitesse de
ruissellement selon la formule de Manning-
Strickler. C'est le rapport de la surface du BV
à la longueur du BV. Cette longueur est
Width m 122
obtenue en considérant la moyenne de
plusieurs parcours hydrauliques sur la surface
du BV.
Elle doit refléter, davantage les écoulements
lents sur les surfaces perméables que les
écoulements rapides sur les surfaces
imperméables. C'est un paramètre de calage.
% Slope % Pente moyenne du bassin versant 0.5
Pourcentage de surface imperméable sur le
% Imperv. % 50
BV
Coefficient de Manning- Strickler pour les
N imperv. s/m1/3 0.01
surfaces imperméables
Coefficient de Manning- Strickler pour les
N-perv. s/m1/3 0.1
surfaces perméables
Hauteur d'eau stockée sur les surfaces
imperméables du BV(dépressions du sol,
Dstore-Imperv. mm 1.3
interception par la végétation, mouillage des
surfaces)
Hauteur d'eau stockée sur les surfaces
Dstore-Perv mm 1.3
perméables du BV
Pourcentage de la surface imperméable
dépourvue de stockage dans les dépressions
%Zero-Imperv % 25
du sol
Modèle d'infiltration : Trois modèles sont
disponibles: GREEN-
Infil. Model - Horton AMPT
- Green Ampt
- courbes SCS
Potentiel de succion ou potentiel matriciel
Suction Head mm 89
capillaire au niveau du front d'humidification
Conductivity mm/h Conductivité hydraulique en milieu saturé 12.7
Déficit Initial : variation de la teneur en eau
du sol (Θf- Θi): eneur en eau finale (à
Initial Deficit fraction 0.26
saturation) - teneur en eau initiale.
4
4- Dans l'onglet Nodes/Links, entrer les valeurs suivantes :
Valeur à
Property Unité SI Signification
Rentrer
Cote inférieure du nœud (cote du radier du
Node Invert m 0
regard de jonction)
Profondeur Maximale du Nœud : distance
Node Max. Depth m verticale entre le radier du regard de 1.2
jonction et le terrain naturel
Surface d'eau stockée au dessus du nœud
après débordement. Si l'option de calcul
correspondante est activée,
Node Ponded Area m² le volume de débordement est stocké puis 0
évacué par le réseau lorsque celui-ci en a
retrouvé la capacité.
Conduit Length m Longueur de la conduite 120
Conduit Geometry Géométrie de la conduite
Nombre de conduites en parallèle
Barrels
Forme de la section CIRCULAR
Shape
Hauteur d'eau maximale dans la section
Max. Depth m (diamètre pour une section circulaire) 0.3
Rugosité de la conduite au sens de 0.01
Conduit Roughness s/m1/3 Manning-Strickler
Unités de débit. Choisir une unité métrique
(CMS :m/s, LPS :L/s, MLD : 1000m3/j)
implique que toutes les autres grandeurs du
logiciel sont exprimées en unités métriques.
Flow Units Les valeurs déjà rentrées ne sont pas CMS
automatiquement converties lorsqu'on
change de système d'unités.
5
Décalage du collecteur par rapport au fond
du regard. La position du radier du
collecteur peut être indiquée sous forme
d'une hauteur au dessus du fond du regard
(DEPTH) = distance entre les points 1 et 2 )
ou sous forme d'une cote absolue
(ELEVATION) = cote absolue du point 1
Link Offsets DEPTH
Modèle de transfert. Trois modèles sont
disponibles :
Steady Flow : Écoulement permanent et
uniforme dans chaque tronçon et à chaque
pas de temps (Translation simple des
Kinematic
Routing Model hydrogrammes)
Wave
Kinematic Wave : Modèle de l'onde
cinématique
Dynamic Wave : Modèle de l'onde
dynamique (Résolution des équations de
Barré de Saint Venant)
5- Cliquer sur OK pour accepter ces choix et fermer la boîte de dialogue. Pour appliquer ces
choix à tous les futurs projets, cocher la case "Save as defaults for all new project".
Dans l'étape suivante, nous allons activer les options d'affichage pour que les symboles et les
noms des objets s'affichent automatiquement sur le plan au fur et à mesure de la création des
objets, et pour que les conduites portent des flèches indiquant le sens de l'écoulement:
1- Sélectionner Tools>>Map Display Options pour faire apparaître la boîte de dialogue des
options d'affichage du plan.
2- Sélectionner la rubrique Subcatchment et choisir Diagonal pour Fill Style. Ceci aura
pour effet de représenter les bassins versants sous forme d'un polygone hachuré
diagonalement.
3- Sélectionner la rubrique Node et choisir une taille de symbole (Node Size) de 5.
4- Sélectionner la rubrique Annotations et cocher les cases SubCatch IDs, Node Ids et
Link Ids pour faire s'afficher sur le plan les noms des bassins versants, des nœuds et des
conduites.
Laisser toutes les autres cases non cochées.
6
5- Enfin, sélectionner la rubrique Flow Arrows et choisir Filled pour Arrow Style. Ceci aura
pour effet de dessiner des flèches pleines.
6- Cliquer sur OK pour accepter ces choix et fermer la boîte de dialogue.
Avant de dessiner les objets sur le plan, il faut spécifier ses dimensions :
1- Sélectionner View>>Dimensions pour faire apparaître la boîte de dialogue des
dimensions du plan.
2- Pour cet exemple, conserver les valeurs par défaut.
Enfin, dans la barre d'état, en bas à gauche de la fenêtre principale, s'assurer que l'option
AutoLength est bien sur Off. (Cette option lorsqu'elle est activée permet de calculer
automatiquement les longueurs de conduites sur un plan établi en coordonnées réelles).
3. Dessin des objets
A présent nous allons ajouter des éléments au plan du projet en commençant par les bassins
versants :
1- Cliquer sur le bouton de la barre à outils des objets pour commencer à dessiner des
bassins versants. (Si cette barre à outils n'est pas visible, sélectionner
View>>Toolsbar>>Object pour la faire apparaître). Le pointeur de la souris prend alors la
forme d'un crayon.
2- Déplacer le pointeur de la souris sur le plan jusqu'à l'emplacement d'un des sommets du
périmètre du bassin versant S1 et cliquer (clic gauche).
3- Faire de même avec les autres sommets et faire un clic droit (ou appuyer sur Entrée) pour
fermer le périmètre et créer ainsi le bassin versant S1. Appuyer sur ESC pour annuler le
dessin d'un bassin versant en cours de saisie. Noter que si nécessaire, il sera possible par la
suite de modifier ou de déplacer le bassin versant.
4- Répéter l'opération pour les bassins versants S2 et S3.
Observer comment les objets sont automatiquement étiquetés à mesure qu'ils sont créés.
La prochaine étape consiste à dessiner les nœuds de jonctions et l'exutoire du réseau.
1- Cliquer sur le bouton de la barre à outils des objets pour commencer à dessiner des
nœuds de jonction.
2- Déplacer le pointeur de la souris sur le plan jusqu'à l'emplacement du nœud de jonction J1
et cliquer (clic gauche). Faire de même pour le nœuds J2 à J4.
7
3- Pour ajouter l'exutoire Out1, cliquer sur le bouton de la barre à outils des objets,
déplacer le pointeur de la souris sur le plan jusqu'à l'emplacement de l'exutoire et cliquer (clic
gauche).
A ce stade, le plan doit ressembler à la figure suivante :
Figure 3: Bassins versants et nœuds de jonction de l'exemple
A présent, nous allons créer les collecteurs qui relient les nœuds de jonction (il faut avoir créé
les nœuds avant de créer les conduites). Commençons par le collecteur C1 qui relie les nœuds
J1 et J2.
1- Cliquer sur le bouton de la barre à outils des objets. Le pointeur de la souris prend
alors la forme d'une croix.
2- Cliquer sur le nœud J1. Le pointeur de la souris prend alors la forme d'un crayon.
3- Déplacer la souris jusqu'au nœud J2. Une ligne se dessine entre les deux nœuds. Cliquer
(clic gauche) pour créer le tronçon. L'opération peut être annulée en cours de saisie en
cliquant sur le bouton droit de la souris ou en appuyant sur la touche ESC.
4- Répéter la procédure pour les collecteurs C2 à C4
Tous les collecteurs ont été représentés par des droites, mais il est possible de les représenter
par des polylignes dont les points intermédiaires sont créés en faisant des clic-gauche.
Pour finir, nous allons ajouter une station de mesure pluviométrique :
1- Cliquer sur le bouton de la barre d'outils des objets,
2- Déplacer la souris sur la carte jusqu'à l'emplacement souhaité pour la station
pluviométrique et cliquer (clic gauche).
8
A ce stade, le dessin de tous les composants du système est terminé. Le dessin obtenu doit
ressembler à celui de la figure 1. Si un objet (station pluviométrique, bassin versant ou nœud
est mal placé, il est possible de le déplacer selon la procédure suivante :
1- Cliquer sur dans la boîte à outils de la carte pour activer l'outil de sélection ,
2 - Cliquer sur l'objet à déplacer,
3- Déplacer l'objet tout en maintenant le bouton de la souris enfoncé puis relâcher le bouton
de la souris.
Pour modifier le contour des bassin versants :
1- Lorsque l'outil de sélection est actif, cliquer sur le centroïde du bassin versant à modifier.
(le centroïde est symbolisé par un carré noir à l'intérieur du bassin versant)
2- Cliquer sur le bouton de la boîte à outils de la carte pour activer l'outil d'édition des
sommets.
3- Cliquer sur un des sommets du contour du bassin versant pour le sélectionner (noter que
celui-ci devient alors noir).
4- Déplacer ce sommet vers sa nouvelle position tout en maintenant le bouton de la souris
enfoncé puis relâcher le bouton de la souris.
5- Si nécessaire, des sommets peuvent être rajoutés ou supprimés sur le contour d'un bassin
versant en faisant un clic droit dessus et en sélectionnant l'action appropriée dans la fenêtre
qui apparaît.
Cette même procédure peut-être utilisée pour modifier la forme d'un collecteur.
[Link] les propriétés des objets
A mesure que les objets sont dessinés, SWMM leur affecte des propriétés par défaut. Pour
modifier le contenu ou la valeur d'une propriété spécifique d'un objet, il faut appeler l'éditeur
de propriétés de cet objet (figure 4). Il y a différentes façons de procéder :
Si l'éditeur de propriétés est déjà visible, il suffit de cliquer sur l'objet sur la carte
ou de le sélectionner dans la liste des données à gauche de la carte. Si l'éditeur de propriétés
n'est pas visible, on peut le faire apparaître grâce à l'une des actions suivantes :
• Double cliquer sur l'objet sur la carte
• Faire un clic droit sur l'objet et choisir l'option « properties » dans la fenêtre qui
apparaît
• Sélectionner l'objet dans la liste des données à gauche de la carte et cliquer sur le
bouton de l'explorateur de données.
9
Figure 4: Editeur de propriétés
Lorsque la fenêtre de l'éditeur de propriétés est active, il suffit de cliquer sur F1 pour obtenir
une description plus détaillée des propriétés listées.
Deux propriétés doivent être renseignées pour les bassins versants : la station de
mesure pluviométrique qui décrit les précipitations sur le bassin versant ainsi que le nœud du
réseau qui reçoit l'eau ruisselant sur le bassin versant. Mais puisque ici tous les bassins
versants utilisent la même station pluviométrique, nous pouvons utiliser une méthode
permettant d'affecter cette propriété à tous les bassins versants à la fois :
1- Dans le menu principal, sélectionner Edit>>Select All,
2- Ensuite sélectionner Edit>>Group Edit pour faire apparaître la boîte de dialogue
d'édition de groupe (figure 5)
3- Choisir Subcatchment pour le type d'objet à sélectionner, Rain Gage comme propriété à
éditer et taper « Pluvio1 » comme nouvelle valeur.
4- Cliquer sur OK pour affecter la station pluviométrique « Pluvio1 » à tous les bassins
versants. Une boîte de dialogue apparaît alors, signalant que les trois bassins versants
ont été modifiés. Choisir « No » pour ne pas continuer l'édition.
Figure 5: Editeur de groupe
10
Comme chaque bassin versant a un nœud exutoire différent, on doit maintenant procéder
pour chaque bassin versant de la façon suivante :
1- Double cliquer sur le bassin versant S1 ou le sélectionner dans la liste de données et
cliquer sur pour faire apparaître l'éditeur de propriétés.
2- Taper « J1 » dans la rubrique Outlet et taper Entrée. Remarquer qu'une ligne en
pointillés s'est dessinée entre le bassin versant et son exutoire.
3- Cliquer sur le bassin versant S2 et entrer « J2 » dans la rubrique Outlet
4- Cliquer sur le bassin versant S3 et entrer « J3 » dans la rubrique Outlet
Nous souhaitons également considérer le fait que le bassin versant S3 est moins urbanisé que
les autres. Sélectionner alors S3 dans l'éditeur de propriétés pour lui affecter un
coefficient d'imperméabilisation (%Imperv) de 25%.
Les cotes du radier (Invert Elevation) des regards ou noeuds de jonction et de l'exutoire
doivent ensuite être renseignées. Comme pour les bassins versants, sélectionner
individuellement les différents nœuds du système et entrer les données ci-dessous :
Cote radier
Nœud
J1 29.26m
J2 27.43m
J3 28.35m
J4 26.82m
Out1 25.91m
Dans notre exemple, un seul des collecteurs a une propriété dont la valeur diffère des valeurs
par défaut. Il s'agit du collecteur C4 dont le diamètre est de 450 mm au lieu de 300 mm. Pour
changer cette valeur, sélectionner le collecteur C4 dans l'éditeur de propriétés et entrer la
valeur 0,450 dans la rubrique Max Depth.
Pour saisir des données pluviométriques, il faut renseigner les propriétés de la station
pluviométrique. Sélectionner l'objet « Pluvio1 » dans la liste des données et éditer ses
propriétés.
Entrer les valeurs ci-dessous :
Unité Valeur à
Property Signification
SI Rentrer
Format de représentation des données
pluviométriques :
INTENSITY : intensité en mm/h
Rain Format INTENSITY
mesurée sur chaque intervalle
de temps
VOLUME : hauteur précipitée en mm
11
mesurée sur chaque intervalle de
temps
CUMULATIVE : cumul de hauteur
précipitée en mm depuis le début de
l'événement pluvieux
h ou Intervalle de temps de mesure du
Rain Interval 01:00
hh:mm pluviomètre
Source des données pluviométriques :
TIMESERIES : données entrées par
Data Source l'utilisateur TIME
FILE : données dans un fichier
externe.
Nom de la série temporelle de
Series Name pluviométrie.
TS1
Double cliquer pour éditer la série
Nom du fichier externe contenant les
File Name données
Numéro d'enregistrement de la station
Station No
Unité de hauteur des données
pluviométriques dans le fichier
Rain Units
externe : IN (pouces) ou MM
(millimètres)
Comme mentionné précédemment, nous allons simuler le fonctionnement du réseau pour une
pluie de projet de 76,2 mm sur 6 heures. Les valeurs horaires de l'intensité de pluie vont
constituer la série temporelle nommée TS1. Pour cela il faut créer un objet « Time Series » et
y rentrer les données.
1- Dans l'explorateur d'objets à gauche de l'écran sélectionner la catégorie « Time Series »
2- Cliquer sur le bouton pour faire apparaître la fenêtre de l'éditeur de série
temporelle(figure 6)
3- Entrer « TS1 » dans le champ « Time Series Name »
4- Entrer les valeurs figurant sur la figure 6 dans les colonnes "Time" et "Value".
5- Cliquer sur le bouton "View", pour visualiser le graphique correspondant puis cliquer sur
OK pour valider.
12
Figure 6: Editeur de séries temporelles (Time Series)
Toutes les données du projet ayant été renseignées, enregistrer les dernières modifications en
cliquant sur le bouton ou en sélectionnant File>>Save.
Les données du projet ont été enregistrées dans un fichier au format texte. Il est
possible de visualiser son contenu en sélectionnant la commande Project>>Details dans le
menu principal.
Pour ouvrir, lors d'une prochaine session de travail, un projet préalablement enregistré, il
suffit de sélectionner la commande File>>Open dans le menu principal.
5. Réaliser une simulation:
5.1 Options de simulation
Avant de simuler le fonctionnement du réseau, il faut fixer un certain nombre de paramètres
qui détermineront la façon dont sera conduite la simulation. Pour cela :
1- Sélectionner la rubrique "Options" dans l'explorateur de données à gauche de l'écran et
cliquer sur le bouton .
2- Dans l'onglet "General" de la fenêtre qui apparaît, vérifier que les cases
"Rainfall/Runoff" et "Flow Routing" de la rubrique "Process Models" sont cochées (pour
réaliser respectivement la transformation pluie/débit sur les bassins versants et pour
modéliser l'écoulement dans le réseau). Dans la rubrique "Routing Model", cocher la
13
case "Kinematic Wave" pour simuler l'écoulement dans le réseau selon le modèle de
transfert de l'onde cinématique. Dans la rubrique "Infiltration Model", la case "Green
Ampt" doit être cochée et dans la rubrique "Miscellaneous", toutes les cases doivent être
décochées.
Figure 7: Options de simulation
Options de la
rubrique Routing Signification
Model
Modèle de la translation simple en régime permanent uniforme
Modèle le plus simple possible qui considère à chaque pas de temps
de calcul, que l'écoulement est permanent et uniforme. Les
hydrogrammes entre l'entrée et la sortie des conduites, ne subissent
alors ni décalage dans le temps, ni amortissement.
Steady Flow
Ne peut représenter ni les effets de stockage de l'eau dans les réseaux,
ni les phénomènes d'influences aval, ni les inversions du sens de
l'écoulement, ni les écoulements en charge. Son utilisation est
restreinte aux réseaux ramifiés (non maillés). Il n'est réellement
utile que pour une analyse préliminaire d'écoulements variant peu
dans le temps.
Modèle de l'onde cinématique
Il résout les deux équations de Barré de Saint Venant en négligeant,
dans l'équation de conservation de laquantité de mouvement, les deux
Kinematic Wave
termes d'inertie et en faisant l' approximation que la pente de la
surface libre est égale à la pente du radier.
Dès que le débit à l'entrée d'un collecteur est supérieur au débit à
pleine section du collecteur, l'excédent est soit perdu, soit stocké au
14
dessus du nœud d'entrée puis réinjecté dans le collecteur lorsque
celui-ci a retrouvé sa capacité d'évacuation.
Permet de représenter des écoulements variant à la fois dans le temps
et dans l'espace, ce qui se traduit par un décalage temporel ainsi
qu'un amortissement des hydrogrammes entre l'entrée et la sortie d'un
collecteur.
Ne prend en compte ni les influences aval, ni les inversions du sens
d'écoulement, ni les écoulements en charge et son utilisation se limite
aux réseaux ramifiés (non maillés). Utilisable avec des pas de temps
de calcul relativement longs, de l'ordre de 5 à 15 minutes
Modèle de l'onde dynamique
Résout les équations complètes de Barré de Saint Venant et produit
donc théoriquement les résultats les plus précis.
Permet de représenter les écoulements en charge dans les collecteurs,
lorsque le débit y excède le débit à pleine section. Lorsqu'en un
nœud, la cote piézométrique dépasse la cote du terrain naturel, il y
a débordement et le débit en excès est alors, soit perdu soit stocké au
dessus du nœud et réinjecté dans le réseau en fonction de ses
capacités d'évacuation.
Permet de représenter les effets de stockage dans le réseau, les
Dynamic Wave phénomènes d'influence aval, l'inversion du sens de l'écoulement et
les écoulements en charge. Il peut être appliqué à n'importe-quel
réseau, y compris les réseaux maillés et/ou comportant des
diffluences.
Modèle à choisir quand on peut s'attendre à des influences aval
significatives, notamment lorsque le réseau comporte des limitations
de débit à l'aval ou des systèmes de régulation tels que déversoirs ou
orifices calibrés.
Requiert cependant des pas de temps de calcul beaucoup plus courts,
de l'ordre de la minute, ou moins.(SWMM réduit automatiquement la
durée du pas de temps choisie par l'utilisateur, afin de garantir la
stabilité numérique du calcul).
Options de la
rubrique Signification
Miscellaneous
Autorise le modèle à prendre en compte le stockage en surface des
Allow Ponding volumes excédentaires lorsque la capacité de transit des collecteurs
est dépassée.
Inclut dans le rapport de simulation, la liste des actions de contrôle
Report Control
ponctuelles (discrete control action) qui se produisent selon les règles
Actions
de contrôle (Control Rules) définies dans projet.
Report Input Inclut dans le rapport de simulation, un résumé des données du
Summary projet.
Utilise, lorsque l'écoulement est permanent, les résultats du dernier
calcul au lieu de procéder à un nouveau au calcul. L'écoulement est
Skip Steady
jugé permanent sur un pas de temps si le débit externe sur chaque
Periods
nœud est inférieur à 15 L/s et si la différence relative entre les débits
entrant et sortant du système n'excède pas 5%.
15
Valeur minimale autorisée pour la pente des conduites (%). Si aucune
Minimum Conduit
valeur n'est imposée, SWMM utilisera la valeur par défaut de
Slope
0.00035.
3- Dans l'onglet "Dates", fixer la fin de la simulation (End Analysis) à 12:00.
4- Dans l'onglet "Time Steps", fixer le pas de temps de calcul (Routing) à 60s.
Onglet
Signification
Time Steps"
Reporting Pas de temps pour la présentation des résultats de calculs
Runoff : Pas de temps de calcul des débits ruisselés à l'exutoire des bassins
DryWeather versants.
RunOff : Pas de temps de calcul de l'accumulation des polluants à la surface
WetWeather des bassins versants en temps sec (modèles de qualité)
Pas de temps de calcul pour la modélisation des écoulements dans le
Routing
réseau.
5- Cliquer sur OK pour fermer la fenêtre des options de simulation.
5.2 Lancer une simulation
Pour lancer une simulation, sélectionner Project>>Run Simulation dans le menu
principal ou cliquer sur le bouton . En cas de problèmes lors de la simulation, un rapport
de simulation apparaîtra à l'écran, décrivant les erreurs rencontrées. Une fois que la
simulation s'est effectuée avec succès, il y a plusieurs manières d'en visualiser les résultats :
5.3 Consulter le rapport de simulation
la simulation. Pour le consulter, sélectionner Report>>Status dans le menu principal.
Le rapport affiché indique principalement :
• Une bonne qualité de la simulation puisque les erreurs sur la conservation de la
masse (Continuity Errors) sont négligeables, pour le modèle de ruissellement
comme pour le modèle de transfert (Runoff Quantity Continuity Error = -0,04 %
et Flow Routing Continuity Error = -0,032%).
• Sur les 76,2 mm de pluie précipitée, 44,5 mm se sont infiltrés et 31,4 mm ont ruisselé.
• La rubrique des débordements aux nœuds (Node Flooding Summary) signale
un débordement au nœud J2.
• La rubrique des surcharges des conduites (Conduit Surcharging Summary) signale
que la conduite C2, juste à l'aval du nœud J2 a été surchargée et donc apparaît
comme sousdimensionnée.
16
5.4 Représenter les résultats sur le plan
Les résultats de simulation (comme certains paramètres descriptifs du réseau tels que les
surfaces de bassins versants, les cotes radier des regards et les diamètres des
collecteurs) peuvent être représentés sur la carte avec des couleurs fonction de leurs
valeurs.
1- Sélectionner l'onglet "Map" dans l'explorateur à gauche de l'écran.
2- Dans la catégorie "Themes", sélectionner les variables à visualiser pour les bassins
versants, les noeuds et les conduites. Par exemple, les débits ruisselés (Runoff) pour les
bassins versants (Subcatchments) et les débits (Flow) dans les collecteurs (Links).
3- Les légendes des couleurs s'affichent alors sur la carte. Pour supprimer ou faire apparaître
ces légendes, sélectionner View>>Legend.
4- Les cadres des légendes sont déplaçables vers une autre position en maintenant le bouton
gauche de la souris enfoncé.
5- Pour modifier une légende (couleurs et amplitude des différentes classes), sélectionner
View>>Legend>>Modify puis la légende à modifier ou faire simplement un clic droit sur le
cadre de légende. Pour voir les valeurs numériques des variables s'afficher sur la carte,
sélectionner Tools>>Map Display Option et sélectionner la rubrique "Annotation".
Cocher alors les cases "Subcatchment Values", "Node Values" et "Link Values" selon les
annotations désirées.
6- Dans l'explorateur à gauche de l'écran, les boutons de la rubrique "Time Period"
permettent de consulter les résultats de la simulation à différents instants.
7- Les boutons "magnétoscope" de la rubrique "Animator" permettent de contrôler une
animation de la carte en fonction du temps.
5.5 Visualiser les résultats sous forme de graphes temporels
Pour créer un graphe temporel à partir d'un résultat de simulation :
1- Sélectionner Report>>Graph>>Time Series ou cliquer sur le bouton dans la
barre d'outils standard.
2- La boîte de dialogue qui apparaît, permet de choisir les variables à représenter en fonction
du temps.
Pour notre exemple, la boîte de dialogue Times Series Plot peut être utilisée pour
représenter le débit dans les collecteurs C1 et C2 :
1- Sélectionner "Links" dans "Object Category"
2- Cocher "Flow" comme variable à représenter
17
3- Cliquer sur la conduite C1 (soit sur la carte soit dans l'explorateur d'objet) et cliquer sur
le bouton pour ajouter cette conduite à la liste des conduites concernées par le graphe.
Faire de même avec la conduite C2.
4- Cliquer sur OK pour créer le graphique.
Une fois le graphique créé, il est possible de :
• Modifier son apparence en sélectionnant Report>>Customize ou en faisant un clic
droit sur le graphique.
• Le copier dans le presse papier et de le coller dans un autre logiciel en
sélectionnant Edit>>Copy ou en cliquant sur le bouton dans la barre d'outils
standard.
• L'imprimer en sélectionnant File>>Print ou File>>Print Preview (utiliser au
préalable File>>Page Setup pour définir les marges, l'orientation de la page, etc.).
5.6 Créer un profil en long
SWMM permet de générer des profils en long montrant comment varie la ligne d'eau dans
une succession de tronçons du réseau. Créons par exemple un tel profil pour les collecteurs
reliant le nœud J1 à l'exutoire général du réseau (Out1).
1- Sélectionner Report>>Graph>>Profile ou cliquer sur le bouton de la barre d'outils
standard.
2- Dans la fenêtre qui apparaît (cf. figure 8), entrez J1 dans le champ "Start Node" pour
indiquer le point de départ du profil, ou bien cliquer sur le nœud J1 sur la carte ou dans
l'explorateur d'objets et cliquer sur le bouton .
Figure 8: Définition du profil en long
3- Faire de même avec le noeud Out1 dans le champ "End Node"
18
4- Cliquer sur "Find Path". Une liste ordonnée des conduites reliant le nœud de départ au
nœud d'arrivée s'affiche alors dans la rubrique "Links in Profile". Il est possible si
nécessaire d'éditer cette liste pour la modifier.
5- Cliquer sur OK pour créer le profil montrant la position de la ligne d'eau à l'instant
indiqué dans l'onglet "Map" de l'explorateur, à gauche de l'écran.
Lorsque que l'on évolue dans le temps en utilisant les boutons de contrôles de l'onglet "Map"
de l'explorateur, ou en lançant une animation avec l"Animator", on voit évoluer la forme et
la position de la ligne d'eau. On observe notamment qu'il y a débordement au noeud J2
pendant une heure, deux heures après le début de la pluie.
L'apparence du profil en long peut être modifiée et il peut être copié et imprimé, selon les
mêmes procédures que pour les graphes temporels.
5.7 Réaliser une simulation avec le modèle de l'onde dynamique
Ce modèle est plus précis que celui de l'onde cinématique mais il requiert un temps de
calcul supérieur car les pas de temps doivent être beaucoup plus courts pour garantir la
stabilité numérique du modèle.
Dans notre exemple, le modèle de l'onde cinématique a montré que la conduite C2 était
surchargée au cours de l'épisode pluvieux. Il est donc probable qu'en se mettant en charge,
elle puisse faire transiter un débit supérieur à son débit à pleine section. Le modèle de l'onde
dynamique est capable de représenter cette mise en charge.
Pour réaliser une simulation selon le modèle de l'onde dynamique :
1- Dans l'explorateur d'objets, sélectionner la catégorie "Options" et cliquer sur le bouton
2- Dans la rubrique "Routing Model" de l'onglet "General", cocher "Dynamic Wave"
3- Dans l'onglet "Dynamic Wave" reproduire les choix de la figure 9 ci-dessous :
19
Figure 9: Options du modèles de l'onde dynamique
Options du
modèle de l'onde Signification
dynamique
Choix du mode de traitement des termes d'inertie de l'équation de
Barré de Saint Venant:
KEEP : les termes d'inertie sont pleinement conservés partout et à
tout instant.
DAMPEN : les termes d'inertie sont réduits lorsque l'écoulement
s'approche du régime critique et sont ignorés lorsque l'écoulement
Inertial Terms
devient torrentiel.
IGNORE : les termes d'inertie sont ignorés partout et à tout instant,
ce qui revient au modèle de l'onde diffusante.
Choix du critère permettant de savoir que l'écoulement devient
torrentiel :
SLOPE : pente de la surface libre supérieure à la pente du radier
Define SuperCritical
FROUDE N° : nombre de Froude supérieur à 1
Flow by
BOTH : l'une ou l'autre des deux conditions précédentes
(recommandé)
Equation utilisée pour calculer la perte de charge unitaire dans les
Force MainEquation écoulements en charge.
20
USE : Utilisation ou non d'un pas de temps variable, calculé
automatiquement.
ADJUSTMENT FACTOR : valeur du facteur d'ajustement (ou
Variable TimStep facteur de sécurité) applicable à ce pas de temps (75%
recommandé).
Le pas de temps variable ne doit pas être inférieur à 0,5s ni excéder
le pas de temps fixe de calcul choisi par l'utilisateur.
Conduit Lengthening
Time Step
Section minimale des nœuds de jonction utilisée pour le calcul des
Minimum Surface
évolutions de la hauteur d'eau dans les regards. La valeur par défaut
Area
est 1,167 m², ce qui correspond à un regard de 4feet de diamètre.
4- Cliquer sur OK et lancer la simulation.
21