Théorie et méthodes de l'hydrogramme unitaire
Théorie et méthodes de l'hydrogramme unitaire
•Sherman - 1932
•Horton - 1933
•Wisler & Brater - 1949 - "l'hydrogramme de surface"
ruissellement résultant d'une pluie relativement courte et intense,
appelé une tempête unitaire.
Le hydrogramme de ruissellement peut être "composé" de ruissellement qui est
généré comme un flux à travers le sol (Black, 1990).
1
Théorie de l'hydrogramme unitaire
2
Hydrogramme Unitaire "Lingo"
Durée
•Lag Time
Temps de concentration
Limbre Montant
Limbre de récession (limbre descendant)
Débit de pointe
3
Représentation graphique
Durée de
excès
précipitation.
Lag time
Temps de
concentration
Écoulement de base
4
Méthodes de développement des UHG
De manière synthétique
–Snyder
–SCS
–Temps-Aire (Clark, 1945)
Distributions "ajustées"
Géomorphologique
5
Hydrogramme unitaire
6
Hydrogramme unitaire dérivé
700.0000
600 000,00
500.0000
Surface
Réponse
400.0000
Débit de base
300.0000
200.0000
100.0000
0,0000
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 6 2 8 4 0 6 2 8 4 0 6 2 8 4 0 6 2 8 4 0 6 2 8
.0 .1 .3 .4 .6 .8 .9 .1 .2 .4 .6 .7 .9 .0 .2 .4 .5 .7 .8 .0 .2 .3 .5 .6
0 0 0 0 0 0 0 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3
7
Hydrogramme unitaire dérivé
700.0000
600.0000 Total
Hydrographe
500.0000
Surface
400.0000 Réponse
300.0000
Débit de base
200.0000
100,0000
0,0000
0,0000 0,5000 1.0000 1.5000 2,0000 2.5000 3.0000 3.5000 4.0000
8
Hydrogramme Unitaire Dérivé
Règles empiriques :
… la tempête devrait être assez uniforme par nature et le
les précipitations excessives devraient être également uniformes tout au long
le bassin. Cela peut nécessiter les conditions initiales tout au long
le bassin doit être spatialement similaire.
… Deuxièmement, la tempête devrait être relativement constante dans le temps,
ce qui signifie qu'il ne doit pas y avoir de pauses ou de périodes sans
précipitation.
… Enfin, la tempête devrait produire au moins un pouce de
précipitations excessives (la zone sous l'hydrogramme après
correction pour le flux de base).
9
Dériver un UHG d'une tempête
bassin versant échantillon = 450 mi2
25000 0,8
0,7
20000
0,6
0,5
15000
0,4
10000
0,3
0,2
5000
0,1
0 0
0 8 4 2 0 8
16 24 32 40 48 56 64 72 80 88 96 0 1 2 2
1 1 1 1
Time (hrs.)
10
Séparation du débit de base
... généralement admis que le point d'inflexion sur le flanc de la récession
d'un hydrogramme est le résultat d'un changement dans le contrôle physique
processus de l'excès de précipitations s'écoulant vers la sortie du bassin.
Dans cet exemple, le débit de base est considéré comme une ligne droite
connecter ce point où l'hydrographe commence à monter rapidement
et le point d'inflexion du côté récession de l'hydrogramme.
le point d'inflexion peut être trouvé en traçant l'hydrogramme en semi-
mode logarithmique avec le flux tracé sur l'échelle logarithmique et en notant le temps
où le côté de la récession s'ajuste à une ligne droite.
11
Diagramme semi-logarithmique
100000
1000
100
10
1
4 9 4 9 4 9 4
29 34 39 44 49 54 59 64 69 74 79 84 89 94 99 0 0 1 1 2 2 3
1 1 1 1 1 1 1
Time (hrs.)
12
Hydrograph & Baseflow
25000
20000
15000
10000
5000
Time (hrs.)
13
Écoulement de base séparé
25000
20000
15000
10000
5000
0
0 7 5 2 9 6 3
14 21 28 35 42 49 56 63 70 77 84 91 98 0 1 1 2 3
1 1 1 1 1
Temps (hrs.)
14
Calculs d'échantillons
•Dans l'exemple présent (pas de temps horaire), les flux sont sommés
et ensuite multiplié par 3600 secondes pour déterminer le volume de
écoulement en pieds cubes. Si souhaité, cette valeur peut ensuite être convertie
en pieds-acre en divisant par 43 560 pieds carrés par acre.
La profondeur des ruissellements directs en pieds se trouve en divisant le total
volume de précipitations excédentaires (maintenant en acres-pieds) par le
zone de bassin versant (450 mi2converti à 288 000 acres.
•Dans cet exemple, le volume de précipitations excédentaires ou direct
Le ruissellement pour la tempête n°1 a été déterminé à 39 692 acres-pieds.
La profondeur du ruissellement direct est trouvée à 0,1378 pieds après division.
par le bassin versant de 288 000 acres.
Enfin, la profondeur du ruissellement direct en pouces est de 0,1378 x 12 = 1,65
pouces.
15
Obtenir les ordonnées UHG
16
UHG final
25000
Hydrographe de la tempête n°1
20000 hydrographe
15000
Tempête n° 1 unité
hydrogramme
10000
Tempête n°1
écoulement de base
5000
Time (hrs.)
17
Déterminer la durée de l'UHG
La durée de l'hydrogramme unitaire dérivé est trouvée par
examinant les précipitations pour l'événement et déterminant que
précipitation qui est en excès.
Cela est généralement accompli en traçant les précipitations dans
forme d'hyetographe et dessiner une ligne horizontale de sorte que le
les précipitations au-dessus de cette ligne sont égales à la profondeur de l'excès
précipitation comme précédemment déterminé.
Cette ligne horizontale est généralement appelée leF-index et est
basé sur l'hypothèse d'une infiltration constante ou uniforme
taux.
L'infiltration uniforme nécessaire pour provoquer 1,65 pouces de
Les précipitations excessives ont été estimées à environ 0,2.
pouces par heure.
18
Estimation des précipitations excessives.
0,8
0,7
0,6
0,5
0,3
0,2
0,1
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
Time (hrs.)
19
Précipitations excessives
0,9
Petites quantités de
0,6
précipitations excessives à
beginning and end may
0,5
être omis.
0,4
0,3
0,2
0,1
0
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
Time (hrs.)
20
Changer la durée
•Très souvent, il sera nécessaire de changer la durée de l'unité
hydrogramme.
•Si les hydrogrammes unitaires doivent être moyennés, alors ils doivent être de la
même durée.
• De plus, la convolution de l'hydrogramme unitaire avec une précipitation
l'événement nécessite que la durée de l'hydrographe unitaire soit égale
au pas de temps des précipitations incrémentales.
La méthode la plus courante pour modifier la durée d'une unité
L'hydreographe est par la méthode de la courbe en S.
60000,00
50000,00
40000,00
30 000,00
20000,00
10000,00
0,00
Time (hrs.)
22
Convertir à une durée d'1 heure
Pour obtenir un hydrogramme unitaire de 1 heure, la courbe S est décalée de 1
heure et la différence entre les deux courbes S retardées est trouvée à
être un hydrographe unitaire d'une heure.
23
Dernière heure UHG
14000,00 60000,00
12000,00
50000,00
Les courbes en S sont
10000,00 retardé d'1 heure
et la différence 40000,00
est trouvé.
8000,00 unité d'une heure
hydrogramme résultant
30000,00
du S en retard -
6000,00 courbes et
multiplier le
différence de 6. 20000,00
4000,00
10000,00
2000,00
0,00 0,00
Time (hrs.)
24
Méthode de raccourci
•Par exemple, si vous aviez un hydrogramme unitaire de deux heures et que vous
voulu le changer en un hydrographe unitaire de quatre heures.
25
Exemple de méthode abrégée
•Tout d'abord, un hydrogramme unitaire de deux heures est donné et un hydrogramme unitaire de quatre heures
un hydrogramme est nécessaire.
Il y a deux possibilités, développer la courbe en S ou puisque elles sont
les multiples utilisent la méthode raccourcie.
Time(hr) Q
0 0
1 2
2 4
3 6
4 10
5 6
6 4
7 3
8 2
9 1
10 0
26
Exemple de méthode abrégée
Time(hr) Q Déspa
l cé
U
HG
0 0
1 2
2 4 0
3 6 2
4 10 4
5 6 6
6 4 10
7 3 6
8 2 4
9 1 3
10 0 2
11 1
12 0
27
Exemple de méthode raccourcie
Time(hr) Q Dp
si a
l cé
U
HG Somme
0 0 0
1 2 2
2 4 0 4
3 6 2 8
4 10 4 14
5 6 6 12
6 4 10 14
7 3 6 9
8 2 4 6
9 1 3 4
10 0 2 2
11 1 1
12 0 0
28
Exemple de méthode abrégée
31
Snyder
•Puisque le débit de pointe et le temps de débit de pointe sont deux des plus importants
paramètres caractérisant un hydrogramme unitaire, la méthode de Snyder
employe des facteurs définissant ces paramètres, qui sont ensuite utilisés dans
la synthèse du graphe unitaire (Snyder, 1938).
•Les paramètres sont Cp, le facteur de débit de pointe, et Ct, le facteur de retard.
•L'hypothèse de base dans cette méthode est que les bassins qui ont
des caractéristiques physiographiques similaires sont situées dans la même zone
aura des valeurs de C similairestet Cp.
•Par conséquent, pour les bassins non mesurés, il est préférable que le bassin soit proche
ou similaire à des bassins jaugeés pour lesquels ces coefficients peuvent être
determined.
32
Relations de base
t DÉCALAGE
tdurée
5,5
t DÉCALAGE
tbase 3
8
640ACp
q sommet
t LAG
33
Forme finale
La forme finale de l'hydrogramme unitaire de Snyder est contrôlée
par les équations pour la largeur à 50 % et 75 % du pic de
le UHG :
34
SCS
Rapports de flux
Messecum
0,8
0,6
0,4
0,2
0
0 0,5 1 1,5 2 2.5 3 3.5 4 4,5 5
T/Tpeak
35
Rapports sans dimension
Rapports de temps Rapports de décharge Rapports de courbe de masse
(t/tp) (q/qp) (Qa/Q)
0 .000 .000
.1 .030 .001
.2 .100 .006
.3 .190 .012
.4 .310 .035
.5 .470 .065
.6 .660 .107
.7 .820 .163
.8 .930 .228
.9 .990 .300
1.0 1.000 .375
1.1 .990 .450
1.2 .930 .522
1.3 .860 .589
1.4 .780 .650
1,5 .680 .700
1.6 .560 .751
1,7 .460 .790
1.8 .390 .822
1.9 .330 .849
2.0 .280 .871
2.2 .207 .908
2.4 .147 .934
2,6 .107 .953
2.8 .077 .967
3.0 .055 .977
3.2 .040 .984
3.4 .029 .989
3.6 .021 .993
3.8 .015 .995
4.0 .011 .997
4,5 .005 .999
5.0 .000 1.000
36
Représentation triangulaire
D SCS sans dimension UHG & Représentation triangulaire
1.2 Excès
Precipitation
Tretard
0,8
Rapports de flux
Messe C
Triangulaire
0,6
Point de
Flexion
Tc
0,4
0,2
0
0,0 Tp 1.0 2.0 3.0 4.0 5.0
Tb
T/Tpointe
37
Représentation Triangulaire
D Représentation UHG sans dimension SCS et triangulaire
Tlatence
Rapports de flux
Cum. Mass
Triangulaire
0,6
Point de
Flexion
Tc
q p Tp q p Trq
0,4
Q= + = p( TpT
+r )
2 2 2 0,2
0
0,0 Tp 1,0 2.0 3.0 4.0 5.0
Tb
2Q
T/Tpointe
qp =
+r
Tp T
654,33 x 2 x A x Q
qp= Le 645,33 est la conversion utilisée pour
TpT+r
livrer 1 pouce de ruissellement (la zone
sous l'hydrographe unitaire) d'1 kilomètre carré
484AQ
qp= mile en 1 heure (3600 secondes).
Tp
38
484 ?
484AQ
qp=
Tp
39
Durée et horaires ?
Encore depuis le triangle
D
Tp = + L
2
L = Lag time
L 0,6*Tc
Tc D 1,7Tp
D
+0.6 T=cTp
2
40
Temps de concentration
Équations de régression
Approche segmentale
41
Une équation de régression
L0,8(S 1) 0,7
Tretard
1900(%Pente) 0,5
42
Approche segmentale
Plus "hydraulique" par nature
•Le paramètre estimé est essentiellement le temps de
concentration ou temps de trajet le plus long dans le bassin.
En général, le temps de trajet le plus long correspond au plus long
chemin de drainage
Le chemin d'écoulement est divisé en segments avec le flux dans chaque segment.
être représenté par un type de régime d'écoulement.
Les représentations de flux les plus courantes sont les flux de surface, les draps, les rigoles et
43
Une approche de base
1
V kS 2
0,25 Forêt avec un sol lourd de débris, prairie de foin (ruissellement de surface)
0,5 Cultivation en culture de couverture ou en semis direct ; en courbes de niveau ou en bandes
cultivé ; forêt (écoulement de surface) McCuen (1989) et SCS
0,7 Pâturage de courte herbe (écoulement de surface) (1972) fournir des valeurs de k
0,9 Rangée cultivée droite (écoulement de surface) pour plusieurs situations d'écoulement
1.0 Presque nu et non cultivé (écoulement de surface) ; cônes alluviaux dans (pente en %)
régions montagneuses de l'ouest
1.5 Voie d'eau herbacée
2.0 Zone pavée (écoulement en feuille); petits ravins de terre
Type d'écoulement K
Petite Affluent - Permanent ou intermittent 2.1
des flux qui apparaissent comme solides ou en pointillés
lignes bleues sur les cartes topographiques de l'USGS.
Voie navigable - Tout itinéraire de circulation terrestre qui 1.2 Sorell & Hamilton, 1991
est un fossé bien défini par l'élévation
contours, mais ce n'est pas une section de ruisseau comme
défini ci-dessus.
Écoulement de feuillets - Tout autre chemin d'écoulement en surface 0,48
qui ne correspond pas à la définition de
une voie navigable.
44
Forme triangulaire
En général, on peut dire que la version triangulaire ne causera pas ou
introduire des différences notables dans la simulation d'un événement de tempête,
particulièrement lorsque l'on s'intéresse au débit de pointe.
Pour les simulations à long terme, l'hydrographie unitaire triangulaire a
un impact potentiel, en raison de la forme du bras de récession.
•Le Corps des ingénieurs de l'armée des États-Unis (HEC 1990) ajuste une unité Clark
hydrogramme pour correspondre aux débits de pointe estimés par l'unité de Snyder
procédure hydrographique.
Il est également possible d'ajuster une fonction synthétique ou mathématique à la
paramètres de débit de pointe et de timing de l'unité hydrographique souhaitée.
• Aron et White (1982) ont ajusté une distribution de probabilité gamma en utilisant
données de débit de pointe et de temps jusqu'au pic.
45
Ajustement d'une distribution Gamma
t a e t b
f(t;a,b)
b a1 (a 1)
500.0000
450.0000
400.0000
350.0000
300.0000
250.0000
200.0000
150.0000
100.0000
50,0000
0,0000
0,0000 1.0000 2.0000 3,0000 4,0000 5,0000 6,0000
46
Temps-Surface
47
Temps-Surface
100%
Temps
Q % Zone de conc.
Temps Temps
48
Temps-Surface
49
Exemple hypothétique
•A 190 mi2 Le bassin versant est divisé en 8 isochrones de temps de trajet.
Le coefficient de routage du réservoir linéaire, R, estimé à 5,5 heures.
Un intervalle de temps de 2,0 heures sera utilisé pour les calculs.
50
Règle générale
51
Analyse des bassins
52
Zone incrémentale
40
35
30
25
20
15
10
0
1 2 3 4 5 6 7 8
Time Increment (hrs)
53
Courbe temporelle de surface cumulative
0
0 20 40 60 80 100 120 140 160 180 200
Time (hrs)
54
Difficulté à obtenir un temps-aire
Courbe?
TAje 1,414T1,5
je pour (0 Ti 0,5)
55
UHG instantané
IUHje cIje (1 c)IUH(je1)
2 t
c
2R t
56
Calculs
Temps Inc. Inc. Inst. IUHG 2 heures
(heures) Zone Traduit UHG Retardé 2 UHG
(mi2) Flow (cfs) heures (cfs)
(1) (2) (3) (4) (5) (6)
0 0 0 0 0
2 14 4 515 1391 0 700
4 44 14 190 5333 1 391 3,360
6 53 17 093 8955 5,333 7 150
8 79 25,478 14043 8 955 11 500
10 0 0 9717 14 043 11 880
12 6724 9,717 8 220
14 4653 6 724 5 690
16 3220 4 653 3 940
18 2228 3,220 2 720
20 1542 2,228 1 890
22 1067 1,542 1 300
24 738 1,067 900
26 510 738 630
28 352 510 430
30 242 352 300
32 168 242 200
34 116 168 140
36 81 116 100
38 55 81 70
40 39 55 50
42 26 39 30
44 19 26 20
46 13 19 20
48 13
57
Zones incrémentales
90
80
70
60
50
40
30
20
10
0
0 2 4 6 8 10
Increments de temps (2 heures)
58
Flux incrémentaux
30000
25000
20000
15000
10000
5000
0
1 2 3 4 5 6
Time Increments (2 hrs)
59
UHG instantané
16000
14000
12000
10000
8000
6000
4000
2000
0
0 10 20 30 40 50 60
Temps (heures)
60
Lag & Average
16000
14000
12000
10000
8000
6000
4000
2000
0
0 10 20 30 40 50 60
Time (hrs)
61
Géomorphologique
62
Classement des ruisseaux de Strahler
1
1
1 1
2
2 1
1 2
63
Concepts de probabilité
64
Équations GIUH
• Basé sur le canal, triangulaire instantané
hydrogramme unitaire
1.31 0,44LWR B
0,55
qp R 0L.43V tp RL0,38
LW Réalité VirtuelleU n
65