REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

MINISTERE DE L’ENSEIGNEMENT SUPERIEUR ET

DE LA RECHERCHE SCIENTIFIQUE
UNIVERSITE MOHAMED BOUDIAF - M’SILA

FACULTE DE TECHNOLOGIE DOMAINE: SCIENCE ET TECHNIQUE

DEPARTEMENT D’HYDRAULIQUE FILIERE : HYDRAULIQUE
N° : ……………………………………….. OPTION :RESOURCES HYDRAULIQUE

Mémoire présenté pour l’obtention

Du diplôme de Master en Resources Hydraulique
Par :
 CHEBABHI RADHWANE
 DERADJI MOUSSA

Intitulédu mémoire :

Etude hydrologique appliquée pour le dimensionnement

des ouvrages hydrauliques sur oued. Cas d'un ouvrage sur
oued sahel

Dirigé par : Mr : SAHLI.Y

Année universitaire : 2019 /2020

 Remerciements

Au moment de clôturer ce mémoire, je remercie en premier lieu dieu de m’avoir donné la

force et la patience de terminer ce travail, je souhaiterais apporter aussi ma vive
reconnaissance à certaines personnes qui m’ont aidé et encouragé au cours de ces deux
dernières années de Master.

Je voudrais remercier très chaleureusement Mr Youcef Sahli, mon promoteur, pour le temps
qu’il a sue me consacrer et les précieux conseils qu’il a pue me donner.

Je ne saurais oublier ma mère, mon père, mon frère et ma sœur.

Enfin, je tiens remercier tous les personnels techniques et administratifs de département

d’hydraulique.
 Didicase
C’est avec un immense plaisir que je dédie ce modeste travail à :

A ma chère belle mère

Aux supports de mon âme, vous mes très chères frères : aed el
baset, Badro et mes sœurs.

Et a mes amies : zakaria dahmani ,badreddine aouina ,bilal mezidi,

, abdelbassit mekafes ,yaakoub hamoma,badre lahmar,basem
haouchat ,
A Tous les enseignants qui m’ont aidé de proche ou de loin pour être
un jour un master hydraulique

A tous les étudiants master hydraulique de l’université de M’sila et

surtout les étudiants de deuxième Année Master « Resources
Hydraulique » promo 2019/2020

Je vous remercie tous la famille chebabhi

 Sommaire
Introduction Générale .................................................................................................................... 2
Chapitre 1 : Présentation de la zone d’étude
1-1- Définitions du bassin versan .................................................................................................................... 4
1-2- ituation géographique.......................................................................... ...................................................... 5
1-3- Conditions climatiques.............................................................................................................................. 7
1-3-1- Température ................................................................................................................. .................. 7
1-4- Le couvert végétal...................................................................................................................................... 7
1-5- La géologie................................................................. ..................................................................... ............8
1-6- Les sols.............................................................................. .......................................................................... 9
1-7- Les Vents................................................................................................... ................................................10
1-8- Précipitations....................................... ..................................................................................................... 11
1-9- L’Humidité Relative................................................................................................................... ..............11
1-10- Morphologie du bassin...................................................................................................................... ...... 11
1-11- Le réseau hydrographique......................................................................... .......................................... .... 11

Chapitre 2 : Etude hydrologique

2-1- Introduction................................................................................................................................................ 14
2-2- Données topographique............................................................................................................................ 14
2-3- Caractéristiques morphologiques du bassin versant d’Oued Sahel .................................................14
2-4- Détermination de l’intensité moyenne des precipitations. .............................................................. 15
2-5- Estimation du débit de crue ..................................... ............................................................................... 18
2-5-1- Formule Rationnelle ......................................................................................................................... 18
2-5-2- Formule Crupédix ........................................................................ ..................................................... 19
2-5-3- Formule de Transition ....................................................................................................................... 19
2-5-4- Formule de Turazza ......................................................................................................... ..................20
2-6- Temps de concentration.......................................................................................................................... 20
2-7- Coefficient de ruissellement « C »......................................................................................................... 21
2-8- Application avec Turazza .................................................................................................................... ...23
2-9- Conclusion............................................................................................................................................... 24

Chapitre 3 : Etude Hydraulique

3-1- Introduction............................................................................................................................... ................26
3-2- Débit d’apport et débit capable de l’ouvrage hydraulique.................................................................... ..26
3-2-1- Estimation du débit de crue (rappel) ........................................................................................... 26
3-2-2- Choix du coefficient de ruissellement........................................................................................... 28
3-2-3- Choix du période de retou............................................................................................................... 28
3-3- Dimensionnement de l’ouvrage hydraulique....................................................................................... 29
3-3-1- Calcul du débit de saturation......................................................................................................... 29
3-3-2- Conception des ouvrages hydrauliques........................................................................................ 30
3-3-3- Protections des ouvrages hydrauliques........................................................................................ 30
3-4- méthode de vérification hydraulique..................................................................................................... 31
3-4-1- méthode de vérification hydraulique............................................................................................ .31
3-5- Application pour notre ouvrage (Oued SAHEL)............................................. .................................. 33
3-6- Conclusion ............................................................................................................................................... 35

Chapitre 4 : Présentation du logiciel HEC-RAS

Introduction................................................................................................................................................. .......37
4-1- Principe d’utilisation de logiciel HEC-RAS ................................................................................ ......37
4-2- Etapes de la modélisatio..................................................... ................................................................ .. 38
4-2-1- Projet ........................................................................... .......................................................... ...........38
4-2-2- Géométrie des sections............................................... ....................................................................39
4-2-3- Débit et conditions aux limites.................................................................................................. .. 43
4-2-4- Simulation hydraulique ............................................. ............................................................... .....45
4-2-5- Visualisation des resultants...................................................................................................... .....45
4-3- Théorie de base de calcul du modèle HEC-RAS........................................................................ ....... 47
4-3-1- Les équations de Saint-Venant 1D ....................................................................................... .......48
4-3-2- Le cas stationnaire (SteadyFlow) ................................................................................................48
4-3-3- Le cas instationnaire (Unsteady Flow)............................................................................. .......... 50
4-4- Conclusion............................................................................................................................................. ..51

Chapitre 5 : Calcul de la PHE (Plus Haute Eau) par HEC-RAS

5.1 Introduction ............................................................................................................................................... 53
5.2 Le logiciel Hec Ras.................................................................................................................................. 53
5.3 Les étapes de développement du model HEC-RAS.......................................................................53
Conclusion générale
Bibliographie
 LISTE DES Tableau
Numéro de Tableau de titre page
la Tableau
Tableau 1.1 Superficie des sous bassins versants de l'oued sahel(A.N.R.H., 1979, 1980) 7
Tableau 1.2 Moyennes mensuelles des températures (max ‘M’ et min ‘m’) dans la 7
région de la basse vallée de la Soummam (période de 1970 - 2011)
Tableau 1.3 Vitesse mensuelle moyenne du vent (1994-2003) 10
Tableau 1.4 Précipitations moyennes annuelles et mensuelles en mm 11
Tableau 1.5 Humidités relatives mensuelles en %(1994-2004) 11
Tableau 2.1 caractéristiques morphologiques du bassin versant 14
Tableau 2.2 Implantation et données pluviométriques de la station de TAZMALT 16
Tableau 2.3 Intensités obtenues en[ mm / heure] de la station de TAZMALT 17
Tableau 2.4.1 : Valeur du Coefficient « C1 » pour diverses pentes 22
Tableau 2.4.2 Valeur du Coefficient « C2 » par rapport à la nature du SolTableau 22
Tableau 2.4.3 Valeur du Coefficient « C3 » relative à la couverture du sol 22
Tableau 2.5 Les caractéristiques des paramètres du bassin versant sont 23
Tableau 2.6 Résultats de calcul du débit de pointe de dimensionnement de l’ouvrage 24
hydraulique
Tableau 3.1 (modèle) : Débit d’apport et débit capable de l’ouvrage hydraulique 31
Tableau 3.2 Résultats de calcul 34
Tableau 4.1 Valeurs estimées pour le coefficient de Manning en fonction de la nature 42
des parois
 LISTE DES FIGURES
Numéro de la Figure de titre page
figure
Figure Ι-1 Distinction entre bassin versant réel et bassin versant topographique 5
Figure I.2 Localisation du bassin versant de l'oued sahel 6
Figure I.3 Les sous bassins versants de l'oued sahel 6
Figure I.4 Carte géologique du bassin versant de l'oued sahel 9
Figure I.5 Carte des sols du bassin versant de l'oued sahe 10
Figure I.6 Carte du réseau hydrographique du bassin de l'oued sahel 12
Figure 4-1 La fenêtre principale de HEC-RAS 37
Figure 4-2 Fenêtre New Project 39
Figure 4-3 Fenêtre Geometric Data 40
Figure 4-4 Fenêtre Cross Section Data 41
Figure 4-5 Fenêtre XS interpolation 43
Figure 4-6 Fenêtre Steady Flow Data 44
Figure 4-7 Fenêtre Reach Bondary Conditions 44
Figure 4-8 Fenêtre Steady Flow Analysis 45
Figure 4-9 Fenêtre de profil de surface de l’eau 46
Figure 4-10 Fenêtre de surface de l’eau 3D 46
Figure 4-11 Résume des erreurs et des avertissements 47
Figure 4-12 Conservation de la charge entre deux sections (source: HEC-RAS) 49
Figure 4-13 Découpage de la section pour le calcul du frottement (source: HEC-RAS). 50
Figure 4-14 Volume de contrôle élémentaire de dérivation de la continuité et 50
d’équations de moment (Inflow : débit entré et Outflow : débit sortie)
Introduction Générale
Le comportement hydrologique des bassins versants se traduit parfois par le phénomène
des crues qui peut être la source de certaines catastrophes appelées : inondations.
Ces événements hydrologiques exceptionnels constituent l’un des plus importants risques
naturels responsable parfois de la perte de vies humaines et de biens matériels.
L’étude hydrologique et hydraulique pour l’implantation d’un ouvrage routier traversant
un oued du bassin versant de sahel est justement pour objectif de prendre en charge ce
risque par quelques solutions qui seront dégagées dans notre étude.
L’objectif principal de notre étude est de déterminer la variation du niveau de remonté
des eaux dans la section ou l’ouvrage sera réalisé sur oued Sahel. Autrement dit, l’objectif
final de notre travail est de déterminer la plus haute eau (PHE) sur cette section de l’oued
à fin de la comparer à la hauteur naturelle des berges correspondant à la même section.
Cet objectif de notre travail, sera fait pas deux manières à savoir : l’utilisation des
formules de calcul avec les donnes hydrologique de la zone (étude hydrologique et
hydraulique) et l’utilisation en deuxième lieu le logiciel HEC-RAS avec ARC GIS. Ces
derniers nous permettent la simulation du niveau d’eau de l‘oued du bassin versant de
Sahel après le passage de la crue.
Ce travail sera reparti sur plusieurs chapitres qui sont présenté dans le sommaire.

