COMPÉTITIVITÉ

ÉCONOMIQUE DU MAROC
CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET
HYDRAULIQUE
APPORTS ET PRINCIPES D’UTILISATION DES OUTILS HEC, HEC-HMS ET HEC-RAS,
PLATEFORME D’INTEGRATION WMS ET OUTIL RIBASIM

FASCICULE 2 : MODELISATION EN HYDROLOGIE – HEC HMS

DATE : JUIN 2013 MEC DOCUMENT 138fr

Ce document a été préparé par DAI pour évaluation par l'Agence des États-Unis
pour le développement international (USAID). 1
COMPÉTITIVITÉ
ÉCONOMIQUE DU MAROC

CONCEPTS DE BASE DE LA MODELISATION

HYDROLOGIQUE ET HYDRAULIQUE

APPORTS ET PRINCIPES D’UTILISATION DES OUTILS HEC, HEC-HMS

ET HEC-RAS, PLATEFORME D’INTEGRATION WMS ET OUTIL RIBASIM

FASCICULE 2 : MODELISATION EN HYDROLOGIE – HEC HMS

FASCICULE 2 : MODELISATION EN HYDROLOGIE – HEC HMS

Soumis à USAID/Maroc, Bureau de croissance économique ‐ Objectif d'aide n° 3 : Réduction des barrières au
commerce et à l'investissement

Par DAI

Numéro de contrat : EEM‐I‐00‐07‐00009‐00 : Ordre de mission : EEM‐I‐07‐07‐00009

L'opinion de l'auteur de ce document engage uniquement la responsabilité de ce

dernier et ne reflète pas nécessairement le point de vue de l'Agence des États‐
Unis pour le développement international (USAID) ni celui du Gouvernement des
États‐Unis
Programme Compétitivité Economique du Maroc

8, rue du Rif

Souissi

10 000 Rabat

Maroc

Tel: (212) 05 37 63 05 59

Fax: (212) 05 37 63 05 61

[email protected]

Http://www.mecprogram.ma
TABLE DES MATIERES
Préambule : .....................................................................................................................................................................................1

1. Introduction ................................................................................................................................................................................1
2. Modélisation des précipitations .............................................................................................................................................2
3. Modélisation des pertes...........................................................................................................................................................8
4. Modélisation du ruissellement direct ................................................................................................................................ 11
5. Modélisation des écoulements fluviaux avec HEC-HMS .............................................................................................. 15
6. Etapes d’une simulation avec HEC-HMS .......................................................................................................................... 19
7. Références bibliographiques ................................................................................................................................................ 19

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE i

TABLE DES ILLUSTRATIONS
Figure 1. : Polygones de Thiessen ................................................................................................................... 3
Figure 2. : Exemple de courbes isohyètes ..................................................................................................... 4
Figure 3. : illustration de la méthode ............................................................................................................. 5
Figure 4. : Facteurs de réduction des hauteurs déduits à partir des courbes HDF ............................ 6
Figure 5.: Exemple de la distribution de la Pluie basée sur la fréquence ............................................... 7
Figure 6. : Averse hypothétique SCS de type I ............................................................................................. 8
Figure 7. : Schéma conceptuel de l’algorithme du modèle SMA .......................................................... 11
Figure 8. : Modélisation d’un modèle de bassin pour le modèle d’onde cinématique ...................... 14

ii COMPETITIVITE ECONOMIQUE DU MAROC

PRÉAMBULE
MODÈLE HYDROLOGIQUE
HEC-HMS
1. INTRODUCTION
Le modèle hydrologique déterministe et conceptuel HEC-HMS (Hydrologic Modeling System) a pendant
longtemps été un modèle essentiellement appliqué à la simulation d'un événement donné. Les versions
récentes permettent de simuler en continu sur une longue période grâce à l'intégration de modules de
calculs des différents processus d'évapotranspiration, d’humidité du sol et d’accumulation/fonte de neige.

HEC-HMS est un modèle distribué qui permet de subdiviser un bassin versant en plusieurs sous-bassins,
qui sont considérés comme ayant chacune des caractéristiques homogènes. Il est particulièrement bien
adapté pour simuler le comportement hydrologique de bassins versants non urbanisés. HEC-HMS
permet également de simuler et d’incorporer des réservoirs et des dérivations.
Le modèle hydrologique HEC-HMS comprend une interface graphique conviviale, des capacités pour la
manipulation des données, résultats et graphiques, un système spécifique de gestion et de stockage de
données (DSS View), ainsi que des possibilités d’affichage et d’impression de résultats, et des outils de
gestion.

Ce logiciel permet de calculer des hydrogrammes de crues pour plusieurs objectifs à savoir les études
de drainage urbain, la prévision des crues, la conception des réservoirs et la réduction des effets des
inondations.

Les modules modèles développés sous HEC-HMS permettent la prise en charge des quatre processus
hydrologiques essentiels : Modèles pour le calcul des précipitations, du volume de ruissellement, du
ruissellement direct et des écoulements souterrains.

