MODELAMIENTO DE LA QUEBRADA “LA HUESO”

Angela Sofia Gaviria

Estudiante ingeniería civil
Universidad de Antioquia
Daniela Vásquez Puerta
Estudiante ingeniería civil
Universidad de Antioquia
Simón Alejandro Gallo
Estudiante ingeniería sanitaria
Universidad de Antioquia

Abstract:

This report studies and analyzes the micro-watershed of the La Hueso stream in Medellín,
Antioquia, with the objective of modeling it, analyzing the potential risks of flooding, and
establishing preventive measures. The HEC-RAS software was used to conduct this stream
study. The area was visited, and photographs were taken, along with a visual inspection of the
different structures present in the section. These collected data are necessary for accurately
modeling the stream in the HEC-RAS software. Different models were created to evaluate the
behavior of the stream under various flow conditions. Special attention was given to the
structures in the area, including box culverts, bridges, and piles, as they can generate scour
effects in the channel. Based on the obtained results, an analysis was conducted for each return
time and flow regime. The passage of floodwater in a selected section was assessed, and the
effects of the flow on the different identified structures in the study area were analyzed.

Keywords: Stream La Hueso, HEC-RAS, flow, flood, box culvert, pile, bridge, scour.

Resumen:
Este informe estudia y analiza la microcuenca del arroyo La Hueso en Medellín, Antioquia,
con el objetivo de modelar, analizar los riesgos potenciales de inundación y establecer medidas
preventivas. Se utilizó el software HEC-RAS para llevar a cabo este estudio del arroyo. Se
visitó el área y se tomaron fotografías, junto con una inspección visual de las diferentes
estructuras presentes en la sección. Estos datos recopilados son necesarios para modelar con
precisión el arroyo en el software HEC-RAS. Se crearon diferentes modelos para evaluar el
comportamiento del arroyo bajo diversas condiciones de flujo. Se prestó especial atención a las
estructuras en el área, incluyendo alcantarillas, puentes y pilotes, ya que pueden generar efectos
de socavación en el canal. Basándonos en los resultados obtenidos, se realizó un análisis para
cada período de retorno y régimen de flujo. Se evaluó el paso del agua de inundación en una
sección seleccionada y se analizaron los efectos del flujo en las diferentes estructuras
identificadas en el área de estudio.

Palabras clave: Quebrada La Hueso, HEC – RAS, Caudal, inundación, flujo, box culvert, pila,
puente, socavación.
Introducción.
En el diseño de canales hay diferentes
factores de gran importancia que van
encaminados a la simulación de un flujo
para determinar su nivel de acuerdo a
variables como la pendiente, forma del
cauce y aquello que pueda tener influencia
en este. Todo esto con el objetivo de
realizar estudios como las probabilidades
de inundación, socavación, y en general, de
toda condición que se pueda presentar con
el tiempo y la influencia climática, con la
finalidad de prevenir diferentes daños tanto Ilustración 1: Tramo de estudio. Fuente: Google
Earth.
para las obras hidráulicas, como para las
La quebrada La Hueso en el tramo de
comunidades aledañas.
estudio se encuentra canalizado con
El modelamiento se realiza mediante el
secciones transversales rectangulares, y con
software HEC-RAS y girará en torno al
secciones trapezoidales en la parte más
estudio de la quebrada La Hueso en
cercana a la desembocadura como se
Medellín, la cual se constituye como la más
observa en las ilustraciones 2 y 3.
importante de la comuna de San Javier,
atravesando sectores como San Cristóbal,
Laureles y Estadio.

Descripción del tramo de estudio.

La quebrada La Hueso se encuentra
localizada en la ciudad de Medellín. Su
nacimiento se encuentra en el barrio San
Javier, ubicado a 2160 msnm, y con una
longitud aproximada de 8,3 Kilómetros
desemboca en el río Medellín a una altura
de 1483 msnm a un costado del centro de
espectáculos La Macarena. La estación de
monitoreo se ubica en el barrio Estadio, en Ilustración 2: Sección trapezoidal aguas abajo de
la zona centro occidental del municipio de la avenida regional. Fuente propia.

Medellín, en la Carrera 65 con Calle 48 A

