Geomorfología Quebrada La Iguaná
Geomorfología Quebrada La Iguaná
Geomorfología Quebrada La Iguaná
Trabajo de Grado
Presentado como requisito para optar el Título
de Maestría en Ciencias de la Tierra
Por:
Director:
Michel Hermelin
Asesores:
Geovany Bedoya S.
Rafael Madrigal
Universidad EAFIT
Maestría en Ciencias de la Tierra
15 de Agosto de 2001
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Tabla de contenido
Página
1 Introducción 1
1.1 Objetivos generales 3
1.2 Objetivos específicos 4
1.3 Descripción de la zona de estudio 4
1.3.1 Zona suroccidental 6
1.3.2 Zona centro - occidental, cuenca de la quebrada La Iguaná 9
1.3.3 Zona centro – oriental 10
1.3.4 Zona noroccidental 11
1.3.5 Cuenca de la quebrada Doña María, corregimiento de San Antonio de 12
Prado
2 Metodología 14
2.1 Revisión de información secundaria 14
2.2 Reconocimiento preliminar de las cuencas 15
2.3 Instrumentación hidrometeorológica 15
2.4 Modelo y análisis cartográfico 15
2.4.1 Generación de zonas homogéneas 16
2.4.2 Determinación en campo de coeficientes de escorrentía 16
2.4.3 Análisis hidrológico 18
2.4.4 Análisis hidráulico 19
2.5 Construcción de mapas de amenazas 19
2.6 Validación del modelo 20
2.7 Conclusiones y recomendaciones 20
3 Revisión de información secundaria 21
3.1 Revisión de material bibliográfico (Tesis y trabajos) 21
3.1.1 Biblioteca Central de la Universidad Nacional de Colombia, Sede Medellín 21
(Efe Gómez)
3.1.2 Centro de Documentación del Postgrado en Aprovechamiento de 22
Recursos Hidráulicos, Facultad de Minas
3.1.3 Biblioteca de la Facultad de Minas 22
3.1.4 Centro de Documentación del Instituto Mi Río 24
3.1.5 Biblioteca Central, Universidad De Antioquia 25
3.1.6 Biblioteca de las Empresas Públicas de Medellín 26
Universidad EAFIT i
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Página
3.1.7 Biblioteca de la Universidad EAFIT (Luis Echavarría Villegas) 27
3.2 Revisión e interpretación de imágenes de satélite y fotografías aéreas 27
3.3 Revisión de cartografía básica y temática 29
Revisión de datos hidrometeorológicos 31
4 Instrumentación 33
5 Análisis Cartográfico 37
5.1 Procesamiento de información base 37
5.1.1 Modelos de elevación digital 37
5.1.2 Generación de divisorias de aguas a partir de los modelos de elevación 41
digital
5.1.3 Generación de redes de drenaje a partir de los modelos de elevación 43
digital
5.1.4 Redes de drenaje con números de orden según Strahler y Shreve 45
5.1.5 Curvas hipsométricas de las cuencas 48
5.1.6 Análisis de pendientes 51
5.1.7 Parámetros morfométricos de las cuencas obtenidos a partir del MDE 54
5.2 Superposición de mapas para la generación de zonas homogéneas 57
5.3 Análisis Espacial 65
6 Determinación de coeficientes de escorrentía 72
6.1 Características del Simulador de Lluvias Portátil 74
6.1.1 Unidad de Alimentación y Control 74
6.1.2 Cabeza del Simulador 75
6.1.3 Trípode de soporte y pantalla contra el viento 76
6.1.4 Unidad de base 77
6.1.5 Sistema para medir la escorrentía 79
6.2 Procedimiento del ensayo de simulación de lluvias 80
6.3 Errores observados en los ensayos de simulación de lluvias 83
6.3.1 Errores debidos al funcionamiento del equipo 84
6.3.2 Errores de lectura 84
6.4 Localización de los ensayos de simulación de lluvias 84
6.5 Presentación, análisis y extrapolación de datos 86
6.5.1 Comportamiento del coeficiente de escorrentía durante un ensayo 86
6.5.2 Intensidad vs Coeficiente de Escorrentía 91
6.6 Criterios de extrapolación a otras zonas homogéneas 104
Universidad EAFIT ii
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Página
6.7 Manejo de resultados en el Sistema de Información Geográfica 104
7 Análisis Hidrológico 107
7.1 Descripción del modelo hidrológico aplicado 107
7.2 Cálculo del tiempo de concentración para las cuencas 109
7.3 Definición de lluvias 111
7.4 Caudales generados a partir del método racional distribuido 117
7.5 Cálculo de caudales a partir de la aplicación de modelos agregados 122
7.5.1 Método de las Ecuaciones de Regionalización de las Características 122
Medias de Smith et al. (1997)
7.5.2 Método Racional Probabilístico 126
7.5.3 Método de Índice de Crecientes 128
7.5.4 Hidrogramas Sintéticos, Modelo del Servicio de Conservación de Suelos 130
de los Estados Unidos
7.5.5 Hidrogramas sintéticos, Modelo de Williams y Hann 137
7.6 Comparación de caudales máximos determinados según diferentes 142
métodos
7.7 Comparación de caudales máximos con otras investigaciones 148
8 Análisis hidráulico 151
8.1 Datos geométricos 151
8.2 Datos de caudal 155
8.3 Datos de salida 156
9 Mapas preliminares de amenaza por inundación 160
10 Análisis del evento del 16 de septiembre de 1988 163
10.1 Descripción del evento 163
10.2 Resultados arrojados por el modelo para un evento similar al del 16 de 168
septiembre de 1988
10.3 Comparación de resultados con los registros del evento 171
11 Aforos en estaciones limnimétricas y limnigráficas 172
11.1 Lista de Estaciones 172
11.2 Observaciones de campo 172
11.3 Realización de Aforos 178
11.4 Resultados de los aforos 180
11.4.1 Altavista Almacén EEPPM 180
11.4.2 Picacha Conquistadores 182
11.4.3 Hueso Metro Floresta 183
Página
11.4.4 Ana Díaz Colegio Jesús María 183
11.4.5 Quintana Carrera 65 184
11.4.6 Iguaná Carrera 74 185
11.4.7 Santa Elena Plaza Minorista 186
12 Conclusiones y recomendaciones 189
13 Bibliografía 194
Anexos
Universidad EAFIT iv
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Lista de tablas
Página
Tabla 1. Cuencas de la zona suroccidental (Tomada de Mi Río, 1995). 7
Tabla 2. Cuencas de la zona noroccidental (Tomada de Mi Río, 1995). 11
Tabla 3. Relación de fotografías aéreas utilizadas en la interpretación. 28
Tabla 4. Estaciones del IDEAM con datos diarios de precipitación. 32
Tabla 5. Localización preliminar de estaciones en las cuencas de estudio. 33
Tabla 6. Localización de las estaciones hidrológicas y meteorológicas 35
instaladas por el SIMPAD en el Valle de Aburrá.
Tabla 7. Intervalos de agrupamiento natural para la clasificación de pendientes, 53
determinados por el método de Quiebres Naturales, para las cuencas del
estudio.
Tabla 8. Parámetros morfométricos de las cuencas obtenidos a partir de los 55
modelos de elevación digital.
Tabla 9. Índices de forma de las cuencas del estudio trabajadas hasta el 56
momento.
Tabla 10. Densidad de drenaje en km/km2 para las cuencas del estudio, con 57
base en la información obtenida de las redes de drenaje generadas a partir de
los MDEs.
Tabla 11. Codificación para el mapa de Formaciones Superficiales. 62
Tabla 12. Codificación para el mapa de Usos del Suelo. 62
Tabla 13. Codificación para el mapa de Pendientes por cuencas. 63
Tabla 14. Zonas homogéneas para la cuenca de la quebrada La Quintana, 66
número de celdas pertenecientes a cada zona, área y porcentaje de la misma
en la cuenca descartando la zona urbana, cuyos códigos comienzan por 4.
Tabla 15. Zonas seleccionadas para ensayos con simulador de lluvias en cada 70
cuenca, área cubierta por bosques y su porcentaje de área sobre el área rural
de la misma.
Tabla 16. Ejemplo de toma de datos en el ensayo de simulación de lluvias de 1 82
hora de duración.
Tabla 17. Ejemplo de resultados para los datos de la tabla 19. 82
Tabla 18. Lecturas de campo y resultados del ensayo 12. 87
Tabla 19. Lecturas de campo y resultados del ensayo 31. 88
Tabla 20. Lecturas de campo y resultados del ensayo 80. 89
Tabla 21. Lecturas de campo y resultados del ensayo 135. 90
Tabla 22. Ecuaciones de regresión del Coeficiente de escorrentía en función de 91
la Intensidad y comparación de los coeficientes de correlación para cada sitio.
Con amarillo se señala la ecuación de regresión escogida para cada sitio.
Universidad EAFIT v
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Página
Tabla 23. Tiempos de concentración en minutos para las cuencas estudio. 111
Tabla 24. Parámetros de las ecuaciones de las curvas IDF para la estaciones 112
ubicadas en la zona de estudio (tomado de comunicación escrita de EEPPM,
Area de Hidrometría e Instrumentación, abril de 2000).
Tabla 25. Intensidades de lluvia para cada cuenca con respecto a diversos 117
períodos de retorno.
Tabla 26. Caudales máximos instantáneos en los puntos más bajos de las 121
quebradas analizadas, para diversos períodos de retorno, obtenidos siguiendo
el método racional distribuido.
Tabla 27. Precipitación promedia multianual de las estaciones utilizadas en los 124
cálculos de caudales máximos de las quebradas del estudio (datos tomados de:
EEPPM, 1990. Anuario Hidrometeorológico).
Tabla 28. Precipitación media anual en cada cuenca estimada por el método de 124
los polígonos de Thiessen y área de influencia de cada estación en cada
cuenca.
Tabla 29. Parámetros morfométricos utilizados y características medias de los 125
caudales máximos instantáneos, según el método de las ecuaciones de
regionalización de Smith et al. (1997).
Tabla 30. Caudales máximos instantáneos calculados con las ecuaciones de 126
regionalización propuestas por Smith et al. (1997).
Tabla 31. Coeficientes de escorrentía probabilísticos para las cuencas de las 127
quebradas de estudio (tomados de las figuras 5.1 a 5.7 de Smith et al., 1997).
Tabla 32. Caudales máximos instantáneos calculados utilizando el método 128
racional agregado con coeficientes de escorrentía probabilísticos.
Tabla 33. Caudales Q2.33 en m3/s para los puntos más bajos de las cuencas de 129
drenaje consideradas en el estudio, considerando su área de drenaje en km2 y
utilizando la ecuación para la cuenca del río Medellín: Q2.33 = -0.00006 A2 +
0.3166 A + 10.789.
Tabla 34. Caudales máximos calculados en m3/s por el método del índice de 130
crecientes y la relación Q/Q2.33 = -0.4461 FX(X)2 + 1.0492 FX(X) + 0.6593, para
los puntos más bajos de las cuencas de drenaje consideradas en el estudio.
Tabla 35. Caudales máximos calculados en m3/s por el método del índice de 130
crecientes y la relación Q/A = -0.1648 FX(X)2 + 0.397 FX(X) + 0.1214, para los
puntos más bajos de las cuencas de drenaje consideradas en el estudio.
Tabla 36. Cálculo del CN para la cuenca de la quebrada La Hueso. 131
Tabla 37. CN para las cuencas de estudio. 132
Tabla 38. Parámetros para la construcción del hidrograma unitario sintético del 133
SCS.
Tabla 39. Cálculo de la precipitación efectiva según el método del SCS para 136
lluvias máximas en la cuenca de la quebrada La Hueso.
Universidad EAFIT vi
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Página
3
Tabla 40. Caudales máximos instantáneos en m /s obtenidos según el método 137
del Hidrograma Sintético, según el SCS.
Tabla 41. Parámetros y tiempos del hidrograma unitario sintético de Williams y 140
Hann calculados para la cuenca de la quebrada La Hueso.
Tabla 42. Caudales máximos para las cuencas del estudio determinados por el 141
método de Williams y Hann.
Tabla 43. Caudales máximos de 100 años de periodo de retorno para cada 142
quebrada, determinados según los diferentes métodos.
Tabla 44. Comparación de caudales máximos determinados por Duque y 149
Nigrinis (1992), Alvarado (1994) y en la presente investigación para la quebrada
Doña María.
Tabla 45. Comparación de caudales máximos determinados por Salazar y 149
Betancur (1991) y en la presente investigación para la quebrada La Picacha.
Tabla 46. Comparación de caudales máximos determinados por Salazar y 150
Betancur (1991) y en la presente investigación para la quebrada Ana Díaz.
Tabla 47. Comparación de caudales máximos determinados por Salazar y 150
Betancur (1991) y en la presente investigación para la quebrada La Hueso.
Tabla 48. Extracto de la tabla de valores generada por el SIG Arcinfo al hacer 152
las secciones transversales de la cuenca de la quebrada La Hueso (ST 01 HU).
Tabla 49. Altura máxima de agua para diferentes períodos de retorno para las 157
secciones transversales definidas en las quebradas Ana Díaz y La Picacha.
Datos arrojados por el programa HEC-RAS.
Tabla 50. Caudales máximos para una intensidad de 110 mm/hr, según el 168
método racional distribuido, para los puntos más bajos de las quebradas.
Tabla 51. Lista de estaciones para medición de nivel de agua en las quebradas, 172
instaladas por el SIMPAD a 27 de junio de 2001.
Lista de figuras
Página
Figura 1. Localización de la ciudad de Medellín en un contexto regional. 5
Figura 2. Localización de las cuencas seleccionadas para el estudio en el Valle 6
de Aburrá.
Figura 3. Esquema de las cuencas de la zona suroccidental en el Valle de 8
Aburrá.
Figura 4. Esquema de la cuenca de la quebrada La Iguaná. 9
Figura 5. Esquema de las cuencas de las quebradas Santa Elena y La Castro. 10
Figura 6. Esquema de las cuencas Minitas y La Quintana, de la zona 11
noroccidental del Valle de Aburrá.
Figura 7. Esquema de la cuenca Doña María, en el Corregimiento de San 13
Antonio de Prado.
Figura 8. Modelo cartográfico del estudio. 17
Figura 9. Distribución de las estaciones pluviográficas, limnimétricas y 36
limnigráficas instaladas por el SIMPAD.
Figura 10. Diagrama de flujo. Generación de modelos de elevación digital, 40
divisorias de agua y redes de drenaje.
Figura 11. Modelo de elevación digital para la cuenca de la quebrada La 41
Castro.
Figura 12. Dirección de flujo: el valor representado es el valor que toma la 42
celda del centro de acuerdo con la diferencia de alturas de las celdas vecinas.
Figura 13. Aplicación de la función flowdirection en un modelo de elevación 42
digital hipotético.
Figura 14. Aplicación de la función flowaccumulation en una matriz de dirección 44
de flujo.
Figura 15. Número de orden según Strahler. 46
Figura 16. Número de orden según Shreve. 46
Figura 17. Números de orden de la red de drenaje de la cuenca de la quebrada 47
La Castro, según Strahler.
Figura 18. Números de orden de la red de drenaje de la cuenca de la quebrada 48
La Castro, según Shreve.
Figura 19. Curva hipsométrica de la cuenca de la quebrada La Castro. 49
Figura 20. Cálculo de la pendiente para cada celda. 51
Figura 21. Histograma de frecuencias o distribución natural de las pendientes 52
en la cuenca de la quebrada La Castro.
Figura 22. Mapa de pendientes clasificadas de la cuenca de la quebrada La 54
Castro.
Página
Figura 23. Mapa de usos del suelo de la cuenca de la quebrada La Castro. 59
Figura 24. Mapa de Formaciones Superficiales de la cuenca de la quebrada La 61
Castro.
Figura 25. Esquema del proceso de rasterización. 63
Figura 26. Esquema de superposición de los mapas en formato raster, 64
mediante una suma sencilla de los valores correspondientes a cada celda para
la obtención del mapa de zonas homogéneas con respecto al coeficiente de
escorrentía.
Figura 27. Mapa de Zonas Homogéneas de un sector de la cuenca de la 65
quebrada Santa Elena. El tamaño de celda es de 5 metros de lado.
Figura 28. Histograma de frecuencias para el porcentaje de área de las zonas 67
homogéneas en la cuenca de la quebrada La Picacha.
Figura 29. Histograma de frecuencias para el porcentaje de área de las zonas 67
homogéneas en la cuenca de la quebrada Altavista.
Figura 30. Histograma de frecuencias para el porcentaje de área de las zonas 67
homogéneas en la cuenca de la quebrada Ana Díaz.
Figura 31. Histograma de frecuencias para el porcentaje de área de las zonas 68
homogéneas en la cuenca de la quebrada La Hueso.
Figura 32. Histograma de frecuencias para el porcentaje de área de las zonas 68
homogéneas en la cuenca de la quebrada Minitas.
Figura 33. Histograma de frecuencias para el porcentaje de área de las zonas 68
homogéneas en la cuenca de la quebrada La Quintana.
Figura 34. Histograma de frecuencias para el porcentaje de área de las zonas 69
homogéneas en la cuenca de la quebrada Santa Elena.
Figura 35. Histograma de frecuencias para el porcentaje de área de las zonas 69
homogéneas en la cuenca de la quebrada La Iguaná.
Figura 36. Histograma de frecuencias para el porcentaje de área de las zonas 69
homogéneas en la cuenca de la quebrada Doña María.
Figura 37. Simulador de lluvias portátil. Diseño inicial publicado por McQueen 73
en 1963 (tomado de McQueen, 1963, U.S. Geological Survey Circular 482).
Figura 38. Simulador de lluvias portátil construido en la Universidad EAFIT 74
(primera versión).
Figura 39. Unidad de alimentación y control del simulador de lluvias portátil. 75
Figura 40. Cabeza del simulador de lluvias portátil. 76
Figura 41. Trípode de soporte y pantalla contra el viento del simulador de 77
lluvias portátil.
Figura 42. Detalle de la unidad de base cilíndrica del simulador de lluvias 78
portátil.
Figura 43. Acople de la unidad base al resto del sistema, con el embudo para 79
evitar pérdidas por desviación de gotas.
Universidad EAFIT ix
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Página
Figura 44. Sistema para medir la escorrentía. 80
Figura 45. Número de ensayos de simulación de lluvias efectuados en cada 85
sitio.
Figura 46. Localización de los sitios en los que se efectuaron ensayos con el 85
simulador de lluvias.
Figura 47. Cambio del coeficiente de escorrentía con respecto al tiempo 87
durante el ensayo 12.
Figura 48. Cambio del coeficiente de escorrentía con respecto al tiempo 88
durante el ensayo 31.
Figura 49. Cambio del coeficiente de escorrentía con respecto al tiempo 89
durante el ensayo 80.
Figura 50. Cambio del coeficiente de escorrentía con respecto al tiempo 90
durante el ensayo 135.
Figura 51. Relación entre la Intensidad de la lluvia y el Coeficiente de 92
Escorrentía en el sitio 01.
Figura 52. Relación entre la Intensidad de la lluvia y el Coeficiente de 93
Escorrentía en el sitio 02.
Figura 53. Relación entre la Intensidad de la lluvia y el Coeficiente de 93
Escorrentía en el sitio 03.
Figura 54. Relación entre la Intensidad de la lluvia y el Coeficiente de 94
Escorrentía en el sitio 04.
Figura 55. Relación entre la Intensidad de la lluvia y el Coeficiente de 94
Escorrentía en el sitio 07.
Figura 56. Relación entre la Intensidad de la lluvia y el Coeficiente de 95
Escorrentía en el sitio 08.
Figura 57. Relación entre la Intensidad de la lluvia y el Coeficiente de 95
Escorrentía en el sitio 09.
Figura 58. Relación entre la Intensidad de la lluvia y el Coeficiente de 96
Escorrentía en el sitio 10.
Figura 59. Relación entre la Intensidad de la lluvia y el Coeficiente de 96
Escorrentía en el sitio 11.
Figura 60. Relación entre la Intensidad de la lluvia y el Coeficiente de 97
Escorrentía en el sitio 12.
Figura 61. Relación entre la Intensidad de la lluvia y el Coeficiente de 97
Escorrentía en el sitio 13.
Figura 62. Relación entre la Intensidad de la lluvia y el Coeficiente de 98
Escorrentía en el sitio 14.
Figura 63. Relación entre la Intensidad de la lluvia y el Coeficiente de 98
Escorrentía en el sitio 15.
Universidad EAFIT x
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Página
Figura 64. Relación entre la Intensidad de la lluvia y el Coeficiente de 99
Escorrentía en el sitio 16.
Figura 65. Relación entre la Intensidad de la lluvia y el Coeficiente de 99
Escorrentía en el sitio 17.
Figura 66. Relación entre la Intensidad de la lluvia y el Coeficiente de 100
Escorrentía en el sitio 18.
Figura 67. Relación entre la Intensidad de la lluvia y el Coeficiente de 100
Escorrentía en el sitio 19.
Figura 68. Relación entre la Intensidad de la lluvia y el Coeficiente de 101
Escorrentía en el sitio 20.
Figura 69. Relación entre la Intensidad de la lluvia y el Coeficiente de 101
Escorrentía en el sitio 21.
Figura 70. Relación entre la Intensidad de la lluvia y el Coeficiente de 102
Escorrentía en el sitio 22.
Figura 71. Relación entre la Intensidad de la lluvia y el Coeficiente de 102
Escorrentía en el sitio 23.
Figura 72. Relación entre la Intensidad de la lluvia y el Coeficiente de 103
Escorrentía en el sitio 24.
Figura 73. Relación entre la Intensidad de la lluvia y el Coeficiente de 103
Escorrentía en el sitio 25.
Figura 74. Mapa de zonas generadoras de escorrentía para la cuenca de la 106
quebrada La Castro.
Figura 75. Curvas I D F para la estación San Cristóbal. 113
Figura 76. Curvas I D F para la estación Villa Hermosa. 113
Figura 77. Curvas I D F para la estación San Antonio de Prado. 114
Figura 78. Curvas I D F para la estación Miguel Aguinaga. 114
Figura 79. Curvas I D F para la estación Aeropuerto Olaya Herrera. 115
Figura 80. Curvas I D F para la estación Vasconia. 115
Figura 81. Ubicación de estaciones con curvas IDF. 117
Figura 82. Localización de secciones transversales en la quebrada Doña María. 118
Figura 83. Localización de secciones transversales en las quebradas Altavista 119
y La Picacha.
Figura 84. Localización de secciones transversales en las quebradas Ana Díaz 119
y La Hueso.
Figura 85. Localización de secciones transversales en la quebrada La Iguaná. 120
Figura 86. Localización de secciones transversales en las quebradas Minitas y 120
La Quintana.
Figura 87. Localización de secciones transversales en las quebradas Santa 121
Elena y La Castro.
Universidad EAFIT xi
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Página
Figura 88. Polígonos de Thiessen formados por las estaciones con 125
precipitación media anual y su área de influencia en cada una de las cuencas.
Figura 89. Hidrograma Unitario para T igual a la duración de la lluvia efectiva. 134
Cuenca de la quebrada La Hueso.
Figura 90. Curvas S para la cuenca de la quebrada La Hueso para lluvias de 135
duración igual a la lluvia efectiva (0.132 horas) y tiempo de concentración.
Figura 91. Hidrograma unitario (SCS) para la cuenca de la quebrada La Hueso, 135
para una lluvia de duración 1.19 horas (igual al tiempo de concentración).
Figura 92. Hidrograma de escorrentía directa para la cuenca de la quebrada La 137
Hueso producido por una lluvia con Tr = 100 años y duración igual al tiempo de
concentración.
Figura 93. Hidrograma unitario sintético de Williams y Hann (tomada de Smith 138
et al., 1997, figura 6.4, página 98).
Figura 94. Relación entre B y n del hidrograma unitario de Williams y Hann 139
(tomada de Smith et al., 1997, figura 6.4, página 100).
Figura 95. Hidrograma Unitario Sintético de Williams y Hann para la cuenca de 141
la quebrada La Hueso y de duración igual al tiempo de concentración de la
cuenca, el círculo rojo indica el punto de caudal máximo.
Figura 96. Comparación de caudales máximos para la quebrada Doña María. 143
Figura 97. Comparación de caudales máximos para la quebrada Altavista. 143
Figura 98. Comparación de caudales máximos para la quebrada La Picacha. 144
Figura 99. Comparación de caudales máximos para la quebrada Ana Díaz. 144
Figura 100. Comparación de caudales máximos para la quebrada La Hueso. 145
Figura 101. Comparación de caudales máximos para la quebrada La Iguaná. 145
Figura 102. Comparación de caudales máximos para la quebrada Minitas. 146
Figura 103. Comparación de caudales máximos para la quebrada La Quintana. 146
Figura 104. Comparación de caudales máximos para la quebrada Santa Elena. 147
Figura 105. Comparación de caudales máximos para la quebrada La Castro. 147
Figura 106. Sección transversal número 01 sobre la quebrada La Hueso, 153
gráfica de MS Excel 97.
Figura 107. Fotografía de la sección transversal ST 01 HU con referencias de 154
distancia.
Figura 108. Sección transversal tomada del DEM (Cauce corregido a partir de 154
la foto), ST 01 HU.
Figura 109. Datos geométricos de la sección transversal ST 01 HU. 155
Página
Figura 111. Sección transversal ST 01 HU de la quebrada La Hueso con los 158
niveles alcanzados debido a caudales máximos instantáneos (líneas
horizontales de color violeta).
Figura 112. Sección longitudinal de la quebrada La Hueso con los niveles 158
normales y los niveles alcanzados debido a caudales máximos instantáneos.
Figura 113. Sección longitudinal de la quebrada La Hueso entre las secciones 159
01 y 04, vista en perspectiva, con los niveles normales y los niveles alcanzados
debido a caudales máximos instantáneos.
Figura 114. Curva de gasto (caudal vs elevación del nivel de agua) para la 159
sección transversal ST 01 HU de la quebrada La Hueso.
Figura 115. Proceso de construcción de las superficies de inundación a partir 160
de los datos de nivel de agua en cada sección transversal determinados en el
programa HEC – RAS.
Figura 116. Resta de matrices para determinar las zonas de inundación. Las 161
celdas con valores negativos corresponden a zonas secas y las celdas con
valores positivos corresponden a zonas que son susceptibles de ser inundadas
en los eventos de caudales máximos.
Figura 117. Zonas de amenaza por inundación de un tramo de la quebrada La 162
Castro (una cuadra arriba de la calle 52). La matriz resultante del proceso está
dibujada con color rojo, la zona de inundación digitalizada se presenta en color
amarillo y corresponde a la zona de amenaza por inundación en este sector.
Figura 118. Fotocopia del pluviograma de la estación San Cristóbal del viernes 163
16 al domingo 19 de septiembre de 1988 (datos suministrados por EEPPM).
Figura 119. Quebrada La Iguaná en el sector de Los Colores. 164
Figura 120. Barrio La Iguaná, aguas arriba de la carrera 65. 164
Figura 121. Barrio La Iguaná en la margen izquierda de la quebrada del mismo 165
nombre.
Figura 122. Sector Plaza de Toros La Macarena, vista hacia el occidente. 166
Figura 123. Puente de la calle 35 con carrera 89 sobre la quebrada Ana Díaz. 166
Figura 124. Garajes en el barrio Laureles que quedaron completamente 167
anegados.
Figura 125. Avenida 33 con calle 65, desborde la quebrada La Picacha. 167
Figura 126. Quebrada La Picacha en el cruce con la avenida 33, zona de 169
inundación, en rojo, para una lluvia de 110 mm/hr y duración 1 hora.
Figura 127. Zona de inundación en el sector Naranjal – Suramericana, en rojo, 170
asociada a la quebrada La Hueso y para una lluvia de 110 mm/hr y duración 1
hora.
Figura 128. Zona de inundación en el sector Estadio, en rojo, asociada a la 171
quebrada La Hueso y para una lluvia de 110 mm/hr y duración 1 hora.
Página
Figura 129. Localización de las estaciones limnimétricas y limnigráficas en el 173
contexto del Valle de Aburrá.
Figura 130. Fotografía del limnímetro de la estación Iguaná Ladrillera San 174
Cristóbal. En ésta puede observarse la pequeña caída de agua que se formó
hacia la parte inferior de la foto y las rocas y el pequeño empozamiento hacia la
parte central y en el borde de la quebrada, evidencias de la inestabilidad de la
sección.
Figura 131. Estación Hueso Metro Suramericana. Obsérvese cómo se modifica 175
la geometría de la sección al acumular las piedras al lado izquierdo del canal,
justo donde se encuentra el limnímetro (adentro del círculo azul).
Figura 132. Vista del limnímetro de la estación Hueso Metro Suramericana. 175
Obsérvese la diferencia entre el nivel del agua en el limnímetro y en el centro
del canal (parte inferior de la foto), en este caso la lectura del limnímetro no
corresponde al nivel real del agua.
Figura 133. Sitio adecuado para la reubicación de la estación Hueso Metro 176
Suramericana, en el sitio indicado por la flecha roja, a aproximadamente 20
metros aguas abajo de la confluencia de las quebradas Hueso y Ana Díaz.
Figura 134. Estación Doña María Vereda La Verde. Las piedras arrastradas por 176
la quebrada formaron una pequeña caída de agua en el sector donde se ubicó
el limnímetro (círculo azul).
Figura 135. Detalle del limnímetro de la estación Doña María Vereda La Verde. 177
Obsérvese la diferencia en el nivel de agua de la quebrada en la zona inferior
de la foto y el señalado en el limnímetro.
Figura 136. Estación Castro Hogar Antioquia. El limnímetro (ver círculo azul) 177
está ubicado justo debajo del puente peatonal, a pocos centímetros de las
piedras que se observan en la foto y que hacen que el flujo de agua sea nulo en
ese lado del canal.
Figura 137. Ubicación de los jalones en una línea perpendicular al cauce, con 178
la cuerda amarrada entre ellos. Estación Iguaná Carrera 74, sitio con limnímetro
y limnígrafo.
Figura 138. Detalle de la cuerda (color amarillo) amarrada entre los jalones con 179
marcas cada 50 cm (negras, ver recuadro azul). A la izquierda puede
observarse el sitio donde está ubicado el limnímetro (barra roja).
Figura 139. Esquema de una sección transversal, dibujada con datos de 179
profundidad cada 50 cm (h1, h2) y los puntos donde se toma el dato de
velocidad (v).
Figura 140. Esquema de medición de la velocidad de la corriente en canales en 180
“V” angostos, x1, x2 y x3 son las distancias a las cuales se mide la profundidad
del cauce, h1 y h2.
Figura 141. Medición de la velocidad de la corriente en el punto medio de los 180
tramos señalados.
Figura 142. Estación Altavista Almacén EEPPM. 181
Página
Figura 143. Estación Altavista Almacén EEPPM, vista desde el puente 181
peatonal.
Figura 144. Estación Picacha Conquistadores, vista desde aguas arriba, 182
margen izquierda.
Figura 145. Estación Picacha Conquistadores, vista desde aguas arriba, 182
margen derecha.
Figura 146. Estación Hueso Metro Floresta. 183
Figura 147. Estación Ana Díaz Colegio Jesús María, visto desde la margen 184
derecha.
Figura 148. Estación Ana Díaz Colegio Jesús María, al fondo se observa el 184
puente peatonal.
Figura 149. Estación Quintana Carrera 65, vista desde la margen izquierda. 185
Figura 150. Detalle de la estación Quintana Carrera 65. 185
Figura 151. Estación Iguaná Carrera 74, al fondo se observa el puente de la 186
carrera 74.
Figura 152. Detalle de la realización del aforo en la estación Iguaná Carrera 74, 186
la medición de la velocidad del agua se hizo cada 50 cm.
Figura 153. Estación Santa Elena Plaza Minorista. El limnímetro se encuentra 187
en la zona demarcada por el círculo azul.
Figura 154. Detalle del limnímetro de la estación Santa Elena Plaza Minorista. 188
Universidad EAFIT xv
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
1 Introducción
De cada una de las cuencas de estudio se presenta un anexo que contiene los mapas
obtenidos tanto en pasos intermedios de la investigación como los resultados finales, esto
es, los mapas preliminares de amenazas, detallando a una escala adecuada las zonas
que poseen amenaza por inundación, superpuesta con la información de viviendas e
infraestructura ubicada en dichas zonas.
Dicha ocupación del territorio se ha llevado a cabo por falta de otras alternativas, en zonas
que corresponden a áreas de inundación de caños y quebradas y en zonas de laderas
Universidad EAFIT 1
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
con fuertes pendientes en las partes altas de las vertientes del Valle de Aburrá. Esta
ocupación desordenada y sin ningún control ha producido un aumento de la vulnerabilidad
de la ciudad ante eventos naturales que pueden calificarse como catastróficos o
desastres, ya que su ocurrencia es responsable de la destrucción de los bienes y servicios
de la comunidad y en algunos casos hasta de la muerte de personas.
El caso de las quebradas del Valle de Aburrá cobra importancia debido al comportamiento
torrencial de las mismas, que históricamente se ha reflejado en la ocurrencia de
crecientes súbitas. Una creciente súbita se presenta cuando el tiempo transcurrido entre
la lluvia y el aumento de nivel en el cauce es mínimo. Este problema ya ha sido tratado
anteriormente, sin que hasta ahora se hayan efectuado labores tendientes a dar una
respuesta definitiva. Tal es el caso del planteamiento efectuado por HIMAT (1990): “En el
caso particular de Medellín, la característica común de todas las quebradas es que una
lluvia intensa puede provocar una crecida repentina. Tales crecidas se desarrollan tan
rápidamente que constituyen un fenómeno de sorpresa para los ocupantes de la zona de
inundación. En el mejor de los casos, el tiempo transcurrido entre la presencia de una
lluvia intensa, y el nivel máximo de crecida, puede ser de media hora como máximo.”
Universidad EAFIT 2
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Dentro del proyecto mencionado se han definido diez quebradas del Valle de Aburrá para
elaborar un estudio piloto que parta de las características físicas del terreno y de la
información hidrometeorológica disponible y permita instrumentar las cuencas de dichas
quebradas con el fin de implementar, ajustar y evaluar una metodología adecuada para la
puesta en marcha y operación de sistemas de alarma temprana y su posterior ampliación
al resto de quebradas del valle.
Las quebradas seleccionadas para este estudio, que poseen mayor cantidad de reportes
de emergencia ante eventos de inundación (recopilados por el SIMPAD), han sido
clasificadas inicialmente según el área de la cuenca (Mi Río, Municipio de Medellín,
PNUD, 1995): Santa Elena (y su afluente La Castro), La Iguaná y Doña María como
Cuencas Superiores o Tipo A (mayores de 20 Km2); Ana Díaz, La Hueso, La Picacha y
Alta Vista como Cuencas Mayores o Tipo B (entre 5 y 20 Km2); y La Minitas y La Quintana
como Cuencas Medianas o Tipo C (entre 1 y 5 Km2).
Universidad EAFIT 3
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 4
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
un estrecho valle encerrado en altas montañas que definen la cuenca del Valle del Río
Medellín, valle que comparte con otros nueve municipios.
Universidad EAFIT 5
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
La totalidad de las quebradas consideradas en el presente estudio son afluentes del río
Medellín (ver figura 2). El recuento de las quebradas se basa en la clasificación por zonas
geográficas realizada por el Instituto Mi Río y el Municipio de Medellín, 1995.
Universidad EAFIT 6
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
La cuenca de la quebrada Altavista, limita por el norte con la cuenca de La Picacha, por el
sur con la cuenca de La Guayabala, por el occidente con la cuenca de la quebrada Doña
María y por el oriente con el río Medellín. Posee en su cabecera dos ramales: quebrada
Barcino, que nace en las estribaciones del Alto El Encanto a 2340 msnm, con una longitud
de 2.7 Km y quebrada Buga, con nacimiento en el cerro El Barcino en la cota 2270 msnm
con una longitud de 2.4 Km; ambas quebradas confluyen en la cota 1720 msnm donde
toman el nombre de quebrada Altavista, y luego de 7.2 Km tributa sus aguas al río
Medellín en la cota 1475 msnm, al pie del cerro Nutibara. Son peculiares los drenajes
largos en su margen derecha, con un patrón de drenaje subparalelo dendrítico controlado
por el diaclasamiento con sentido NS y NE. En su margen izquierda, la quebrada Santa
Rita es el único drenaje largo.
Limita al norte con la cuenca de la quebrada Malpaso, al noroccidente con la Cuchilla Las
Baldías, divisoria de aguas con la cuenca de la quebrada El Hato, al occidente con la
divisoria de aguas de la cordillera El Frisol y la quebrada La Frisola hasta el Cerro del
Padre Amaya, al suroccidente con el Alto El Astillero y la Cuchilla Los Arrayanes, por el
sur con la cuenca de la quebrada La Hueso y por el oriente con el río Medellín. Nace en
la cota 2555 msnm en las estribaciones del Alto de las Repetidoras, después de recorrer
en dirección SE un tramo aproximado de 5.6 Km cambia a dirección WE hasta su
desembocadura en el río Medellín, en la cota 1457 msnm. Posee una longitud total de
15.29 Km y un área de 51.22 Km2, de los cuales 6.35 Km2 corresponden a la zona urbana
y 44.87 Km2 a la zona rural. La figura 4 muestra un esquema de la cuenca de la quebrada
La Iguaná.
Universidad EAFIT 9
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Limita al norte con las cuencas de las quebradas El Ahorcado, El Molino y parte de
Piedras Blancas, por el oriente con la cuchilla La Gulupera, al Sur con el municipio de
Envigado y las cuencas de las quebradas La Presidente y La Poblada y al occidente con
el río Medellín. El área total de la cuenca es de 45.61 Km2 de los cuales el 74.95%
pertenece a zona rural y el 25.05% a zona urbana, el estudio hidrológico de la cuenca se
efectuó hasta la entrada del box-coulvert, unos metros arriba del teatro Pablo Tobón
Uribe, con un área de captación de 39.46 Km2. Nace en el Cerro Espíritu Santo en la cota
2620 msnm y luego de un recorrido de 14.43 Km desemboca en el río Medellín en la cota
1457 msnm. La longitud del cauce hasta la entrada del box-coulvert es de 12.08 Km a
una altura de 1505 msnm. La figura 5 muestra un esquema de la cuenca de la quebrada
Santa Elena.
Universidad EAFIT 10
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Este afluente de la quebrada Santa Elena presenta una longitud de 4.24 Km. Su cuenca
hace parte de la vereda Piedras Blancas y los barrios Las Estancias y Villa Tina,
comprendiendo un área de 4.22 Km2, de los cuales 4.06 Km2 pertenecen a zona rural y
0.16 Km2 a zona urbana. Nace en la cota 2500 msnm y desemboca en la cota 1630
msnm.
Esta zona comprende las cuencas de las quebradas localizadas en la parte noroccidental
de Medellín, en las comunas 5,6 y 7 limitada por el norte con la cuenca de la quebrada La
Loca del municipio de Bello, por el sur y suroccidente con la cuenca de la quebrada La
Iguaná, por el occidente con el Alto de Yolombó y por el oriente con el río Medellín. El
área total de la zona corresponde a 19.5 Km2 subdividida en quince cuencas, dos de las
cuales son objeto de estudio: La Quintana y Minitas, esquematizadas en la figura 6. En la
tabla 2 se reseñan las principales características de las cuencas en mención.
Universidad EAFIT 11
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Limita al sur y occidente con las cuencas de las quebradas Malpaso y El Malpaso, la
cuenca de la quebrada La Cantera al norte y el río Medellín por el oriente. Nace en la
cota 2207 msnm y desemboca en la cota 1445 msnm. La cuenca se caracteriza por
poseer un patrón de drenaje subparalelo con tendencia noreste y su forma rectangular –
oblonga.
Limita por el norte con la cuenca de la quebrada San Francisco, por el oriente con las
cuencas de las quebradas La Hueso, La Picacha, Altavista, Jabalcona, por el sur con la
cuenca de la quebrada Grande y por el occidente con las cuencas de las quebradas la
Frisola y La Sucia. Hacia el sur del Cerro del Padre Amaya nacen dos afluentes de la
quebrada Doña María: quebrada La Frisola en la cota 2905 msnm y la quebrada Los
Vallados en la cota 2590 msnm y a partir de la confluencia de estas en la cota 2240
msnm, toma el nombre de quebrada Doña María. Desemboca en el río Medellín en la
cota 1522 msnm luego de un recorrido de 22.1 Km. El área total de la cuenca es de 72.6
Km2, con un total de 48 microcuencas en la margen derecha y 43 microcuencas en la
margen izquierda. La figura 7 muestra el mapa de la cuenca de la quebrada Doña María.
Universidad EAFIT 12
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 13
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
2 Metodología
Se revisaron varias tesis y trabajos que abordaban la misma temática del trabajo, sirvieron
inicialmente para reconocer las características propias de cada cuenca y posteriormente
para comparar los resultados obtenidos en la investigación actual. De cada trabajo se
reseñaron sus planteamientos y enfoques.
Universidad EAFIT 14
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Se revisaron e interpretaron las fotografías aéreas del área del Valle de Aburrá, tomadas
en 1998, en escalas aproximadas 1:10.000 para el área urbana y 1:20.000 para el área
rural y se revisó la imagen de satélite de agosto de 1996 del Valle de Aburrá disponible en
el Laboratorio de Imágenes de la universidad EAFIT. Las fotos aéreas se utilizaron
inicialmente para revisar las unidades geomorfológicas y formaciones superficiales y
familiarizarse con cada una de las cuencas y posteriormente para definir los usos del
suelo de las cuencas.
El siguiente paso consistió en diseñar la metodología con que se iba a abordar los análisis
en cada una de las cuencas. Ubicadas las fuentes de información y teniendo en cuenta la
Universidad EAFIT 15
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Esta etapa de campo tuvo como objeto final determinar la relación entre la intensidad de
la lluvia y el coeficiente de escorrentía para zonas que se supusieron con comportamiento
hidrológico similar. Para lograr esto se construyó y utilizó un simulador de lluvias portátil
siguiendo el diseño propuesto por McQueen (USGS, 1963). Durante el período de
campaña se hizo simulación de lluvias en 23 sitios diferentes, de 4 a 11 ensayos por sitio,
para un total de 159 ensayos.
Universidad EAFIT 16
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Simulación
de lluvias 2
geoproce rast, reclas
usc_cuenca CAMPO
Funciones
3 Flujo 2
Rast,inter (*) llr_cuenca zhCe_cuenca
pg_cuenca lluvia real
Conocida la relación entre una lluvia y la escorrentía que se genera para cada celda que
compone una zona homogénea en particular y teniendo para cada cuenca el modelo de
elevación digital, se aplicaron funciones de flujo propias del SIG ArcInfo para determinar
primero la dirección de flujo para cada celda del modelo y luego la acumulación de flujo en
cada celda. Si se supone una lluvia con intensidad constante durante un tiempo igual al
tiempo de concentración de la cuenca, la matriz de acumulación de flujo representa para
cada celda el caudal generado por dicha lluvia. En la investigación actual se tuvieron en
cuenta las lluvias extremas de acuerdo con curvas de Intensidad – Duración – Frecuencia
de estaciones con registros históricos en la zona de estudio, definiendo una estación foco
para cada cuenca y suponiendo que la lluvia presentaba una distribución espacial
homogénea.
El cálculo del caudal en cada punto de la cuenca surgido a partir de una lluvia en
particular se determina con la aplicación de dos funciones propias de ArcInfo llamadas
flowdirection y flowaccumulation ejecutadas en el modelo digital de elevación de la cuenca
(dem_cuenca). La función flowdirection permite definir la dirección de flujo para cada
celda del modelo digital de elevación, mientras que la función flowaccumulation calcula el
caudal acumulado en cada celda de la cuenca determinado por la suma de las celdas que
fluyen hacia esta.
Q= C * I * A
Donde Q: Caudal en el punto de captación
C: Coeficiente de Escorrentía
I: Intensidad de lluvia
A: Área de captación
El método utilizado se denominó “método racional distribuido”, esto es, a cada celda se le
aplica el método racional y se suman los aportes de caudal de cada celda para encontrar
el caudal en cada celda de la matriz.
Obtenidos los caudales máximos instantáneos generados por una lluvia se procedió a
compararlos con caudales arrojados por otros métodos ampliamente utilizados en nuestro
medio. Además se compararon los valores obtenidos en la investigación con los
resultados publicados en trabajos anteriores. Ya que se supone una lluvia con distribución
espacial y temporal homogénea, y que se determinaron los caudales en cada celda
utilizando un modelo lineal sencillo, el modelo hidrológico utilizado es apenas una
aproximación hacia la determinación de caudales máximos y en ningún caso debe verse
Universidad EAFIT 18
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Las secciones transversales a los cauces de las quebradas se tomaron a partir de los
modelos de elevación digital con la utilización del SIG ArcInfo. En la mayor parte de las
secciones no se detallaba la forma exacta del cauce, por lo cual fue necesario corregir la
forma del cauce a partir de medidas efectuadas en campo y de fotos tomadas
perpendiculares a las secciones y con referencias de distancia que posibilitaban la
corrección de la forma.
Además de la corrección de la forma del canal, las fotos tomadas en cada una de las
secciones fueron útiles para la determinación de los números de Manning que
caracterizan la rugosidad de las llanuras de inundación y el canal en las secciones. Estos
valores se definieron a partir de tablas que relacionan las características de la sección con
coeficientes de rugosidad. Otros datos que el programa requiere son las distancias entre
secciones tanto por el canal como por las bancas y las abscisas que definen el canal para
cada sección transversal.
Del modelamiento hidráulico se obtienen los niveles de agua producidos por un caudal
máximo derivado de una lluvia crítica para cada sección transversal. Para definir zonas
de amenazas es necesario construir una superficie de inundación partiendo de perfiles
con un nivel de agua particular y de un modelo de elevación digital. Se interpolan los
niveles de agua de los perfiles, obteniendo una superficie que representa el nivel del agua
que produce una lluvia. A dicha superficie se le resta el modelo de elevación digital,
obteniendo una nueva superficie en la cual los valores positivos representan zonas de
inundación y los valores negativos zonas hasta donde no llega el agua generada por una
Universidad EAFIT 19
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
creciente. Por último se comparan las zonas de inundación con las observaciones
efectuadas en campo para ajustar las zonas de amenazas y al superponer la cobertura de
vías y viviendas se determinan zonas de riesgo.
Como no se disponen de datos medidos de lluvias y caudal no fue posible efectuar una
validación sistemática del modelo. Con el modelamiento de una lluvia extrema que
ocasionó desbordamientos en las quebradas de la zona occidental del Valle de Aburrá en
septiembre de 1988 y de la cual se cuenta con un registro fotográfico tomado de la prensa
local, fue posible una validación inicial del modelo seguido en la presente investigación.
Universidad EAFIT 20
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
De los perfiles descritos en este estudio sólo dos están dentro del área de las cuencas de
las quebradas del presente proyecto, uno dentro de la cuenca de la quebrada Doña
María, en la vía de San Antonio de Prado a Heliconia, a 1.8 Km del primero; y otro en la
cuenca de la quebrada La Iguaná, en el sitio Boquerón. Los perfiles pertenecen a dos
conjuntos diferentes (Abejorrál y Boquerón, respectivamente), dentro de la misma
asociación (Tequendamita).
Universidad EAFIT 21
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Salazar, Jose & Betancur Jaime, 1991. Estudio hidrológico e hidráulico de las
quebradas La Hueso, La Ana Díaz, La Picacha. Tesis de grado. Universidad
Nacional de Colombia, Sede Medellín, Facultad de Minas.
Se compararon los caudales máximos obtenidos en dicho trabajo con los arrojados en la
actual investigación.
Angel, Jorge & Duque Teresa, 1990. Zonificación de riesgo para la cuenca de la
Quebrada Iguaná. Tesis de grado. Universidad Nacional de Colombia, Sede
Medellín, Facultad de Minas.
T551.3/A54
Universidad EAFIT 23
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Smith, R.; Vélez, M. V.; Acosta, Z. L.; Blandón, J. R.; Posada, J.E.; Botero, A.L.;
Campo, J.M.; Montoya, R.D.; Múnera, J.C., Muñoz, S.I.; & Ruiz, C.D., 1997.
Hidrología de Antioquia. Ed. Imprenta Departamental de Antioquia. 176 p, anexos.
De los temas que fueron útiles en el trabajo actual, se encuentran los capítulos de
geomorfología y de hidrología. En este último se realiza un estudio de niveles de
inundación y caudales de diseño de estructuras por el método de regionalización.
Además presenta resultados de caudales de un estudio desarrollado por ACP Ltda,
utilizando hidrogramas unitarios sintéticos. Estos resultados se compararon con los
obtenidos en la actual investigación.
Universidad EAFIT 24
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 25
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Se adoptó una frecuencia de 100 años para determinar el caudal máximo teniendo en
cuenta obras hidráulicas de regulación y canalización. Dicho análisis se efectuó a partir
de datos pluviográficos de las estaciones de EPM en Caldas, San Cristóbal, San Antonio
de Prado, Planta de Filtros, Chorrillos, Mazo y Las Palmas (registros desde 1948). Se
construyó una hidrógrafa unitaria sintética para la cuenca del río Medellín, puesto que no
se tenían datos de limnígrafos.
Aguilar, Parra & Sañudo, Jairo, 1997. Plan de ajustes para el ordenamiento
ambiental de la cuenca de la quebrada La Iguaná (desde su nacimiento hasta el
sitio de captación por EEPPMM). Tesis de especialización en Ingeniería Ambiental,
Universidad Pontificia Bolivariana, Medellín.
333.91715 / A 283
Este libro contiene información climática e hidrológica de las estaciones de las Empresas
Públicas de Medellín ubicadas en el Valle de Aburrá y sus alrededores. Se ha utilizado
para el análisis de información secundaria para encontrar caudales máximos instantáneos
en algunas de las quebradas de la zona de estudio.
Contiene información básica de cada una de las quebradas de estudio tales como área,
perímetro, longitud del cauce, pendientes, etc. Además se reseñan las características
geológicas, geomorfológicas, de suelos y sus usos, entre otras. Este trabajo fue útil en la
etapa inicial de la actual investigación para reconocer las cuencas de estudio.
Se revisaron e interpretaron las fotografías aéreas del área del Valle de Aburrá, tomadas
en 1998, en escalas aproximadas 1:10.000 para el área urbana y 1:20.000 para el área
Universidad EAFIT 27
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
digital, se decidió trabajar con esta última información en lugar de trabajar con la imagen
de satélite antes mencionada.
Con respecto a los usos del suelo, luego de acceder al mapa suministrado por Planeación
Municipal y encontrar que la información cartográfica no se pudo procesar, se averiguó la
procedencia de dicho mapa y se encontró que éste se había levantado en el año de 1987
con una verificación en campo tendiente a la ubicación de zonas de protección para el
Plan de Ordenamiento Territorial, en 1999.
Posteriormente se procedió a una verificación de los usos del suelo, comparando las
unidades presentes en el mapa en cuestión con las fotografías aéreas de 1998
pertenecientes a Planeación Municipal, en dos zonas escogidas al azar, una en el área de
Santa Elena y otra en el área de San Cristóbal, y no se encontró una correspondencia
adecuada.
Teniendo en cuenta que el mapa de usos del suelo es de suma importancia como insumo
en el modelo que se pretende implementar en este estudio, y que si éste modelo parte de
información errónea va a producir resultados erróneos, se decidió generar un mapa de
usos del suelo con la fuente de información más actualizada que se tiene y que posee la
escala más cercana a la escogida para este estudio: las fotografías aéreas que se
encuentran en Planeación Municipal, tomadas en 1998, con escalas aproximadas de
1:10000 en la zona urbana y de 1:20000 en la zona rural.
Universidad EAFIT 29
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Planchas urbanas:
Dicha cartografía se encuentra en escala 1:2000 con curvas de nivel cada dos metros,
cubriendo la totalidad de la zona urbana correspondiente al municipio de Medellín en las
cuencas del estudio. La lista de planchas es (75 planchas):
145, 151, 152, 153, 154, 155, 156, 160, 161, 162, 163, 164, 165, 169, 170, 171, 172, 173,
174, 178, 179, 180, 181, 182, 187, 188, 189, 190, 191, 192, 193, 194, 196, 197, 198, 199,
200, 201, 202, 203, 205, 206, 207, 208, 209, 210, 211, 212, 214, 215, 216, 217, 218, 219,
220, 221, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 237B, 238, 239,
240, 247, 248, 249.
Planchas rurales:
La cartografía rural se encuentra en escala 1:5000 con curvas de nivel cada cinco metros.
Se cubrió la totalidad de la zona rural de las cuencas del estudio. La lista de planchas
rurales es (20 en total):
Revisión de la información.
Cartografía base: la primera labor efectuada fue la revisión de cada uno de los archivos
correspondiente a cada categoría mencionada anteriormente de cada plancha (570
archivos en total), de la cartografía del proyecto SIGAME. Esta revisión fue necesaria
puesto que la cartografía entregada no había pasado por la interventoría final. Se
identificaron errores de edición y de relación de atributos, estos últimos solo se
manifestaban en los Sistemas de Información Geográfica utilizados en el proyecto
Universidad EAFIT 30
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Luego se unieron las planchas de cobertura de topografía e hidrografía tanto rural como
urbana por cuencas, se construyeron polígonos para cada cuenca que cubrieran una zona
algo mayor que la divisoria de aguas y se cortaron las planchas anteriormente pegadas
teniendo como límite los polígonos generados en la fase anterior, para cada cuenca.
Con respecto a la cartografía de Usos del Suelo, aunque se encontró que fue publicada
en el POT de Medellín como si estuviera actualizada a 1999, la correspondencia con las
fotografías aéreas de 1998 es muy pobre, y por lo tanto se determinó que su calidad no es
buena. Anteriormente se mencionaron las dificultades encontradas y las decisiones que
se tomaron al respecto (ver numeral 3.2, Revisión e interpretación de imágenes de
satélite y fotografías aéreas).
Universidad EAFIT 32
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
4 INSTRUMENTACIÓN
En el mes de mayo de 2000 se seleccionaron los sitios para la ubicación de las estaciones
tal como se especifica en la tabla 5.
Universidad EAFIT 34
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 35
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGIO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Como consecuencia de lo anterior, no fue posible efectuar ningún tipo de validación del
modelo desarrollado en el trabajo, además de la imposibilidad de implementar un sistema
de alarmas, puesto que no se conocen los umbrales de lluvia (intensidad y duración) que
produzcan inundaciones.
Universidad EAFIT 36
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
5 Análisis Cartográfico
En este capítulo se explica detalladamente el procesamiento de la información
cartográfica digital y los productos que se obtuvieron a partir del mismo. El análisis
cartográfico consta de 3 etapas diferentes: procesamiento de información base para
análisis morfométrico de las cuencas, superposición de mapas para la obtención de zonas
homogéneas y análisis espacial para la selección de las zonas más importantes en cada
cuenca para su posterior caracterización.
Con la cartografía base del proyecto SIGAME, topografía e hidrografía en escalas 1: 2000
para la zona urbana y 1: 5000 para la zona rural, se elaboraron modelos de elevación
digitales para las cuencas de las quebradas Santa Elena, La Castro, Minitas, La Quintana,
La Iguaná, La Hueso, Ana Díaz, La Picacha y Altavista.
Para la cuenca de la quebrada Doña María la cartografía del proyecto SIGAME no cubre
la parte baja por ser jurisdicción del municipio de Itagüí. Esta parte está cubierta por la
cartografía del proyecto SIGMA de Empresas Públicas de Medellín, pero no fue posible
tener acceso a ella por la negativa de EEPPM a la solicitud del SIMPAD de utilizarla. Por
lo tanto se decidió trabajar con cartografía base disponible en Planeación Municipal para
esta cuenca, con topografía e hidrografía en escala 1: 10000, menos detallada y menos
actualizada que la de los proyectos SIGAME y SIGMA.
Los MEDs son la base para el análisis hidrológico que se efectuó en la presente
investigación, pues a partir de éstos se generó la red de drenaje, la divisoria de aguas de
la cuenca y mediante un procesamiento de la información, las curvas hipsométricas de las
cuencas y los mapas de pendientes de las mismas y números de orden de las corrientes
según los métodos de Shreve y de Strahler. Además es posible extraer del MDE
secciones o perfiles de líneas rectas o irregulares, tales como las quebradas o las vías de
transporte o secciones transversales a las mismas. Otra aplicación importante de los
MEDs es la generación de matrices de dirección y acumulación de flujo que permiten
definir caudales instantáneos de todos los puntos de las cuencas a partir de una lluvia
espacializada y de un mapa de zonas homogéneas en cuanto a la relación precipitación –
escorrentía.
Universidad EAFIT 37
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
1. En primer lugar se exportaron los archivos recibidos del proyecto SIGAME (top.dgn)
que contenía las curvas de nivel en un formato que fuera posible leer en ArcInfo, esto
es el formato de transferencia dxf (top.dxf).
3. Por otro lado se construyeron polígonos por cuenca cuyos límites estuvieran por fuera
de la divisoria de aguas de cada cuenca, para esto se tuvo en cuenta la red de
drenaje de la cartografía suministrada por el proyecto SIGAME. Dicha cobertura se
nombró cdiv_cuenca.
Como se mencionó con anterioridad, los datos de entrada para la función topogrid
utilizados en la actual investigación, fueron la cobertura con curvas de nivel
(ctop_cuenca), la cobertura con la hidrografía (crds_cuenca) y la cobertura con el límite
aproximado de la cuenca (cdiv_cuenca).
Universidad EAFIT 38
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
asignarle la dirección de flujo al MDE, el sistema escoge las celdas con el valor de
elevación más bajo y el canal que genera en la red de drenaje es dibujado exactamente
por el medio de dichas celdas, por lo tanto, sería dibujado a 2.5 m del canal real.
5. Luego de este paso se obtuvieron MED para las cuencas los cuales presentaban
depresiones (sinks) debido a pequeños errores en la información básica (curvas de
nivel), a depresiones naturales del terreno y a errores propios del método de
interpolación. Es necesario llenar estos huecos o depresiones para poder aplicar las
funciones de flujo y así generar la red de drenaje y la divisoria de aguas de cada
cuenca, así como para aplicar el modelo hidrológico propuesto en la actual
investigación.
6. Para hallar las depresiones del MED el paso inicial es aplicarle la función flowdirection
que define la dirección de flujo para cada celda del MDE por comparación del valor de
la altura de dicha celda con sus ocho vecinas, asignando dicha dirección hacia la
celda que presenta una diferencia de altura mayor. Si una celda no posee una vecina
con un valor menor de altura, entonces se trata de una depresión y el flujo entonces
se concentrará en ella. Para los análisis hidrológicos es necesario que todas las
celdas fluyan hacia un mismo punto en una cuenca en particular: la desembocadura o
punto de cota más baja en toda la cuenca, es decir, el MED solo puede tener una
depresión que coincida con dicho punto.
Universidad EAFIT 39
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
A los modelos de elevación digital generados en este punto se les aplica funciones para
construir la red de drenaje y la divisoria de aguas. El siguiente diagrama de flujo
esquematiza el proceso de generación de un modelo de elevación digital sin depresiones
y a partir de éste construir la divisoria de aguas de la cuenca y la red de drenaje.
Universidad EAFIT 40
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 41
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Figura 12. Dirección de flujo: el valor representado es el valor que toma la celda del
centro de acuerdo con la diferencia de alturas de las celdas vecinas.
Por ejemplo si la dirección del flujo de la celda en cuestión es hacia la esquina superior
izquierda, entonces en la matriz resultante dicha celda tendrá un valor de 32.
Para hallar la distancia entre las celdas se multiplica la longitud de cada celda por 1 entre
celdas ortogonales y por 1.414 entre celdas diagonales.
Si todas las celdas vecinas presentan una elevación mayor que la celda analizada, esta
se considera una depresión (dirección de flujo indefinida) y debe ser llenada con el valor
más bajo de sus vecinas asignándole una dirección de flujo siguiendo esta celda. Si dos
celdas fluyen una hacia la otra, se consideran como una depresión y se les asigna una
dirección de flujo indefinida. Un ejemplo que representa la aplicación de la función
flowdirection en una matriz de elevación o MED es el siguiente:
flowdirection
gdem_cuenc gfd_cuenca
a
Figura 13. Aplicación de la función flowdirection en un modelo de elevación digital
hipotético.
Universidad EAFIT 42
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Por último se cambia la estructura de la matriz de raster a vector para construir coberturas
de polígonos que representan la cuenca o área de captación, esto se logra con la
aplicación del comando gridpoly del sistema ArcInfo. La cobertura generada se denominó
cwshed_cuenca.
Para las cuencas del estudio, se consideró la cuenca de drenaje hasta un punto de interés
que coincide con el punto seleccionado para la instrumentación limnigráfica o limnimétrica
en cada una de ellas, el cual no siempre coincide con la desembocadura de la corriente
principal. De las cuencas del estudio, sólo coincide el punto de interés con la
desembocadura de la cuenca en las quebradas La Picacha, La Hueso, La Iguaná y La
Castro.
Para el caso de la cuenca de la quebrada Altavista, este punto está a la altura de la calle
65, en el cruce con la calle 30. Aguas debajo de dicho punto la quebrada recibe las aguas
de la quebrada La Guayabala, por lo tanto, si se considera su desembocadura en el río
Medellín, se estaría teniendo en cuenta las redes de drenaje de las dos cuencas en lugar
de solamente la primera.
Para la cuenca de la Ana Díaz, se escogió el punto donde ésta le cae a la quebrada La
Hueso para así poder considerar el área de drenaje de las dos cuencas por separado.
En la cuenca de la quebrada Santa Elena, se consideró el punto más bajo como aquél en
el que la quebrada entra al box coulvert de la avenida La Playa, ya que a partir de ese
punto y hasta cerca de su desembocadura en el río Medellín el cauce está completamente
cubierto.
Para generar una red de drenajes a partir de un modelo de elevación digital se deben
efectuar los siguientes pasos: aplicar la función flowdirection al modelo de elevación digital
(ver numeral anterior); a la matriz generada se le aplica la función flowaccumulation y a la
matriz arrojada por este proceso se le define el umbral de flujo, el cual corresponde al
número de celdas que drenan a un punto en particular a partir del cual se presenta un
Universidad EAFIT 43
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
El comando flowaccumulation genera una matriz de flujo acumulado para cada celda,
calculando el total de celdas que fluyen hacia cada celda de la matriz. Los parámetros de
entrada son la matriz de dirección de flujo y una matriz de peso que representa el peso
asignado a cada celda de la matriz. Si no se especifica la matriz de peso, el sistema le
aplica un peso de uno a cada celda y así la matriz resultante será el número total de
celdas que fluyen hacia cada celda.
Las celdas con una gran acumulación de flujo son celdas con flujo concentrado y pueden
ser utilizadas para definir las corrientes permanentes, mientras que las celdas con valor
de acumulación igual a cero son máximos topográficos que representan cuchillas o topes
de colinas.
Si se asigna una matriz de peso que represente la cantidad de escorrentía generada por
una celda en una lluvia en particular, entonces la matriz resultante será el caudal
instantáneo que presenta cada celda en el evento analizado.
flowaccumulation
gfd_cuenca gfa_cuenca
Definición del valor umbral a partir del cual se presenta un drenaje permanente.
Universidad EAFIT 44
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Un recurso usual es generar varias redes de drenaje con diferentes valores umbrales
teniendo en cuenta condiciones similares para la totalidad de la cuenca. Posteriormente
se efectúa una verificación visual con una red de drenaje generada a partir de una
fotointerpretación en la escala de trabajo y se define cual de los valores umbrales se
acerca más a la red de drenaje real.
En la actual investigación se definió el valor umbral en 500 celdas por comparación con la
cobertura de hidrografía de la cartografía suministrada por planeación municipal, que se
encuentra en escala 1:10000, igual a la escala de trabajo del actual proyecto. Este valor
corresponde a un área de captación de 12500 m2, que sería el área necesaria para que se
presente una corriente cartografiable a dicha escala.
El comando utilizado para el proceso descrito se ejecuta en el módulo Grid del sistema
ArcInfo y su sintaxis es:
Lo que significa que a las celdas de la matriz de acumulación que presentan un valor
mayor de 500 se les asigna el valor de 1, mientras que a las demás se les asigna el rótulo
de NODATA.
Teniendo la red drenaje generada con la aplicación del proceso descrito en el numeral
anterior se aplica la función streamorder del sistema ArcInfo para encontrar los números
de orden de las corrientes según los métodos propuestos por Strahler (1957) y Shreve
(1966).
El orden de las corrientes identifica y clasifica los tipos de corriente basados en el número
de corrientes tributarias. En los dos métodos propuestos se asigna un número de orden
de 1 a las corrientes que no poseen tributarios. En el método de Strahler el orden de la
corriente se incrementa cuando dos corrientes del mismo orden se interceptan, así la
intersección de dos corrientes de orden 1 genera una corriente de orden dos, mientras la
intersección de dos corrientes de orden 2 genera una corriente de orden 3. La
Universidad EAFIT 45
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
El método de Shreve cuenta todas las uniones entre la red de drenaje. Al igual que el
método de Strahler se le asigna un número de orden 1 a las corrientes sin tributarios, la
diferencia es que todas las intersecciones cuentan para aumentar el número de orden de
las corrientes, así la unión de una corriente de orden 2 con una corriente de orden 3
genera una corriente de orden 5, es decir, es un método aditivo. Debido a esto, el número
de orden de Shreve se toma más como una descripción de la magnitud más que el orden
o jerarquía de la corriente. La magnitud de una corriente según este método es el número
de intersecciones que convergen en esta.
Universidad EAFIT 46
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Las siguientes figuras presentan la red de drenaje con el número de orden según Strahler
y Shreve generadas a partir de los modelos de elevación digital para la cuenca de la
quebrada La Castro. Los números de orden del resto de cuencas del estudio se presentan
en sus respectivos anexos.
Universidad EAFIT 47
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 48
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
% del área acumulada, pues toda el área de la cuenca tiene elevaciones que están por
encima de este último.
Para la construcción de las curvas, se tomó el modelo digital de elevación de cada cuenca
en el sistema ArcInfo y se clasificó en intervalos de elevación de 50 m a los cuales se les
asignó el valor de su altura media, por ejemplo, al intervalo que está comprendido entre
1650 y 1700 msnm, se le asignó el valor de su altura media, 1675 msnm. A este nuevo
modelo de elevación clasificado se le construyó una tabla de valores de atributos que
indica cuántas celdas pertenecen a cada intervalo. La tabla VAT se exportó a un archivo
ascii y éste fue importado desde el programa Microsoft Excel, en el cual se procedió a
ampliar la tabla calculando el área de cada intervalo (multiplicando el número de celdas
por 25 m2 que es el área que representa cada celda), el porcentaje de área de cada
intervalo con respecto al total de la cuenca y el porcentaje de área acumulada de cada
intervalo desde el intervalo más alto al más bajo. Por último se graficó la columna de
elevaciones, en las ordenadas, contra la columna de porcentaje de área acumulada, en
las abscisas.
Universidad EAFIT 49
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
La curva hipsométrica de una cuenca puede considerarse como una especie de perfil
generalizado de la cuenca de drenaje (Remeneiras, 1974). Según Monsalve (1995), una
curva hipsométrica cóncava representa una cuenca con valles amplios y filos escarpados
en la parte alta de la cuenca, mientras que una convexa representa una cuenca con valles
profundos y zonas planas hacia la parte alta. En la realidad la interpretación no es tan
fácil, ya que lo que se tiene es una situación más compleja en el relieve de cada cuenca.
En el caso de la cuenca Ana Díaz, las pendientes fuertes predominan también en la parte
alta de la cuenca, con un cambio que puede ubicarse alrededor de los 2450 msnm,
pasando a pendientes moderadas que predominan hasta cerca de los 1700 msnm, punto
en el que se presenta un nuevo cambio de pendiente y que puede asimilarse al inicio de
la zona de pendientes bajas en la cuenca, que es mucho más pequeña, en porcentaje de
área, que las de las otras cuencas vecinas y que coincidiría con la zona de la planicie
aluvial del río Medellín. Esto se debe en parte a que aquí se considera la cuenca de la
quebrada Ana Díaz hasta el punto de su desembocadura en la quebrada La Hueso, a la
altura del colegio Jesús María.
Universidad EAFIT 50
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
pendientes bajas a moderadas y debajo de los 1900 msnm las laderas de la cuenca
presentan pendientes mucho más fuertes, lo cual es evidencia de la incisión de la
quebrada en este tramo y el encañonamiento que se observa en la zona cercana a su
desembocadura, en los barrios Caicedo y Villa Liliam.
donde:
y = diferencia de alturas entre celdas
x = distancia entre celdas (tamaño de celda x 1 para celdas ortogonales ó
tamaño de celda x 1.414 para celdas diagonales)
θ = pendiente en ángulo
θ = tan-1 (y/x)
El proceso de cálculo utilizado por el sistema ArcInfo tiene en cuenta una matriz vecina de
3x3 para cada celda y por medio de una técnica denominada el máximo promedio, calcula
la pendiente promedia máxima para cada celda.
Universidad EAFIT 51
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 52
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Cuencas
Zona Suroccidental La Iguaná y Doña María Santa Elena Zona Noroccidental
0° a 10° 0° a 10° 0° a 13° 0° a 10°
11° a 22° 11° a 20° 14° a 22° 11° a 20°
22° a 32° 21° a 31° 23° a 32°
21° a 52°
33° a 60° 32° a 61° 32° a 61°
Para esta clasificación se agruparon las cuencas de la zona suroccidental, puesto que
presentan características litológicas y de uso del suelo similares, y se definió la misma
clasificación para las cuatro cuencas. En el caso de la cuencas Santa Elena y La Castro,
se determinaron los intervalos para toda la cuenca Santa Elena y se trabajó la subcuenca
de la quebrada La Castro con los mismos intervalos. En el caso de La Iguaná y la Doña
María, los intervalos resultantes fueron muy similares, por lo que se decidió unificarlos
para estas cuencas.
Con los intervalos de pendientes por cuencas definidos en el numeral anterior se genera
una nueva matriz que presenta las clases espacializadas. Esto se obtiene reclasificando
la matriz de pendientes con los intervalos definidos, los cuales se especifican en una tabla
denominada REMAP TABLE, que es un archivo de texto que el sistema ArcInfo toma
como parámetro de entrada. Estas nuevas matrices son las que se utilizan para la
generación de zonas con coeficiente de escorrentía homogéneo, con la suposición de que
en los intervalos definidos la relación precipitación – escorrentía presenta un
comportamiento similar. A continuación se presenta el mapa con los intervalos de
pendientes obtenidos con la aplicación del proceso descrito para la cuenca de la
quebrada La Castro, para el resto de cuencas del estudio este mapa se presenta en sus
respectivos anexos.
Universidad EAFIT 53
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 54
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Parámetros
Cuencas:
morfométricos:
La Santa La Doña
La Castro La Hueso Ana Díaz Altavista Minitas Quintana
Iguaná Elena Picacha María
Area de la cuenca
(km2):
51.22 39.46 4.22 12.38 9.17 10.65 13.71 1.28 2.44 72.6
Longitud del cauce
principal (km):
15.29 12.08 4.238 6.767 9.275 10.878 10.549 4.187 4.787 22.1
Longitud del cauce
principal hasta la 16.37 12.27 4.423 7.043 9.479 11.028 10.725 4.502 5.035 22.2
divisoria (km):
Cota mayor del cauce
(msnm):
2555 2620 2500 2132 2485 2424 2227 2146 2206 2906
Cota mayor de la
cuenca (msnm):
3137 2800 2560 2346 2661 2540 2411 2223 2328 3081
Cota menor (msnm): 1457 1505 1630 1478 1480 1464 1475 1446 1444 1522
Perímetro de la cuenca
(km):
50.85 40.6 14.48 22.03 26.56 27.72 28.67 11.96 13.81 60.3
Diferencia de cotas en
el cauce principal (m):
1098 1115 870 654 1005 960 752 700 762 1384
Pendiente del cauce
(%) (dif. entre ptos 7.18 9.23 20.53 9.66 10.84 8.83 7.13 16.72 15.92 6.26
extremos):
Pendiente media de la
cuenca (%):
39.03 33.17 44.41 33.53 44.26 33.40 35.10 25.25 23.82 40.40
Pendiente media de la
cuenca (grados):
21.32 18.35 23.94 18.53 23.87 18.47 19.40 14.17 13.40 22
Índice de Gravelius.
Kc = 0.28P / A1/2
Donde P es el perímetro de la cuenca, en Km, y A, el área de drenaje en Km2. Mientras
más irregular sea la cuenca de drenaje, es decir, de forma más diferente a un círculo,
mayor será su índice de Gravelius. Las cuencas de forma aproximadamente circular
tienen coeficientes de Gravelius mínimos, con valores cercanos a 1, y tienen mayor
tendencia a presentar inundaciones, ya que el tiempo de concentración tiende a ser
menor en cuencas de forma circular que en cuencas de forma elongada o de forma
rectangular. En la tabla 9 se presenta el Índice de Gravelius para las cuencas del estudio.
Universidad EAFIT 55
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Factor de forma.
El factor de forma relaciona el ancho medio de la cuenca con la longitud axial de la misma
(Monsalve, 1995). Ésta última se mide siguiendo la longitud del cauce principal, corriente
más larga, desde la desembocadura hasta la divisoria de aguas de la cuenca. El ancho
medio de la cuenca se obtiene de dividir el área de la cuenca por su longitud axial.
Kf = A / L2
En donde A es el área de la cuenca en Km2 y L la longitud axial de la cuenca en Km. En
una cuenca con un factor de forma bajo es probable que se presenten menos
inundaciones que en una cuenca del mismo tamaño pero con un factor de forma más alto.
La importancia de calcular estos factores para las cuencas del estudio, en particular al
comparar las de área similar, radica en que éstos pueden dar una idea preliminar de
cuáles de ellas merecen más atención en cuanto a las labores de prevención de
inundaciones, ya que es probable que su ocurrencia sea mayor que en las otras. En la
tabla 9 se presenta el factor de forma de las cuencas del estudio.
Tabla 9. Índices de forma de las cuencas del estudio trabajadas hasta el momento.
Índices de Cuencas
Santa La Doña
forma La Iguaná
Elena
La Castro La Hueso Ana Díaz
Picacha
Altavista Mnitas Quintana
María
Indice de
Gravelius 1.99 1.81 1.97 1.75 2.46 2.38 2.17 2.96 2.48 1.98
Factor de
forma 0.18 0.24 0.18 0.22 0.09 0.08 0.11 0.05 0.08 0.14
Densidad de drenaje.
Universidad EAFIT 56
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Tabla 10. Densidad de drenaje en km/km2 para las cuencas del estudio, con base en la
información obtenida de las redes de drenaje generadas a partir de los MDEs.
Cuencas
Doña
Iguaná St Elena Castro Hueso Ana Díaz Picacha Altavista Minitas Quintana
María
6.76 6.48 5.57 5.16 5.64 6.45 6.22 8.77 8.26 5.23
Teniendo en cuenta estas ideas, cabe comparar la densidad de drenaje de las cuencas de
la zona suroccidental de la ciudad, en las cuales la litología y el relieve son similares. En
la tabla se observa una densidad de drenaje un poco menor en las cuencas de las
quebradas Ana Díaz y La Hueso que en las cuencas de las quebradas Altavista y La
Picacha, lo que podría estar indicando leves diferencias en la pluviosidad de ambas
zonas, con presencia de clima más seco al norte que al sur. Esta hipótesis se puede
comprobar luego de obtener datos a partir de la instrumentación efectuada en la presente
investigación. Una observación de la vegetación presente en las cuencas de las
quebradas de la zona suroccidental que están ubicadas hacia el norte, es decir, La Hueso
y la Ana Díaz, permite inferir que la pluviosidad puede variar tal como se expuso
anteriormente.
Mapas de Pendientes.
Los usos del suelo tienen una gran importancia en el comportamiento hidrológico de una
zona, puesto que la cobertura vegetal determina en gran medida la relación entre la
precipitación y la escorrentía generada en un evento de lluvia. En la actual investigación
se hizo una delimitación cartográfica de usos del suelo teniendo en cuenta las siguientes
categorías:
Universidad EAFIT 57
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Suelo desnudo: Zonas degradadas por minería, zonas afectadas por movimientos en
masa que no presenten cobertura vegetal. Son zonas con altos coeficientes de
escorrentía y de gran aporte de sedimentos a las corrientes.
El mapa temático de usos del suelo de la zona de estudio no se pudo tomar de ninguna
de las fuentes existentes, lo que hizo necesaria la fotointerpretación de las fotografías
aéreas de 1998 existentes en la Planoteca del Municipio de Medellín, que fueron la base
para la fotorestitución cartográfica digital utilizada en los proyectos SIGMA y SIGAME.
Luego de acceder al mapa suministrado por Planeación Municipal y encontrar que la
información cartográfica no se dejaba procesar, se averiguó la procedencia de dicho mapa
y se encontró que éste se había levantado en el año de 1987 con una verificación en
campo tendiente a la ubicación de zonas de protección para el Plan de Ordenamiento
Territorial, en 1999.
Posteriormente se procedió a una verificación de los usos del suelo, comparando las
unidades presentes en el mapa en cuestión con las fotografías aéreas de 1998, en dos
zonas escogidas al azar, una en el área de Santa Elena y otra en el área de San
Cristóbal, y no se encontró una correspondencia adecuada. Debido a lo anterior se
decidió generar un mapa de usos del suelo con la fuente de información más actualizada
que se tiene y que posee la escala más cercana a la escogida para este estudio: las
fotografías aéreas que se encuentran en Planeación Municipal, tomadas en 1998, con
escalas aproximadas de 1:10000 en la zona urbana y de 1:20000 en la zona rural.
Universidad EAFIT 58
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
útiles para ubicar bien las zonas y encontrar la correspondencia precisa entre los mapas y
las fotografías aéreas, de manera que la información fuera lo más precisa posible y
pudiera producirse rápidamente. Una vez estaban delimitadas las zonas en los nuevos
mapas, se procedió a la digitalización de algunas líneas que cerraban los polígonos y a
editar las coberturas. La edición se inició con la generación de una cobertura con
topología de polígonos por cuenca a partir de las polilíneas organizadas en el paso
anterior y finalmente con la inclusión de un campo en la tabla de atributos en donde se
especificaba el uso del suelo de acuerdo con un código establecido para cada uso el cual
se reseñará posteriormente. En la figura 23 se presenta el mapa de usos del suelo para la
cuenca de la quebrada La Castro, para el resto de cuencas los mapas de usos del suelo
se presentan en los respectivos anexos.
Universidad EAFIT 59
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 60
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Una vez obtenidos los mapas de Formaciones Superficiales, Usos del Suelo y Pendientes
se realizó una codificación de cada una de las categorías de los mapas. Para este
proceso se añadió un campo numérico a la tabla de atributos de cada uno de los mapas
para cada cuenca, y se asignaron los códigos de acuerdo a la categoría de la variable
representada, tal como se ilustra en las tablas 11, 12 y 13.
Universidad EAFIT 61
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 62
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Una vez obtenidos los mapas en formato raster, se procedió a realizar una superposición
simple, mediante la suma de mapas, para identificar las zonas que tienen las mismas
características en cuanto a tipo de formación superficial, tipo de uso del suelo e intervalo
de pendiente. Sumando los códigos de cada uno de los mapas se obtiene el código de la
zona homogénea en cada celda. Por ejemplo, para una celda perteneciente a la cuenca
de la quebrada Santa Elena puede tenerse:
Universidad EAFIT 63
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
+ +
Formaciones
Uso del Suelo Pendientes Superficiales
Zonas
Homogéneas
Figura 26. Esquema de superposición de los mapas en formato raster, mediante una
suma sencilla de los valores correspondientes a cada celda para la obtención del mapa de
zonas homogéneas con respecto al coeficiente de escorrentía.
Universidad EAFIT 64
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Luego de hacer la superposición para encontrar las zonas en cada cuenca, se hizo un
análisis de la importancia relativa según el área, o porcentaje de cubrimiento, de cada una
de esas zonas con respecto a la cuenca. Con este análisis se determinaron las zonas más
importantes según su área dentro de cada cuenca, hacia las cuales se orientaron los
ensayos con el simulador de lluvias para determinar el coeficiente de escorrentía.
Universidad EAFIT 65
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Para definir las zonas en las que se efectuaron los ensayos, se escogió un valor umbral
del 5% del área rural, es decir, se decidió efectuar ensayos en aquellas zonas que
tuvieran un área acumulada sobre la parte rural de la cuenca igual o mayor al 5 % (p. ej:
ver tabla 14, valores escritos en rojo oscuro). De dichas zonas se descartaron las zonas
de bosque, ya que por la naturaleza de los ensayos, los datos obtenidos en ensayos
realizados en bosque están sobredimensionados, puesto que no se tiene en cuenta el
porcentaje de captación y almacenamiento de la lluvia por parte de la vegetación del
bosque.
Con base en las tablas como la mostrada anteriormente se construyeron los histogramas
de frecuencias que se presentan en las siguientes figuras y se seleccionaron las zonas en
las cuales se efectuaron los ensayos. La línea roja señalada en los histogramas de
frecuencia indica el límite del 5 %.
Universidad EAFIT 66
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 67
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 68
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 69
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Tabla 15. Zonas seleccionadas para ensayos con simulador de lluvias en cada cuenca,
área cubierta por bosques y su porcentaje de área sobre el área rural de la misma.
Universidad EAFIT 70
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Lo anterior se hizo con el fin de tener un valor estimado para comparar entre pastos y uso
silvopastoril para la misma pendiente y misma formación superficial, para tener un
indicativo de cuál es el orden del coeficiente de escorrentía en los bosques medido con el
simulador de lluvias y para tener un valor del coeficiente de escorrentía en zonas con
suelo desnudo y comparar con zonas de pastos en el mismo intervalo de pendiente y para
la misma formación superficial.
Universidad EAFIT 71
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
El objeto final de los pasos seguidos es determinar la relación entre la intensidad de una
lluvia y el coeficiente de escorrentía para cada celda perteneciente a una zona
homogénea con relación a su comportamiento hidrológico, caracterizado por el coeficiente
de escorrentía. Esto es, la obtención de mapas que presenten el coeficiente de
escorrentía espacializado para cada cuenca de estudio.
Las zonas homogéneas surgen a partir de la superposición de mapas de usos del suelo,
formaciones superficiales y pendientes, teniendo en cuenta que es la vegetación
(representada en el mapa de usos del suelo), los primeros niveles del suelo
(representados por el mapa de formaciones superficiales) y la pendiente los parámetros
que más influyen en el comportamiento hidrológico de una zona.
El resultado final de los ensayos para determinar los coeficientes de escorrentía en zonas
homogéneas es el mapa de zonas productoras de escorrentía por cuenca, insumo
fundamental en la toma de decisiones para planificar sistemas de alarma y acciones
tendientes a la regulación del coeficiente de escorrentía para reducir el riesgo por
inundación cerca de los cauces de las quebradas (cambios en los usos del suelo).
Universidad EAFIT 72
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Figura 37. Simulador de lluvias portátil. Diseño inicial publicado por McQueen en 1963
(tomado de McQueen, 1963, U.S. Geological Survey Circular 482).
Universidad EAFIT 73
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 74
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Tubo de Mariotte
Probeta o cilindro
graduado
Manguera de
alimentación
El tamaño de las gotas está controlado por el tamaño del filamento introducido en cada
aguja, a mayor diámetro de filamento, menor tamaño de gota y por ende, menor
intensidad de la lluvia aplicada.
Universidad EAFIT 75
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Círculo de acrílico
Agujas con
filamentos de
nylon (amarillos
y rojos)
Manguera de
alimentación
Tornillo o
pasador
Universidad EAFIT 76
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Pantalla contra
el viento
Trípode de
topografía
Figura 41. Trípode de soporte y pantalla contra el viento del simulador de lluvias portátil.
En el presente estudio se trabajó con una unidad de base cilíndrica de acero inoxidable
(ver figuras 42 y 43), cuya función es únicamente asegurar que el flujo de agua que se da
en la parcela sea vertical, es decir, que no hay pérdidas de agua por flujo lateral hacia el
exterior de la parcela, y así permitir medir el agua en exceso que se succiona de la
parcela con el sistema para medir la escorrentía.
Universidad EAFIT 77
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Figura 42. Detalle de la unidad de base cilíndrica del simulador de lluvias portátil.
Esta unidad se entierra en el suelo y delimita la parcela del ensayo, a una profundidad
que varía entre los 5 y 15 centímetros dependiendo del tipo de suelo, hasta que se pueda
asegurar que no hay pérdidas de agua por flujo lateral a través de los límites de la
parcela. De esta forma se asegura que el agua de la lluvia que no se infiltra en el suelo se
acumule adentro de la parcela y que pueda ser succionada y medida en el sistema para
medir la escorrentía.
Universidad EAFIT 78
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Escudo para
salpicadura
(embudo)
Unidad de
base cilíndrica
Figura 43. Acople de la unidad base al resto del sistema, con el embudo para evitar
pérdidas por desviación de gotas.
Universidad EAFIT 79
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Mangueras
Probeta de
almacenamiento
Pera de
succión
Universidad EAFIT 80
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Es necesario asegurarse de que todos los filamentos del simulador están generando
gotas antes de comenzar el ensayo, para lo cual se hace un chequeo rápido y se induce
la generación de gotas en todos los filamentos removiendo las burbujas de aire en los
conductos (agujas de jeringa) e introduciendo agua a presión con una pequeña
manguera.
Existe una relación directa no lineal entre la altura de la columna de agua dentro del
simulador y la intensidad de la lluvia simulada, pero dicha relación se ve afectada, aunque
en forma mínima, por los cambios que induce el viento en la presión dentro del tubo de
Mariotte y por los cambios en la presión atmosférica durante el día.
Toma de datos:
Por último, se conecta la manguera de alimentación a la cabeza del simulador, se espera
a que se igualen las presiones de la unidad de alimentación y control y de la cabeza del
simulador, y una vez éstas igualen comienza la toma de datos.
Para la toma de datos se espera que el nivel del agua llegue a 1000 ml dentro de la
unidad de almacenamiento para tomar las condiciones iniciales del ensayo. Una vez se
alcanza este nivel, se lleva el cronómetro a ceros, se anota la hora del ensayo y la presión
atmosférica, y se anota el nivel del agua en la probeta de almacenamiento de escorrentía.
Cada 5 minutos se lee el nivel del agua dentro de la unidad de alimentación y control, el
nivel del agua dentro de la probeta de almacenamiento de escorrentía y se chequea el
nivel del agua dentro de la cabeza del simulador para asegurarse de que el agua que está
llegando a la parcela en forma de lluvia es la misma cantidad de agua que está pasando
de la unidad de almacenamiento y control a la cabeza del simulador. Si esto no sucede,
se debe corregir el volumen de lluvia con la diferencia de niveles dentro de la cabeza del
simulador. En la tabla 16 se presenta un ejemplo de la toma de datos.
Universidad EAFIT 81
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Lecturas de campo
0 1000 120
5 900 210
10 825 278
15 765 337
20 715 385
25 675 418
30 620 460
35 575 495
40 540 525
45 500 560
50 460 590
55 420 620
60 380 645
0 1000 120
5 900 210 60 60 54 54 0.9 0.9
10 825 278 45 52.5 40.8 47.4 0.90666667 0.90285714
15 765 337 36 47 35.4 43.4 0.98333333 0.92340425
20 715 385 30 42.75 28.8 39.75 0.96 0.92982456
25 675 418 24 39 19.8 35.76 0.825 0.91692307
30 620 460 33 38 25.2 34 0.76363636 0.89473684
35 575 495 27 36.4285714 21 32.1428571 0.77777778 0.88235294
40 540 525 21 34.5 18 30.375 0.85714286 0.88043478
Universidad EAFIT 82
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
La intensidad de la lluvia simulada, que es la cantidad de lluvia por unidad de área por
unidad de tiempo, se calcula de la siguiente forma:
Las unidades de volumen son mililitros, las de área centímetros cuadrados y las de tiempo
minutos, para convertir ml / cm2 min a mm / hr, se tiene en cuenta lo siguiente:
teniendo en cuenta que esta lluvia se aplica sobre una parcela de 200 cm2, la intensidad
de la lluvia unitaria tenida en cuenta anteriormente, es
I = 0.3 cm /hr,
Ya que 1 cm = 10 mm,
I = 3 mm / hr,
I = 3 v / t (mm / hr)
Son varias las fuentes de error observadas durante el proceso de simulación de lluvias.
Estos errores que se reflejan en datos inconsistentes en la relación entre la intensidad de
la lluvia aplicada a la parcela y el coeficiente de escorrentía generado, se pueden
clasificar en dos tipos: los debidos al funcionamiento del equipo y los debidos a errores de
lectura.
Universidad EAFIT 83
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
La manguera que alimenta la probeta receptora del agua en exceso (aquella que no se
infiltra) dejaba de recoger agua. Esto se debía a que la entrada de la misma se obstruía
con partículas de suelo o se desacomodaba, quedando en algún punto de la parcela
diferente al punto más bajo (que es en donde se concentra el agua que no se infiltra).
La escogencia de las zonas homogéneas a las cuales se les efectuó ensayos con el
simulador de lluvia surge a partir del análisis espacial, detallado en el capítulo 5 (análisis
cartográfico).
Universidad EAFIT 84
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
La ubicación de los sitios en los que se llevaron a cabo los ensayos se muestra en la
figura 46.
Figura 46. Localización de los sitios en los que se efectuaron ensayos con el simulador de
lluvias.
Universidad EAFIT 85
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Al iniciar cada ensayo el suelo se encuentra saturado: Si se considera que las lluvias
máximas que generan caudales máximos se presentan durante un período invernal,
cuando el suelo se encuentra saturado, los ensayos de simulación de lluvias aplicados en
la actual investigación, tendiente a encontrar la relación entre la intensidad de la lluvia y el
coeficiente de escorrentía, se acercan a las condiciones reales si se prehumedece el
suelo hasta saturarlo.
Las tablas y figuras muestran además que la tendencia de cada ensayo si presenta una
correspondencia entre la intensidad de la lluvia y el coeficiente de escorrentía, por lo cual
es posible encontrar una relación entre estos dos parámetros.
Universidad EAFIT 86
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
0 1000 75
5 960 82 24 24 4.2 4.2 0.175 0.175
10 920 90 24 24 4.8 4.5 0.2 0.1875
15 885 95 21 23 3 4 0.142857143 0.173913043
20 850 103 21 22.5 4.8 4.2 0.228571429 0.186666667
25 815 108 21 22.2 3 3.96 0.142857143 0.178378378
30 780 115 21 22 4.2 4 0.2 0.181818182
35 745 120 21 21.85714286 3 3.857142857 0.142857143 0.176470588
40 710 125 21 21.75 3 3.75 0.142857143 0.172413793
45 670 135 24 22 6 4 0.25 0.181818182
50 630 142 24 22.2 4.2 4.02 0.175 0.181081081
55 580 155 30 22.90909091 7.8 4.363636364 0.26 0.19047619
60 538 162 25.2 23.1 4.2 4.35 0.166666667 0.188311688
65 490 170 28.8 23.53846154 4.8 4.384615385 0.166666667 0.18627451
70 447 180 25.8 23.7 6 4.5 0.23255814 0.189873418
75 400 188 28.2 24 4.8 4.52 0.170212766 0.188333333
Figura 47. Cambio del coeficiente de escorrentía con respecto al tiempo durante el
ensayo 12.
Universidad EAFIT 87
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
0 1000 95
5 920 120 48 48 15 15 0.3125 0.3125
10 830 145 54 51 15 15 0.277777778 0.294117647
15 740 165 54 52 12 14 0.222222222 0.269230769
20 650 185 54 52.5 12 13.5 0.222222222 0.257142857
25 560 212 54 52.8 16.2 14.04 0.3 0.265909091
30 470 235 54 53 13.8 14 0.255555556 0.264150943
35 410 257 36 50.57142857 13.2 13.88571429 0.366666667 0.274576271
40 330 280 48 50.25 13.8 13.875 0.2875 0.276119403
42.3 290 290 48 50.11764706 12 13.76470588 0.25 0.274647887
Figura 48. Cambio del coeficiente de escorrentía con respecto al tiempo durante el
ensayo 31.
Universidad EAFIT 88
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
0 1000 90
5 940 127 36 36 22.2 22.2 0.616666667 0.616666667
10 880 163 36 36 21.6 21.9 0.6 0.608333333
15 840 195 24 32 19.2 21 0.8 0.65625
20 800 225 24 30 18 20.25 0.75 0.675
25 750 257 30 30 19.2 20.04 0.64 0.668
30 700 285 30 30 16.8 19.5 0.56 0.65
35 660 310 24 29.14285714 15 18.85714286 0.625 0.647058824
40 615 340 27 28.875 18 18.75 0.666666667 0.649350649
45 575 370 24 28.33333333 18 18.66666667 0.75 0.658823529
50 540 395 21 27.6 15 18.3 0.714285714 0.663043478
Figura 49. Cambio del coeficiente de escorrentía con respecto al tiempo durante el
ensayo 80.
Universidad EAFIT 89
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
0 1000 68
5 920 85 48 48 10.2 10.2 0.2125 0.2125
10 840 100 48 48 9 9.6 0.1875 0.2
15 760 115 48 48 9 9.4 0.1875 0.195833333
20 690 125 42 46.5 6 8.55 0.142857143 0.183870968
25 615 140 45 46.2 9 8.64 0.2 0.187012987
30 540 150 45 46 6 8.2 0.133333333 0.17826087
35 470 155 42 45.42857143 3 7.457142857 0.071428571 0.164150943
40 400 170 42 45 9 7.65 0.214285714 0.17
45 330 180 42 44.66666667 6 7.466666667 0.142857143 0.167164179
50 265 195 39 44.1 9 7.62 0.230769231 0.172789116
55 195 205 42 43.90909091 6 7.472727273 0.142857143 0.170186335
60 130 215 39 43.5 6 7.35 0.153846154 0.168965517
Figura 50. Cambio del coeficiente de escorrentía con respecto al tiempo durante el
ensayo 135.
Universidad EAFIT 90
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Universidad EAFIT 91
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Cuando el coeficiente de correlación, R2, es muy bajo, como en los sitios 9, 10, 12, 19, 23
y 25, indica que no existe relación entre los datos de Coeficiente de Escorrentía e
Intensidad de la lluvia simulada. En estos casos la ecuación de regresión tiende a ser una
línea recta, o casi recta, paralela al eje x, lo que indica que el coeficiente de escorrentía
tiende en esos sitios a un valor constante independiente de la intensidad de la lluvia
simulada. A continuación se presentan las gráficas que relacionan las intensidades de
lluvia y el coeficiente de escorrentía por sitio.
Sitio 01:
Fecha de los ensayos: Febrero 9 de 2001.
Localización: Divisoria de aguas entre las cuencas de la Picacha y Doña María (Carretera
Aguas Frías- San Antonio de Prado.
Anotaciones del lugar: Pastos sobre saprolito del Stock de Altavista, pendiente 0o.
Código de la zona homogénea: 201008
Universidad EAFIT 92
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Sitio 02:
Fecha de los ensayos: Febrero 12 de 2001.
Localización: Divisoria de aguas entre las cuencas de la Picacha y Altavista (Carretera
Aguas Frías – San Antonio de Prado.
Anotaciones del lugar: Pastos sobre saprolito del Stock de Altavista, pendiente 0o.
Código de la zona homogénea: 201008.
Sitio 03.
Fecha de los ensayos: Febrero 14 de 2001.
Localización: Divisoria de aguas entre las cuencas de la Picacha y Altavista (Carretera
Aguas Frías - San Antonio de Prado).
Anotaciones del lugar: Pastos sobre saprolito del Stock de Altavista, pendiente de 22o a
26o
Código de la zona homogénea: 203008
Universidad EAFIT 93
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Sitio 04:
Fecha de los ensayos: Febrero 16 de 2001.
Localización: Divisoria de aguas entre las cuencas de la Picacha y Altavista (Carretera
Aguas Frías - San Antonio de Prado)
Anotaciones del lugar: Pastos sobre saprolito del Stock de Altavista, pendiente de 15o
Código de la zona homogénea: 202008
Sitio 07:
Fecha de los ensayos: Febrero 21 de 2001.
Localización: Barrio Pajarito, cerca de la escuela Alfonso Upegui Orozco.
Anotaciones del lugar: Pastos sobre suelo residual de depósito de flujo de lodos y
escombros, pendiente de 16o.
Código de la zona homogénea: 202014
Universidad EAFIT 94
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Sitio 08.
Fecha de los ensayos: Febrero 23 de 2001
Localización: Margen norte de la carretera a San Pedro, un km arriba de las partidas (500
metros antes del peaje)
Anotaciones del lugar: Pastos sobre suelo residual de Anfibolita, pendiente de 15o.
Código de la zona homogénea: 202006
Sitio 09.
Fecha de los ensayos: Febrero 27 de 2001.
Localización: Vía a San Pedro, frente a la escuela Santa Angela.
Anotaciones del lugar: Pastos sobre depósito coluvial, pendiente de 8o.
Código de la zona homogénea: 201015
Universidad EAFIT 95
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Sitio 10.
Fecha de los ensayos: Febrero 28 de 2001.
Localización: Vía a San Pedro, frente a la escuela Santa Angela.
Anotaciones del lugar: Pastos sobre depósito coluvial, pendiente de 16o.
Código de la zona homogénea: 202015
Sitio 11.
Fecha de los ensayos: Marzo 02 de 2001
Localización: Santa Elena, vereda El Plan
Anotaciones del lugar: Pastos sobre suelo derivado de cenizas volcánicas, pendiente de
8o.
Código de la zona homogénea: 201019
Universidad EAFIT 96
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Sitio 12.
Fecha de los ensayos: Marzo 05 de 2001
Localización: Divisoria de aguas entre las cuencas de la Picacha y Altavista (Carretera
Aguas Frías - San Antonio de Prado)
Anotaciones del lugar: Pastos sobre saprolito del Stock de Altavista, pendiente de 35o.
Código de la zona homogénea: 204008
Sitio 13.
Fecha de los ensayos: Marzo 07 de 2001.
Localización: Barrio Pajarito, cerca de la escuela Alfonso Upegui Orozco y al tanque de
EEPPM.
Anotaciones del lugar: Pastos sobre suelo residual de depósito de flujo de lodos y
escombros, pendiente de 25o.
Código de la zona homogénea: 203015
Universidad EAFIT 97
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Sitio 14.
Fecha de los ensayos: Marzo 09 de 2001.
Localización: Vereda Yolombó, entrada a la derecha 1 km después de las partidas para
San Pedro.
Anotaciones del lugar: Pastos sobre suelo residual de anfibolita, pendiente de 25o.
Código de la zona homogénea: 203006
Sitio 15.
Fecha de los ensayos: Marzo 12 de 2001.
Localización: Santa Elena, Vía a Piedras Blancas.
Anotaciones del lugar: Pastos sobre suelo derivado de cenizas volcánicas, pendiente de
8o.
Código de la zona homogénea: 201006
Universidad EAFIT 98
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Sitio 16.
Fecha de los ensayos: Marzo 14 de 2001.
Localización: Cuenca Santa Elena, via a Santa Elena en el estadero el Pingüino,
carretera al occidente. Cien metros antes de la entrada a Hogares Claret.
Anotaciones del lugar: Pastos sobre suelo residual de Dunita, pendiente de 23o.
Código de la zona homogénea: 203005
Sitio 17.
Fecha de los ensayos: Marzo 16 de 2001.
Localización: Corregimiento Santa Elena, vereda El Plan, al frente de la finca El Palmar.
Anotaciones del lugar: Pastos sobre suelo derivado de ceniza volcánica, pendiente de
10o.
Código de la zona homogénea: 201019.
Universidad EAFIT 99
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Sitio 18.
Fecha de los ensayos: Marzo 20 de 2001.
Localización: Vía a Piedras Blancas, corregimiento de Santa Elena. 20 metros después
de la tienda Pescadero.
Anotaciones del lugar: Pastos sobre derivado de cenizas volcánicas, pendiente de 10o.
Código de la zona homogénea: 201019.
Sitio 19.
Fecha de los ensayos: Marzo 23 de 2001.
Localización: Vía a Boquerón, Vereda San José de la Montaña, 30 metros arriba del
estadero Manantiales.
Anotaciones del lugar: Pastos en suelo residual de Anfibolita, pendiente de 8o.
Código de la zona homogénea: 201006
Sitio 20.
Fecha de los ensayos: Marzo 26 de 2001.
Localización: Vereda La Montañita, Corregimiento de San Antonio de Prado.
Anotaciones del lugar: Pastos sobre suelo residual de Gabro de Romeral, pendiente de
18o.
Código de la zona homogénea: 202002
Sitio 21.
Fecha de los ensayos: Marzo 27 de 2001.
Localización: Vereda La Verde, Corregimiento de San Antonio de Prado, al oriente de la
finca Hiesen y Ofelia.
Anotaciones del lugar: Pastos y árboles - pinos y eucaliptos cada 10 metros aprox.- (uso
silvopastoril) sobre suelo residual derivado de depósito de flujo, pendiente de 8o.
Código de la zona homogénea: 201014.
Sitio 22.
Fecha de los ensayos: Marzo 29 de 2001.
Localización: Vereda La Verde, Sector los Tanques, costado sur de la finca Hiesen y
Ofelia, Corregimiento de San Antonio de Prado.
Anotaciones del lugar: Pastos sobre suelo residual derivado de depósito de flujo,
pendiente de 8o.
Código de la zona homogénea: 201014.
Sitio 23.
Fecha de los ensayos: 2 de abril de 2001.
Localización: Vereda Yarumalito (por la finca Hoyo Rico), San Antonio de Prado.
Anotaciones del lugar: Bosque de pino pátula sobre suelo residual de esquistos,
pendiente 32o.
Código de la zona homogénea: 104009.
Sitio 24.
Fecha de los ensayos: 3 de abril de 2001
Localización: Corregimiento de Altavista.
Anotaciones del lugar: Suelo desnudo en Stock de Altavista, pendiente 18o.
Código de la zona homogénea: 502008
Sitio 25.
Fecha de los ensayos: 5 de abril de 2001.
Localización: Cruce vía al parque Piedras Blancas y la quebrada Piedras Blancas, cien
metros aguas abajo de la quebrada, 20 m margen izquierda.
Anotaciones del lugar: Bosque muy húmedo natural sobre ceniza volcánica (anfibolita),
pendiente de 10o. Arboles entre 10 y 12 metros de alto, varios doseles, mucha hojarasca.
Código de la zona homogénea: 101019.
Las zonas homogéneas que se tuvieron en cuenta para los ensayos fueron las que
representaban mayor proporción areal en las cuencas de estudio, tal como se menciona
en el numeral 5.3. En las otras zonas homogéneas que se encuentran sobre pastos y
cultivos se determinó el coeficiente de escorrentía partiendo de dos criterios de
extrapolación:
Los cultivos presentan un coeficiente de escorrentía con valores intermedios entre los
obtenidos en pastos y en suelos desnudos, debido a que en su mayoría los cultivos
poseen partes del suelo protegidas por la especie vegetal que se está cultivando y partes
del suelo que se encuentran desnudas.
Una zona homogénea que comparta el mismo uso del suelo y la misma formación
superficial que otra, presenta un coeficiente de escorrentía mayor en cuanto su pendiente
sea mayor.
Para las zonas homogéneas que se encuentran sin cobertura vegetal (suelo desnudo) se
utilizó la misma línea de tendencia obtenida en el ensayo efectuado en el Stock de
Altavista, cambiando sólo el intercepto.
Los resultados obtenidos hasta el momento son líneas de tendencia que relacionan la
intensidad de una lluvia con el coeficiente de escorrentía para cada una de las zonas
homogéneas que posee cada cuenca. Además se tienen intensidades de lluvia tomadas
de curvas de intensidad – duración – frecuencia de estaciones ubicadas en la zona de
estudio. Anteriormente se definieron las estaciones pluviográficas foco de cada cuenca.
Se definen las intensidades de las lluvias para cada cuenca de acuerdo con los diferentes
períodos de retorno y el tiempo de concentración de la misma (teniendo en cuenta las
curvas I-D-F de la estación foco de la cuenca).
Se construyen una tabla en Excel con las zonas homogéneas para cada cuenca en la cual
se diseña una rutina de cálculo con la fórmula de la línea de tendencia cuyo parámetro es
la intensidad para cada período de retorno. El resultado de la rutina de cálculo es el
coeficiente de escorrentía asociado a cada período de retorno para cada zona
homogénea.
Se construye una nueva tabla para cada período de retorno y se guarda como archivo de
texto. Dicha tabla contiene para cada zona homogénea un valor correspondiente al
coeficiente de escorrentía relacionado con la intensidad de la lluvia de acuerdo con el
período de retorno de la misma.
Se reclasifica la matriz que contiene las zonas homogéneas, asignándole a cada celda el
valor del coeficiente de escorrentía para cada lluvia. Este paso se desarrolla en el SIG
ArcInfo.
7 Análisis Hidrológico
El objeto de este capítulo es detallar la manera como se calcularon los caudales máximos
instantáneos producidos por lluvias críticas, determinadas por análisis de Intensidad –
Duración – Frecuencia de estaciones ubicadas en la zona de estudio. Se parte del
supuesto que lluvias críticas generan caudales máximos. Además de la obtención de
caudales máximos instantáneos siguiendo el modelo diseñado en la actual investigación,
se aplicaron modelos hidrológicos agregados con amplia utilización en nuestro medio para
poder comparar los resultados obtenidos. Por último se comparan los resultados con los
publicados en diversos trabajos y tesis en donde se presentan caudales máximos en la
desembocadura de las quebradas.
Conocida la relación entre una lluvia y la escorrentía que se genera para cada celda que
compone una zona homogénea en particular y teniendo para cada cuenca el modelo de
elevación digital, se aplican funciones de flujo propias del SIG ArcInfo para determinar
primero la dirección de flujo para cada celda del modelo y luego la acumulación de flujo en
cada celda. Si se supone una lluvia con intensidad constante durante un tiempo igual al
tiempo de concentración de la cuenca, la matriz de acumulación de flujo representa para
cada celda el caudal generado por dicha lluvia. En la investigación actual se tuvieron en
cuenta las lluvias extremas de acuerdo con curvas de Intensidad – Duración – Frecuencia
de estaciones con registros históricos en la zona de estudio, definiendo una estación foco
para cada cuenca y suponiendo que la lluvia presentaba una distribución espacial
homogénea. El coeficiente de escorrentía se determinó a partir de ensayos de campo con
la utilización de un simulador de lluvias portátil.
Para conocer el caudal que produce cada celda en una lluvia en particular se aplica el
método racional para cada celda multiplicando la matriz que contiene los coeficientes de
escorrentía por la matriz que contiene los valores de intensidad de lluvia. Es necesario
además multiplicar este resultado por el área de cada celda.
En notación del módulo Grid del sistema ArcInfo, la sintaxis del comando para la
multiplicación de matrices que representa la multiplicación es:
Grid: gqs_cu = 25 * gce_cu * gitr_cu / 36000000
Donde gqs_cu : matriz que tiene los valores de caudal de salida para cada celda.
gce_cu: matriz con los coeficientes de escorrentía para cada celda.
gtr_cu: matriz con los valores de intensidad de lluvia para cada celda.
Cabe anotar que cada lluvia analizada corresponde a una lluvia crítica determinada por
curvas IDF. Los coeficientes de escorrentía se encuentran relacionados con la intensidad
de cada lluvia de acuerdo con el análisis reseñado en el capítulo anterior. Es así como
cada lluvia con un período de retorno específico, presenta una matriz diferente de
intensidad de lluvia y de coeficiente de escorrentía.
Conocido el caudal de salida de cada celda producido por una lluvia en particular se
procede a calcular el caudal acumulado en cada celda, que siguiendo unas suposiciones
hidrológicas, representa el caudal máximo instantáneo producido por dicha lluvia. Estas
suposiciones son:
Aplicando la función flowaccumulation del módulo grid del sistema ArcInfo a la matriz que
representa el caudal de salida de cada celda se obtiene el caudal máximo instantáneo
para cada celda, la notación es:
Grid: gqa_cu = flowaccumulation (gfd_cu, gqs_cu)
El tiempo de concentración se ha definido como el tiempo que tarda una gota de agua en
llegar desde el punto más alejado de la cuenca al sitio de interés (Smith et al., 1997).
Éste es uno de los factores más importantes en el cálculo de caudales con base en
relaciones precipitación – escorrentía, pues la duración de la lluvia crítica que
posiblemente ocasionará una inundación debe suponerse igual a él, ya que para
duraciones iguales al tiempo de concentración de lluvias distribuidas sobre el área de la
cuenca toda la cuenca aporta al caudal en el punto analizado (desembocadura, estación
hidrológica, zonas críticas, etc.) para lluvias de duraciones menores al tiempo de
concentración no toda el área de la cuenca contribuye al mismo tiempo al caudal en la
salida y para duraciones mayores al mismo no hay incremento en el área contribuyente
sino un aumento en la duración de la inundación.
Teniendo en cuenta su definición, resulta claro que debe existir una relación entre el
tiempo de concentración y las características morfométricas de la cuenca, en la que se
cumple que a mayor área, mayor tiempo de concentración; a mayor longitud del cauce,
mayor tiempo de concentración; a mayor pendiente, menor tiempo de concentración, etc.
En el presente trabajo se parte de la definición física del tiempo de concentración como se
expuso anteriormente y se retoma el trabajo de Campo y Múnera (1997) en el que se
compararon las diferentes relaciones propuestas en la literatura y se determinó cuáles de
ellas se comportan mejor para el caso de las cuencas de Antioquia.
Ecuaciones recomendadas.
Según Campo y Múnera (1997), las relaciones entre el tiempo de concentración y los
parámetros morfométricos de la cuenca que mejor resultados presentan en el caso de 43
cuencas de Antioquia, aunque fueron desarrolladas para cuencas en ambientes y climas
diferentes al local, son las 4 siguientes:
Donde:
A: área de la cuenca en kilómetros cuadrados,
L: longitud del cauce principal en kilómetros,
So: diferencia de cotas entre los puntos extremos de la corriente en m/m.
AE: diferencia de cotas en kilómetros en el cauce principal desde el punto más alto al más
bajo.
Promedio aritmético 116.0 97.6 36.3 71.2 62.5 74.4 82.2 32.8 39.4 142.6
Promedio aritmético
desechando valores 117.2 96.8 36.9 71.4 64.3 75.9 84.2 32.8 38.4 143.1
extremos
Al hacer esto se está suponiendo que las estaciones mencionadas son representativas de
toda el área de las cuencas en cuestión, lo cual se aleja un poco de la realidad y es una
fuente de error en los cálculos que se deriven de estos datos. El conocimiento del
comportamiento meteorológico e hidrológico de cada cuenca se mejorará a medida que
los instrumentos instalados por el SIMPAD produzcan series de datos continuas y
confiables sobre la lluvia en las cuencas y el caudal de las quebradas.
I = C / (h + Tc)m
Donde:
C ,h y m son parámetros de la ecuación, y
Tc: tiempo de concentración en minutos.
Los parámetros de las ecuaciones para los diferentes periodos de retorno se presentan en
la tabla 24. Las curvas IDF de cada una de las estaciones mencionadas se presentan en
las figuras 75 a 80.
Tabla 24. Parámetros de las ecuaciones de las curvas IDF para la estaciones ubicadas
en la zona de estudio (tomado de comunicación escrita de EEPPM, Area de Hidrometría e
Instrumentación, abril de 2000).
Para la cuenca de la quebrada Doña María se definió la estación San Antonio de Prado
como la estación foco y para la cuenca de la quebrada La Iguaná se tuvo en cuenta la
estación San Cristóbal, ya que estas estaciones se encuentran ubicadas dentro de las
cuencas.
Las cuencas de las quebradas Altavista, La Picacha y Ana Díaz presentan a la estación
de San Antonio de Prado como estación foco y La Hueso a la estación San Cristóbal.
Para las cuencas de la zona suroccidental de la ciudad en donde no se ubican estaciones
se tenía la posibilidad de utilizar alguna de las dos estaciones mencionadas como
estación foco, al hallar las intensidades de lluvia para cada período de retorno de acuerdo
con los diversos tiempos de concentración de las cuencas se observó que los resultados
obtenidos eran muy similares teniendo en cuenta cualquiera de estas dos estaciones para
las cuencas en cuestión, comprobando que los registros históricos de lluvias máximas
para las estaciones de San Antonio de Prado y San Cristóbal son similares y pueden
representar el comportamiento temporal de lluvias máximas para las cuencas de la zona
suroccidental.
Para las cuencas de las quebradas Minitas y Quintana se tomó como estación foco a la
estación San Cristóbal. La estación de Villa Hermosa se definió como la estación foco de
la cuenca de la quebrada La Castro y la estación Vasconia como la estación foco de la
cuenca de la quebrada Santa Elena.
En la figura 81 se observa la ubicación de las estaciones con curvas IDF y los límites de
las cuencas.
Tabla 25. Intensidades de lluvia para cada cuenca con respecto a diversos períodos de
retorno.
zonas que pudieran presentar riesgos de inundación, ubicadas en zonas bajas con la
presencia de llanuras de inundación y de depósitos aluvio – torrenciales. Para cada
sección se definió la celda que surge de la intersección de la sección transversal y el
curso de la quebrada y se extrajo el valor de caudal de dicha celda de la matriz de
acumulación. Posteriormente se efectuó el análisis hidráulico de cada quebrada con la
utilización del programa HEC-RAS, el cual se explica en el próximo capítulo.
Tabla 26. Caudales máximos instantáneos en los puntos más bajos de las quebradas
analizadas, para diversos períodos de retorno, obtenidos siguiendo el método racional
distribuido.
Caudales máximos instantáneos (m3/s) para los puntos más bajos analizados de las quebradas
Tr Doña La La La Santa La
Altavista Ana Díaz La Hueso Minitas
María Picacha Iguaná Quintana Elena Castro
2.33 años 66.49 39.51 37.85 21.90 94.43 93.98 19.83 5.79 91.30 19.69
5 años 79.62 45.02 44.42 27.33 115.39 119.14 26.12 7.70 107.31 23.87
10 años 99.88 52.52 51.41 30.71 128.01 139.59 29.49 8.53 123.59 28.31
25 años 121.56 62.36 62.18 37.08 151.43 168.24 38.12 10.98 150.61 35.04
50 años 137.65 70.95 70.31 43.65 175.66 205.16 42.75 12.70 181.13 40.38
100 años 155.34 79.90 78.84 49.84 197.85 216.30 51.17 14.95 196.62 45.40
En la quebrada Santa Elena la sección transversal más baja se ubica en la entrada del
box-coulvert, aguas arriba del teatro Pablo Tobón Uribe, justo cuando se canaliza para
formar la avenida La Playa, en la carrera 38 entre las calles 51 y 54 (Caracas). La
sección más baja de la quebrada La Minitas está ubicada en la entrada del box-coulvert
en el barrio San Martín de Porres en la carrera 76 entre las calles 98 y 99. La sección
más baja de la quebrada Altavista se encuentra ubicada en la calle 30 una cuadra aguas
abajo de la carrera 65, antes de que la quebrada La Guayabala desemboque a la
quebrada Altavista.
Los datos de precipitación se tomaron de los análisis de lluvias máximas que desarrolló
Empresas Públicas de Medellín y que se presentan en las curvas de Intensidad –
Duración – Frecuencia de las figuras 75 a 80. Como hasta la finalización de la
investigación no se produjeron datos con los instrumentos de medición adquiridos, la
única otra alternativa posible fue utilizar las curvas IDF de las estaciones analizadas por
las Empresas Públicas de Medellín para los análisis actuales. A partir del procesamiento
de información cartográfica basado en la cartografía digital del proyecto SIGAME se
generaron Modelos de Elevación Digital para cada cuenca con los cuales se obtuvo los
parámetros morfométricos de las cuencas de una manera precisa.
Para el cálculo de caudales máximos con base en los parámetros morfométricos de las
cuencas se utilizaron cinco métodos diferentes: el método de las ecuaciones de
regionalización, el método racional probabilístico, el método del índice de crecientes, el
hidrograma unitario sintético del SCS y el hidrograma unitario sintético de Williams y
Hann. A continuación se presenta una breve reseña de cada uno de los métodos y una
comparación de los resultados obtenidos para cada una de las quebradas de estudio.
Este método infiere algunos aspectos de la respuesta hidrológica de una cuenca con base
en las características geomorfológicas. La cuenca se ve entonces como un filtro que
transforma la precipitación en escorrentía. Las características físicas de una cuenca se
deben a la interacción de los procesos exógenos o morfogenéticos con la geología de la
cuenca, que es un componente constante en el tiempo. A su vez, los procesos
morfogenéticos están gobernados en gran medida por la disponibilidad de agua y su
distribución temporal, por lo que se puede afirmar que las características morfométricas
de la cuenca son en parte el resultado de la interacción entre el agua y la geología de la
cuenca y por lo tanto que existe una relación determinística entre las características
morfométricas de la cuenca y las variables que describen el comportamiento hidrológico
de la misma (Smith et al., 1997).
El criterio utilizado para el caso de Antioquia por Smith et al. (1997) para la definición de
zonas de comportamiento hidrológico similar fue el de ubicación geográfica apoyado en el
mapa de isoyetas de precipitación media total multianual y el parámetro climático fue la
precipitación media total multianual, que fue hallada con base en 61 estaciones
pluviográficas y 23 estaciones pluviométricas del departamento de Antioquia y se utilizó el
método de interpolación de Krigging para la elaboración del mapa de isoyetas. Luego se
hizo ajuste a las distribuciones Gumbel y Log Normal y se encontró que se debe utilizar la
distribución de Gumbel para mayor facilidad de cálculo.
QTR = +K
Donde:
QTR : caudal máximo para un periodo de retorno TR,
: media estimada de los caudales máximos instantáneos,
: desviación estándar estimada de los caudales máximos instantáneos, y
K : factor de frecuencia que depende de la distribución y del periodo de retorno.
cada cuenca se presenta en la tabla 28. En la figura 88 se observan las estaciones con
precipitación media anual tenidas en cuenta en el presente análisis, los polígonos de
Thiessen formados y su área de influencia en cada una de las cuencas de estudio.
Tabla 27. Precipitación promedia multianual de las estaciones utilizadas en los cálculos
de caudales máximos de las quebradas del estudio (datos tomados de: EEPPM, 1990.
Anuario Hidrometeorológico).
Precipitación media
Estación
anual (mm)
San Antonio de Prado 2080.5
San Cristóbal 1703.0
Vasconia 2235.4
Miguel de Aguinaga 1397.4
Villa Hermosa 1526.5
Aeropuerto Olaya Herrera 1500.0
Tabla 28. Precipitación media anual en cada cuenca estimada por el método de los
polígonos de Thiessen y área de influencia de cada estación en cada cuenca.
2
Área de influencia de cada estación, en Km Precipitación
Cuenca San Ant. San Olaya Miguel Villa media anual
Vasconia Área total (mm)
De Prado Cristóbal Herrera Aguinaga Hermosa
Doña María 55.48 14.78 2.29 72.55 1985
Altavista 5.10 8.16 13.26 1723
La Picacha 1.27 3.82 4.98 0.59 10.66 1636
Ana Díaz 8.03 0.56 0.60 9.19 1673
La Hueso 10.27 0.99 11.26 1676
La Iguaná 48.45 2.70 51.15 1687
Minitas 0.20 1.08 1445
Quintana 0.84 1.59 1503
Santa Elena 17.72 0.39 21.34 39.45 1844
La Castro 4.22 4.22 1527
Teniendo en cuenta los datos de la tabla anterior y los parámetros morfométricos de cada
cuenca, se hallaron las características medias de los caudales máximos instantáneos para
las quebradas del estudio, en los puntos más bajos de las cuencas de drenaje
consideradas en el presente estudio, según el método de las ecuaciones de
regionalización de Smith et al. (1997) y los parámetros morfométricos utilizados en el
cálculo, los cuales se presentan en la tabla 29.
Figura 88. Polígonos de Thiessen formados por las estaciones con precipitación media
anual y su área de influencia en cada una de las cuencas.
Precipitación
área perímetro Pendiente Longitud de la Media Desviación
Cuenca (km2) (km) prom (%) corr. Ppal. (km)
media anual
Estándar
(mm)
Doña María 72.60 60.30 40.40 22.10 1985 36.63 97.07
Altavista 13.71 28.67 35.10 10.55 1723 4.75 32.88
Picacha 10.65 27.72 33.40 10.88 1636 2.79 33.22
Ana Díaz 9.17 26.56 44.26 9.28 1673 3.48 49.44
Hueso 22.24 30.73 36.01 8.18 1676 11.64 25.97
Iguaná 51.22 50.85 39.03 15.29 1687 24.50 71.22
Santa Elena 39.46 40.60 33.17 12.08 1844 18.24 31.04
La Castro 4.22 14.48 44.41 4.24 1527 2.68 20.29
Minitas 1.28 11.96 25.25 4.19 1445 0.14 5.97
Quintana 2.44 13.81 23.82 4.79 1503 0.36 5.86
El método racional fue desarrollado inicialmente por Mulvaney en 1851 y fue mejorado
después por Kuicling en 1888 y por Lloyd Davis en 1906. Aunque aparentemente es muy
sencillo, la aplicación del método se complica al estimar los coeficientes de escorrentía
correctamente. El concepto en el que se basa el método racional es que si una lluvia con
una intensidad determinada empieza en forma instantánea y continúa en forma indefinida,
el aumento de escorrentía continuará hasta que llegue el tiempo de concentración, en el
cual toda la cuenca está contribuyendo al caudal en la salida.
La entrada de agua a la cuenca está dada por el producto del área de captación por la
intensidad de la lluvia para el tiempo de concentración de la cuenca, y la relación entre
ésta y el caudal pico está dada por el coeficiente de escorrentía, que es el factor que
indica cuánta agua de esa que está entrando a la cuenca se convierte en escorrentía
(Chow, Maidment & Mays, 1994).
La fórmula general del método es:
Q = C * I * A / 3.6
Donde:
Q: caudal pico en m3/s asociado a un periodo de retorno dado, TR,
C: coeficiente de escorrentía
I: intensidad en mm/h para TR, y un tiempo de concentración Tc,
A: área de la cuenca en km2.
En este método los efectos de la lluvia y del tamaño de la cuenca son considerados
explícitamente y otros procesos, tales como infiltración, textura y estructura del suelo,
usos, etc., son considerados implícitamente en el tiempo de concentración necesario para
encontrar la intensidad de la lluvia crítica y en el coeficiente de escorrentía. No se tiene en
cuenta el almacenamiento temporal ni las variaciones espacio – temporales de la lluvia,
por lo que se recomienda aplicar este método a cuencas pequeñas, de menos de 50 km2.
determina el coeficiente de escorrentía con base en los usos del suelo y la pendiente del
terreno.
De esos mapas (figuras 5.1 a 5.7 de Smith et al., 1997) se extrajeron los valores para los
coeficientes de escorrentía de las cuencas de las quebradas objeto de estudio, los cuales
se presentan en la tabla 31.
Tabla 31. Coeficientes de escorrentía probabilísticos para las cuencas de las quebradas
de estudio (tomados de las figuras 5.1 a 5.7 de Smith et al., 1997).
Smith et al. (1997) aplicaron este método para el Departamento de Antioquia: dividieron el
Departamento en 15 zonas y obtuvieron los índices para cada una de las estaciones a
partir de los eventos generados por el modelamiento probabilístico de las series de
caudales máximos instantáneos. Calcularon la Curva Regional de Frecuencia para cada
una de las zonas del Departamento teniendo en cuenta los parámetros de escala, Q 2.33, y
área de drenaje para periodos de tiempo de base común y más largo de registros. Por
último, hicieron ajustes lineales, polinómicos, exponenciales y potenciales para cada una
de las regiones y para todo el Departamento. Las ecuaciones que cabe aplicar en este
trabajo se presentan a continuación:
Para la cuenca del río Medellín, una vez se calcula Q2.33, se procede a hallar los caudales
máximos asociados a otros periodos de retorno utilizando cualquiera de las dos
ecuaciones de las curvas regionales de frecuencia que se presentan a continuación:
Todas las cuencas del estudio pertenecen a la cuenca del Río Medellín, la cual presenta
índices característicos de 1.53 para la distribución Log Normal, 1.50 para la distribución
Gumbel por el Método de los Momentos y 1.48 para la distribución Gumbel por el método
de Máxima Verosimilitud.
Utilizando la ecuación para la cuenca del río Medellín para hallar Q2.33 con base en el área
de drenaje, se tiene (ver tabla 33):
Tabla 33. Caudales Q2.33 en m3/s para los puntos más bajos de las cuencas de drenaje
consideradas en el estudio, considerando su área de drenaje en km2 y utilizando la
ecuación para la cuenca del río Medellín: Q2.33 = -0.00006 A2 + 0.3166 A + 10.789.
Q2.33 en m3/s
Doña
Altavista Picacha Ana Díaz Hueso Iguaná Minitas Quintana St Elena Castro
María
33.46 15.12 14.15 13.69 17.80 26.85 11.19 11.56 23.19 12.12
En la tabla 34 se presentan los resultados para los caudales máximos calculados con
base en la ecuación de Q/Q2.33, mientras que en la tabla 35 se presentan los resultados
para los caudales máximos calculados en la ecuación de Q/A.
Tabla 34. Caudales máximos calculados en m3/s por el método del índice de crecientes y
la relación Q/Q2.33 = -0.4461 FX(X)2 + 1.0492 FX(X) + 0.6593, para los puntos más bajos de
las cuencas de drenaje consideradas en el estudio.
Tabla 35. Caudales máximos calculados en m3/s por el método del índice de crecientes y
la relación Q/A = -0.1648 FX(X)2 + 0.397 FX(X) + 0.1214, para los puntos más bajos de las
cuencas de drenaje consideradas en el estudio.
Los hidrogramas unitarios son modelos que suponen que el efecto de la cuenca que
transforma precipitación en escorrentía se puede representar por una función lineal de tipo
convolutiva y se utilizan para enfrentar el problema de información escasa. Un
hidrograma unitario se puede definir como el hidrograma de escorrentía directa que
resulta de una precipitación efectiva de profundidad unitaria, uniformemente distribuida
sobre la cuenca y de duración específica (Smith et al., 1997).
El cálculo del número de curva en cada cuenca se requiere para calcular la precipitación
efectiva, que es la parte de la precipitación que produce escorrentía superficial directa. El
método propuesto por el SCS se basa en el número de curva (CN), que relaciona las
características de los suelos presentes en una cuenca con sus condiciones hidrológicas,
expresada en valores de CN tabulados. Según las condiciones de humedad antecedente
se definen tres grupos: suelos secos, suelos intermedios y suelos húmedos. En el trabajo
actual se escogió los números de curva que representan la condición de suelos húmedos,
suponiendo que los eventos de crecientes se presentan en períodos invernales en donde
los suelos se encuentran en condiciones de humedad alta.
Según el método del SCS se define el CN para toda la cuenca y se propone que este
valor se obtenga como un promedio ponderado de las áreas con CN diferente. En el
presente trabajo se utilizó la información de zonas homogéneas para encontrar el CN,
relacionando cada zona homogénea con un valor de CN y encontrando su proporción
areal de acuerdo con el número de celdas que corresponden a dicha zona homogénea.
La tabla 36 muestra la forma como se obtuvo el CN de la cuenca de la quebrada La
Hueso.
Doña
Altavista Picacha Ana Díaz Hueso Iguaná Minitas Quintana St Elena Castro
María
84.57 84.16 84.94 81.32 86.26 84.73 94.10 91.81 84.68 86.07
Obtención del Hidrograma Unitario para una lluvia de duración igual a la lluvia
efectiva (según SCS).
Donde:
TR = Tiempo de rezago (tiempo en horas desde el centroide del hietograma de la
precipitación efectiva hasta el caudal pico del hidrograma unitario)
LC = Longitud del canal principal en pies
S = Factor de retención o almacenamiento en pulgadas
SC = Pendiente promedio de la cuenca en porcentaje
S = ( 1000 / CN ) - 10
Donde:
S = Factor de retención o almacenamiento en pulgadas
CN = Número de Curva
El tiempo de concentración, según el SCS, se halla a partir del tiempo de rezago según la
siguiente relación:
TC = 5/3 TR
Donde:
TC = Tiempo de concentración en horas
TR = Tiempo de rezago en horas
T = 0.133 TC
Donde:
T = Duración de la lluvia efectiva en horas
TP = (T/2) + TR
Donde:
TP = Tiempo al pico en horas
T = Duración de la lluvia efectiva en horas
TR = Tiempo de rezago en horas
UP = (2.08*AC) / (TP/10)
Donde:
UP = Caudal pico en metros cúbicos por segundo por milímetro
AC = Área de la cuenca en kilómetros cuadrados
TP = Tiempo al pico en horas
Tabla 38. Parámetros para la construcción del hidrograma unitario sintético del SCS.
La Santa La Doña
Variables La Castro Ana Díaz Altavista Minitas Quintana La Hueso
Iguaná Elena Picacha María
CN (III) 84.73 84.68 86.07 81.32 84.94 84.16 94.10 91.81 84.57 86.26
Lc 50164 39633 13904 30430 35689 34610 13737 15705 72507 26837
Sc 39.03 33.17 44.41 44.26 33.40 35.10 25.25 23.82 40.40 36.01
Ac 51.22 39.46 4.22 9.17 10.65 13.71 1.28 2.44 72.6 22.24
S 1.80 1.81 1.62 2.30 1.77 1.88 0.63 0.89 1.82 1.59
TR 1.00 0.90 0.32 0.70 0.82 0.80 0.30 0.38 1.32 0.60
Tc 1.66 1.50 0.53 1.17 1.36 1.33 0.50 0.64 2.21 0.99
T 0.22 0.20 0.07 0.16 0.18 0.18 0.07 0.08 0.29 0.13
Tp 1.11 1.00 0.36 0.78 0.91 0.89 0.33 0.43 1.47 0.66
Up 9.61 8.23 2.47 2.44 2.45 3.22 0.80 1.19 10.26 6.98
Smm 45.78 45.94 41.09 58.33 45.04 47.79 15.92 22.66 46.35 40.47
Tm 13.3 11.9 4.3 9.36 10.8 10.6 4.0 5.1 17.6 7.9
Figura 89. Hidrograma Unitario para T igual a la duración de la lluvia efectiva. Cuenca de
la quebrada La Hueso.
Obtención del Hidrograma Unitario para una lluvia de duración igual al tiempo
de concentración.
Ahora es necesario construir el hidrograma unitario para una lluvia de duración igual al
tiempo de concentración de la cuenca, para esto, se construye la curva S (suma de
hidrogramas unitarios hasta que el caudal se vuelva constante) para la lluvia efectiva y se
le resta una nueva curva S que inicia en el tiempo igual al tiempo de concentración de la
cuenca. Cabe anotar que el tiempo de concentración de la cuenca es el obtenido
utilizando las expresiones que lo relacionan con los parámetros morfométricos de la
cuenca, que para el caso de la quebrada La Hueso es de 1.19 horas. La figura 90
muestra las curvas S construidas para la cuenca de la quebrada la Hueso.
Figura 90. Curvas S para la cuenca de la quebrada La Hueso para lluvias de duración
igual a la lluvia efectiva (0.132 horas) y tiempo de concentración (1.19 horas).
Figura 91. Hidrograma unitario (SCS) para la cuenca de la quebrada La Hueso, para una
lluvia de duración 1.19 horas (igual al tiempo de concentración).
el presente trabajo se calculó la precipitación efectiva para las lluvias máximas tomadas
de las curvas IDF reseñadas con anterioridad.
Pe = (P – 0.2*S)2 / (P + 0.8*S)
Donde:
Pe = Precipitación efectiva en mm
P = Precipitación total en mm
S = Factor de retención o almacenamiento en mm
Tabla 39. Cálculo de la precipitación efectiva según el método del SCS para lluvias
máximas en la cuenca de la quebrada La Hueso.
Este mismo procedimiento se aplicó a todas las cuencas de estudio y se calcularon los
caudales máximos instantáneos para los diversos períodos de retorno, los cuales se
presentan en la tabla 40.
Tabla 40. Caudales máximos instantáneos en m3/s obtenidos según el método del
Hidrograma Sintético, según el SCS.
Las ecuaciones utilizadas en este método para determinar la forma del hidrograma
unitario sintético (ver figura 93) son:
Figura 93. Hidrograma unitario sintético de Williams y Hann (tomada de Smith et al.,
1997, figura 6.4, página 98).
t1 = t0 + 2K
Figura 94. Relación entre B y n del hidrograma unitario de Williams y Hann (tomada de
Smith et al., 1997, figura 6.4, página 100).
Universidad EAFIT 139
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
El caudal pico del hidrograma unitario, Up en pies cúbicos por segundo, se calcula con la
siguiente expresión:
Up = B * Ac / Tp
1. Se calculan K y Tp.
3. Se calcula t0 y t1.
4. Se calcula Up.
6. Se obtienen las ordenadas del hidrograma unitario sintético, Ut, reemplazando cada
abscisa t, en las ecuaciones del hidrograma, para los diferentes t. El hidrograma así
construido tiene unidades de pies3/s/pul para los caudales y de horas para los
tiempos.
Tabla 41. Parámetros y tiempos del hidrograma unitario sintético de Williams y Hann
calculados para la cuenca de la quebrada La Hueso.
Parámetro valor
K 0.6831
Tp 1.0427
n 5.9089
B 490
Up (ft3/s.in) 4082.3701
t0 (hr) 1.5134
t1 (hr) 2.8795
U0 (ft3/s.in) 2772.0604
U1 (ft3/s.in) 375.1576
Para la determinación de los hidrogramas instantáneos de Williams y Hann para una lluvia
en particular, se sigue el mismo procedimiento, explicado en el punto 7.5.4, de
convolución de la precipitación efectiva con el hidrograma unitario. Ya que la precipitación
efectiva se supuso de intensidad constante, el procedimiento se reduce a multiplicar las
ordenadas del hidrograma unitario instantáneo por la precipitación efectiva para hallar el
hidrograma instantáneo de escorrentía directa para esa lluvia.
Tabla 42. Caudales máximos para las cuencas del estudio determinados por el método
de Williams y Hann.
Tr
Caudales máximos en m3/s
(años) La Santa La La La Doña La Ana
Altavista Minitas Quintana
Iguana Elena Castro Picacha Maria Hueso Díaz
2.33 82.90 92.51 16.84 28.54 28.85 10.53 16.39 77.62 75.15 21.91
5 112.68 113.23 21.86 34.71 35.40 13.81 21.62 99.76 99.13 27.64
10 134.15 132.88 26.76 41.68 42.86 15.55 24.29 127.43 112.07 33.73
25 161.23 166.33 34.87 51.99 53.58 19.66 31.01 156.57 141.47 43.44
50 203.51 199.69 40.84 59.72 61.50 22.08 35.17 178.95 167.92 50.70
100 217.96 212.75 47.71 67.08 69.40 25.93 41.44 203.05 192.75 57.89
Tabla 43. Caudales máximos de 100 años de periodo de retorno para cada quebrada,
determinados según los diferentes métodos.
Las figuras se construyeron como diagramas de barras, graficando juntos los resultados
de los cálculos de caudales máximos asociados al periodo de retorno de 100 años. En las
leyendas, M Rac D: Método Racional Distribuido; Reg C M: Método de Regionalización de
las Características Medias de Smith et al., (1997); M Rac. P: Método Racional
Probabilístico; I C (Q/Q2.33): Método de Índice de Crecientes, utilizando la ecuación de
Q/Q2.33; I C (Q/A): Método de Índice de Crecientes, utilizando la ecuación de Q/A; SCS:
Hidrograma Sintético según el SCS y W&H: Hidrograma Sintético según Williams y Hann.
Del análisis de las figuras se puede concluir que los únicos métodos cuyos resultados
guardan algún tipo de consistencia entre sí son el Método Racional Distribuido y los
hidrogramas sintéticos, según el SCS y según Williams y Hann, ya que las diferencias
entre estos tienden a ser pequeñas en todos los casos excepto para las cuencas más
pequeñas – Minitas y La Quintana – y para la cuenca de la quebrada Doña María, cuya
diferencia es de más del 25 %.
Los métodos del SCS y Williams y Hann son ampliamente utilizados en el medio para el
diseño de estructuras hidráulicas, y con respecto a ellos, el Racional Distribuido ofrece la
ventaja de poder determinar el caudal máximo en cualquier punto de la cuenca de una
forma más eficiente, cualidad que fue aprovechada para la construcción de los mapas de
amenaza por inundación.
En este punto se presenta una comparación de los caudales máximos para algunas de las
quebradas del estudio, calculados con el método racional distribuido en la presente
investigación y con otros métodos en trabajos de investigación diferentes. El objetivo es
poner de manifiesto la variabilidad en cuanto a resultados que presentan las diferentes
metodologías utilizadas hasta el momento y de las cuales prácticamente ninguna ha
validado sus resultados debido a la escasez de información.
unitario sintético de William y Hann, los publicados por Alvarado (1994) en el “Plan de
ordenamiento territorial rural de la cuenca de la quebrada Doña María” que fueron
calculados por las firmas METROMED, siguiendo el método de regionalización, y ACP
Ltda, que hace un promedio de los resultados obtenidos por los métodos de hidrogramas
sintéticos de Snyder, Williams y Hann y SCS, y los obtenidos con el método racional
distribuido en la presente investigación.
Tabla 44. Comparación de caudales máximos determinados por Duque y Nigrinis (1992),
Alvarado (1994) y en la presente investigación para la quebrada Doña María.
Quebrada La Picacha:
Presente
Salazar y Betancur (1991)
Tr (años) Investigación
Williams y Snyder - Snyder -
H U Geom XSRAIN Rac Dist
Hann AEI Integral
5 45.8 19.9 59.9 86.4 54.6 44.4
10 56.6 29.4 70.8 108.1 69.5 51.4
25 68.2 41.2 81.9 131.7 83.3 62.2
50 81.7 56.4 95.0 159.4 99.1 70.3
100 91.1 67.8 103.7 178.9 110.1 78.8
Presente
Salazar y Betancur (1991)
Tr (años) Investigación
Williams y Snyder - Snyder -
H U Geom XSRAIN Rac Dist
Hann AEI Integral
5 52.7 27.0 61.9 95.7 70.2 27.3
10 62.1 36.4 70.3 113.6 82.4 30.7
25 80.8 56.4 86.7 149.2 105.9 37.1
50 92.7 70.2 97.1 172.1 120.6 43.7
100 104.1 84.2 106.9 194.4 134.8 49.8
Quebrada La Hueso:
Presente
Salazar y Betancur (1991)
Tr (años) Investigación
Williams y Snyder - Snyder -
H U Geom XSRAIN Rac Dist
Hann AEI Integral
5 112.5 109.2 111.7 193.9 159.8 115.4
10 138.6 141.0 132.3 239.3 193.1 128.0
25 173.6 181.5 159.8 301.0 246.1 151.4
50 193.8 210.7 175.3 336.7 273.3 175.7
100 225.4 237.2 199.7 392.6 279.8 197.9
En las tablas anteriores puede observarse que para las quebradas Doña María y Ana
Díaz, los valores obtenidos por el método racional distribuido en la presente investigación
son considerablemente más bajos que los obtenidos por otros autores, mientras que para
las quebradas La Hueso y La Picacha, los valores guardan cierta similitud.
Vale la pena reiterar que hasta el momento no es posible afirmar cuál de los métodos
representa mejor la realidad de cada cuenca, ya que la falta de datos de caudal reales y
confiables no permite hacer una evaluación ni validación de ningún método. Afirmar que
uno es mejor que otro sería una simple especulación. Hasta el momento, la ventaja de
haber aplicado el método distribuido consiste en tener disponibilidad de los valores
caudales máximos en cada punto de la cuenca de una manera más eficiente y no sólo
para ciertos puntos escogidos para hacer los cálculos.
8 Análisis hidráulico
En el presente capítulo se describe la forma como se obtuvieron los niveles de flujo en las
diferentes quebradas a partir del modelamiento hidráulico en algunos sectores de las
quebradas, utilizando el programa HEC-RAS. Se parte de cortes topográficos de las
secciones tomados inicialmente de los modelos de elevación digital y corregidos en
campo. Cada corte está definido geométricamente por su distancia a la sección
transversal aguas abajo, las abscisas que definen el cauce e hidráulicamente por los
valores de rugosidad, representados por números de Manning, tanto en el cauce como en
las bancas. El modelamiento hidráulico parte además de la información de caudal
obtenida en el análisis hidrológico. Es así como de cada quebrada se obtuvieron seis
perfiles hidráulicos, correspondientes a períodos de retorno de 2.33, 5, 10, 25, 50 y 100
años. Para determinar las condiciones de frontera, necesarias para que el programa inicie
el proceso de modelación, se supuso un régimen mixto, esto es que el flujo en cada
sección podía ser tanto subcrítico como supercrítico. Además para las secciones
extremas se supuso un flujo con profundidad normal, para lo cual el programa requiere las
pendientes del cauce.
Definidas las secciones transversales a las cuales se les iba a aplicar el modelamiento
hidráulico tendiente a encontrar los niveles de agua producidos por un caudal, se requiere
inicialmente entrar los datos geométricos al programa HEC-RAS.
Las secciones transversales a los cauces de las quebradas se tomaron a partir de los
modelos de elevación digital con la utilización del SIG ArcInfo. Se crearon coberturas con
la localización de las secciones transversales que se querían generar para cada cuenca,
las cuales pueden observarse en las figuras 82 - 87 del capítulo 7. Para generar las
secciones transversales con datos de elevación y distancia, en las corrientes principales
de cada quebrada, se utilizó el comando surfaceprofile del módulo MapTools del SIG
ArcInfo.
En dicho módulo se abre una vista nueva y se carga la superficie de elevación (MDE) en
el sistema y se le da el nombre de la cobertura que tiene las líneas por donde se quiere
El sistema genera una tabla que da la ubicación cartesiana de cada punto de muestreo,
un identificador para la superficie sobre la que se hizo el perfil, un identificador para cada
línea o sección, la elevación del punto y la distancia a partir del nodo inicial de la línea (ver
tabla 48). Dicha tabla puede graficarse luego o exportarse, desde el módulo Tables del
sistema, a un archivo ascii para su manipulación posterior en cualquier otro programa que
maneje tablas de valores, tal como Microsoft Excel o Lotus entre otros.
Tabla 48. Extracto de la tabla de valores generada por el SIG Arcinfo al hacer las
secciones transversales de la cuenca de la quebrada La Hueso (ST 01 HU).
En la mayor parte de las secciones no se detallaba la forma exacta del cauce, por lo cual
fue necesario corregir su forma a partir de medidas efectuadas en campo y de fotos
tomadas perpendiculares a las secciones y con referencias de distancia que posibilitaban
su corrección.
Posteriormente se corrige manualmente la forma del cauce ajustando los valores de altura
de las abscisas a partir de la foto. La figura 107 muestra la fotografía digital de la sección
transversal ST 01 HU con las referencias de distancia.
Figura 108. Sección transversal tomada del DEM (Cauce corregido a partir de la foto), ST
01 HU.
Además de la corrección de la forma del canal, las fotos tomadas en cada una de las
secciones fueron útiles para la determinación de los números de Manning que
caracterizan la rugosidad de las llanuras de inundación y el canal en las secciones. Estos
valores se definieron a partir de tablas que relacionan las características de la sección con
coeficientes de rugosidad. Otros datos que el programa requiere son las distancias entre
secciones tanto por el canal como por las bancas y las abscisas que definen el canal para
Distancia hasta
la próxima
sección aguas
abajo
Coeficientes
Geometría de la de Manning
sección
Abscisas que
definen el cauce
Con los datos de nivel que arroja el programa para cada sección transversal y para
diferentes períodos de retorno se construyen superficies de inundación por interpolación
de secciones tal como se reseñará en el próximo capítulo.
Tabla 49. Altura máxima de agua para diferentes períodos de retorno para las secciones
transversales definidas en las quebradas Ana Díaz y La Picacha. Datos arrojados por el
programa HEC-RAS.
Figura 112. Sección longitudinal de la quebrada La Hueso con los niveles normales y los
niveles alcanzados debido a caudales máximos instantáneos.
Figura 113. Sección longitudinal de la quebrada La Hueso entre las secciones 01 y 04,
vista en perspectiva, con los niveles normales y los niveles alcanzados debido a caudales
máximos instantáneos.
Figura 114. Curva de gasto (caudal vs elevación del nivel de agua) para la sección
transversal ST 01 HU de la quebrada La Hueso.
Universidad EAFIT 159
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
En este capítulo se presenta el proceso final, utilizado para la obtención de los mapas de
amenaza por inundación de cada una de las quebradas del estudio. Una vez se
obtuvieron los caudales máximos en cada una de las secciones transversales de las
quebradas se procedió a determinar el nivel de dicho caudal en cada sección utilizando el
programa HEC – RAS, tal como se explicó en el capítulo 8.
A partir de dicha cobertura se construyó una red triangular irregular (TIN, por sus iniciales
en inglés) haciendo interpolación lineal de los valores de nivel de agua, para así obtener
una superficie de nivel de agua asociada al caudal máximo determinado. Dicha superficie
se denominó tni100_cuenca.
Posteriormente se pasó la superficie del formato TIN a grid (o raster, matriz de valores) y
se denominó gni100_cuenca, con tamaño de celda igual al del modelo de elevación
digital (ver figura 115).
Figura 115. Proceso de construcción de las superficies de inundación a partir de los datos
de nivel de agua en cada sección transversal determinados en el programa HEC – RAS.
Para hallar las zonas inundables se restó una superficie que representa el modelo de
elevación digital de la superficie que representa el nivel de agua, gni100_cuenca, y el
mapa resultante se denominó gzi100_cuenca (ver figura 116).
Figura 116. Resta de matrices para determinar las zonas de inundación. Las celdas con
valores negativos corresponden a zonas secas y las celdas con valores positivos
corresponden a zonas que son susceptibles de ser inundadas en los eventos de caudales
máximos.
Una vez obtenidas las zonas inundables asociadas a cada quebrada se compararon los
resultados para los diferentes periodos de retorno y se decidió presentar sólo los
resultados del proceso para los caudales máximos de 100 años de periodo de retorno, ya
que las zonas de inundación son muy pequeñas y las diferencias entre estas son apenas
de algunos metros.
Luego de obtener las matrices que representan las zonas de inundación se hizo un
recorrido de campo para verificar los resultados obtenidos. En dicho recorrido se hizo una
cartografía de llanuras de inundación y de depósitos aluvio – torrenciales en las
márgenes de cada quebrada, en las zonas en las que todavía fue posible reconocerlos –
zona rural y suburbana -, en escala 1: 2000.
Los mapas preliminares de amenaza por inundación presentan estas zonas así como las
cartografiadas en campo. En la figura 117 se observa un tramo de la quebrada La Castro,
con la matriz resultante del proceso y la cobertura de polígonos que finalmente representa
las zonas de amenaza por inundación corregidas.
Los mapas de amenaza por inundación de cada cuenca se presentan en los respectivos
anexos en escala 1: 10000 en una vista general de la quebrada, y se resaltan las zonas
de amenaza en recuadros a escala más detallada, ya que en escala 1:10000 no se
visualizan bien las zonas.
Entre las 11 p.m. del viernes 16 de septiembre de 1988 y la 1 a.m. del sábado 17, se
presentó una lluvia intensa generalizada en el sector occidental de la ciudad. Dicha lluvia
ocasionó crecientes súbitas en las quebradas La Iguaná, La Hueso, Ana Díaz y La
Picacha, las cuales se desbordaron en algunos tramos, inundando parte de los barrios El
Salado, La América, La Floresta, Calasanz, El Estadio, El Corazón, Santa Teresita,
Laureles, San Joaquín, Conquistadores, San Javier, 20 de Julio, Blanquizal, La Iguaná,
Fuente Clara, Moravia y San Isidro entre otros y dejaron un saldo de 18 personas muertas
y 100 familias damnificadas (El Colombiano, 18 de septiembre de 1988).
Figura 118. Fotocopia del pluviograma de la estación San Cristóbal del viernes 16 al
domingo 19 de septiembre de 1988 (datos suministrados por EEPPM).
Universidad EAFIT 163
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
Figura 121. Barrio La Iguaná en la margen izquierda de la quebrada del mismo nombre.
La quebrada La Iguaná, a la altura de Los Colores, se llevó parte de la calzada (ver figura
119). A la altura del barrio La Iguaná (figuras 120 y 121), la quebrada arrasó varias
viviendas que se habían construido prácticamente sobre el cauce de la misma, en su
margen izquierda y afectó a todo el barrio La Iguaná.
Dicha corriente también se desbordó a la altura de los barrios Blanquizal y Fuente Clara,
pero de estos sitios no se encontró registro fotográfico.
Figura 123. Puente de la calle 35 con carrera 89 sobre la quebrada Ana Díaz.
Para evaluar dicho evento con el modelo planteado en el presente estudio, se supone que
la tormenta en cuestión tuvo una duración igual o mayor que el tiempo de concentración
de las cuencas, y que los suelos estaban saturados.
Dicho evento ocurrió en la mitad de la segunda temporada invernal del año, de una rápida
revisión de los pluviogramas y de la prensa de los días anteriores, se puede afirmar que
se cumplió la hipótesis de que los suelos estaban saturados. Por otro lado, la duración de
la lluvia fue de una hora aproximadamente, por lo que para las cuencas de las quebradas
La Hueso, Ana Díaz y La Picacha se cumple la hipótesis de que la lluvia tuvo una
duración al menos tan larga como el tiempo de concentración de las cuencas, a saber de
aproximadamente 1 hora; para la cuenca de la quebrada La Iguaná la duración de la lluvia
fue aproximadamente la mitad del tiempo de concentración de la cuenca, que es de
aproximadamente 2 horas, por lo que para esta cuenca no se presentan resultados del
modelo.
La intensidad de la tormenta que se trabajó en el modelo fue de 110 mm/hr; los caudales
máximos obtenidos en las quebradas, luego de desarrollar todo el proceso con el método
racional distribuido, se presentan en la tabla 50.
Tabla 50. Caudales máximos para una intensidad de 110 mm/hr, según el método
racional distribuido, para los puntos más bajos de las quebradas.
Figura 126. Quebrada La Picacha en el cruce con la avenida 33, zona de inundación, en
rojo, para una lluvia de 110 mm/hr y duración 1 hora.
Debido a que el modelo solo tiene en cuenta la cantidad de agua generada por una lluvia
crítica, no es posible modelar el transporte de sedimentos, que en ocasiones son los
causantes de desbordamientos producidos por taponamientos del cauce. Tal es el caso
de lo ocurrido en el puente de la carrera 86 con calle 35 (figura 123) que pasa sobre la
quebrada Ana Díaz. En este sitio el desborde de la quebrada fue producido por el
taponamiento de la sección, situación imposible de modelar teniendo en cuenta solo el
agua generada por la lluvia.
Universidad EAFIT 171
Maestría en Ciencias de la Tierra
ESTUDIO HIDROMETEOROLÓGICO Y GEOMORFOLÓGICO DE 10 QUEBRADAS DEL VALLE DE ABURRÁ
El objeto del presente capítulo es presentar los resultados de los aforos efectuados en
varias de las estaciones limnigráficas y limnimétricas instaladas durante la investigación,
se detalla el procedimiento seguido durante los aforos, los cuales deben realizarse
regularmente para determinar curvas de gasto (relación entre el nivel de la quebrada y el
caudal que pasa por la sección). Además se propone mejorar algunas estaciones en las
cuales su ubicación no permite realizar aforos confiables debido a la irregularidad e
inestabilidad del canal.
Tabla 51. Lista de estaciones para medición de nivel de agua en las quebradas,
instaladas por el SIMPAD a 27 de junio de 2001.
Todas las estaciones relacionadas en la tabla fueron visitadas en campo, los días 27 de
junio y 3 de julio del año en curso. De las visitas de campo se concluyó lo siguiente:
Las estaciones Doña María Bomba gasolina (LG) y Santa Elena Box Colvert (LG) deben
ser reubicadas según criterio del Ing. Civil Rafael Madrigal, por lo tanto no se hicieron
medidas de caudal en dichas estaciones.
La estación Hueso Metro Suramericana (ver figuras 131 y 132), debe ser limpiada
periódicamente o reubicada, ya que en este punto la sección de la quebrada cambia
constantemente con la depositación del material que arrastra la quebrada y no es posible
garantizar una relación confiable entre el nivel medido en el limnímetro y el caudal de la
quebrada. El mejor punto para reubicar esta estación es 20 metros aguas debajo de la
confluencia de las quebradas La Hueso y Ana Díaz, cerca al puente del patinódromo (ver
figura 5) ya que el canal en ese punto es trapezoidal y no se presenta depositación de
ningún tipo.
Las estaciones Doña María Pte Vereda La Verde (ver figuras 133 y 134) y Castro Hogar
Antioquia (ver figura 136) deben ser limpiadas antes de hacer las mediciones de caudal vs
nivel de agua, ya que las piedras que arrastra la quebrada cambian sustancialmente la
geometría de la sección y de esta forma no se podría encontrar una curva de gastos para
dichas estaciones en la que se pueda relacionar caudal vs nivel de agua con un grado de
confianza aceptable.
Figura 130. Fotografía del limnímetro de la estación Iguaná Ladrillera San Cristóbal. En
ésta puede observarse la pequeña caída de agua que se formó hacia la parte inferior de
la foto y las rocas y el pequeño empozamiento hacia la parte central y en el borde de la
quebrada, evidencias de la inestabilidad de la sección.
Figura 132. Vista del limnímetro de la estación Hueso Metro Suramericana. Obsérvese la
diferencia entre el nivel del agua en el limnímetro y en el centro del canal (parte inferior de
la foto), en este caso la lectura del limnímetro no corresponde al nivel real del agua.
Figura 133. Sitio adecuado para la reubicación de la estación Hueso Metro Suramericana,
en el sitio indicado por la flecha roja, a aproximadamente 20 metros aguas abajo de la
confluencia de las quebradas Hueso y Ana Díaz.
Figura 134. Estación Doña María Vereda La Verde. Las piedras arrastradas por la
quebrada formaron una pequeña caída de agua en el sector donde se ubicó el limnímetro
(círculo azul).
Figura 135. Detalle del limnímetro de la estación Doña María Vereda La Verde.
Obsérvese la diferencia en el nivel de agua de la quebrada en la zona inferior de la foto y
el señalado en el limnímetro.
Figura 136. Estación Castro Hogar Antioquia. El limnímetro (ver círculo azul) está ubicado
justo debajo del puente peatonal, a pocos centímetros de las piedras que se observan en
la foto y que hacen que el flujo de agua sea nulo en ese lado del canal.
En primer lugar se levanta la sección transversal detallada del sitio en donde se encuentra
ubicado el limnímetro o limnígrafo, tomando datos de profundidad del cauce cada 50 cm.
Para esto se ubican dos jalones a ambos lados del cauce, y se amarra una cuerda entre
ellos (ver figura 137).
Figura 137. Ubicación de los jalones en una línea perpendicular al cauce, con la cuerda
amarrada entre ellos. Estación Iguaná Carrera 74, sitio con limnímetro y limnígrafo.
Luego se marca la cuerda con cinta negra cada 50 cm (ver figura 138), para determinar
los puntos en los que se va a medir la profundidad del cauce.
Las mediciones de profundidad del cauce se toman con una mira de topografía
convencional en los sitios señalados en la cuerda.
Figura 138. Detalle de la cuerda (color amarillo) amarrada entre los jalones con marcas
cada 50 cm (negras, ver recuadro azul). A la izquierda puede observarse el sitio donde
está ubicado el limnímetro (barra roja).
Figura 139. Esquema de una sección transversal, dibujada con datos de profundidad
cada 50 cm (h1, h2) y los puntos donde se toma el dato de velocidad (v).
En los canales en “V” de las quebradas del estudio, ya que definen cauces muy angostos
(de aproximadamente 1 a 1.5 metros de ancho), se tomaron datos de velocidad en tres
tramos, tal como se ilustra en la figura 11, en los puntos v1, v2 y v3.
Para determinar el caudal de agua que pasa a través del área delimitada por cada
trapecio, se hallan las áreas de los mismos delimitados por las profundidades medidas (h1
y h2 en las figuras 138 y 140) y se multiplica por el valor de la velocidad medido dentro de
cada trapecio. Por último, se suman todos los caudales medidos en la sección para
determinar el caudal total.
Esta estación está ubicada en la calle 30, al frente de la entrada al Almacén General de
EEPPM (ver figuras 142 y 143).
Figura 143. Estación Altavista Almacén EEPPM, vista desde el puente peatonal.
Esta estación está ubicada en el puente peatonal a 100 metros de la carrera 65 hacia la
desembocadura de la quebrada, en el barrio Conquistadores (ver figuras 144 y 145).
Figura 144. Estación Picacha Conquistadores, vista desde aguas arriba, margen
izquierda.
Figura 145. Estación Picacha Conquistadores, vista desde aguas arriba, margen derecha.
Esta estación está ubicada en la estación Floresta del Metro, a pocos metros aguas abajo
de la carrera 80 (ver figura 146).
Está ubicada cerca al puente peatonal al frente del Colegio Jesús María (ver figuras 147 y
148).
Figura 147. Estación Ana Díaz Colegio Jesús María, visto desde la margen derecha.
Figura 148. Estación Ana Díaz Colegio Jesús María, al fondo se observa el puente
peatonal.
Figura 149. Estación Quintana Carrera 65, vista desde la margen izquierda.
Está ubicada a aproximadamente 50 metros del puente de la carrera 74 (ver figuras 151 y
152). En esta estación de aforo se encuentran instalados un limnímetro y un limnígrafo.
Figura 151. Estación Iguaná Carrera 74, al fondo se observa el puente de la carrera 74.
Figura 152. Detalle de la realización del aforo en la estación Iguaná Carrera 74, la
medición de la velocidad del agua se hizo cada 50 cm.
Esta estación está ubicada en la mitad del tramo de la quebrada Santa Elena
comprendido entre la glorieta de Fatelares y la Autopista, al lado de la Plaza Minorista (ver
figura 153).
Figura 153. Estación Santa Elena Plaza Minorista. El limnímetro se encuentra en la zona
demarcada por el círculo azul.
El acceso a esta parte de la quebrada para la realización del aforo es muy difícil, por lo
que se recomienda mejorarlo mediante la instalación de barras de metal o madera, a
modo de escalas, aprovechando los orificios que se encuentran en el muro de concreto
(ver figura 154).
Orificios en el muro
de concreto
Figura 154. Detalle del limnímetro de la estación Santa Elena Plaza Minorista.
12 Conclusiones y recomendaciones.
Los Modelos de Elevación Digital (MDEs) son la base para el análisis hidrológico que
se efectuó, pues a partir de éstos fue posible generar la red de drenaje, la divisoria de
aguas de la cuenca y mediante un procesamiento de la información, las curvas
hipsométricas de las cuencas y los mapas de pendientes de las mismas. Otra
aplicación importante de los MEDs es la generación de matrices de dirección y
acumulación de flujo que permitieron definir caudales instantáneos de todos los puntos
de las cuencas a partir de lluvias espacializadas y zonas homogéneas en cuanto a la
relación precipitación – escorrentía.
Las curvas hipsométricas de cada cuenca del presente estudio son un aporte nuevo al
conocimiento del entorno físico de la ciudad. De éstas es posible extraer información
valiosa acerca de la morfología de las cuencas, y su importancia y validez radica en
que son una forma objetiva de presentación de la información de la distribución del
relieve de la cuenca, que no está sesgada por las ideas a priori que pueda tener el
investigador.
Para las cuencas de las quebradas del estudio, los tiempos de concentración varían
entre media hora (Minitas) y 2 horas 20 minutos (Doña María) aproximadamente, lo
que indica que el tiempo de respuesta ante eventos en estas cuencas es muy bajo, y
por lo tanto el esfuerzo para disminuir los desastres en ellas debe enfocarse en las
labores de prevención.
Los caudales calculados con el método racional presentan un aumento menos abrupto
que los calculados con el método de regionalización de las características medias a
medida que aumenta el periodo de retorno, pero hasta el momento no es posible
determinar cuál de estos métodos se acerca más a la realidad de cada quebrada,
pues no existen datos de caudal que sirvan para validar los resultados. Una vez se
tengan series de datos largas y confiables de estas quebradas, será posible
determinar cuáles métodos se ajustan más a la realidad de cada cuenca.
Del análisis de las figuras en las que se compara los caudales máximos para 100 años
de periodo de retorno se puede concluir que los únicos métodos cuyos resultados
guardan algún tipo de consistencia entre sí son el Método Racional Distribuido y los
hidrogramas sintéticos, según el SCS y según Williams y Hann, ya que las diferencias
entre estos tienden a ser pequeñas en todos los casos excepto para las cuencas más
pequeñas – Minitas y La Quintana – y para la cuenca de la quebrada Doña María,
cuya diferencia es aproximadamente del 25 %. Es importante aclarar que no se está
afirmando aquí que el Método Racional Distribuido de estimaciones más cercanas a la
realidad que los otros métodos, sino que los métodos mencionados dan resultados
similares entre sí.
Los métodos del SCS y Williams y Hann son ampliamente utilizados en el medio para
el diseño de estructuras hidráulicas, y con respecto a ellos, el Racional Distribuido
ofrece la ventaja de poder determinar el caudal máximo en cualquier punto de la
cuenca, cualidad que fue aprovechada para la construcción de los mapas de amenaza
por inundación.
pero en casi todas las quebradas se observa que son radicalmente diferentes de los
obtenidos con los demás métodos. Al comparar los resultados para las quebradas Ana
Díaz y La Hueso se observa una inconsistencia: el caudal calculado para la quebrada
Ana Díaz es más alto que para la quebrada La Hueso, lo cual no es posible en la
realidad, ya que la primera es afluente de la segunda. Por lo tanto, se recomienda no
usar este método para el cálculo de caudales máximos en el Valle de Aburrá.
Los mapas presentados fueron el resultado de aplicar una lluvia máxima de periodo de
retorno de 100 años, de acuerdo con los datos disponibles de las curvas IDF de las
estaciones San Cristóbal, San Antonio de Prado, Villa Hermosa, Olaya Herrera y
Miguel de Aguinaga.
Para que el modelamiento hidráulico sea más refinado es necesario levantar más
secciones transversales al cauce de las quebradas y entrarle datos de caudal y
niveles de agua a medida que se avance en el monitoreo hidrológico de las cuencas.
13 Bibliografía
Angel, J. & Duque, T., 1990. Zonificación de riesgo para la cuenca de la Quebrada
Iguaná. Tesis de grado. Universidad Nacional de Colombia, Sede Medellín, Facultad
de Minas.
Chow, V. T.; Maidment, D. R. & Mays, L. W. 1988. Applied Hydrology. McGraw Hill
Book Co. New York. 574 p.
Chow, V. T., 1982. Hidráulica de los Canales Abiertos. Editorial Diana, México. 620
p.
Duque, J. & Nigrinis, O., 1992. Estudio hidrológico e hidráulico de la quebrada Doña
María. Tesis de Grado, Universidad Nacional de Colombia, Sede Medellín, Facultad
de Minas.
Harden, C. 1988. Meoscale estimation of soli erosion in the Rio Ambato drainage,
Ecuadorian Sierra. En: Mountain Research and Development, Vol. 8. No. 4, pp. 331 –
341.
Harden, C., 1992. A new look at soil erosion processes on hillslopes in highland
Ecuador. En: Erosion, Debris Flows and Environment in Mountain Regions. IAHS Publ.
No. 209. pp. 77 – 85.
Hutchinson, M. F., 1989. A new procedure for gridding elevation and stream line data
with automatic removal of spurious pits. Journal of Hydrology:106, 211-232
Salazar, J. & Betancur, J., 1991. Estudio hidrológico e hidráulico de las quebradas
La Hueso, La Ana Díaz, La Picacha. Tesis de grado. Universidad Nacional de
Colombia, Sede Medellín, Facultad de Minas.
Salinas, I. C., González, R., Vargas, J. P., Jiménez, J. L., Chica, G & Bustamante,
M., 1990. Las Formaciones Geológicas Superficiales de Medellín. En: AGID Report
No. 13, Medellín, p. 503 – 512.
Smith, R.; Vélez, M. V.; Acosta, Z. L.; Blandón, J. R.; Posada, J.E.; Botero, A.L.;
Campo, J.M.; Montoya, R.D.; Múnera, J.C., Muñoz, S.I.; & Ruiz, C.D., 1997.
Hidrología de Antioquia. Ed. Imprenta Departamental de Antioquia. 176 p, anexos.
Toro, G. & Velásquez, A., 1984. Estudio geomorfológico del Valle de Aburrá desde
Caldas hasta las quebradas Santa Elena e Iguaná. Proyecto de grado. Facultad de
Minas, Universidad Nacional de Colombia. Medellín.
Anexos.
Se presentan a continuación los resultados obtenidos para cada cuenca de los análisis
cartográficos, los mapas temáticos trabajados para la modelación hidrológica y los mapas
preliminares de amenazas.
Además se incluye una breve reseña de las características principales de cada cuenca,
profundizando en los usos actuales del suelo. Un análisis de la distribución espacial de
los usos del suelo es importante puesto que el cambio en dichos usos permitirá en un
futuro la disminución de las zonas de amenaza por inundación, específicamente, el
cambio de las amplias zonas que en la actualidad se encuentran en pastos por bosques.
Dicho cambio disminuirá los coeficientes de escorrentía con la consecuente disminución
de caudales máximos instantáneos.
Adicionalmente se incluye un anexo digital con el registro fotográfico de las cuencas, las
secciones transversales con sus respectivas fotografías, los mapas temáticos y de
amenazas en archivo digital, y los resultados de los ensayos efectuados con el simulador
de lluvias.
Limita por el norte con la cuenca de la quebrada San Francisco, por el oriente con las
cuencas de las quebradas La Hueso, La Picacha, Altavista, Jabalcona, por el sur con la
cuenca de la quebrada Grande y por el occidente con las cuencas de las quebradas la
Frisola y La Sucia. Hacia el sur del Cerro del Padre Amaya nacen dos afluentes de la
quebrada Doña María: quebrada La Frisola en la cota 2905 msnm y la quebrada Los
Vallados en la cota 2590 msnm y a partir de la confluencia de estas en la cota 2240
msnm, toma el nombre de quebrada Doña María. Desemboca en el río Medellín en la
cota 1522 msnm luego de un recorrido de 22.1 Km. El área total de la cuenca es de 72.6
Km2, con un total de 48 microcuencas en la margen derecha y 43 microcuencas en la
margen izquierda.
Distribución espacial de los usos del suelo actuales. Cuenca de la quebrada Doña
María.
La cuenca de la quebrada Altavista, limita por el norte con la cuenca de La Picacha, por el
sur con la cuenca de La Guayabala, por el occidente con la cuenca de la quebrada Doña
María y por el oriente con el río Medellín. Posee en su cabecera dos ramales: quebrada
Barcino, que nace en las estribaciones del Alto El Encanto a 2340 msnm, con una longitud
de 2.7 Km y quebrada Buga, con nacimiento en el cerro El Barcino en la cota 2270 msnm
con una longitud de 2.4 Km; ambas quebradas confluyen en la cota 1720 msnm donde
toman el nombre de quebrada Altavista, y luego de 7.2 Km tributa sus aguas al río
Medellín en la cota 1475 msnm, al pie del cerro Nutibara.
Son peculiares los drenajes largos en su margen derecha, con un patrón de drenaje
subparalelo dendrítico controlado por el diaclasamiento con sentido NS y NE. En su
margen izquierda, la quebrada Santa Rita es el único drenaje largo.
La parte baja comprende la zona plana correspondiente a la llanura aluvial, donde las
pendientes varían entre 3% y 15%.
Se observan zonas en rastrojo alto y bajo a lo largo del cauce, café asociado con plátano
como sombrío y a plena exposición en sentido de la pendiente y hortalizas y pastos de
corte en áreas pequeñas. En su parte alta se han desarrollado pequeños núcleos
urbanos, correspondientes a viviendas de interés social que están en proceso de
densificación. En la parte media predominan las canteras, las ladrilleras y los tejares, que
deterioran las terrazas debido a la explotación de material con técnicas artesanales y
rudimentarias.
Los procesos erosivos predominantes son los deslizamientos en roca saprolitizada con
formación de surcos y cárcavas y erosión fluvial por socavación en las márgenes.
En la parte alta el cañón es profundo en forma de "V" con vertientes empinadas que
evidencian pequeños derrumbes y degradación del suelo, con sobrepastoreo intenso en
las laderas de la margen derecha; otros usos son las plantaciones de café sin sombrío y
asociadas con plátano. En la cabecera y a lo largo del lecho se observan rastrojos altos y
explotación de material pétreo en el lecho. Como procesos erosivos se mencionan surcos
y cárcavas en cortes de talud y pequeñas corrientes de agua y socavación lateral del
cauce.
En su parte alta presenta un cañón profundo en forma de "V", cauce rectilíneo desde el
nacimiento hasta la confluencia con la quebrada Romedala, que sugiere un control
estructural, posiblemente un lineamiento el cual se prolongaría en la cuenca de la
quebrada La Hueso a lo largo de su afluente El Salado. El cauce se encuentra sobre
buena cobertura vegetal de regeneración espontánea, el lecho es pedregoso con aguas
limpias. Se observa sobrepastoreo intenso en laderas y plantaciones de ciprés sobre
cuchillas y laderas.
Distribución espacial de los usos del suelo actuales. Cuenca de la quebrada Ana
Díaz.
Su parte alta se caracteriza por poseer un cauce rectilíneo de cañón poco profundo, el
cual se va profundizando en forma de “V”, volviéndose además sinuoso a medida que
desciende. Presenta vertientes escarpadas con pendientes de hasta el 60% que se
suavizan en las demás zonas hasta el 30%. Se anota la presencia de socavación en las
márgenes desestabilizando taludes. Los retiros de la quebrada se encuentran cubiertos
por vegetación de rastrojo altos y bajos.
La parte media que comprende las cotas 1700 y 1500 msnm posee pendientes
variables entre el 30 y el 15%. Las vertientes se han desarrollado sobre el Batolito
de Altavista altamente fracturado y meteorizado, lo que permite su utilización
como material pétreo. Se reportan movimientos en masa debido a la alta
meteorización y sobresaturación de agua en los suelos, reptación, deslizamientos
antiguos y recientes causados principalmente por acción antrópica (canteras y
banqueos). Se encuentra modificada debido a la urbanización intensa que
incluso ha invadido sus retiros en ciertas zonas.
En la parte baja se localizan los depósitos aluviales y posee rangos de pendiente entre
3% y 15%, allí el desarrollo ha transformado por completo los cauces naturales por medio
de canalizaciones. Se reseñan dos eventos recientes con altos efectos negativos: el 20
de octubre de 1984, cuando se presentó un fuerte aguacero con una duración de 90
minutos y una precipitación de 56 mm, según los datos de la estación pluviométrica de
San Cristóbal perteneciente a Empresas Públicas de Medellín. La gran cantidad de
escorrentía generada sobrepasó la capacidad de las obras de canalización, causando
desbordamiento y daños en los barrios Santa Lucía, Antonio Nariño y El Socorro (Correa,
Montoya y otros, 1984, citados por Mi Río, Municipio de Medellín, PNUD, 1995).
Posteriormente, hubo otro evento de grandes proporciones donde se desbordaron varias
quebradas de la zona suroccidental en el año de 1988.
Limita al norte con la cuenca de la quebrada Malpaso, al noroccidente con la Cuchilla Las
Baldías, divisoria de aguas con la cuenca de la quebrada El Hato, al occidente con la
divisoria de aguas de la cordillera El Frisol y la quebrada La Frisola hasta el Cerro del
Padre Amaya, al suroccidente con el Alto El Astillero y la Cuchilla Los Arrayanes, por el
sur con la cuenca de la quebrada La Hueso y por el oriente con el río Medellín. Nace en
la cota 2555 msnm en las estribaciones del Alto de las Repetidoras, después de recorrer
en dirección SE un tramo aproximado de 5.6 Km cambia a dirección WE hasta su
desembocadura en el río Medellín, en la cota 1457 msnm. Posee una longitud total de
15.29 Km y un área de 51.22 Km2, de los cuales 6.35 Km2 corresponden a la zona urbana
y 44.87 Km2 a la zona rural.
Su parte alta se encuentra por encima de la cota 1820 msnm, se caracteriza por estar
constituida en la vertiente sur por esquistos en su parte más alta con alto grado de
fracturamiento y meteorización, con zonas escarpadas y con alta susceptibilidad a la
erosión donde son frecuentes los movimientos en masa. Se destacan los escarpes del
Cerro del Padre Amaya y el Boquerón. El resto de la vertiente sur se encuentra sobre
saprolito del Batolito de Altavista, con dedicación al pastoreo, cultivos de flores y
hortalizas. Por su parte, la vertiente norte está desarrollada sobre saprolito de anfibolita
en superficies inclinadas y onduladas, como el Alto de Yolombó en la Cuchilla Las
Baldías. Se reportan numerosos deslizamientos superficiales y profundos, flujos, caídas
de roca, reptación, surcos y cárcavas.
La parte media se ubica entre las cotas 1515 msnm y 1820 msnm. La vertiente norte está
desarrollada sobre extensos depósitos de vertiente compuestos en esencia por anfibolita
saprolitizada, presentando superficies cóncavas y convexas, debido a intercalaciones de
flujos de lodos y escombros. La vertiente sur se desarrolla sobre un cuerpo neisico
altamente meteorizado.
La parte baja se define desde la cota 1515 msnm hasta la desembocadura. Esta zona
completamente urbanizada, corresponde a un abanico aluvial en donde se destacan dos
cerros denominados Indural y El Volador. Sobre esta zona se han presentado eventos de
inundación de alta magnitud.
Limita al sur y occidente con las cuencas de las quebradas Malpaso y El Malpaso, la
cuenca de la quebrada La Cantera al norte y el río Medellín por el oriente. Nace en la
cota 2207 msnm y desemboca en la cota 1445 msnm. La cuenca se caracteriza por
poseer un patrón de drenaje subparalelo con tendencia noreste y su forma rectangular –
oblonga.
La cuenca tiene una forma alargada que evidencia una expresión morfológica en una
serie de ondulaciones que constituyen una pendiente homogénea debido al material de
depósitos de flujos saprolitizados que conforman la zona. La mayor parte de su recorrido
lo hace en zona urbana por lo cual se encuentra altamente intervenida.
La parte alta se caracteriza por la presencia de lecho pedregoso en los tramos naturales
con inestabilidad y socavación. En la parte media su cauce ha sido modificado y
transformado por medio de coberturas y en la parte baja se encuentra un tramo natural
con retiros conservados en pastos y rastrojos bajos, de cauce estrecho y aguas muy
contaminadas.
Limita al norte con las cuencas de las quebradas El Ahorcado, El Molino y parte de
Piedras Blancas, por el oriente con la cuchilla La Gulupera, al Sur con el municipio de
Envigado y las cuencas de las quebradas La Presidente y La Poblada y al occidente con
el río Medellín. El área total de la cuenca es de 45.61 Km2 de los cuales el 74.95%
pertenece a zona rural y el 25.05% a zona urbana, el estudio hidrológico de la cuenca se
efectuó hasta la entrada del box-coulvert, uno metros arriba del teatro Pablo Tobón Uribe,
con un área de captación de 39.46 Km2. Nace en el Cerro Espíritu Santo en la cota 2620
msnm y luego de un recorrido de 14.43 Km desemboca en el río Medellín en la cota 1457
msnm. La longitud del cauce hasta la entrada del box-coulvert es de 12.08 Km a una
altura de 1505 msnm.
La parte alta entre las cotas 2720 y 2250 corresponde a una antigua superficie de erosión
y posee pendientes entre el 20% y el 30%. Las rocas presentan una cobertura de cenizas
volcánicas de un poco más de un metro de espesor con excepción de áreas mínimas
afectadas por erosión intensa debido a procesos antrópicos. La topografía es suave y los
cañones de la quebrada son uniformes. En general, la tasa de erosión es baja (James,
1988).
La parte media alta corresponde a la franja ubicada entre las cotas 2520 msnm y 1850
msnm, caracterizada por ser la parte más abrupta de la cuenca. Esta zona escarpada
presenta planos de diaclasamiento muy marcado sobre la vertiente norte que favorecen
los desplazamientos de bloques de roca, la anfibolita se presenta como una roca fresca
recubierta por cuñas de depósitos de ladera. La vertiente sur presenta serpentinita de
grano fino y fracturada. Las altas pendientes y precipitación favorecen los eventos
torrenciales y facilitan los movimientos en masa tales como los deslizamientos
rotacionales y las caídas de bloques de roca.
La parte media abarca desde la cota 1850 msnm hasta la 1540 msnm con un rango de
pendiente entre 30% y 40%. Sobre las laderas se localizan numerosos depósitos de flujos
de lodo y escombros meteorizados que llegan hasta los depósitos aluviales y torrenciales.
Se reporta reptación y movimientos de masa debido a la filtración de aguas subterráneas.
La parte baja comprende desde la cota 1540 msnm hasta la desembocadura de la Santa
Elena en el río Medellín. Existen depósitos torrenciales, terrazas y conos de deyección.
Esta zona completamente urbanizada a perdido su expresión geomorfológica debido a la
acción antrópica.
Distribución espacial de los usos del suelo actuales. Cuenca de la quebrada Santa
Elena.
Este afluente de la quebrada Santa Elena presenta una longitud de 4.24 Km. Su cuenca
hace parte de la vereda Piedras Blancas y los barrios Las Estancias y Villa Tina,
comprendiendo un área de 4.22 Km2, de los cuales 4.06 Km2 pertenecen a zona rural y
0.16 Km2 a zona urbana. Nace en la cota 2500 msnm y desemboca en la cota 1630
msnm.
En su parte alta presenta un cañón amplio con pendientes mayores al 60% sobre roca
fracturada, el lecho es pedregoso con rocas angulares. En este sector predomina la
vegetación arbustiva de rastrojo bajo con tramos reforestados. En la parte media el cañón
es profundo e incisado de taludes casi verticales, los cuales forman una garganta
angosta. Sus retiros se encuentran con frecuencia invadidos por asentamientos
subnormales. En la parte baja el cañón se amplia un poco entre cerros de poca altura, al
mismo tiempo que desarrolla pequeñas terrazas donde se estableció el barrio Las
Estancias.