Le modèle hydrologique HEC-HMS a été développé par le Hydrologic Engineering Center (HEC) de
l’U.S. Army Corps of Engineers (USACE).

Une description détaillée de HEC-HMS et sa documentation complète peuvent être consultées et

téléchargés gratuitement sur le site Internet de l’USACE pour Windows, à
http://www.hec.usace.army.mil/software/hec-hms/download.html.

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE 1

2. MODÉLISATION DES PRÉCIPITATIONS

2.1. Définition. Méthodes basées sur les coefficients de pondération

2.2. Précipitations moyennes

A partir des mesures pluviométriques effectuées sur le bassin versant, on obtient la moyenne surfacique
des précipitations (Pmoy). Elle est calculée à partir de la moyenne arithmétique après avoir attribué un
coefficient de pondération pour chaque station pluviométrique:

 
  w  p t 
i 1
Pmoy  i T

w i
i

Où, wi : Coefficient de pondération affecté au pluviomètre i ;

Pi (t) : hauteur des précipitations mesurée au pluviomètre i au temps t.

Le poids affecté à un pluviomètre est fonction de la surface qu’il couvre. Les méthodes utilisées pour
calculer le coefficient de pondération sont :

Moyenne arithmétique

Cette méthode permet de calculer un coefficient égal à l’inverse du nombre de stations. C’est le même
coefficient qui sera appliqué à l’ensemble des stations de mesure se trouvant sur le bassin versant ou à
ses alentours.

1
wi 
N
Avec, N : Nombre de stations pluviométriques.

Polygones de Thiessen

Cette méthode permet de déterminer le coefficient de pondération en utilisant les polygones de

Thiessen (figure IV.1). Le coefficient de pondération est calculé par la formule ci- dessous :

aij
wi 
Aj

avec, aij : Surface d’intersection du " polygone j " et du " sous bassin i "
Ai : Surface totale du sous bassin i.

2 COMPETITIVITE ECONOMIQUE DU MAROC

Figure 1. : Polygones de Thiessen

Les isohyètes

Une isohyète est une ligne dont la pluviosité est la même sur tous ses points. Dans ce cas les coefficients
de pondération peuvent être calculés de la manière suivante :

bij
wj 
Ai

Avec, b ij : Surface délimitée par le sous-bassin i et les deux isohyètes j et j+1 ;

A i : Surface totale du sous-bassin i.

La pluie moyenne sur le sous-bassin i est calculée par la suite en utilisant la formule suivante :
K
Pmoy   w j Pj
j 1

Avec

h j  h j 1
Pj 
2

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE 3

 Avec, P j : Précipitation moyenne entre deux isohyètes j et j+1 ;
K : nombre total d'isohyètes.

Figure 2. : Exemple de courbes isohyètes

2.3. Répartition temporelle

Après avoir calculé la moyenne totale des précipitations, il faut ensuite trouver la répartition temporelle
des précipitations à partir de la répartition typique de l'évènement pluvieux notée Ptyp(t):

 
 Ptyp t  
Pmoy  P
  Ptyp t   moy
 t 

 w p t 
j j

où : Ptyp t   j

wj
j

Avec wj : coefficients de pondération affectés au pluviomètre j.

2.4. Inverse du carré de la distance

C’est une autre alternative à la méthode précédente qui comprend en fait deux étapes (calcul de Pmoy
puis la répartition temporelle). Cette méthode permet la réalisation directe de l’hyétogramme en faisant
intervenir la notion de nœuds entre pluviomètres. Les nœuds sont positionnés sur le bassin versant de
sorte à avoir une répartition spatiale adéquate de la précipitation. HEC-HMS définit alors la répartition
temporelle des pluies en un nœud par pondération des données pluviométriques donnée par la relation
suivante :

4 COMPETITIVITE ECONOMIQUE DU MAROC

 k
Pnode t    wi Pi
i

Où : wi : le coefficient de pondération affecté à la station i ;

K : nombre de stations.

Pour calculer les coefficients de pondération, HEC-HMS construit un repère orthogonal basé sur un
nœud (figure IV.3), puis il détermine pour chaque quadrangle délimité par deux axes une station (la plus
proche du nœud). Le facteur de pondération est calculé par la relation suivante :
1
d 2j
wj 
 1 2
dj
Où, 1 / dj2 : Inverse du carré de la distance de la station j au « nœud P ».

Figure 3. : Illustration de la méthode

Une fois la série Pnode (t) est établie pour tous les nœuds, la moyenne temporelle des précipitations est
calculée par la relation suivante :

P loy (t ) 
 w t P t 
node node

w node

Où : wnode est le poids affecté pour chaque nœud. Si un seul nœud est utilisé dans le bassin versant, wnode
sera égal à 1.00.

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE 5

 2.5. Evénements pluvieux hypothétiques
