RAPPORT DE STAGE :

ANALYSE HYDROLOGIQUE ET
MODELISATION HYDRAULIQUE DU
SOUS-BASSIN VERSANT D’OURIKA

2e année de génie civil, option hydraulique

Institut spécialisé des travaux publics Marrakech ISTPM

Réalisé par : -SLIMAN CHOUKOUK Encadrée par : Mme RAJI LLAH ASMA
-HAMZA ENOUIRI Fonction : INGENIEUR Au sein
du Service de Gestion du
Domaine Public Hydraulique

Soutenu : le 21 Mai 2025

LISTE DES FIGURES : ................................................................................................................................ 4
LISTE DES TABLEAUX : ............................................................................................................................. 5
REMERCIEMENT : .................................................................................................................................... 6
INTRODUCTION : ..................................................................................................................................... 7
I. CONTEXTE INSTITUTIONNEL : L’ABHT ET LA GESTION DE L’EAU : .................................................. 9
II. OBJECTIFS DE L’ABHT :..................................................................................................................... 9
IV. SITUATION GEOGRAPHIQUE DU BASSIN VERSANT TENSIFT. ................................................... 11
I. SITUATION GEOGRAPHIQUE DU SOUS BASSIN VERSANT D’OURIKA A LA STATION AGHBALOU 13
II. LA DELIMITATION DU BASSIN VERSANT D’OURIKA AGHBALOU ................................................. 14
1. Présentation du logiciel Arc GIS : Fonctionnalités et Composants.......................................... 14
2. Les cordonnées Lambert et de l’exutoire du bassin versant ................................................... 14
3. Le réseau hydrographique :...................................................................................................... 15
III. CADRE GEOLOGIQUE : .............................................................................................................. 16
V. CARTOGRAPHIE DE LA COUVERTURE VEGETALE ET DE FACTEUR D'EROSION HYDRIQUE DU
BASSIN VERSANT :................................................................................................................................. 17
1. La végétation : ........................................................................................................................... 17
2. L’érosion du bassin versant ....................................................................................................... 18
IV. LE REGIME PLUVIOMETRIQUE ET HYDROMETRIQUE DU BASSIN VERSANT D’OURIKA A LA
STATION AGHBALOU : ........................................................................................................................... 19
1. Le régime pluviométrique ........................................................................................................ 19
a. Précipitations annuelles (mm) : ........................................................................................... 19
b. Précipitations moyennes saisonnières ................................................................................. 20
c. Précipitations moyennes mensuelles :................................................................................. 21
d. Précipitations journaliers maximales : ................................................................................. 22
2. Le régime hydrométrique ......................................................................................................... 23
a. Débits moyens annuels......................................................................................................... 23
b. Débits moyens saisonniers : ................................................................................................. 24
c. Débits moyens mensuels : .................................................................................................... 24
I. CARACTERISTIQUES GEOMETRIQUES :......................................................................................... 27
1. Surface (A) : .............................................................................................................................. 27
2. Périmètre (P) :........................................................................................................................... 27
3. Forme : ...................................................................................................................................... 27
1. La longueur de la cour d’eau principale (Li) :........................................................................... 28

1
 2. Le rectangle équivalent :........................................................................................................... 28
II. CARACTERISTIQUES TOPOGRAPHIQUES ...................................................................................... 29
1. La pente : ................................................................................................................................... 29
a. La pente moyenne : .............................................................................................................. 29
b. La carte de pente : ................................................................................................................ 29
2. Hypsométrie :............................................................................................................................ 30
a. La carte hypsométrique ........................................................................................................ 30
b. La courbe Hypsométrique. ................................................................................................... 31
c. Indice de pente IG : ............................................................................................................... 33
III. CALCUL LE TEMPS DE CONCENTRATION................................................................................... 35
1. Formule espagnole : ................................................................................................................. 35
2. Formule de Californienne : ....................................................................................................... 35
3. Formule de Passini : .................................................................................................................. 35
4. Formule de Kripich: ................................................................................................................... 36
5. Formule de Turazza : ................................................................................................................. 36
6. Formule de Ventura .................................................................................................................. 36
V. CALCUL DES DEBITS DE POINTE : .................................................................................................. 38
1. Ajustement statistique :............................................................................................................ 38
a. Présentation du logiciel Hyfran-plus : .................................................................................. 38
b. Lois d’ajustement : ................................................................................................................ 38
c. L’analyse des graphes :.......................................................................................................... 38
d. L’analyse des critères statistiques : ...................................................................................... 40
2. Transformation pluie – débit : .................................................................................................. 42
a. Présentation du modèle utilisé : .......................................................................................... 42
b. Paramétrage du modèle ....................................................................................................... 43
c. Paramètres clés de sortie : ................................................................................................... 43
1. Méthodes empiriques :............................................................................................................. 45
a. Débits par Fuller 2 : ............................................................................................................... 45
b. Débits par Mallet Gauthier : ................................................................................................. 46
c. Débits par Hazan Lazareviç : ................................................................................................. 46
II. MODELISATION DES ECOULEMENTS HYDRAULIQUES : ............................................................... 49
1. Intégration des outils HEC-RAS, HEC-GeoRAS et ArcGIS dans la modélisation hydraulique : 49
a. Utilité d’ArcGIS : .................................................................................................................... 49
b. Utilité de HEC-GeoRAS : ........................................................................................................ 49
c. Utilité de HEC-RAS : .............................................................................................................. 49
2. Construction de la géométrie hydraulique à l’aide de HEC-GeoRAS : ..................................... 49

2
 a. Préparation du TIN : .............................................................................................................. 50
b. L’extraction de la géométrie hydraulique :........................................................................... 50
3. Simulation : ............................................................................................................................... 52
4. Résultats : .................................................................................................................................. 55
ANEXES .................................................................................................................................................. 57
CONCLUSION : ....................................................................................................................................... 61

3
 LISTE DES FIGURES :

Figure 1: organigramme du L'ABHT ....................................................................................................... 10

Figure 2 : situation géographique du bassin versant Tensift ................................................................. 11
Figure 3 : situation géographique du sous bassin versant Ourika ......................................................... 13
Figure 4 : réseau hydrographique du sous bassin versant Ourika ......................................................... 15
Figure 5 : Carte de vegetation ............................................................................................................... 17
Figure 6 : carte de facteur d’erosion...................................................................................................... 18
Figure 7 : précipitations annuelles en mm ............................................................................................ 19
Figure 8 : Précipitations moyennes saisonnières en mm ...................................................................... 20
Figure 9 : Précipitations moyennes mensuelles en mm ........................................................................ 21
Figure 10 : Précipitations journaliers maximales en mm ...................................................................... 22
Figure 11 : Débits moyens annuels en m3/s ......................................................................................... 23
Figure 12 : Débits moyens saisonniers en m3/s .................................................................................... 24
Figure 13 : Débits moyens mensuels en m3/s ....................................................................................... 24
Figure 14 : carte de pente...................................................................................................................... 29
Figure 15 : la carte hypsométrique........................................................................................................ 30
Figure 16 : Courbe hypsométrique du sous bassin versant d'Ourika .................................................... 32
Figure 17 : Ajustement graphique des débits maximaux selon les lois statistiques.............................. 40
Figure 18 : Hydrogramme simulé du sous-bassin ourika – Précipitation, pertes et débit de
ruissellement ......................................................................................................................................... 44
Figure 19 : débit de pointe retenue....................................................................................................... 44
Figure 20 : Model de terrain sous forme de triangles irregiliers ........................................................... 50
Figure 21 : segmen trepresnte le flowpath ........................................................................................... 50
Figure 22 : Extraction des sections transversales à partir du TIN à l’aide de HEC-GeoRAS ................... 51
Figure 23 : importation des fichiers dans hec ras .................................................................................. 52
Figure 24 : profile en travers de la simulation ....................................................................................... 53
Figure 25 : Extraction des sections transversales à partir du TIN à l’aide de HEC-GeoRAS ................... 53
Figure 26 : insertion le debit.................................................................................................................. 54
Figure 27 : insertion du pente amont et pente aval du bassin .............................................................. 55
Figure 28 : Modélisation de l'écoulement sur le tronçon étudié .......................................................... 55
Figure 29 : modélisation hydrologique sur hyfran plus pour obtient XT pour la modélisation
hydraulique ............................................................................................................................................ 60

4
 LISTE DES TABLEAUX :