(Pichincha) a un costado de la estación
Suramericana del Metro.
El tramo de estudio se localiza desde Av. 33
#66b-12 hasta la desembocadura en el río
Medellín por la avenida Regional, dicho
tramo tiene una longitud aproximada de
1.045 Kilómetros. Ilustración 3: Sección rectangular entre la calle 48
y 47D. Fuente propia.
Además, en el tramo de interés, se También es importante mencionar que en la
encuentran diferentes obras físicas como sección analizada, las pendientes de los
pilas de puentes y box culvert, que márgenes del canal de la corriente de agua
interfieren directamente con el cauce de la están habitadas por individuos que, en
quebrada. A continuación, se presenta un repetidas ocasiones, arrojan residuos
inventario de algunas de las obras en este sólidos al río.
tramo:
● Una pila ubicada a escasos metros
de la Av. Regional que se encuentra
en el centro del cauce de la
quebrada. Ilustración 4.

Ilustración 6: Desechos arrojados a la quebrada.

Fuente propia.

HEC-RAS
Inicialmente se trazó el cauce del río a partir
del boceto obtenido en Google Earth para
dar una representación cercana del área de
estudio.
En la imagen (7) se muestra el modelo
digital desarrollado en HEC-RAS con las
secciones transversales agregadas a lo largo
de la quebrada.

Ilustración 4: Pila cerca de la avenida Regional.

Fuente propia.

● Puente abajo de la estación

Suramericana del metro.

Ilustración 5: Puente. Fuente propia. Ilustración 7: Sección del tramo de estudio.

Elaboración propia.
Modelado del tramo en HEC-RAS.
Para realizar la edición de las secciones se
debe hacer desde la opción “View/Edit
geometric date” , en la cual se
desplegará una nueva ventana, en esta se
accede al comando “ Edit and/ or create
cross section” , al hacer clic en dicha
opción se desplegará una ventana con toda
la información de la sección transversal,
entre las cuales se encuentra “Station” y
“Elevation”, cada dato de estos dos
Ilustración 8: Sección trapezoidal. Elaboración
representa un punto que define la sección propia.
transversal, entonces modificando estos dos
valores se puede cambiar la ubicación de
cada punto y así cambiar la sección Tabla 2: Medidas de sección transversal
rectangular.
transversal. También es de anotar que entre
Sección transversal (6 a 16)
los datos que se pueden modificar en esta
ventana se encuentra el n de Manning, el Forma del canal Rectangular

cual es de 0.015 para todos los lados tanto Ancho (B) 13 m

en la sección rectangular como la Altura 2,4 m
trapezoidal. Diferencia de elevación 0,3 m
Pendiente 0,006 m
Se tomaron distancias entre secciones cada
50 m iniciando en estadio con un canal
trapezoidal, posteriormente, alrededor de
unos 200 m el canal toma una forma de
sección rectangular con las medidas
especificadas en las siguientes tablas.

Tabla 1: Medidas de sección transversal

trapezoidal.
Sección transversal (1 a 5)
Forma del canal Trapezoidal
Ancho (B) 16 m
Ancho (b) 5m
Altura 2,2 m
Ilustración 9: Sección rectangular. Elaboración
Diferencia de elevación 0,2 m propia.
Pendiente 0,004 m
 Tabla 3: Medidas de sección transversal A la hora realizar de las obras hidráulicas
rectangular. que se desean incluir en el modelo, se debe
Sección transversal (17 a 20) seleccionar el ícono en el editor de
Forma del canal Rectangular geometría, para así crear las pilas y los
Ancho (B) 12 m culverts necesarios y realizar su pertinente
Altura 1,4 m modelado de flujo para diferentes
condiciones.
Diferencia de elevación 0,5 m
Pendiente 0,01 m

Ilustración 10: Sección rectangular. Elaboración

propia.

Es importante resaltar que, para las

secciones con forma rectangular, se tienen
los mismos valores del n de Manning
correspondientes a un valor de 0,015, ya
que este tipo de obra tiene el mismo
material tanto en los lados como en la base
del canal. Esto se observa en la Ilustración
11.

Ilustración 12: Pila en la sección rectangular y

Ilustración 11: n de Manning para las secciones. trapezoidal respectivamente.
 Después de tener las secciones
transversales y haber agregado la pila y el
box culvert, se procede a simular el
comportamiento de la quebrada con un
flujo estacionario, el cual se realiza desde el
comando “View/Edit steady flow data”
En el cual se ingresan el dato de caudal para
un determinado tiempo de retorno, estos
datos fueron tomados de: (La red de
monitoreo ambiental en la cuenca
hidrográfica del río Medellín). Lo anterior
Ilustración 13: Culvert entre suramericana y se puede observar en la ilustración 15.
desembocadura.

