Hidraulica e Hidrologia - Capitulo 6 - Q Max

DISEO Y PUESTA EN MARCHA DE LA RED DE MONITOREO AMBIENTAL EN LA CUENCA HIDROGRFICA DEL
RO MEDELLN EN JURISDICCION DEL AREA METROPOLITANA

UNIVERSIDAD DE ANTIOQUIA - UNIVERSIDAD PONTIFICIA BOLIVARIANA UNIVERSIDAD DE MEDELLN UNIVERSIDAD NACIONAL
METODOLOGAS
(MXIMOS)
PARA
LA
ESTIMACIN
DE
CAUDALES
EXTREMOS
6.1 MTODO GRADEX

El mtodo GRADEX (gradiente de valores extremos) se emplea para determinar los
caudales mximos asociados a determinado perodo de retorno, Tr, a partir de la
distribucin de frecuencias de los valores mximos anuales de precipitacin. El mtodo
supone que en condiciones de saturacin de una cuenca, la distribucin de probabilidad
de la precipitacin mxima anual es paralela a la distribucin de probabilidad del caudal
mximo.
La aplicacin del mtodo requiere informacin anual de precipitacin mxima y caudal.
En las siguientes secciones se describe la metodologa empleada en su aplicacin.
6.1.1. Valores mximos de precipitacin
Adicional a las hiptesis presentadas al inicio del captulo, el mtodo supone que la
precipitacin mxima que cae sobre una cuenca se distribuye probabilsticamente de
forma exponencial. En este tipo de distribucin pueden mencionarse la Gumbel y la
LogNormal, ampliamente utilizadas en nuestro medio en el estudio de valores extremos.
De esta forma, la precipitacin mxima correspondiente a un perodo de retorno, Tr,
puede determinarse segn la ecuacin 6.1 para la distribucin Gumbel.
1
P = + ln ln(1 )
Tr
(6.1)
donde:
= 0.577
(6.2)
(6.2)
P es la precipitacin mxima sobre la cuenca, y corresponden a la media y desviacin

estndar de sta, estimadas con la serie de datos disponible.
La pendiente de la distribucin definida por las ecuacin (6.1), , es denominada Grdex
de precipitacin, y es uno de los parmetros que finalmente permite estimar las avenidas
mximas.
REA METROPOLITANA DEL VALLE DE ABURRA
CONVENIO 366 / 2003
6-1

6.1.2. Saturacin
Los caudales mximos que se obtienen a partir del mtodo GRADEX se asocian a lluvias
de Tc (tiempo de concentracin) horas de duracin, que son las que hacen que se
cumplan las condiciones de saturacin en una cuenca cuando estn asociadas a perodos
de retorno de 10 aos, segn Guillot y Duband (1967). En Colombia, sin embargo, se ha
encontrado que las condiciones de saturacin se dan para perodos de retorno menores
(Smith & Vlez, 1997) que se encuentran entre 1.5 y 2.5 aos.
Una vez se presenta la saturacin, cualquier exceso de lluvia se transforma en
escorrenta, de modo que si se presenta una lluvia de duracin Tc con un perodo de
retorno Tr puede hacerse el siguiente balance volumtrico para dicho exceso:
V = 1000 P A = 3600 Q tc
(6.3)
donde V (m) es el volumen de agua que ingresa a la cuenca en forma de precipitacin,

A (Km) es el rea de sta, Tc (h) es la duracin de la lluvia, Q (m/s) es el incremento en
el caudal producido por el incremento P (mm) asociado a una lluvia de perodo retorno
Tr, donde Tr es mayor que el perodo de retorno de la lluvia que ocasiona la saturacin de
la cuenca. De acuerdo con la ecuacin 6.1, P viene dado por:
T
P = ln r
T0
(6.4)
donde T0 corresponde al perodo de retorno de saturacin.

De esta manera, el incremento de caudal Q puede expresarse como:
Q tc =
A Tr
ln
3.6t c T0
(6.5)
El Grdex de caudal, q, expresado en m/s se obtiene a partir del Grdex de

precipitacin, , a partir de la expresin:
q =
A
3.6tc
(6.6)
Obtenidos el caudal pivote Qtc (T0 aos), como es conocido en la literatura, y el Grdex de
caudal q se pueden extrapolar los caudales para diferentes perodos de retorno, T, como:
T
Q tc (T ) = Q tc (T0 ) + q Ln
T0
(6.7)
CONVENIO 366 / 2003
6-2

6.1.3. Obtencin del caudal pivote

Como ya se mencion, la aplicacin del mtodo requiere la determinacin del caudal
correspondiente a T0 aos de perodo de retorno, estado en cual la cuenca se satura. La
metodologa que mayor acogida ha tenido para ste propsito, plantea la obtencin de
dicho caudal a partir de la consideracin de banca llena. Los investigadores han
encontrado resultados consistentes para el perodo de retorno del caudal a partir de ste
mtodo: para Wolman y Leopold (1957), el perodo de retorno del caudal a seccin llena
vara entre 1 y 2 aos, Dury et al (1973) postulan un perodo de 1.58 aos (mejor ajuste
en Norteamrica), Nixon (1959) con datos del Reino Unido propone 1.67 aos. Segn los
resultados para el departamento de Antioquia de Ramrez y Giraldo (2001), trabajar con
caudales pivotes de perodos de retorno iguales a 2.33 y 1.3 aos provoca errores
inferiores al 25%.
6.1.4. Caudales mximos en cuencas no instrumentadas
La determinacin de caudales mximos en cuencas que no poseen informacin de
precipitacin y/o de caudales requieren la aplicacin de tcnicas de regionalizacin como
las que se han venido presentando en las restantes metodologas. Tanto el caudal pivote
como el Grdex de precipitacin deben extrapolarse a stas.
6.1.5. Resultados
Ya que se han presentado las bases tericas en los cuales se fundamenta el mtodo,
pueden estimarse los parmetros necesarios para su aplicacin. Los resultados
obtenidos provienen del anlisis y procesamiento de informacin de precipitacin y caudal
suministrada al proyecto, as como de la implementacin del software HidroSIG
desarrollado en el Posgrado en Aprovechamiento de Recursos Hidrulicos de la
Universidad Nacional de Colombia, Sede Medelln.
6.1.5.1 Grdex de Precipitacin
De acuerdo con la metodologa expuesta en las secciones anteriores, se procedi a
determinar el Grdex de Precipitacin de 39 estaciones que ya han sido filtradas en los
anlisis previos de homogeneidad. El proceso consisti en determinar inicialmente, la
bondad de ajuste que posee la distribucin tipo Gumbel sobre la serie de mximos
anuales de cada una de las estaciones seleccionadas dentro de la regin de estudio. La
Figura 6.1 muestra la localizacin espacial de stas estaciones.
Una vez evaluada la bondad de ajuste, la cual se hizo a partir de las pruebas estadsticas
Smirnov Kolmogorov y Cramer Von Mises e implementando el software MPH desarrollado
en el Posgrado en Aprovechamiento de Recursos Hidrulicos de la Universidad Nacional,
se estimaron los parmetros de la distribucin por mtodos como el de mxima
verosimilitud y de los momentos insesgados.
CONVENIO 366 / 2003
6-3

Figura 6.1
Estaciones de precipitacin utilizadas en la metodologa gradex
De las 39 estaciones de precipitacin seleccionadas para aplicar la metodologa Grdex,

slo 2 (estaciones El Buey y Ro Abajo) no fueron ajustables por medio de una
distribucin Gumbel, por lo cual fueron desechadas del anlisis.
En la Tabla 6.1 se presentan los resultados obtenidos despus de someter las series de
precipitacin mxima anual a las pruebas de ajuste en el MPH.
Tabla 6.1 Parmetros de ajuste de series de precipitacin mxima anual
PARMETROS DE AJUSTE
MOMENTOS
ESTACIN
NOMBRE
MXIMA VEROSIMILITUD
LOCALIZACION ESCALA LOCALIZACION ESCALA
2308021
FE_LA
50.49
9.11
50.58
8.54
2308022
SEVERA_LA
49.52
9.84
49.18
10.82
2308024
VASCONIA
59.03
9.45
58.88
9.86
2308026
MOSCA_LA
51.36
13.32
50.70
18.40
2308027
RIONEGRO_LA_MACARE
48.48
9.33
48.66
8.54
2701034
MAZO
51.80
9.23
51.69
9.21
2701035
CHORRILLOS
51.12
9.51
51.14
9.46
2701036
CALDAS
57.42
8.77
56.76
11.30
2701037
FABRICATO
55.67
13.54
55.64
13.34
CONVENIO 366 / 2003
6-4

PARMETROS DE AJUSTE
MOMENTOS
MXIMA VEROSIMILITUD
2701038
SAN_ANTONIO_DE_PRA
55.50
10.27
55.63
2701045
VILLAHERMOSA_PLANT
47.07
7.98
47.10
10.09
7.85
2701046
SAN_CRISTOBAL
52.18
10.57
51.69
12.31
2701047
MEDELLIN_MIGUEL_AG
44.27
7.23
44.11
7.65
2701051
BOTON_EL
47.23
7.63
47.25
7.57
2701056
GIRARDOTA
46.58
8.84
46.21
10.03
2701057
BARBOSA
59.84
13.64
59.64
14.31
2701062
BOQUERON
48.04
11.92
47.45
14.24
2701065
DON_MATIAS
57.49
11.50
58.67
8.58
2701066
GABINO
73.57
13.88
72.70
17.74
12.29
2701076
NIQUIA
53.79
13.90
54.16
2701078
BOQUERON
71.18
16.05
71.40
15.38
2701080
MESETA_LA_SN_PEDRO
35.93
20.15
36.25
19.81
2701081
STA_HELENA
57.42
13.36
56.96
14.82
2701082
CUCHILLA_LA
52.18
18.31
53.97
13.46
2701093
AYURA
48.42
8.62
48.57
8.14
2701111
ASTILLEROS
60.73
16.88
61.59
13.63
25.40
2701112
IGUANA_LA
69.33
31.51
71.30
2701114
CUCARACHO
53.56
7.85
53.18
9.19
2701115
ASTILLERO
53.30
18.06
55.04
12.77
2701122
CONVENTO_EL
53.85
7.96
53.72
8.41
2701507
APTO_OLAYA_HERRERA
47.35
10.67
48.13
8.81
2701509
TULIO_OSPINA
52.31
10.81
51.40
15.64
2701512
PIEDRAS_BLANCAS
43.15
26.34
45.40
20.75
2701515
PROGRESO_HDA_EL
51.72
15.37
51.43
16.25
2701526
SALADA_LA
57.62
17.85
57.67
18.68
2701106
PLANTA_MANANTIALES
49.29
8.75
48.93
9.44
2308023
PALMAS_LAS
53.44
12.44
53.79
11.48
Una vez estimado el parmetro (gradex de precipitacin) en cada una de las estaciones
de precipitacin (parmetro de escala de la distribucin Gumbel), se aborda la idea de
regionalizar el parmetro de forma que puedan extrapolarse los resultados obtenidos a
cuencas sin instrumentacin. Para tal fin, se han incorporado puntos de apoyo extrados
del trabajo Adaptacin del mtodo Grdex en cuencas Antioqueas de la Universidad
Nacional de Colombia (Figura 6.2), que permiten densificar la nube de puntos dentro de la
regin de estudio y darle mayor credibilidad a las tcnicas de interpolacin que se han
efectuado para obtener el mapa de grdex que se presenta en la Figura 6.3. Vale la pena
mencionar que los resultados posteriores a la interpolacin, han sido evaluados teniendo
en cuenta el comportamiento regional que posee el parmetro .
CONVENIO 366 / 2003
6-5

Figura 6.2
Estaciones auxiliares de precipitacin utilizadas en la metodologa grdex
Figura 6.3
Mapa de Grdex de precipitacin
CONVENIO 366 / 2003
6-6

A partir del mapa presentado en la figura anterior, puede determinarse el valor medio del
parmetro p al interior de cada una de las cuencas definidas por las 21 estaciones de
monitoreo. De esta forma, la variabilidad espacial de p es tenida en cuenta en nico valor
con el cual puede agregarse el grdex de caudal correspondiente a la cuenca. En la
Tabla 6.2 se presenta el valor medio del gradex de precipitacin en las cuencas definidas
por las estaciones de monitoreo.
Tabla 6.2 Gradex promedio de precipitacin
ESTACIN
GRADEX PRECIPITACIN
PROMEDIO (mm)
Alto San Miguel

La Primavera
Ancn Sur
Antes San Fernando
Despus San Fernando
Puente Guayaquil
Aula Ambiental
Puente Acevedo
Puente Machado
Copacabana
Metromezclas
Puente Girardota
Parque de Las Aguas
Hatillo
Papelsa
Popalito
Pradera
Eade
Puente Gabino
Quebrada Doa Mara
Quebrada La Garca
16.23
15.73
11.73
11.57
11.21
11.00
11.43
11.27
12.11
11.93
11.85
11.72
11.56
11.54
11.60
11.79
11.90
11.98
11.98
12.37
15.42
6.1.5.2 Gradex de caudal