2
Chapitre 1 : Présentation de la zone d’étude

Introduction:

L'étude détaillée du bassin versant d'un cours d'eau est indispensable avant le
dimensionnement des ouvrages hydraulique situant dans ce bassin versant, ou étude
d'impact sur l’environnemental. L'objet de ce chapitre est :
- Une connaissance générale des principales caractéristiques du bassin versant de l'oued
Sahel. Cette partie sera donc avant tout descriptive et qualitative. Dans un premier temps,
nous allons décrire globalement le bassin versant, en le situant géographiquement, en
étudiant sa topologie et sa nature géologique. Ensuite, nous allons chercher à étudier plus
précisément le cours d'eau qui sera traversé par notre ouvrage.

1-1. Définitions du bassin versant:

Le bassin versant est l'ensemble des territoires superficiels et souterrains qui se drainent
vers le même exutoire.
On distingue habituellement:

- Le bassin fonctionnel : portion de l'espace géographique dont les eaux se

drainent vers un même cours d'eau.
- Le bassin topographique : portion de l'espace géographique dont les
pentes convergent vers un même cours d'eau. Il peut être identique au bassin
fonctionnel ou différent selon la géologie ou la perméabilité du terrain.
Toutefois, la délimitation topographique nécessaire à la détermination en surface du bassin
versant naturel n'est pas suffisante. Lorsqu'un sol perméable recouvre un substratum
imperméable, la division des eaux selon la topographie ne correspond pas toujours à la
ligne de partage effective des eaux souterraines (voir Figure 1-2). Le bassin versant est
alors différent du bassin versant délimité strictement par la topographie. Il est appelé dans
ce cas bassin versant réel.

4
Chapitre 1 : Présentation de la zone d’étude

Figure 1-1: Distinction entre bassin versant réel et bassin versant topographique.

1-2.Situation géographique

Le bassin versant de l'oued sahel est situé dans la partie Nord-est du pays entre 3° 60’ et
4° 70’ de longitude Est et entre 36° 00’ et 36° 50’ de latitude Nord.

Le bassin versant de l'oued sahel s’étendant selon un axe Nord -ouest. Il se compose de
région:
les plateaux de Bouira. Il est limité:

- Au Nord : par les montagnes de la grande Kabylie (massif de Djurdjura).

- Au Sud : par les montagnes de Bibans et de Mansourah, au Sud-est par le contrefort des
monts du Hodna.

- Et à l’Ouest : il est limité par les cours de l’Isser, et du Sébaou. IL présente une forme
très irrégulière.

5
Chapitre 1 : Présentation de la zone d’étude

Figure 1-2. Localisation du bassin versant de l'oued sahel

Le bassin versant de l'oued sahel est subdivisé en 05 sous bassins versants (Fig. 1.3) et il
2
couvre une superficie de 3750 km .

Figure I.3. Les sous bassins versants de l'oued sahel (extrait de la carte hydro-
climatologique de surveillance de la qualité de l’eau dressée par l’ANRH, 2005)

6
 Chapitre 1 : Présentation de la zone d’étude

Tableau 1.1. Superficie des sous bassins versants de l'oued sahel

(A.N.R.H., 1979, 1980)

Sous bassin versant Code Superficie (km²)

Oued Eddous 1501 903.000
Oued Zaiane 1502 1012.00
Oued Sahel amont 1503 430.000
Oued Azerou 1504 1084.00
Oued Sahel Aval 1505 326.000

1.3- Conditions climatiques :

1.3.1-Température :
La température, second facture consécutif du climat, influe sur le développement de la
végétation et le fonctionnement du système radiculaire et elle affecte les relations
physiologiques entre les organes souterrains et les organes aériens, d’après (GRECO, 1966).
Les températures moyennes mensuelles relevées sur 41 ans de 1970 à 2011 sont représentées
par le tableau suivant:
Tableau 1-2 : Moyennes mensuelles des températures (max ‘M’ et min ‘m’) dans la
région de la basse vallée de la Soummam (période de 1970 - 2011)

Mois Jan. Fév. Mars Avr. Mai Juin Juil. Aout Sep. Oct. Nov. Déc.
M(C°) 16.47 17.07 18.61 20.26 22.86 26.37 29.47 30.16 28.14 25.28 20.85 17.57
m(C°) 7.50 7.67 9.04 10.83 13.96 17.57 20.37 21.18 19.04 15.73 11.71 8.70
Tmoy(C°) 11.99 12.37 13.82 15.54 18.41 21.97 24.92 25.67 23.59 20.50 16.28 13.13

m : température moyenne minimales du mois

M : température moyenne maximales du mois
Tmoy : Température moyenne mensuelle

1.4.Le couvert végétal

La végétation du bassin versant de l'oued sahel (Soummam) est de type méditerranéen et
est très variée. Les régions alluvionnaires du haut sahel sont occupées par des chênaies. La
basse et moyenne sahel sont occupés par des oliveraies et des vignobles. Dans la région de la
moyenne sahel, le couvert végétal est moins dense.

Dans la plaine, la végétation est dense mais essentiellement temporaire ; elle est formée
par les grands et formidables champs de cultures maraîchère diverses.

7
 Chapitre 1 : Présentation de la zone d’étude

Les forêts qui couvraient la région fournissaient le bois nécessaire à une industrie du bois
florissante mais malheureusement ce capital tend à disparaître sous les feux qui ravagent
chaque année des milliers d’hectares. A cela s’ajoute l’absence d’une politique claire de
reboisement et de lutte contre les incendies.

1.5.La géologie
On distingue au niveau du bassin versant de l'oued sahel différentes formations
géologiques (Fig. I.4).

- Le Quaternaire : il s’agit des alluvions le long des Oueds et au niveau des hautes plaines.

- Le Pliocène Continental : le coteau du sahel est couronné par des graviers siliceux et
fluviaux. Ailleurs, il est conglomératique sinon caillouteux.

- Le Miocène : il s’agit essentiellement du Miocène inférieur.

- L’Oligocène : il est généralement argilo-gréseux. Il affleure en abondance à partir de

Tazmalt et de part et d’autre de la vallée de la sahel, en allant vers l’ouest. Il est présent au
niveau de la chaîne de Gouraya et dans la région de Béni-Mansour.

- L’Eocène : il est complètement marneux, parfois calcaro-marneux. On le trouve surtout au

niveau des hautes plaines.

- Le Crétacé : se présente, en général, sous forme de flysch. Il est étendu, continue depuis
la chaîne des Bibans jusqu’aux montagnes de la petite Kabylie. Le faciès est marneux au
Crétacé supérieur (Sénonien). Il est schisteux, sous forme de bancs de plusieurs centaines de
mètres d’épaisseur avec quelques intercalations de grès au crétacé moyen (Albien) et
n’apparaît que

Dans la région de Béni-Mansour et la chaîne des Bibans. Il est grésopélitique, parfois

carbonaté-gréseux au crétacé inférieur.

- Le Jurassique : il affleure surtout au niveau de la chaîne calcaire et le long de sa continuité

vers l’Est. Il est principalement calcaro-dolomitique.

- Le trias : il est de faciès gypso-salin et apparaît le long de contacts anormaux (contact des
nappes de charriage) dans le Tell méridional et est présenté par des conglomérats et des
grés rougeâtres dans Tell septentrional.

8
 Chapitre 1 : Présentation de la zone d’étude

Figure 1.4. Carte géologique du bassin versant de l'oued sahel

1.6.Les sols
La majorité des sols dans le bassin de l'oued sahel sont des sols calcaires (riches en
calcaire). Leur complexe absorbant est modifié par l’eau. Ces sols ont généralement une
texture légère et sont donc perméabl
Au nord du bassin, on rencontre des sols insaturés qui ne contiennent pas de calcaire dont
l’argile peut être plus abondante en surface qu’en profondeur ; leurs roches mères sont
généralement imperméables ou donnent des produits de décomposition imperméables (Fig.
I.5). Le long des Oueds, les sols sont des dépôts alluviaux appelés sols alluviaux à profil
indifférencié ou présentent une simple accumulation de matières organique en surface. En
plus on rencontre aussi des sols calciques

9
 Chapitre 1 : Présentation de la zone d’étude

Figure 1.5. Carte des sols du bassin versant de l'oued sahel

1.7.Les Vents :
Dans la vallée de la Soummam, les vents dominats sont les vents du nord-ouest (vents
marins qui s’engouffrent facilement dans la vallée de la Soummam). En hiver, ils apportent
les nuages, la pluie. En été, ce sont les vents sud-est qui sont beaucoup mois importants et
surtout le vent sud, le sirocco, vent sec d’intensité variable ayant une action funeste sur la
végétation.
Les moyennes mensuelles des vitesses du vent moyen sur dix ans sont données par le Tableau
suivant :
Tableau 1-3: Vitesse mensuelle moyenne du vent (1994-2003) [4]
Mois Jan Fév. Mars Avr. Mai Juin Juil. Aout Sep Oct. Nov. Déc.
V(m/s) 4,7 4,28 3,92 3,52 3,58 3,47 3,63 3,38 3,54 4,13 4,4 4,91

La vitesse du vent la plus faible a été enregistrée le mois d’aout de l’ordre de 3,38 m/s
contrairement à la période hivernale où les vents sont froids et violents surtout dans le mois de
décembre et janvier dont la vitesse du vent varie entre 4,91 et 4,7 m/s,

10
 Chapitre 1 : Présentation de la zone d’étude

1.8-Précipitations :

Les traitements statistiques des séries d’observations de ces stations donnent les
résultats suivant :

Tableau 1-4: Précipitations moyennes annuelles et mensuelles en mm

Station
Sep Oct. Nov. Déc. Jan Fév. Mars Avr. Mai Juin Juil. Aout Année
Mois
Akbou 29 40 48 83 74 67 54 56 37 13 5 7 513
Tazmalt 9,8 34,4 39,4 55,7 50,4 49,3 44,8 43,2 19,4 6,2 3,2 9,5 365,3

1.9-L’Humidité Relative
Le tableau suivant illustre les données relatives aux humidités moyennes
mensuelles (moyennes mensuelles et annuelle de l’humidité) :

Tableau 1-5: Humidités relatives mensuelles en %(1994-2004)

Mois Jan Fév. Mars Avr Mai Juin Juil. Aout Sep Oct. Nov. Déc. Moy.
Humidité
74,9 77,0 76,8 78,20 78,5 77,0 72,7 75,5 75,4 74,2 76,30 75,50 76,0
(%)

Les pourcentages d’humidité varient entre 72 et 78 % au cours de l’année, La moyenne