On peut créer avec le logiciel HEC-HMS des évènements pluvieux hypothétiques, c'est-à-dire réaliser
une simulation à partir de données qui ne sont pas issues de relevés réels. Cela peut s'avérer utile en
particulier lorsqu'on veut connaître les volumes d'écoulements pour dimensionner des ouvrages d’art.
On peut créer trois types d’évènements pluvieux standards :

 Pluie hypothétique basée sur la fréquence ;

 Pluie de projet standard : cette méthode fait intervenir des paramètres définis uniquement
pour le territoire américain par certains organismes ;
 Pluie hypothétique dont la distribution est définie par l'utilisateur.

a. Pluie hypothétique basée sur la fréquence

L’objectif de cette méthode est de définir un événement pour lequel la hauteur et la durée des
précipitations sont déterminées pour une probabilité de dépassement donnée. Pour définir cet
évènement avec HEC-HMS l’utilisateur doit spécifier la hauteur de précipitation en chaque pas de temps
pour une probabilité de dépassement choisi. Les hauteurs des précipitations peuvent être déduites à
partir des courbes HDF (Hauteur-Durée-Fréquence). HEC-HMS applique un coefficient de correction de
surface aux hauteurs données par les courbes HDF. Généralement, la distribution de l’intensité des
précipitations est non uniforme sur le bassin versant, la hauteur moyenne sur la surface du bassin est
inférieure à la hauteur donnée par les courbes HDF. Pour remédier à ce problème, le U.S. Weather
Bureau, a définit des facteurs qui permettent d’ajuster les hauteurs données par les courbes HDF à la
hauteur moyenne. Ces facteurs, exprimés en pourcentage, sont fonction de la surface du bassin et de la
durée de la pluie (figure IV.4).

Figure 4. : Facteurs de réduction des hauteurs déduits à partir des courbes HDF

6 COMPETITIVITE ECONOMIQUE DU MAROC

 HEC-HMS fait une interpolation linéaire pour calculer les hauteurs d’eau en pas de temps régulier.
Ensuite, il détermine un histogramme à partir des valeurs incrémentales des hauteurs de précipitations
par la méthode dite en « Bloc ». Cette méthode positionne la hauteur maximale en milieu de
l’hyétogramme. Les blocs restants sont arrangés alors dans l'ordre décroissant, alternativement avant et
après le bloc central. Cette méthode est illustrée sur la figure 5.

Figure 5.: Exemple de la distribution de la Pluie basée sur la fréquence

b. Pluie hypothétique dont la distribution est définie par l’utilisateur

Cette option donne à l’utilisateur la possibilité de définir la hauteur et la distribution de la pluie

hypothétique. On peut aussi donner la quantité totale d'eau tombée sur le bassin versant et la fraction
de cette quantité à chaque instant. Parmi les pluies hypothétiques les plus utilisées dans les simulations,
figurent les pluies SCS. La figure IV.6 illustre le cas d’une pluie SCS de type I.

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE 7

Figure 6. : Averse hypothétique SCS de type I

3. MODÉLISATION DES PERTES

HEC-HMS calcule les volumes d'écoulements en soustrayant aux précipitations les quantités d'eau qui
sont stockées, infiltrées ou évaporées sur le bassin versant. L’interception, l’infiltration, le stockage et
l’évaporation sont représentés par HEC comme étant des « pertes ».

3.1. Concepts de base

Les surfaces d'un bassin versant sont classées en deux catégories :

1. Surfaces directement connectées et imperméables, où l'écoulement est direct et se fait sans

pertes. Dans ce cas on utilise le modèle « sans pertes » ;

2. Surfaces perméables soumises à des pertes décrites par les différents modèles suivants :
 Modèle de perte initiale et à taux constant ;
 Modèle à déficit et à taux de perte constant ;
 Modèle basé sur le Curve Number (CN) ;
 Modèle de Green et Ampt.

Pour chacun de ces modèles, les pertes sont évaluées pour chaque intervalle de temps et soustraites à la
moyenne surfacique des précipitations de cet intervalle. La quantité d’eau restante désigne l'excès de

8 COMPETITIVITE ECONOMIQUE DU MAROC

précipitation. Cette quantité est considérée uniforme sur tout le bassin versant et représente le volume
d'écoulement de surface.

3.2. Le modèle de pertes initiales et à taux constant

Ce modèle considère que le potentiel du taux de pertes maximum, noté fc, est constant, et inclut le taux
de pertes initiales Ia qui représente l'interception et le stockage dans les dépressions de surfaces.
L’interception est une conséquence de l’absorption de la pluie par le couvert végétal et le stockage de
surface est la conséquence de la topographie du bassin versant : l’eau stockée dans les dépressions de
surface sera soit évaporée soit infiltrée. Tant que Ia n'est pas atteint, il n'y aura pas de ruissellement.

Ceci est traduit par le schéma suivant :