Como se present en el numeral 6.1.2, el exceso de escorrenta a la salida de una cuenca
puede estimarse a partir del rea de la misma, de la pendiente de la distribucin de
probabilidad de los mximos anuales de precipitacin, esto es, el Grdex de precipitacin,
y del tiempo de concentracin al interior de la cuenca. Teniendo en cuenta la resolucin
diaria de la informacin utilizada en la aplicacin del mtodo Grdex, debe suponerse que
el tiempo de concentracin en toda la regin de estudio es del mismo orden. Este
supuesto puede ocasionar tanto subestimacin como sobreestimacin de caudales, segn
el tamao y caractersticas morfomtricas de las subcuencas que se estudien.
Ahora bien, ya que se cuenta con valores medios del gradex de precipitacin en las
cuencas de inters, sobre el eje del ro Medelln y algunos de sus afluentes, se ha
estimado el correspondiente gradex de caudal, q, empleando la expresin 6.8, donde A
CONVENIO 366 / 2003
6-7

(Km2) corresponde al rea de la cuenca, p (mm) al valor medio del gradex de

precipitacin y tc= 24 horas.
q =
A p
3.6tc
(6.8)
En Tabla 6.3 se presenta el valor medio de ste parmetro en la regin de estudio.

Tabla 6.3 Gradex promedio de caudal
ESTACIN
GRADEX DE CAUDAL
3
(m /s)
Alto San Miguel

La Primavera
Ancn Sur
Antes San Fernando
Despus San Fernando
Puente Guayaquil
Aula Ambiental
Puente Acevedo
Puente Machado
Copacabana
Metromezclas
Puente Girardota
Parque de Las Aguas
Hatillo
Papelsa
Popalito
Pradera
Eade
Puente Gabino
Quebrada Doa Mara
Quebrada La Garca
2.8
9.1
18.1
31.7
35.5
42.5
62.8
66.6
89.6
98.2
101.3
106.2
112.5
115.2
138.3
151.5
160.8
167.5
276.9
10.8
14.7
6.1.5.3 Caudal pivote

Como ltimo parmetro necesario para determinar avenidas mximas se encuentra el
caudal pivote o de saturacin Qtc (To), el cual, al igual que la informacin de precipitacin
empleada, corresponde a un nivel promedio diario. El caudal pivote utilizado en la
aplicacin del mtodo, corresponde a la media de los caudales mximos anuales ya que
es una buena aproximacin al caudal correspondiente a un perodo de retorno de 2.33
aos o caudal a banca llena (Smith y Vlez, 199)4 como es conocido en la literatura.
Al igual que los dems parmetros presentados, es necesario proponer una expresin
regional para Q2.33 que permita obtener caudales mximos para diferentes perodos de
retorno.
CONVENIO 366 / 2003
6-8

Para regionalizar el caudal pivote se obtuvo la media de los caudales mximos anuales en
las estaciones de caudal localizadas sobre el eje del ro Medelln, excluyendo aquellas
con series demasiado cortas. Posteriormente, se ajust este parmetro con la
correspondiente rea de drenaje definida por cada estacin. De esta forma se puede
obtener el parmetro Q2.33 en cuencas instrumentadas o no instrumentadas a partir de la
expresin (6.9):
Q2.33 = 0.1308 A1.005
(6.9)
6.1.5.4 Avenidas mximas

Una vez determinados los parmetros requeridos en la aplicacin del mtodo Grdex, se
procedi a determinar los caudales mximos diarios en las estaciones de monitoreo
preestablecidas en el proyecto. En la Tabla 6.4 se presentan los caudales para diferentes
perodos de retorno.
Tabla 6.4 Caudales mximos para diferentes perodos de retorno
3
Caudal mximo (m /s)

ESTACIN
Perodo de retorno (aos)

2.33
10
25
50
100
Alto San Miguel
5.39
7.53
9.47
12.04
13.98
15.92
La Primavera
17.54
24.47
30.77
39.09
45.38
51.68
Ancn Sur
45.77
59.58
72.11
88.68
101.22
113.75
Antes San Fernando
80.36
104.61
126.61
155.70
177.71
199.71
Despus San Fernando
92.54
119.68
144.32
176.89
201.52
226.16
Puente Guayaquil
112.16
144.58
174.02
212.93
242.36
271.80
Aula Ambiental
158.20
206.12
249.62
307.13
350.62
394.12
Puente Acevedo
170.09
220.98
267.17
328.23
374.43
420.62
Puente Machado
211.73
280.16
342.27
424.38
486.50
548.61
Copacabana
234.82
309.77
377.80
467.74
535.78
603.81
Metromezclas
243.72
321.05
391.25
484.05
554.25
624.45
Puente Girardota
258.06
339.15
412.77
510.09
583.70
657.32
Parque de Las Aguas
276.65
362.54
440.52
543.59
621.57
699.54
Hatillo
283.72
371.71
451.59
557.18
637.06
716.94
Papelsa
337.42
443.05
538.93
665.69
761.58
857.47
Popalito
362.92
478.59
583.60
722.41
827.42
932.42
Pradera
381.16
503.93
615.37
762.70
874.14
985.59
Eade
394.20
522.12
638.25
791.76
907.89
1024.01
Puente Gabino
643.85
855.28
1047.21
1300.92
1492.85
1684.78
Quebrada Doa Mara
26.18
34.40
41.86
51.72
59.18
66.64
Quebrada La Garca
28.66
39.90
50.10
63.59
73.79
83.99
CONVENIO 366 / 2003
6-9

Las metodologas empleadas en la estimacin de caudales mximos no incluyen los

efectos de las principales descargas de aguas residuales que se presentan a lo largo del
ro y que fueron analizadas en el Captulo 4, ya que no representan una cantidad
importante frente a los caudales que naturalmente se presentan en las zonas sur y centro
de la cuenca, incluso para los menores perodos de retorno. Sin embargo, en la zona
norte, los caudales que aporta la central de la Tasajera al ro constituyen porcentajes
superiores al 15 % respecto a los valores mximos diarios estimados a partir del mtodo
Gradex y para perodos de retorno bajos. A pesar de stos rdenes de magnitud no
puede tipificarse este tipo de descarga razn por la cual se ha aplicado las metodologas
para mximos bajo el supuesto del no funcionamiento de la central.
6.2 REGIONALIZACIN DE CARACTERSTICAS MEDIAS
De igual forma que la regionalizacin de las caractersticas medias para caudales
mnimos, pueden estimarse caudales mximos para asociados a cualquier periodo de
retorno, relacionando los parmetros estadsticos de los caudales extremos con
caractersticas geomorfolgicas, climticas y topogrficas de cada subcuenca.
La ecuacin de Ven te Chow, permite obtener caudales mximos asociados a cualquier
perodo de retorno as:
QTr = Q max + K Q max
(6.10)
Donde QTr es el caudal mximo para un periodo de retorno Tr; Q max es la media
estimada de los caudales mximos instantneos anuales; Q max es la desviacin estndar
de los caudales mximos instantneos anuales y K es el factor de frecuencia que
depende de la distribucin, del periodo de retorno y del tipo de caudal extremo.
El factor de frecuencia K de caudales mximos utilizando una distribucin de dos
parmetros, en este caso la distribucin Gumbel (Chow, 1953), se puede calcular as:
K=
6
Tr

0.5772 + ln ln

Tr 1
(6.11)
Utilizando la metodologa propuesta por UNALMED-UPME (2000), la media y la

desviacin estndar de los caudales mnimos esta dada por:
donde
max = K (P E ) A
(6.12)
max = K (P E ) A
(6.13)
K y K son el factor de frecuencia para la media y la desviacin estndar,
respectivamente; P, la precipitacin media sobre la cuenca; E, la evapotranspiracin; A, el

CONVENIO 366 / 2003
6-10

rea de la cuenca y , , y son constantes asociadas a la media y a la

desviacin estndar.
Los valores de dichas constantes se estimaron mediante
procedimientos de regresin mltiple, utilizando como variables independientes los
parmetros climticos (precipitacin y evapotranspiracin) y geomorfolgicos de la cuenca
(el rea de drenaje, en este caso) y como variables dependientes los valores de max y
max calculados para cada estacin de caudal localizada al interior de la cuenca del ro
Medelln.
Las estaciones de medicin con informacin de caudales mximos instantneos
mensuales en la cuenca del ro Medelln utilizadas para elaborar las regresiones son:
Tabla 6.5 Estaciones con informacin de caudales mximos instantneos sobre el ro
Medelln
CODIGO
NOMBRE
2701733
SALADA_LA_RMS_11
2701727
ANCON_SUR_RMS_17
2701734
MACHADO_RMS.12
2701781
ANCON_NORTE_RMS_20
2701803
GIRARDOTA
2701738
HATILLO_EL_RMS.13
2701735
YARUMITO_RMS_14
2701736
GABINO_RMS.15
Para el clculo de max y max se eligieron los periodos que se consideren tienen
informacin confiable, adems se consideraron los periodos en las estaciones de la parte
baja de la cuenca que se encuentran influenciados por la descarga de la central La
Tasajera, pues se considera que en eventos extremos mximos esta planta de
generacin no realiza descargas al ro.
La informacin con la que se realizaron las regresiones mltiples se muestra a
continuacin:
Tabla 6.6 Informacin empleada para la obtencin de las ecuaciones de regionalizacin
de caudales mximos
max
max
P
E
A
(mm/ao) (mm/ao) (km2)
ESTACIN
PERIODO
SALADA_LA_RMS_11
ANCON_SUR_RMS_17
MACHADO_RMS.12
ANCON_NORTE_RMS_20
GIRARDOTA
HATILLO_EL_RMS.13
YARUMITO_RMS_14
81-86
80-95
80-83
89-97
95-03
81-03
81-03
8.5422
4.5819
69.4036 47.2235
219.0800 62.4881
418.1125 98.8338
334.0667 55.3267
271.5078 103.8690
321.7652 72.1432
2186.89
2159.5
1905.58
1872.67
1863.86
1858.55
1928.95
1019.6
1016.6
981.59
976.3
975.09
974.64
983.47
43.77
121.93
639.24
738.49
783.42
862.79
1099.1
GABINO__RMS.15
72-03
452.6875 124.7329
2084.30
996.37
2541.4
CONVENIO 366 / 2003
6-11

Las variables P y E representan la precipitacin y la evaporacin media en cada cuenca

respectivamente. El valor de la evaporacin (E) es el ledo del mapa de evaporacin de
Cenicaf.
Una vez obtenida la media y desviacin estndar de los caudales mximos en cada
estacin, se realizan las regresiones mltiples, que tuvieron como resultado las siguientes
ecuaciones:
max = 0.38546 A 0.97650

max = 0.64579 A0.70852
R = 0.946395
(6.14)
R = 0.880422
(6.15)
No se encontr relacin de estos parmetros con la lmina de escorrenta estimada (P-E).

El coeficiente R presentado es el valor del coeficiente de correlacin mltiple.
Los caudales mximos para diferentes perodos de retorno, estimados por el Mtodo de
Regionalizacin de caractersticas medias, en los puntos de aforo y muestreo de calidad
de aguas en el ro Medelln se presentan en la Tabla 6.7.
Tabla 6.7 Caudales mximos estimados en los sitios de aforo, segn el mtodo de
Regionalizacin de Caractersticas Medias.
CAUDAL MXIMO (m3/s)
PERODO DE RETORNO (Aos)

2.33 5
10
25
50
100
ESTACIN
AREA
(km2)
max
max
Alto San Miguel

La Primavera
Ancn Sur
Antes de San Fernando
Despus de San Fernando
Puente Guayaquil
Aula Ambiental
Puente Acevedo
Puente Machado
Copacabana
Metromezclas
Puente Girardota
Parque de Las Aguas
Hatillo
Papelsa
Popalito
Pradera
Eade
Puente Gabino
Quebrada Doa Maria
14.91
49.88
133.20
237.08
273.94
333.54
474.38
510.91
639.35
710.85
738.43
782.96
840.77
862.79
1030.37
1110.18
1167.36
1208.27
2541.40
75.17
5.39
17.54
45.77
80.36
92.54
112.16
158.20
170.09
211.73
234.82
243.72
258.06
276.65
283.72
337.42
362.92
381.16
394.20
814.77
26.18
4.38
10.31
20.67
31.10
34.45
39.61
50.84
53.58
62.81
67.71
69.56
72.50
76.26
77.67
88.08
92.86
96.22
98.60
166.98
13.78
5.4
17.6
45.8
80.4
92.6
112.2
158.3
170.1
211.8
234.9
243.8
258.1
276.7
283.8
337.5
363.0
381.3
394.3
814.9
26.2
8.5
25.0
60.6
102.7
117.3
140.7
194.8
208.6
256.9
283.5
293.8
310.2
331.5
339.6
400.8
429.7
450.4
465.1
934.9
36.1
11.1
31.0
72.7
120.9
137.5
163.8
224.5
240.0
293.7
323.2
334.5
352.6
376.1
385.0
452.3
484.1
506.7
522.8
1032.6
44.2
14.3
38.6
88.0
143.9
163.0
193.1
262.1
279.6
340.1
373.2
385.9
406.2
432.5
442.5
517.4
552.7
577.8
595.7
1156.0
54.3
16.7
44.3
99.4
161.0
181.9
214.8
290.0
309.0
374.5
410.3
424.0
446.0
474.3
485.1
565.7
603.6
630.6
649.8
1247.6
61.9
19.1
49.9
110.6
177.9
200.6
236.4
317.7
338.2
408.7
447.2
461.9
485.5
515.8
527.3
613.7
654.2
683.0
703.5
1338.5
69.4
Quebrada La Garca
82.47
28.66
14.72
28.7
39.2
47.9
58.7
66.8
74.8
CONVENIO 366 / 2003
6-12

