11728, 1226 PM Hydrologie générale : Analyse tréquentille Corrigé de l'exercice "Etimation des débits de crue pour différents temps de retour par la méthode statistique — Application au bassin versant de la Mentue a Yvonand (VD, Suisse)" object de cet exercice est d'estimer les débits de pointes (déblts maximaux) correspondants & un certain temps de retour, cesta-dlre a une certaine probabilié d apparition donnée. Le corigé se trouve aussi dans la fchior Excel suivant 1 Méthode & appliquer : ajustement statistique d'une sére de données -Gumbel analyse fréquentielle dune longue série de débits maximaux permet destimerle temps de retour dune valeur particuliére, Cette prédiction repose surla definition et la mise en ceuvre d'un modele fréquentiel qui est une équation décrvant (modélisant) le comportement statistique dun processus. Ces modéles décrivent la probabilté d'appariion d'un événement de valeur donnée. C’est du choix du modele fréquentiel (et plus partculigrement de son type) que dépendra la validté des résultats de analyse fréquentelle, Un modéle fréquentialtrés souvent uilisé 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éparttion de la li de Gumbel F(x) s'exprime de la maniére suivante x-a FQ@)- oxf- {4 (1) avec a variable réuite suivante: w = 5% (ay 01a et b sont les paramétres du moddle de Gumbel. La distnbution s'écit alors de la maniére suivante : (x) = exp (~ exp(—te)) (9) et x = - In (- In F(x)}) (4). ‘avantage dutliser la variable réduite est que Vexoression dun quantile est alors linéaire : x, ~ a+ Buty En conséquence, és lors que les points de la série & ajuster peuvent are reporiés dans un systéme d'axes x- u il est possible 'ajuster une Groite qui passe le mieux par ces points t d’en déduire les deux paramatres a ot b dela lo. Il existe différentes méthodes d ajustement ‘méthode graphique (ajustement a oxi ou a aide d une regression statistque), méthode des moments ect En pratique, i sagit essentillement 'estimer la probabilté de non dépassement F(xi) quil convient qattribuer & chaque valeur xi Il existe de ombreuses formules c'estimation de la fonction de répartition & Taide de la fréquence empirique, Elles reposent toutes sur un tide la série par valeurs croissantes permettant d associer & chaque valour son rang r. Des simulations ont montré que pour la loi de Gumbel, faut uiliser Ta fréquence empirique de Hazen 05 F © ca st ng done sro do domes classé par valu cssriog, noe a tile de chan, x] vow de rng Rappoons encore quo le tops de retour T dun éénerrent est fi comm stant ware de a frequen capparion de Féénemert Sit 1 (6) =F) 2 Démarche et résultats Etape 1 : Préparation de la série de données des débits de point. * Tier las valeurs dans Fordre rolssant * Attribuer un rang & chaque valeur Etape 2 : Calcul de la fréquence empiriaue 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, x) de la série & luster (igure 1). hitpssocho2 op chie-rologiefchapitres/annexes/AnalFrequ_rep html us11728, 1226 PM Hydrologie générale : Analyse tréquentille — z ee § g po y= 7076 EEN i al g s & a oo 0 4 2 3 4 6 variable réduite de Gumbel Etape 5 : Ajustement d'une relation linalre de type aux couples (ul, x!) figure 1) et en déduire es deux paramatres a et b) ‘Avec un ajustement de type graphique (acl), on a alors une estimation des paramatres a et b :a voir aussi avec la méthode des moments, 5 5a 7.98 Etape 6 : Utlisation ou modéle statistque pour estimer des débits de pointe de différents temps de retour T. Par exemple pour 100 ans, on suit es étapes suivantes ‘Calcul de la fréquence de non-dépassement daprés la relation (6) : F =0.99 ‘Calcul de la variable réduite de Gumbel correspondante d apres la relation (4) u= 4.60: *Caleul du quantile correspondant d'aprés la relation linéaire (avec a et b fours par'étape 5 précédente ) 100 = 25.5+4.60 . 7.98 = 62.2 mils (On a de méme pour les autres temps de retour 05-37.5 m/s 220 =49.2 mis (050=56.6 ms ‘Autre méthode : Méthode des moments. La méthode des moments consiste a égaler les moments des échantilons avec les moments théoriques de la loi ‘moments les paramétres a et b sont calculés daprés les formules p-B5 “7 avec y = 0.5772 (constante Euler), a= g-by «r évarttype des valeurs composant 'échantilon. H:: moyenne de Féchantillon. Das ores posible estimor las bits dont la représontaton graphique est une droite dbquation Gratin avec: u variable ite (t équation (4) (on obtent (05=36.7 m/s (020 =47.3 m/s (050-54 mils (2100= 59 ms [Link]-rologiefchapitres/annexes/AnalFrequ_rep ml Parla méthode des11728, 1226 PM Hydrologie générale : Analyse tréquentille Remarque: Dans lexemple ci-dessus, las deux méthodes donnent des résultats tr8s praches lune de autre. La méthode des moments est rettement plus rapige & appliquer, elle présente cependant un désavantage par rapport la méthode graphique. L'ajustement graphique permet en effet de repérerd'éventuels points qui ne sont pas bien alignés et de ne pas en tenie compte, On pourrait également voir sla série Comportalt une « rupture » cest-d-die un changement de pente et donc un changement des paramétres de la lo statistique. De maniére générale, 'ajustement manuel donne souvent beaucoup dinformations surla série étudiée, [Link]-rologie/chapitres/annexes/AnalFrequ_rep html11728, 1226 PM Hydrologie générale : Analyse tréquentielle [Link]-drologiefchapitreslannexesiAnalFrequ_rep html 48