minimale a été enregistrée le mois de juillet avec une valeur de 72,7% et la maximale au mois
de Mai avec une valeur de 78,5 %.
1.10-Morphologie du bassin
Le bassin versant de la Soummam se compose de dix sous bassins correspondant aux
unités hydrologiques certaines sont contrôlées par des stations hydrométriques. Il a une
forme irrégulière, étiré dans la direction Est-Ouest et accède a la méditerrané par le golf de
Bejaia. Le bassin versant de la Soummam est limité par les chaînes du Djurdjura, à l'Est
par les montagnes de la petite Kabylie et au Sud par les montagnes des Bibans et de
Mansourah. Son orographie est très prononcé, aux environs de Bouira jusqu'à Akbou, la
chaîne de Djurdjura n'est a aucun endroit inférieur a 1500 m. C'est dans cette partie que se
situe le plus haut sommet D'Akbou en direction Nord-Est de la mer, la hauteur diminue
graduellement.
Nous constatons ainsi, la présence d'un système de barrières naturelles représentée par
les reliefs. Des chaînes de montagnes qui se dressent face aux vents pluviaux du littoral, il est
en effet un bassin assez bien arrosé

11
 Chapitre 1 : Présentation de la zone d’étude

1.11. Le réseau hydrographique

Le bassin versant de l'oued sahel possède un réseau hydrographique très développé. une
seule rivière et leurs affluents constituent son réseau hydrographique principal : le Sahel,
affluent venant du Nord-Ouest. (Figure 1.6)

Figure 1-6. Carte du réseau hydrographique du bassin de l'oued sahel

Conclusion :
Le bassin versant de Sahel, avec sa superficie de 3750 Km2, et un périmètre de 320 Km, a
une forme allongée, ce qui lui donne un temps de concentration assez long, il est situe
dans la zone Nord-est de l'Algérie qui est considérée comme la région la plus arrosée en
Algérie; donc il est soumise à des fortes probabilités d'inondation

12
 Chapitre 2 : Etude hydrologique

2.1 Introduction
L’étude hydrologique consiste à déterminer le débit de crue de l’oued ou écoulement
qui est destiné à être franchit par un tracé routier. Par conséquent, l’écoulement doit être
rétabli par ouvrage hydraulique (Pont, dalots ou buses) pour notre cas, il s’agit d’un pont.
Les méthodes de calcul des débits de ruissellement utilisent généralement un modèle
statique de la relation “pluie - ruissellement” (Q = f(I)) et leur application présente la
difficulté d’apprécier les caractéristiques physiques et géomorphologiques du bassin versant
tels que forme, pente, superficie, longueur du cours d’eau principal, temps de parcours;
coefficient de ruissellement, etc.

2.2. Données topographiques :

Les supports topographiques qu’on a utilisés dans ces types sont les cartes d’états
majors qui sont scannées et numérisées.
Le bassin d’apports d’eaux pluviales a été délimité en fonction de la structure des talwegs et
des lignes de crêtes sur ces cartes, leurs surfaces sont déterminées à l’aide du logiciel
Autocad. La carte de ce bassin versant est jointe ci après à ce mémoire.
2. 3 Caractéristiques morphologiques du bassin versant d’Oued Sahel :
La détermination des caractéristiques du bassin versant consiste à calculer, ou évaluer, les
paramètres explicatifs, indispensables à l’évaluation des débits de crues. Les paramètres
physiques peuvent être calculés, sans grandes difficultés, avec une précision satisfaisante.

Les paramètres suivants sont évalués et calculés.

- Surface
- Périmètre
- Longueur de l’oued
- Altitude max
- Altitude min
- Longueur du rectangle équivalent
Les caractéristiques des paramètres du bassin versant sont données dans le tableau suivant :
Tableau 2.1 : caractéristiques morphologiques du bassin versant

S P L Lr Hmax Hmin H Ig
Kc TC Ds I(%)
(km2) (km) (km) (km) (m) (m) .moy(m) (m/Km)

3800 291,75 101,87 1,44 111,92 15,39 2305 200 1252,5 10,33 2105 2,07

14
 Chapitre 2 : Etude hydrologique

Avec :
S (km2) : Surface (km2)
P (km) : Périmètre (km)
L (km) : Longueur du cours d'eau principal
Kc : Indice de compacité
Lr (km) : Longueur du rectangle équivalent
TC : Temps de concentration en heures
H.max(m) : Altitude max
H.min(m) : Altitude min
H.moy(m) : Altitude moyenne
Ig (m/Km) : Indice de pente globale
Ds : dénivelée spécifique
I (%) : pente du bassin versant

2.4. Détermination de l’intensité moyenne des précipitations :

o Données pluviométriques
Nous avons recueillis toutes les renseignements pluviométriques disponibles au
niveau de l’Agence Nationale des Ressources Hydrauliques (A.N.R.H.) Bir-Mourad-Rais
(Alger) qui couvrent notre zone d’étude.
Les données hydrologiques existent souvent sous forme de courbes Intensité-Durée-
Fréquence (IDF) au niveau de certaines régions et sous forme d’un listing de pluies
journalières mensuelles annuelles maximales dans d’autres régions.
Les précipitations
Les données pluviométriques se présentent sous forme de pluies journalières
mensuelles annuelles maximales, existants au niveau du bassin versant étudié, concerne
la station de TAZMALT. Cette station est située à l’intérieur du bassin versant concerné
par l’étude dont ses informations sont :

15
 Chapitre 2 : Etude hydrologique

Tableau 2.2 Implantation et données pluviométriques de la station de TAZMALT

Nom de la station Code X Y Altitude Nombre Pjmoy CV

(m) d’années (mm) (mm)
Lambert Lambert

TAZMALT 150501 652.75 344.25 324 39 450 0.4

o Intensités
La détermination de l’intensité de pluie relative à ces stations comprend différentes étapes de
calcul qui sont :

a. Hauteur de la pluie journalière maximale annuelle

Pj 
Pjmoy
C 1
2

 exp u LnCv2  1 
v

Pjmoy : pluie journalière moyenne (mm).

Cv : Coefficient de variation.
U : Variable de Gauss.
Ln : Log. Nepérien.
b. Hauteur de la pluie de durée ‘t = tc’
b
t 
P(t )  1.13. p j .  c 
 24 
Pj : Hauteur de la pluie journalière maximale (mm).
b : Exposant climatique.
tc : Temps de concentration (heure).
c. L’intensité horaire
Pt 
i
tc
i : Intensité de la pluie (mm/h).
tc : Temps de concentration (heure).
P(t) : Hauteur de la pluie de durée tc (mm).

16
 Chapitre 2 : Etude hydrologique

 Courbe IDF (Intensité – Durée – Fréquence)

o Pluie moyenne dans la région :
La pluie annuelle moyenne dans la région est comprise entre 400 et 500mm d’ après
la carte des pluies de l’Algérie du nord de l’ANRH, nous prenons en considération la valeur
de :
Pmoy = 450 mm.
o Pluies de courtes durées :
La relation Intensité - Durée - Fréquence recherchée est de la forme :
I(t,T) = a(T)/tb
I(t,T) : Intensité de pluie en mm/h
T : période de retour
T : temps en heures.
a et b : paramètres de Montana.
b = 0,65

I(t,T) = a (T)/t0.7

T 2 5 10 20 50 100

a 11,7 16,2 19,2 22,0 25,7 28,4

o Intensités obtenues en[ mm / heure] de la station de TAZMALT

I(t,T) 2 5 10 20 50 100
15 mn 28,9 39,9 47,2 54,1 63,2 70,0
30 mn 18,4 25,4 30,1 34,5 40,3 44,6
1h 11,7 16,2 19,2 22,0 25,7 28,4
2h 7,5 10,3 12,2 14,0 16,4 18,1
3h 5,7 7,9 9,4 10,8 12,6 13,9
6h 3,7 5,1 6,0 6,9 8,0 8,9
12h 2,3 3,2 3,8 4,4 5,1 5,7
24h 1,5 2,1 2,4 2,8 3,3 3,6

17
 Chapitre 2 : Etude hydrologique
Courbes Intensité-Durée-Fréquence / Station de TAZMALT

1000.0

100.0
Intensité en mm/h

10.0

1.0
1 10 100
Durée en heures

Série1 Série2 Periode

Série3 de Retour
Série4 Série5 Série6

2.5. Estimation du débit de crue :

Le débit de projet correspond au débit de pointe pour une période de retour donnée,
dimensionnant l’ouvrage hydraulique. Ce débit est calculé par les plusieurs méthodes
méthodes (en fonction de la surface du bassin versant et sa situation géographique):
Rationnelle, Crupédix, transition et Turazza.

Pour notre bassin versant d’une surface qui dépasse largement les 200 km², nous avons
opté pour la formule de Turazza.

2.5.1 Formule Rationnelle : Pour les bassins dont la surface ne dépasse pas 10km2, nous
adoptons cette formule :

Q(T )  K .C(T ) .I (T ) . A

Dont : Q(T) : Débit maximum d’eau pluviale de période de retour (m3/s).

C(T) : Coefficient de ruissellement pondéré pour la période de retour T.

I(T) : Intensité moyenne de pluie en mm/h pour la période de retour T
pendant
le temps de concentration tc.

18
 Chapitre 2 : Etude hydrologique

K : Coefficient de conversion des unités.

A : Aire du bassin d’apport (km²).

2.5.2 Formule Crupédix : Pour les bassins versants de surfaces supérieures à 50km2, la
formule de Crupédix sera retenue.
2
P 
Q10  R .  10  . A0.8
 80 
(2.6)

Dont : Q10 : Débit décennal, en (m3/s).

R: Coefficient régional traduisant l’aptitude de ruissellement :

R=0.2 pour des terrains perméables
R=1.5 à 1.8 pour des terrains imperméables
R=1 pour des terrains intermédiaires
P10 : Pluie journalière de période de retour décennale en (mm).
A : Aire du bassin d’apport (km²).
Le débit centennal est obtenu à partir de la corrélation :
 
Q100  b . Q10
à priori : 1.4  b  4


Jusqu’à 20km2, b est déterminé en appliquant la formule rationnelle (calcul de Q10 et Q100
comme si la formule rationnelle était applicable).


Au-delà de 20km2, b est déterminé à partir des données provenant des cours d’eau jaugés


sur des bassins versants représentatifs à proximité du projet. A défaut, b =2 au minimum.

2.5.3 Formule de Transition : Entre les deux méthodes précédentes, une formule de
transition est appliquée pour des surfaces comprises entre 10km2 et 50km2.

 50  A   A  10 
QTransition     QRationnelle     QCrupédix
 40   40 

19
 Chapitre 2 : Etude hydrologique

2.5.4 Formule de Turazza :