Los caudales mximos calculados mediante esta metodologa, se consideran que no

estn afectados notablemente por el influjo de las descargas de aguas residuales de los
municipios del Valle del Aburr, pues en eventos extremos mximos el aporte de estas
aguas es mnimo respecto al caudal natural en el ro.
6.3 ESTIMACIN DE CAUDALES MXIMOS PARA DIFERENTES PERODOS DE
RETORNO MEDIANTE MODELACIN HIDROLGICA UTILIZANDO HEC-HMS
Para la estimacin de la creciente de diseo correspondiente a tormentas estimadas para
diferentes periodos de retorno se emplearon metodologas de clculo semi - agregado de
crecientes; en las cuales la cuenca de estudio puede ser subdividida en una serie de
subcuencas y canales interconectados entre si, permitiendo considerar de una forma un
poco ms detallada la variabilidad espacio temporal de las caractersticas
geomorfolgicas y climticas de la cuenca de estudio. Para tal fin, se utiliza el software
HEC-HMS del cuerpo de ingenieros de la armada de los Estados Unidos (U.S. Army
Corps of Engineers, 2001), que integra varias metodologas de clculo de crecientes para
el diseo hidrolgico.
El modelo implementado para la cuenca del Ro Medelln permite realizar el trnsito de la
creciente en toda la cuenca mediante la representacin del sistema de drenaje como la
unin de las subcuencas principales que drenan al Ro Medelln y los tramos de canales
existentes entre ellas; permitiendo obtener los caudales mximos en cada subcuenca, en
los canales y en diferentes puntos de control dentro de la red de drenaje.
Para la determinacin de los parmetros geomorfolgicos requeridos por el modelo se
empleo la cartografa disponible y el modelo de elevacin digital con una resolucin de
pixel de 90 x 90 m de la NASA (Disponible de forma gratuita en
http://ham.eit.com/SRTMEarthThreeArcSecond). Para el procesamiento de dicha
informacin en SIG se empleo el software IDRISI Kilimanjaro versin 14. En la Figura 6.4
se muestra un esquema discretizado para las principales subcuencas que drenan al Ro
Medelln y el modelo de elevacin digital para toda la cuenca. En el Anexo B se presentan
los parmetros obtenidos para cada cuenca.
Para la modelacin en HEC-HMS fueron consideradas otras cuencas ms pequeas
teniendo en cuenta el siguiente criterio de agregacin: a) Zonas altas y parte media (zona
urbana) reas pequeas entre 0 y 3 Km2. b) Zonas Bajas: reas entre 3 y 10 Km2 En el
anexo C se presenta la lista de subcuencas agregadas y en el Anexo D se presenta los
puntos de control seleccionados.
Para la estimacin de caudales mximos segn los diferentes componentes empleados
para la modelacin deben definirse los mtodos y parmetros a emplear en los siguientes
procesos:
Cuencas
Mtodos de transformacin
Modelos de Perdida
Mtodo para estimar flujo base
CONVENIO 366 / 2003
6-13

Canales
Mtodo de transito de crecientes
Modelo meteorolgico
En este se definen el mtodo deseado para estimar la lluvia total de diseo.
Con base a la informacin disponible se presenta a continuacin una descripcin de los
diferentes mtodos de transformacin, perdidas, flujo base y transito hidrolgico
empleados en la modelacin de los caudales mximos en la cuenca del Ro Medelln.
Figura 6.4
6.3.1.
Cuencas principales que drenan al ro Medelln.
Modelos de prdida
Para simular las perdidas de la precipitacin, esto es la cantidad de agua que se pierde
por intercepcin en la vegetacin, almacenamiento en depsitos del suelo y por
infiltracin, el HEC-HMS permite emplear diferentes mtodos tales como el nmero de
curva de la S.C.S (Soil Conservation Service), el mtodo del ndice (U.S. Army Corps of
Engineers, 2001; Chow V.T., 1994), mtodo de Green Ampt (1991) y el SMA (.Soil
Moisture Accounting Loss Method), entre otros.
El primero fue desarrollado por el Servicio de Conservacin de suelos de los Estados
Unidos (SCS) y relaciona la lluvia y la retencin en el suelo con la lluvia efectiva o
escorrenta directa. El SCS defini un nmero adimensional de curva CN, luego de
realizar un estudio sobre muchas cuencas en los Estados Unidos, este valor depende del
tipo y el uso de suelo y de las condiciones de humedad antecedente. Por otro lado, el
mtodo del ndice utiliza la prdida inicial y una tasa constante de prdidas en la cuenca
luego de que ha sucedido la lluvia. Estos valores son obtenidos a partir de los
hietogramas de precipitacin y de calibraciones previas.
CONVENIO 366 / 2003
6-14

6.3.1.1 Modelo Green-Ampt

Mediante el anlisis del flujo del agua a travs de una columna de suelo de rea
transversal unitaria A y de longitud L Green-Ampt (1911) combina la ley de Darcy para
obtener una ecuacin de para estimar la tasa de infiltracin f(t) como una funcin de
variables como la succin del suelo que a la vez es una funcin del contenido de
humedad del suelo, el dficit de humedad estimado como la diferencia entre la porosidad
del suelo y el contenido de humedad en un tiempo dado y la conductividad hidrulica.
Diferentes autores han realizado estudios para determinar algunas de las variables ms
importantes del modelo Greem - Ampt tales como Brooks y Corey (1964), Brakensiek,
Engleman y Rawls (1981) encontrando por ejemplo valores tpicos de porosidad,
conductividad hidrulica y cabeza de succin de algunos suelos Kensiek y Miller (1983),
como se muestra en la Tabla 6.8.
Tabla 6.8 Parmetros para el mtodo de Green - Ampt para algunos tipos de suelo
Fuente HEC-HMS
Clase de Textura
Arena
Marga arenosa
Arena Margosa
Marga
Marga limosa
Marga arenoarcillosa
Marga arcillosa
Marga limoarcillosa
Arcilla arenosa
Arcilla limosa
Arcilla
Porosidad
(cu.cm/cm)
0.437
0.437
0.453
0.463
0.501
0.398
0.464
0.471
0.430
0.479
0.475
Conductividad
Hidrulica
Saturada
(cm/hr)
21.00
6.11
2.59
1.32
0.68
0.43
0.23
0.15
0.12
0.09
0.06
Cabeza de
Succin de
Frente hmedo
(cm)
10.6
14.2
22.2
31.5
40.4
44.9
44.6
58.1
63.6
64.7
71.4
6.3.1.2 Modelo SMA (.Soil Moisture Accounting Loss Method)

El modelo SMA estima el exceso de precipitacin considerando la intercepcin de la lluvia
la cual corresponde a la parte que es retenida por los rboles, arbustos, hierbas y que no
alcanzan el perfil del suelo, el almacenamiento superficial debido a las depresiones del
terreno la cual ocurre cuando se excede la capacidad de infiltracin, la Infiltracin, el
proceso de percolacin y la evapotranspiracin mediante un modelo de almacenamiento
de cinco capas en cada una de las cuales se consideran las entradas y salidas de
acuerdo al tipo de proceso involucrado (balance hdrico). Algunos autores como Dunne
(1978), Caulder (1980), Zinke (1967), Crawford (1966) han realizado estudios para tratar
de determinar el porcentaje de intercepcin para diferentes especies vegetales y otros
autores como Smith (2003), Flemmint et al (2004) han estudiado la variabilidad espacio
temporal de las variables del modelo SMA mediante utilizacin de herramientas SIG.
CONVENIO 366 / 2003
6-15

A pesar que los modelos SMA y Greem Ampt presentan una sustentacin fsica un poco
mas detallada respecto a los modelos del nmero de curva de la SCS y el mtodo del
ndice , requieren tambin una cantidad de parmetros mayor que muchas veces no
estn disponibles debiendo ser obtenidos a partir de valores calibrados en sitios con
caractersticas geomorfolgicas y climticas diferentes a las nuestras.
6.3.1.3 Mtodo del ndice
Teniendo en cuanta lo anterior para la cuenca del Ro Medelln se empleo el mtodo del
ndice (Initial/Constant) basado en los valores presentados por Montoya, J.J. y Patio,
J.E., (2001) en su trabajo Estudio Hidrometeorolgico y Geomorfolgico de diez
quebradas del Valle de Aburr. En este trabajo se presentan relaciones del coeficiente de
escorrenta para varias cuencas, y se utilizan las estimaciones de este coeficiente con los
datos de campo para calcular la tasa de infiltracin (ndice 10 mm/h). En cuanto a la
perdida inicial, esta se supone igual a la obtenida en el modelo del nmero de curva
presentado por el Servicio de Conservacin de Suelos de Estados Unidos (S.C.S). Para
estudios posteriores se recomienda la estimacin de dichos parmetros de forma
espacialmente distribuida a partir de los datos de formaciones superficiales, coberturas
vegetales y los valores presentados por la S.C.S como se muestra en la Tabla 6.9.
Tabla 6.9 Valores de tasas de infiltracin segn grupos de la S.C.S. SCS, 1986; Skaggs
and Khaleel, 1982).
Grupo
de suelo
A
B
C
Descripcin
Arenas profundas, loess profundos, limos
agregados.
Loess superficial, Marga arenosa.
Arcillas
Margosas,
Margas
arenosas
superficiales, Suelos con bajo contenido de
materia orgnica y suelos con alto contenido
de arcillas.
Suelos que se ablandan significativamente
cuando estn hmedos, Arcillas de alta
plasticidad y ciertos suelos salinos 0.00-0.05
Rangos de
Tasas de prdida
[mm/h]
7.62 11.43
3.8 7.62
1.27 3.81
0 1.27
6.3.1.4 Mtodo Initial/Constant (HEC-HMS)

Los parmetros requeridos por el mtodo Initial/Constant del HEC-HMS son la tasa de
infiltracin constante, la prdida inicial y el coeficiente de impermeabilidad, el cual
corresponde al porcentaje total del rea de cada subcuenca que es totalmente
impermeable.
En este modelo el flujo de escorrenta es generado nicamente cuando la cantidad de
precipitacin supera las prdidas iniciales de precipitacin, y adems, la tasa constante de
CONVENIO 366 / 2003
6-16

infiltracin es inferior a la tasa de precipitacin promedia en el rea de la subcuenca. El

modelo matemtico se presenta a continuacin.
0
si pi < I a
pet = pt f c si pi > I a y pt > f c
si pi > I a y pt < f c
0
(6.16)
Donde,
Pet
Pt
Pi
Fc
Ia
:Excesos de precipitacin
:Precipitacin promedia en el rea de la subcuenca
:Precipitacin promedia para el intervalo i
:Tasa de infiltracin
:Perdida Inicial de precipitacin
Para la determinacin del coeficiente de impermeabilidad se emplea la informacin

cartogrfica digital, implementando las herramientas del SIG se obtienen las reas
urbanas para todas las subcuencas de la zona de estudio. Los valores obtenidos fueron
multiplicados por los coeficientes de impermeabilidad para techos o para zonas comercial
o industrial (segn la zona de ubicacin) propuestos por Obras Pblicas como se muestra
en la Tabla 6.10.
Tabla 6.10 Coeficientes de impermeabilidad por zona urbana Segn obras publicas.
Impermeabilidad
Caracterstica
[I]
1.Techos
0.9
2.Comercial o industrial
0.9
3.Residencial con casas contiguas, predominio de zonas impermeables.
0.75
4.Residencial mutifamiliares
impermeables.
0.75
con
bloques
contiguos
zonas
5.Residencial unifamiliares con casas contiguas, predominio de jardines.
0.55
6.Residencial con casas con jardines
0.45
7.Residencial con predominio de zonas verdes, cementerios tipo jardn.
0.30
8.Laderas desprotegidas con vegetacin
0.6
9.Laderas protegidas con vegetacin.
0.3
A continuacin se presenta una descripcin del mtodo de la S.C.S y la metodologa

empleada para la determinacin del nmero de curva (CN).
CONVENIO 366 / 2003
6-17

6.3.1.5 Mtodo de la Soil Conservation Service (S.C.S)

El enfoque del Numero de Curva (CN) es el mtodo ms ampliamente usado en los
Estados Unidos. Fue desarrollado en los aos 70 por el Servicio de Conservacin de
Suelos de los Estados Unidos (SCS) para clasificar la escorrenta potencial para
diferentes combinaciones de suelo y cobertura. El mtodo relaciona la escorrenta
superficial (R) con la precipitacin total (P) y con la capacidad de almacenamiento (Smax)
mediante la siguiente ecuacin emprica:
R=
( P 0.2S max ) 2
, Para P>0.2 Smx
( P + 0.8S max )
(6.17)
donde: La precipitacin (en pulgadas) es medida, y Smax es calculada como una funcin
del numero de curva (CN) donde:
S max =
1000
10
CN
(6.18)
El intervalo del numero de curva va desde 0 a 100 y es funcin de la textura del suelo, la
cobertura del suelo y los antecedentes de humedad del suelo. Entre menor sea el
nmero de curva, ms bajo es el potencial de escorrenta, las condiciones del suelo estn
asignadas en cuatro grupos hidrolgicos con base a la textura y estructura del suelo como
se presento en la Tabla 6.9. Con base a la humedad antecedente se definen tres grupo
AMC, como se muestra en la Tabla 6.11.
Tabla 6.11 Clasificacin de clases antecedentes de humedad (AMC) para el mtodo de la
S.CS. Fuente:Soil Conservation Service.
Grupo
AMC
Ia
IIb
IIIc
Lluvia antecedente total