Para el culvert representado hace referencia

al culvert de la Carrera 65, de
Ilustración 15: Caudales máximos instantáneos.
aproximadamente 13 metros de ancho. En Fuente: Red de Monitoreo en la cuenca
la representación se tuvo en cuenta que el hidrográfica del rio Medellín.
pilar se encontraba en el cauce del río así
que este se someterá posteriormente a También es necesario ingresar condiciones
análisis de socavación. Para construir la de frontera, para este caso la profundidad
carretera para el culvert se ingresaron los crítica, como se observa en los perfiles de
siguientes datos mostrados en la ilustración flujo para diferentes tiempos de retorno.
14, incluyendo posición y ancho de la obra
hidráulica.

Ilustración 16.1: Perfil de flujo para tiempo de

retorno de 2.33 años. Elaboración propia

Ilustración 14: Datos para la carretera del culvert.

Elaboración propia.
Además de los datos mostrados
anteriormente se debe ingresar en la opción
“Htab Param.” la altura máxima que puede
alcanzar el agua, en este caso se ingresó la
cota correspondiente al gálibo del puente,
para así posteriormente realizar la
Ilustración 16.2: Perfil de flujo para tiempo de
representación del puente de la Carrera 65
retorno de 5 años. Elaboración propia
que pasa por la quebrada La Hueso.
 Luego de ingresar todos los datos
anteriores, se procede a modelar el flujo
estacionario, para esto se debe crear un
plan, para esto se ingresa a la pestaña
“Perform a steady flow simulation” , al
ingresar en dicha opción, aparecerá una
nueva ventana.
Al final de la modelación se puede observar
Ilustración 16.3: Perfil de flujo para tiempo de el nivel del agua en las diferentes secciones
retorno de 10 años. Elaboración propia para los periodos de retorno ingresados,
como se observa en la ilustración 17.

Ilustración 16.4: Perfil de flujo para tiempo de Ilustración 17: Flujo estacionario en una sección
retorno de 25 años. Elaboración propia trapezoidal. Elaboración propia.

Además, se puede observar el

comportamiento del flujo en 3D, esto
permitiendo dar una mejor idea de cómo es
el comportamiento a lo largo de la zona de
interés y así poder observar las zonas más
críticas, como se muestra a continuación:

Ilustración 16.5: Perfil de flujo para tiempo de

retorno de 50 años. Elaboración propia

Ilustración 16.6: Perfil de flujo para tiempo de

retorno de 100 años. Elaboración propia.

Ilustración 18: Modelo 3D para flujo estacionario.

Elaboración propia.
 Ilustración 19: Tabla de datos para el Tr de 100
años. Elaboración propia.
Ilustración 21: Datos para Flow Hydrograph.
Elaboración propia.
Posteriormente para realizar la simulación
donde se presenta un flujo inestable y Luego por efectos del estudio de la
posteriormente un desbordamiento con modelación se decidió tomar datos con una
respecto a las condiciones iniciales del variación un poco uniforme y así obtener la
canal, se debe ingresar los datos de un siguiente gráfica de distribución:
histograma de frecuencias aguas arriba, y
tomamos el flujo aguas abajo como una
profundidad normal, datos los cuales serán
mostrados a continuación:

Ilustración 22: Curva de distribución para los

datos tomados. Elaboración propia.
Una vez ingresados los datos se procede a
Ilustración 20:: Datos para flujo inestable. correr la simulación propuesta donde
Elaboración propia. inicialmente presenta errores a la hora de
A continuación, se ingresan los datos con hacer correr el programa, unos con fácil
una fecha fija donde se supone que empieza solución y otros más complejos, los errores
la simulación con intervalos de 10 minutos más comunes y fáciles de solucionar son los
entre cada dato, como se muestra en la asociados a el cambio de región y zona
siguiente ilustración: horaria, ya está se debe modificar desde la
configuración y cambiar la zona horaria a
Hora Central (Estados Unidos Canadá), y
en región se debe seleccionar Estados
Unidos.

Por otro lado, se presentan errores más

complejos y difíciles de hallar, donde se
deben modificar factores de la geometría en
la sección que presenta problemas, muchos
de estos se solucionan como se menciona toman datos del tamaño y tipo de grano
anteriormente, configurando el “HTab presente en el sedimento cargado por la
param” de las obras hidráulicas, y de las quebrada.
secciones en general, en la opción de editar

