Science & Mathematics > Earth Sciences">
Informe Hidrologico e Hidraulico Canal La Sierrita
Informe Hidrologico e Hidraulico Canal La Sierrita
Informe Hidrologico e Hidraulico Canal La Sierrita
Preparado por
TABLA DE CONTENIDO
Pág.
TABLA DE CONTENIDO 2
LISTA DE TABLAS 4
LISTA DE FIGURAS 5
LISTA DE ANEXOS 6
1. GENERALIDADES 7
2. OBJETIVOS 8
3. ESTUDIO HIDROLÓGICO 9
4.2.1 Parteaguas 10
4.2.2 Área A 10
4.2.3 Perímetro P 10
4.2.4 Coeficiente de Compacidad Kc 11
4.2.5 Factor de forma Kf 11
5. DISEÑO HIDRÁULICO 29
8. CONCLUSIONES 36
9. RECOMENDACIONES 37
BIBLIOGRAFÍA 38
ANEXOS 39
LISTA DE TABLAS
Pág.
LISTA DE FIGURAS
Pág.
LISTA DE ANEXOS
Pág.
1. GENERALIDADES
En este informe se presentarán los resultados del estudio hidrológico y diseño hidráulico de 1465
metros de canalización de una corriente tributaria a la quebrada El Tupe en el municipio de
Chiriguaná – Cesar, localizado en el corregimiento la Sierra entre El Cruce y Rincón Hondo.
El estudio consiste en el análisis hidrológico de las áreas aferentes y diseño hidráulico de las obras
de conducción de aguas lluvias que afectan el sitio en consideración. Para el análisis hidrológico
se revisaron los datos climatológicos disponibles en la zona, se calculó la precipitación de diseño
considerando los tiempos de retorno recomendados por el Reglamento Técnico del Sector de
Agua Potable y Saneamiento Básico RAS y con el uso de métodos desarrollados para modelar el
proceso de transformación lluvia – escorrentía, se calculó el caudal máximo probable que
circularía por el tramo de corriente a canalizar. Para el diseño hidráulico se tuvo en cuenta la
geomorfología de la zona, las condiciones actuales de la corriente y estructuras existentes y se
realizó el cálculo hidráulico con los caudales obtenidos en el análisis hidrológico, finalmente se
obtuvieron las secciones hidráulicas recomendadas en cada caso.
Para lograr un diseño económico se subdividió el área de la cuenca en áreas más pequeñas con el
ánimo de obtener caudales en función de la longitud de la canalización planteada, de este modo
se obtuvo un primer tramo de 745 m con ancho de base de 4.00 m y altura variable de 0.70 m
hasta 1.1 m, y un segundo tramo de 740 m con ancho de base de 5.00 m y altura variable de 1.10
m hasta 1.80 m, el cual entrega sus aguas al boxcolvert existente de 5.00 x 1.70 el cual cumple los
requerimientos hidráulicos del proyecto y por tanto no requiere ser remplazado desde el punto
de vista hidráulico.
2. OBJETIVOS
Pronosticar los caudales máximos probables para el chequeo y diseño de las obras hidráulicas
presentes en el sitio del proyecto.
Ejecutar el chequeo y diseño de las secciones hidráulicas y cotas máximas que se requieren
para evacuar satisfactoriamente el caudal máximo probable de las estructuras requeridas en
el proyecto.
3. ESTUDIO HIDROLÓGICO
La hoya hidrográfica que aporta escorrentía al sitio de interés del proyecto abarca un área total de
0.79 Km2. Se ubica en la zona central del departamento del Cesar a los 9° 22’ 9.557’’ de latitud
norte y 75° 31’ 28.346’’ de longitud oeste, ver figura 1, es surcada por la corriente de interés en este
estudio. Para efectos de seleccionar secciones del canal acorde a las condiciones del lugar, la hoya
hidrográfica fue subdividida en 6 unidades más pequeñas tal como se ve en la figura 2 y en la
tabla 2.
4.2.1 Parteaguas. Conocido también como divisoria de aguas, se designa como la línea que separa
las precipitaciones que caen en hoyas hidrográficas inmediatamente vecinas, y que encaminan la
escorrentía resultante para uno u otro sistema fluvial [2].
4.2.2 Área A. Se define como la proyección horizontal de la superficie de drenaje que va dirigida
a un cauce y que está delimitada por la divisoria de aguas o parteaguas [2].
4.2.3 Perímetro P. Es la longitud sobre el plano horizontal, que recorre la línea divisoria de aguas,
dicho de otro modo, corresponde a la longitud del polígono cuya proyección horizontal, contiene
el área de la hoya hidrográfica.
=
2 ⁄
es el perímetro de la cuenca [Km]
es el área de drenaje de la hoya [Km2]
Cuanto más irregular sea la cuenca mayor será su coeficiente de compacidad. Una hoya circular
posee el coeficiente mínimo, igual a uno. Hay mayor tendencia a las crecientes en la medida en
que este número sea próximo a la unidad [2].
Valores de Kc Forma
1.00 - 1.25 Redonda a oval redonda
1.25 - 1.50 Oval redonda a oval oblonga
1.50 -1.75 Oval oblonga a rectangular oblonga
4.2.5 Factor de forma Kf. Se determina a partir del cociente entre el área de la cuenca A, y el
cuadrado de la longitud hidráulica de la hoya hidrográfica.
=
es el área de drenaje de la hoya [Km2]
es la longitud hidráulica o curso de agua más largo de la hoya hidrográfica [Km]
Si > 1 la hoya es achatada y tiende a formar fácilmente grandes crecidas; un < 1 indica que
la hoya es alargada, menos sujeta a crecientes que una de la misma área y mayor factor de forma;
por otro lado si = 1 la hoya tiene forma de un cuadrado [2].
A continuación se muestran las propiedades morfométricas de las subhoyas del área de estudio.
En cuanto al relieve de una hoya hidrográfica, se pueden definir varios parámetros topográficos
que controlan la velocidad con que se da la escorrentía superficial y por lo tanto el tiempo de
concentración. Entre estos parámetros se tiene la elevación máxima y mínima, la pendiente media,
elevación medía. Como se observa a continuación, la zona de estudio es plana.
Tabla 3. Parámetros del relieve de las subhoyas hidrográficas del área de estudio
Subhoya Elevación máxima [m] Elevación mínima [m] Pendiente media [%]
Sierra 1 57.00 44.67 2.04
Sierra 2 57.00 46.62 2.04
Sierra 3 57.00 47.44 2.36
Sierra 4 53.00 48.00 2.02
Sierra 5 53.00 48.20 2.43
Sierra 6 47.00 45.10 1.95
La red de drenaje se puede definir como la trayectoria, disposición o arreglo de los cauces y lechos
por donde de manera superficial y aparente corre el agua excedente, producto de la precipitación
hacia un depósito natural o artificial [3].
4.4.1 Cauce principal. Es la corriente que pasa por la salida de la misma. Toda cuenca tiene una
y solo una corriente principal.
4.4.2 Pendiente del cauce principal. La pendiente del cauce principal se calcula a partir del
cociente entre el desnivel entre los puntos definidos por las cotas que marcan la desembocadura
y el nacimiento del cauce principal y la longitud del mismo. Este parámetro es indispensable para
determinar el tiempo de concentración de una cuenca e influye directamente en la velocidad del
agua que fluye por el cauce principal [4].
4.4.3 Cota mayor y menor del cauce principal. Corresponde a la elevación más alta y más baja
del cauce principal respectivamente.
4.4.4 Densidad de drenaje Dd. Corresponde a la relación entre la longitud de drenaje y el área
de la hoya hidrográfica, se calcula con la siguiente expresión.
4.4.5 Longitud hidráulica Lm. Corresponde a la mayor trayectoria de flujo desde el parte aguas
de la hoya hasta el punto de salida de la misma.
A continuación se muestras los parámetros de drenaje de las subhoyas del área estudiada.
La referencia [5] lo define como el tiempo para el cual toda la hoya hidrográfica empieza a
contribuir, es decir, el tiempo necesario para que toda la hoya aporte eficazmente a la generación
de flujo en el sitio de desagüe.
Una hipótesis ampliamente usada en hidrología es suponer que la duración de la lluvia de diseño
coincide con el tiempo de concentración de la hoya hidrográfica, con lo cual se maximiza el caudal
de diseño. Sin embargo, [6] establece una diferencia importante al expresar que el tiempo de
concentración es la fase del aguacero que determina el máximo caudal, y esta fase puede ser igual
o no a la duración del aguacero, ya que la duración de la lluvia es ajena al tamaño de la cuenca y
es solo función de variables climáticas.
Para su cálculo se pueden emplear diferentes fórmulas que se relacionan con otros parámetros
propios de la cuenca. Para la estimación del tiempo de concentración usualmente se recomienda
emplear el promedio obtenido de varias ecuaciones empíricas disponibles en la literatura
científica, sin embargo, en este estudio se usó la mediana de los datos debido a que esta no es
susceptible a los valores extremos.
A continuación se muestran las ecuaciones usadas para el cálculo del tiempo de concentración en
este proyecto.
.
= 0.95
Bransby Williams.
= 0.605 . .
×
Kirpich.
.
= 0.0663 .
.
= 0.3 .
Pasini.
( × ) .
= 0.108 .
.
= 0.871
.
= 0.191 .
California Culvert
Bransby Williams
Mediana [horas]
Corps Engineers
Reclamation
Bureau of
[minutos]
Mediana
Subhoya
practice
Kirpich
Témez
Pasini
Sierra 1 0.68 1.44 1.04 0.75 2.45 0.68 0.38 0.75 44.77
Sierra 2 0.59 1.26 0.94 0.64 2.17 0.59 0.34 0.64 38.23
Sierra 3 0.29 0.78 0.67 0.37 1.64 0.29 0.22 0.37 22.31
Sierra 4 0.20 0.53 0.50 0.23 1.19 0.20 0.16 0.23 14.07
Sierra 5 0.10 0.24 0.17 0.09 0.31 0.10 0.08 0.10 6.28
Sierra 6 0.13 0.23 0.16 0.08 0.22 0.13 0.07 0.13 7.54
∑ 280562.165
= = = 0.3551
∑ 790141.0511
Para el análisis del comportamiento climatológico en la zona se analizaron los datos de la estación
meteorológica Apto las Flores que pertenece a la red de estaciones del IDEAM. En la tabla 2, se
muestran los valores medios mensuales de precipitación, temperatura medía máxima,
temperatura media mínima y temperatura media que fueron extraídas de las series de tiempo de
35 años suministrados por el IDEAM (ver anexos).
Fuente. Instituto de hidrología, Meteorología y Estudios Ambientales, modificado por los autores.
En la tabla anterior se aprecia que la temperatura media anual es de 28.8 °C, y que oscila entre
una media mínima de 21.6 °C y una media máxima de 36.1 °C, siendo los meses de enero, febrero,
marzo y abril los que reportan las temperaturas más altas, durante los meses de mayo, junio, julio
y agosto, se estabiliza cerca de la media con la llegada de las lluvias y finalmente en septiembre,
octubre, noviembre y diciembre se aprecian las temperaturas más bajas, ver figura 5.
La precipitación media anual en la zona de estudio es del orden de los 1849 mm, siendo abril,
mayo y junio, agosto, septiembre, octubre y noviembre los meses más lluviosos, y julio el mes en
que se presenta una pequeña disminución de la precipitación, lo cual le confiere un
comportamiento bimodal. El periodo de escasez de lluvias lo conforman los meses de diciembre,
enero, febrero y marzo, ver figura 4.
4.8.1 Periodo o tiempo de retorno. Según [5] el periodo de retorno puede definirse como el
intervalo de recurrencia promedio entre eventos que igualan o exceden una magnitud especifica.
La probabilidad de ocurrencia del evento ≥ en cualquier observación puede relacionarse con
el periodo de retorno de la siguiente manera.
1 1
= =
( ≥ ) 1− ( ≤ )
[6] establece en el numeral D 4.3.2 que se deben adoptar los siguientes periodos de retorno para
el cálculo de las lluvias de diseño.
Fuente. Reglamento técnico del sector de agua potable y saneamiento básico Título D.
4.8.2 Precipitación máxima probable. [5] la define como la mayor profundidad de precipitación
estimada analíticamente para una duración dada que sea físicamente posible y que caracterice
razonablemente una región geográfica particular en un determinado periodo del año.
Dicho procedimiento se hace con el uso de distribuciones de probabilidad [5]. En este caso
particular se utilizó la distribución de probabilidad de valores extremos tipo I o distribución de
Gumbel, la cual es una distribución de probabilidad utilizada usualmente para valores máximos
aleatorios extraídos de poblaciones suficientemente grandes [7], en hidrología se recomienda un
mínimo de 25 años de datos. Se trabajó con los datos de precipitación máxima en 24 horas de la
estación pluviométrica Rincon Hondo por ser la mas cercana a la zona de estudio (ver anexos).
( )= −∞ ≤ ≤ ∞
√
= ; = ̅ − 0.5772
−
=
1
= − ln ln
( )
1 1
= =
1− ( ≤ ) 1− ( )
De donde
−1
( )=
Por lo cual
= − ln ln
−1
= +
Fuente. Instituto de hidrología, Meteorología y Estudios Ambientales, modificado por los autores.
Tabla 10. Parámetros para el cálculo de la distribución de Gumbel, estación Rincón Hondo
Parametro Valor
No. Datos 35
Media 109.8571
Desviación estándar 31.6052
Coeficiente de asimetría 0.0534
α 24.6425
u 95.6335
Nivel de confianza β 90.00%
Nivel de significancia α 5.00%
Probabilidad acumulada 95.00%
Zα 1.6450
Tabla 11. Distribución de Gumbel para P máx. en 24 horas, estación Rincón Hondo
Como se observa en la tabla anterior, los datos cumplen con la prueba de bondad y ajuste K – S y
por tal motivo pueden ser usados bajo la distribución de valores extremos tipo I (Gumbel)
4.8.4 Curvas Intensidad – Duración – Frecuencia. [7] las define como arreglos en los cuales se
presentan las lluvias (estimadas como intensidad de precipitación) contra su duración y el periodo
de retorno. Se obtienen a partir de datos históricos de precipitación de una estación provista de
pluviógrafo. Cuando el nivel de complejidad del sistema sea bajo, [6] propone el uso del método
simplificado para la obtención de las curvas IDF. El procedimiento consiste en usar la siguiente
expresión
=
60
Donde:
Corregimiento de
La Sierra – Cesár
Región a b c d
R1 0.94 0.18 0.66 0.83
R2 24.86 0.22 0.50 0.10
R3 13.92 0.19 0.58 0.20
R4 5.53 0.17 0.63 0.42
En la figura 15 se observa que la zona de estudio se encuentra en la región R1, por lo tanto, los
parámetros de regresión y datos para determinar las curvas IDF son los siguientes:
REGION a b c d M*
R1 0.94 0.18 0.66 0.83 109.857
* M corresponde a la media de la muestra presentada en la tabla 10
Aplicando lo anterior se obtiene la siguiente expresión para calcular las curvas IDF de la zona de
estudio
. .
0.94 109.857
= .
60
.
692.779.12
= .
Donde:
La referencia [6] indica que para obtener los caudales de diseño de sistemas de recolección y
evacuación de aguas lluvias de áreas relativamente pequeñas se puede usar el método racional
siempre y cuando la hoya hidrográfica sea menor a 80 ha.
El área de interés en este estudio es de 79 ha, y fue subdividida como se vio en la figura 2 para
efectos de escoger una sección relativamente pequeña en el inicio del canal e ir aumentándola
gradualmente. Por lo anteriormente expuesto se utilizó el método racional para el cálculo de los
caudales de diseño. El procedimiento consistió en suponer que la duración de la lluvia es igual al
tiempo de concentración de cada hoya hidrográfica, luego de eso se es escogió un periodo de
retorno de diseño de 25 según el criterio de los autores a pesar de que [6] recomienda 10, lo
anterior a raíz de las fuertes variaciones climáticas que se han presentado recientemente en todo
el territorio nacional, con esos datos se obtuvo la intensidad de la lluvia para cada hoya
hidrográfica a partir de las curvas IDF de la figura 16, finalmente se aplicó la siguiente ecuación
correspondiente al método racional.
= 0.278 × ×
Donde
4.9.3 Caudales hidrológicos resultantes. Los caudales máximos probables obtenidos con el
método racional aplicando el procedimiento ya descrito, constituyen la base para el diseño
hidráulico que será abordado en el siguiente capítulo. A continuación los resultados obtenidos
del modelo hidrológico.
Se aclara que el abscisado del canal se realizó dejando la abscisa 0+000 en la sección aguas abajo
del boxculvert existente y la abscisa 1+485 en la sección aguas arriba del canal, lo anterior debido
a que ese es el modo de funcionamiento del modelo hidráulico HEC –RAS el cual fue usado en el
presente diseño. También se aclara que la canalización será realizada desde el 0+020 hasta el 1+485
con lo cual se tienes 1465 m de canalización. El 0+020 es el punto de entrada del boxculver existente
en la vía Chiriguana – Rincón hondo.
5. DISEÑO HIDRÁULICO
El diseño hidráulico consiste en proponer una sección hidráulica con base en las condiciones
encontradas en el sitio del proyecto y luego hacer el chequeo correspondiente con el uso de un
modelo hidráulico. En este caso particular se utilizó el modelo HEC-RAS, desarrollado por cuerpo
de ingenieros del ejército de los Estados Unidos, ver [8], el cual permite realizar cálculos de flujo
uniforme unidimensional, flujo inestable unidimensional y bidimensional, transporte de
sedimentos en lecho móvil y modelos de temperatura y calidad del agua.
Para un análisis hidráulico básico se requiere de la topografía del canal (secciones transversales,
perfil longitudinal), de la rugosidad de la superficie del terreno (coeficiente de rugosidad de
Manning) y las condiciones de frontera que se encargan de definir el comportamiento del modelo
en la primera sección (aguas arriba) y en la última (aguas abajo). Cuando se modela con flujo
uniforme estas condiciones pueden ser: un nivel de aguas en alguna sección de control, el tirante
critico definido por las condiciones de la sección hidráulica, el tirante normal que está en función
de la sección transversal y la pendiente del canal y, finalmente una curva de calibración de
caudales si se tiene.
5.1.2. Condiciones de frontera. Las condiciones de frontera se definieron tanto aguas arriba como
aguas abajo para realizar un análisis de flujo en régimen mixto (normal y crítico). Aguas arriba se
usó la condición del tirante normal con una pendiente de 0.002 m/m y en la sección aguas abajo
se usó el tirante normal con una pendiente de 0.002 m/m, ambos datos extraídos del diseño
preliminar de la rasante de la base o solera del canal.
5.1.3 Situación actual del sitio del proyecto. La longitud total del canal es de 1465 metros. Para
efectos de crear un adecuado modelo hidráulico se incluyó un tramo de 20 metros aguas abajo del
sitio donde comenzara la canalización y a ese punto se le asigno como nombre “K0+000”, ver
figura 10. A lo largo del tributario a canalizar existen 3 estructuras hidráulicas, la primera es un
boxculvert de 15 metros de longitud y sección de 5.00 m x 1.70 m que recibe las aguas del canal
del proyecto la cual será conservada en caso de cumplir con los requisitos hidráulicos, las dos
siguientes son alcantarillas circulares de 38” de diámetro localizadas en las abscisas K0+310 y
K0+990 del modelo hidráulico las cuales deben ser reemplazadas debido a que su sección es
insuficiente para los caudales del proyecto y las pendientes que se deberían usar para no
demolerlas serian inferiores al 0.05% lo cual traería problemas de colmatación y tirantes altos. A
continuación se muestran las secciones propuestas para realizar el chequeo hidráulico del canal.
Sección hidráulica
Tramo del canal
Inicial Final
K1+485 hasta K0+760 4.00 x 0.70 4.00 x 1.10
K0+760 hasta K0+740 Transición de 4.00 a 5.00
K0+740 hasta K0+020 5.00x1.10 5.00 x 1.80
5.1.4 Chequeo del tramo K1+485 hasta K1+130. El caudal para el chequeo de estas secciones es
de 3.68 m3/S. A continuación se muestran los gráficos resultantes del modelo en los dos tramos.
Como se observa en la imagen anterior, la sección con ancho de base de 4.00 es satisfactoria para
drenar el caudal de diseño preservando un buen margen de borde libre.
5.1.5 Chequeo del tramo K1+130 hasta K0+740. El caudal para el chequeo de este tramo es de 7.57
m3/S. A continuación se muestran los gráficos resultantes del modelo.
Como se observa en la gráfica anterior, el boxculvert es satisfactorio para evacuar las aguas lluvias.
5.1.6 Chequeo del tramo K0+740 hasta K0+290. El caudal para el chequeo de este tramo es de
12.54 m3/S. A continuación se muestran los gráficos resultantes del modelo.
Al igual que el tramo anterior, en este también existe una alcantarilla circular de 38” de diámetro
sobre una de las vías de la zona de estudio, por tal motivo se planteó la construcción de un
boxculvert de 6.00 m x 1.80 m para evacuar las aguas de drenaje remplazando la estructura
existente, a continuación se muestran la solución planteada.
5.1.7 Chequeo del tramo K0+290 hasta K0+000. El caudal para el chequeo de este tramo es de
16.12 m3/S. A continuación se muestran los gráficos resultantes del modelo.
La sección K0+020 del modelo hidráulico coincide con la entrega de la canalización proyectada al
boxculvert existente. A continuación se muestran los gráficos resultantes para el chequeo de dicha
estructura.
Como se observa, la sección existente evacua sin problemas el caudal de diseño por lo cual no
debe ser remplazada desde el punto de vista hidráulico.
5.1.8 Chequeo de canal tributario en la abscisa K0+025. Este tramo corresponde a 20 metros de
canal que entrega sus aguas en la abscisa K0+025 de la canalización principal, el caudal de diseño
es de 0.64 m3/S. A pesar de que el caudal es pequeño, este tramo de canal estará sujeto a sobre
elevaciones causadas por las profundidades de flujo que se presenten en el canal principal objeto
de este estudio. A continuación se muestran las gráficas del modelo hidráulico para este tramo
que inicia con una sección de 3.00 m x 1.00 m y finaliza con 4.00 m x 1.00 m.
Figura 19. Tributario de la abscisa K0+025, perfil con secciones inicial y final.
5.1.9 Chequeo de velocidades de flujo. Con el ánimo de evitar daños en la estructura, las
velocidades de flujo en el canal deben ser inferiores a 8 m/S si el revestimiento se realiza en
concreto de 210 kg/cm2 [7], a continuación se muestran la velocidades de flujo encontradas en el
proyecto.
Como se observa, en el caso critico la velocidad máxima es de 4.0 m/S lo cual es satisfactorio.
8. CONCLUSIONES
Luego de realizar el análisis hidrológico y diseño hidráulico del canal en consideración se puede
concluir lo siguiente:
La hoya hidrográfica estudiada posee un relieve que varía de plano a ligeramente plano con
pendientes menores al 7%, su forma es redonda a oval redonda con coeficientes de
compacidad cercanos a la unidad, lo que sugiere susceptibilidad a crecientes.
El clima de la zona es cálido seco, la temperatura media anual es de 28.8 °C, y oscila entre una
media mínima de 21.6 °C y una media máxima de 36.1 °C. La precipitación media anual en la
es del orden de los 1849 mm y tiene un comportamiento bimodal con los picos más altos en
los meses de mayo y octubre.
Se logró hacer una selección adecuada de las secciones y obras hidráulicas requeridas a lo
largo del canal en estudio, lo anterior se hizo con base a las recomendaciones del titulo D del
reglamento técnico del sector de agua potable y saneamiento básico RAS. Obteniendo los
siguientes resultados para canal y boxculvers
9. RECOMENDACIONES
BIBLIOGRAFÍA
[1] National Aeronautics and Space Administration NASA, «Shuttle Radar Topography Mission
(SRTM),» [En línea]. Available: https://lta.cr.usgs.gov/SRTM1Arc. [Último acceso: 15 8 2016].
[2] G. Monsalve Saenz, Hidrología en la ingeniería, Bogotá: Escuela Colombiana de Ingeniería,
1999.
[3] W. Gámez Morales, Texto Básico de Hidrología, Managua: Universidad Agraria, 2009.
[4] G. A. Barros Cantillo, Curso de Hidrología, Sincelejo: Universidad de Sucre, 1994.
[5] V. T. Chow, D. R. Maidment y L. W. Mays, Hidrología Aplicada, Bogotá: McGRAW - HILL
Interamericana S.A., 1994.
[6] Ministerio de Vivienda, Reglamento Técnico del Sector de Agua Potable y Saneamiento
Básico:Título D. Sistemas de recolección y evacuación de aguas residuales domésticas y aguas
lluvias., Bogotá: Viceministerio de Aguas y Saneamiento Básico., 2012.
[7] Ministerio de Transporte, Instituto Nacional de Vías, Manual de drenaje de carreteras, Bogotá,
2009.
[8] US Army Corps of Engineers, Institute for Water Resources, Hydrologic Engineering Center.,
HEC - RAS, River Analysis System, Davis, California , 2016.
ANEXOS