de cinco das (pulg)
Estacin Inactiva
Estacin de Crecimiento
Menor que 0.5
0.5 a 1.1
Sobre 1.1
Menor que 1.4

1.4 a 2.1
Sobre 2.1
Los antecedentes de humedad afectan la escorrenta potencial. Entre ms hmedas sean

las condiciones antecedentes ms alto es el potencial de escorrenta.
Debido a que los datos de capacidad de infiltracin generalmente no estn disponibles
para todos los suelos, la textura y usos del suelo es usada para asignar a cada tipo de
suelo un grupo valor de CN. En la Tabla 6.12 se presentan los valores de CN para
diferentes usos y tratamientos del suelo.
a
b
c
Suelos secos
Suelos intermedios
Suelos humedos
CONVENIO 366 / 2003
6-18

Tabla 6.12 Valores de CN para AMC II. Para pendiente del terreno mayores que el
1%.Fuente: Vlez 2000
Uso de la tierra y
cobertura
Tratamiento
del suelo
Tipo de suelo
Sin cultivo
Surcos rectos
A
77
Cultivos en surcos
Surcos rectos
Contorneo
Terrazas
72
70
66
81
79
74
88
84
80
91
88
82
Cereales
Surcos rectos
Contorneo
Terrazas
Terrazas
65
63
61
59
76
74
72
70
84
82
79
78
88
85
82
81
Surcos rectos
Contorneo
Terrazas
66
64
63
77
75
73
85
83
80
89
85
83
---------Contorneo
68
47
79
67
86
81
89
88
56
46
36
26
15
75
68
60
52
44
86
78
70
62
54
91
84
77
69
61
72
74
82
84
87
90
89
92
Leguminosas o
praderas con rotacin
Pastizales
Bosques
naturales
Muy ralo
Ralo
Normal
Espeso
Muy espeso
-------------------
-------------------
Caminos
De terracera
De superficie dura
B
86
C
91
D
94
Estimacin del nmero de curva (CN)

Los nmeros de curva promedio para cada una de las 137 subcuencas se obtuvieron
mediante SIG a partir de las imgenes en formato raster de cobertura vegetal, mapa de
pendientes y tipos de suelo para la condicin de humedad antecedente B (segn la Tabla
6.12). Por ejemplo, para las cuencas altas con bajo componente urbano como La Gmez
parte alta, se utiliz un nmero de curva relativamente bajo debido a que esta presenta
poca densidad poblacional y por ende baja componente urbana (viviendas, vas, placas
polideportivas, etc.). Para las cuencas en la zona urbana de los diferentes municipios las
cuales presentan alta intervencin como por ejemplo quebradas como la Iguana, La
Hueso, La Garca, La Ayura, La Poblada etc, se seleccionaron nmeros de curva
CONVENIO 366 / 2003
6-19

relativamente altos lo que indica una mayor escorrenta superficial directa. En la Figura
6.5 se presentan los valores obtenidos.
Figura 6.5
Valores medios de CN para las cuencas tributarias al ro Medelln
Para la determinacin de las perdidas iniciales se considero que aquellas cuencas menos
antropizadas tienen mayores prdidas iniciales que las cuencas que se encuentran ms
intervenidas.
6.3.2.
Medicin de la Escorrenta Superficial Directa (Mtodos de transformacin)
Mediante esta aplicacin se gener un modelo lluvia escorrenta acoplado para toda la
cuenca, donde se incluyeron las subcuencas tributarias del ro Medelln. A continuacin se
exponen en orden metodolgico los pasos para el clculo de los eventos de crecida en los
puntos seleccionados para la cuenca, como tambin aquellas cuencas que no poseen un
nombre definido.
Como mtodo de transformacin de lluvia se empleo el mtodo de las hidrgrafas
unitarias sintticas, el cual permite construir el hidrograma de escorrenta superficial
directa para una lluvia de duracin y profundidad unitaria para una cuenca sin
informacin. En gran parte de los mtodos existentes para construir hidrgrafas unitarias
sintticas el hidrograma es definido a partir de los componentes principales que describen
la forma de este como son el tiempo de rezago, el tiempo de concentracin, el tiempo
base, el tiempo al pico, etc. Estas variables son calculadas como una funcin de
parmetros morfomtricos de la cuenca tales como el rea de la cuenca, pendiente de la
cuenca, longitud del cauce principal, Cota mxima, cota mnima, distancia al centroide etc;
variando la forma de calculo entre las diferentes metodologas.
Existen otras metodologas que permiten realizar un balance de masas de forma
espacialmente distribuida a nivel de celda o pixel, sin embargo dichas metodologas
CONVENIO 366 / 2003
6-20

requieren de mucha informacin espacial que no esta disponible, debido a la magnitud

total de la cuenca del Ro Medelln y la poca informacin para calibrar y validar dichos
modelos.
Es importante anotar que las subcuencas que drenan al Ro Medelln, pocas de ellas
poseen registros de caudal a resolucin diaria, lo cual impide la calibracin y validacin de
los parmetros de los modelos.
Como mtodos de transformacin el HEC- HMS ofrece diferentes posibilidades dentro de
las cuales, entre otros se encuentran el hidrograma unitario de Clark (1945) el cual
conceptualiza la cuenca como un modelo de embalse lineal en el cual se representa el
impacto agregado en la produccin de la escorrenta superficial directa del
almacenamiento en el suelo, en la superficie y en los canales. Para esto se considera el
tiempo de viaje del agua desde los puntos ms lejanos de la cuenca hasta la salida de la
misma empleando un modelo de canal lineal (Dooge 1959); esto es obtenido a travs del
histograma de rea - tiempo el cual representa al rea contribuyente al flujo a la salida de
la cuenca como una funcin del tiempo.
Otro de los mtodos de transformacin contenidos dentro del HEC-HMS corresponde al
mtodo de onda cinemtica en el cual se conceptualiza la cuenca en un modelo de
elementos finitos por dos superficies planas representando cada una de las vertientes de
la cuenca y un canal central que representa el cauce principal. Para el flujo en ladera y en
el canal central el modelo emplea las ecuaciones bsicas de flujo en canales, la ecuacin
de momentum y la ecuacin de continuidad para generar las ecuaciones iterativas del
modelo.
Las otras opciones incluidas en el HEC-HMS permiten la utilizacin de los muy conocidos
y empleados en nuestro medio hidrogramas unitarios de Snyder y S.C.S, as como la
posibilidad emplear un hidrograma unitario deseado por el usuario (Posibilidad de
emplear hidrogramas unitarios obtenidos a partir de datos simultneos de lluvia y
caudales).
Para la cuenca del Ro Medelln se emplearon tres de las metodologas de transformacin
mencionadas anteriormente, las cuales se presentan a continuacin de forma mas
detallada.
6.3.2.1 Modelo Hidrograma Unitario de Snyder*
El modelo de Snyder deriva un Hidrograma Unitario a partir de algunas caractersticas
fsicas de la cuenca, para ser aplicado en las cuencas donde no se tienen registros
histricos de caudal. Este modelo considera los siguientes componentes del hidrograma.
Caudal Pico
El modelo propone calcular el caudal pico por milla cuadrada, uP, mediante la siguiente
expresion :
*
Fuente:Velez y Smith (1997)
CONVENIO 366 / 2003
6-21

up = Cp
640
TR
(6.19)
Donde uP es el caudal pico del Hidrograma Unitario por unidad de rea en pie/s.mi , CP
es un coeficiente, T es la duracin de la precipitacin efectiva en horas, TR es el tiempo de
rezago y tS es la duracin de la lluvia efectiva en horas. El coeficiente CP depende de la
topografa de la cuenca variando entre 0.5 y 0.8, por ejemplo, para cuencas de alta
pendiente se recomienda empleara valores cercanos a 0.8 (menor atenuacin del
hidrograma de respuesta). Este coeficiente puede ser obtenido si se tiene informacin
simultanea de precipitacin y caudal de la cuenca, lo cual solo es posible en cuencas
instrumentadas.
Cuando el hidrograma unitario sinttico a construir corresponda a una precipitacin
efectiva cuya duracin coincide con la ecuacin 6.22, el caudal pico del hidrograma
unitario por unidad de rea puede calcularse con la ecuacin 6.16; en caso contrario si la
duracin de la lluvia es diferente puede calcularse como.
up = Cp
640
[TR + (T t s ) / 4]
(6.20)
Una vez obtenido el caudal pico por unidad de rea de la cuenca, el caudal pico total se
obtiene como:
UP = uP AC
(6.21)
En donde AC es el rea de la cuenca en mi.y UP es el caudal pico del Hidrograma Unitario

Sinttico en pie/s/pul.
Tiempo al pico
El tiempo puede calcularse como:
TP = T/2 + TR (6.22)
En donde TP es el tiempo al pico en horas, T es la duracin de la lluvia en horas y TR es el
tiempo de rezago en horas.
Tiempo de rezago
El modelo emplea como definicin de tiempo de rezago (TR), el tiempo comprendido entre
el centroide del hietograma de precipitacin efectiva y el pico del hidrograma de
escorrenta directa correspondiente. Se asume que el rezago es constante para una
cuenca, ya que depende de algunas caractersticas fsicas y no est determinado por el
tipo de lluvia o sus variaciones. El tiempo de rezago se calcula mediante la siguiente
ecuacin de regresin:
TR =
L0c.8 ( S + 1) 0.7
1900 * S c0.5
(6.23)
CONVENIO 366 / 2003
6-22

En donde TR es el tiempo de rezago en horas, Lc es la longitud del canal principal en pies,

S es el factor de retencin o almacenamiento calculado en trminos del nmero de curva
CN y Sc es la pendiente de la cuenca en porcentaje. Para las cuencas en donde se
obtienen tiempos de rezago mayores que el tiempo de concentracin se recomienda
utilizar la siguiente expresin
TR = 0.6TC
(6.24)
Duracin de la lluvia seleccionada por Snyder

Snyder consider lluvias que estuvieran de acuerdo con el tamao de la cuenca,
definiendo para ellas una duracin de 1/5.5 veces el tiempo de rezago de la cuenca, es
decir:
T
ts = R
5 .5
(6.25)
Donde ts es la duracin de la lluvia efectiva en horas.

Duracin de la escorrenta superficial
Para el clculo del tiempo base del Hidrograma Unitario se propone la siguiente ecuacin:
tb = 3 + 3
TR
(6.26)
24
En donde tb es el tiempo base en das con TR en horas. Esta ecuacin da una estimacin
razonable del tiempo base para cuencas grandes, pero produce valores excesivamente
altos para cuencas pequeas. En estas el tiempo base puede calcularse, en forma
aproximada, como 3 a 5 veces el tiempo al pico.
Esquematizacin de la hidrgrafa
El Cuerpo de Ingenieros de los EEUU introdujo despus de los estudios de Snyder dos
ecuaciones adicionales a este modelo, con el objeto de obtener otros cuatro puntos del
Hidrograma Unitario Sinttico de Snyder que facilitan su definicin. Estas ecuaciones
son:
W50 =
770
u 1p.08
W75 =
440
u 1p.08
(6.27)
En donde: uP es el caudal pico por unidad de rea en pie/s mi, W50 es el intervalo de
tiempo en horas, correspondiente al 50% del caudal pico y W75 es el intervalo de tiempo
en horas, correspondiente al 75% del caudal pico.
Como se puede ver en la Figura 6.6 cada intervalo de tiempo se ubica en la curva de tal
forma que su tercera parte quede a la izquierda de la vertical que pasa por el pico del
Hidrograma Unitario y las otras dos terceras partes a la derecha de dicha lnea (Snyder,
1938).
CONVENIO 366 / 2003
6-23

Figura 6.6
Esquematizacin del Hidrograma Unitario Sinttico de Snyder.
6.3.2.2 Modelo Hidrograma Unitario del S.C.S*

El Servicio de Conservacin de Suelos de los Estados Unidos. (Soil Conservation Service,
SCS) desarroll un Hidrograma Unitario Adimensional a partir de una serie de hidrgrafas
observadas, correspondientes a cuencas de muy diversos tamaos y ubicadas en
distintos sitios de los Estados Unidos.
Relacin entre tiempos caractersticos.
El Hidrograma Unitario adimensional curvilneo del S.C.S puede ser representado por un
Hidrograma Unitario Triangular equivalente, con las mismas unidades de tiempo y caudal,
teniendo por consiguiente el mismo porcentaje del volumen en el lado creciente del
Hidrograma como se observa en la Figura 6.7.
Figura 6.7
*
Hidrograma Unitario Adimensional del S.C.S.
Fuent:Velez y Smith (1997)
CONVENIO 366 / 2003
6-24

Caudal pico
El caudal pico en pies cbicos por segundo por pulgada se puede calcular como:
Up =
484 Ac
Tp
(6.28)
En donde UP es el caudal pico correspondiente a un Hidrograma Unitario, AC es el rea de