Tableau 1 : Cordonnées d’exutoire (station aghbalou) .......................................................................... 14

Tableau 2: Répartition de la superficie selon l’élévation des courbes de niveau .................................. 31
Tableau 3 : résumé des caractéristiques topographiques ..................................................................... 34
Tableau 4 : résultats de calcul de temps de concentration en h ........................................................... 37
Tableau 5 : PJ max retenue, BIC, AIC, pour T=100 ans .......................................................................... 41
Tableau 6 : Pj max retenue (hyfran plus) ............................................................................................... 41
Tableau 7 : paramétrage du model........................................................................................................ 43
Tableau 8 : paramétres regionaux a,N de fuller 2 ................................................................................. 45
Tableau 9 : paramètres a,N utilisées ..................................................................................................... 45
Tableau 10 : débit calculé par Fuller 2 ................................................................................................... 46
Tableau 11 : paramétres utilisées de mallet gautier ............................................................................. 46
Tableau 12 : débit calcule par mallet gautier ........................................................................................ 46
Tableau 13 : paramètres régionaux d’hazan Lazarevic .......................................................................... 47
Tableau 14 : paramètres d’hazan Lazarevic utilises ............................................................................... 48
Tableau 15 : débit calcule par hazan Lazarevic ...................................................................................... 48
Tableau 16 : Les résultats de la simulation sous format du tableau...................................................... 56
Tableau 17 : données pluviometriques (station aghbalou) ................................................................... 58
Tableau 18 : donnees hydrometriques .................................................................................................. 59

5
 REMERCIEMENT :

Au terme de notre stage de fin de formation en option Hydraulique à l’Institut

Spécialisé des Travaux Publics de Marrakech (ISTPM), nous souhaitons exprimer nos plus
sincères remerciements à toutes les personnes qui ont contribué, de près ou de loin, à la
réussite de cette expérience professionnelle.
Nous tenons à adresser notre profonde gratitude à Madame RAJI ILAH Asma, notre
encadrante pédagogique, pour son accompagnement rigoureux, sa disponibilité et la qualité
de ses conseils tout au long de la réalisation de ce stage, mené en collaboration avec l’Agence
du Bassin Hydraulique du Tensift (ABH Tensift).
Nos remerciements s’adressent également à Madame INTISSAR, pour son aide
précieuse dans la recherche et l’obtention de ce stage, ainsi que pour son engagement dans
l’encadrement et le suivi des étudiants de notre filière.
Nous remercions chaleureusement l’ensemble du personnel de l’ABH Tensift, pour leur
accueil, leur disponibilité et leur collaboration durant toute la période du stage. Leur
encadrement et leurs conseils techniques nous ont permis d’enrichir nos connaissances et de
mieux appréhender les réalités du terrain.
Nous exprimons également notre reconnaissance à l’ensemble du corps enseignant de
l’ISTPM Marrakech, pour la qualité de leur formation, ainsi qu’à l’administration, représentée
par Monsieur BOUTAGLAY, directeur de l’institut, pour son soutien indéfectible, sa
disponibilité et sa gestion humaine, qui ont grandement contribué à créer un environnement
propice à notre apprentissage et à notre réussite.
À toutes ces personnes, nous adressons nos remerciements les plus sincères.
Enfin, nous tenons particulièrement à remercier nos familles pour leur encouragement
constant et leur patience, sans lesquels ce travail n’aurait pas été possible, ainsi que toutes
les personnes qui nous ont soutenus tout au long de ce parcours.

6
 INTRODUCTION :
L’eau constitue l’une des ressources naturelles les plus importantes, autour de laquelle
s’articulent toutes les activités humaines. Toutefois, cette ressource essentielle fait face à de
nombreux défis, notamment des problèmes de gestion, de pénurie et d’inondations,
exacerbés par les changements climatiques.
Au Maroc, les phénomènes hydrologiques extrêmes ont longtemps retenu l’attention
des chercheurs et des autorités. En effet, le pays subit régulièrement des crues et inondations
dans diverses régions, comme en témoignent les événements survenus dans la vallée du Ziz
en 1965, à Ourika en 1995, dans la région d’El Hajeb en 1997, ou encore à Taroudant et
Errachidia en 2019. Ces crues, bien que considérées comme exceptionnelles, se révèlent
violentes et destructrices, impactant les vies humaines, les infrastructures, les terrains
agricoles et les plaines en aval. Ainsi, leur prévision et leur prévention demeurent des enjeux
majeurs.
Notre zone d’étude se situe dans le Haut Atlas, plus précisément au niveau du bassin
versant d’Ourika à la station AGHBALOU. Cette région montagneuse, caractérisée par un climat
semi-aride à subhumide, est particulièrement exposée aux risques hydrologiques.
Ce travail s’inscrit dans le cadre d’un stage professionnel réalisé en collaboration avec
l’Agence du Bassin Hydraulique du Tensift (ABHT). Son objectif principal est la modélisation
hydrologique du bassin versant d’Ourika à la station AGHBALOU, en utilisant une approche
combinant les logiciels ArcGIS et HYFRAN Plus.
La première partie de cette étude est consacrée à la présentation du contexte général
de la zone d’étude, incluant ses caractéristiques géographiques, géologiques et climatiques,
ainsi que la délimitation du bassin versant et une analyse pluviométrique et hydrométrique.
La deuxième partie aborde l’étude hydrologique proprement dite, en examinant les
caractéristiques géométriques et topographiques du bassin, le calcul du temps de
concentration et l’estimation des débits de pointe.

7
8
 I. CONTEXTE INSTITUTIONNEL : L’ABHT ET LA GESTION DE
L’EAU :
L’Agence du Bassin Hydraulique de Tensift (ABHT) est instaurée par le décret N° 2-00-479 du
14 Novembre 2000. Elle a pour mission d’évaluer, de planifier, de développer et de gérer les
ressources en eau au niveau de sa zone d’action. Elle traduit la politique de décentralisation
de la gestion de l’eau induite par le nouveau concept de la gestion intégrée des ressources en
eau à l’échelle du bassin hydraulique, et constitue ainsi un organisme fédérateur des acteurs
concernés par la gestion de l’eau au niveau régional.
L’Agence du bassin est chargée des missions suivantes :

✓ Elaborer un plan directeur d’aménagement intégré des ressources en eau relevant de sa

zone d’action.
✓ Veiller à l’application du plan directeur d’aménagement intégré des ressources en eau à
l’intérieur de sa zone d’action.
✓ Assurer l’approvisionnement en eau en cas de pénurie d’eau déclarée.
✓ Gérer et contrôler l’utilisation des ressources en eau mobilisées.
✓ Réaliser les infrastructures nécessaires à la prévention et à la lutte contre les
inondations.
✓ Réaliser toutes les mesures piézométriques et de jaugeages ainsi que les études
hydrologiques, hydrogéologiques, de planification et de gestion de l’eau tant au plan
quantitatif que qualitatif.
✓ Fournir toute aide financière et toute prestation de service, notamment d’assistance
technique, aux personnes publiques ou privées qui en feraient la demande, soit prévenir
la pollution des ressources en eau, soit en vue d’un aménagement ou d’une utilisation du
domaine public hydraulique.

II. OBJECTIFS DE L’ABHT :

✓ Protection et conservation quantitative et qualitative des ressources en eau.

✓ Utilisation rationnelle et optimale de l’eau.
✓ Gestion intégrée, décentralisée et forcée des ressources en eau et du domaine public
hydraulique.
✓ Planification cohérente et obligatoire à l’échelle du bassin hydraulique.

9
III. ORGANIGRAMME DU L’ABHT :

Figure 1: organigramme du L'ABHT

10
 IV. SITUATION GEOGRAPHIQUE DU BASSIN VERSANT TENSIFT.
Le bassin Tensift est une aire située au centre ouest du Maroc, s’étendant de la ligne de partage
des eaux de surface d’Oum R’bia au nord et à l’est, aux crêtes du haut atlas au sud, l’océan
atlantique à l’ouest. Il occupe 3% de la superficie totale du pays.

Figure 2 : situation géographique du bassin versant Tensift

11
12
 I. SITUATION GEOGRAPHIQUE DU SOUS BASSIN VERSANT
