ENAC/ISTE/HYDRAM Logo optimisé par

HYDROTHEQUE : base de données d’exercices en Hydrologie

J.-D.Bonjour, SI-DGR
13.4.93

Cours : Hydrologie Appliquée / Thématique : Etude des Crues

ÉCOLE POLYTECHNIQUE
FÉDÉRALE DE LAUSANNE

Exercice n° HA 0808 - Corrigé

Prédétermination des débits de crue de fréquence rare par la méthode du
GRADEX – Application au bassin versant du Mouhoun à Boromo
(Burkina-Faso)

Données de l’exercice :
L’exercice porte sur le bassin versant du Mouhoun à Boromo dont les caractéristiques sont
présentées dans le tableau 1-énoncé. Les données concernant les précipitations maximales journalières
annuelles enregistrées de 1930 à 1975 se trouvent dans le tableau 2-énoncé. Les données nécessaires à
la réalisation de cet exercice sont aussi regroupées dans le fichier Excel « HA0808_enonce.xls » ou
dans le fichier Excel « HA0808_feuilledecalcul.xls » qui est à compléter. Les résultats sont
disponibles sur le fichier Excel « HA0808_corrige.xls ».

Question 1 : Calcul de la valeur du « GRADEX » des pluies journalières

observées
La loi de Gumbel est souvent utilisée pour ajuster les séries de pluies maximales et les débits
correspondants. Dans ce cas et dans ce cas uniquement, le caractère exponentiel de cette distribution
est décrit par la pente de la droite d’ajustement des pluies observées. La pente de cette droite est le
GRAdient de cette distribution Exponentielle, d’où le nom de la méthode GRADEX.
Chercher le « GRADEX » des pluies journalières observées consiste donc à procéder à
l’ajustement des valeurs maximales annuelles selon une distribution de Gumbel pour en déduire une
estimation du paramètre b1 de la droite d’ajustement x q = a + bu q (avec u : variable réduite de
Gumbel).

 Méthode à appliquer : ajustement statistique d’une série de données

Un modèle fréquentiel très souvent utilisé pour décrire le comportement statistique des valeurs
extrêmes est la distribution statistique de Gumbel (loi double exponentielle ou loi de Gumbel). La
fonction de répartition de la loi de Gumbel F(x) s’exprime de la manière suivante :
⎛ ⎛ x − a ⎞⎞ x−a
F ( x) = exp⎜⎜ − exp⎜ − ⎟ ⎟⎟ (1) avec la variable réduite suivante : u = (2)
⎝ ⎝ b ⎠⎠ b
où a et b sont les paramètres du modèle de Gumbel.
La distribution s’écrit alors de la manière suivante :
F ( x) = exp(− exp(−u ) ) (3) et u = − ln(− ln(F ( x) )) . (4)
L’avantage d’utiliser la variable réduite est que l’expression d’un quantile est alors linéaire
( x q = a + bu q ).

En conséquence, dès lors que les points de la série à ajuster peuvent être reportés dans un système
d’axes x − u , il est possible d’ajuster une droite qui passe le mieux par ces points et d’en déduire les
deux paramètres a et b de la loi. L’estimation des paramètres a et b de l’ajustement peut se faire
graphiquement (ajustement à l’œil ou à l’aide d’une régression statistique), ou selon une méthode
mathématique comme celle des moments.

Mise à jour le 14.04.2004 HA 0808 - Page 1

En pratique il s’agit essentiellement d’estimer la probabilité de non dépassement F(xi) qu’il convient
d’attribuer à chaque valeur xi. Il existe de nombreuses formules d’estimation de la fonction de
répartition Fˆ ( x ) à l’aide de la fréquence empirique. Elles reposent toutes sur un tri de la série par
valeurs croissantes permettant d’associer à chaque valeur son rang r. Des simulations ont montré que
pour la loi de Gumbel, il faut utiliser la fréquence empirique de Hazen :
r − 0.5
( )
F x[r ] =
n
(5)

où r est le rang dans la série de données classée par valeurs croissantes, n est la taille de l’échantillon,
x[r] la valeur de rang r.
Rappelons encore que le temps de retour T d'un événement est défini comme étant l'inverse de la
fréquence d'apparition de l'événement. Soit :
1
T=
1 − FQ ( xQ )
. (6)
A l’aide de l’ajustement, il est alors possible d’estimer la pluie maximale pour un temps de retour
donné.
 Démarche et résultats
Etape 1 : Préparation de la série de données des précipitations maximales journalières annuelles de
pointe.
ƒ Trier les valeurs dans l’ordre croissant.
ƒ Attribuer un rang à chaque valeur.
Etape 2 : Calcul de la fréquence empirique pour chaque rang (Hazen, équation (5)).
Etape 3 : Calcul de la variable réduite « u » du Gumbel (équation (4)).
Etape 4 : Représentation graphique des couples (ui, xi) de la série à ajuster (figure 1).
90

80

70
Lame précipitée [mm/24h]

60

50

40

30

20

10

0
-2 -1 0 1 2 3 4 5
Variable réduite de Gumbel u [-]

précipitations observées Droite d'ajustement des pluies

Figure 1. Ajustement graphique du modèle (calcul des paramètres « a » et « b » de la droite

d’ajustement de Gumbel par la méthode graphique)

Mise à jour le 14.04.2004 HA 0808 - Page 2

 Etape 5 : Ajustement d’une relation linaire de type x q = a + bu q aux couples (ui, xi) (figure 1).