la cuenca en millas cuadradas y TP es el tiempo al pico en horas.
Obtencin del tiempo de concentracin segn el SCS
El mtodo del SCS utiliza el tiempo de concentracin, el cual define como el tiempo que
demora la escorrenta en llegar desde el punto ms lejano de la cuenca hasta el punto de
aforo, o de manera equivalente, el tiempo que transcurre desde el final de la lluvia efectiva
hasta el punto de inflexin de la rama decreciente del Hidrograma Unitario. Con base en
esto, el SCS propone una relacin promedio entre el tiempo de rezago, TR y el tiempo de
concentracin, TC, como:
TC = 5/3 TR
(6.29)
El tiempo de rezago, TR, definido como el tiempo en horas desde el centroide del
hietograma de la precipitacin efectiva hasta el caudal pico del hidrograma unitario, puede
calcularse como:
L0c.8 (S + 1)
TR =
1900 S c0.5
0.7
(6.30)
En donde LC es la longitud del canal principal en pies, SC es la pendiente promedio de la

cuenca en porcentaje y S es el factor de retencin o almacenamiento en pulgadas.
El factor de retencin o almacenamiento S se obtiene del mtodo de la S.C.S para
perdidas mediante la siguiente expresin:
S=
1000
10
CN
(6.31)
Donde CN es el Nmero de Curva.

Duracin del Hidrograma Unitario, T
El tiempo al pico, TP est dado como:
Tp =
T
+ TR
2
(6.32)
CONVENIO 366 / 2003
6-25

En donde T es la duracin de la lluvia efectiva en horas y TR es el tiempo de rezago en

horas.
La relacin promedio entre el rezago y el tiempo de concentracin, puede escribirse de la
siguiente manera:
TR = 0.6 TC
(6.33)
Los valores del tiempo de rezago segn esta expresin son los mismos obtenidos para el
Hidrograma Unitario de Snyder.
Se puede demostrar que la relacin entre la duracin T de la precipitacin efectiva a la
que se le va a construir el hidrograma unitario y el tiempo de concentracin est dada
como:
T = 0.133 TC (6.34)
6.3.2.3 Modelo Hidrgrafa Sinttica Unitaria de Clark
Este mtodo puede ser aplicado a cuencas no instrumentadas con geomorfologa, suelos,
cobertura vegetal o usos del suelo, y caractersticas climticas similares a las de las
cuencas instrumentadas. El proceso de atenuacin y translacin dominan el movimiento
del flujo gradiente abajo a travs de la cuenca en respuesta a la gravedad. La atenuacin
resulta a partir de la fuerzas friccinales y de los efectos de almacenamiento del canal que
resisten el flujo.
Partiendo de la ecuacin de continuidad
ds
= I t Ot (6.35)
dt
Donde,
ds/dt
It
Ot
:Tasa de cambio del almacenamiento en un tiempo t

:Flujo de entrada promedio en el tiempo t
: Flujo de salida promedio en el tiempo t
La atenuacin del flujo puede ser representado con tanque lineal para el cual el
almacenamiento esta relacionado con el flujo de salida de la siguiente manera:
St = ROt
(6.36)
Donde,
St
R
:Almacenamiento de la cuenca en el tiempo t

: Coeficiente de almacenamiento de la cuenca.
CONVENIO 366 / 2003
6-26

Ot
:Flujo de salida de la cuenca en el tiempo t.
Combinando y solucionando las ecuaciones 6.16 y 6.17 y empleando un modelo en

diferencias finitas aproximado se obtiene la siguiente expresin.
Ot =C A It +CBOt 1
(6.37)
Donde CA y CB son coeficientes que pueden ser obtenidos mediante la siguiente

expresin:
CA =
t
R + 0.5t
(6.38)
CB = 1 C A
El flujo promedio de salida en el instante t se calcula como:
Ot =
Ot 1 Ot
2
(6.39)
Clark (1945) observ que la translacin del flujo a travs de la cuenca puede ser descrito
mediante la curva de tiempo-rea, la cual expresa la curva de la fraccin del rea de
cuenca que contribuye a la escorrenta de la cuenca como una funcin del tiempo desde
que comienza la precipitacin efectiva. La precipitacin efectiva es aquella precipitacin
que ni es retenida en la superficie ni es infiltrada dentro del suelo. La curva tiempo rea es
delimitada en el tiempo por el tiempo de concentracin (Tc) de la cuenca.
El HEC-HMS presenta una curva tiempo - rea tpica obtenida a partir de estudios
realizados en muchas cuencas de Estados Unidos, la cual esta dada por la siguiente
expresin.
1.5
t
t
1.414 para t c
2
At
tc
=
A
t
tc
1 1.4141 t para t 2
c
(6.40)
Donde,
At
A
tc
:Area acumulada contribuyente en el tiempo t

:Area total de la cuenca
:Tiempo de concentracin de la cuenca
La multiplicacin del rea contribuyente en el tiempo t multiplicada por la intensidad de la

lluvia en el instante t produce el caudal de salida Qt.
CONVENIO 366 / 2003
6-27

Numerosos investigadores han encontrado que no es necesario determinar la curva

tiempo-rea para la cuenca. Por ejemplo, Turner y Burdoin (1941) encontraron una
hidrgrafa con resultados relativamente aceptables obtenidas cuando la curva de tiempo rea fue sustituida por la formas geomtricas simples.
En Ilinois, el programa de modelacin hidrolgica, HEC-1 (este software es la versin
anterior del HEC-HMS), es utilizado tpicamente para calcular la hidrgrafa unitaria
sinttica de Clark. Tc y R son los parmetros requeridos para la modelacin en el HEC-1.
El Tc de Clark es ligeramente diferente a la definicin tpica del tiempo de concentracin
aplicado en otras metodologas, tal como el mtodo racional. La definicin tpica del
tiempo de concentracin es: tiempo de viaje para la primera gota de la precipitacin
efectiva en el punto ms distante de la cuenca para alcanzar la salida de la cuenca.
En la hidrografa de Clark, Tc es definido como el tiempo desde el fin de la precipitacin
efectiva hasta el punto de inflexin de la hidrgrafa de escorrenta. El punto de inflexin
de la hidrografa de escorrenta corresponde al tiempo cuando el flujo sobre el terreno
hacia los canales cesa desde el drenaje de los canales de almacenamiento. Por ende, el
Tc de Clark es el tiempo de viaje requerido para que la ltima gota agua de la
precipitacin efectiva en el punto ms distante hidrulicamente en la cuenca alcance la
red de canales. En la mayora de las aplicaciones del HEC-1, Tc es determinado a partir
de valores calibrados con precipitaciones medidas y datos de escorrenta superficial
directa.
Debido a la no presencia de cuencas instrumentadas dentro de zona de estudio, se hace
necesario recurrir a la aplicacin de metodologas regionalizadas en otros pases,
especficamente Estados Unidos, con el propsito de poder obtener los parmetros
requeridos para la modelacin hidrolgica con el hidrograma de Clark.
El coeficiente de atenuacin, el cual representa el efecto de almacenamiento de un canal,
es calculado a partir de un hidrograma observado en la cuenca. Como se muestra en la
Figura 6.8.
Figura 6.8
Calculo del coeficiente de atenuacin de almacenamiento (R)
CONVENIO 366 / 2003
6-28

El valor del coeficiente de alamacenamiento (R) se calcula como:
R=
Q(t )dt
poi
Q poi
(6.41)
La metodologa regionalizada con la cual se estimaron los parmetros R y Tc, fue

implementada en el estado de Illinois, Estados Unidos. Esta metodologa utiliz el anlisis
de regresin mltiple lineal para correlacionar el logaritmo de la media de Tc y la media de
R para cada cuenca con los logaritmos del rea de la cuenca, pendiente y longitud del
canal. Las ecuaciones regionalizadas para Tc y R se expresan a continuacin:
TC = 1.54 L0.875 S 0.181 (6.42)

R = 16.4 L0.342 S 0.790
(6.43)
Donde:
L es la longitud del cauce expresado en millas.
S es la pendiente del canal expresado en ft / milla.
Tanto R como Tc estn expresados en horas.
Empleando esta metodologa se obtuvieron valores demasiado altos para el coeficiente de
almacenamiento, por lo cual se estimo un valor promedio para toda la cuenca trabajando
con una longitud promedio de ladera de 100 m y una velocidad promedio en la misma de
de 14 cm/s, de donde se obtiene un coeficiente de almacenamiento constante de R =0.2
h.
En el Anexo E se presenta un resumen de otros parmetros obtenidos para los diferentes
mtodos de transformacin empleados y para cada cuenca, y en el Anexo F se presentan
los datos empleados como datos de entrada al modelo de HEC-HMS para cada una de
las cuencas agregadas y los diferentes mtodos de transformacin.
6.3.3. Mtodos de Transito en canales
Segn Chow (1994) El transito de canales es un procediendo para determinar el tiempo y
la magnitud del hidrograma de caudal en un punto de un curso de agua utilizando
hidrogramas conocidos o supuestos en uno o mas puntos aguas arriba. Si el flujo es una
creciente, el procedimiento se conoce especficamente como transito de crecientes.
Las metodologas disponibles para el transito de crecientes en canales en el HEC-HMS
son:
CONVENIO 366 / 2003
6-29

6.3.3.1 Lag
Este mtodo supone que el hidrograma de salida del canal es exactamente igual al de
entrada del canal pero retardado. El parmetro requerido es simplemente el retardo lag
en minutos , esto implica tener un conocimiento de las caractersticas de atenuacin del
cauce en cuanto a efectos de traslacin de la onda de creciente. Este mtodo es
principalmente empleado en canales de drenaje urbano (Pilgrim y Cordery, 1993).
Matemticamente las ordenadas del hidrograma aguas abajo en el canal son calculadas
como:
I t
Ot =
I t lag
t < lag
t lag
(6.44)
Donde,
Ot
It
Lag
:Ordenadas del hidrograma de entrada en el tiempo t

:Ordenadas del hidrograma de salida en el tiempo t
:tiempo en el cual las ordenadas del hidrograma de entrada son retardadas
6.3.3.2 Muskingum:
El mtodo de Muskingum es un mtodo de transito hidrolgico que se usa para manejar
relaciones caudal - almacenamiento variables. Inicialmente el esquema se desarrollo
como un sistema de flujo uniforme que esquematiza el almacenamiento volumtrico de
creciente en un tramo de canal o de un ro como la combinacin del almacenamiento por
cua y prisma, como se observa en la Figura 6.9.
qwdd
Figura 6.9
Muskingum. Transito hidrolgico de crecientes en canales.
El esquema de la figura 1.6 se representa matemticamente mediante la siguiente

expresin:
S = KQ + KX ( I Q) (6.45)
CONVENIO 366 / 2003
6-30

Donde,
KQ
KX(I-Q)
: Representa el almacenamiento por cua

: Representa el almacenamiento por prisma
Reordenando la anterior expresin se obtiene la funcin de almacenamiento para el

mtodo de Muskingum:
S = K [XI + (1 X )Q ]
(6.46)
Donde,
K
X
: Coeficiente de proporcionalidad
: Es el factor de ponderacin
El valor de X varia entre 0 y 5 pero para corrientes naturales se encuentra entre 0.y 0.3.
Para la solucin de la funcin de almacenamiento de Muskingum con cambio de
almacenamiento en el tiempo, pueden emplearse tcnicas de diferencias finitas, por
medio de las cuales se llega al siguiente sistema simplificado.
Q j+1 = C1 I j+1 + C 2 I j + C 3 Q j
(6.47)
Donde,
C1 =
t 2KX
2K (1 X ) + t
C2 =
t + 2KX
2K (1 X ) + t
C3 =
2K (1 X ) t
2K (1 X ) + t
(6.48)
Los parmetros K y X pueden ser calibrados cuando se tienen datos del hidrograma de
entrada y de salida observados en el tramo de canal. Esto obviamente hace compleja su
utilizacin para efectos del modelo implementado.
El mtodo de Muskingum originalmente aplicado es criticado por que incluye parmetros
que no tienen fundamentos fsicos y por consiguiente sus valores son difciles de obtener.
Ademas , el modelo se basa en supuestos que no se cumplen en canales naturales
(manual HEC-HMS, 1998).
CONVENIO 366 / 2003
6-31

6.3.3.3 Muskingum Cunge y Muskingum Cunge 8 puntos

El mtodo de Muskingum Cunge o Muskingum Cunge 8 puntos propuesto por Cunge
(1969), combina el mtodo de onda cinemtica y el mtodo de Muskingum tradicional
incluyendo flujos no uniformes y no permanentes en la solucin de la ecuacin 31, de ah
que el caudal en x = (i+1)x y t=(j+1)t (Discretizacin en malla de elementos finitos)
Qi j++11 = C1Qi j +1 + C2Qi j + C3Qi j+1
(6.49)
En esta ultima propuesta a la solucin del transito de crecientes, Cunge demuestra,

despus de aplicar ciertas restricciones (K y t constantes) que la ecuacin 6.46 es una
solucin a la ecuacin de onda cinemtica y que adems esta es un caso especial de la
ecuacin de difusin siempre y cuando:
K=
X=
x
ck
1
Q
1
2 Bc k S 0 x
(6.50)
(6.51)
Donde,
Ck
B
So
: Celeridad de la onda correspondiente a Q y B

: Ancho de la superficie del agua
: Pendiente de la lnea de Energa
El mtodo de Muskingum Cunge se lleva a cabo resolviendo la ecuacin 6.49. Los