D’OURIKA A LA STATION AGHBALOU
L’oued d’Ourika fait partie du grand bassin du Tensift, il se situe à 31° et 31°20’ de latitude
Nord et à 7°30’ et 8° de longitude. Son bassin versant à Aghbalou se situe dans le Haut Atlas
occidental au centre de la province d’Al Haouz à 35 km de Marrakech. Il se trouve entre les
bassins versants de la Ghiraya et du N’fis. Le bassin est un hydro système bien individualisé. Il
se répartit principalement sur trois Communes rurales : Ourika, Setti Fadma, et l’Oukaïmden .

Figure 3 : situation géographique du sous bassin versant Ourika

13
 II. LA DELIMITATION DU BASSIN VERSANT D’OURIKA AGHBALOU
1. Présentation du logiciel Arc GIS : Fonctionnalités et Composants
Arc GIS est la suite logicielle SIG produite par ESRI (Environmental Systems Research Institute).
C’est un système complet qui permet de collecter, organiser, gérer, analyser, communiquer
et diffuser des informations géographiques.
Arc GIS comprend une suite d’applications intégrées les unes aux autres
❖ Arc Catalog facilite l’organisation et la gestion de toutes les données SIG (cartes, jeux
de données, modèles…)
❖ Arc Globe et Arc Scene sont spécialisées dans l'analyse et la visualisation en trois
dimensions.
❖ Arc ToolBox est Composée d’un ensemble d’outils de transformation des données.
❖ Arc Map est l'application principale dans Arc GIS. Il permet l'édition, la modification
de données géographiques, l'analyse, la conception et l'édition de cartes. Il représente
les informations géographiques sous forme d'un ensemble de couches et les autres
éléments sous forme d'une vue de carte. Il y a deux vues cartographiques principales
dans ArcMap : le mode donné et le mode mise en page.
- Le mode donné présente une fenêtre géographique, ou bloc cartographique pour afficher et
utiliser les informations géographiques sous forme d'une série de couches.
- Le mode Mise en page affiche une page pour organiser les éléments cartographiques (par
exemple, le bloc de données, une barre d'échelle et un titre de carte) pour l'impression de la
carte.

2. Les cordonnées Lambert et de l’exutoire du bassin versant

Types
Cordonnées Lambert zone 1 X= 276459.00 Y= 82616.00
Géographiques LAT= 31.314079 LONG= -7.74648

Tableau 1 : Cordonnées d’exutoire (station aghbalou)

14
 3. Le réseau hydrographique :

²
Ourika Watershed Hydrographic Map

Legend
! station aghbalou
ourika watershed
Main watercourse

0 2 4 8
Miles

Figure 4 : réseau hydrographique du sous bassin versant Ourika

Le réseau hydrographique est l’ensemble des ruisseaux de drainage permanents ou

temporaires, par ou s’écoulent les eaux de ruissellement ou celles restituées par les nappes
souterraines sous forme de sources.

15
III. CADRE GEOLOGIQUE :
Le bassin versant d’Ourika à la station AGHBALOU montre deux grands types de faciès :
Terrains perméables Terrains semi-perméable Terrains imperméable
6% 39 % 55
❖ La partie montagne :
• Située à des altitudes supérieures à 2000 m.
• Constituée par :
-Des roches magmatiques
-Des roches métamorphiques (Gneiss, migmatites)
-Des roches plutoniques (Granites, Granodiorites)
-Des roches volcaniques (Andésites, Rhyolite)
❖ La partie plaine :
• Située à des altitudes inférieures à 2000 m.
• Composée par des dépôts permo-triasiques et quaternaires plus tendre.
• Comprendre au Nord un faciès formé de Conglomérats, Grés et Silites. Au sud un
facies formé essentiellement de Silites argileuses et localement de Grés massif.

La lithologie du bassin versant est donc dans l’ensemble assez imperméable selon les
chiffres de l’Agence du Bassin Hydraulique Tensift. (Saidi et al. 2006 ; Fniguire et Boukrim,
2007)

16
 V. CARTOGRAPHIE DE LA COUVERTURE VEGETALE ET DE
FACTEUR D'EROSION HYDRIQUE DU BASSIN VERSANT :
1. La végétation :
Le bassin présente une couverture végétale insuffisante, ce qui influence directement la
dynamique des écoulements en les accélérant. En effet, une végétation dense joue un rôle
essentiel dans la régulation hydrologique : elle limite les ruissellements rapides et réduit
l’intensité des crues en favorisant l’infiltration de l’eau. Or, la végétation actuelle se concentre
principalement dans quelques zones boisées, dominées par des chênes verts, ainsi que dans
des vergers de pommiers et de noyers situés en fond de vallée.

Station aghbalou

Figure 5 : Carte de vegetation

17
 2. L’érosion du bassin versant
La carte du facteur d'érosion révèle une forte variabilité spatiale, avec des valeurs allant de
0.55 (faible) à 0.9 (élevé). Les zones à haut risque (en rouge) correspondent aux pentes
abruptes et aux sols peu végétalisés, favorisant le ruissellement. Les secteurs en vert, moins
érosifs, coïncident avec les forêts et vergers. Cette analyse identifie les priorités pour la
conservation des sols, notamment via la restauration du couvert végétal dans les zones
critiques (ex : T1 à T10).

aghbalou

Figure 6 : carte de facteur d’erosion

18
IV. LE REGIME PLUVIOMETRIQUE ET HYDROMETRIQUE DU
BASSIN VERSANT D’OURIKA A LA STATION AGHBALOU :
1. Le régime pluviométrique
Les précipitations sont essentiellement mesurées à la station D’AGHBALOU, localisée à
l’exutoire du bassin versant de l’Ourika.
a. Précipitations annuelles (mm) :

PRECIPITATIONS ANNUELLE
1200.0

1000.0

800.0

600.0

400.0

200.0

117,7

0.0

Figure 7 : précipitations annuelles en mm

La courbe montre une forte irrégularité intra annuelle des pluies avec une valeur maximale
(1053,5 mm) en 1970 et une valeur minimale (117,7 mm) en 1992, ce qui traduit
l’appartenance de la zone de l’étude au climat semi-aride.

19
 b. Précipitations moyennes saisonnières

PLUIES MOYENNES SAISONNIEES (mm)

250.00
208.30
200.00
166.7

150.00 128.30

100.00

50.00 28.7

0.00
AUTOMNE HIVER PRINTEMPS ETE

Figure 8 : Précipitations moyennes saisonnières en mm

L’histogramme ci-dessous illustre une augmentation progressive des précipitations, débutant

en automne et atteignant un pic au printemps. En revanche, la saison estivale enregistre les
niveaux les plus faibles. La répartition saisonnière des pluies peut être structurée comme
suit :
- Une période fortement pluvieuse : correspondant aux saisons d’hiver et de
printemps, avec des précipitations variantes entre 166,7 mm et 208,3 mm .
- Une période modérément pluvieuse : représentée par l’automne, avec un cumul de
128,3 mm .
- Une période sèche : observée durant l’été, avec un minimum de 28,7 mm.

20
 c. Précipitations moyennes mensuelles :

PLUIE MOYENNE MENSUELLE ( mm)

90

80

70

60

50

40

30

20

10

Figure 9 : Précipitations moyennes mensuelles en mm

La courbe des précipitations mensuelles moyennes à Ourika révèle une moyenne annuelle
d’environ 44,3 mm. On observe une augmentation des hauteurs de pluie dès le mois de
septembre, atteignant un premier pic en novembre. Une légère baisse survient en décembre,
suivie d’une reprise jusqu’au mois de mars, où les précipitations atteignent 82,1 mm Ensuite,
les valeurs diminuent progressivement jusqu’à atteindre un minimum en juillet (5,4 mm),
avant de remonter à nouveau. Cette évolution traduit une alternance marquée entre des mois
humides et des mois secs tout au long de l’année.

21
 d. Précipitations journaliers maximales :

[Link] (mm)
100.0

90.0

80.0

70.0

60.0

50.0

40.0

30.0

20.0

10.0

0.0

Figure 10 : Précipitations journaliers maximales en mm

Le graphique montre une forte variabilité annuelle des précipitations journalières maximales
entre 1968 et 2021.
-On observe plusieurs pics importants, notamment en 1968, 1990 ou 2005 (le plus élevé, >90
mm), témoignant d’épisodes pluvieux extrêmes.
Certaines années comme 1992, 2000 ou 2011 enregistrent des valeurs très faibles.
Globalement, il n’y a pas de tendance claire (hausse ou baisse), mais la fréquence
d'événements extrêmes semble plus marquée après 2000, ce qui pourrait être lié au
changement climatique.

22
2. Le régime hydrométrique
a. Débits moyens annuels

Débits moyens annuels (m³/s)

12

10

Figure 11 : Débits moyens annuels en m3/s

Les débits moyens annuels présentent des variations interannuelles notables, marquées
par l’alternance d’années sèches et humides. La station D’AGHBALOU a enregistré des débits
particulièrement élevés durant l’année hydrologique 1970/1971, avec un maximum atteint de
10,4 m³/s.

23
 b. Débits moyens saisonniers :

DEBITS MOYENNES SAISONNIERS (m³/s)

8.00

7.00

6.00

5.00

4.00

3.00

2.00

1.00

0.00
AUTOMNE HIVER PRINTEMPS ETE

Figure 12 : Débits moyens saisonniers en m3/s

L’histogramme ci-dessous présente les débits moyens saisonniers de l’Ourika à la station

D’AGHBALOU sur la période allant de 1970 à 2017. Les faibles débits sont observés durant les
saisons d’été, d’automne et d’hiver, tandis que les débits les plus élevés sont enregistrés au
printemps
c. Débits moyens mensuels :

D É BI T MOY E N ME N SUE LLES (M³/ S)

10

Figure 13 : Débits moyens mensuels en m3/s

24
La courbe montre que les débits moyens mensuels de l’Ourika sont élevés en mars,
avril et mai, tandis que les autres mois présentent des débits nettement plus faibles.

25
26
 I. CARACTERISTIQUES GEOMETRIQUES :

1. Surface (A) :
Il s'agit de la surface en plan du bassin versant. Elle est généralement déterminée soit à
l'aide d'outils informatiques, soit de manière approchée suivant un principe de quadrillage .

S = 503 Km2

2. Périmètre (P) :
Le périmètre est la longueur (km) de la ligne de contour du bassin. La mesure se faite à
l’aide d’un curvimètre ou par des logiciels.

P = 122 Km

3. Forme :
L’indice de Gravelius (ou indice de compacité) est un indicateur clé en hydrologie pour
évaluer la forme d’un bassin versant et son influence sur le comportement hydrologique.
L’indice de compacité défini comme le rapport du périmètre du bassin au périmètre du cercle
ayant la même surface :

KG = 0,28*P/ξ𝑺 = 𝟏, 𝟓𝟐
Avec :
• KG : Indice de compacité de Gravelius
• A : Surface du bassin versant (km²)
• P : Périmètre du bassin (km)
On a KG = 1,52 on peut conclure que le bassin d’étude présente une forme allongée

27
 1. La longueur de la cour d’eau principale (Li) :
La longueur de la cour d’eau principale renseigne sur le temps de transit de l’eau pour
atteindre l’exutoire.

Li = 50,76 km

2. Le rectangle équivalent :
Est une Représentation simplifiée d'un bassin sous forme rectangulaire, conservant sa
surface, son périmètre et son indice de Gravelius. Permet d'analyser et comparer les temps de
concentration et les écoulements de différents bassins de manière standardisée.
Utilité :
- Comparaison objective entre bassins
- Modélisation hydrologique plus aisée
- Estimation des temps de réponse aux précipitations

• Longueur : L=0.25*P*[1+ ((1-1.12/KC) ^2) ^1/2] = 51,5

• Largueur :
L=0.25*P*[1- ((1-1.12/KC) ^2) ^1/2] = 9,76

28
 II. CARACTERISTIQUES TOPOGRAPHIQUES

1. La pente :
a. La pente moyenne :
La pente moyenne est une caractéristique importante dans le ruissèlement. Elle donne une
bonne indication sur le temps de parcours du ruissellement direct et donc sur le temps de
concentration. Elle influence directement le débit de pointe lors d'une averse.

Im = ΔH /L = 5,93

Avec :
• Im : La pente moyenne
• ΔH : Altitude max – Altitude min [m].
• L : Longueur du cours d’eau principal en [m].

b. La carte de pente :

Figure 14 : carte de pente

29
 2. Hypsométrie :
a. La carte hypsométrique
La carte hypsométrique donne une idée sur les proportions des tranches d’altitudes. Les
études Une carte hypsométrique en hydrologie représente la relation entre l’altitude et la
surface cumulée du bassin versant. Elle montre comment la surface du terrain se répartit selon
les altitudes, ce qui aide à analyser le relief, l’érosion potentielle, et le comportement
hydrologique du bassin.
1176000 ,000000 1184000 ,000000 1192000 ,000000 1200000,000000 1208000,000000 1216000,000000
,000000

3490000,000000
Hypsometric map of the Ourika watershed

²
3490000

,000000

3483000,000000
3483000

,000000

3476000,000000
3476000

Legend
Elevations
,000000

3469000,000000
998 - 1 470
3469000

1 470 - 1 778
1 778 - 2 061
2 061 - 2 357
,000000

2 357 - 2 640 3462000,000000

3462000

2 640 - 2 924
2 924 - 3 208
3 208 - 3 491
,000000

3455000,000000

0 2,5 5 10 3 491 - 4 012

Km
3455000

1176000 ,000000 1184000 ,000000 Figure1192000

15 : la carte hypsométrique
,000000
1200000,000000 1208000,000000 1216000,000000

La carte hypsométrique fournit une vue synthétique sur la pente du bassin. Donc sur le
relief cette courbe représente la répartition de la surface du bassin versant en fonction de son
altitude. Il est compris entre 981 m et 4011 m.

30
b. La courbe Hypsométrique.
La courbe hypsométrique est construite en calculant les surfaces situées entre les
courbes de niveau, puis en exprimant ces surfaces en pourcentage de la surface totale du
bassin. On représente ensuite un graphique avec les altitudes en ordonnée et les surfaces
cumulées correspondantes en abscisse.

Altitudes Surface
En Km2 En % Surface cumule en %
998 - 1 470 33 6,57 100,00

1 470 - 1 778 56 11,16 93,43

1 778 - 2 061 70 13,94 82,27

2 061 - 2 357 72 14,34 68,33

2 357 - 2 640 69 13,75 53,98

2 640 - 2 924 67 13,35 40,24

2 924 - 3 208 55 10,96 26,89

3 208- 3 491 47 9,36 15,94

3 491 - 4 012 33 6,57 6,57

Tableau 2: Répartition de la superficie selon l’élévation des courbes de niveau

31
 Surfaces cumulées [%] COURBE HYPSOMETRIQUE
Min : 998
100

90

80

70

60
Médiane : 2400
50

40

30

20

10 Altitudes
Max : 4012
0
1217

3289
975
1097

1340
1462
1583
1706
1828
1948
2071
2193
2313
2436
2558
2678
2801
2923
3044
3167

3409
3532
3654
3774
3897
Figure 16 : Courbe hypsométrique du sous bassin versant d'Ourika

La courbe montre une élévation rapide des altitudes vers l’amont, elle donne une idée sur le
profil du bassin versant. Par ailleurs, l’examen de la courbe permet de déduire :
• L'altitude maximum : 4012 m (correspondant au sommet du bassin versant)
• L'altitude minimale : 981 m
• L'altitude médiane correspond à une surface cumulée de 50% : 2400 m
-L'altitude moyenne correspond à l’altitude lue au point d’abscisse 50 % de la surface totale
du bassin, sur la courbe hypsométrique. Elle est donnée par l’équation suivante.

32
 c. Indice de pente IG :
Cet indice permet de déterminer l’importance du relief sur le bassin comme étant le rapport
entre le dénivellement utile (DU) et la longueur L du rectangle équivalent cet indice
caractérise le relief du bassin.

Ig = DU/Léq = H5% - H95%/Léq

• DU : Dénivèlement utile
• Léq : longueur de rectangle équivalent
• H5% : Altitude correspondant à 5% de la surface total du bassin.
• H95% : Altitude correspondant à 95% de la surface du bassin
• D’après la courbe hypsométrique : H5% = 3747 m ; H95%= 1247 et on a Léq = 51,5 km

Donc : Ig = 48,54

33
 Paramètres Valeur

Surface 503 km²

Périmètre 122 km

Indice de compacité 1.52

Longueur du cours d’eau principale 50,76 km

Altitude maximale 4011 m

Altitude minimale 981 m

Altitude médiane 2445 m

Altitude moyenne 2160 m

Dénivelée 1464 m

Longueur du rectangle équivalant 51.5 km

Largeur du rectangle équivalant 9.76 km

Pente moyenne 4.57 %

Tableau 3 : résumé des caractéristiques topographiques

34
III. CALCUL LE TEMPS DE CONCENTRATION

Le temps de concentration (Tc) est le maximum de durée nécessaire à une goutte d'eau pour
parcourir le chemin hydrologique entre un point situé le plus loin possible de l’exutoire pour
atteindre ce dernier. Ce paramètre est déterminé par des formules empiriques.

1. Formule espagnole :

Tc=60*0,3*((L/P^0,25) ^0,77) = 10,63

Avec
Tc : Temps de concentration en h
L : longueur du cours d’eau principal en km
P : Pente moyenne du bassin en m/m

2. Formule de Californienne :

Tc=60*0,1452*(L/P^0,5) ^0,77 = 8,86

Avec :
Tc : Temps de concentration en h
L : longueur du cours d’eau principal en km
P : Pente moyenne du bassin en m/m

3. Formule de Passini :

Tc=64,8*(S*L) ^0,333*P^-0,5 =13,01

Avec :
Tc : Temps de concentration en (h)
S : Superficie du bassin versant en (km²)
L : longueur du cours d’eau principal en km
P : Pente moyenne du bassin en %

35
 4. Formule de Kripich:

Tc=0,0195*L^0,77*P^-0,385 = 4,05

Avec :
Tc : Temps de concentration en (h)
L : longueur du cours d’eau principal en m
P : Pente moyenne du bassin en m/m

5. Formule de Turazza :

Tc=60*0,108(S*L) ^0,333*P^-0,5 =13,01

Avec :
Tc : Temps de concentration en (h)
S : Superficie du bassin versant en (km²)
L : longueur du cours d’eau principal en km
P : Pente moyenne du bassin en m/m

6. Formule de Ventura

Tc=76,3*(S/P)^0,5 = 11,71

Avec :
Tc : Temps de concentration en (h)
S : Superficie du bassin versant en (km²)
L : longueur du cours d’eau principal en km
P : Pente moyenne du bassin en %

36
 Formule Temps de concentration

Formule de Kripich 4,05

Formule de Turazza 11,71

Formule de Passini 13,01

Formule espagnole 10,63

Formule de Californienne 8,86

Formule de ventura 11,71

Tableau 4 : résultats de calcul de temps de concentration en h

`
Le temps de concentration retenu est celui de la formule de Kirpich (4,05 h), conformément
aux recommandations du Guide d’assainissement routier 2022.

37
V. CALCUL DES DEBITS DE POINTE :
1. Ajustement statistique :
La méthode utilisée se base sur les données des débits instantanés maximaux annuels de
l’oued Ourika. Grâce au logiciel Hyfran Plus, plusieurs lois statistiques ont été testées (comme
Gumbel, Log-Normal, etc.) pour ajuster les données.
Ensuite, on a comparé les résultats de chaque loi en utilisant deux critères importants : le
critère d'information d'Akaike (AIC) et le critère bayésien (BIC). Ces critères permettent de
choisir la meilleure loi en donnant un score : plus la valeur est faible, meilleur est
l’ajustement.
La loi qui a obtenu les valeurs les plus basses en AIC et BIC a été retenue comme la plus
représentative pour modéliser les débits extrêmes. Cette méthode permet d’avoir des
résultats fiables pour mieux comprendre le comportement de l’oued et pour bien
dimensionner les ouvrages hydrauliques.

a. Présentation du logiciel Hyfran-plus :

Le logiciel Hyfran-Plus permet d’ajuster un large éventail de distributions statistiques à
une série de données hydrologiques, à condition que celles-ci respectent certaines hypothèses
fondamentales : l’indépendance, l’homogénéité et la stationnarité.
Pour faciliter le choix de la distribution la plus adaptée, un Système d’Aide à la Décision (SAD)
intégré au logiciel oriente l’utilisateur vers la classe de lois la plus appropriée, notamment
pour l’estimation des quantiles associés à des périodes de retour élevées, utiles dans la
gestion des risques hydrologiques.
Dans notre cas, la courbe de distribution expérimentale issue des valeurs maximales
d’intensité a été confrontée à plusieurs lois statistiques. L’objectif était de déterminer les
quantiles de crue avec la plus grande précision possible, en s'appuyant sur les résultats des
tests d’ajustement.
b. Lois d’ajustement :
Il existe plusieurs lois (modèles) d’ajustement citant à titre d’exemple :
Loi fweibul
Loi des extrêmes généralisée (GEV)
Loi log- normale (Gualton)
Loi de Gumbel
Loi gamma

c. L’analyse des graphes :

Les figures ci-dessous illustrent les graphiques des résultats obtenus, en comparant les
probabilités empiriques aux probabilités théoriques pour différentes lois statistiques : Gamma,
Log-Normale à 3 paramètres, Weibull, Gamma généralisée, fweibul et Gumbel.

38
À la lecture de ces graphiques, on observe que les lois Weibull et Gamma offrent les meilleurs
ajustements aux données observées du débit instantané maximal, traduisant une bonne
adéquation entre le modèle statistique et le comportement réel de la variable hydrologique
étudiée.

39
 d. L’analyse des critères statistiques :
Le tableau ci-dessous synthétise les valeurs des critères d’information Akaike (AIC) et
bayésien (BIC).

Figure 17 : Ajustement graphique des débits maximaux selon les lois statistiques

40
Sur la base des résultats numériques, c’est la loi de Weibull qui minimise le plus ces critères,
ce qui en fait le modèle le plus adapté pour représenter le débit instantané maximal.

Modèle XT BIC AIC

Webull (Method of moments) 82,594 394,741 390,999

GEV (Method of moments) 84,619 398,122 392,508

Gamma (Method of moments) 90,132 395,56 391,817

Lognormal (Maximum 97,259 397,316 393,573

Likelihood)
Gumbel (Method of moments) 96,14 402,633 398,891

Tableau 5 : PJ max retenue, BIC, AIC, pour T=100 ans

Le modèle de Weibull a été retenu car il présente les valeurs les plus faibles de BIC (394,741)
et d’AIC (390,999), indiquant ainsi le meilleur ajustement aux données parmi les distributions
testées. Ce critère statistique confirme sa pertinence pour modéliser les événements
extrêmes analysés.

Débit de pointe retenue :

[Link]

T (ans) Gumbel lognormal gev GAMMA weibull

2 50,6 51 52.6 51,7 53,4

5 62,8 64,4 64,8 64 64,8

10 70,9 72,8 71,1 71,1 70,4

20 78,6 80,5 76,1 77,4 74,8

50 88,6 90,2 81,4 84,9 79,5

100 96,1 97,3 84,6 90,1 82,594

Tableau 6 : Pj max retenue (hyfran plus)

41
2. Transformation pluie – débit :
La transformation pluie-débit est une étape clé de la modélisation hydrologique. Elle
permet de convertir la pluie efficace, issue d’un événement pluvieux, en débit à l’exutoire du
bassin versant. Cette opération repose sur un modèle hydrologique, souvent conceptuel ou
déterministe, capable de représenter les processus de ruissellement et de concentration. Elle
offre une vision précise de la réponse du bassin, utile pour la gestion des ressources en eau, la
prévention des crues et l’aménagement du territoire.

a. Présentation du modèle utilisé :

Pour simuler la réponse hydrologique du bassin versant de l’Ourika, nous avons utilisé le
modèle SCS (Soil Conservation Service). Ce modèle, développé par le service de conservation
des sols des États-Unis, repose sur le concept de pluie efficace, calculée à partir des
précipitations totales en tenant compte des pertes initiales (infiltration, interception, etc.).
Le modèle SCS utilise le paramètre Curve Number (CN), qui reflète les caractéristiques du
sol, de l’occupation du sol et de l’état hydrologique antérieur. Une fois la pluie efficace estimée,
la méthode permet de générer un hydrogramme unitaire basé sur une forme standardisée et
symétrique.
Ce modèle est largement utilisé pour sa simplicité, sa rapidité d’exécution et sa capacité à
fournir des résultats fiables, notamment pour l’estimation des crues en milieu rural ou semi-
urbanisé.

42
 b. Paramétrage du modèle
Les principaux paramètres du modèle SCS ont été déterminés à partir des données
topographiques, pédologiques et d’occupation du sol du bassin versant de l’Ourika. Le tableau
suivant présente les paramètres utilisés :

Paramètre Description Méthode d’obtention Valeur retenue

Superficie du Surface totale du bassin Analyse SIG / carte 503 km²

bassin (A) versant topographique

Curve Number Indice d’humidité du sol Selon la carte 90

(CN) et d’occupation du sol d’occupation du sol

Pluie nette (Pe) Pluie efficace après 52,81

pertes initiales Pe = (P - Ia) ² / (P - Ia + S)

Pertes initiales (Ia) Infiltration, interception, Ia = 0.2 × S 5,644

stockage de surface S=(25400-254CN)/CN

Temps de Temps nécessaire à l’eau 243,04

Formule empirique (
concentration (Tc) pour atteindre l’exutoire
Kirpich )

Lag time Temps entre pluie 145,82

Lag time = Tc (min) * 0,6
maximale et débit
maximal.

Tableau 7 : paramétrage du model

c. Paramètres clés de sortie :

Cette section présente les résultats principaux issus de la simulation hydrologique du
bassin versant de l’Ourika, réalisée avec le modèle HEC-HMS. Les paramètres hydrologiques
estimés comprennent la superficie drainée, le débit de pointe, l'heure du pic de crue et le
volume ruisselé. Les captures d'écran ci-dessous illustrent ces valeurs, permettant une analyse
précise du comportement hydrologique du bassin sous différents scénarios de précipitation.

43
 Figure 19 : débit de pointe retenue.

La simulation "Run 2" pour le sous-bassin "bv1" fournit un débit de pointe de 1564,1 m³/s,
survenant le 10 mai 2025 à 11h05. Le volume total de précipitation enregistré est de 82,59
mm, dont 26,29 mm ont été perdus (infiltration, évaporation). Le ruissellement direct atteint
53,35 mm, représentant l’ensemble du volume d’écoulement de surface. Le débit de base est
négligeable, avec une valeur de seulement 0,01 mm.

Figure 18 : Hydrogramme simulé du sous-bassin ourika – Précipitation, pertes et débit de

ruissellement

44
Le graphique ci-dessus représente les résultats de la simulation hydrologique effectuée sous
HEC-HMS pour le sous-bassin "bv1" du bassin versant de l’Ourika, correspondant à la
simulation "Run 2". La partie supérieure du diagramme illustre l'événement de précipitation
ainsi que les pertes hydrologiques (notamment l'infiltration). La partie inférieure montre
l’hydrographe de ruissellement résultant.
-On observe un pic de débit atteint vers 14h00 le 10 mai 2025, avec un débit de pointe estimé
à environ 1531,7 m³/s. Ce pic résulte de l’intensité de la pluie et des caractéristiques du bassin.
Le débit diminue ensuite progressivement, traduisant l’évacuation de l’eau de ruissellement.
Ce graphique met en évidence la dynamique de réponse du bassin versant à un événement
pluvieux intense, permettant ainsi d’évaluer les risques de crue et de mieux planifier les
mesures de gestion des eaux.
1. Méthodes empiriques :
a. Débits par Fuller 2 :

Q = (1 + a ·log T) ·(S^0,8 + 8/3· s^0,5) ·4 · N/300

Avec :
QT : débit de pointe de période de retour en m3/s
T : Période de retour en ans
S : superficie du bassin versant en km²
a : Coefficient empirique sans dimension, fonction des caractéristiques climatiques du bassin.

Paramètre Valeur Région

0.7 à 0.8 Humides
a
0.8 à 2 Arides
3 à 3.5 Sahariens
80 Plaines
N
85 Reliefs accidentés
100 Montagnes
Tableau 8 : paramétres regionaux a,N de fuller 2

Paramètres utilisés
S 503,13
N 85
a 2

Tableau 9 : paramètres a,N utilisées

45
 periode de retour 5 10 20 50 100
debit calcule par fuller 2 556,59 696,35 836,1 1120,82 1160,57

Tableau 10 : débit calculé par Fuller 2

b. Débits par Mallet Gauthier :

Q = 2 · k ·log(1 + aH) · S ·(1 + 4 ·log T - log S)^0,5 / L^0,5

Avec
Q : Débit de pointe pour la période de retour T (m3/s)
H : Pluie moyenne annuelle en m
T : La période de retour en ans
S : Superficie du bassin versant en km²
L : Longueur du cours d’eau principal en km
a : Coefficient variant de 20 à 30. On prend 20 pour le Maroc
K : Paramètre régional variant de 0.5 à 5. On prend K= 2 pour le bassin versant étudié.

S 503,13
L 50,76
a 20
K 2
H 0,511
Tableau 11 : paramétres utilisées de mallet gautier

periode de retour 5 10 20 50 100

debit calcule par mallet gautier 310,25 449,64 555,07 669,43 744,35

Tableau 12 : débit calcule par mallet gautier

c. Débits par Hazan Lazareviç :

La formule d’Hazan Lazareviç donne le débit millénial en fonction de la superficie du bassin
versant. Sa formule est comme suit :

46
 Q(1000) = K1 · S^k2

Avec :
S : surface du bassin versant en km².
T : période de retour en ans.
K1 et K2 : paramètres dont les valeurs sont résumées dans les tableaux ci-dessous :

Zone Pluviométrie (mm) X Y

Rif Central 1000 - 1300 15.53 0.776

Rif Occidental 800 - 1000 9.78 0.793

Rif Oriental 600 - 800 7.58 0.808

Moyen Atlas 700 - 900 14.94 0.636

Saharien
Moyen Atlas 500 - 700 13.51 0.613
Saharien
Moyen Atlas 400 - 500 13.47 0.587
Saharien
Haut Atlas Saharien 200 - 400 9.38 0.742

Tableau 13 : paramètres régionaux d’hazan Lazarevic

On s’intéresse de la zone du Haut Atlas Saharien

-Pour calculer les débits on utilise la formule initiale de Fuller

Q(T) = Q(1000) ·(1 + a ·log T)/(1 + a ·log 1000)

Avec :
Q (T) : Débit de pointe de période de retour en (m3/s)
Q (1000) : La formule d’Hazan Lazareviç
a : Coefficient régional, égale 2
T : La période de retour

47
 Paramètres utilisés

X 13,47
Y 0,587
a 2

Tableau 14 : paramètres d’hazan Lazarevic utilises

période de retour 5 10 20 50 100

débit calcule par Mallet 622,77 667,44 712,06 771,08 815,73

Gautier

Tableau 15 : débit calcule par hazan Lazarevic

48
 II. MODELISATION DES ECOULEMENTS HYDRAULIQUES :

La simulation hydraulique a pour vocation de représenter de manière précise le comportement

de l’écoulement de l’eau le long du lit d’un cours d’eau, en tenant compte de la topographie,
des conditions hydrologiques et des caractéristiques géométriques de la zone étudiée. Elle
permet d’analyser les régimes d’écoulement, d’identifier les zones susceptibles d’être
inondées, de vérifier la capacité hydraulique des sections, et d’assister à la conception ou à
l’optimisation d’ouvrages hydrauliques.
Dans le cadre de cette étude, la simulation a pour objectif de produire un modèle fiable du
comportement hydraulique du tronçon étudié, en se basant sur un MNT détaillé et une
géométrie extraite via HEC-GeoRAS. Ce modèle permet d’évaluer différents scénarios
d’écoulement, qu’ils soient en conditions ordinaires ou extrêmes, afin d’anticiper les impacts
hydrauliques et d’appuyer les décisions d’aménagement.

1. Intégration des outils HEC-RAS, HEC-GeoRAS et ArcGIS dans la

modélisation hydraulique :
a. Utilité d’ArcGIS :
ArcGIS est un Système d’Information Géographique (SIG) utilisé pour la gestion, la visualisation
et l’analyse des données spatiales. Dans cette mission, ArcGIS permet :
L’affichage et le traitement du TIN (modèle de terrain).
La création et l’édition des entités géométriques (chenal, berges, lignes transversales, etc.).
L’environnement de travail pour HEC-GeoRAS.
b. Utilité de HEC-GeoRAS :
HEC-GeoRAS est une extension d’ArcGIS qui sert de passerelle entre les données géospatiales
et le modèle hydraulique. Son utilité principale est :
L’extraction automatique de la géométrie hydraulique (profils transversaux, lit du cours d’eau,
berges).
La structuration des données spatiales au format requis par HEC-RAS.
La génération de fichiers d’export vers HEC-RAS (.sdf, .xml).
c. Utilité de HEC-RAS :
HEC-RAS (Hydrologic Engineering Center – River Analysis System) est le logiciel de simulation
hydraulique unidimensionnelle. Il est utilisé pour :
Importer la géométrie extraite par HEC-GeoRAS.
Paramétrer les conditions hydrauliques (débits, conditions aux limites).
Simuler les écoulements (permanents ou transitoires).
Analyser les hauteurs d’eau, vitesses d’écoulement et zones inondables.

2. Construction de la géométrie hydraulique à l’aide de HEC-GeoRAS :

Dans cette phase, la géométrie du cours d’eau est extraite à partir du Modèle Numérique de
Terrain (MNT) en utilisant l’extension HEC-GeoRAS sous ArcGIS. Cette étape permet de
générer les entités géométriques nécessaires à la construction du modèle hydraulique dans
HEC-RAS.

49
 a. Préparation du TIN :
La préparation du TIN consiste à générer un modèle de terrain sous forme de triangles
irréguliers à partir de points altimétriques ou courbes de niveau, permettant une
représentation précise du relief pour l'extraction des entités géométriques en modélisation
hydraulique.

Figure 20 : Model de terrain sous forme de triangles irregiliers

Après l’activation de l’extension HEC-GeoRAS dans ArcGIS, le TIN est converti en raster pour
permettre l’extraction précise des entités géométriques nécessaires à la modélisation
hydraulique, telles que le chenal, les berges et les profils transversaux. Ces entités sont ensuite
structurées et exportées vers HEC-RAS. Ce dernier permet de simuler l’écoulement de l’eau,
d’analyser les hauteurs d’eau et d’identifier les zones potentiellement inondables.
b. L’extraction de la géométrie hydraulique :
- On commence par définir un segment représentant le flowpath, c’est-à-dire l’axe
d’écoulement principal de la river.

Figure 21 : segmen trepresnte le flowpath

50
Ensuite :
- Dessin des berges (banks)
Identifiez et dessinez les left bank (berge gauche) et right bank (berge droite) le long
du tracé du flowpath pour délimiter le lit du cours d’eau.
- Création des lignes de coupe transversales (XS cutlines)
Tracez des XS cutlines (profils en travers) perpendiculairement au flowpath. Ces lignes
permettent de définir les sections du lit de la rivière à des intervalles réguliers pour
analyser les variations de terrain.
- Amélioration de la géométrie
Si nécessaire, améliorez la géométrie en ajustant les cutlines pour mieux refléter les
changements de pente et de largeur de la rivière. Vous pouvez également affiner le
tracé du flowpath pour qu’il suive plus précisément le terrain réel.
- Création des zones de calcul (floodplain et autres zones)
Définissez les zones d’écoulement, y compris les floodplains (zones inondables), en
fonction des caractéristiques topographiques extraites du TIN. Cette étape est cruciale
pour la modélisation des inondations.
- Vérification de la topologie et de l’orientation
Assurez-vous que toutes les entités (river, banks, flowpath, cutlines) sont
correctement orientées, avec une direction fluide du tracé de l’amont vers l’aval.
- Génération de la géométrie finale
HEC-GeoRAS génère automatiquement la géométrie hydraulique à partir des entités
créées. Cette géométrie comprend les informations nécessaires pour simuler
l'écoulement de l’eau.
•

Figure 22 : Extraction des sections transversales à partir du TIN à l’aide de HEC-GeoRAS

51
 River : C’est le tracé principal du cours d'eau dessiné manuellement.
• Flowpaths : Représente la direction de l'écoulement dans le cours d’eau.
• Banks : Définissent les berges gauche et droite du lit du cours.
• XT Cutlines : Lignes transversales perpendiculaires au flowpath pour créer des
sections.
3. Simulation :

- L’importation des fichiers de géométrie dans HEC-RAS se fait à l’aide de l’outil « Import
Geometry Data ». Pour ce faire, on sélectionne le format « GIS Format », qui correspond à
l’extension des fichiers exportés par ArcGIS via HEC-GeoRAS. Cette étape permet de charger
toutes les données géométriques cruciales, telles que les sections transversales, les profils en
long et les structures, dans le modèle HEC-RAS.

Figure 23 : importation des fichiers dans hec ras

La partie de la rivière qui sera étudiée est représentée d’un trait bleu fin sur cette figure ainsi que
les points berges sont représentées en points rouges et les profils en travers en ligne vert.

52
Figure 25 : Extraction des sections transversales à partir du TIN à l’aide de HEC-GeoRAS

Figure 24 : profile en travers de la simulation

53
L’image représente la définition d’un profil en travers hydraulique dans HEC-RAS. Les
coordonnées de station et d’élévation permettent de modéliser la géométrie du terrain. Les
berges sont délimitées par les points de type Bank Sta, essentiels pour les calculs
d’écoulement.

Dans ce cas on utilise les Débits de crues (m3/s) qui correspondent aux périodes de retour 100
ans pour notre zone d’étude (1670,087 m3/s).

Figure 26 : insertion le debit

Pour compléter notre fichier d’écoulement, il faut donner des conditions aux limites (Pente en
amant et en aval, une valeur d’un seuil s’il existe).
Le logiciel HEC-RAS nous donne différentes formes de restitution des résultats, sous forme des
tableaux récapitulant l’ensemble des données entrées et des paramètres hydrauliques issus
de la simulation au niveau de chaque point (cross section) du tronçon étudié, ainsi que des
résultats graphiques tel que le profil en long et les profiles en travers…

54
 Figure 27 : insertion du pente amont et pente aval du bassin

4. Résultats :
La figure ci-dessous présente le résultat de la simulation hydraulique réalisée avec HEC-RAS,
correspondant à un débit centennal (Q100ans). On y observe la surface libre de l’écoulement
(en bleu), la topographie du lit du cours d’eau (ligne noire « Ground ») ainsi que les points de
berge (Bank Sta), représentés par les lignes rouges latérales.
Les niveaux d’eau sont représentés sur chaque section transversale, ce qui permet d’analyser
la variation de la ligne d’eau le long du tronçon modélisé. Cette visualisation est essentielle

Figure 28 : Modélisation de l'écoulement sur le tronçon étudié

55
pour vérifier les hauteurs d’eau atteintes, évaluer les zones susceptibles d’être inondées, et
ajuster le modèle si nécessaire.

River Q Total
Min W.S. Crit E.G. Vel Chl Flow Top Froude
Sta (m³/s)Ch Elv Elev WS Elv Slope (m/s) Area Width # Chl
(m) (m) (m) (m²) (m)
561.134 1650.00 723.82 724.86 727.77 0.000988 1.40 1176.00 100.00 0.13
474.404 1670.00 724.30 726.58 736.90 0.00121 1.48 1126.00 100.00 0.14
360.396 1622.00 725.85 726.91 736.93 0.00124 1.51 1118.69 100.00 0.14
359.300 1633.00 726.35 727.64 736.94 0.00163 1.63 994.15 100.00 0.17
304.629 1650.00 726.41 727.90 736.94 0.00191 1.74 976.37 100.00 0.18
230.367 1623.00 727.16 728.38 736.94 0.00238 2.18 805.79 100.00 0.24
203.345 1645.00 727.45 728.62 736.92 0.00325 2.08 805.19 100.00 0.23
148.375 1700.00 728.75 730.04 736.91 0.00452 2.38 683.12 100.00 0.32
118.375 1750.00 729.00 730.59 736.91 0.00513 3.12 535.01 100.00 0.43
116.943 1752.00 729.44 735.79 736.43 0.01747 3.52 497.37 100.00 0.50
0 1709.00 730.45 734.39 735.95 0.007954 5.54 308.63 100.00 1.01

Tableau 16 : Les résultats de la simulation sous format du tableau

Le tableau ci-dessus présente les résultats de la simulation hydraulique obtenus à l’aide

du logiciel HEC-RAS pour un débit centennal (Q100 = 1700 m³/s en moyenne). Plusieurs
paramètres hydrauliques importants sont extraits pour chaque station le long de la rivière
ourika.
On observe une élévation progressive du niveau d’eau (W.S. Elev) depuis l’amont vers
l’aval, ce qui est conforme à la pente naturelle du terrain. Le pente énergétique (E.G. Slope)
varie de 0.000988 à 0.01747, indiquant une augmentation du gradient hydraulique dans les
sections les plus basses, en particulier à l’approche de l’exutoire (station 0), où la vitesse (5.54
m/s) et le nombre de Froude (1.01) traduisent un écoulement potentiellement torrentiel.
Le nombre de Froude reste généralement inférieur à 1, indiquant un écoulement
subcritique (tranquille), sauf dans la dernière section où il atteint 1.01, ce qui marque le
passage à un écoulement critique ou supercritique. La largeur de surface (Top Width) est
maintenue constante à 100 m, ce qui suppose un canal de largeur uniforme ou une hypothèse
simplifiée.
Les vitesses d’écoulement augmentent progressivement (de 1.40 m/s à 5.54 m/s),
corrélées à la diminution de la surface de la section d’écoulement (Flow Area), qui passe de
1176 m² à 308.63 m². Cela reflète un effet de resserrement du lit ou un changement
topographique important.

56
ANEXES

57
ANNEE DATE P. MAX P. NB_MOIS NB_J_AVEC0 NB_J_SANS0 Y
ANNUELLE
1968 19/05/1969 73.5 320.6 5 153 26 1.4
1969 08/04/1970 69.1 363.0 12 365 35 1.3
1970 01/01/1971 63.0 1053.5 12 365 89 1.2
1971 17/01/1972 59.5 802.0 12 366 82 1.2
1972 27/11/1972 35.1 563.3 12 365 71 0.7
1973 10/04/1974 48.6 921.4 12 365 72 0.9
1974 18/04/1975 46.7 499.1 12 365 56 0.9
1975 18/03/1976 32.0 583.8 12 366 68 0.6
1976 19/01/1977 48.9 489.0 12 365 50 0.9
1977 22/11/1977 56.7 647.9 12 365 68 1.1
1978 20/12/1978 37.8 334.4 12 365 52 0.7
1979 09/10/1979 65.2 646.5 12 366 68 1.3
1980 16/03/1981 39.8 407.7 12 365 61 0.8
1981 26/04/1982 66.7 567.7 12 365 59 1.3
1982 21/12/1982 32.2 314.9 12 365 42 0.6
1983 21/03/1984 46.4 376.3 12 366 57 0.9
1984 03/05/1985 55.4 551.8 12 365 67 1.1
1985 28/01/1986 36.0 467.5 12 365 65 0.7
1986 11/02/1987 57.1 313.8 12 365 43 1.1
1987 02/11/1987 42.5 596.1 12 366 75 0.8
1988 27/03/1989 67.2 680.9 12 365 75 1.3
1989 26/12/1989 38.3 421.1 12 365 71 0.7
1990 13/03/1991 71.6 617.1 12 365 76 1.4
1991 13/06/1992 41.5 487.2 12 366 76 0.8
1992 13/03/1993 22.7 276.9 12 365 48 0.4
1993 06/03/1994 64.9 647.8 12 365 70 1.3
1994 10/02/1995 67.8 359.4 12 365 50 1.3
1995 09/01/1996 41.4 725.1 12 366 94 0.8
1996 19/04/1997 39.3 515.2 12 365 67 0.8
1997 27/09/1997 58.4 560.3 12 365 61 1.1
1998 04/10/1998 67.9 552.4 12 365 49 1.3
1999 23/05/2000 36.5 449.4 12 366 54 0.7
2000 12/10/2000 41.9 379.6 12 365 44 0.8
2001 11/04/2002 48.2 481.9 12 365 65 0.9
2002 15/11/2002 49.4 406.1 12 365 72 1.0
2003 16/11/2003 64.1 704.7 12 366 72 1.2
2004 29/10/2004 34.3 303.7 10 304 44 0.7
2005 24/04/2006 86.5 574.1 12 365 66 1.7
2006 11/04/2007 62.5 552.3 12 365 60 1.2
2007 21/03/2008 55.3 418.4 12 366 50 1.1
2008 20/01/2009 60.1 686.3 12 365 91 1.2
2009 16/08/2010 67.8 591.0 12 365 77 1.3
2010 29/04/2011 53.5 679.1 12 365 75 1.0
2011 02/04/2012 32.5 459.8 11 335 52 0.6
2012 05/04/2013 55.2 458.8 12 365 48 1.1
2013 29/01/2014 42.8 455.5 12 365 52 0.8
2014 28/11/2014 67.2 780.8 12 365 62 1.3
2015 15/02/2016 36.0 280.8 12 365 38 0.7
2016 22/11/2016 51.6 369.7 12 365 44 1.0
2017 07/05/2018 38.7 571.1 12 365 63 0.8
2018 18/09/2018 72.4 654.4 12 365 58 1.4
2019 20/03/2020 55.4 440.9 12 366 45 1.1
2020 14/04/2021 41.6 412.2 12 365 43 0.8
2021 04/04/2022 35.1 214.2 12 365 30 0.7
2022 11/10/2022 55.4 117.7 3 91 11 1.1
Tableau 17 : données pluviometriques (station aghbalou)

58
Année Débits moy. j (m³/s) Débit insta (m³/s) Débits moy. j
(mm)
1988-89 268.6 823 46.05
1999-00 250.14 762 42.88
1993-94 138 226 23.06
2014-15 134.84 343.3 23.12
2012-13 132.52 402 22.7
1986-87 109.06 250 18.7
2010-11 89.34 339.1 15.32
1995-96 68.61 165 11.76
2009-10 65.07 234.4 11.15
2008-09 60.79 267 10.42
1970-71 56.06 117 9.61
1973-74 47 77.5 8.06
1994-95 45.84 163 7.86
2000-01 42.56 150 7.82
1992-93 40.4 74.5 7.58
1971-72 40.7 55.2 6.98
2015-16 35.16 200 6.03
1981-82 32.9 93.3 5.64
1996-97 31.87 65 5.46
1972-73 29.7 72 5.09
1987-88 26.19 651 4.49
1975-76 25.6 101 4.39
1984-85 24.6 40.1 4.22
1977-78 23.6 95.6 4.05
2005-06 21.15 208 3.63
2007-08 20.23 558 3.47
2006-07 19.03 358.5 3.26
2011-12 18.79 206 3.22
2003-04 18.4 139 3.16
1976-77 17.9 52.6 3.07
1983-84 17.8 37.4 3.05
1974-75 17 23.4 2.91
1997-98 15.66 8.2 2.68
1982-83 14.6 24.4 2.5
1998-99 11.65 31.4 2
1978-79 10.73 18.4 1.84
2002-03 10.57 11 1.81
2000-01 9.93 26.4 1.64
2004-05 9.58 19.4 1.64
1980-81 7.94 8.1 1.36
1985-86 5.72 50.5 0.98
2013-14 4.18 6.7 0.72
2016-17 2.88 6 0.49
Tableau 18 : donnees hydrometriques

59
Figure 29 : modélisation hydrologique sur hyfran plus pour obtient XT pour la modélisation
hydraulique

60
 CONCLUSION :

Dans le cadre de notre étude hydrologique, nous avons pu analyser le fonctionnement du

bassin versant de l’Ourika grâce à une approche combinée utilisant les logiciels ArcGIS et
HYFRAN+. Cette étude, centrée sur la station d’Aghbalou dans le Haut Atlas occidental, nous
a permis d’aborder les aspects géographiques, climatiques, géologiques, ainsi que l’analyse
pluviométrique, hydrométrique et les caractéristiques physiques du bassin.
Les résultats montrent que la situation géo-climatique particulière du bassin favorise une forte
vitesse d’écoulement et des débits de pointe très élevés. Face à cela, plusieurs mesures
structurelles et non structurelles ont été mises en œuvre dans la vallée de l’Ourika :
aménagement de seuils pour atténuer les crues, stabilisation des thalwegs, construction de
murs de soutènement, ouvrages d’art et d’assainissement routier, ainsi que l’aménagement
d’une piste en crête et de zones de refuge.
Un système de prévision et d’alerte a également été installé en amont pour prévenir les
populations en cas de crue. Le choix de l’Ourika comme bassin pilote montre l’importance
stratégique de cette zone. Nous espérons que ce dispositif pourra être étendu aux bassins
voisins afin de renforcer la résilience de l’ensemble du territoire face aux risques
hydrologiques.

61
62