Si  P  I  P
1 et 0

Si  P  I et P  f
i t c  Pet  Pt  f c

Si  P  I et P  f
i t c  Pet  0
Où : Pt est la moyenne surfacique des précipitations au temps t
Pet le ruissellement au temps t donné par :

N.B : en général on ne déterminera pas directement les pertes initiales et le taux de pertes mais on
procèdera plutôt à un calage du modèle à partir de données réelles.

3.3. Le modèle de Curve Number (CN)

Ce modèle estime l'excès de précipitations comme une fonction des précipitations cumulées, de la
couverture et de l'humidité initiale du sol à partir de l’équation suivante :

Pe 
 p  I 2
p  I  S

Où, Pe : L’excès de précipitation ;

P : est le total des précipitations accumulées au temps t ;
Ia : Les pertes initiales ;
S : est le potentiel maximum de rétention.

On a par ailleurs la relation empirique ci-dessous, qui a été développée à partir d’une analyse menée sur
un ensemble de petits bassins versants expérimentaux qui donne Ia en fonction de S:

On obtient donc :

Pe 
 p  0.2S 2
p  0.8S

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE 9

Le potentiel de rétention maximum, S, et les caractéristiques du bassin versant sont reliées par
l’intermédiaire du numéro de courbe CN (Curve Number) par :

25400  254CN
S
CN

CN peut être estimé comme une fonction du type du sol, d’occupation du sol et des conditions
hydriques précédentes du bassin.

3.4. La méthode CN en grille

Une autre alternative de la méthode de CN qui peut être utilisée par le HEC-HMS est celle de CN en
grille, avec cette option, chaque cellule doit inclure : sa localisation, sa distance de l’exutoire, ses
dimensions et son CN (qui est donné en fonction du type de sol). Ainsi, HEC calcule l’excès des
précipitations pour chaque cellule indépendamment en utilisant l’équation (IV.14), ensuite, il calcule
l’excès total au niveau de l’exutoire en utilisant la méthode de ModClark (présentée ultérieurement). Ce
qui veut dire que la méthode de CN en grille n’est compatible qu’avec cette méthode de calcul.

3.5. Le modèle de Green & Ampt

Il est l’un des plus utilisé, c’est un modèle conceptuel d’infiltration et de précipitation. Sa fonction de
production est basée sur la loi de Darcy généralisée aux milieux non saturés combinée avec le principe
de conservation de la masse.

3.6. Le modèle continu SMA (Sol-Moisture Accounting)

(Prise en compte de l’humidité du sol)

Le modèle SMA représente le bassin versant avec une série de couches de stockage interconnectées
entre elles. La figure IV.7 représente la façon dont les couches de stockages sont connectées. Chacune
est représentée par le taux d’apport, de sortie et de sa capacité de stockage qui commande la quantité
d’eau perdue ou supplémentaire.

10 COMPETITIVITE ECONOMIQUE DU MAROC

Figure 7. : Schéma conceptuel de l’algorithme du modèle SMA

Le modèle SMA calcule l’écoulement à l’entrée et à la sortie de chaque couche de stockage.

3.7. Applications et limitations des modèles de ruissellement avec HEC-HMS

Le choix du modèle et l’estimation des différents paramètres sont une étape critique dans le
développement des données avec le logiciel HEC-HMS. On ne peut pas connecter n’importe quel
modèle de pertes avec n’importe quel modèle de transfert. Par exemple, la méthode de pertes en grille
n’est utilisée qu’avec le modèle de transfert de « ModClark » (méthode de Clark modifiée). Cependant,
ceux-ci sont seulement des directives et devraient être complétées par la connaissance et l'expérience
de l’utilisateur.

4. MODÉLISATION DU RUISSELLEMENT DIRECT

Dans cette partie, on présente les modèles qui simulent le processus du ruissellement direct des
précipitations excédentaires sur un bassin versant. HEC-HMS se rapporte à ce processus comme
transformation de l'excès de précipitation pour chacun des sous bassins en un écoulement à son
exutoire.

4.1. Les modèles liés à la méthode de l’hydrogamme unitaire (HU)

Ces modèles donnent une relation empirique entre l'excès de précipitations et le ruissellement direct. La
méthode de l'hydrogramme unitaire vise à déterminer l'hydrogramme de ruissellement direct à l'exutoire
d'un bassin versant à partir des hyétogrammes de l'averse correspondante reçue par ce même bassin.

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE 11

L'hydrogramme unitaire donne le débit de ruissellement par unité de hauteur d'eau tombée sur le bassin
versant. Cette méthode repose donc principalement sur l'hypothèse de linéarité entre l'excès de
précipitations et le ruissellement. On a donc l'équation de convolution suivante :
nm
Qn  P U
m 1
m n  m 1

où, Qn est le débit de ruissellement au temps n. ∆ t,