Pour les bassins versants de superficie supérieure à 200 km2 les temps de concentrations sont
supérieure à 12 heures, la pluie continue de durée égale au temps de concentration présente
une probabilité de retour quasi-nulle.
Le ruissellement est provoqué par une succession d’averses de plus en moins grande intensité.
Les débits de crues sont alors fonction des rythmes de succession des averses continues.
Le calcul de débit est obtenu en majorant ; par un coefficient variant de 1.5 à 4, le débit
moyen de la crue définie par la hauteur d’eau ruisselée sur le bassin versant, divisée par le
temps dit de concentration Tc.

Q : débit du bassin versant (m3/s) ;

P (t) : pluie tombée pour une durée égale au temps de concentration du bassin versant en
(mm) ;
Cr : coefficient de ruissellement ;
A : superficie du bassin versant en (km2).
Ka : coefficient d’abattement égale du rapport de pluie moyenne sur la pluie maximale et
varient de 0.5 à 0.9 ;
Tc : temps de concentration détermine par la formule de GIONDOTTI en heures.
Temps de concentration
Pour calculer le débit maximal (débit de projet), le temps de concentration (tc) est un facteur
prépondérant, ainsi la durée de pluie utilisée sera égale à (tc) du bassin. Le temps de
concentration du bassin versant est défini comme le temps nécessaire à une particule d’eau de
parcourir le plus éloigné du bassin jusqu’à l’exutoire de ce bassin. Différentes types de
formules expérimentales sont utilisées pour calculer (tc), tel que les formules de Ventura,
Passini, Giandotti, Pezolli, Horton, Kirpich (1940) et Setra, etc.

 Pour S < 5 km2, on applique la formule de Ventura

S
t C10  0,127
P
Avec : tc10 : Temps de concentration en heure.
S : Surface du bassin versant en km2.

20
 Chapitre 2 : Etude hydrologique

P : Pente moyenne du bassin versant en m/m.

 Pour 5 < S < 25 km2, on applique la formule de Passini

3
Sx L
t C10  0,11
P

Avec : L : Longueur du bassin versant en km.

S : Surface du bassin versant en km2.
P : Pente moyenne du bassin versant en m/m.
 Pour S 25 km , il est calculé à l’aide de la formule de Giandotti
2

4 S  1,5 L
t C10 
0,8 H

Avec : S : Surface du bassin versant en km2.

L : Longueur du bassin versant en km.
H : Altitude moyenne du bassin versant.
Pour une période de retour supérieur à 10ans : le temps de concentration est déterminé par
l’expression suivante :
0.23
 PT   P0 
tcT  tc10   
 10 0 
P _ P

Avec : Tct : Temps de concentration pour une période de retour T en mn

Tc10 : Temps de concentration décennal en mn.
P(T): Pluie journalière de période de retour T en mm.
P10 : Pluie journalière décennale, en mm.
P0 : Rétention initial, en mm.
2.6 Coefficient de ruissellement « C »
Le coefficient de ruissellement « C10 », qui est un indicatif du pourcentage de l’eau ruisselée
sur un sol saturé par les précipitations antérieures par rapport à la totalité de l’averse, sera
estimé à partir de la méthode de Kenessey (Hongrie) où "C10" est la somme de trois
coefficients partiels :
C10  C1  C 2  C 3

C1 : dépend de la pente "p" du thalweg principal estimée comme étant égal au rapport de la
dénivelée du talweg (falaises exclues) par sa longueur développée. Ces informations sont

21
 Chapitre 2 : Etude hydrologique

obtenues généralement à partir des cartes topographiques disponibles. Les valeurs de « C1 »

sont présentées dans le Tableau 2.4.1 ci-dessous.
C2 : dépend de la perméabilité du sol, déterminée à partir des informations concernant la
nature du sol et ses caractéristiques. Les valeurs du coefficient « C2 » sont présentées dans le
Tableau 2.4.2 ci-après.

C3 : dépend de la couverture végétale du bassin versant. Ces informations seront obtenues à partir d’une
inspection visuelle des bassins versants, des photos prises du site et des photos aériennes. Les valeurs du
coefficient « C3 » sont présentées dans le Tableau 2.4.3 ci-après

Tableau 2.4.1 : Valeur du Coefficient « C1 » pour diverses pentes

Pente (P) (%) Coefficient (C1)

P ≤ 3.5 0,01 – 0,05
3,5 ≤ P ≤ 11 0,06 – 0,10
11 < P ≤ 35 0,12 – 0,20
35 < P 0,22 – 0,30

Tableau 2.4.2 : Valeur du Coefficient « C2 » par rapport à la nature du Sol

Nature du Sol Coefficient (C2)

Imperméable 0,22 – 0,30
Peu perméable 0,10 – 0,20
Perméable 0,06 – 0,10
Très perméable 0,06 – 0,05

Tableau 2.4.3 : Valeur du Coefficient « C3 » relative à la couverture du sol

Couverture Végétale Coefficient (C3)

Rocheux 0,22 – 0,30
Prairie 0,17 – 0,25
Labours – Champs 0,07 – 0,15
Forêts et Terrains Sableux 0,03 – 0,05
Remarque :
Après analyse de la topographie du site, le couvert végétal ainsi que la pente moyenne du
terrain, le coefficient de ruissellement C = C1 + C2 + C3 = 0.20 + 0.25 + 0.25 = 0.70
C = 0.7

22
 Chapitre 2 : Etude hydrologique

Pour une période de retour supérieur à 10ans : le coefficient de ruissellement est déterminé
par l’expression suivante :

 p 
CT  0.8  1  0  
 P 
 T  

Avec : P(T): Pluie journalière de période de retour T en mm.

P0 : Rétention initial, en mm, donnée par les formules suivantes :

 C 
P0  1  10   P10
Pour C100.8,  0.8 

Pour C100.8 ; C(T)=C(10)

2.7 Application avec Turazza :
Application de la formule de Turazza pour le bassin versant d’oued Sahel (calcul du débit max à
son exutoire)
Caractéristiques morphologiques du bassin versant de Oued Sahel :
La détermination des caractéristiques du bassin versant consiste à calculer, ou évaluer, les
paramètres explicatifs, indispensables à l’évaluation des débits de crues. Les paramètres
physiques peuvent être calculés, sans grandes difficultés, avec une précision satisfaisante.

Les paramètres suivants sont évalués et calculés .

- Surface
- Périmètre
- Longueur de l’oued
- Altitude max
- Altitude min
- Longueur du rectangle équivalent
Tableau 2.5 : Les caractéristiques des paramètres du bassin versant sont données par le
tableau suivant :

S P L Lr Hmax Hmin H.moy Ig I

(km2) (km) (km) Kc (km) TC (m) (m) (m) (m/Km) Ds (%)
3800 291,75 101,87 1,44 111,92 15,39 2305 200 1252,5 10,33 2105 2,07

Le débit est calculé pour une période de retour de 100 ans.

23
 Chapitre 2 : Etude hydrologique

Tableau 2.2 : Résultats de calcul du débit de pointe de dimensionnement de l’ouvrage

hydraulique

S Hmin Hmax L I Q
Cr
(km2) (m) (m) (km) (%) (m3/s)
Q100 =
3800 200 2305 101,87 0.7 2,07
1350

2.8 Conclusion :

La méthode d’étude adoptée est française. Elle est utilisée pour les zones
méditerranéennes.
Pour utiliser ces formules en dehors de la zone méditerranéenne, il suffit de changer les
indices climatiques des stations utilisées qui change par apport à celles des zones
méditerranéennes.

24
Chapitre 3 : Etude hydraulique

3.1 Introduction :
L’étude hydraulique consiste à dimensionner les ouvrages hydrauliques nécessaires pour
faire traverser les écoulements et les oueds en dehors des zones urbaines et dimensionner
les ouvrages de drainage en zones urbaines après avoir déterminé le débit de crue de chaque
oued ou écoulement consterné.
L’objectif principal de cette partie d’études est de déterminer la plus haute eau (PHE) d’un
oued qui sera traversé par un ouvrage qui est un pont au niveau de oued Sahel, en utilisant
les formules empiriques qui tiennent compte des donnes hydrologiques de la région et en
utilisant les résultats de l’étude hydrologique (calcul du débit du bassin versant).
3.2 Débit d’apport et débit capable de l’ouvrage hydraulique

3.2.1Estimation du débit de crue (rappel) :

Le débit de projet correspond au débit de pointe pour une période de retour donnée,
dimensionnant l’ouvrage hydraulique. Ce débit est calculé par les plusieurs méthodes
méthodes (en fonction de la surface du bassin versant et sa situation géographique):
Rationnelle, Crupédix, transition et Turazza.

Pour notre bassin versant d’une surface qui dépasse largement les 200 km², nous avons opté
pour la formule de Turazza expliqué en chapitre 2.

a. Formule Rationnelle : Pour les bassins dont la surface ne dépasse pas 10km2, nous
adoptons cette formule :

Q(T )  K .C(T ) .I (T ) . A

Dont : Q(T) : Débit maximum d’eau pluviale de période de retour (m3/s).

C(T) : Coefficient de ruissellement pondéré pour la période de retour T.

I(T) : Intensité moyenne de pluie en mm/h pour la période de retour T
pendant
le temps de concentration tc.
K : Coefficient de conversion des unités.
A : Aire du bassin d’apport (km²).
b. Formule Crupédix : Pour les bassins versants de surfaces supérieures à 50km2, la
formule de Crupédix sera retenue.

2
P 
Q10  R .  10  . A0.8
 80 

26
Chapitre 3 : Etude hydraulique

Dont : Q10 : Débit décennal, en (m3/s).

R: Coefficient régional traduisant l’aptitude de ruissellement :

R=0.2 pour des terrains perméables
R=1.5 à 1.8 pour des terrains imperméables
R=1 pour des terrains intermédiaires
P10 : Pluie journalière de période de retour décennale en (mm).
A : Aire du bassin d’apport (km²).
Le débit centennal est obtenu à partir de la corrélation :
 
Q100  b . Q10 à priori : 1.4  b  4

Jusqu’à 20km2, b est déterminé en appliquant la formule rationnelle (calcul de Q10 et Q100
comme si la formule rationnelle était applicable).

Au-delà de 20km2, b est déterminé à partir des données provenant des cours d’eau jaugés sur

des bassins versants représentatifs à proximité du projet. A défaut, b =2 au minimum.
c. Formule de Transition : Entre les deux méthodes précédentes, une formule de
transition est appliquée pour des surfaces comprises entre 10km2 et 50km2.

 50  A   A  10 
QTransition     QRationnelle     QCrupédix
 40   40 
d. Formule de Turazza :

Pour les bassins versants de superficie supérieure à 200 km2 les temps de concentrations sont
supérieure à 12 heures, la pluie continue de durée égale au temps de concentration présente
une probabilité de retour quasi-nulle.
Le ruissellement est provoqué par une succession d’averses de plus en moins grande intensité.
Les débits de crues sont alors fonction des rythmes de succession des averses continues.
Le calcul de débit est obtenu en majorant ; par un coefficient variant de 1.5 à 4, le débit
moyen de la crue définie par la hauteur d’eau ruisselée sur le bassin versant, divisée par le
temps dit de concentration Tc.

