Modelisation Hydraulique
Modelisation Hydraulique
Modelisation Hydraulique
MODELISATION HYDRAULIQUE
Email: [email protected]
1
OBJECTIFS
• Cerner la notion de modélisation hydraulique;
• Indiquer les principes de base de la modélisation
hydraulique
e f g
5
1. Introduction
❑ Les 4 fonctions principales d'un modèle
• On représente le modèle par une fonction de
transformation ‘’f’’ produisant une sortie ‘’s’’
en fonction d’une entrée ‘’e’’ et des variables
de contrôle ‘’c’’.
6
1. Introduction
❑ Les 4 fonctions principales d'un modèle
• Fonction cognitive : le modèle sert à représenter les relations qui existent entre
variables d'entrée et variables de sorties du système;
• Fonction prévisionnelle : le modèle sert à prévoir comment évolueront les
variables de sortie du système;
7
1. Introduction
❑ Les 4 fonctions principales d'un modèle
• Fonction décisionnelle : le modèle sert à déterminer comment agir sur les
variables d'entrée pour atteindre les objectifs que l'on s'est fixé sur les variables de
sortie;
• Fonction normative : le modèle sert à représenter les relations souhaitables entre
les variables d'entrée et les variables de sortie.
8
1. Introduction
• Modélisation: schématisation, une simplification des
processus réels.
9
1. Introduction
• Relation entre hydraulique et l’hydrologie:
❑l’hydrologue caractérise les débits qui doivent transiter dans les
cours d’eau,
❑ l’hydraulicien en déduit les cotes
atteintes par l’eau dans le cours d’eau
étudié, ainsi que les vitesses
d’écoulement.
11
1. Introduction
• Etape d’une modélisation hydraulique:
– Phase d’identification: Permettant de caractériser comment le
modèle se comporte (calage). Des données d’entrée et des données
de sortie correspondantes, on peut donc caler les différents
paramètres du modèle (les coefficients de Strickler)
– Utilisation en prédiction: Connaissant par exemple la géométrie
d’un cours d’eau, le débit en tête et les caractéristiques hydrauliques
(pertes de charge, prélèvements), on peut alors prédire le niveau
atteint par l’eau, ainsi que les vitesses d’écoulement.
– Utilisation en détection: Connaissant les lois de comportement du
système et les effets observés, on peut remonter aux causes
génératrices des conditions.
12
1. Introduction
• Modélisation mono(uni)dimensionnelle (1-D) :
– Considération des variables d’écoulement selon un axe
longitudinal.
– Résolution des équations de Bernoulli généralisées (prenant
en compte les pertes de charge) en régime permanent, et de
Barré de Saint-Venant en régime transitoire.
14
1. Introduction (Applications des modèles hydrauliques)
• La modélisation hydraulique apparaît partout où un écoulement
varie de façon spatiale ou temporelle.
• Multiples implications et application:
– études d’écoulements de rivières, protection contre les crues,
– transport solide, gestion de systèmes d’irrigation,
– transferts de polluants,
– assainissement, collecte d’eaux pluviales…
• Le principe de base des modèles d’hydraulique 1-D actuels est de
fournir au chercheur, à l’ingénieur, des outils solides dans le but
d'améliorer la compréhension des phénomènes, le
dimensionnement et la gestion des ouvrages hydrauliques.
15
1. Introduction (Applications des modèles hydrauliques)
• Gestion de canaux d’irrigation
– Des projets d’irrigation nécessaires pour répondre à l’augmentation de la
demande en nourriture résultant de la croissance de la population alors
que, dans le même temps, la compétition pour l’eau devient plus serrée.
26
3. Calage des modèles hydrauliques
• Le calage de ces modèles consiste en pratique à établir la valeur:
- de la rugosité des lits mineurs et majeurs (coefficients de Strickler) et
- des échanges entre les différents lits (bras) et entre le lit et les casiers
d'inondation,
- les paramètres liés aux ouvrages singuliers : loi de seuil pour des
vannages ou des déversoirs, loi de perte de charge pour les ponts et
autres ouvrages de franchissement.
• Cette quantification doit ensuite être analysée à l'aune des incertitudes sur les
données observées (débits, hauteurs) auxquelles on a comparé les résultats du
modèle au cours du calage.
31
6. UNE FORTE DIVERSITE DE LOGICIELS PROPOSES
Exemples de modèles hydrauliques
• LISFLOOD-FP qui modélise l’écoulement dans le lit
mineur avec un modèle d’onde de crues 1D et dans le
lit majeur avec l’équation de Manning.
36
6. UNE FORTE DIVERSITE DE LOGICIELS PROPOSES
• Les critères de choix d'un modèle adapte
– Du fait de la très forte diversité des modèles proposés, il n'est pas facile
pour un utilisateur de discerner le modèle qui conviendra le mieux au type
d'utilisation auquel il le destine.
– Avant de choisir un modèle, l'utilisateur doit être capable de définir
précisément le type de problèmes qu'il se propose de résoudre.
– En fonction du type de système à modéliser et des objectifs de l’étude,
l’utilisateur devra chercher l’outil générique capable de modéliser les
phénomènes à étudier.
37
38
7. CONTEXTES D'UTILISATION DE MODELES D'HYDRAULIQUE 1 -D
• UTILISATION POUR LA RECHERCHE
– Le but premier de ces modèles est rarement la commercialisation.
– outils de recherche permettant d’acquérir une connaissance en profondeur
du comportement hydraulique des ensembles étudiés.
– HEC-RAS : L’US Army Corps of Engineers
C’est actuellement le logiciel 1 -D le plus utilisé au monde (disponible en accès
libre sur Internet, avec un guide d’utilisation complet, un forum de discussion
permettant aux utilisateurs de poser des questions, et d’échanger des
informations sur l’utilisation du logiciel.
– WOLF1D: modèle d’écoulement 1 -D en réseau hydrographique quelconque
– SIC (Simulation de Canaux d’Irrigation): Cemagef de Montpellier
39
7. CONTEXTES D'UTILISATION DE MODELES D'HYDRAULIQUE 1 -D
• UTILISATION PAR LES BUREAUX D'ETUDE
– Les logiciels Carima et Canoe, développés par le bureau d’étude français
SOGREAH
– SIC pour la régulation des canaux d'irrigation,
– CARIMA pour la modélisation de plaines d'inondation et de rivières à forte
influence aval,
– HEC RAS pour les régimes mixtes (fluvial/torrentiel) en rivière
40
7. CONTEXTES D'UTILISATION DE MODELES D'HYDRAULIQUE 1 -D
• UTILISATION DANS L'ENSEIGNEMENT
Le plus important n'est pas l'utilisation du logiciel lui-même, mais
l'apprentissage de la démarche générale de modélisation, depuis la collecte et
l’entrée de données dans le logiciel, qui représentent souvent une part
importante du travail de modélisation, jusqu'au calcul de tracé des lignes d'eau.
Les critères de choix des logiciels utilisés sont donc la simplicité, ainsi que
l'aspect pédagogique
– Canal 9, développé par le Cemagref, offre ainsi l'avantage de présenter
clairement les calculs réalisés pas à pas, et propose des interfaces
graphiques didactiques.
– SIC du Cemagref peut être également utilisé pour d'autres modélisations
plus élaborées (propagations de crues, écoulements de rivières…)
– HEC RAS de l'Usace peut etre utilisé pour des études d'aménagements de
rivières en régime permanent
41
7. CONTEXTES D'UTILISATION DE MODELES D'HYDRAULIQUE 1 -D
• UTILISATION DANS L'ENSEIGNEMENT
– HEC RAS de l'Usace peut être utilisé pour des études d'aménagements de
rivières en régime permanent
– Possibilité de couplage de HEC-RAS au logiciel de SIG ArcView via
l'extension HECGeoRAS, ce qui permet/
• d'extraire des profils en travers du modèle numérique de terrain afin de
construire le modèle hydraulique,
• après calcul des lignes d’eau, de visualiser sous SIG l’extension d’une
inondation
• de produire des vues en plan d'inondation à partir de la ligne d'eau calculée
par HEC-RAS.
– Le choix de ce logiciel pour l'enseignement s'appuie sur les critères de
rapidité, fiabilité et facilité d'utilisation, ainsi que ses possibilités variées de
visualisations.
42
43
Présentation de HEC RAS
8. Présentation de HEC RAS
3.1 HEC-RAS
57
Modélisation hydraulique (Application)
Application 1:
- Section du canal supposée trapézoïdale et l’écoulement est supposé
permanant.
- Longueur du canal : 800m
- Nombre initial de cross sections (profile de la section ou profile en travers): 9
et uniformément repartie à 100 m de distance le long du canal
- Coefficient de Manning sur les berges (gauche et droite) du canal est de 0.042
et dans le canal est de 0.035
- Débits s’écoulant dans le canal : 65m3/s
- L’amont du canal est plus élevé que l’aval de 4m
58
Modélisation hydraulique (Application)
Application 1:
Profile de la section trapézoïdale
Stations Elevation
0 7
30 6
50 4
70 4
90 6
120 7
59
Modélisation hydraulique (Application)
Application 1: Etapes
- Créer un nouveau projet (Canal_Trapezoidal) dans un dossier nommé
‘Traperzoidal’ a new project in a folder (File>New project> donner le nom du
projet tout en sélectionnant le dossier qui doit le contenir)
- Utiliser l’unité de mesure internationale (Options> Unit system> System
International >OK>Yes)
- Fournir les infos sur la géométrie du canal
Edit>Geometric data
- Cliquer sur ‘River Reach’ pour tracer la rivière ou le canal. Apres avoir tracé,
double cliquer puis enregistrer sous un nom (Riviere, Principale)
- Cliquer sur ‘Cross Section’ pour entrer les informations sur les profiles en travers
(profiles des sections) du canal
60
Modélisation hydraulique (Application)
Application 1: Etapes
Edit>Geometric data
- Cliquer sur ‘Cross Section’ pour entrer les informations sur les profiles en
travers (profiles des sections) du canal
- Option> ‘Add a new cross-section’ pour activer un point d’entrer des données de profile en
travers d’une station > entrer le numéro 1 correspondant à la première station la plus en aval
du canal (point vers ou les eaux s’écoulent)
- Fournir les données de ‘cross section coordinates’ à partir du tableau précédant. Les
informations doivent être relevées de la gauche vers la droite du canal. Station=> Abscisses
cumulatives
- Fournir les informations sur les coefficients de Mannings ‘Mannings Values’
- ‘Main Channel Bank stations’: Left Bank: Abscisse de la première station la plus à gauche du
canal (0); Right Bank: Abscisse de la dernière station la plus à droite du canal (120)
- Cliquer sur ‘Apply Data’, la figure du canal apparaitra à droite. 61
Modélisation hydraulique (Application)
Application 1: Etapes
Edit>Geometric data
- Cliquer sur ‘Cross Section’ pour entrer les informations sur les profiles en
travers (profiles des sections) du canal
- Option> ‘Copy current cross-section’ pour dupliquer les info précédemment entrées > ‘River
Sta’ égale 9, la station la plus en amont du canal > OK
- Elévation de la station amont doit être plus élevée que la station avale, d’où la nécessité
d’ajuster l’élévation: Option> Adjust the elevation> 4 pour ajouter 4m de plus que l’amont
- ‘Downstream Reach Length’= 800, 800, 800 correspond à la distance à laquelle se trouve le
prochain cross section. Nous n’avons à fournir que deux ici.
- Cliquer sur ‘Apply Data’ pour enregistrer dans le projet
- Exit
62
Modélisation hydraulique (Application)
Application 1: Etapes
Edit>Geometric data
- Tools>XS interpolation> Within a Reach : pour interpoler les données de cross
section précédemment entrées
- S’il y a plusieurs River ou Reach, en choisir ceux qui seront modifiés
- Upstream Riv Sta: choisir 9; Downstream Riv Sta: choisir 1
- Maximum Distance…: 100
- Decimal places in interpolation: 0.00
- Cliquer sur interpolate XS > Done
- File> Save Geometric Data> pour enregistrer les données sur la géométrie du
canal. Nom possible: Trapesoide_Geometrie > OK
- Les valeurs interpolées sont avec *
- Fermer la fenêtre de Geometric data 63
Modélisation hydraulique (Application)
Application 1: Etapes
Edit>Steady Flow data
- PF= Débit dans le canal> Entrer 65
- Cliquer sur ‘Reach Boundary Conditions’ > Downstream devrait être
automatiquement sélectionné en bleu en bas
- Cliquer sur ‘Normal Depth’: Ajouter le Slope: 0.005 (4m/800m )> OK>OK
- Cliquer sur Apply Data
- Cliquer sur ‘ File> Save flow data’ > Nom: Steady_flow > OK
- Close Steady Flow data
64
Modélisation hydraulique (Application)
Application 1: Etapes
Run>Steady Flow Analysis
- Short ID: Plan01
- Cliquer sur Compute > Close
- File> Save Plan
- Close
65
Modélisation hydraulique (Application)
Application 1: Etapes
View>Cross sections
- Energy Profile: EG PF
- Ground
- Critical Profle (Crt PF )
View> X-Y-Y perspective plots
- Azimuth Angle: 10; Rotation Angle: 60;
View> Rating Curve
- Pour trace la vitesse d’ecoulement à travers une section donnée: run>Steady Flow
Analysis>Option>Flow Distribution Locations> LOB=ROB=Channel=10 en haut comme en bas
>OK > Save Plan> Compute pour excecuter le modele encore.
- Cliquer sur Cross Section> Option>Velocity Distribution> Ticker Volocity Distribution> OK, On voit
apparaitre le profile de distribution des vitesses le long d’une section;
- Cliquer sur l’avant 2e derniere icone pour voir le tableau associé a la simulation
- General Profile Plot of the channel (16e icone)
66
Modélisation hydraulique (Application)
Application 2: Etapes
- On désire implanter un pont sur une branche d’une rivière
où 20 cross-sections ont été mesurées. L’ image dans le
dossier ‘HEC-RAS_Data’ fourni un schéma des opérations.
Nous avons donc à modéliser le comportement de cette
rivière dans HEC-RAS afin d’avoir les caractéristiques
nécessitant le dimensionnement de cet ouvrage sur la
rivière.
-
-
67
Modélisation hydraulique (Application)
Application 2: Etapes
- Comme précédemment, ouvrir un projet, tracer une
rivière puis donner les noms qui conviennent:
- Dans ‘Cross-Section’, Utiliser les données de ‘CS-UCS-
17.xls’ pour remplir Cross Section Coordinate. Ceci
concerne la station la plus en amont du canal
-
-
68
Modélisation hydraulique (Application)
Application 2: Etapes
69
Modélisation hydraulique (Application)
Application 2: Etapes
Cross-Section
- Utiliser les données de ‘CS-DCS-3.xls’ pour remplir Cross Section Coordinate. Ceci
concerne la station la plus en aval du canal (voir figure suivante)
- Choisir la valeur de ‘River Sta’ correspondant à 20 puis
- Option>copy current cross section. ‘River Sta’ serait égal maintenant à 19.
Répéter la même opération 17 fois tout en réduisant le numéro à chaque fois.
- Modifier à nouveau les valeurs de Cross Section Coordinate tout en
respectant l’ordre d’apparition, et le Main Channel Bank Station (left and rigth)
- Remember to click on Apply everytime.
- Close
70
71
Modélisation hydraulique (Application)
Application 2: Etapes
72
Modélisation hydraulique (Application)
Application 2: Etapes
Edit>Geometric Data
- File>Save geometry : Original_geometry
- Interpoler la géométrie par pas de 100m (Tools>XS
interpolation>Within Reach)
- Enregistrer la géométrie comme:
Interpolated_geometry> close
-
- 73
Modélisation hydraulique (Application)
Application 2: Etapes
Edit>Steady Flow Data
- PF (Profile flow rate): 5500 (Periode de retour de 100ans)
- Add flow change condition: pour plusieurs rivières qui se rencontrent. Ce n’est
pas le cas ici.
- Option>Edit profile name> 100ans
- Cliquer sur Reach Boundary Condition> Known WS> 7.85m>OK>OK
- File > Save the Flow Data> nom:100ansInondation> OK
- Apply > Close
-
74
Modélisation hydraulique (Application)
Application 2: Etapes
Run>Steady Flow Analysis
- Sélectionner la géométrie pour l’analyse
- Short ID: 01
- Option> Save Plan> Plan
- Compute pour faire tourner le modèle
- Close
75
Modélisation hydraulique (Application)
Application 2: Etapes
Analyse des resultats:
- Geometric Data> Cross Section pour voir les niveaux dans les différentes sections du
canal en modifiant River Sta
- Close Geometric Data and Cross Section
- Cliquer sur l’avant 2e derniere icone pour voir le tableau associé a la simulation
- General Profile Plot of the channel (16e icone)
- Pour trace la vitesse d’ecoulement à travers une section donnée: run>Steady Flow
Analysis>Option>Flow Distribution Locations> LOB=ROB=Channel=10 en haut
comme en bas >OK > Save Plan> Compute pour excecuter le modele encore.
- Fermer les deux fenetres
-
76
Modélisation hydraulique (Application)
Application 2: Etapes
Analyse des resultats:
- Cliquer sur Cross Section> Option>Velocity Distribution> Ticker Volocity
Distribution> OK, On voit apparaitre le profile de distribution des vitesses le
long d’une section;
- Exporter les figures: File>Copy Plot to Clipboard; on peut coller dans word ou
paint puis fermer le cross section fenetre
- View> X-Y-Z Perspective Plots (Classic)>
77
78
Modélisation hydraulique (Application)
Application 3:
Se rappeler pour changer le format de la date en "Anglais (Etat Unis)" dans les paramètres régionnaux.
79
Application 3
Différentes étapes énumérées ci-après :
- Create a new project (File>New Project> chose the
containing folder and name the project)
- Set unit to international system (Option>Unit system>
System International)
Open RAS Mapper (GIS Tools)
- Right click on Terrain and Create a new terrain > Select no for
projection > click + > add the raster file> select Yes> Create
terrain>Close
- Double click on terrain> Zoom to layer to show the map> save
80
2. Modélisation hydraulique (Application)
Edit>Geometric data
- right click to set the geometric extent
- Click on 2D flow area, draw the extent of the area
- Double click to stop and name the file: 2DflowArea
- Click on the map and select “Edit 2D flow Area”
- Select “generate computational points on regular…”
- Provide dx and dy values > select “Generate points in 2D…”
- Select ‘force mesh recomposition’
- Save the file as “Geometric”
81
2. Modélisation hydraulique (Application)
Edit>Geometric data
- If the DEM is lost, Go to RASmapper, save the terrain map and add the
saved file through “Map Layers” and close RASmapper
- In ‘Geometric data’, Add the created map by selecting layer to view in
background (it is the 10th icon right, just before description)
Defining the boundary condition ( Clip on SA/2D Area BC)
- Clip on SA/2D Area BC
- Draw the line on the upper BC and save it
- Draw the line on the Dowstream BC and save it
- Click on the map and select “Edit 2D flow Area” then ‘force mesh…’
- Save the Geometric data and close the window
82
2. Modélisation hydraulique (Application)
Edit>Unsteady flow data
- Select Upper stream> Click on flow hydrograph> Copy
the hydrograph and paste in the column name “Flow
- Select ‘fix start time’ and provide date and Add start
time >Time:00:00>Data interval: 1 day
- Provide the Energy slope ‘EG slope…’ (0.05)> Click OK
- Select now downstream> Normal depth> 0.01
- Save the data as UnsteadyFlowData
- Initial conditions are optional
83
2. Modélisation hydraulique (Application)
Run>Unsteady flow analysis
- Save a plan as PlanA > ID=1
- Programs to run: select the four (Geometric
preprocessing , Unsteady flow, … , Flood plain
- Simulation time in line with the data provided, Time:
00:00
- C omputational interval: 6h for all for example
- Save the plan and close the window
84
2. Modélisation hydraulique (Application)
RASMapper
- Geometry>2D flow Area (Right click) > compute 2D flow
Area Hydraulic table
- Save the file and close the window
85
2. Modélisation hydraulique (Application)
Run>Unsteady flow analysis
- Options>Calculation Options>2D Flow Options>
- 6. Equation set> Select following the preference
(Diffusive Wave)
- 7. Initial condition (Warm Up): 48h
- Click OK > Save the plan> Compute
- Close after computation
86
2. Modélisation hydraulique (Application)
RASMapper to view the results
87