Avec un ajustement de type graphique (à l’œil), on a alors une estimation des paramètres a1 et b1
(a1 = 40.5 et b1 = 9.9 ). La pente b1 de cette droite n’est autre que le GRAdient Exponentielle des
pluies.

Question 2 : Estimation des débits de pointe de temps de retour 20, 50, 100 et
500 ans par la méthode du GRADEX

 Méthode à appliquer : Méthode du GRADEX

La méthode du GRADEX a été conçue, mise au point, testée et validée par le groupe de recherche
d’EDF 1 (Guillot and Duband, 1967). Elle a pour but de rechercher les débits maximaux de crues pour
des fréquences d’apparition rares à très rares ( temps de retour de plus de 100 ans). Elle s’applique
notamment lorsque l’on dispose d’une longue série de pluie et d’une courte série de débit (env.10 ans)
sur le bassin, et spécialement dans les régions de montagne.
Le principe de la méthode est de considérer que la droite de distribution des débits est parallèle à
partir du seuil, correspondant au temps de retour 10 ans, à celle des pluies (figure ci-contre).
Cette méthode s’appuie sur différentes hypothèses :
1. Les débits maximums recherchés sont provoqués uniquement par des pluies maximales,
uniformément réparties sur le bassin.
2. Les pluies maximales et les débits correspondants (débits maximaux) suivent une même loi de
distribution statistique, dite des "extrêmes" en raison de la nature du phénomène recherché (crues
rares).
3. Ceci exprime surtout le fait qu’à partir d’une certaine valeur de pluie, le comportement
asymptotique du débit sera identique à celui des pluies. Selon les auteurs de cette méthode, ce seuil
représente le taux de saturation du bassin qui est atteint après un événement pluviométrique qui
provoque un débit décennal (T = 10 ans).

L’application de la méthode du GRADEX

)
4h
4h

implique plusieurs contraintes :

p l dr es

r2
s) e
es la d
2

uie oit

Pluies
su
t d à ion
urs

Débits [mm/24h.]
en llèle olat

1. La durée des pluies considérées doit

es
lui

t em ara ap

strictement correspondre à celle des débits

p

jus (p xt r
es

tg(α )
d'a bit s e d'e

(même Δt et en général 24 h). Elle est

td

α
en

dé roit
t em

conditionnée par le temps de concentration des

D
jus

Q10
eaux du bassin au point d’intérêt.
d'a
e
oit
Dr

2. Les unités des pluies et des débits doivent être

débits mesurés identiques, si l’on procède à l’application de
cette méthode en utilisant la loi de Gumbel (en
0 mm/24h).
0 Variable réduite de Gumbel u = -ln(-ln(F)) 2.25
3. Les limites d’application de cette méthode sont
Fréquence F 0.9 conditionnées par des temps de concentration
tc variant de 1 heure à 4 jours. La méthode ne
Temps de retour T = 1(1-F) [an]
10 peut donc s’appliquer qu’à des bassins versants
de 5000 km2 au maximum.
Il convient encore de signaler que les résultats obtenus par extrapolation sont des débits maximaux
moyens qui résultent de pluies maximales moyennes. Il s’agit donc de multiplier ces valeurs de débits

1
Electricité De France

Mise à jour le 14.04.2004 HA 0808 - Page 3

par le coefficient de pointe moyen (coefficient de pointe=débit de pointe/débit moyen) pour obtenir le
débit maximum instantané.

 Résultats :
L’ajustement des précipitations selon la loi de Gumbel étant fait, on peut alors tracer la fonction de
répartition de la loi de distribution des débits moyens journaliers maximaux, comme la parallèle à la
droite des pluies (b2 = b1 = 9.9) et passant par le point pivot correspondant au débit moyen journalier
de temps de retour 10 ans. Il ne faut pas oublier de convertir les valeurs de débit journalier (en m3/s)
en lame écoulée (en mm) afin de pouvoir appliquer la méthode du GRADEX (ici,
Qp(10) = 150 m3/s = 43.2 mm/24h).
En appliquant la loi d’ajustement extrapolée des débits pour les temps de retour 20, 50, 100 et 500
et en faisant la conversion adéquate pour avoir des valeurs en m3/s, on obtient des valeurs de débits
moyens journaliers. Il faut encore calculer les débits de pointe en multipliant les valeurs de débits
trouvées précédemment par le coefficient de pointe.
Les résultats sont présentés dans la figure 1 et le tableau 1 ci-dessous :
90

80
Lame précipitée / écoulée [mm/24h]

70

60

50

40

30

20

10

0
-2 -1 0 1 2 3 4 5 6 7
Variable réduite de Gumbel u [-]

précipitations observées Droite d'ajustement des pluies Droite d'ajustement des débits

Figure 1. Ajustement graphique pour les pluies et droite d’extrapolation pour les débits.

Tableau 1. Débits moyens maximums journaliers et débit de pointe

de temps de retour 10, 20, 50, 100 et 500 ans

T F(x) u QT QT Qp
[ans] [-] [-] [mm/24h] [m3/s] [m3/s]

10 0.900 2.25 43.2 150 240

20 0.950 2.97 50.3 175 280
50 0.980 3.90 59.6 207 331
100 0.990 4.60 66.5 231 370
500 0.998 6.21 82.6 287 459

Mise à jour le 14.04.2004 HA 0808 - Page 4