27
Chapitre 3 : Etude hydraulique

Q : débit du bassin versant (m3/s) ;

P (t) : pluie tombée pour une durée égale au temps de concentration du bassin versant en
(mm) ;
Cr : coefficient de ruissellement ;
A : superficie du bassin versant en (km2).
Ka : coefficient d’abattement égale du rapport de pluie moyenne sur la pluie maximale et
varient de 0.5 à 0.9 ;
Tc : temps de concentration détermine par la formule de GIONDOTTI en heures.

3.2.2. Choix du coefficient de ruissellement

Le coefficient de ruissellement est défini comme le rapport de la hauteur d’eau ruisselée à la
hauteur d’eau précipitée.
Il est choisit compte tenu de la nature du bassin versant en particulier de sa couverture du type
de sol, de la pente, ainsi qu’en fonction du temps de retour de l’événement.
Le coefficient de ruissellement dépend des caractéristiques et conditions d’infiltration du sol,
l’intensité de précipitation, de la proximité de la nappe phréatique, du degré de compactage du
sol, de sa porosité, du couvert végétal, de la typologie d’occupation, de la pente du bassin, de
la récurrence, etc….
Le coefficient de ruissellement est estimé par la méthode préconisée de KENESSEY qui est la
somme de trois (03) coefficients partiels :
- C1 : dépendant de la pente du bassin versant.
- C2 : dépendant de la nature du sol (perméabilité).
- C3 : dépendant de la couverture végétale du bassin versant.

3.2.3 Choix du période de retour T

Pour chaque type d’infrastructure, les conditions d’écoulement et l’effet d’une crue
exceptionnelle doivent être appréciés.
Pour les buses en à la période de retour de T= 10 ans
Pour les dalots en à la période de retour de T= 50 ans
Pour les ponts en à la période de retour de T= 100 ans

28
Chapitre 3 : Etude hydraulique

3.3 Dimensionnement de l’ouvrage hydraulique

3.3.1 Calcul du débit de saturation
Le calcul des ouvrages ne peut se faire qu’après avoir disposé des contraintes de l’écoulement
naturel en aval recueillies jusqu’à au moins 100 m de l’ouvrage hydraulique.
Par ailleurs, les ouvertures des ouvrages hydrauliques de rétablissement sont généralement
inférieures à la section courante du ruisseau ou du thalweg pour des raisons de coûts.
Ce rétrécissement hydraulique n’est pas neutre pour son fonctionnement et notamment en
période d’évacuation du débit de crue.
L’écoulem.1ent est entièrement conditionné par les caractéristiques de l’oued.
Le débit est calculé par la formule de MANNING - STRICKLER.

Sm : section mouillée en (m2) ;

Rh : rayon Hydraulique en (m) ;

Pm : périmètre mouillé en (m) ;

I : pente moyenne de l’oued ;
K : coefficient de rugosité de Strickler.
Voici quelques ordres de grandeur du coefficient de Strickler

29
Chapitre 3 : Etude hydraulique

Nous avons adopté la valeur 30 pour les ponts

70 pour les dalots
80 pour les buses

De façon à contrôler les effets de l’érosion, la vitesse maximale de l’écoulement devra être
limitée à 4.5 m/s.
Pour la hauteur utile, si la hauteur totale

- Pour Ht  2,5 m  hauteur utile = 80 % hauteur totale.

- Pour Ht > 2,5 m  hauteur utile = hauteur totale - 0,50 m.

3.3.2 Conception des ouvrages hydrauliques

On distingue généralement 5 familles d’ouvrages :
Les buses circulaires, les dalots, les buses arches, les ouvrages à voûte cintrée, et les ouvrages
d’art.
Dans la mesure du possible, les produits industrialisés seront à rechercher plutôt que des
ouvrages coulés en place plus coûteux. Les ouvrages en béton armé, sous réserve de
dispositions constructives soignées, présentent d’excellentes garanties de solidité et de
longévité.
3.3.3 Protections des ouvrages hydrauliques
Il peut être envisagé de caler le radier de l’ouvrage hydraulique à au moins 0,30 m sous le
fond du lit du cours d’eau pour permettre la reconstitution d’un fond naturel dans l’ouvrage.
La surélévation du niveau amont des écoulements et l’accroissement des vitesses en sortie
d’ouvrage nécessitent le plus souvent des protections en amont et en aval des ouvrages.
Toute rectification du tracé nécessitera :
 la continuité de l’écoulement hydraulique ;
 la protection efficace des berges aux changements de direction par des techniques