Pm est l’excès de précipitation entre m ∆ t et (m+1) ∆ t (en hauteur d’eau),
Un-m+1 est la valeur de l’hydrogramme unité au temps (n-m+1) ∆ t .

4.2. Modèles d’hydrogrammes unitaires synthétiques

4.2.1. Le modèle de SNYDER

Ce modèle permet d’estimer l’hydrogramme unitaire à partir des caractéristiques du bassin versant. Il est
basé sur trois paramètres : l'amplitude du pic, la base du temps total et le décalage entre le maximum de
l’hyétogramme de la pluie excédentaire et le pic de l'hydrogramme. Ce modèle est fondé par ailleurs sur
un événement standard tel que le décalage temporel (t p) qui est déterminé par la relation suivante :
t p  5 .5 t r
où : tr est la durée de l’évènement.

Si le décalage temporel est remarquablement différent du temps calculé par la relation précédente, on
utilise la relation suivante :
tr  t R
t pR  t p 
4
où, tR : Le décalage temporel de l’hydrographe désiré (introduit par l’utilisateur) ;
tpR : Le décalage temporel de l’hydrogramme unitaire (introduit par l’utilisateur).

Pour ce type d’événement, Snyder a établi une relation entre le décalage temporel et le pic de
l'hydrogramme unitaire.
Up Cp
C
Ad t
où, Up : Le pic de l’hydrogramme unitaire,
Ad : La surface drainée,
Cp : Coefficient du pic de l’hydrogramme unitaire,
C : Constante de conversion = 2.08 pour le Système international d'unités (SI),
t : Le temps de décalage. Il est égal à t p pour des événements standard et à tpR pour d’autres
durées.

12 COMPETITIVITE ECONOMIQUE DU MAROC

 4.2.2. Le modèle du SCS (Soil Conservation Service)
Ce modèle repose sur l'hydrogramme unitaire normalisé (qui est la moyenne de nombreux
hydrogrammes unitaires calculés pour différents bassins versant). Cet hydrogramme normalisé
représente le débit d'écoulement, Ut, comme une fraction du débit maximal, Up, et, Tp, l'instant du pic. On
a par ailleurs les relations empiriques suivantes :

Ad
U p  2.08
Tp
t
Tp   0.6Tc
2
où, ∆t : La durée de l’excès de précipitation ; Tc : le temps de concentration du bassin. Ainsi, il suffit de
connaître le temps de concentration pour remonter à Tp et Up et d’obtenir ainsi l'hydrogramme unitaire
désiré par simple multiplication de l'hydrographe unitaire normalisé.

4.2.3. Le modèle de CLARK (ModClark)

Il s'agit d'un modèle quasi-conceptuel puisqu'il prend en compte les deux processus intervenant dans la
transformation d'un excès de précipitation en ruissellement, à savoir le mouvement de l'eau de son
origine jusqu'à l'exutoire du bassin versant et l'atténuation de cette quantité d'eau par stockage lors de
son parcours. Le modèle de Clark modifié, fait intervenir un maillage du bassin versant qui permet de
mieux modéliser les temps de parcours.

Une fois que l’excès de précipitation est déterminé pour chaque grille du maillage du bassin versant, le
temps de transfert de l’excès de précipitation de la grille à l’exutoire est calculé par la relation suivante :

Ad
U p  2.08
Tp

t
Tp   0.6Tc
2

où, tcell : Temps de transfert de l’excès de précipitations de la grille ;

tc : Temps de concentration du bassin versant ;
dcell : Distance entre la grille et l’exutoire du bassin versant ;
dmax : Le plus long parcours du bassin versant.

4.2.4. Modèle de l'onde cinématique

Il s'agit d'un modèle conceptuel qui représente les versants du bassin comme des canaux à surface libre
dont le débit entrant correspond à l’excès de précipitations. Il résout donc les équations (de continuité
et de quantité du mouvement) qui régissent les écoulements non stationnaires en eau peu profonde
pour tracer l'hydrogramme de ruissellement.

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE 13

 Le bassin versant étant découpé en deux plans séparés par le cours d'eau dans lequel les écoulements se
déversent (figure IV.8). Le modèle de l'onde cinématique représente les écoulements de surface sur ces
plans.

Figure 8. : Modélisation d’un modèle de bassin pour le modèle d’onde cinématique

Sans rentrer dans la théorie complète, et sous quelques hypothèses simplificatrices, l’équation de
mouvement et celle de continuité donnent :
AT A
  m 1 T  q
t x

avec

S 0. 5
  1.486
N

5
m
3N
où , AT : La section transversale ;
S : est la pente,
N : est un coefficient de rugosité tabulé;
q : est le débit d’apport latéral par unité de longueur du canal.

HEC-HMS représente les écoulements de surface sur un canal rectangulaire de largeur unité. Cette
méthode est aussi utilisée pour les écoulements fluviaux. Dans les deux cas, la résolution des équations
se fait par la méthode des différences finies pour assurer la précision et la stabilité. Pour utiliser la