geometría se selecciona el icono y allí

se debe variar los puntos y el incremento,
de manera que el flujo quede por debajo de
estas.
Luego de verificar los posibles errores se
procede a correr la simulación de flujo
inestable hasta que no se muestre ningún
error como se muestra en la siguiente
ilustración:

Ilustración 23: Verificación de simulación de flujo

Ilustración 24: Socavación. Fuente: elaboración
inestable. Elaboración propia.
propia.

Por último, se debe determinar la

Análisis de resultados.
socavación presente en el flujo. Esta se
Para las condiciones de flujo estacionario,
presenta en la pila del puente de la carrera
se deduce que no se presenta
65 y en los bordes del canal. Para hallar la
desbordamientos para caudales de periodos
socavación, se tomó el caudal para un
de retorno de 2.33 y 5 años, esto es bueno
periodo de ocurrencia de 2.33 años, es
ya que indica que las obras de canalización
decir, de 11.67 𝑚3 /𝑠; esto debido a que la
en el tramo de estudio se hicieron pensando
pila no es alcanzada por el flujo. Para este
en poder soportar altos caudales que
cálculo, en el software se utilizó el
corresponden también a altos periodos de
comando “Perform hydraulic design retorno.
computations”. También es claro las diferencias de
velocidad entre diferentes tramos de la zona
La socavación se hace presente igualmente de interés, en las cuales se observa una
en los bordes del canal, debido a que se mayor velocidad aguas arriba y a medida
encuentra canalizado en el tramo estudiado. que avanza el flujo se va reduciendo esta
También, debido a la canalización, no se velocidad. Esto es coherente con la realidad
ya que a lo largo de este tramo se cotas del terreno, ya que las
encuentran diferentes disipadores de secciones transversales no
energía, los cuales hacen que se reduzca la coinciden con las del canal ya que
velocidad del flujo en la quebrada. este se encuentra canalizado, por lo
que toca modificar ciertas secciones
En el flujo inestable se presentan unas para estas se adecúen más a la
condiciones donde hay un desbordamiento realidad del canal.
del flujo en el terreno, condición la cual ● El N de Manning que se supone
hace que el programa presente una cantidad corresponde a secciones de la
considerable de errores, que deben ser quebrada que se encuentran
corregidos para poder realizar el modelo de canalizadas, por lo que se conoce un
simulación. En el tramo que seleccionamos, aproximado de la superficie.
normalmente las condiciones de flujo son ● El canal presenta una geometría
estables ya que aguas arriba hay un homogénea, además de un flujo que
vertedero que sirve como disipador de se distribuye a lo ancho del canal, ya
energía, sin embargo, con las características que no se presentan barras
de la simulación hacen que en nuestro hidráulicas que eviten el flujo del
tramo se presente un desbordamiento, agua.
hecho que en la realidad sería desastroso ya ● El inicio del tramo presenta flujo
que en las inmediaciones de la quebrada supercrítico, cuando se encuentra
canalizada, se encuentran un gran número con pilas cambia a subcrítico debido
de viviendas. a que estas estructuras disipan la
energía del flujo.
En lo que respecta a socavación, esta existe
en el pilar del puente de la carrera 65 y en
los bordes de la quebrada, ya que se Referencias
encuentra canalizada. ● https://es.scribd.com/document/475
588116/Quebrada-La-Hueso#
Conclusiones. ● https://www.metropol.gov.co/ambi
● El tramo de estudio podría soportar ental/recurso-hidrico/PORH/fase-
una caudal correspondiente a un 4-formulacion/programa-de-
periodo de retorno de 5 años, esto, monitoreo.pdf
aunque puede dar tranquilidad, no ● https://silo.tips/download/red-de-
es algo seguro, ya que los caudales monitoreo-ambiental-en-la-cuenca-
de periodos de retorno se basan en hidrografica-del-rio-aburra-
probabilidades, por lo cual podría medellin-en
ser posible que en cualquier ● https://www.corantioquia.gov.co/ci
momento aumente tanto el caudal adoc/AGUA/AIRNR_CN_1658_1
que cause un desbordamiento de la 999.pdf
quebrada.
● A la hora de importar el modelo de
elevación digital, para así obtener
una aproximación de la geometría,
solo es útil para tener una idea de las