coeficientes se calculan utilizando las ecuaciones 6.50 y 6.51 en conjunto con las
ecuaciones 6.48; para una discretizacin dada del espacio y el tiempo.
Una de las principales ventajas de los mtodos de Muskingun- Cunge y Muskingun Cunge 8 puntos es que no requieren la calibracin de los parmetros K y X, pues stos
son calculados durante la simulacin.
Otros mtodos son el de la Onda Cinemtica y el Mtodo del Pulso Modificado.
Se propone emplear el mtodo de Muskingum Cunge puntos para el cual se requiere de
la siguiente informacin.
Longitud de todos los tramos de canal acuerdo a los nodos definidos en el modelo
principal (Estos fueron obtenidos a partir de la red de drenaje generada en SIG ,
empleando el Software Autocad landDesktop).
Pendiente promedio de cada canal: Fueron obtenidas a partir de las curvas de nivel en
forma digital las cuales fueron obtenidas su vez del Modelo de Elevacin Digital
empleando el modulo GIS Analyisis Surface Analysis Contour de IDRISI.
CONVENIO 366 / 2003
6-32

Ancho del fondo del canal: Estos fueron seleccionados de forma aproximada de acuerdo
a la ubicacin del mismo mediante datos de campo ( si se conocen datos de la seccin
transversal puede emplearse el mtodo de Muskingum Cunge 8 puntos).
Coeficiente de rugosidad de Manning: Para la estimacin del nmero de Manning se
tomaron valores presentados en la bibliografa teniendo en cuenta las caractersticas del
material del fondo el cual presenta alto contenido de canto rodado y piedras en las partes
altas y medias de la cuenca (valores entre 0.05 a 0.06) para los canales ubicados en la
parte baja de la cuenca se tomaron valores mas altos (0.06 a 0.08) teniendo en cuenta la
existencia de grandes rocas en algunas zonas, cambios bruscos en el alineamiento, poco
mantenimiento, posos y vegetacin espesa en las mrgenes. Para algunos de los tramos
en la parte canalizada entre Ancon Sur y Acevedo se validaron los valores obtenidos
mediante el software HEC-RAS (ver Captulo 7) con base a los valores obtenidos de la
campaa de aforos. Dada la importancia y sensibilidad del modelo ante dicho parmetro
se recomienda la medicin mas detallada del coeficiente de rugosidad de Manning
mediante muestreos de campo (El mtodo Muskingum Cunge 8 puntos permite considerar
distintos valores de Manning para el fondo y las bancas).
Pendiente del talud: Para esto se asumieron pendientes constante de los taludes de 45
(1H:1V)
En el Anexo G se presentan un resumen de los datos de entrada al modelo para el
mtodo de Musking Cunge para cada uno de los canales definidos en el modelo.
6.3.4. Modelo Meteorolgicos (Definicin de Escenarios de lluvia)
Se definen unos eventos de lluvias mximas posibles que puedan ocurrir sobre la cuenca
y que se traduzcan en la produccin de la escorrenta. Debido a la falta de informacin de
distribucin de lluvias mximas en escala horaria sobre el rea de estudio, se utiliza para
la generacin de hietogramas de precipitacin mediante el mtodo del hietograma de Huff
(Chow V.T., 1994) y el mtodo del Hietograma de Chicago. Este ultimo mtodo utiliza las
curvas Intensidad-Frecuencia-Duracin (IDF) para determinar la secuencia temporal de
precipitacin sobre un rea dada. La ventaja de utilizar este mtodo es que proporciona
una distribucin temporal de precipitacin extractando de las curvas IDF las intensidades
para cualquier duracin dada y para un perodo de retorno dado. A continuacin se
presenta la informacin y los hietogramas de precipitacin tpica para las tormentas
calculadas por los metodos anteriores.
La tormenta de diseo se aplica con la finalidad de poder predecir un evento de
precipitacin para diferentes periodos de retorno, de modo que se pueda tener una
aproximacin de la lluvia que caera en determinado lugar, de las diferentes cuencas que
componen toda la red del Ro Medelln.
La estimacin de las tormentas de diseo parte inicialmente de una zonificacin que se
lleva a cabo en todas las subcuencas que son tributarias de la red de drenaje del Ro
CONVENIO 366 / 2003
6-33

Medelln. Esta zonificacin se realiza con base en los mapas de precipitacin total
multianual de toda el rea de estudio, y asimismo de las caractersticas topogrficas de la
zona, es decir, empleando la cartografa y el Modelo de Elevacion Digital (DEM)
disponible se identifican las principales barreras topogrficas que puedan incidir sobre la
distribucin de la lluvia, como se muestra en la Figura 6.10(a)
La zonificacin se realiz en un sistema de informacin geogrfica, SIG, en el cual se
puede observar toda el rea de drenaje del Ro Medelln, dividido por subcuencas y cada
una de stas con sus atributos especficos, como pendiente, usos del suelo, entre otras.
El SIG permite observar los diferentes rasgos del relieve para un mejor entendimiento del
comportamiento de las lluvias en la zona. En la Figura 6.7(b) se presenta la zonificacin
de tormentas para la cuenca del Ro Medelln.
Figura 6.10
(a) caracterizacin topogrfica (b) Zonificacin de las tormentas de diseo.
Una vez definida la zonificacin de toda el rea tributante, se procedi a estimar el tiempo
de concentracin para cada cuenca (ver Anexo E) y luego con base a stos, se procedi a
determinar el tiempo de concentracin media para cada de las zonas que forman el rea
tributante, el resultado de esta metodologa es una lluvia de corta duracin de alta
intensidad la cual se aplica como lluvia de diseo en el HEC-HMS para cada zona.
Existen diversas metodologas para determinar el tiempo de concentracin de una cuenca
a partir de los parmetros morfomtricos, bsicamente de la pendiente y longitud del
canal y de la cuenca, todas ellas determinadas a partir de ajustes empricos de registros
hidrolgicos. La duracin de la lluvia se hace igual al tiempo de concentracin de la
cuenca, puesto que es, para esta duracin, cuando la totalidad de la cuenca est
aportando al proceso de escorrenta, por esto, es de esperarse que se presenten los
caudales mximos.
En la literatura existen mltiples expresiones para el clculo del tiempo de concentracin:
propuestas por Temez, William, Kirpich, Johnstone y Cross, California Coulverts Practice,
CONVENIO 366 / 2003
6-34

Giandotti, Ventura -Heron, Brausby-William. Debido a las diferentes formas en las fueron
concebidas estas expresiones, la variabilidad en los resultados de una a otra es bastante
alta por lo cual se hace necesario escoger el tiempo de concentracin de forma apropiada
descartando aquellos mtodos que no se ajusten a las condiciones morfomtricas locales
de la cuenca.
Las expresiones usadas para el clculo de los tiempos de concentracin (Tc) sugeridas
por Smith y Vlez (1997) son:
Tmez(1978)
L
Tc = 0.3 0.25
So
0.75
(6.52)
Tc: Tiempo de concentracin en horas, L: Longitud del cauce principal en kilmetros, So:
Diferencia de cotas sobre L en porcentaje.
Williams (1922)
L A 0 .4
Tc =
DS 0o.2
(6.53)
A: rea de la cuenca en millas cuadradas, L: distancia en lnea recta desde el sitio de

inters al punto mas alto en millas, So: diferencia de cotas entre los puntos ms extremos
divida por L en porcentaje, d: dimetro de una cuenca circular con rea A en millas.
Kirpich (1990)
Tc = 0.066
S
0.77
(6.54)
L: longitud desde la estacin de aforo hasta la divisoria siguiendo en cauce principal en

kilmetros, So: diferencia de cotas entre los puntos extremos de la corriente en m/m.
Johnstone y Cross (1949)
Tc = 5
S
0 .5
(6.55)
L: longitud del cauce principal en millas, So: pendiente del canal en pies/milla.
California Culverts Practice (1942)
L
Tc = 0.01
H
0.385
(6.56)
L: longitud del canal principal en kilmetros, H: diferencia de cotas entre el punto de

inters y la divisoria en metros.
CONVENIO 366 / 2003
6-35

Giandiotti (1990)
Tc =
4 A + 1 .5 L
25.3 LS o
(6.57)
A: rea de la cuenca en kilmetros cuadrados, L: Longitud del cauce principal en

kilmetros, So: diferencia de cotas entre puntos extremos de la corriente en m/m.
Snyder
( )
Tc = C t L L
0.3
(6.58)
L: longitud de la cuenca en millas, L : distancia desde el punto de inters al centro de

gravedad de la cuenca en millas, C: es una constante (1.2 en reas montaosas, 0.72 en
zonas de pie de ladera y 0.35 en valles).
Luego de aplicar estas expresiones para la estimacin de Tc, se adopta tomar la media
del tiempo de concentracin para cada una de las zonas. Esto se hace con el propsito de
no generar tendencias muy marcadas en el valor de Tc, lo cual conllevara a simulaciones
hidrolgicas diferentes a las reales.
En la Tabla 6.13 se presenta las estaciones seleccionadas para la determinacin de la
elaboracin de las IDF.
Tabla 6.13 Estaciones con IDF utilizadas para el diseo de las tormentas
CDIGO
NOMBRE
2308024 Vasconia
2701034 Mazo
2701035 Chorrillos
2701036 Caldas
2701037 Fabricato
2701038 San Antonio de Prado
2701045 Villahermosa Planta
2701046 San Cristbal
2701057 Barbosa
2701066 Gabino
2701093 Ayur
2701507 Apto Olaya Herrera
2701053 Alto San Andrs
MUNICIPIO
Rionegro
Medelln
Medelln
Caldas
Bello
Medelln
Medelln
Medelln
Barbosa
Santo Domingo
Envigado
Medelln
Girardota
COTA
844895
841940
842280
828568
831500
824640
837340
827560
861540
875150
835380
831858
849230
1178000
1184200
1188220
1160510
1195470
1175253
1183990
1186530
1204180
1216360
1173830
1179127
1203540
2510
2455
2353
1875
2422
2000
1690
1890
1290
1080
1770
1490
2240
En las zonas donde se presenta ms de una estacin pluviomtrica, se le asigna a cada

estacin un determinado peso utilizando la tcnica de polgonos de Thiessen. Las
estaciones y los valores de las IDF se presentan en la Tabla 6.14.
Para determinar la intensidad de lluvia asociada a diferentes perodos de retorno se
utilizan las curvas IDF (Intensidad-Frecuencia-Duracin), que relacionan la duracin de la
lluvia con la intensidad de sta para diferentes perodos de retorno. Las curvas IDF estn
representadas por una ecuacin de la siguiente forma para las zonas 1, 2, 3, 4 y 5:
CONVENIO 366 / 2003
6-36

I (mm / h) = C (h + D) m
Donde: c, h y m son parmetros constantes de la estacin; I es la intensidad (mm/h); TR
es el perodo de retorno de inters (aos); D es la duracin de la lluvia (minutos).
En las estaciones donde no se cuenta con un registro necesario para determinar la curva
IDF, se implementa la metodologa propuesta por Wilches (1999), la cual permite
determinar curvas IDF a partir de un periodo corto de registros. Para la zona 6 se utiliz
la siguiente expresin para calcular la intensidad de la Lluvia:
46.2 43.05
I(mm/h) = K 0.75
(6.59)
d
d
Donde: k es un parmetro calculado por la metodologa de Wilches para cada perodo de
Retorno y d es la duracin de la lluvia en minutos.
Tabla 6.14 Parmetros obtenidos de las curvas IDF para cada una de las zonas
Zona 1: Sector Caldas
Tr
(Aos)
2.33
5
10
25
50
100
c
162369
30819.79
13123.79
9805.85
9515.12
8923.1
Estacin
Caldas
h
58
34
24
18
16
14
m
-1.72029
-1.41291
-1.22745
-1.15086
-1.1295
-1.10273
Zona 2: Sector Ancn sur

Tr
(Aos)
2.33
5
10
25
50
100
Estacin
San Antonio de Prado
c
h
m
c
71663.74
50
-1.5677
5612.910
21973.75
32
-1.34471
7081.760
10254.29
22
-1.18512 10204.550
11629.01
20
-1.19038 14533.170
11495.52
18
-1.17556 16923.510
10953.22
16
-1.15445 20165.180
Ayur
h
22.000
18.000
18.000
18.000
18.000
18.000
m
-1.127
-1.155
-1.211
-1.262
-1.258
-1.276
Olaya Herrera
h
12.00
10.00
8.00
6.00
6.00
10.00
m
-0.97
-0.89
-0.83
-0.78
-0.76
-0.86
Zona 3: Sector Medelln

Tr
(Aos)
2.33
5
10
25
50
100
Estacin
C
28902.26
6183.28
2527.91
3817.56
2298.3
4531.96
San Cristbal
h
56
28
14
14
8
12
m
-1.376
-1.082
-0.901
-0.964
-0.839
-0.964
c
1678.78
1434.62
1231.60
1057.15
1113.70
1929.04
CONVENIO 366 / 2003
6-37

Tr
Aos
2.33
5
10
25
50
100
c
14613.09
7283.46
3854.3
3722.39
3389.62
3716.04
Villa Hermosa
h
30
22
16
14
12
12
m
-1.32762
-1.16758
-1.01342
-0.97874
-0.94236
-0.94303
c
5106.01
3669.03
3613.09
3878.45
2510.39
3246.96
Vasconia
h
16
10
8
6
2
2
m
-1.07642
-0.99514
-0.9707
-0.96827
-0.85871
-0.90548
Zona 4: Sector Bello