14 COMPETITIVITE ECONOMIQUE DU MAROC

 méthode de l’onde cinématique dans la modélisation du ruissellement, le bassin versant est décomposé
en divers éléments :
- Les plans de ruissellement de surface ;
- Les canaux de collection ;
- Le cours d'eau principal.

Les informations nécessaires à connaître pour chacun de ces éléments sont représentées sur le tableau
ci-dessous

Plan de ruissellement Canaux Cours d’eau

de surface de collection principal

Longueur caractéristique Aire drainée par le canal Pente du canal

Pente représentative Longueur caractéristique Forme de la section e dimensions
Coefficient de rugosité Pente du canal Coefficient de Manning
Aire représentée par le plan Forme de la section et dimensions
Paramètre de perte Coefficient de Manning

Tableau 1. : Informations à connaitre pour chaque élément du bassin

5. MODÉLISATION DES ÉCOULEMENTS FLUVIAUX AVEC HEC-

HMS
Les modèles d'écoulements fluviaux inclus dans le HEC-HMS permettent de calculer un hydrogramme
en aval du bassin versant, connaissant l'hydrogramme amont. Tous ces modèles utilisent les équations de
continuité et de quantité de mouvement. On trouvera plus de détails à la section 3.

1) Equation de continuité

Elle s’écrit sous la forme suivante :

V y y
A  VB  B  q
x x t
où, B : Largeur de la surface libre ;
q : Le débit unitaire, calculé par unité de longueur du canal ;
A : L’aire de la section du cours d’eau ;

2) Equation de quantité de mouvement.

L’équation du mouvement s’écrit comme suit :

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE 15

 y V V 1 V
S f  S0   
x g x g t
où, Sf : est le gradient hydraulique ;
S0 : est le gradient au fond du canal ;
V : vitesse d’écoulement ;
y : hauteur hydraulique ;
x : la longueur du chemin de l’écoulement ;

Les principales hypothèses retenues lors de l’élaboration des ces équations sont :

 La vitesse est constante sur la surface et elle reste horizontale sur chaque section du cours
d’eau ;

 La pression reste hydrostatique en chaque point du cours d’eau ;

 Aucune circulation latérale et secondaire ne se produit ;

 Le lit du canal est fixe ; l’érosion et le dépôt ne changent pas la forme de la section du canal ;

 La résistance à l'écoulement peut être décrite par des formules empiriques, telles que l'équation
de Manning ou de Chézy.

5.1. MODÈLE DE PULS (MODIFIED PULS)

Ce modèle est basé sur une approximation en différences finies de l'équation de continuité, et sur une
représentation empirique de l'équation de quantité de mouvement. On considère qu'il n'y a pas
d'apports latéraux. L'équation de continuité devient alors :

On utilise alors la méthode des différences finies en écrivant que :

où It et Ot sont les débits moyens respectivement amont et aval pendant une période ∆ t,

∆ St, étant la variation du stock d'eau durant le même pas de temps. On obtient alors l'équation suivante
:

 S t Ot   I t 1  I t   S t 1 Ot 1 
     
 t 2   2   t 2 

Ainsi, au temps t, on connaît les termes de droite, et l'on cherche à déterminer les termes de gauche.
On a donc deux inconnues St et Ot, et l'on doit donc connaître une relation supplémentaire entre ces
deux variables. Une fois cette relation déterminée, on résout le nouveau système d'équations par
récurrence. En ce qui concerne la relation liant St à Ot, on peut la déterminer par des observations

16 COMPETITIVITE ECONOMIQUE DU MAROC

historiques des valeurs des débits et du volume d'eau stocké au cours du temps.

5.2. Modèle de Muskingum

Ce modèle utilise comme le modèle précédent une approximation simple en différences finies de
l'équation de continuité :

 I t 1  I t   Ot 1  Ot   St  St 1 
     
 2   2   t 
On écrit ensuite que le volume d'eau stocké est la somme d'un volume stocké constant et d'une
variation du stock :
S t  T p Ot  T p X I t  Ot   T p  XI t  1  X Ot 
avec, T p : Temps de parcours et X paramètre permettant une pondération (0<X<0.5). Ainsi, si le
stockage d'eau dans la rivière est contrôlé par les conditions aval, on met X = 0 (On retrouve alors le
modèle de réservoir linéaire présenté plus haut). Au contraire, on prendra X = 0.5, pour donner un
poids similaire aux débits entrants et sortants.

5.3. Modèle du décalage (lag model)

Ce modèle est le plus simple de tous et est extrêmement utilisé surtout pour modéliser le drainage
urbain. On considère que l'hydrogramme aval est simplement l'hydrogramme amont décalé dans le temps
d'une certaine durée (lag) : La valeur de ce décalage peut être estimée comme la durée entre les deux
pics des deux hydrographes (ou leurs centres de gravité respectifs).