pérennes relevant prioritairement du génie végétal.
Les techniques de renforcement par enrochements et gabions devront être réservées aux
sections fortement sollicitées par la vitesse de l’écoulement, si les enjeux sont importants en
terme de sécurité des personnes et des biens à fortes valeurs ajoutées ;
Les écoulements en pente importante p = 4% posent des problèmes spécifiques (détermination
de la hauteur d’eau amont, vitesse dans les ouvrages…

30
 Chapitre 3 : Etude hydraulique

Tableau 3.1 (modèle) : Débit d’apport et débit capable de l’ouvrage hydraulique

Vitesse
Surface du Débit de Pente
Débit capable Max de Type
B.V crue ouvrage Dimension
(m3/s) L’ouvrage d’ouvrage
(Km2) (m3/s) (%)
(m/s)

3800 1350 1 1350 4.5 ponts

3.4 CALCUL HYDRAULIQUE DES OUVRAGES TRAVERSANTS

Les ouvrages hydrauliques sont dimensionnés par la méthode citée ci-dessous avec la valeur
Q100, pour ce débit le dimensionnement devra respecter la condition suivante :
Ham
 1 .2
H
Avec : Ham : hauteur d’eau amont, H : hauteur de l’ouvrage

Il est souhaitable d’avoir une revanche de :

- R0.3m pour des OH où H  1.5m,
- R0.5m pour des OH où 1.5H2m
- 0.7mR1.5m au-delà.
3.4.1 méthode de vérification hydraulique

1. Le calcul de la section est effectué par la méthode de Manning Strickler :

Qs  K . S . Rh
2/3
. I 1/ 2
Avec :
I : Pente de pose de l’ouvrage (m/m).
S : Section mouillée (m²).
Rh : Rayon hydraulique moyen (m) avec :
Section mouillée
Rh 
Perimètre mouillé
Avec : K : Coefficient de Manning Strickler donné par :
 K = 60 pour les structures en béton coulées sur place
 K = 70 pour les structures préfabriquées en B.A. et pour les collecteurs circulaires.
2. Le calcul du remous amont est effectué grâce au théorème de Bernoulli :

31
 Chapitre 3 : Etude hydraulique

On peut obtenir par la formule Manning Strickler la hauteur normale (Hn) d’eau dans
l’ouvrage ainsi que la hauteur critique (Hc). En fonction de leur rapport, le régime
d’écoulement peut être fluvial ou torrentiel.
Si Hn  Hc, le régime d’eau est fluvial. On peut calculer Ham comme suit :
2
Ham  Hn 
VN
2g
1  Ke

Si Hn  Hc , le régime d’eau est torrentiel ; la hauteur d’eau à l’entrée de l’ouvrage est critique
.
Selon le niveau d’eau et le régime d’écoulement en aval de l’ouvrage il peut y avoir formation
d’un ressaut nécessitant une protection contre les affouillements (lit à enrocher par exemple).
La formule simplifiée fournissant la hauteur d’eau amont par excès s’écrit :

VC2
Ham  Hc 
2g
1  Ke
Dont : Ke : Coefficient de perte de charge à l’entrée de l’ouvrage hydraulique (fonction du
type de l’entée)
 Ke=0.5 pour l’extrémité avec mur de tête et murs en aile
 Ke=0.7 pour l’extrémité taillée en sifflet
 Ke=0.9 pour l’extrémité en saillie
VN : Vitesse d’écoulement correspondant à la hauteur d’eau normale dans l’ouvrage

VN  Q , SN est calculé par la formule de Manning Strickler

SN
VC : Vitesse d’écoulement correspondant à la hauteur d’eau critique dans l’ouvrage

VC  Q , SC est calculé par la formule de Manning Strickler

SC
3. V : Vitesse d’écoulement dans l’ouvrage hydraulique

V Q
SN
Dont : Q : Débit de charge (m3/s)
SN : Surface mouillée (m2)

32
Chapitre 3 : Etude hydraulique

3.5 Application pour notre ouvrage (Oued SAHEL)

Remarque : le résultat du débit de crue de notre projet est crue de projet pour une crue
centennale de : Q100 = 1350 m3/s (voir chapitre 2 : Etude hydrologique)

a. Détermination de la cote des Plus Hautes Eaux :(P.H.E)

Pour un ouvrage hydraulique de type pont et pour une finition correcte de la ligne rouge, il
faut déterminer la côte correspondant au débit de crue centennal. Pour ce faire un calcul
dérivant de la formule de Manning Strickler V = Kst x J1/2 x Rh2/3, et Q= Vx S permet
d’obtenir avec une approximation acceptable, le niveau d’une crue, connaissant le débit
hydrologique Q de la crue centennale et les caractéristiques géométriques du cours d’eau
(section mouillée S et périmètre mouillé P).

La formule la plus utilisée est celle de Manning Strickler.

V = Kst x Rh2/3 x I 1/2 et Q=VxS

Avec:
V: Vitesse moyenne d’écoulement (m/s)
R: Est le rayon hydraulique (en m)
J : La pente de l’écoulement (m/m)
Kst:= 1/n est le coefficient de Rugosité, et n varie en fonction de plusieurs facteurs :
- Nature des parois (terre, rocher, gravier fin, gravier grossier)
- Influence des variations, forme et de dimension de la section mouillée
- Influence des obstructions
- Influence de la végétation
- Importance des méandres

Dans le cas du franchissement d’Oued Sahel, nous prenons une valeur du coefficient de
rugosité à Kst = 30.
b. Calcul de la pente longitudinale J du lit l’oued :
On a pris au niveau du lit d’oued six (06) pentes J1, J2, J3, J4, J5 et J6
J1 = 168,63-168.28 / 38,40 = 0.009 m/m
J2 = 168,28-168,14 / 32,52 = 0.004 m/m
J3 = 168,14-167,72 /28,63 = 0,014 m/m
J4 = 166,49-166,02 / 41,59 = 0.011 m/m
J5 = 166,49-165.71 / 93,80 = 0.008 m/m

33
Chapitre 3 : Etude hydraulique

J6 = 166,28-165,63 / 71,89 = 0.009 m/m

J moyen = (J1+J2+J3+J4 +J5+J6)/ 6 = 0.009 m/m
Calculs des vitesses et débits au site de franchissement sur Oued Sahel (pour K = 30 et J =
0,009 m/m)

c. Résultats de calcul

Côte (m) Hauteur Pm (m) Sm Rh Debit Vitesse

3
d’eau(m) (m2) (m) (m /s) d’écoulement
(m/s)
168,24 0,50 83,672 27,752 0,332 37,826 1,363
168,74 1,00 98,640 71,485 0,725 164,130 2,296
169,24 1,50 112,171 119.621 1,066 364,844 3.050
169,74 2,00 123.480 172.119 1,394 611,022 3,550
170,24 2,50 133.480 224.716 1,684 904,931 4,027
170.74 3.00 143.480 277.283 1,933 1224,276 4,415
170,94 3,20 147.480 297.693 2,019 1353,132 4,545

d. Choix du TIRANT D’AIR

Un cours d’eau charrie souvent détritus, corps flottants, branchages, etc… Au passage sous un
pont ces derniers s’accrochent et bouchent peu à peu la section d’écoulement mettant ainsi
l’ouvrage en risque d’être en charge .C’est pourquoi, il est obligatoire de prévoir un tirant
d’air pour diminuer ce risque d’obstruction partielle ou totale du pont .Pour des ponts de
grandes longueurs le tirant d’air varie de 1,00 m à 2,50 m .. Le tirant d’air considéré est de H
= 1.20 m

Le gabarit de l’ouvrage (ouverture verticale) est de : 3.20 + 1.20 = 4.40m.

e. Détermination de la cote finale de la ligne rouge de l’ouvrage (pont):

Si nous prenons la hauteur de construction de l’ouvrage h = 1,60 m, la côte finie (de

roulement) est de : 170,94 + 1.20 + 1,60 = 173,74 m

34
Chapitre 3 : Etude hydraulique

3.6 Conclusion :
Le dimensionnement d’un ouvrage hydraulique en général et de drainage en particulier est
une grande responsabilité pour l’hydraulicien dont l’erreur n’est pas permise afin d’éviter des
défaillances après réalisation. Pour cela le choix de la méthode et les formules adéquates
doivent se faire convenablement, sans oublier le choix des données à utiliser. Pour notre cas,
le calcul des plus hautes eaux (PHE) nous permet de poser l’ouvrage de traversé sans aucun
risque de débordement de l’eau ou de déstabilisation de l’ouvrage.
Selon le résultat trouvé, l’ouvrage de traversé (pont) doit être au dessus d’une hauteur de
4.40m à partir du lit de l’oued traverse.

35
Chapitre 4 : Présentation des logiciels HEC-RAS

Introduction :
HEC-RAS est un logiciel intégré pour l'analyse hydraulique qui permet de simuler les
écoulements à surface libre. Il a été conçu par le Hydrologic Engineering Center du U.S Army
Corps of Engineers. Il s'agit d'une nouvelle version d'un modèle hydraulique auparavant
nommé HEC-2, qui comporte maintenant une interface graphique permettant d'éditer, de
modifier et de visualiser les données d'entrée, de même qu'observer les résultats obtenus. Il est
présentement utilisé dans plusieurs firmes d'ingénierie et organismes gouvernementaux. La
version 5.0.6 de HEC-RAS permet de simuler les écoulements permanents et non permanents
4-1-Principe d’utilisation de logiciel HEC-RAS :
HEC-RAS est un logiciel intégré pour l’analyse hydraulique qui permet de simuler
les écoulements à surface libre. Il a été conçu par le Hydrologic Engineering Center de l’US
Army Corps of Enginners.
La version 4.0 Beta de HEC-RAS que nous utiliserons dans notre travail, permet
de simuler les écoulements à surface libre permanents et non permanents, le transport
de sédiments, la qualité des eaux et certaines fonctions facilitant la Conception d'ouvrages
hydrauliques.

Pour démarrer HEC-RAS, double-cliquez sur l’icône placé sur le bureau, ou bien
allez dans le menu Démarrer et choisissez le programme HEC-RAS4.0 Beta. Fait important
à noter, HEC-RAS utilise comme symbole décimal le point, et non la virgule utilisée
habituellement dans notre système d’unité. Un message d’erreur peut apparaître lors
du démarrage si le symbole décimal spécifié pour votre ordinateur n’est pas le point.
Vous pouvez changer ce symbole dans le Panneau de configuration de Windows.
Suite au démarrage de HEC-RAS vous devez obtenir la fenêtre principale de HEC-RAS qui
illu

Figure 4-1 : La fenêtre principale de HEC-RAS.

37
Chapitre 4 : Présentation des logiciels HEC-RAS

La barre de menu comprend toutes les fonctions disponibles de HEC-RAS. Le menu

File permet d’ouvrir, de créer et de sauvegarder un projet (avec extension.prj). Les autres
options permettent de modifier le titre du projet ou d’effacer tous les fichiers reliés à un
projet. Cette dernière option est utilisée avec précaution.
Outre le menu File, le menu Edit est également important et permet de spécifier
les données de base décrivant le système hydrographique devant être modélisé ainsi
que les conditions du débit dans ce système.
La principale étape de la création d’un projet de modélisation avec HEC-RAS est
de définir la géométrie du cours d’eau, au moyen de sections transversales. Cette étape est
réalisée en choisissant l’option Geometric Data…Toutes les données reliées à la géométrie
du cours d’eau sont sauvegardées dans les fichiers dont l’extension est .G** ou ** désigne
les chiffres représentant une numérotation séquentielle.
Une autre option disponible dans le menu Edit est Steady Flow Data…Cette option
permet de spécifier la ou les valeurs de débit dans le cours d’eau qui devront être modélisés,
ainsi que les conditions limites de l’écoulement. Les informations qui sont spécifiées avec
cette option sont contenues dans les fichiers .F**.
4-2-Etapes de la modélisation :
4-2-1-Projet :
Dans le menu File, choisir l’option New Projet pour obtenir une fenêtre similaire à
celle illustrée plus bas (Figure V-2).
La première étape est de choisir le dossier de travail où le projet sera sauvegardé.
Un nouveau répertoire peut d’ailleurs être crée avec le bouton Create Folder. Le titre du
projet sera par la suite spécifié dans la ligne haute et un nom de fichier doit aussi être
indiqué dans la case suivante, tout en conservant l’extension .prj. Il ne faut pas utiliser
les accents pour un nom de fichier. Après avoir appuyé sur OK, un message apparaît pour
confirmer les informations soumises. Appuyez à nouveau sur OK et les données
seront sauvegardées dans le fichier et le répertoire indiqué.
Avant d’entrer les informations reliées à la géométrie et aux débits, il faut spécifierle
système d’unité utilisé. Dans le menu Option, allez à Units System et choisissez System
International (Metric System).

38
Chapitre 4 : Présentation des logiciels HEC-RAS

Figure 4-2: Fenêtre New Project.

4-2-2-Géométrie des sections :

Après avoir défini les principaux paramètres du projet, la deuxième étape est de définir
les caractéristiques géométriques du système modélisé. Sélectionnée Edit / Geometric Data….
et la fenêtre Base Geometric Data apparaît (Figure V-3).
Pour faciliter le tracé du tronçon de rivière d’étude, on peut ajouter une image de fond
à cette fenêtre. Appuyez sur le bouton Add /Edit background pictures…et appuyez sur Add
afin de sélectionner le fichier image.
Après avoir appuyé sur Close, un message apparaît car la taille de l’image est plus
grande que la taille de l’écran. Il est possible d’ajuster l’échelle en choisissant dans le menu
View l’option Full Plot. Pour sauvegarder les étapes accomplies, allez à File / Save Geometric
Data as… Sauvegardez les données de la géométrie sous un fichier. Pour représenter le
tronçon à étudier, cliquez sur le bouton River Reach. Le curseur se transforme en crayon
et vous pouvez alors dessiner le tronçon de rivière à l’aide d’une suite de points qui vont de
l’amont vers l’aval, en suivant le centre de la rivière sur l’image de référence. Vous double-
cliquez pour indiquer la fin du tronçon. A ce moment, une fenêtre apparaît vous demandant
d’indiquer le nom de la rivière et le nom du tronçon.

39
Chapitre 4 : Présentation des logiciels HEC-RAS

Figure 4-3 : Fenêtre Geometric Data.

Les autres caractéristiques géométriques nécessaires à cette étape sont les

sections transversales à différents endroits de la rivière. Les sections transversales sont
représentées au moyen de points représentant des coordonnées X-Y, où X est la distance
par rapport à un point de référence arbitraire placé sur la rive et Y est l’élévation du
fond de la rivière.
Pour entrer les données, vous devez cliquer sur le bouton Cross-Section dans la partie
gauche pour obtenir une fenêtre similaire (Figure V-4). Dans le menu Option, choisissez Add
a new Cross Section…Il vous est alors demandé de définir la station (River Station) de cette
section transversale qui est en fait un identificateur numérique. HEC-RAS place sur un
tronçon de rivière les stations en ordre décroissant de la partie amont ver la partie aval. La
signification de chacun des termes que l’on y retrouve est décrite ci-dessous :

40
Chapitre 4 : Présentation des logiciels HEC-RAS

Figure 4-4 : Fenêtre Cross Section Data.

 River : Nom de la rivière sur laquelle la nouvelle section sera ajoutée ;

 Reach : Nom du tronçon de la rivière sur laquelle la nouvelle section sera ajoutée ;
 River station : Identification numérique de la section transversale ;
 Description : Commentaire de l’utilisation sur cette section ;
 Cross-section coordinates : Coordonnées relatives dans le plan X-Y des points
définissant la section transversale ;
 Downstream Reach Lengths : Distance en mètres jusqu’à la prochaine section
transversale située en aval. LOB signifie la partie gauche de la plaine inondable et
ROB sa partie droite, alors que Channel désigne le lit principal de la rivière ;
 Manning’s n values : Coefficient de Manning de chaque portion de la section
transversale ;

41
Chapitre 4 : Présentation des logiciels HEC-RAS

Tableau 4-1 : Valeurs estimées pour le coefficient de Manning en fonction de la nature

des parois.

n (Mannig) K (Strickler)
Nature des parois
(s/m1/3) (m1/3/s)
Béton lisse 0.011 à 0.013 77 à 91
Béton brut 0.013 à 0.016 62 à 77
Canal en terre, non enherbé 0.017 60
Canal en terre, enherbé 0.02 50
Rivière de plaine, sans végétation 0.025 à 0.029 35 à 40
arbustive
Rivière de plaine, large, végétation 0.033 30
peu dense
Rivière à berges étroites très 0.067 à 0.1 10 à 15
végétalisées
Lit majeur en prairie 0.033 à 0.05 20 à 30
Lit majeur en vigne 0.067 à 0.1 10 à 15
Lit majeur urbanisé 0.067 à 0.1 10 à 15
Lit majeur en forêt >0.1 <10

 Main channel bank stations : Coordonnées, dans le plan X seulement, des limites
gauches et droite du lit principal de la rivière. Les valeurs fournies doivent
correspondre à une valeur déjà présente dans la partie Cross-section coordinates.
Après avoir entré les données, cliquez sur Apply Data. Vous pouvez alors visualiser
la section transversale que vous venez de définir dans la partie adjacente de la fenêtre. Les
autres sections seront rentrées de la même manière.
Pour augmenter la stabilité des calculs de niveau par le modèle. Il est recommandé
d’avoir une distance raisonnable entre deux sections transversales. En assumant que la pente
et les propriétés des sections varient de façon linéaire entre deux sections consécutives,
la fonction XS interpolation dans le menu Tools permet d’ajouter par interpolation de
nouvelles sections entre deux sections existantes (Figure V-5).

42
Chapitre 4 : Présentation des logiciels HEC-RAS

Figure 4-5 : Fenêtre XS interpolation.

Pour terminer, vous pouvez vérifier les données spécifiées avec le menu Tables. Les
mêmes paramètres pour toutes les sections y sont regroupés à l’intérieur d’un même tableau et
peuvent y être modifiés. Avant de quitter la fenêtre Geometric Data, sauvegardez les
informations que vous avez rentrées.
4-2-3-Débit et conditions aux limites :
L’étape suivante de la modélisation hydraulique sur HEC-RAS est de spécifier les
débits utilisés pour calculer les profils d’écoulement. Dans la fenêtre principale, cliquez sur le
bouton Steady Flow Data (Figure V-6). Pour entrer les valeurs des débits, indiquées d’abord
Number of Profiles. Par la suite, dans le menu Options, allez à Edit Profil Names…et
indiquée les noms de profils appropriés. Les valeurs de débit qui sont entrées aux cases
correspondantes sont représentatives de la section amont et sont considérés valides sur tout le
tronçon de rivière. HEC-RAS permet toutefois de représenter des changements de débits aux
sections transversales, lorsqu’un affluent important entraîne un changement de débit dans le
tronçon.

43
Chapitre 4 : Présentation des logiciels HEC-RAS

Figure 4-6: Fenêtre Steady Flow Data.

Lorsque les débits correspondants aux différents profils devant être calculés sont
spécifiés, il faut par la suite définir les conditions limites de l’écoulement en cliquant sur
le bouton Reach Bondary Conditions (Figure 4-7). Les conditions limites sont nécessaires
pour calculer la hauteur d’eau initiale aux extrémités de chaque tronçon.
Pour un écoulement fluvial, seulement les conditions à l’aval sont nécessaires, tandis
que pour un écoulement torrentiel, les conditions à l’amont seulement sont nécessaires. Selon
le régime d’écoulement modélisé, une seule des deux conditions limites peut être indiquée.
On a le choix entre trois conditions aux limites : hauteur critique, hauteur normale et hauteur
connue. Avant de passer à l’étape suivante, sauvegardez ces données avec File / Save Flow
Data as…Entrez le titre et quittez cette fenêtre pour revenir au menu principal de HEC-RAS.

Figure 4-7: Fenêtre Reach Bondary Conditions.

44
Chapitre 4 : Présentation des logiciels HEC-RAS

4-2-4-Simulation hydraulique :
La dernière étape nécessaire à la modélisation avec HEC-RAS est de définir le plan
utilisé. Cliquez sur le bouton Perform a Steady Flow Simulation pour obtenir la fenêtre
présentée plus bas. Le Geometry File et Steady Flow File que vous avez créés y sont indiqués.
Dans le menu File, choisissez New Plan. Lorsque requis, entrez le titre et l’identificateur
indiqué sur la figure III-10.
Choisissez le régime d’écoulement et appuyez sur Compute pour débuter la simulation.
Une fenêtre montrant la progression de la simulation s’ouvrira et les calculs s’effectueront
(FigureV-8). Lorsque le programme a terminé, vous pouvez fermer la fenêtre Hydraulic
Computation, ainsi que la fenêtre Steady Flow Analysis.

Figure 4-8: Fenêtre Steady Flow Analysis.

4-2-5-Visualisation des résultats :

A la suite de la simulation les résultats sont automatiquement sauvegardés dans un
fichier. Un des résultats intéressant à consulter est la vue en profil du tronçon simulé. Cette
fonction, View Profiles, est située dans la barre de boutons. Allez également dans le menu
Option de cette fenêtre. Vous y trouverez les options d’affichage, telles que Zoom et Pan. Il y
a aussi d’autres options qui permettent d’afficher les résultats d’un ou plusieurs plans,
d’ajouter ou supprimer l’affichage de certains profiles, de changer le tronçon dont les résultats
sont actuellement affichés. Ce menu Options est d’ailleurs similaire pour différentes fenêtre
graphiques permettant de visualiser les résultats (Figure 4-9).

45
Chapitre 4 : Présentation des logiciels HEC-RAS

Figure 4-9. Fenêtre de profil de surface de l’eau.

Une autre option graphique intéressante est accessible via le bouton View 3D multiple
cross section plot…. Les sections transversales de début et de fin peuvent être modifiées pour
afficher qu’une partie du tronçon à l’étude. Rotation et Azimuth permettent quant à eux
de modifier l’angle de vue. Pour ne voir qu’une seule section. Deux autres boutons, View
Detailed Output at XS et View Output at Multiple…, permettent de voir les résultants sous
forme tabulaire (Figure 4-10).

Figure 4-10 : Fenêtre de surface de l’eau 3D.

46
Chapitre 4 : Présentation des logiciels HEC-RAS

Finalement, comme on peut voir dans le tableau des résultats des sections transversales,
il y a plusieurs variables qui sont calculées par HEC-RAS en plus des niveaux d’eau. Dans
certains cas, on s’intéressera aux vitesses afin de quantifier par exemple les possibilités
d’érosion.
Dans la fenêtre principale de l’HEC-RAS, sélectionnez erreur avertir les notes
sommaires de vue tirer vers le bas menu est appuyez sur fermer lorsque terminer (Figure 4-
11).

Figure 4-11 : Résume des erreurs et des avertissements.

4-3-Théorie de base de calcul du modèle HEC-RAS :
L’objectif principal de ce programme est tout à fait simple, il est conçu pour exécuter
le calcul hydraulique unidimensionnel pour tous les réseaux de canaux naturels et artificiels
par le calcul des hauteurs de la surface de l’eau en toute section d’intérêt pour un ensemble
de données d’écoulement en régime permanent et par des hydrogrammes de propagation
des crues pour la simulation de l’écoulement en régime non permanent.
HEC-RAS est un logiciel basé sur le modèle 1D de Saint Venant permanent ou non
permanent de calcul de la ligne d’eau en graduellement varié. Il résout « l’équation de
l’énergie unidimensionnelle », les pertes étant évaluées par la formule de frottement au fond
de Manning-Strickler et par des formules de contraction/expansion de l’écoulement. Pour les
situations rapidement variées telles que les ressauts hydrauliques, les écoulements à proximité

47
Chapitre 4 : Présentation des logiciels HEC-RAS

des ponts, et les confluences de rivière, l’équation de l’énergie est remplacée par l’équation
de quantité de mouvement.
4-3-1-Les équations de Saint-Venant 1D :
Le logiciel HEC-RAS utilise les équations de Saint-Venant 1D pour relier les hauteurs
d'eau et les débits. Ces équations se déduisent des équations de Navier-Stokes grâce à des
simplifications liées au modèle de la rivière.
4-3-2-Le cas stationnaire (SteadyFlow) :
Comme nous l'avons dit dans la partie précédente, le maillage consiste en une
succession de sections. Le logiciel calcul la hauteur d'une section à partir de la hauteur d'eau
de la section précédente à l'aide de l'équation de conservation de l'énergie (conservation de la
charge).

