Remerciez-le!

Remerciez @Admin pour avoir partagé cet document gratuitement, de la manière la plus simple, en partageant sur les réseaux sociaux.

ENAC/ISTE/HYDRAM HYDROTHEQUE : base de données d’exercices en Hydrologie Cours

ENAC/ISTE/HYDRAM HYDROTHEQUE : base de données d’exercices en Hydrologie Cours : Hydrologie Générale / Thématique : Étude des crues ÉCOLE POLYTECHNIQUE FÉDÉRALE DE LAUSANNE Logo optimisé par J.-D.Bonjour, SI-DGR 13.4.93 Exercice no HA 0810 - Corrigé Régionalisation des débits de pointe pour une région hydrologique relativement homogène Données de l’exercice Le fichier Excel « HA0810_enonce.xls » donne les principales caractéristiques hydrologiques et physiographiques des quatre groupes de bassins versants (G0, G1, G2 et G3). L’exercice porte sur trois bassins versants : La Mentue à Yvonnand (surface : 105 km2, altitude moyenne : 679 m). Le Lyssbach (surface : 22.5 km2, altitude moyenne : 616 m). Un bassin versant hypothétique non jaugé (surface : 100 km2, altitude moyenne : 700 m). Le fichier « HA0810_enonce.xls » donne les séries de crues annuelles des deux premiers bassins. Le fichier « L-moments-théoriques » donne quelques informations sur les lois théoriques (Gumbel, GEV, lognormale, etc.). Il existe aussi dans le fichier Excel « HA0810_feuilledecalcul.xls » des feuilles de calcul compléter. Les résultats sont disponibles sur le fichier Excel « HA0810_corrige.xls ». Question 1 : Méthode de « Index Flood » Rappelons que trois étapes sont nécessaires pour aboutir à une estimation des débits de pointe selon la méthode index flood: Délimitation de régions (sur la base de critères géographiques, physiographiques, climatiques, etc.) et vérification de leur homogénéité en vérifiant que certaines statistiques des débits de pointe annuels observés (coefficients de variation et d'asymétrie) peuvent être considérées comme les mêmes pour tous les bassins de la région. Pour chacune de ces régions, ajustement d’une loi statistique à la distribution normalisée des débits de pointe annuels (i.e. débit de pointe divisé par le débit de pointe annuel moyen). Cette distribution normalisée, représentative de la région est aussi appelée "courbe de croissance". Enfin, détermination de l'expression permettant d'obtenir pour chacun des sites le débit de pointe annuel moyen (facteur d’échelle) à partir de certaines caractéristiques physiographiques ou climatiques des bassins versants. Cette expression peut être obtenue notamment par régression simple ou multiple. La distribution des débits de pointe recherchée au site i est alors le produit du débit de pointe annuel moyen du site i et de la courbe de croissance de la région correspondante. Question 2 : Evaluation de l’homogénéité et choix d’une loi théorique sur la base du diagramme des L-moments. Les diagrammes des L-moments (ou plus exactement les diagrammes des L-coefficients de variation t2, d’asymétrie t3 et d’aplatissement t4) sont similaires aux diagrammes des moments classiques, mais ont l’avantage d’être moins sensibles aux fluctuations d’échantillonnage. Ils permettent d’évaluer le degré d’homogénéité d’une région et de choisir une distribution théorique pour construire une courbe de croissance. Mise à jour le 14.04.2004 HA 0810 - Page 1 On admettra que la région formée par les bassins du groupe G0 est homogène. Autrement dit, si ces bassins disposaient chacun d’une série de mesures infiniment longue, ils seraient tous représentés par le même point dans ces diagrammes. La dispersion des points observée est donc attribuée aux fluctuations d’échantillonnage. Les diagrammes des L-moments construits pour les groupes G0, G1, G2 et G3 sont présentés ci- dessous. 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 -0.2 0.0 0.2 0.4 t3 t2 G0 G1 G2 G3 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 -0.2 0.0 0.2 0.4 t3 t4 G0 G1 G2 G3 Sur la base de ces diagrammes, on peut faire quelques remarques : a) Groupe G1 : La dispersion des points du groupe G1 (bassins disposant de moins de 10 ans de mesure) dans les deux diagrammes ci-dessus est plus grande que celle du groupe G0. Les causes principales de cette dispersion sont les fluctuations d’échantillonnage importantes dues à des séries de mesures très courtes sur les bassins du groupe G1. L’acquisition de nouvelles données aux stations du groupe G1 aura pour effet de réduire cette dispersion. On ne peut donc pas conclure, dans ce cas, que la dispersion est due à des réponses hydrologiques particulières, et donc à un phénomène d’hétérogénéité du groupe G1. De même, il n’est pas possible d’affirmer que l’incorporation des bassins du groupe G1 dans le groupe G0 augmente l’hétérogénéité de ce dernier. b) Groupe G2 : Pour les bassins du groupe G2 (bassins préalpins) on observe un décalage du nuage de points vers des valeurs de t3 et t4 plus grandes (les valeurs de t2 restant sensiblement identiques). Contrairement à celles du groupe G1, les séries de mesures du groupe G2 ne sont pas significativement plus courtes que celle du groupe G0. L’incorporation des bassins du groupe G2 au groupe G0 a donc pour effet d’augmenter l’hétérogénéité (augmentation de la taille du nuage de points non imputable aux fluctuations d’échantillonnage). Cependant les valeurs du groupe G2 sont peu dispersées (de même ordre de grandeur que celles du groupe G0). Compte tenu des données à disposition, on peut s’attendre à ce que les bassins du groupe G2 appartiennent à une autre région homogène (une vérification formelle de l’homogénéité peut se faire au moyen de tests statistiques comme ceux de Chowdhury ou de Hosking et Wallis). Les différences au niveau des L-coefficients entre le groupe G0 et G2 sont le reflet des différences au niveau des réactions hydrologiques entre les bassins versants du Plateau et ceux des Préalpes. En effet, les bassins préalpins sont caractérisés par : (1) la quasi-imperméabilité des surfaces, (2) des pentes élevées (diminution du temps de concentration) et (3) des températures plus basses que sur le plateau (précipitations plus souvent sous forme de neige en hiver). Les crues des bassins préalpins se produisent ainsi généralement en été à la suite de précipitations convectives. En revanche, sur le plateau, les crues annuelles sont essentiellement provoquées par des précipitations hivernales de type frontal. Mise à jour le 14.04.2004 HA 0810 - Page 2 c) Groupe G3 : Enfin, pour les bassins du groupe G3 (petits bassins versants du plateau) on observe que les points sont plus dispersés que ceux du groupe G0, sans pour autant que l’on puisse y observer un décalage dans l’un ou l’autre des axes t2, t3 ou t4 (contrairement aux bassins du groupe G2). Par ailleurs la dispersion plus élevée n’est pas imputable à des fluctuations d’échantillonnage importantes (les séries de mesures du groupe G3 ne sont pas significativement plus courtes que celle du groupe G0). Les bassins du groupe G3 ne constituent donc apparemment pas une région homogène et, l’inclusion de ces bassins dans le groupe G0 va augmenter l’hétérogénéité de ce dernier. Ceci s’explique par le fait que la variabilité des caractéristiques (et par conséquent des réponses hydrologiques) des petits bassins (moins de 10 km2) est en général beaucoup plus grande que pour les grands bassins versants. d) Lois théoriques pour le groupe G0 : La figure ci-dessous montre le diagramme des L-coefficients d’asymétrie (t3) et d’aplatissement (t4) pour le groupe G0. On y a également représenté les valeurs théoriques de certaines lois de distributions utilisées pour l’analyse fréquentielle (Gumbel, GEV, Exponentielle, etc.). On y observe que les distributions de Gumbel (à deux paramètres) et GEV (à trois paramètres) constituent de bons candidats pour modéliser la courbe de croissance des bassins du groupe G0. Le choix entre la loi de Gumbel et la loi GEV revient à décider si l’on fixe le paramètre k à 0 ou non. Une vérification formelle de l’adéquation peut se faire au moyen de tests statistiques comme ceux de Chowdhury ou de Hosking et Wallis ou par le test de Kolmogorov-Smirnov (présenté par exemple dans [Meylan et Musy, 1998]). 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 -0.2 0.0 0.2 0.4 t3 t4 G0 normale exp gumbel uniforme GEV GPA logN P3 Question 3. Construction d'une courbe de croissance Pour construire la courbe de croissance des bassins du groupe G0, une des approches possibles est de calculer des valeurs centrales des L-coefficients de variation et de déterminer les paramètres de la loi théorique retenue (cf. cours). Pour calculer ces valeurs centrales et afin de tenir compte de l'effet de l'incertitude associé au nombre d'années de mesures à disposition, la moyenne pondérée (où les poids sont inversement proportionnels à la longueur des séries observées) est préférable à la moyenne arithmétique simple. On a ainsi pour le groupe G0: t2 t3 t4 Valeurs pondérées = 0.221 0.171 0.174 On obtient les paramètres suivants : Pour une loi de Gumbel : k =0, b = t2 / ln(2) = 0,319, et a = 1 - 0,5772 . b = 0.815 Pour une loi de distribution GEV les paramètres estimés sont très peu différents d’une distribution de Gumbel qui est donc tout à fait appropriée dans ce cas. Mise à jour le 14.04.2004 HA 0810 - Page 3 La figure ci-dessous donne la courbe de croissance. Cette courbe permet de lire le facteur de croissance, qui est le rapport entre le débit de pointe de temps de retour T et le débit de pointe annuel moyen. Pour un temps de retour 100 ans, il uploads/Geographie/ ha0810-corrige 1 .pdf