5.4. Modèle de l’onde cinématique

Ce modèle est basé sur une approche différences finies de l’équation de continuité et une approche
simplifiée de l’équation du mouvement. Cette approche est décrite en détail dans la section précédente.
Les informations nécessaires à ce modèle sont :

 Le type de la section du cours d’eau : trapézoïdale, rectangulaire ou circulaire ;

 Les dimensions principales : la base de la section du canal, le diamètre de la conduite ;
 Pente des berges ;
 La longueur du bief ;
 La pente de ligne d’énergie ;
 Coefficient de la rugosité de Manning (n).

5.5. Modèle de Muskingum-Cunge

Le modèle Muskingum est facile à utiliser mais emploie des paramètres qui n'ont pas de signification
physique et qui sont donc difficiles à évaluer. De plus ce modèle est basé sur des approximations qui
s'avèrent souvent fausses. Le modèle de Muskingum-Cunge évite ces erreurs. Il est basé sur l'équation de
continuité incluant un débit latéral et sur la forme diffusive de l'équation de quantité de mouvement :

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE 17

 A Q
 q
t X

Y
S f  S0 
X

En combinant ces deux équations, on obtient alors l'équation de diffusion–convection de Miller et Cunge,
1975.

Q Q  2Q
c   c q1
t X X 2

Avec c la célérité la diffusivité hydraulique définies par :

Q
c
X

Q

2 B S0

L'approximation en différences finies des dérivées partielles permet d'obtenir l'équation suivante:

Qt  C1 I t 1  C2 I t  C3 Ot 1  C4 q1 x

t t
 2X 2 1  X  
C1  K C3  K
t t
 2 1  X   2 1  X 
K K
t t
 2X 2
C2  K C4  K
t t
 2 1  X   2 1  X 
K K

x
K
c
1 Q 
X  1  
2  B S0 c x 

Les coefficients, Q, B et les quatre paramètres Ci sont des variables. Ils sont recalculés à chaque pas de
temps et d'espace par HEC-HMS.

18 COMPETITIVITE ECONOMIQUE DU MAROC

6. ETAPES D’UNE SIMULATION AVEC HEC-HMS
Pour effectuer une première simulation, il est nécessaire de créer une base de données pour
caractériser le bassin versant et déterminer les averses qui seront employées comme données
d’entrée du modèle. Les principales étapes de la simulation sont :
 Démarrage d’un nouveau projet
 Définition du modèle de bassin versant
 Entrée des données météorologiques
 Définition du modèle de précipitation
 Définition des paramètres de contrôle de la simulation
 Création et exécution d’une simulation
 Visualisation des résultats

7. RÉFÉRENCES BIBLIOGRAPHIQUES
Aspects Hydrologiques, hydrauliques, logiciels

Vieux, B.E., 2004. Distributed Hydrologic Modeling using GIS. Kluwer Academic Publishers, Dordrecht,
The Netherlands.

Davie, T., 2008. Fundamentals of Hydrology, Second edition, Routledge, Taylor & Francis Group, 200p.

Ambroise, B., 1991. Hydrologie des petits bassins versants ruraux en milieu tempéré, processus et
modèles. Séminaire du Conseil Scientifique du département « Science du Sol » de l’Institut National
de Recherches Agronomiques de Dijon, 26-27 mars 1991.

Ambroise, B., 1999. La dynamique du cycle de l’eau dans un bassin versant. Processus, Facteurs,
Modèles. Editions H.G.A., Bucarest.

Hydrologic Engineering Center. 2003. Geospatial Hydrologic Modeling Extension HEC- GeoHMS
User’s Manual Version 1.1. Davis, CA: U.S. Army Corps of Engineers.

Hydrologic Engineering Center. 2006. GageInterp A Program for Creating a Sequence of

HEC-DSS Grids from Time-Series Measurements Version 1.4. Davis, CA: U.S. Army Corps of
Engineers.

Hydrologic Engineering Center. 2006. HEC Data Storage System Visual Utility Engine (DSSVue) User’s
Manual Version 1.2. Davis, CA: U.S. Army Corps of Engineers.

Hydrologic Engineering Center. 2008. Hydrologic Modeling System HEC-HMS User’s Manual Version
3.2. Davis, CA: U.S. Army Corps of Engineers.

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE 19

Panigrahi, B.K., Harris, D.J., and C.S. Pathak, 2012. Hydrology and Hydraulics: Watershed Management
and Flood Assessment, Short Course offered during 5th International Perspective on Water
Resources & the Environment January 4-7, 2012, Marrakech, Morocco

Beven, K.J., 2002. Rainfall-Runoff Modelling The Primer, John Wiley & Sons Ltd, 359 p.

Singh V.P., 1995. Watershed modeling. In: Singh V.P., ed. Computer models of watershed
Davie, T. (2008) : FUNDAMENTALS OF HYDROLOGY, Second edition