Avec :

Z1, Z2 : élévation du fond du canal ;

Y1, Y2 : hauteur normale de la section ;
V1, V2 : vitesse moyenne d’écoulement ;
α1, α2: coefficients de pondération de vitesse ;
g : accélération de gravité ;
he : pertes de charge.

Le coefficient de pertes de charge he se calcule à l'aide de la formule suivante:

Avec L la distance entre les deux profils, Sf la pente de la ligne d'énergie (pente de
frottement) et C : le coefficient d'expansion ou de contraction.

Pour calculer le débit on utilise l'équation de Manning :

Q=K.Sf1/2 avec K = 1/n*A*R2/3

48
Chapitre 4 : Présentation des logiciels HEC-RAS

Avec :
K : flux pour la subdivision.
n : coefficient de frottement de Manning.
A : section mouillée.
R : rayon hydraulique.
Sf : pertes de charges entre les deux sections.
On voit donc que le débit sur chaque section est calculé à partir de la hauteur d'eau
puisque le rayon hydraulique, la section mouillée ainsi que les pertes de charge dépendent de
la hauteur d'eau. On peut visualiser cette conservation de la charge à l'aide de la figure 4-12.

Figure 4-12 : Conservation de la charge entre deux sections (source: HEC-RAS).

Pour obtenir des résultats satisfaisants il est nécessaire d'évaluer correctement les pertes
de charges, et donc le frottement. Pour cela HEC-RAS découpe chaque section en plusieurs
domaines verticaux et évalue le frottement sur chacun de ces domaines. En effet le frottement
dépend de la vitesse et on se doute bien que la vitesse n'est pas la même sur les bords, dans le
lit majeur ou dans le lit mineur. Le découpage se fait comme sur la figure 4-13 qui suit:

49
Chapitre 4 : Présentation des logiciels HEC-RAS

Figure 4-13 : Découpage de la section pour le calcul du frottement (source: HEC-RAS).

Pour faire les calculs lorsqu'il n'y a pas de stationnarité une première approche pourrait
être de calculer les débits et les hauteurs pour chaque temps, et donc chaque débit d'entrée, à
l'aide de l'étude stationnaire. Cependant une étude plus rigoureuse consiste à faire le calcul à
l'aide d'équations qui prennent en compte le non stationnarité.
4-3-3-Le cas instationnaire (Unsteady Flow) :
Les lois physiques qui régissent l’écoulement d’eau dans un canal sont :
Le principe de la conservation de la masse et le principe de la conservation de la quantité de
mouvement. Ces lois sont exprimées mathématiquement sous forme d’équations
différentielles partielles, qui sont désigné ci-après sous le nom des équations de continuité et
de la quantité de mouvement.

Figure 4-14 : Volume de contrôle élémentaire de dérivation de la continuité et

