Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Importer

Bienvenue sur Scribd !

  • 267 vues

    HA0801 Corrige

    Transféré par

    jolegende
    Ce résumé présente les étapes clés de l'estimation du débit de crue de temps de retour 30 ans pour un bassin versant à l'aide de trois séries de données. La méthode statistique de Gumbel est appliquée à chaque série pour ajuster les paramètres du modèle aux données.

    Droits d'auteur :

    Attribution Non-Commercial (BY-NC)
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd

    HA0801 Corrige

    Transféré par

    jolegende
    267 vues4 pages
    Ce résumé présente les étapes clés de l'estimation du débit de crue de temps de retour 30 ans pour un bassin versant à l'aide de trois séries de données. La méthode statistique de Gumbel est appliquée à chaque série pour ajuster les paramètres du modèle aux données.

    Description originale:

    methodfe du gradex

    Copyright

    © Attribution Non-Commercial (BY-NC)

    Formats disponibles

    PDF, TXT ou lisez en ligne sur Scribd

    Avez-vous trouvé ce document utile ?

    Ce contenu est-il inapproprié ?

    Ce résumé présente les étapes clés de l'estimation du débit de crue de temps de retour 30 ans pour un bassin versant à l'aide de trois séries de données. La méthode statistique de Gumbel est appliquée à chaque série pour ajuster les paramètres du modèle aux données.

    Droits d'auteur :

    Attribution Non-Commercial (BY-NC)
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    267 vues4 pages

    HA0801 Corrige

    Transféré par

    jolegende
    Ce résumé présente les étapes clés de l'estimation du débit de crue de temps de retour 30 ans pour un bassin versant à l'aide de trois séries de données. La méthode statistique de Gumbel est appliquée à chaque série pour ajuster les paramètres du modèle aux données.

    Droits d'auteur :

    Attribution Non-Commercial (BY-NC)
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    Vous êtes sur la page 1sur 4

    ENAC/ISTE/HYDRAM HYDROTHEQUE : base de donnes dexercices en Hydrologie Cours : Hydrologie Applique / Thmatique : Etude des Crues

    Logo optimis par J.-D.Bonjour, SI-DGR 13.4.93

    COLE POLYTECHNIQUE FDRALE DE LAUSANNE

    Exercice n HA 0801 - Corrig


    Estimation du dbit de crue de temps de retour 30 ans par la mthode statistique partir de 3 sries (annuelle, tronque et gonfle) Application au bassin versant bassin versant du Parimbot (FR, Suisse)
    Donnes de lexercice :
    Lexercice porte sur le bassin versant du Parimbot (6.75 km2). Les donnes ncessaires la ralisation de cet exercice se trouvent dans le tableau 1 (il sagit des dbits instantans enregistrs sur plusieurs annes et plus grand que 0.35 m3/s). Ces donnes sont aussi regroupes dans le fichier Excel HA0801_enonce.xls ou dans le fichier HA0801_feuilledecalcul.xls complter. Les rsultats sont aussi disponibles sur le fichier Excel HA0801_corrige.xls .

    Question 1 : Construction des sries de maxima annuels, gonfle et tronque


    A partir de la srie de dbits de lnonc il est facile de construire les trois sries de dbits demandes. On obtient : Une srie de 20 dbits de pointe annuels. Une valeur de dbit de pointe a t supprime du fait se sa trop grande faiblesse. Une srie gonfle comprenant 39 valeurs de dbit choisis parmi les deux plus grands dbits de pointe pour chaque anne ( srie gonfle ). Une srie de 48 dbits de pointe ayant une valeur suprieure 6.00 m3/s ( srie tronque ).

    Question 2 : Estimation du dbit de pointe de temps de retour 30 ans par la mthode statistique Ajustement de Gumbel
     Mthode appliquer : ajustement statistique dune srie de donnes distribution de Gumbel Lobjectif de cet exercice est destimer les dbits de pointes (dbits maximaux) correspondants un certain temps de retour, cest--dire une certaine probabilit dapparition donne. L'analyse frquentielle d'une longue srie de dbits maximaux permet destimer le temps de retour d'une valeur particulire. Cette prdiction repose sur la dfinition et la mise en uvre dun modle frquentiel qui est une quation dcrivant (modlisant) le comportement statistique dun processus. Ces modles dcrivent la probabilit dapparition dun vnement de valeur donne. Cest du choix du modle frquentiel (et plus particulirement de son type) que dpendra la validit des rsultats de lanalyse frquentielle. Un modle frquentiel trs souvent utilis pour dcrire le comportement statistique des valeurs extrmes est la distribution statistique de Gumbel (loi double exponentielle ou loi de Gumbel). La fonction de rpartition de la loi de Gumbel F(x) sexprime de la manire suivante :

    xa x a F ( x) = exp exp b (1) avec la variable rduite suivante : u = b


    o a et b sont les paramtres du modle de Gumbel. La distribution scrit alors de la manire suivante :

    (2)

    F ( x) = exp( exp(u ) ) (3)


    Mise jour le 08.04.2004

    et u = ln( ln(F ( x) )) .

    (4)

    HA 0801 - Page 1

    Lavantage dutiliser la variable rduite est que lexpression dun quantile est alors linaire ( x q = a + bu q ). En consquence, ds lors que les points de la srie ajuster peuvent tre reports dans un systme daxes x u , il est possible dajuster une droite qui passe le mieux par ces points et den dduire les deux paramtres a et b de la loi. Lestimation des paramtres a et b de lajustement peut se faire graphiquement (ajustement lil ou laide dune rgression statistique), ou selon une mthode mathmatique comme celle des moments (cf. ci-dessous). En pratique il sagit essentiellement destimer la probabilit de non dpassement F(xi) quil convient dattribuer chaque valeur xi. Il existe de nombreuses formules destimation de la fonction de ( x ) laide de la frquence empirique. Elles reposent toutes sur un tri de la srie par rpartition F valeurs croissantes permettant dassocier chaque valeur son rang r. Des simulations ont montr que pour la loi de Gumbel, il faut utiliser la frquence empirique de Hazen :

    F x[r ] =

    ( )

    r 0.5 n

    (5)

    o r est le rang dans la srie de donnes classe par valeurs croissantes, n est la taille de lchantillon, x[r] la valeur de rang r. Rappelons encore que le temps de retour T d'un vnement est dfini comme tant l'inverse de la frquence d'apparition de l'vnement. Soit : 1 T= 1 FQ ( xQ ) . (6) A laide de lajustement, il est alors possible destimer le dbit de pointe pour un temps de retour donn.  Mthode appliquer : Mthode des moments

    La mthode des moments consiste galer les moments des chantillons avec les moments thoriques de la loi 1 . Par la mthode des moments les paramtres a et b sont calculs daprs les formules :

    6 b = avec = 0.5772 (constante dEuler). (7) a . b =


    avec

    : cart-type des valeurs composant lchantillon. : moyenne de lchantillon.

    Ds lors il est possible destimer les dbits dont la reprsentation graphique est une droite dquation :

    u =a +b Q
    avec : u: variable rduite (cf. quation (4)).

    (8)

    Autrement dit : la mthode des moments consiste galer les caractristiques de la distribution (empirique) des chantillons avec les caractristiques thoriques de la loi. Les caractristiques utilises pour dcrire une distribution sont les moments dont les plus connus sont la moyenne et la variance (la moyenne dune distribution est appele premier moment, la variance est le deuxime moment). Les moments dordre plus levs sont utiles pour caractriser dautres aspects de la distribution. Le troisime moment est par exemple li lasymtrie ou la dyssimtrie.
    Mise jour le 08.04.2004

    HA 0801 - Page 2

     Dmarche et rsultats Etape 1 : Prparation de la srie de donnes des dbits de pointe. Trier les valeurs dans lordre croissant.

    Attribuer un rang chaque valeur.

    Etape 2 : Calcul de la frquence empirique pour chaque rang (Hazen, quation (5)). Etape 3 : Calcul de la variable rduite u du Gumbel (quation (4)). Etape 4 : Reprsentation graphique des couples (ui, xi) de la srie ajuster (figure 1).
    Ajustement des dbits maxima annuels selon une distribution de Gumbel et intervalle de confiance
    14

    12

    10

    dbit [m3/s]

    0 -2 -1 0 1 variable rduite de Gumbel u [-] 2 3 4

    dbit observ

    dbit estim

    intervalle de confiance 90%

    Figure 1. Ajustement graphique du modle (calcul des paramtres a et b de la droite dajustement de Gumbel par la mthode des moments)

    Etape 5 : Ajustement dune relation linaire de type x q = a + bu q aux couples (ui, xi) (figure 1).

    Avec un ajustement par la mthode des moments, on a alors pour ce premier ajustement une estimation des paramtres a1 et b1 : a1 = 167 et b1 = 38.2 Vrifier de manire visuelle lajustement. On peut remarquer que lajustement est mdiocre pour les grandes valeurs de u. A ce stade, il serait ncessaire de vrifier statistiquement que les valeurs observes sont estimes de manire satisfaisante laide de tests dajustement appropris.
    Etape 6 : Calcul de lintervalle de confiance laide de la formulation de lnonc (voir figure). Etape 7 : Utilisation du modle statistique pour estimer des dbits de pointe de diffrents temps de retour T. Pour T=30 ans, on suit les tapes suivantes : Calcul de la frquence de non-dpassement daprs la relation (6) :

    F (Q p (T )) = 1

    1 1 = 1 = 0,97 T 30

    Calcul de la variable rduite de Gumbel correspondante daprs la relation (4) :

    u = ln ln ( F (Q p (T )) ) = ln ( ln ( 0,97 ) ) = 3,38

    Mise jour le 08.04.2004

    HA 0801 - Page 3

    Calcul du quantile correspondant daprs la relation linaire (avec a et b fournis par ltape 5 prcdente ) :

    Qp (30) = a1 + b1 u30 = 8 + 1 0.97 = 11.2

    m3 /s

    Question 3 : Application aux sries tronque et gonfle


     Mthode appliquer : : ajustement statistique dune srie de donnes et corrections frquentielles (voir rappels-nonc).

    Lorsque la srie ajuster est gonfle ou tronque, la dmarche dcrites ci-dessus peut tre applique mais en apportant des corrections aux formulations des probabilits cumules de nondpassement et aux paramtres de la droite dajustement (voir rappels-nonc).
     Rsultats

    Les dbits de pointe de temps de retour 30 ans ainsi que les paramtres a et b pour les ajustements prcdents se trouvent dans le tableau 1. Les rsultats sont reprsents dans la figure 2. On peut constater que les valeurs de dbit sont peu diffrentes. Tableau 1 : Dbits de pointe de temps de retour 30 ans et paramtres a et b des ajustements des trois sries de dbit.
    Srie maxima annuels Srie gonfle Srie tronque

    a b Qp [m3/s]

    8 1 11,2

    8,1 0,9 11,3

    8 1,2 11,9

    Comparaison des ajustement des dbits


    14

    12

    10

    dbit [m3/s]

    0 -3 -2 -1 0 1 2 3 4 variable rduite de Gumbel u [-]

    dbit estim - maxima annuels

    dbit estim - srie gonfle

    dbit estim - srie tronque

    Figure 2: Comparaison des ajustements laide de la srie des maxima annuels, de la srie gonfle et tronque

    Mise jour le 08.04.2004

    HA 0801 - Page 4

    Vous aimerez peut-être aussi