Tr
(Aos)
2.33
5
10
25
50
100
Tr
(Aos)
2.33
5
10
25
50
100
Estacin
c
6321.12
3531.49
3526.21
4240.8
4270.88
4974.99
Fabricato
h
38
22
18
18
16
16
m
c
-1.09507 123854.2
-0.97203 11521.81
-0.95604 7171.86
-0.95738 5153.19
-0.93943 5410.29
-0.95062 4444.21
c
125463.1
37374.23
22743.7
10713.66
10908.04
10676.29
Estacin
Mazo
h
66
42
32
22
20
18
m
-1.66488
-1.44639
-1.34665
-1.17491
-1.16256
-1.14436
Chorrillos
h
68
34
26
18
16
12
m
-1.66941
-1.23658
-1.13258
-1.05198
-1.04871
-0.99994
Zona 5: Sector Ancn Norte

Tr
(Aos)
2.33
5
10
25
50
100
c
2603.4
2103.46
2047.97
2106.21
1977.72
2224.58
Estacin
San Andrs
h
20
14
12
10
8
8
m
-0.95913
-0.9066
-0.88165
-0.85802
-0.82659
-0.83197
Zona 6: Sector Norte

Tr
(Aos)
2.33
5
10
25
50
100
Estacin
Barbosa
Gabino
k
k
26.465842
32.010
31.274333
36.947
34.936882
40.635
39.316355
44.974
42.433225
48.021
45.447276
50.937
CONVENIO 366 / 2003
6-38

La magnitud de la precipitacin se determina a partir de la intensidad de lluvia para

diferentes perodos de retorno mediante la siguiente expresin, donde i es la intensidad
en mm/h, d es la duracin de la lluvia en minutos.
P = i*
d
60
(6.60)
Las intensidades y precipitaciones asociadas a los tiempos de concentracin de cada

zona para diferentes perodos de retorno, se aprecian en la Tabla 6.15.
Tabla 6.15 Intensidades y Precipitaciones para diferentes perodos de retorno.
Tr(aos)
Zona
2,33
5
10
25
50
100
Tr(aos)
Zona
2,33
5
10
25
50
100
Tr(aos)
Zona
2,33
5
10
25
50
100
I(mm/h)
59,90
69,11
78,13
89,84
98,89
107,49
P(mm)
40,93
47,22
53,39
61,39
67,58
73,45
Zona
I(mm/h)
47,00
54,66
60,29
68,94
75,40
83,28
P(mm)
35,25
41,00
45,22
51,70
56,55
62,46
Zona
I(mm/h)
62,29
69,51
77,52
90,85
99,98
110,29
P(mm)
30,11
33,60
37,47
43,91
48,32
53,31
Zona
I(mm/h)
60,29
70,44
79,14
91,53
103,76
113,67
P(mm)
38,18
44,61
50,12
57,97
65,72
71,99
I(mm/h)
53,65
62,88
70,85
83,54
93,07
102,32
P(mm)
33,98
39,82
44,87
52,91
58,95
64,80
I(mm/h)
62,21
72,58
80,40
89,68
96,24
102,55
P(mm)
32,14
37,50
41,54
46,33
49,72
52,98
Para la distribucin temporal de la lluvia se utiliz el diagrama elaborado por Huff, 1967,
para lluvia con una probabilidad de excedencia del 50%. As, para incrementos en el
tiempo de 10%, se tiene el comportamiento de las lluvias mostrado en la Tabla 6.16.
Tabla 6.16 Porcentajes de distribucin temporal de la lluvia para las estaciones empleadas
(probabilidad de excedencia de 50%).
%T
%P
0
0
10
18
20
51
30
60
40
79.5
50
85
60
89
70
92
80
95
90
98
100
100
CONVENIO 366 / 2003
6-39

HIETOGRAMA TOTAL ZONA SUR T01
HIETOGRAMA TOTAL ZONA ANCN SUR T02
25
25
Lmina (mm)
Tr = 2.33 aos
15
Tr = 5 aos
10
Tr = 10 aos
Lmina (mm)
20
20
Tr = 5 aos
Tr = 10 aos
10
Tr = 25 aos
5
Tr = 25 aos
Tr = 2.33 aos
15
Tr = 50 aos
Tr = 50 aos
Tr = 100 aos
Tr = 2.33 aos
Tr = 100 aos
Tr = 2.33 aos
12
16
21
25
29
33
37
11
15
19
23
27
30
34
38
41
Tiempo acumulado (min)
HIETOGRAMA TOTAL ZONA MEDELLN T03
HIETOGRAMA TOTAL ZONA BELLO T04

25
25
20
15
Tr = 2.33 aos
Tr = 5 aos
10
Tr = 10 aos
Tr = 25 aos
Lmina (mm)
Lmina (mm)
20
Tr = 2.33 aos
15
Tr = 5 aos
Tr = 10 aos
10
Tr = 25 aos
Tr = 50 aos
Tr = 50 aos
Tr = 100 aos
Tr = 2.33 aos
14
18
23
27
32
36
41
Tr =aos
100
Tr = 2.33
45
11
15
19
23
27
30
34
aos
38
HIETOGRAMA TOTAL ANCON NORTE T05
HIETOGRAMA TOTAL NORTE T06
18
18
16
16
14
12
Tr = 2.33 aos
10
8
Tr = 5 aos
Tr = 10 aos
Tr = 25 aos
Tr = 50 aos
Tr = 100 aos
Tr = 2.33 aos
3
12
15
17
20
23
26
29
Lmina (mm)
Lmina (mm)
14
12
10
Tr = 2.33 aos
Tr = 5 aos
Tr = 10 aos
Tr = 25 aos
Tr = 50 aos
= aos
100
Tr =Tr
2.33
0
3
12
16
19
22
25
28
aos
31
Figura 6.11 Hietogramas de precipitacin por zonas para tormentas con Tr = 2.33, 5,
10, 25, 50 y 100 aos
Para el mtodo del Hietograma de Chicago Se utiliz la informacin de las curva IDF para
las estaciones existentes en cada zona. Se seleccion una duracin de tormenta de 3
horas por ser un valor acorde con las caractersticas de lluvias de gran intensidad en la
regin de Antioquia y a que las valores obtenidos para este tipo de lluvias al principio y al
final de la misma son pequeos respecto a los valores obtenidos en la parte central del
periodo de ocurrencia de la tormenta (Valores cercano a 1.5 h). En la Figura 6.12 se
presentan los resultados obtenidos para las lluvias de Chicago.
Se recomienda la modelacin de la cuenca empleando una condicin aun mas critica para
una lluvia de Chicago con duracin igual al tiempo de concentracin de la cuenca (8 horas
aproximadamente).
CONVENIO 366 / 2003
6-40

Zona 1
Precipitacin (mm)
Precipitacin (mm)
10,00
9,00
8,00
7,00
6,00
5,00
4,00
3,00
2,00
1,00
0,00
30,00
L m in a ( m m )
L m in a ( m m )
25,00
20,00
15,00
10,00
5,00
4
12
20
28
36
44
52
60
68
76
84
92
100
108
116
124
132
140
148
156
164
172
180
52
60
68
76
84
92
10
0
10
8
11
6
12
4
13
2
14
0
14
8
15
6
16
4
17
2
18
0
4
12
20
28
36
44
0,00
Tiempo (min)
Tiempo (min)
Tr 2.33 aos
Tr 100 aos
Zona 2
Precipitacin (mm)
Precipitacin (mm)
30,00
10,00
9,00
8,00
7,00
6,00
5,00
4,00
3,00
2,00
1,00
0,00
L m in a ( m m )
L m in a ( m m )
25,00
20,00
15,00
10,00
5,00
4
12
20
28
36
44
52
60
68
76
84
92
100
108
116
124
132
140
148
156
164
172
180
68
76
84
92
10
0
10
8
11
6
12
4
13
2
14
0
14
8
15
6
16
4
17
2
18
0
4
12
20
28
36
44
52
60
0,00
Tiempo (min)
Tiempo (min)
Tr 2.33 aos
Tr 100 aos
Zona 3
Precipitacin (mm)
Precipitacin (mm)
30,00
10,00
25,00
6,00
4,00
20,00
15,00
10,00
2,00
5,00
0,00
0,00
Tiempo (min)
Tr 2.33 aos
4
12
20
28
36
44
52
60
68
76
84
92
100
108
116
124
132
140
148
156
164
172
180
L m in a (m m )
8,00
4
12
20
28
36
44
52
60
68
76
84
92
100
108
116
124
132
140
148
156
164
172
180
L m in a ( m m )
12,00
Tiempo (min)
Tr 100 aos
CONVENIO 366 / 2003
6-41

Zona 4
Precipitacin (mm)
Precipitacin (mm)
25,00
8,00
20,00
6,00
L m in a ( m m )
L m in a ( m m )
7,00
5,00
15,00
4,00
10,00
3,00
2,00
5,00
1,00
0,00
4
12
20
28
36
44
52
60
68
76
84
92
100
108
116
124
132
140
148
156
164
172
180
4
12
20
28
36
44
52
60
68
76
84
92
100
108
116
124
132
140
148
156
164
172
180
0,00
Tiempo (min)
Tiempo (min)
Tr 2.33 aos
Tr 100 aos
Zona 5
Precipitacin (mm)
Precipitacin (mm)
9,00
8,00
7,00
6,00
5,00
4,00
3,00
2,00
1,00
0,00
4
12
20
28
36
44
52
60
68
76
84
92
100
108
116
124
132
140
148
156
164
172
180
4
12
20
28
36
44
52
60
68
76
84
92
10
0
10
8
11
6
12
4
13
2
14
0
14
8
15
6
16
4
17
2
18
0
L m in a ( m m )
L m in a ( m m )
20,00
18,00
16,00
14,00
12,00
10,00
8,00
6,00
4,00
2,00
0,00
Tiempo (min)
Tiempo (min)
Tr 2.33 aos
Tr 100 aos
Zona 6
Precipitacin (mm)
12,00
10,00
8,00
6,00
4,00
2,00
84
92
10
0
10
8
11
6
12
4
13
2
14
0
14
8
15
6
16
4
17
2
18
0
4
12
20
28
36
44
52
60
68
76
0,00
Tiempo (min)
4
12
20
28
36
44
52
60
68
76
84
92
100
108
116
124
132
140
148
156
164
172
180
L m in a ( m m )
L m in a ( m m )
Precipitacin (mm)
20,00
18,00
16,00
14,00
12,00
10,00
8,00
6,00
4,00
2,00
0,00
Tiempo (min)
Tr 2.33 aos
Tr 100 aos
Figura 6.12 Hietogramas de precipitacin de Chicago por zonas para periodos de
retorno de Tr = 2.33 y 100 aos
CONVENIO 366 / 2003
6-42

6.3.5. Resultados
Para la conformacin general del sistema de drenaje de la cuenca del Ro Medelln se
realiz una agrupacin de subcuencas y canales teniendo en cuenta los criterios
mencionados anteriormente.
El modelamiento se inicia desde la parte alta de la cuenca en proceso secuencial, pues
los hidrogramas de salida de las cuencas son los hidrogramas de entrada para los tramos
de canal y as los hidrogramas de salida de los tramos superiores son los hidrogramas de
entrada de los inmediatamente inferiores. En la Figura 6.13 se presenta la
esquematizacin en HEC-HMS de una zona de las zonas que drenan al ro Medelln, se
puede ver que la discretizacin comienza desde las zonas ms altas, de modo que la
direccin del flujo quede definida desde las partes altas hacia las ms bajas.
Q.La Iguana
Q. Santa Elena
Q.El Hueso
Q. El Indio
Q.La Picacha
Q.La Poblada
Q.Altavista
Q.Doa Maria
Q.LaAguacatala
Q.La Ayura
Q.La Grande
Q.La Doctora
Figura 6.13 Representacin de la parte alta y baja de la cuenca del Ro Medelln para
realizar la modelacin usando HEC-HMS.
Definicin de Escenarios de Modelacin
Basados en el estudio de la variabiliad espacial de las tormentas en la cuenca del ro
Medelln y con las lluvias de diseo obtenidas para las diferentes zonas, los diferentes
periodos de retorno y los dos mtodos empleados se definieron los siguientes escenarios.
Escenario 1: Este escenario presenta un aguacero con un periodo de retorno de 2.33
aos en la parte alta de la cuenca en la zona 1 (Sur Caldas) y ausencia de lluvia en el
resto de la cuenca.
Escenario 2: Este escenario presenta el aguacero con un periodo de retorno de 2.33
aos en la media alta de la cuenca en la zona 2 (Ancn Sur) y ausencia de lluvia en el
resto de la cuenca.
CONVENIO 366 / 2003
6-43


aos en la parte media de la cuenca en la zona 3 (Medelln) y ausencia de lluvia en el
resto de la cuenca.
aos en la parte media de la cuenca en la zona 4 (Bello) y ausencia de lluvia en el resto
de la cuenca.
aos en la parte media baja de la cuenca en la zona 5 (Ancon Norte) y ausencia de lluvia
en el resto de la cuenca.
aos en la parte baja de la cuenca en la zona 6 (Norte) y ausencia de lluvia en el resto de
la cuenca.
Escenario 7: Este escenario presenta un aguacero con un periodo de retorno de 100
aos en la parte alta de la cuenca en la zona 1 (Sur Caldas) y un aguacero con un periodo
de retorno de 2.33 en el resto de la cuenca.
Escenario 8: Este escenario presenta el aguacero con un periodo de retorno de 100 aos
en la parte media alta de la cuenca en la zona 2 (Ancn Sur) y un aguacero con un
periodo de retorno de 2.33 en el resto de la cuenca.
Escenario 9: Este escenario presenta el aguacero con un periodo de retorno de 100 aos
en la parte media de la cuenca en la zona 3 (Medelln) y un aguacero con un periodo de
retorno de 2.33 en el resto de la cuenca.
Escenario 10: Este escenario presenta el aguacero con un periodo de retorno de 100
aos en la parte media de la cuenca en la zona 4 (Bello) y un aguacero con un periodo de
aos en la parte media baja de la cuenca en la zona 5 (Ancon Norte) y un aguacero con
un periodo de retorno de 2.33 en el resto de la cuenca.
aos en la parte baja de la cuenca en la zona 6 (Norte) y un aguacero con un periodo de
CONVENIO 366 / 2003
6-44