d’équations de moment (Inflow : débit entré et Outflow : débit sortie).
On obtient l'équation de conservation de la masse:
∂A /∂t + ∂Q/∂x −ql = 0
Avec A la section, Q le débit et ql le débit entrant latéral par unité de longueur.
L'équation de conservation de la quantité de mouvement s'écrit:
∂Q/∂t + ∂QV/∂x +gA(∂z/∂x+Sf) = 0
Avec ∂z/∂x la pente de la hauteur d'eau, Sf la pente de la ligne d'énergie (pente de frottement),
A la section et V la vitesse moyenne.
Pour résoudre ces équations, le logiciel utilise la théorie des différences finies qui est
à la fois pratique et simple pour ce cas 1D. Il utilise un schéma implicite qui permet d'utiliser
un pas de temps beaucoup plus important que dans le cas du schéma explicite.

50
Chapitre 4 : Présentation des logiciels HEC-RAS

4-4 Conclusion :
Dans ce chapitre nous avons présenté le mode de travail de logiciels HEC-RAS avec des
exemples de simulation d’écoulement permanent et non permanent. La méthode utilisée par le
logiciel HEC-RAS apporte des réponses face au risque d’inondations, Des renseignements
complémentaires sont également fournis qui permettent d’estimer le débit maximal et la
vitesse maximale de l’onde ainsi la hauteur d’eau maximale pour procéder à des résultats
cartographiques facilement lisibles

51
Chapitre :05 Calcul de la PHE (Plus Haute Eau) par HEC-RAS

5.1 Introduction : L’objectif de ce chapitre est de déterminer la variation du niveau de

remonté des eaux dans la section ou l’ouvrage sera réalisé sur oued Sahel appelée la plus
haute eau (PHE). Ce travail se fera en utilisant le logiciel HEC-RAS avec ARC GIS. Ces
derniers nous permettent la simulation du niveau d’eau de l‘oued du bassin versant de Sahel
après le passage de la crue.
5.2. Le logiciel Hec Ras :

HEC-RAS est un logiciel intégré pour l’analyse hydraulique qui permet de simuler les
écoulements à surface libre. Il a été conçu par le Hydrologic Engineering Center de l’US
Army Corps of Enginners.

La version 5.0 Beta de HEC-RAS que nous utiliserons dans ce travail, permet de simuler les
écoulements à surface libre permanents et non permanents, le transport de sédiments la qualité
des eaux et certaines fonctions facilitant la conception d'ouvrages hydrauliques

5.3. Les étapes de développement du model HEC-RAS:

a. Ouvrir HEC-RAS

b Entrer les informations du projet, Sélectionnez File, New Project

Créer un répertoire de travail : Create Folder

Ajouter Nom du fichier : File Name
Ajouter titre du projet: Title
Cliquer ok
Ajouter la description du projet : Description
Enregistrer : File, Save Project

53
Chapitre :05 Calcul de la PHE (Plus Haute Eau) par HEC-RAS

c. Choisir système d’unités :

Options, Unit system, System International, Ok

Set as defaults for new projects

d. Introduire les données géométriques:

Sélectionner Edit, Geometric Data

-Ajouter:
Reaches : tronçon
Cross Sections : Sections transversales
-Dessinez le tronçon en sélectionnant le bouton Reach (cliquez une fois)

54
Chapitre :05 Calcul de la PHE (Plus Haute Eau) par HEC-RAS

-Déplacez la souris sur l'écran pour dessiner le tronçon

-Définir le tronçon en procédant d'amont en aval

-Avec le curseur dans fenêtre de dessin, cliquez une fois sur le bouton de la souris pour
commencer

-Déplacez le curseur, cliquez sur le bouton de la souris une fois pour ajouter un point ou
changer de direction
-Au point Aval, double clic sur la souris pour finir.

- Enter le nom de la rivière River Name et le nom du tronçon Reach Name

Reaches

Cross
Sections

55
Chapitre :05 Calcul de la PHE (Plus Haute Eau) par HEC-RAS

e. Ajouter les sections transversales du tronçon choisi :

-Sélectionner le bouton Cross Section (un seul clic)

-Sélectionner: Options, Add a new Cross Section

-Enter le numéro d’emplacement de la station, saisir (1.0) dans River Sta

-La station doit être une valeur numérique et sont définis comme suit : Station inférieure aval,
stations plus élevées en amont.

Station Elevation
5 171
70 170.10
95 169.94
110 168.28
125 167.98
155 167.84
180 167.74
195 168.34
210 169.63
260 169.57
325 170.91

56
Chapitre :05 Calcul de la PHE (Plus Haute Eau) par HEC-RAS

Downstream Reach Lengths : Distance en mètres jusqu’à la prochaine section transversale

située en aval.
LOB signifie la partie gauche de la plaine inondable et ROB sa partie droite, alors que
Channel désigne le lit principal de la rivière ;

Ajouter n la valeur de Manning :

LOB Channel ROB
0.03 0.03 0.03

f. Répéter pour les sections transversales supplémentaires

Edit, Copy, Options Add a New Cross Section

Name 01
Downstream Reach Lengths:
Exit Cross Section Data Editor

g. Répéter pour les sections transversales supplémentaires

Tools, XS Interpolation, Within a Reach
Maximum distance:
Interpolate XSs, Close
Exit Cross Section Data Editor

57
Chapitre :05 Calcul de la PHE (Plus Haute Eau) par HEC-RAS

Les sections transversales interpolée ont "*" par le nom de la station. On peut changer ce nom
comme suit

Options, Rename River Station

Enregistrer les données géométriques :

Save Geometric Data, donner le non du fichier

h. les données de la simulation :

Démarrez Edit - Steady Flow Data

Entrez le débit de la rivière à la section transversale amont : 1353.132 m³/s

58
Chapitre :05 Calcul de la PHE (Plus Haute Eau) par HEC-RAS

Enter les conditions aux limites : Reach Boundary Conditions- Critical Depth

SAVE STEADY FLOW DATA, file name automatically assigned

i. Simulation

Run – Steady Flow Analysis

59
Chapitre :05 Calcul de la PHE (Plus Haute Eau) par HEC-RAS

Select Geometry File

Select Steady Flow File
Flow Regime: Subcritical
Save Plan: donner nom fichier

COMPUTE
Steady Flow Analysis
10. Affichages des résultats:

View, X-Y-Z Perspective Plot

View, Water Surface Profile

print using File, Print from screen or copy to Windows clipboard

60
Chapitre :05 Calcul de la PHE (Plus Haute Eau) par HEC-RAS

61
Chapitre :05 Calcul de la PHE (Plus Haute Eau) par HEC-RAS

62
Conclusion Générale

Conclusion générale

Apres une étude hydrologique et hydraulique comme première méthode, et étude de

simulation par le logiciel HEC RAS, nous avons calculé la plus haute eau sur un tronçon de
l’oued sahel qui sera traversé par un ouvrage de franchissement.
Avec les différentes étapes de notre étude, nous avons identifié le niveau de la remonté des
eaux en période de crues pour différentes périodes de retour de 100ans .
La période de retour est celle qui convienne au dimensionnement des grands ouvrages
hydrauliques.
D’après les résultats de simulation par HEC RAS et de calculs par des formules, le résultats
de la PHE (plus haute eau) trouvé par deux méthodes sont proches et le résultat est
satisfaisant.
Enfin, nous espérons avoir contribué avec ce modeste travail à donner un aperçu général
sur les méthodes de calcul de remonté des eaux dans un tronçon et section d’un oued en
périodes

64
Bibliographie
- zouaoui abdelhamid. «etude faisabilite et d'execution d'une retenue collinaire sur oude
elkhaloua (m'sila) ».
- MokhtariChoayb, « modélisation Pluie débit par logiciel HEC-HMS de sous bassin
versant l'oued SAHEL » Mémoire Master, Université Mohamed Boudiaf M’sila Juin
2016
- ould fah Mohamed Mahmoud : Etude hydrologique du bassin versant de la
Soummam
- Mahfoud Baaidja: Evaluation des pertes en sols et leur identification dans le
bassin versant d’Oued Azzerou (BBA )
- Guide technique Assainissement routier - Sétra - Service d'Études techniques des routes
et autoroutes, octobre 2006.
- Manuel Pratique D’hydrologie Par Bénina Touaibia Mars 2004
- L’hydrologie de surface par sari AHMED 2002
- Ouvrages d’assainissement - Volume1：Mise en œuvre- Ouvrages annexes Normes
applicables au fascicule 70（AFNOR）
- Ouvrages d’assainissement - Volume2：Canalisations d’assainissement Normes
applicables au fascicule 70（AFNOR）
- Guide de conception des réseaux d’assainissement de la plate-forme autoroutière et
des rétablissements des écoulements naturels.
- Cours .hydrologie ben kadja.
- Merrine –soikaichbach ;livre ,hydrologie et hydrolique ,2003-2004.

66
Study summary
After a hydrological and hydraulic study as the first method, and a simulation
study by HEC RAS software, we calculated the highest water on a section of the
Sahel wadi which will be crossed by a crossing structure.
With the different stages of our study, we identified the level of upwelling
during flooding for different return periods of 100 years.
The return period is that which is suitable for the sizing of large hydraulic
structures.
According to the results of simulation by HEC RAS and calculations by
formulas, the results of PHE (highest water) found by two methods are close and
the result is satisfactory.
Résumé
Apres une étude hydrologique et hydraulique comme première méthode, et
étude de simulation par le logiciel HEC RAS, nous avons calculé la plus haute
eau sur un tronçon de l’oued sahel qui sera traversé par un ouvrage de
franchissement.
Avec les différentes étapes de notre étude, nous avons identifié le niveau de la
remonté des eaux en période de crues pour différentes périodes de retour de
100ans .
La période de retour est celle qui convienne au dimensionnement des grands
ouvrages hydrauliques.
D’après les résultats de simulation par HEC RAS et de calculs par des
formules, le résultats de la PHE (plus haute eau) trouvé par deux méthodes sont
proches et le résultat est satisfaisant.
‫الملخص‬
‫ قمنا‬، HEC RAS ‫ ودراسة محاكاة بواسطة برنامج‬، ‫بعد دراسة هيدرولوجية وهيدروليكية كأول طريقة‬
.‫بحساب أعلى مستوى للمياه في جزء من وادي الساحل الذي سيتم عبوره بواسطة هيكل معبر‬
‫ حددنا مستوى ارتفاع المياه أثناء الفيضانات لفترات عودة مختلفة‬، ‫من خالل المراحل المختلفة لدراستنا‬
.‫ عام‬011 ‫تبلغ‬
.‫فترة اإلرجاع هي تلك المناسبة لتحجيم الهياكل الهيدروليكية الكبيرة‬
‫ (أعلى‬PHE ‫ فإن نتائج‬، ‫ والحسابات بواسطة الصيغ‬HEC RAS ‫وفًقا لنتائج المحاكاة بواسطة‬
.‫مستوى للمياه) التي تم العثور عليها بطريقتين متقاربة والنتيجة مرضية‬