Routledge Fundamentals of Physical Geography, Routledge Taylor and francis Group

Vieux, B.E., 2004. Distributed Hydrologic Modeling using GIS. Kluwer Academic Publishers, Dordrecht,
The Netherlands.

Singh, V.P. and Frevert, D.K., editors, 2002a. Mathematical Models of Large Watershed Hydrology.
Water Resources Publications, Highlands Ranch, CO.

William A. Scharffenberg and Matthew J. Fleming (2010) : Hydrologic Modeling System HEC-HMS,
User's Manual, Version 3.5 August 2010, US Army Corps of Engineers, Hydrologic Engineering Center

Ancey, C. (2012) : Notes de cours Hydraulique à surface libre, version 4.2 du 2 mai 2012,
Laboratoire hydraulique environnementale (LHE), École Polytechnique Fédérale de Lausanne

Chow, V. T., 1959, Open-Channel Hydraulics , McGraw-Hill Book Company, New York, N. Y.

French, R. H., 1985, Open-Channel Hydraulics , McGraw-Hill Book Co., New York, NY.

Henderson, F.M., 1966, Open-Channel Flow , Macmillan Co., London, UK.

Abbott, M. B., 1979, Computational Hydraulics: Elements of the Theory of Free Surface Flows, Pitman,
London.

Abbott, M. B., and Basco, D. R., 1990, Computational Fluid Dynamics, John Wiley, New York, NY.

Basco, D. R., 1983, Computation of Rapidly Varied, Unsteady Free-Surface Flow, U. S. Geological Survey
Report WRI 83-4284.

Bell, S. W., Elliot, R. C., and Chaudhry, M. H., 1992, “Experimental results of two-dimensional dam-break
flows,” Jour. Hyd. Research, vol. 30, no. 2, pp.225-252.

Boussinesq, J., 1877, “Essais sur la theorie des eaux courantes,” Memoires pre- sertes par divers
Savants a l’Academie des Sciences de l’Institut de France, vol. 23, pp. 1-680; vol. 24, 1-64.

20 COMPETITIVITE ECONOMIQUE DU MAROC

Chow, V. T., 1959, Open-Channel Hydraulics, McGraw-Hill, New York, NY. Cunge, J. A., Holly, F. M., and
Verwey, A., 1980, Practical Aspects of Computational River Hydraulics, Pitman, London.

DHI. MIKE 21 Flow Model: Hydrodynamic Module Scientific Documentation. MIKE by DHI, 2009.

Guide de laboratoire HEC-RAS (2009): Exemple pour la rivière du Loup, HEC-RAS Version 3.1.3,
CTN-762 Ressources hydriques, Ecole supérieure de Technologie, Montréal,

Gary W. Brunner (2010) ; HEC-RAS River Analysis System, Hydraulic Reference Manual, Version 4.1,
January 2010 Approved for Public Release. US Army Corps of Engineers, Hydrologic Engineering
Center

Ouazar, D. (2006) : Ecoulements en régime permanent et transitoire, Notes de cours, Ecole Mohammadia
d’Ingénieurs, Rabat.
Ouazar, D. (2010) : Modèles numériques pour l’ingénieur, Notes de cours, Ecole Mohammadia
d’Ingénieurs, Rabat
Ouazar, D. (2001) : Modélisation des problèmes d’hydrologie et d’hydraulique, Notes de cours, Ecole
Mohammadia d’Ingénieurs, Rabat
Manuel d’utilisation du logiciel WMS, Aquaveo

Projets de fin d’études :

Bennajah, A., et E. Messka, 2011. Modélisation de Prévision des crues par utilisation de la Plateforme
WMS, Cas d’étude Bassin de l’Ouergha, Mémoire de Projet de Fin d’Etudes, Ecole Mohammadia
d’Ingénieurs, Université Mohamed V Agdal, Rabat

AIT ZEOUAY, H. 2011. Modélisation de prévision des crues par utilisation de la plateforme ATHYS,
Cas d’étude Bassin de l’Ouergha, Mémoire de Projet de Fin d‟Etudes, Ecole Mohammadia
d’Ingénieurs, Université Mohamed V Agdal, Rabat
ESSADAOUI, M., 2012). Modélisation de Prévision des crues par utilisation de la Plateforme WMS,
Cas d’étude Bassin de l’Ouergha, Mémoire de Projet de Fin d‟Etudes, Ecole Mohammadia
d’Ingénieurs, Université Mohamed V Agdal, Rabat
Sites internet
http://www.aquaveo.com/wms-intro
http://www.xmswiki.com/xms/WMS:Quick_Tour_-_Layout
http://www.xpsoftware.com/products/watershed-modeling-system-wms/
http://srtm.csi.cgiar.org

CONCEPTS DE BASE DE LA MODELISATION HYDROLOGIQUE ET HYDRAULIQUE 21

22 COMPETITIVITE ECONOMIQUE DU MAROC