Obtención Del Hidrograma Unitario
Obtención Del Hidrograma Unitario
Obtención Del Hidrograma Unitario
Universidad de sucre
Facultad de ingeniería
Hidrología aplicada
Sincelejo, Sucre
28 de septiembre, 2017
Hidrograma unitario
“Si fuera posible que se produjeran dos lluvias idénticas sobre una cuenca
esperarse que los hidrogramas correspondientes a las dos lluvias también fueran iguales.
Esta es la base del concepto de hidrograma unitario. En la realidad es muy difícil que
Para la obtención del hidrograma unitario se tendrán en cuenta los siguientes caudales
Tiempo
(horas) Q(m3/seg)
0 0
0,5 12,1
1 54,5
1,5 150
2 258
2,5 300,9
3 222,9
3,5 111,1
4 52,3
4,5 39,7
5 23,5
5,5 3,9
6 0
Escurrimiento directot
350
300
250
200
150
100
50
0
0 2 4 6 8
Hietograma de excesos
60
50
I(mm/hora)
40
30
20
10
0
0.5 1 1.5
∆t=0,5h
Tres de los métodos para calcular el HU para una cuenca son: el método matricial,
método de Ven Te Chow y método de Collins. Para los tres métodos el número de
𝑈 =𝑁−𝑀+1
N=11, M=3
𝑈 =𝑁−𝑀+1
𝑈 = 11 − 3 + 1
𝑈=9
sabiendo que el HU debe tener 9 ordenadas y que se cuenta con 3 precipitaciones se procede
𝑄1 = 𝑃1 𝑈1
𝑄2 = 𝑃2 𝑈1 + 𝑃1 𝑈2
𝑄3 = 𝑃3 𝑈1 + 𝑃2 𝑈2 + 𝑃1 𝑈3
𝑄4 = 𝑃4 𝑈1 + 𝑃3 𝑈2 + 𝑃2 𝑈3 + 𝑃1 𝑈4
𝑄5 = 𝑃5 𝑈1 + 𝑃4 𝑈2 + 𝑃3 𝑈3 + 𝑃2 𝑈4 + 𝑃3 𝑈5
𝑄6 = 𝑃6 𝑈1 + 𝑃5 𝑈2 + 𝑃4 𝑈3 + 𝑃3 𝑈4 + 𝑃2 𝑈5 + 𝑃1 𝑈6
𝑄7 = 𝑃7 𝑈1 +𝑃6 𝑈2 + 𝑃5 𝑈3 + 𝑃4 𝑈4 + 𝑃3 𝑈5 + 𝑃2 𝑈6 + 𝑃1 𝑈7
𝑄8 = 𝑃8 𝑈1 + 𝑃7 𝑈2 + 𝑃6 𝑈3 + 𝑃5 𝑈4 + 𝑃4 𝑈5 + 𝑃3 𝑈6 + 𝑃2 𝑈7 + 𝑃1 𝑈8
𝑄9 = 𝑃9 𝑈1 + 𝑃8 𝑈2 + 𝑃7 𝑈3 + 𝑃6 𝑈4 + 𝑃5 𝑈5 +𝑃4 𝑈6 + 𝑃3 𝑈7 + 𝑃2 𝑈8 + 𝑃1 𝑈9
Pero como solo se presentan 3 precipitaciones, todos los temrinos que contengan
𝑄1 = 𝑃1 𝑈1
𝑄2 = 𝑃2 𝑈1 + 𝑃1 𝑈2
𝑄3 = 𝑃3 𝑈1 + 𝑃2 𝑈2 + 𝑃1 𝑈3
𝑄4 = 𝑃3 𝑈2 + 𝑃2 𝑈3 + 𝑃1 𝑈4
𝑄5 = 𝑃3 𝑈3 + 𝑃2 𝑈4 + 𝑃1 𝑈5
𝑄6 = 𝑃3 𝑈4 + 𝑃2 𝑈5 + 𝑃1 𝑈6
𝑄7 = 𝑃3 𝑈5 + 𝑃2 𝑈6 + 𝑃1 𝑈7
𝑄8 = 𝑃3 𝑈6 + 𝑃2 𝑈7 + 𝑃1 𝑈8
𝑄9 = 𝑃3 𝑈7 + 𝑃2 𝑈8 + 𝑃1 𝑈9
Como se desea conocer los valores de U, y se conocen los de Q y P, para cada caso se
presenta una ecuación con una variable, que puede ser resulta despejando, a continuación se
presentan dichas ecuaciones y una tabla donde se resumen los resultados obtenidos al
𝑄1
𝑈1 =
𝑃1
𝑄2 − 𝑃2 𝑈1
𝑈2 =
𝑃1
𝑄3 − 𝑃3 𝑈1 − 𝑃2 𝑈2
𝑈3 =
𝑃1
𝑄4 − 𝑃3 𝑈2 − 𝑃2 𝑈3
𝑈4 =
𝑃1
𝑄5 − 𝑃3 𝑈3 − 𝑃2 𝑈4
𝑈5 =
𝑃1
𝑄6 − 𝑃3 𝑈4 − 𝑃2 𝑈5
𝑈6 =
𝑃1
𝑄7 − 𝑃3 𝑈5 − 𝑃2 𝑈6
𝑈7 =
𝑃1
𝑄8 − 𝑃3 𝑈6 − 𝑃2 𝑈7
𝑈8 =
𝑃1
𝑄9 − 𝑃3 𝑈7 − 𝑃2 𝑈8
𝑈9 =
𝑃1
2.5
u(m3/seg mm
1.5
0.5
0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Tiempo (horas)
Método Matricial
𝑄1 = 𝑃1 𝑈1
𝑄2 = 𝑃2 𝑈1 + 𝑃1 𝑈2
𝑄3 = 𝑃3 𝑈1 + 𝑃2 𝑈2 + 𝑃1 𝑈3
𝑄4 = 𝑃3 𝑈2 + 𝑃2 𝑈3 + 𝑃1 𝑈4
𝑄5 = 𝑃3 𝑈3 + 𝑃2 𝑈4 + 𝑃1 𝑈5
𝑄6 = 𝑃3 𝑈4 + 𝑃2 𝑈5 + 𝑃1 𝑈6
𝑄7 = 𝑃3 𝑈5 + 𝑃2 𝑈6 + 𝑃1 𝑈7
𝑄8 = 𝑃3 𝑈6 + 𝑃2 𝑈7 + 𝑃1 𝑈8
𝑄9 = 𝑃3 𝑈7 + 𝑃2 𝑈8 + 𝑃1 𝑈9
𝑄10 = 𝑃3 𝑈8 + 𝑃2 𝑈9
𝑄11 = 𝑃3 𝑈9
Este sistema debe ser organizado en forma vectorial, para lo cual se obtienen las
Para solucionar este sistema se debe aceptar cierta aproximación en las componentes,
𝑃𝑇 𝑃{𝑈} = 𝑃𝑇 {𝑄}
La matriz P es:
P
26,95 0 0 0 0 0 0 0 0
49,05 26,95 0 0 0 0 0 0 0
45,95 49,05 26,95 0 0 0 0 0 0
0 45,95 49,05 26,95 0 0 0 0 0
0 0 45,95 49,05 26,95 0 0 0 0
0 0 0 45,95 49,05 26,95 0 0 0
0 0 0 0 45,95 49,05 26,95 0 0
0 0 0 0 0 45,95 49,05 26,95 0
0 0 0 0 0 0 45,95 49,05 26,95
0 0 0 0 0 0 0 45,95 49,05
0 0 0 0 0 0 0 0 45,95
Y su transpuesta es:
PT
26,95 49,05 45,95 0 0 0 0 0 0 0 0
0 26,95 49,05 45,95 0 0 0 0 0 0 0
0 0 26,95 49,05 45,95 0 0 0 0 0 0
0 0 0 26,95 49,05 45,95 0 0 0 0 0
0 0 0 0 26,95 49,05 45,95 0 0 0 0
0 0 0 0 0 26,95 49,05 45,95 0 0 0
0 0 0 0 0 0 26,95 49,05 45,95 0 0
0 0 0 0 0 0 0 26,95 49,05 45,95 0
0 0 0 0 0 0 0 0 26,95 49,05 45,95
PT x P
5243,6 3575,7 1238,4 0 0 0 0 0 0
3575,7 5243,6 3575,7 1238,4 0 0 0 0 0
1238,4 3575,7 5243,6 3575,7 1238,4 0 0 0 0
0 1238,4 3575,7 5243,6 3575,7 1238,4 0 0 0
0 0 1238,4 3575,7 5243,6 3575,7 1238,4 0 0
0 0 0 1238,4 3575,7 5243,6 3575,7 1238,4 0
0 0 0 0 1238,4 3575,7 5243,6 3575,7 1238,4
0 0 0 0 0 1238,4 3575,7 5243,6 3575,7
0 0 0 0 0 0 1238,4 3575,7 5243,6
Q
12,1
54,5
150
258
300,9
222,9
111,1
52,3
39,7
23,5
3,9
Al multiplicar PT y Q se obtiene
PT{Q}
9891,82
20681,38
30523,76
31954,50
24147,55
13859,80
7383,68
4436,60
2401,80
U = {P 𝑇 𝑄}{𝑃𝑇 P}−1
(PT P)-1
4,71E-04 -5,00E-04 2,58E-04 1,27E-05 -1,57E-04 1,53E-04 -7,36E-05 4,13E-06 1,46E-05
-5,00E-04 1,00E-03 -7,74E-04 2,47E-04 1,74E-04 -3,14E-04 2,31E-04 -8,60E-05 4,13E-06
2,58E-04 -7,74E-04 1,14E-03 -7,67E-04 1,60E-04 2,60E-04 -3,55E-04 2,31E-04 -7,36E-05
1,27E-05 2,47E-04 -7,67E-04 1,13E-03 -7,47E-04 1,39E-04 2,60E-04 -3,14E-04 1,53E-04
-1,57E-04 1,74E-04 1,60E-04 -7,47E-04 1,13E-03 -7,47E-04 1,60E-04 1,74E-04 -1,57E-04
1,53E-04 -3,14E-04 2,60E-04 1,39E-04 -7,47E-04 1,13E-03 -7,67E-04 2,47E-04 1,27E-05
-7,36E-05 2,31E-04 -3,55E-04 2,60E-04 1,60E-04 -7,67E-04 1,14E-03 -7,74E-04 2,58E-04
4,13E-06 -8,60E-05 2,31E-04 -3,14E-04 1,74E-04 2,47E-04 -7,74E-04 1,00E-03 -5,00E-04
1,46E-05 4,13E-06 -7,36E-05 1,53E-04 -1,57E-04 1,27E-05 2,58E-04 -5,00E-04 4,71E-04
U
0,4495
1,2044
2,6069
2,7757
1,6684
0,5017
0,3650
0,4209
0,0848
i Ui Tiempo
1 0,4495 0,5
2 1,2044 1
3 2,6069 1,5
4 2,7757 2
5 1,6684 2,5
6 0,5017 3
7 0,3650 3,5
8 0,4209 4
9 0,0848 4,5
HIDROGRAMA UNITARIO
(MATRICIAL)
3.5000
3.0000
u(m3/seg mm
2.5000
2.0000
1.5000
1.0000
0.5000
0.0000
0 1 2 3 4 5
Tiempo (horas)
Método de Collins
Escurrimiento directo
6
4
Q(m3/s)
0
0 1 2 3 4 5 6 7 8
Tiempo (h)
Hietograma de excesos
9
8
7
6
P(mm)
5
4
3
2
1
0
0 1 2 3 4
∆t=1h
Para este caso U=5, y este debe ser el número de ordenadas del hidrograma unitario.
Para la aplicación de este método se hace uso de una tabla iterativa con 8 columnas, a
menores.
precipitación de excesos mayor (en valor absoluto). Esta columna representa las
i Ui Tiempo
1 0,1226 1
2 0,2302 2
3 0,2957 3
4 0,3560 4
5 0,2015 5
HIDROGRAMA UNITARIO
(COLLINS)
0.4
0.35
0.3
u(m3/seg mm
0.25
0.2
0.15
0.1
0.05
0
0 1 2 3 4 5 6
Tiempo (horas)