En la Figura 6.14 se presenta un ejemplo de los resultados grficos presentados por el

software HEC-HMS para el punto de control correspondiente a Puente Gabino y el
escenario 9.
Figura 6.14 Presentacin grfica de resultados Estacin de control Puente Gabino

(Software HEC-HMS).
En la Tabla 6.17 y Tabla 6.18 se presenta un resumen de los resultados empleando los
mtodos de transformacin de lluvia de Snyder, Clark y SCS en los puntos de control
correspondientes a Doa Mara, Metromezclas y Puente Gabino y para las lluvias de
diseo de corta duracin obtenidas a partir del tiempo de contraccin ms alto para cada
zona y mediante el mtodo de Chicago respectivamente.
CONVENIO 366 / 2003
6-45

Tabla 6.17 Valores de caudales Mximos obtenidos para diferentes puntos de control con
mtodo de lluvia de corta duracin (m3/s).
CONVENIO 366 / 2003
6-46

Tabla 6.18 Valores de caudales Mximos obtenidos para diferentes puntos de control en
mtodo de lluvia de Chicago (m3/s).
Puede observarse como para los escenarios del 1 al 6 en los cuales se considera un
aguacero de 2.33 aos de periodo de retorno correspondiente a una lluvia convectiva de
corta duracin en cada una de las zonas, los caudales mximos obtenidos a lo largo de la
cuenca representan adecuadamente la concentracin de la escorrenta en la zona de
influencia del aguacero; atenuando la onda de creciente a lo largo de los canales aguas
abajo. Por ejemplo en el escenario 1 para el cual se tiene el aguacero de 2.33 aos de
periodo retorno en la parte alta de la cuenca en la zona da Caldas (Zona 1) el caudal
obtenido en el punto de control Doa Maria (Zona 2) fue de 492.97 m3/s y aguas abajo en
el punto de control correspondiente a Metromezclas (zona 4 -Sin Lluvia) se presenta un
caudal mximo de 132.23 m3/s atenuando la creciente a lo largo del canal hasta dicho
punto. En puente Gabino el caudal obtenido corresponde al caudal base seleccionado
para el modelo en dicho punto de 87.70 m3/s.
Si se observan los resultados obtenidos para los escenarios 1 y 2 en el mismo punto de
control se observa como para los mtodos de Snyder y S.C.S se presenta un valor un
poco menor para el escenario 2 en el cual se presenta una tormenta de 2.33 aos de
periodo de retorno correspondiente a la zona de ubicacin de dicho punto de control, esto
puede deberse a que el punto de control Metromezclas se encuentra en la parte media de
CONVENIO 366 / 2003
6-47

la zona 2 presentando por lo tanto un rea contribuyente al caudal mucho menor que para
el escenario 1 en el cual se tiene toda el rea de la zona 1aguas arriba contribuyendo a
dicho caudal .
Para los escenarios correspondientes a tormentas con periodos de retorno de 100 aos
en cada zona y de periodo de retorno de 2.33 en el resto (Escenario 7 a 12) se observa
tambin la influencia de la cercana de las tormentas ms intensas en los caudales
mximos obtenidos y de las zonas urbanas en el aumento de la escorrenta superficial
directa.
6.4 ANLISIS COMPARATIVO
El mtodo GRADEX, aplicado en la estimacin de caudales mximos, posee en su
esencia hiptesis que an son cuestionables por investigadores del rea y que hacen del
mtodo slo una alternativa que permite validar o calibrar otras metodologas aplicadas
como regionalizacin de caractersticas medias o trnsito Hidrulico (HEC HMS). Sin
embargo, los caudales mximos obtenidos por el mtodo de regionalizacin de
caractersticas medias y por el mtodo Gradex son similares para periodos de retorno
bajos y para cuencas pequeas, las diferencias entre ambas metodologas son ms
notorias cuando se estiman caudales mximos para periodos de retorno altos y para
cuencas grandes. La Figura 6.15 a Figura 6.20 muestran los caudales estimados por
ambas metodologas.
ANLISIS COMPARATIVO DE CAUDALES MXIMOS (Tr =2.33)
900
800
700
Q (m3/s)
600
500
Regionalizacin
400
Gradex
300
200
100
0
E1 E2 E3 E4 E5 E6 E7 E8 E9 E10 E11 E12 E13 E14 E15 E16 E17 E18 E19 E20 E21
Estaciones
Figura 6.15
Caudales mximos, Tr = 2.33 aos.
CONVENIO 366 / 2003
6-48

ANLISIS COMPARATIVO DE CAUDALES MXIMOS (Tr =5)

1000
900
800
Q (m3/s)
700
600
Regionalizacin
500
Gradex
400
300
200
100
0
Estaciones
Figura 6.16
Caudales mximos, Tr = 5 aos.

1200
1000
Q (m3/s)
800
Regionalizacin
600
Gradex
400
200
0
Estaciones
Figura 6.17
CONVENIO 366 / 2003
6-49


1400
1200
Q (m3/s)
1000
800
Regionalizacin
600
Gradex
400
200
0
Estaciones
Figura 6.18

1600
1400
Q (m3/s)
1200
1000
Regionalizacin
800
Gradex
600
400
200
0
Estaciones
Figura 6.19
CONVENIO 366 / 2003
6-50


1800
1600
1400
Q (m3/s)
1200
1000
Regionalizacin
Gradex
800
600
400
200
0
Estaciones
Figura 6.20
Con respecto al modelamiento distribuido de caudales se puede decir lo siguiente:

El modelo implementado permite observar la influencia de la cercana de las tormentas
ms intensas y de las zonas con alta componente urbana en los caudales mximos
obtenidos. Permite el planteamiento de escenarios de lluvia mas ajustados a las
condiciones reales para una cuenca de tal magnitud y caractersticas climticas como la
del rio Medelln (Lluvias convectivas concentradas en determinada zona de gran
intensidad y corta duracin), as como la inclusin de escenarios que combinen lluvias de
diferentes periodos de retorno y diferentes duraciones totales de la tormenta.
Mediante el modelo es posible considerar de forma un poco mas detallada la influencia en
la variabilidad espacial de variables geomorfolgicas como pendiente de la cuenca,
longitud del cauce principal, usos y tipos del suelo en la concentracin de la escorrentia y
la forma de la respuesta hidrolgica de la cuenca ; as como observar de forma
aproximada la influencia de los cambios en las caractersticas del cauce principal como
pendiente, rugosidad, ancho y forma de la seccin en el transito de la onda de creciente.
El modelo permite realizar de forma aproximada un anlisis de sensibilidad de los
parmetros empleados para la simulacin as como la posible calibracin de algunos de
ellos, teniendo en cuenta obviamente las limitaciones propias de los mtodos de
transformacin, perdidas y transito de crecientes empleadas. Desde el punto de vista
prctico presenta la ventaja de permitir la estimacin de caudales mximos aproximados
en diferentes puntos de la red principal a lo largo del ro Medelln.
Otra de las ventajas del modelo es que permite de forma fcil la Implementacion de
descargas puntuales a lo largo del cauce o de hidrogramas que representen los aportes
de alcantarillados de aguas lluvias durante una tormenta determinada mejorando as la
modelacin de crecientes en cuencas con alta componente urbana.
CONVENIO 366 / 2003
6-51

Como recomendaciones a futuro se pretende en primera instacia, la implementacin de

estudios para la estimacin valores ms precisos de la rugosidad a lo largo de canal de
manera que se consideren aspectos como presencia de vegetacin, variaciones en la
seccin transversal a lo largo del canal, alineamiento, profundidad de flujo, etc. de manera
que permitan mejorar los valores empleados durante la simulacion. Adicionalmente
estudios que permitan determinar parmetros asociados con las perdidas por infiltracin
de forma mas precisa considerando la influencia de variables como pendiente,usos del
suelo, intensidad de la lluvia.
Se recomienda emplear el HEC-HMS como herramienta para el apoyo de toma de
decisiones sin olvidar los posibles errores asociados a los diferentes mtodos empleados
en la modelacin de la cuenca del rio Medelln (mtodos de transformacin, perdidas y
transito de crecientes). Es de vital importancia avanzar hacia la Implementacin de
modelos espacialmente distribuidos que consideren de forma mas detallada la fsica del
fenmeno Lluvia Escorrentia y el transito de crecientes en canales.
CONVENIO 366 / 2003
6-52

Hidraulica e Hidrologia - Capitulo 6 - Q Max

Cargado por

Copyright:

Formatos disponibles

Hidraulica e Hidrologia - Capitulo 6 - Q Max

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Hidraulica e Hidrologia - Capitulo 6 - Q Max

Cargado por

Copyright:

Formatos disponibles

DISEO Y PUESTA EN MARCHA DE LA RED DE MONITOREO AMBIENTAL EN LA CUENCA HIDROGRFICA DEL

RO MEDELLN EN JURISDICCION DEL AREA METROPOLITANA

6.1 MTODO GRADEX

P es la precipitacin mxima sobre la cuenca, y corresponden a la media y desviacin

REA METROPOLITANA DEL VALLE DE ABURRA

CONVENIO 366 / 2003

DISEO Y PUESTA EN MARCHA DE LA RED DE MONITOREO AMBIENTAL EN LA CUENCA HIDROGRFICA DEL

donde V (m) es el volumen de agua que ingresa a la cuenca en forma de precipitacin,

donde T0 corresponde al perodo de retorno de saturacin.

El Grdex de caudal, q, expresado en m/s se obtiene a partir del Grdex de

CONVENIO 366 / 2003

DISEO Y PUESTA EN MARCHA DE LA RED DE MONITOREO AMBIENTAL EN LA CUENCA HIDROGRFICA DEL

6.1.3. Obtencin del caudal pivote

REA METROPOLITANA DEL VALLE DE ABURRA

CONVENIO 366 / 2003

DISEO Y PUESTA EN MARCHA DE LA RED DE MONITOREO AMBIENTAL EN LA CUENCA HIDROGRFICA DEL

Estaciones de precipitacin utilizadas en la metodologa gradex

De las 39 estaciones de precipitacin seleccionadas para aplicar la metodologa Grdex,

LOCALIZACION ESCALA LOCALIZACION ESCALA

REA METROPOLITANA DEL VALLE DE ABURRA

CONVENIO 366 / 2003

DISEO Y PUESTA EN MARCHA DE LA RED DE MONITOREO AMBIENTAL EN LA CUENCA HIDROGRFICA DEL

REA METROPOLITANA DEL VALLE DE ABURRA

CONVENIO 366 / 2003

DISEO Y PUESTA EN MARCHA DE LA RED DE MONITOREO AMBIENTAL EN LA CUENCA HIDROGRFICA DEL

Estaciones auxiliares de precipitacin utilizadas en la metodologa grdex

Mapa de Grdex de precipitacin

REA METROPOLITANA DEL VALLE DE ABURRA

CONVENIO 366 / 2003

DISEO Y PUESTA EN MARCHA DE LA RED DE MONITOREO AMBIENTAL EN LA CUENCA HIDROGRFICA DEL

Alto San Miguel

6.1.5.2 Gradex de caudal

CONVENIO 366 / 2003

DISEO Y PUESTA EN MARCHA DE LA RED DE MONITOREO AMBIENTAL EN LA CUENCA HIDROGRFICA DEL

(Km2) corresponde al rea de la cuenca, p (mm) al valor medio del gradex de

En Tabla 6.3 se presenta el valor medio de ste parmetro en la regin de estudio.

Alto San Miguel

6.1.5.3 Caudal pivote

CONVENIO 366 / 2003

DISEO Y PUESTA EN MARCHA DE LA RED DE MONITOREO AMBIENTAL EN LA CUENCA HIDROGRFICA DEL

Q2.33 = 0.1308 A1.005

6.1.5.4 Avenidas mximas

Caudal mximo (m /s)

Perodo de retorno (aos)

Alto San Miguel

Antes San Fernando

Despus San Fernando

Parque de Las Aguas

Quebrada Doa Mara

REA METROPOLITANA DEL VALLE DE ABURRA

CONVENIO 366 / 2003

DISEO Y PUESTA EN MARCHA DE LA RED DE MONITOREO AMBIENTAL EN LA CUENCA HIDROGRFICA DEL

Las metodologas empleadas en la estimacin de caudales mximos no incluyen los

QTr = Q max + K Q max

Utilizando la metodologa propuesta por UNALMED-UPME (2000), la media y la

K y K son el factor de frecuencia para la media y la desviacin estndar,

respectivamente; P, la precipitacin media sobre la cuenca; E, la evapotranspiracin; A, el

CONVENIO 366 / 2003