Obtención Del Hidrograma Unitario

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 14

Obtención del Hidrograma Unitario

Universidad de sucre

Facultad de ingeniería

Departamento de ingeniería civil

Hidrología aplicada

Presentado a Luis Cerpa Reyes

Sincelejo, Sucre

28 de septiembre, 2017
Hidrograma unitario

“Si fuera posible que se produjeran dos lluvias idénticas sobre una cuenca

hidrográfica cuyas condiciones antes de la precipitación también fueran idénticas, sería de

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

ocurran lluvias idénticas; estas pueden variar su duración; el volumen precipitado; su

distribución espacial; su intensidad.”

Para la obtención del hidrograma unitario se tendrán en cuenta los siguientes caudales

de escurrimiento directo y precipitaciones uniformes:

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

Tiempo (horas) Pe (mm)


0 0
0,5 26,95
1 49,05
1,5 45,95
2 0

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

coordenadas del HU debe ser:

𝑈 =𝑁−𝑀+1

Donde U es el número de ordenadas del HU, M el número de ordenadas de Pe y N el

número de ordenadas de escurrimiento directo. Para este caso, se tiene que:

N=11, M=3
𝑈 =𝑁−𝑀+1

𝑈 = 11 − 3 + 1

𝑈=9

Por tanto, el hidrograma unitario debe tener 9 ordenadas.

Método de Ven Te Chow

Para esta metodología se emplea el procedimiento de eliminación Gaussiana,

sabiendo que el HU debe tener 9 ordenadas y que se cuenta con 3 precipitaciones se procede

a determinar la ecuación para cada una de las ordenadas.

𝑄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

precipitaciones mayores a la 3 se cancelarán, resultando las siguientes ecuaciones:

𝑄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

reemplazar los caudales y precipitaciones necesarias

𝑄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

Haciendo los reemplazos respectivos se obtienen los siguientes resultados.


i Ui Tiempo
1 0,44897959 0,5
2 1,20510393 1
3 2,60701446 1,5
4 2,77370741 2
5 1,67188634 2,5
6 0,49877995 3
7 0,36406926 3,5
8 0,42758681 4
9 0,07413302 4,5

HIDROGRAMA UNITARIO (V.T.C. HOW)


3.5

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

Las ecuaciones para tener en cuenta en este método son:

𝑄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

matrices P, U (de incógnitas) y Q.

Para solucionar este sistema se debe aceptar cierta aproximación en las componentes,

esto se mejora multiplicando U por la transpuesta de P, es decir PT.

La ecuación de matrices que resulta para obtener los valores de U es:

𝑃𝑇 𝑃{𝑈} = 𝑃𝑇 {𝑄}

Para resolver el sistema se debe obtener PTP y PTU

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

Ahora se multiplican las matrices para obtener PTP:

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

Se obtiene la matriz PT{Q}, siendo Q:

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

La solución está dada por

U = {P 𝑇 𝑄}{𝑃𝑇 P}−1

La matriz {𝑃𝑇 P}−1 es:

(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

Con la multiplicación de las matrices anteriormente obtenidas se obtiene la matriz de

solución U que es:

U
0,4495
1,2044
2,6069
2,7757
1,6684
0,5017
0,3650
0,4209
0,0848

Es decir, los valores para el hidrograma unitario son:

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

Con estos datos se obtiene el siguiente hidrograma unitario:

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

Para este método se hace uso de los siguientes datos:

Tiempo (horas) Caudal (m3/s)


1 1.6
2 3.0
3 4.2
4 5.3
5 3.5
6 2.0
7 1.1

Escurrimiento directo
6

4
Q(m3/s)

0
0 1 2 3 4 5 6 7 8
Tiempo (h)

Tiempo (horas) Precipitación (mm)


0 0
1 5
2 8
3 3
4 0

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

continuación, se presenta el proceso que se aplicada a cada una de ellas:

- Columna 1: Caudales de exceso (m3/s)


- Columna 2: Se asigna un U supuesto (m3/s mm)

- Columna 3 y 4: Se multiplican las ordenas del HU por las precipitaciones

menores.

- Columna 5: Suma de las columnas 3 y 4.

- Columna 6: Se restan las columnas 5 y 1.

- Columna 7: Se dividen cada una de las ordenadas de la columna 6 entre la

precipitación de excesos mayor (en valor absoluto). Esta columna representa las

ordenadas del HU siempre y cuando los valores obtenidos no difieran mucho de

los supuestos, de no ser así se prosigue con la columna 8.

- Columna 8: Promedio entre la columna 2 y 7, este sería el nuevo HU supuesto y

para la siguiente iteración seria la columna 2.

Se continua el proceso iterativo hasta obtener valores similares en la columna 2 y 7.

Para el paso Pe1=5mm, Pe2=8mm y Pe3=3mm.

HU HIDROGRAMA HIDROGRAMA HIDROGRAMA HU


Q QED SUPUESTO UxPe1 UxPe3 SUMA RESTA UNITARIO PROMEDIO
1 1,6 0,4 2 - 2 0,4 0,05 0,225
2 3 1 5 - 5 2 0,25 0,625
3 4,2 1,4 7 1,2 8,2 4 0,5 0,95
4 5,3 1,2 6 3 9 3,7 0,4625 0,83125
5 3,5 0,5 2,5 4,2 6,7 3,2 0,4 0,45
6 2 - - 3,6 3,6 1,6 0,2 -
7 1,1 - - 1,5 1,5 0,4 0,05 -
1 1,6 0,225 1,125 - 1,125 -0,475 0,059375 0,142187
2 3 0,625 3,125 - 3,125 0,125 0,015625 0,320312
3 4,2 0,95 4,75 0,675 5,425 1,225 0,153125 0,551562
4 5,3 0,83125 4,15625 1,875 6,03125 0,73125 0,09140625 0,461328
5 3,5 0,45 2,25 2,85 5,1 1,6 0,2 0,325
6 2 - - 2,49375 2,49375 0,49375 0,06171875 -
7 1,1 - - 1,35 1,35 0,25 0,03125 -
1 1,6 0,082812 0,41406 - 0,4140625 -1,1859375 0,14824218 0,234277
2 3 0,320312 1,60156 - 1,6015625 -1,3984375 0,17480468 0,363183
3 4,2 0,551562 2,75781 0,24843 3,00625 -1,19375 0,14921875 0,305273
4 5,3 0,461328 2,30664 0,96093 3,26757812 -2,03242187 0,25405273 0,289526
5 3,5 0,325 1,625 1,65468 3,2796875 -0,2203125 0,02753906 0,013769
6 2 - - 1,38398 1,38398437 -0,61601562 0,07700195 -
7 1,1 - - 0,975 0,975 -0,125 0,015625 -
1 1,6 0,234277 1,17138 - 1,17138671 -0,42861328 0,05357666 0,143927
2 3 0,363183 1,81591 - 1,81591796 -1,18408203 0,14801025 0,255596
3 4,2 0,305273 1,52636 0,70283 2,22919921 -1,9708007 0,24635009 0,275811
4 5,3 0,289526 1,44763 1,08955 2,53718261 -2,76281738 0,34535217 0,317439
5 3,5 0,013769 0,06884 0,91582 0,98466796 -2,51533203 0,31441650 0,164093
6 2 - - 0,86857 0,86857910 -1,13142089 0,14142761 -
7 1,1 - - 0,04130 0,04130859 -1,05869140 0,13233642 -
1 1,6 0,143927 0,71963 - 0,71963501 -0,88036499 0,11004562 0,126986
2 3 0,255596 1,27798 - 1,27798461 -1,72201538 0,21525192 0,235424
3 4,2 0,275811 1,37905 0,43178 1,81083984 -2,38916015 0,29864502 0,287228
4 5,3 0,317439 1,58719 0,76679 2,35398712 -2,94601287 0,36825161 0,342845
5 3,5 0,164093 0,82046 0,82743 1,64790039 -1,85209960 0,23151245 0,197802
6 2 - - 0,95231 0,95231781 -1,04768219 0,13096027 -
7 1,1 - - 0,49227 0,49227905 -0,60772094 0,07596511 -
1 1,6 0,126986 0,63493 - 0,63493156 -0,96506843 0,12063355 0,123809
2 3 0,23542 1,17712 - 1,17712211 -1,82287788 0,22785973 0,231642
3 4,2 0,287228 1,43614 0,38095 1,81710090 -2,38289909 0,29786238 0,292545
4 5,3 0,342845 1,71422 0,70627 2,42050046 -2,87949953 0,35993744 0,351391
5 3,5 0,197802 0,98901 0,86168 1,85069885 -1,64930114 0,20616264 0,201982
6 2 - - 1,02853 1,02853632 -0,97146368 0,12143296 -
7 1,1 - - 0,59340 0,59340820 -0,50659179 0,06332397 -
1 1,6 0,123809 0,61904 - 0,61904966 -0,98095033 0,12261879 0,123214
2 3 0,231642 1,15821 - 1,15821039 -1,84178960 0,2302237 0,230932
3 4,2 0,292545 1,46272 0,37142 1,83415675 -2,36584324 0,29573040 0,294137
4 5,3 0,351391 1,75695 0,69492 2,45188344 -2,84811655 0,35601457 0,353703
5 3,5 0,201982 1,00991344 0,87763 1,88754961 -1,61245038 0,20155629 0,201769
6 2 - - 1,05417 1,05417432 -0,94582567 0,11822821 -
7 1,1 - - 0,60594 0,60594806 -0,49405193 0,06175649 -

De lo anterior se obtiene que los valores para el hidrograma unitario son:

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)

También podría gustarte