Proyecto 560457
Proyecto 560457
Proyecto 560457
FACULTAD DE INGENIERÍA
FACULTAD DE INGENIERÍA
______________________________________
______________________________________
______________________________________
______________________________________
Presidente del Jurado
______________________________________
Jurado
______________________________________
Jurado
Primero que todo a Dios por protegerme durante todo el camino, a mis padres y hermana
por guiarme y darme la fortaleza, y a todas las personas que de una u otra forma han contribuido
para la culminación de este objetivo. Este trabajo es fruto de esfuerzo y dedicación.
TABLA DE CONTENIDO
INTRODUCCIÓN ................................................................................................................................ 19
3 METODOLOGÍA........................................................................................................................ 44
6 CARACTERIZACIÓN CLIMATOLÓGICA............................................................................. 88
FIGURA 2-1. LOCALIZACIÓN DEL MUNICIPIO DE SÁCAMA Y HATO COROZAL EN EL DEPARTAMENTO DEL CASANARE 34
FIGURA 4-2. MAPA DE ELEVACIÓN MEDIA DE LA CUENCA EN ESTUDIO. QUEBRADA CASIRBA .............................. 54
FIGURA 5-3. DISTRIBUCIÓN TEMPORAL DE PRECIPITACIÓN TOTAL - ESTACIÓN PTE QUEMADO ............................ 59
FIGURA 5-6. DISTRIBUCIÓN TEMPORAL DE NO. DÍAS CON PRECIPITACIÓN - ESTACIÓN LA CABUYA ...................... 61
FIGURA 5-7. DISTRIBUCIÓN TEMPORAL DE NO. DÍAS CON PRECIPITACIÓN - ESTACIÓN PTE QUEMADO .................. 62
FIGURA 5-8. DISTRIBUCIÓN TEMPORAL DE NO. DÍAS CON PRECIPITACIÓN - ESTACIÓN LA AGUADA...................... 62
FIGURA 5-9. DISTRIBUCIÓN TEMPORAL DE NO. DÍAS CON PRECIPITACIÓN - ESTACIÓN LA AGUADA...................... 63
FIGURA 5-10. POLÍGONOS DE THIESSEN SOBRE LA CUENCA EN ESTUDIO. QUEBRADA CASIRBA ............................ 64
FIGURA 5-11. DISTRIBUCIÓN TEMPORAL DE PRECIPITACIÓN MÁXIMA EN 24 HR- ESTACIÓN LA CABUYA ................ 65
FIGURA 5-12. DISTRIBUCIÓN TEMPORAL DE PRECIPITACIÓN MÁXIMA EN 24 HR- ESTACIÓN PTE QUEMADO ............ 66
FIGURA 5-13. PRECIPITACIONES MÁXIMAS EN 24 HORAS AJUSTADAS A DIFERENTES FUNCIONES DE DISTRIBUCIÓN DE
PROBABILIDAD – ESTACIÓN LA CABUYA ......................................................................................... 70
FIGURA 5-14. PRECIPITACIONES MÁXIMAS EN 24 HORAS AJUSTADAS A DIFERENTES FUNCIONES DE DISTRIBUCIÓN DE
PROBABILIDAD – ESTACIÓN PTE QUEMADO...................................................................................... 70
FIGURA 5-15. LOCALIZACIÓN DE LAS ESTACIONES E ISOYETAS EN LA CUENCA EN ESTUDIO. QUEBRADA CASIRBA ... 72
FIGURA 5-16. DISTRIBUCIÓN DE CAUDALES MÁXIMOS DIARIOS. ESTACIÓN PTE CARRETERA ............................... 74
FIGURA 5-17. DISTRIBUCIÓN ESPACIAL DE LA PRECIPITACIÓN PARA LA CUENCA DE LA ESTACIÓN CERCANA (PTE
CARRETERA) Y DE LA CUENCA EN ESTUDIO (QDA. CASIRBA) ............................................................... 75
FIGURA 5-18. CURVAS IDF. ESTACIÓN LA CABUYA................................................................................... 78
FIGURA 5-19. ANÁLISIS DE CN PARA LA CUENCA EN ESTUDIO. QDA CASIRBA ................................................. 81
FIGURA 7-2. MODELO METODOLOGÍA DE SCS, CON CNII GRUPO B DE 57,4 Y TR DE 2,33 AÑOS. ....................... 103
FIGURA 7-3. MODELO METODOLOGÍA DE BA, CON CNII GRUPO B DE 57,4 Y TR DE 2,33 AÑOS. ......................... 104
FIGURA 7-4. MODELO METODOLOGÍA DE SCS, CON CNII GRUPO C DE 71,2 Y TR DE 2,33 AÑOS. ....................... 105
FIGURA 7-5. MODELO METODOLOGÍA DE BA, CON CNII GRUPO C DE 71,2 Y TR DE 2,33 AÑOS. ......................... 106
FIGURA 7-6. MODELO METODOLOGÍA DE SCS, CON CNII GRUPO B DE 57,4 Y TR DE 50 AÑOS. .......................... 107
FIGURA 7-7. MODELO METODOLOGÍA DE BA, CON CNII GRUPO B DE 57,4 Y TR DE 50 AÑOS. ........................... 108
FIGURA 7-8. MODELO METODOLOGÍA DE SCS, CON CNII GRUPO C DE 71,2 Y TR DE 50 AÑOS. .......................... 109
FIGURA 7-9. MODELO METODOLOGÍA DE BA, CON CNII GRUPO C DE 71,2 Y TR DE 50 AÑOS. ........................... 110
FIGURA 7-10. MODELO METODOLOGÍA DE SCS, CON CNII GRUPO B DE 57,4 Y TR DE 50 AÑOS. ........................ 111
FIGURA 7-11. MODELO METODOLOGÍA DE BA, CON CNII GRUPO B DE 57,4 Y TR DE 50 AÑOS. ......................... 112
FIGURA 7-12. MODELO METODOLOGÍA DE SCS, CON CNII GRUPO C DE 71,2 Y TR DE 100 AÑOS. ...................... 113
FIGURA 7-13. MODELO METODOLOGÍA DE BA, CON CNII GRUPO C DE 71,2 Y TR DE 100 AÑOS. ........................ 114
FIGURA 7-14. MODELO METODOLOGÍA DE SCS, CON CNIII Y TR DE 2.33 AÑOS. ............................................ 116
FIGURA 7-15. MODELO METODOLOGÍA DE SCS, CON CNIII Y TR DE 50 AÑOS. .............................................. 117
FIGURA 7-16. MODELO METODOLOGÍA DE SCS, CON CNIII Y TR DE 100 AÑOS. ............................................. 118
FIGURA 8-1. PLANTA DEL LEVANTAMIENTO TOPOBATIMÉTRICO DE LA ZONA EN ESTUDIO ................................. 119
FIGURA 9-3. VARIACIÓN DE NIVELES DEL RÍO CASANARE. INFORME HIDROLÓGICO DIARIO NO. 134, 2016........... 134
LISTA DE TABLAS
TABLA 2-4. PRECIPITACIÓN ACUMULADA PARA TRES NIVELES DE CONDICIÓN DE HUMEDAD ANTECEDENTE ........... 33
TABLA 4-6. RANGOS APROXIMADOS DE LA PENDIENTE MEDIA DEL CAUCE PRINCIPAL ........................................ 53
TABLA 5-2. DISTRIBUCIÓN TEMPORAL DE PRECIPITACIÓN TOTAL (MM). ESTACIÓN LA CABUYA .......................... 58
TABLA 5-3. DISTRIBUCIÓN TEMPORAL DE PRECIPITACIÓN TOTAL (MM). ESTACIÓN PTE QUEMADO ....................... 59
TABLA 5-4. DISTRIBUCIÓN TEMPORAL DE PRECIPITACIÓN TOTAL (MM). ESTACIÓN LA AGUADA .......................... 60
TABLA 5-5. DISTRIBUCIÓN TEMPORAL DE PRECIPITACIÓN TOTAL (MM). ESTACIÓN LA AGUADA .......................... 60
TABLA 5-6. DISTRIBUCIÓN TEMPORAL DE NO. DÍAS CON PRECIPITACIÓN (DÍAS). ESTACIÓN LA CABUYA ............... 61
TABLA 5-7. DISTRIBUCIÓN TEMPORAL DE NO. DÍAS CON PRECIPITACIÓN (DÍAS). ESTACIÓN PTE QUEMADO ........... 62
TABLA 5-8. DISTRIBUCIÓN TEMPORAL DE NO. DÍAS CON PRECIPITACIÓN (DÍAS). ESTACIÓN LA AGUADA .............. 63
TABLA 5-9. DISTRIBUCIÓN TEMPORAL DE NO. DÍAS CON PRECIPITACIÓN (DÍAS). ESTACIÓN LA AGUADA .............. 63
TABLA 5-11. DISTRIBUCIÓN TEMPORAL DE PRECIPITACIÓN MÁXIMA EN 24 HR (MM). ESTACIÓN LA CABUYA ......... 65
TABLA 5-12. DISTRIBUCIÓN TEMPORAL DE PRECIPITACIÓN MÁXIMA EN 24 HR (MM). ESTACIÓN PTE QUEMADO ...... 66
TABLA 5-15. DISTRIBUCIÓN DE PRECIPITACIÓN MÁXIMA EN 24 HORAS CON DISTRIBUCIÓN GUMBEL .................... 71
TABLA 5-17. ESTACIONES BASE PARA ISOYETAS TOTALES ANUALES (MM) ...................................................... 72
TABLA 5-18. CAUDALES BASE MÁXIMOS DIARIOS PARA LA CUENCA QDA CASIRBA........................................... 75
TABLA 5-20. DATOS DE LLUVIA DE LA ESTACIÓN EN ESTUDIO PARA ELABORAR LAS CURVAS IDF ........................ 77
TABLA 5-22. CLASIFICACIÓN DE TIPOS DE SUELO PRESENTES EN LA CUENCA DE LA QDA CASIRBA ....................... 79
TABLA 5-23. CLASIFICACIÓN DE USO DE SUELO PRESENTES EN LA CUENCA DE LA QDA CASIRBA ......................... 79
TABLA 6-3. DISTRIBUCIÓN TEMPORAL DEL BRILLO SOLAR (HORAS). ESTACIÓN PAZ DE ARIPORO......................... 91
TABLA 6-4. DISTRIBUCIÓN TEMPORAL DE EVAPORACIÓN (MM). ESTACIÓN PAZ DE ARIPORO .............................. 92
Dicho análisis de crecientes y niveles se asocian a los periodos de retorno de 2,33 (que
simula condiciones medias) 50 y 100 años, estimados mediante el uso de metodologías directas
recomendadas por el “Manual de drenaje de carreteras” (Instituto Nacional de vías, 2009) como
el de la SCS Soil Conservation Service (Oficina de Conservación de Suelos de los Estados
Unidos) conocido como el Hidrograma Unitario y de modelos de cómputo como los software
HEC-HMS (desarrollado por el cuerpo de Ingenieros de Estados Unidos) y ArGIS como sistema
de información geográfica SIG.
Este estudio se llevará a cabo con el fin de definir la cota mínima que debe tener la
estructura (puente sobre la Quebrada Casirba) en el sitio de cruce sobre la vía y garantizar la
seguridad del mismo. Para ello, será necesario realizar los análisis de parámetros hidrológicos en
los que se contemplará la estimación de caudales máximos, precipitación, caracterización
morfométrica detallada de la cuenca Quebrada Casirba y el análisis climatológico con los
parámetros de brillo solar, evaporación, humedad relativa y temperatura a nivel mensual
multianual de la zona de estudio.
17
ABSTRACT
The present work is based on the duration of the growth and the maximum levels of the
Quebrada Casirba in the crossing site on the La Cabuya - Sácama road in the department of
Casanare, framed within the line of investigation Sanitation of communities.
This analysis of floods and levels are associated to the return periods of 2.33 (simulating
average conditions) 50 and 100 years, estimated through the use of direct methodologies
recommended by the “Manual de drenaje de carreteras” (Instituto Nacional de vías, 2009) such
as the SCS Soil Conservation Service, known as the Unitary Hydrogram and computer models
such as the HEC-HMS software (developed by the United States Engineers) and ArGIS as a
geographic information system SIG.
This study refers to a rope in order to define the minimum height that the structure must
have (bridge over the Quebrada Casirba) at the crossing point on the road and guarantee the
safety of it. For this, it will be necessary to carry out the analysis of the hydrological
parameters the estimation of maximum flows, precipitation, detailed morphometric
characterization of the Quebrada Casirba basin and the climatological analysis with the
parameters of solar brightness, evaporation, relative humidity and temperature at the multi-year
monthly level of the study area.
18
INTRODUCCIÓN
Para llevar a cabo el estudio será necesario acceder a información secundaria como el
Plan de desarrollo: Unidos construyendo el cambio con un gobierno para todos (Casanare A. M.,
2016-2019) y el Plan vial departamental de Casanare (públicas, 2010-2019), los cuales recopilan
y presentan de manera global las generalidades de la zona de estudio en el departamento de
Casanare donde se localiza el puente objeto de este proyecto.
La estructura del informe estará basada en: Generalidades del sitio de estudio,
antecedentes, marco de referencia, metodologías, análisis hidroclimatológicos, modelación
hidráulica, conclusiones y recomendaciones. Este documento hace parte de los requisitos de la
Universidad Católica de Colombia para otorgar el título de Especialista en Recursos Hídricos.
19
1 GENERALIDADES DEL TRABAJO DE GRADO
La vía La Cabuya – Sácama es considerada como una vía principal del departamento del
Casanare en la cordillera oriental, la cual de acuerdo con el “Manual de diseño geométrico de
vías” (Ospina, 2002) está catalogada como de orden secundario al unir cabeceras municipales
entre sí. La vía conecta los municipios de Sácama, Hato Corozal, Paz de Ariporo, Pore y Yopal.
De acuerdo con los habitantes de la zona y con los medios de comunicación como la
emisora la voz del Cinaruco 105 AM y su portal web (Cinaruco, 2016), el puente localizado
sobre la Quebrada Casirba y el cruce de la vía mencionada, se encuentra en riesgo de colapsar
dado el deterioro que presenta en sus cimientos y accesos. Este puente es transitado por
vehículos pesados como tractomulas con restricción inferior a las 20 toneladas.
Por lo anterior, este estudio permitirá conocer la cota mínima que debe tener el puente a partir de
los análisis hidroclimatológicos donde se estimará el caudal y los niveles máximos de la
20
corriente con ayuda de curvas de calibración de unas secciones topobatimétricas en el sitio de
análisis construidas a partir de la ecuación de Manning para flujo uniforme.
¿Cuál debe ser la cota mínima de elevación del puente para que la hidráulica de flujo de
la Quebrada Casirba no presente alteraciones ante eventos de crecientes para periodos de retorno
2,33 (condiciones medias) 50 y de 100 años?.
1.3 JUSTIFICACIÓN
Además, los resultados que se obtendrán de dicho análisis son básicos y pueden llegar a
ser usados para que otras especialidades puedan definir a futuro la geometría de la estructura de
puente requerido y poder desarrollar un diseño para la proyección a una futura licitación de la
construcción del mismo, dado las inversiones que ha realizado el Ministerio de Hacienda años
anteriores (Hacienda, 2014) y que en la actualidad (2019) no se ha visto plasmada ninguna
intervención.
21
Finalmente, se espera brindar un beneficio a la comunidad cercana al área de estudio que
requiere de la comunicación terrestre entre los municipios de Sácama, Hato Corozal, Paz de
Ariporo, Pore y Yopal.
1.4 OBJETIVOS
Estimar y analizar las crecientes y los niveles máximos en la Quebrada Casirba sobre la
vía La Cabuya – Sácama en el departamento del Casanare, mediante el software HEC-HMS para
periodos de retorno de 2,33 50 y 100 años.
22
2 MARCOS DE REFERENCIA
Es la superficie del terreno que contribuye al escurrimiento, hacia el cauce principal y sus
tributarios, delimitada por la divisoria de aguas o línea imaginaria que se ubica en los puntos de
mayor elevación topográfica.
2.1.1.2 Orientación
23
La orientación de la cuenca con respecto al movimiento del sol y a la interposición de las
cadenas montañosas con respecto a las corrientes de aire, es un factor que influye en su
comportamiento hidro-meteorológico.
Este coeficiente es propuesto por Gravelius, compara la forma de la cuenca con la de una
circunferencia, cuyo círculo inscrito tiene la misma área de la cuenca en estudio. Se define como
la razón entre el perímetro de la cuenca que es la misma longitud del parteaguas que la encierra y
el perímetro de la circunferencia. Este coeficiente adimensional, independiente del área estudiada
tiene por definición un valor de uno para cuencas imaginarias de forma exactamente circular.
Nunca los valores del coeficiente de compacidad serán inferiores a uno. El grado de
aproximación de este índice a la unidad indicará la tendencia a concentrar fuertes volúmenes de
aguas de escurrimiento, siendo más acentuado cuanto más cercano a uno sea, es decir mayor
concentración de agua y se calcula mediante la siguiente expresión:
P
K C 0,282
A0 , 5
Dónde:
Kc: coeficiente de compacidad
P: Perímetro de la cuenca (longitud de la línea parteaguas)
A: área de la cuenca
24
Cuanto más irregular sea la cuenca, mayor será su coeficiente de compacidad. Este
coeficiente define la forma de la cuenca, respecto a la similaridad con formas redondas, dentro
de rangos que se muestran a continuación (ver Tabla 2-1), determinados por la FAO (Estudio
FAO Riego y Drenaje, 1985):
Es el tiempo que toma el agua superficial en llegar desde el punto más alejado de la
cuenca hasta el sitio de interés. Para la cuenca en estudio, el tiempo de concentración se calculará
a partir de las expresiones propuestas en la Tabla 2-2.
25
Fórmula Ecuación Parámetros
Tc = Tiempo de concentración de la cuenca (h)
4 A 1.5L
Giandotti Tc
L = Longitud del cauce principal (km)
25,3 L S S = Pendiente media del cauce principal (m/m)
A = Área de la cuenca (km2)
Tc = Tiempo de concentración de la cuenca (h)
L
0,5
Johnstone y
Tc 2,6 0,5 L = Longitud del cauce principal (km)
S
Cross
S = Pendiente media del cauce principal (m/km)
Tc = Tiempo de concentración de la cuenca (h)
L3
0 , 385
26
Se consideran corrientes de primer orden, aquellas corrientes fuertes, portadoras de aguas
de nacimientos y que no tienen afluentes. Cuando dos corrientes de orden uno se unen, resulta
una corriente de orden dos. De manera general, cuando dos corrientes de orden i se unen, resulta
una corriente de orden i+1. Cuando una corriente se une con otra de orden mayor, resulta una
corriente que conserva el mayor orden.
Este índice relaciona la longitud de la red de drenaje y el área de la cuenca sobre la cual
drenan las corrientes hídricas.
Longitudcorrientes (km)
Dd
Areacuenca (km 2 )
Con el fin de catalogar una cuenca como bien o mal drenada, analizando su densidad de
drenaje, se puede considerar que valores de densidad de drenaje próximo a 0.5 km/km 2 o
mayores indican la eficiencia de la red de drenaje.
Es la distancia medida desde el punto más alejado de la cuenca hasta el sitio donde se
ubicará la obra de drenaje. Para el caso de estudio será la Quebrada Casirba.
27
Tabla 2-3. Categorías de clasificación del índice de sinuosidad
Tipo de canal Índice
Rectilíneo < 1,05
Sinuoso 1,05 - 1,5
Meandrizante > 1,5
Entrelazado > 1,8
Anastomosado >2
2.1.3.1 Elevación
Es uno de los parámetros más determinantes de la oferta hídrica y del movimiento del
agua a lo largo de la cuenca. De ella dependen en gran medida la cobertura vegetal, la biota, el
clima, el tipo y uso del suelo y otras características fisiográficas de un territorio.
La elevación media de una cuenca, se calcula como el promedio ponderado de las áreas
entre curvas de nivel y el área total, este método se conoce como área-elevación. La ecuación
para el cálculo de la elevación media corresponde a la expresión:
n
Ai ei
Em i 1
AT
Donde
Ai: Es el área i-ésima entre curvas de nivel
Ei: Es la elevación media entre curvas de nivel
AT: Es el área total de la cuenca
28
2.1.3.3 Pendiente:
2
S
Li
1
Si
Dónde:
S: Pendiente media del cauce principal (m/m)
L: Longitud del tramo (m) entre las curvas de nivel
Si: Pendiente media del tramo (m/m)
Son cuerpos de agua cerrados que permanecen en un mismo lugar sin correr ni fluir como
ciénagas, lagos, lagunas, humedales o pantanos
Entendiéndose como Sistema lótico todas las masas de agua que se mueven
continuamente en una misma dirección. Ejemplo: ríos, caños y quebradas.
2.1.4.3 Caudal
El caudal de un río, es la cantidad de agua que fluye a través de una sección transversal y
se expresa en volumen por unidad de tiempo. El caudal en un tiempo dado puede medirse por
diferentes métodos y la elección del método depende de las condiciones del sitio de estudio.
29
2.1.5 Variables y parámetros para modelación hidrológica
Para el desarrollo del estudio se realizará la modelación hidrológica a través del software
HEC-HMS desarrollado por Hydrologic Engineering Center (HEC) del United States Corps of
Engineers, el cual ha diseñado los modelos hidráulicos e hidrológicos con mayor reconocimiento
y aceptación a nivel internacional. El HEC-HMS (Hydrologic Model Sistem versión 4.0)
(Engineers, 2010), fue diseñado para simular la escorrentía superficial de respuesta de una
cuenca a la precipitación mediante la representación de la cuenca como un sistema
interconectado de componentes hidrológicos e hidráulicos, cada uno de las cuales refleja un
aspecto del proceso precipitación-escorrentía a partir de parámetros físicos e hidrológicos.
2.1.5.1 Precipitación
30
2.1.5.2 Distribución espacial de la precipitación
Los eventos de lluvia que van a ser simulados deben ser distribuidos tanto espacialmente
como en forma temporal. La distribución temporal es expresada generalmente en forma gráfica o
analítica mediante hietogramas. Para conseguir esta distribución se deben analizar los
pluviogramas de los eventos extremos presentados identificando en el tiempo los incrementos de
lluvia. Se debe agrupar los eventos de acuerdo a la duración, pues las lluvias se distribuyen en
forma diferente según sea su duración.
El grupo hidrológico de suelos describe el tipo de suelo. El uso y el tratamiento del suelo
describen el tipo y la condición de la cubierta vegetal. La condición hidrológica se refiere a la
capacidad de la superficie de la hoya hidrográfica para aumentar o impedir la escorrentía directa.
La condición de humedad antecedente tiene en cuenta la historia reciente de la precipitación y,
consecuentemente, es una medida de la cantidad almacenada por la hoya.
31
Clasificación hidrológica de los suelos
Los grupos hidrológicos en que se pueden dividir los suelos son utilizados en el
planteamiento de cuencas para la estimación de la escorrentía superficial a partir de la
precipitación. Las propiedades de los suelos “desnudos”, luego de un humedecimiento
prolongado, que son consideradas para estimar la tasa mínima de tración, son: profundidad del
nivel freático en época de invierno, infiltración y permeabilidad del suelo luego de un
humedecimiento prolongado, y la profundidad hasta el estrato de permeabilidad muy baja. La
influencia de la cobertura vegetal es tratada independientemente. Los suelos han sido
clasificados en cuatro grupos A, B, C, y D, de acuerdo con el potencial de escurrimiento.
Grupo A. (Bajo potencial de escorrentía). Suelos que tienen alta tasa de infiltración
incluso cuando estén muy húmedos. Consisten en arenas o gravas profundas, bien a
excesivamente drenadas. Estos suelos tienen una alta tasa de transmisión de agua.
Grupo D. (Alto potencial de escorrentía). Suelos con infiltración muy lenta cuando están
muy húmedos. Consisten en suelos arcillosos con alto potencial de expansión; suelos con nivel
freático alto permanente; suelos con estrato arcilloso superficial; suelos con infiltración muy
32
lenta debido a sales o álcalis y suelos poco profundos sobre material casi impermeable. Estos
suelos tienen una tasa de transmisión de agua muy lenta.
El uso del suelo pertenece a la cobertura de la hoya, incluyendo todo tipo de vegetación,
humus vegetal, tierras en descanso (suelo limpio), así como usos no agrícolas, como superficies
de agua (lagos, ciénagas u otras), superficies impermeables (caminos, techos, etc.) y áreas
urbanas. El método del número de curva de escorrentía distingue entre suelos cultivados, prado
y bosques.
El método del número de curva de escorrentía tiene tres niveles de humedad antecedente,
dependiendo de la precipitación total en los cinco días previos a la tormenta que se analiza. La
condición de humedad antecedente seca (AMC I) tienen el menor potencial de escorrentía, con
los suelos estando lo suficientemente secos para un arado satisfactorio o para que una siembra se
lleve a cabo. La condición de humedad antecedente promedio (AMC II) tiene un potencial de
escorrentía promedio. La condición de humedad antecedente húmeda (AMC III) tiene el mayor
potencial de escorrentía, con la hoya hidrográfica prácticamente saturada de precipitaciones
anteriores. La condición de humedad antecedente puede ser determinada a partir de la
información presentada en la Tabla 2-4.
Tabla 2-4. Precipitación acumulada para tres niveles de condición de humedad antecedente
Condición de humedad Precipitación acumulada de los 5 días
Antecedente (AMC) previos al evento en consideración (mm)
I 0 – 36
II 36,1 – 52,5
III Más de 52,5
33
Para determinar el número de curva de escorrentía CN teniendo en cuenta los criterios
ateriormente descritos, dirigirse al Manual de drenajes de carreteras (Instituto Nacional de vías,
2009, págs. 2-85) en las Tabla 2.14, Tabla 2.15, Tabla 2.16 y Tabla 2.17.
Se debe tener en cuenta que de acuerdo con el método del SCS, se deben ajustar las
características de la cuenca a los requerimientos, condiciones y clasificación básica del método,
usando como referencia las definiciones y conceptos del documento "Computer Assisted
Floodplain Hydrologic & Hydraulics" (Hoggan, 1990).
La Quebrada Casirba se encuentra entre los municipios de Sácama y Hato Corrozal, como
se muestra en la Figura 2-1 dos de los 19 municipios en los que se divide el departamento del
Casanare. Sobre este cuerpo de agua se localiza una estructura puente, el cual permite la
comunicación terrestre entre los municipios de Sácama, Hato Corozal, Paz de Ariporo, Pore y
Yopal en el departamento.
Municipio donde se
localiza la Quebrada
en estudio
34
Fuente: (Tierra Colombiana, Municipios de Casanare , 2016)
Por otro lado, la cuenca de la Quebrada Casirba (cuenca en estudio), como se presenta en
la Figura 2-2 se localiza noroccidente del departamento de Casanare entre los municipios de
Sácama y Hato Corrozal y cuenta con un área de 31,87 km2. La Quebrada Casirba tiene una
longitud aproximada de 9,62 km hasta desembocar en la margen derecha del Rio Casanare. La
altura máxima de su límite hidrográfico es de 1500 msnm y su cota mínima de 600 msnm.
35
Figura 2-2. Localización de la cuenca en estudio. Quebrada Casirba
36
Figura 2-3. Estado de la estructura sobre la Quebrada Casirba
2.3.1 Sácama
De acuerdo con los registros del DANE reportados en el Plan de Desarrollo del municipio
de Sácama (Casanare A. M., 2016-2019), este municipio alberga el 0,6% de la población del
departamento de Casanare, cuenta con una población aproximada de 2001 habitantes, de los
37
cuales el 29% equivalente a 587 personas se encuentran en el área rural y el 71% equivalente a
1414 se encuentran en el área urbana.
En cuanto al sistema de salud, en el municipio la IPS está capacitada para hacer atención
de primer nivel con servicios médicos, odontológicos y de laboratorios básicos, hospitalización y
pequeñas intervenciones quirúrgicas. Carece de médico especializado, por lo cual muchos
pacientes deben ser remitidos a Yopal.
Buena parte del empleo se genera a través de la contratación que realiza el municipio,
para la prestación de los servicios establecidos por la ley así como en las obras de infraestructura
realizadas tanto en la cabecera municipal como en la zona rural y los trabajadores de
mantenimiento de carreteras financiados por el INVIAS. El porcentaje de la población
38
beneficiada por la generación de empleo por parte de la contratación se encuentra entre el 20% y
el 30% del total de habitantes del municipio, dependiendo de la actividad.
De acuerdo con el Plan de Desarrollo del Municipio de Hato Corozal (Casanare A. M.,
2012-2015) y los registros del DANE, Este municipio tiene una población a 2012 de 11.541, de
los cuales 6.103 son hombres y 5.408 mujeres equivalente al 53% y 47% de la población
respectivamente. La distribución de población por áreas urbana y rural del municipio de Hato
Corozal es para el año 2012 de 41.54% y 58.46% respectivamente, sin embargo se puede
observar que entre los años 2012 y 2015 hay una tendencia de aumento de la población urbana,
al pasar del 41.54% al 43.47% lo cual equivale a un aumento absoluto de 492 habitantes durante
el periodo de análisis, que al compararlo con la población rural se presenta una disminución en
similares proporciones al pasar 58.46% y 53.43% en el mismo periodo, pese a lo anterior la
mayor concentración poblacional sigue siendo predominantemente rural.
El municipio de Hato Corozal cuenta con una amplia oferta ambiental representada en la
abundancia del recurso hídrico, la flora y la fauna gracias a su localización en la zona de
transición de la Selva Amazónica y La Cordillera Oriental, cuenta con características físico-
bióticas especiales que le permiten el desarrollo de potenciales forestales, pecuarios y energéticas
en su territorio. Predominan principalmente dos modos de producción agrícola, la producción
agrícola de subsistencia y la mecanizada, la primera se desarrolla en la mayoría de fincas de la
39
zona de sabana, utiliza mano de obra no calificada, ausencia total de fertilizantes y asistencia
téCNIIica. En cuanto a la producción en forma comercial esta se desarrolla en la franja de
piedemonte pero en términos generales esta está poco desarrollada y se encuentra representada
por productos como arroz, plátano, maíz, yuca y cacao este último desarrollado a través de
asociaciones de productores con el apoyo del estado pero que realmente no representa unas altas
producciones.
En cuanto al sistema de salud, El municipio cuenta con una oferta de servicios de salud
de primer nivel de atención, con la red pública la ESE Departamental Red Salud Casanare en el
casco urbano y un puesto de salud en el resguardo indígena de Caño Mochuelo (San José del
Ariporo). De igual manera se cuenta con una Red privada que funciona simultáneamente para
40
aumento de la oferta de servicios, con 18 prestadores de servicios de salud; generándose una
dinámica de prestación de servicios con oportunidad en el casco urbano, presentándose
dificultades en su cobertura en zonas rurales apartadas del municipio y en el resguardo indígena.
Por otra parte, las actividades económicas del municipio corresponden a agricultura y
ganadería con restricciones de uso en el suelo correspondientes a agroforestales, agropecuarios,
ganaderos, agrícolas y mineros plasmados en el ACUERDO No. 500.02-004 de 2016.
2.3.4 Pore
En cuanto al sector de salud, el municipio cuenta con un centro de salud que atiende los
servicios de primer nivel (consulta médica general, laboratorio, odontología y urgencias). En
cuanto a su infraestructura, el Hospital de Pore, requiere mantenimiento y dotación
especialmente en la parte eléctrica.
Por otra parte, la condición económica está definida por las actividades agropecuarias
principalmente, sin embargo durante los últimos años se ha posicionado el sector petrolero como
el valor agregado de mayor importancia de los diferentes renglones de la economía municipal.
Los establecimientos comerciales, también son una fuente de empleo de los habitantes como
venta de lubricantes y combustibles, cafeterías, talleres de lavado, restaurantes, droguerías y
estaciones de servicio. En cuanto al sector industrial, este se encuentra poco desarrollado, solo
cuenta con 3 fábricas dedicadas a la producción y comercialización de quesos 1 empresa
empacadora de café, 2 fábricas de prefabricados y 1 taller artesanal y 1 taller de ornamentación.
41
Sin embargo, el municipio cuenta con una localización favorable para la comercialización de
productos, dado que se encuentra sobre la vía marginal de la selva.
2.3.5 Yopal
Aun así, Yopal mantiene unas características de ciudad oferente de trabajo, vivienda y un
alto nivel de asistencialismo en los servicios básicos a la población vulnerable, lo cual hace que
42
absorba una gran cantidad de población con características de formación e ingreso por debajo de
la media nacional, incluso de la media Regional. Además esta atracción también se ejerce sobre
los mismos habitantes del sector rural del municipio y del departamento que han engrosado los
cinturones de población demandante, al no contar con infraestructura básica adecuada cerca de
sus comunidades de origen. Esto nos da a entender que subsiste la inequidad en la inversión
social entre la ciudad y el campo.
43
3 METODOLOGÍA
Esta tesis está dirigida dentro de un enfoque cuantitativo y descriptivo, debido a que se va
a centrar en la obtención de resultados a partir de análisis y modelaciones en la cuenca de estudio
sobre la Quebrada Casirba. Estos datos son susceptibles a estudios estadísticos los cuales
permitirán entregar resultados óptimos y reales de los caudales máximos y niveles del cuerpo de
agua analizado y servir a futuros estudios y disciplinas para conocer el estado de la estructura
que allí se encuentra.
44
3.3 FASES DEL TRABAJO
45
3.3.3 Caracterización morfométrica de la cuenca de estudio
Elevación
Elevación media
Pendiente
Sistemas lénticos
Sistemas lóticos
Variables y parámetros para modelación hidrológica
46
3.3.5 Caracterización Climatológica
47
máxima en 24 horas utilizando las curvas regionalizadas de Intensidad – Duración –
Frecuencia regionalizadas para Colombia (Rodrigo Vargas, 1998).
48
4 CARACTERIZACIÓN MORFOMÉTRICA
49
Tabla 4-1. Caracterización parámetros morfométricos
Área P Coef. Coef. Forma
No. Cuenca
km² km compacidad Oval redonda a oval oblonga
1 Qda. Casirba 31,866665 31,35 1,57 Oval redonda a oval oblonga
50
4.1.2 Tiempo de Concentración
Entendiéndose que es el tiempo que el agua superficial en llegar desde el punto más
alejado de la cuenca hasta el sitio de interés y de acuerdo con los métodos presentados en el ítem
2.1.1.5 Tiempo de concentración. Para la cuenca en estudio (ver Tabla 4-3) el tiempo de
concentración promedio adoptado es de 171 minutos.
Qda. Casirba
Área (km2) 31,87
Longitud del cauce (m) 9616
Longitud del cauce (km) 9,62
H2 (msnm) 1300
H1 (msnm) 600
diferencia de cotas 700
Pendiente Cauce m/m 0,038
Diámetro cuenca (km) 6,37
Kirpich (min) 79,6
Témez (min) 77,8
Williams 176,4
Bransby – Williams 190,5
Giandotti (min) 144,2
Johnstone y Cross 194,2
SCS – Ranser 62,3
Ventura – Heras 76,3
Ven Te Chow 197,7
Cuerpo de Ingenieros 174,2
Promedio 137,3
Desv est 56,6
Prom + desv 194,0
Prom - desv 80,7
Prom (min) 171
Para este estudio se optó por el método de Horton, el cual considera corrientes de primer
orden (Orden 1), aquellas corrientes fuertes, portadoras de aguas de nacimientos y que no tienen
afluentes, como es el caso de la Quebrada Casirba (ver Figura 4-1).
51
4.1.4 Densidad de drenaje (Dd), Índice de Sinuosidad
El cual, teniendo en cuenta la metodología desarrollada por Morisawa (ver Tabla 4-5), la
Quebrada Casirba es de tipo Meándrico y moderadamente drenada.
Como se observa en la Figura 4-1, la Quebrada Casirba, tiene una longitud de 9,6 km
desde su nacimiento hasta su desembocadura en la margen derecha del Río Casanare y
una pendiente de 0,038 m/m, lo cual de acuerdo a la clasificación de pendientes para
cauces (Parámetros geomorfológicos de cuencas hidrográficas, 2016) (ver Tabla 4-6) se
clasifica como una pendiente del cauce Suave.
52
Tabla 4-6. Rangos aproximados de la pendiente media del cauce principal
53
Figura 4-2. Mapa de elevación media de la cuenca en estudio. Quebrada Casirba
54
Figura 4-3. Mapa de pendiente media de la cuenca en estudio. Quebrada Casirba
55
5 CARACTERIZACIÓN HIDROLÓGICA
Sistemas Lénticos
Son cuerpos de agua cerrados que permanecen en un mismo lugar sin correr ni fluir como
ciénagas, lagos, lagunas, humedales o pantanos. Estos sistemas no se presentan cerca de la zona
de estudio.
Sistemas Lóticos
Son todas las masas de agua que se mueven continuamente en una misma dirección.
Ejemplo: ríos, caños y quebradas. La Quebrada Casirba y sus afluentes (sin nombre) clasifican
dentro de este sistema, al igual que el Río Casanare.
Este análisis de precipitación se realizó para las estaciones cercanas a la zona de estudio,
activas y operadas por el IDEAM (IDEAM, Instituto de Hidrología, Meteorología y Estudios
Ambientales, 2014). Las características de las estaciones se presentan en la Tabla 5-1 y su
localización en la Figura 5-1.
56
Tabla 5-1. Características de las estaciones hidrológicas
Código Estación Tipo msnm X Y
36020020 La Cabuya PM 575 1228871,83 1169824,27
36020010 Pte Quemado PM 1020 1208905,62 1165358,61
36010020 La Aguada PG 500 1229638,2 1144606,41
36010030 Sta Rita PG 350 1248488,51 1148469,63
Dónde:
PG: Pluviográfica, PM: Pluviométrica.
57
5.1.1.1 Precipitación total
Estación La Cabuya
58
Figura 5-3. Distribución temporal de precipitación total - Estación Pte Quemado
Tabla 5-3. Distribución temporal de precipitación total (mm). Estación Pte Quemado
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC ANUAL
59,82 73,28 134,40 303,25 460,39 420,62 423,53 374,65 298,95 382,67 236,25 82,56 3250,37
Estación La Aguada
59
Figura 5-4. Distribución temporal de precipitación total - Estación La Aguada
60
5.1.1.2 Precipitación No. de días
Estación La Cabuya
Esta estación registró un promedio de días con precipitación de 15 días para el periodo
1984-2016, lo cual indica que en un año promedio se registran 185 días con precipitación. El mes
de mayor número de días de lluvias es julio y su distribución se presenta en la Figura 5-6.
Figura 5-6. Distribución temporal de No. Días con precipitación - Estación La Cabuya
Tabla 5-6. Distribución temporal de No. Días con precipitación (días). Estación La Cabuya
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC ANUAL
5,03 4,17 9,45 17,13 22,35 23,00 23,14 20,89 17,26 18,66 15,57 8,13 184,78
Esta estación registró un promedio de días con precipitación de 16 días para el periodo
1984-2015, lo cual indica que en un año promedio se registran 191 días con precipitación. El mes
de mayor número de días de lluvias es junio y su distribución se presenta en la Figura 5-7.
61
Figura 5-7. Distribución temporal de No. Días con precipitación - Estación Pte Quemado
Tabla 5-7. Distribución temporal de No. Días con precipitación (días). Estación Pte Quemado
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC ANUAL
7,03 6,53 11,13 17,56 22,28 22,34 22,03 21,83 17,46 19,88 15,86 7,30 191,24
Estación La Aguada
Esta estación registró un promedio de días con precipitación de 13 días para el periodo
1984-2015, lo cual indica que en un año promedio se registran 157 días con precipitación. El mes
de mayor número de días de lluvias es junio y su distribución se presenta en la Figura 5-8.
Figura 5-8. Distribución temporal de No. Días con precipitación - Estación La Aguada
62
Tabla 5-8. Distribución temporal de No. Días con precipitación (días). Estación La Aguada
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC ANUAL
2,06 3,06 8,21 17,05 19,17 20,79 19,61 17,16 16,47 15,22 11,21 7,17 157,17
Esta estación registró un promedio de días con precipitación de 9 días para el periodo
1995-2016, lo cual indica que en un año promedio se registran 118 días con precipitación. El mes
de mayor número de días de lluvias es julio y su distribución se presenta en la Figura 5-9.
Figura 5-9. Distribución temporal de No. Días con precipitación - Estación La Aguada
Tabla 5-9. Distribución temporal de No. Días con precipitación (días). Estación La Aguada
ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC ANUAL
1,05 3,44 5,63 12,45 15,16 15,63 15,79 13,11 11,47 12,21 8,84 3,68 118,47
63
Thiessen (ver Figura 5-10) que tiene en cuenta todos los registros de cada una de las estaciones
y el factor de representatividad de cada una de las estaciones sobre las cuencas en estudio.
De acuerdo a lo anterior, las estaciones que más influencia tienen sobre la cuenca en
estudio son la estación La Cabuya y estación Pte Quemado. En la Tabla 5-10, se presenta el área
total y el porcentaje de área que representa cada estación sobre la cuenca Quebrada Casirba.
64
Tabla 5-10. Áreas sobre la cuenca de estudio. Qda Casirba
Área de la cuenca
Nombre Área sobre la cuenca km2
%
31,87
1,00
Cuenca Qda. Casirba Área total cuenca
Estación La Cabuya 30,98 0,97
Estación Pte Quemado 0,89 0,03
Estación La Cabuya
65
Estación Pte Quemado
Posteriormente, se realizó para cada estación de influencia una prueba estadística de datos
dudosos método del Water Resources Council de los Estados Unidos, recomendado para detectar
los puntos que se alejan de la tendencia de los valores máximos por encima o por debajo y
facilitar la toma de decisiones de retención o eliminación de datos que pueden afectar
significativamente la magnitud de los parámetros estadísticos. Esta metodología determina los
66
valores de los límites superior e inferior a partir de los cuales se eliminan los valores que no se
consideran estadísticamente aceptables. Para estas estaciones, todos los datos se consideran
aceptables, excepto para los valores de los años 2002 y 2003 de la estación Pte Quemado, por lo
cual se retiran de la serie a analizar al no representar al comportamiento de las precipitaciones en
el sitio de cada estación.
67
Tabla 5-13. Prueba de datos dudosos. Estación La Cabuya
Número Datos 32 Desviación Estandar Sy 0,0836
Kn 2,5910 Coeficiente Asimetría Cs 0,1750
Valor Máximo 139,50
Valor Mínimo 79,20
yH 2,2309 yL 1,7979
QH 170 QL 63
n ( yi y )3
Fórmulas: n
( yi y) 2
1 n
Sy Cs i 1
n 1 i 1 ( n 1)( n 2 ) S y
3
yH y Kn S y yL y Kn S y
Q H 10 y H QL 10 y L
68
Tabla 5-14. Prueba de datos dudosos. Estación Pte Quemado
Número Datos 32 Desviación Estandar S y 0,3032
Kn 2,5910 Coeficiente Asimetría Cs -3,1856
Valor Máximo 138,40
Valor Mínimo 5,80
yH 2,6513 yL 1,0800
QH 448 QL 12
n ( yi y )3
Fórmulas: n
( yi y) 2
1 n
Sy Cs i 1
n 1 i 1 ( n 1)( n 2 ) S y
3
yH y Kn S y yL y KnS y
Q H 10 y H QL 10 y L
69
Figura 5-13. Precipitaciones máximas en 24 horas ajustadas a diferentes funciones de
distribución de probabilidad – Estación La Cabuya
70
Utilizando como criterios el resultado gráfico de los ajustes estadísticos presentados
anteriormente y tomando el mejor valor de Chi Cuadrado, la distribución de probabilidad que
mejor se ajusta es Gumbel para las dos estaciones y su distribución para los periodos de retorno
2,33 – 50 y 10 años se presenta en la Tabla 5-15.
71
En la Tabla 5-17, se presentan las coordenadas y los valores de precipitación anual de las
estaciones empleadas y su localización en la Figura 5-15.
72
Como se puede observar en la Figura 5-15, la precipitación media anual disminuye en
sentido oeste este variando desde los 3100 mm hasta los 1800 mm para la zona de estudio. Se
observa que las precipitaciones más altas se presentan en las estaciones Pte Quemado. La
estación más cercana a la zona de estudio es La Cabuya y presenta una precipitación total
multianual de 2870 mm.
Debido a que sobre la Quebrada Casirba o cerca a ésta no hay estaciones de caudales que
permitan determinarlos de manera directa, se utiliza una metodología indirecta o método regional
de Transposición que utiliza relaciones precipitación - Área – Caudal, mediante la siguiente
expresión:
Q estación Qx
A estación * Pestación A x * Px
Dónde:
Qestación: Caudal medio para cada período de retorno en la estación cercana
Q x : Caudal medio para cada Tr (años) en el sitio de proyecto (cuenca Qda. Casirba)
Se tiene que la estación más cercana que registra parámetros de caudal es la estación
Limnimétrica 36017030 Pte Carretera operada por el IDEAM (IDEAM, Instituto de Hidrología,
Meteorología y Estudios Ambientales, 2014), la cual tiene un periodo de registro desde 1997
hasta 2015. Presenta un comportamiento monomodal donde los meses con caudal por encima
del promedio 9,73 m³/s son de mayo a octubre. Posteriormente, se seleccionaron los caudales
73
máximos diarios para cada año y se realizó un ajuste de distribución para estimar los caudales
máximos diarios esperados para los periodos de retorno 2,33 50 y 100 años para diferentes
distribuciones de probabilidad, entre ellas Normal, Gumbel, Pearson, Log-Normal, Log Pearson
y EV3. Utilizando como criterio el resultado grafico de los ajustes estadísticos se seleccionaron
los valores obtenidos mediante el ajuste de Gumbel.
A partir del análisis de distribución espacial de la precipitación (Isoyetas) y con ayuda del
software ArcGIS (ver Figura 5-17) se calculó la Precipitación media para la cuenca de la
estación cercana (Pte Carretera) y de la cuenca en estudio (Qda. Casirba).
74
Figura 5-17. Distribución espacial de la precipitación para la cuenca de la estación
cercana (Pte Carretera) y de la cuenca en estudio (Qda. Casirba)
Tabla 5-18. Caudales base máximos diarios para la cuenca Qda Casirba
75
5.1.4 Cálculo de Curvas IDF (Intensidad de la Precipitación).
Las curvas intensidad duración frecuencia (curvas IDF) constituyen una forma de
sintetizar la información de precipitación para las tormentas de corta duración, y en la mayoría
de los casos, es el insumo básico para la estimación de tormentas o eventos de diseño en la
modelación hidrológica.
Para el desarrollo de este análisis se utilizó la metodología de Vargas (Vargas & Díaz
Granados, 1998), formulada a partir del análisis de más de 200 curvas IDF y registro
pluviométrico concurrente en Colombia, la cual propone la generación sintética de curvas IDF
basada en la regionalización de los coeficientes para la siguiente ecuación:
76
zona, dado que es la que mejor representa los registros de lluvia de la cuenca en estudio. La
región Orinoquía sobrevalora los registros.
Tabla 5-20. Datos de lluvia de la Estación en estudio para elaborar las curvas IDF
Número de
Precipitación P Máxima
días de
Estación Periodo total anual en 24h (M Región
precipitación
(PT en mm) en mm)
(N)
La Cabuya 1986-2016 2870,1 105,3 184,8 Andina
De acuerdo con la metodología presentada, en la Figura 5-18, se presenta las curvas IDF
obtenidas para la estación La Cabuya y el resumen de datos en la Tabla 5-21.
77
Figura 5-18. Curvas IDF. Estación La Cabuya
Para realizar este análisis, se usó el mapa de cobertura y tierra (2010-2014) del IDEAM
(IDEAM, Instituto de Hidrología, Meteorología y Estudios Ambientales, 2014), y el mapa de
suelos de territorio colombiano a escala 1:100.000 para el departamento de Casanare (IGAC,
2019). (Ver Figura 5-19).
78
Tabla 5-22. Clasificación de tipos de suelo presentes en la cuenca de la Qda Casirba
Área
Litología Características
(km²)
Materiales aluviales
Bien drenados, profundos, de texturas moderadamente gruesas y
heterométricos con 0.114174
medianamente finas, muy fuertemente ácidos y fertilidad muy baja
matriz gruesa
Bien drenados, muy superficiales a moderadamente profundos, de texturas
Arcillolitas, lodolitas y
moderadamente gruesas a finas, muy fuerte a extremadamente ácidos y 0.057093
areniscas
fertilidad baja y muy baja
Bien drenados, muy superficiales a moderadamente profundos, de texturas
Arcillolitas, lodolitas y
moderadamente gruesas a finas, muy fuerte a extremadamente ácidos y 0.918872
areniscas
fertilidad baja y muy baja
Sedimentos aluviales Excesivamente drenados, muy superficiales, texturas gruesas, alta saturación
0.422683
heterométricos de aluminio, muy fuertemente ácidos, fertilidad muy baja
Sedimentos aluviales Bien drenados, profundos y superficiales, texturas medianamente finas,
0.015265
finos fuertemente ácidos, fertilidad moderada y baja
Bien drenados, muy superficiales a moderadamente profundos, de texturas
Arcillolitas, lodolitas y
moderadamente gruesas a finas, muy fuerte a extremadamente ácidos y 1.191338
areniscas
fertilidad baja y muy baja
Bien drenados, muy superficiales a moderadamente profundos, de texturas
Arcillolitas, lodolitas y
moderadamente gruesas a finas, muy fuerte a extremadamente ácidos y 7.610572
areniscas
fertilidad baja y muy baja
Arcillolitas No suelos, en pendientes fuertes y erosión severa y muy severa 4.476455
Bien drenados, muy superficiales a moderadamente profundos, de texturas
Arcillolitas, lodolitas y
moderadamente gruesas a finas, muy fuerte a extremadamente ácidos y 16.060844
areniscas
fertilidad baja y muy baja
En este ítem se incluye todo tipo de vegetación de la superficie en estudio. Los resultados
se presentan en la Tabla 5-23, encontrando que en su mayor porcentaje de área se presenta
pastos y bosques.
79
Tipo de Suelo Área (km²)
Pastos enmalezados 0.250854
Bosque de galería y ripario 0.286889
Bosque fragmentado 0.157787
Mosaico de pastos con espacios naturales 1.480414
Bosque fragmentado 0.570865
Bosque denso 26.259472
Teniendo en cuenta que se considera que la cuenca en estudio pertenece al grupo de suelo
B y C, se realiza un promedio aritmético de los valores teniendo en cuenta la clasificación del
uso del suelo y su porcentaje de área dentro de la cuenca y los valores establecidos en el Manual
de drenaje de carreteras (Instituto Nacional de vías, 2009). En la Tabla 5-24, se presenta el CNII
promedio adoptado para la cuenca en estudio.
El CNIII, es considerado para condiciones extremas es decir que tiene el mayor potencial
de escorrentía, de acuerdo con la metodología se adopta cuando la precipitación es mayor a 52
mm. Dado que las precipitaciones para los diferentes periodos de retorno son superiores a este
valor se adopta un valor de CNIII de 86.
80
Figura 5-19. Análisis de CN para la cuenca en estudio. Qda Casirba
81
directamente proporcionales a la cantidad de escorrentía directa, y v) el hidrograma resultante de
un exceso de lluvia dado refleja las características no cambiantes de una cuenca.
d d
tp t l 0 ,6 t c
2 2
tl
Donde es el tiempo de rezago de la cuenca que corresponde a una fracción (60%) del
tc
tiempo de concentración , definido en el capítulo de morfología, y d es la duración de exceso
de precipitación unitaria. Para finalizar la descripción geométrica del hidrograma, el ancho de la
base se define como:
8
T bt tp
3
82
5.1.6.1 Tormenta de diseño
83
5.1.6.2 Distribución espacial de precipitación y reducción por área
Debido a que los registros de precipitación máxima en las cuencas son puntuales por la
localización de las estaciones, son máximos en 24 horas y no corresponden a duraciones
inferiores de tormentas, se determinó la precipitación máxima promedio sobre las cuencas
mediante el uso de diferentes metodologías.
Para determinar el factor de reducción por área, el manual de drenajes del INVIAS
recomienda el uso de tres ecuaciones y el uso del estudio realizado por INGETEC en 2005 para
la Empresa de Acueducto de Bogotá denominado Estudio y revisión de las curvas IDF y del
Análisis Espacial de las Tormentas, curvas PADF para la Sabana de Bogotá – Referencia 2.14
del manual de drenajes.
Las ecuaciones recomendadas por el INVIAS para determinar el factor de reducción por
área son las siguientes:
Ecuación de Fhrüling:
84
Dónde:
A: Área de drenaje de la Cuenca en km2.
Tc: tiempo de concentración de la Cuenca en horas (h)
En cuanto al uso del estudio realizado por INGETEC en 2005 para el Análisis Espacial de
las Tormentas, se presentan a continuación las curvas PADF recomendadas por el INVÍAS, las
cuales presentan las curvas de reducción de la lluvia puntual por espacialidad en la cuenca.
85
duración. Para el caso de la cuenca en estudio, la precipitación correspondiente a un área de
31,86 km2 y una duración igual al tiempo de concentración de aproximadamente 3 horas es de 58
mm y la precipitación para un área de 1 km2 para la misma duración es de 75 mm, por tanto el
factor de reducción por área será el resultado de dividir 58 mm / 75 mm = 0,77.
A continuación se presenta el resumen de los resultados obtenidos a partir del uso de las
ecuaciones mencionadas anteriormente y del uso de las curvas PADF.
Características Ecuación
Institution of Factor de
Cuenca Datos Curvas
Área (km ) 2
Tc (h) Fhruling Civil reducción
Norteamérica PADF
Engineers
Qda
31,86 2,85 0,59 0,92 0,82 0,77 0,78
Casirba
Los resultados muestran que con la ecuación del Institution of Civil Engineers de
Inglaterra se presentan los mayores factores de reducción y con la ecuación de Fhrüling los
menores. Los valores obtenidos con las curvas PADF son intermedios entre los obtenidos con la
ecuación de los datos de Norteamérica y con la ecuación de Fhrüling.
Este método es una forma simple para desarrollar en hietograma de diseño utilizando una
curva de intensidad-duración-frecuencia. El hietograma de diseño producido por este método
especifica la profundidad de precipitación que ocurre en intervalos de tiempo sucesivos de
duración sobre una duración total de tiempo.
Para dar una secuencia más crítica del aguacero, como es usual en este tipo de análisis
por medio del método de Bloques Alternos, los incrementos de lluvia de los hietogramas
mencionados previamente se arreglan de la siguiente manera: el valor más bajo se coloca en el
86
primer lugar, el segundo valor en orden creciente se coloca en último lugar, el tercer valor en tal
orden se ubica en segundo lugar, el cuarto valor en el penúltimo lugar, y así sucesivamente.
87
6 CARACTERIZACIÓN CLIMATOLÓGICA
El clima está definido por las condiciones meteorológicas registradas durante largos
períodos de tiempo y se describe a partir de variables atmosféricas como la temperatura, el brillo
solar, la humedad relativa, el viento y otras más.
Dónde:
CP: Climatológica Principal, CO: Climatológica ordinaria
88
Figura 6-1. Localización de las estaciones climatológicas de la cuenca en estudio.
Quebrada Casirba
89
Tabla 6-2. Características de las estaciones climatológicas
ESTACIONES IDEAM PERIODO
PARÁMETROS
36015010_Paz Ariporo 35235020_Tamara COMÚN
Brillo Solar Hr 1996-2017 - -
Evaporación mms 1996-2016 1996-2016 -
Humedad Relativa % 1996-2017 1996-2017 1996-2017
Nubosidad Octas 1996-2016 1997-2016 1997-2016
Punto de Rocío °C 1996-2017 1996-2017 1996-2017
Temperatura media °C 1996-2017 1996-2017 1996-2017
Temperatura máxima °C 2006-2017 2006-2017 2006-2017
Temperatura mínima °C 1998-2017 2006-2017 2006-2017
Tensión Vapor Mb 1996-2017 1996-2017 1996-2017
Velocidad del viento m/s 1996-2012 - -
A continuación se realiza el análisis de cada uno de los parámetros registrados por las
estaciones Paz de Ariporo y Tamara.
90
Figura 6-2. Distribución temporal del brillo solar – Estación Paz de Ariporo
Tabla 6-3. Distribución temporal del brillo solar (horas). Estación Paz de Ariporo
ESTACION ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC PROM
PAZ DE
254.03 200.16 160.70 132.32 146.30 142.18 151.68 168.98 185.84 197.87 205.23 239.12 182.03
ARIPORO
6.2.2 Evaporación
91
Figura 6-3. Distribución temporal de evaporación – Estación Paz de Ariporo
Para este parámetro se analizan las dos estaciones (Paz de Ariporo y Tamara) como se
presenta en la Figura 6-4 para el periodo común 1996-2017. En la Tabla 6-5, se presentan los
valores asociados a cada estación.
92
humedad relativa con un valor de 85,33 % y febrero el mes de menor humedad relativa con un
valor de 72,53 %.
PAZ DE ARIPORO 62.62 60.05 61.91 73.73 79.86 82.82 82.45 80.50 78.41 77.82 77.50 70.95 74.05
TAMARA 73.94 72.53 74.12 80.06 83.19 84.43 85.33 84.67 82.43 81.29 79.86 77.86 79.97
6.2.4 Nubosidad
Para este parámetro se analizan las dos estaciones (Paz de Ariporo y Tamara) como se
presenta en la Figura 6-5 para el periodo común 1997-2016. En la Tabla 6-6, se presentan los
valores asociados a cada estación.
93
N es igual 0 el cielo está despejado.
N es mayor o igual que 1 y menor o igual que 3 el cielo es poco nuboso.
N es mayor o igual que 4 y menor o igual que 5 el cielo es nuboso.
N es mayor o igual que 6 y menor o igual que 7 el cielo es muy nuboso.
N es igual a 8 el cielo es cubierto.
La estación Paz de Ariporo, para el periodo en estudio abril es el mes con mayor
nubosidad con un valor de 5,7 octas y enero el de menor precipitación con un valor de 3,5 octas,
siendo su promedio anual 5 octas y de acuerdo a los rangos, se clasifica como una zona nubosa.
La estación Tamara, para el periodo en estudio junio y julio son los meses con mayor
nubosidad con un valor de 5,9 octas y enero el de menor precipitación con un valor de 3,2 octas,
siendo su promedio anual 5 octas y de acuerdo a los rangos, se clasifica como una zona nubosa.
94
6.2.5 Punto de Rocío
Para este parámetro se analizan las dos estaciones (Paz de Ariporo y Tamara) como se
presenta en la Figura 6-7 para el periodo común 1996-2017. En la Tabla 6-8, se presentan los
valores asociados a cada estación.
95
La temperatura media en la estación Paz de Ariporo es de 27,09 °C. Los periodos
diciembre a abril presentan valores por encima del promedio, siendo marzo el mes con mayor
temperatura registrada con un valor de 28,82 °C y julio el mes con menor temperatura con un
valor de 25,57 °C.
96
6.2.7 Temperatura Máxima
Para este parámetro se analizan las dos estaciones (Paz de Ariporo y Tamara) como se
presenta en la Figura 6-8 para el periodo común 2006-2017. En la Tabla 6-9, se presentan los
valores asociados a cada estación.
La temperatura máxima en la estación Paz de Ariporo es de 28,91 °C. Los periodos enero
a abril presentan valores por encima del promedio, siendo marzo el mes con mayor temperatura
máxima registrada con un valor de 35,98 °C y junio el mes con menor temperatura máxima con
un valor de 32,45°C.
La temperatura máxima en la estación Tamara es de 22,32 °C. Los periodos enero a abril
presentan valores por encima del promedio, siendo febrero el mes con mayor temperatura
máxima registrada con un valor de 30,72°C y julio el mes con menor temperatura máxima con
un valor de 27,57°C.
97
Tabla 6-9. Distribución mensual de la Temperatura Máxima (°C).
ESTACION ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC PROM
PAZ DE ARIPORO 34.53 35.53 35.98 35.52 33.63 32.45 32.62 33.12 33.42 33.58 33.54 33.80 33.98
TAMARA 29.81 30.72 30.23 29.97 28.55 27.90 27.57 28.09 28.92 28.50 28.31 28.35 28.91
Para este parámetro se analizan las dos estaciones (Paz de Ariporo y Tamara) como se
presenta en la Figura 6-9 para el periodo común 2006-2017. En la Tabla 6-10, se presentan los
valores asociados a cada estación.
La temperatura mínima en la estación Paz de Ariporo es de 28,91 °C. Los periodos enero
a mayo presentan valores por encima del promedio, siendo febrero el mes con mayor
temperatura mínima registrada con un valor de 20,92 °C y septiembre el mes con menor
temperatura mínima con un valor de 20,22°C.
La temperatura mínima en la estación Tamara es de 22,32 °C. Los periodos enero a abril
y octubre presentan valores por encima del promedio, siendo febrero el mes con mayor
temperatura mínima registrada con un valor de 17,94°C y agosto el mes con menor temperatura
mínima con un valor de 16,92°C.
98
Tabla 6-10. Distribución mensual de la Temperatura Mínima (°C).
ESTACION ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC PROM
PAZ DE ARIPORO 20.60 20.92 21.73 21.52 21.00 20.67 20.38 20.23 20.22 20.63 20.76 20.82 20.79
TAMARA 16.86 17.94 18.36 17.82 17.47 16.95 16.88 16.92 16.69 17.10 17.53 17.07 17.30
99
6.2.10 Velocidad del viento
Para este parámetro se utilizó solo la información suministrada por la estación Paz de
Ariporo, dado que es la estación más cercana al área de estudio con registros de velocidad.
La estación Paz de Ariporo cuenta con registros en el periodo 1996-2012 reportando una
velocidad promedio de viento de 3,36 m/s, como se muestra en la Figura 6-11. Los vientos más
fuertes se presentan en enero, febrero, marzo y diciembre, siendo febrero el mes de mayor
registro con un valor de 4.70 m/s y el más bajo en julio con 2,76 m/s.
100
7 MODELACIÓN HIDROLÓGICA
Esta modelación se realizó a partir del método del Hidrograma Unitario de la Oficina de
Conservación de Suelos de los Estados Unidos de América (Ven Te Chow D. M., 1994) y del
método de bloques alternos empleando el paquete de cómputo HEC-HMS o Sistema de
Modelación Hidrológica desarrollado por el Centro de Ingeniería Hidrológica del Cuerpo de
Ingenieros de los Estados Unidos junto con la herramienta de información geográfica ArcGis
(ver Figura 7-1) y los parámetros descritos en el ítem 5. CARACTERIZACIÓN
HIDROLÓGICA. Para esto, se presenta en la
CNIII (C) 86
Q (m³/s)
5,31
Tr 2,33 (años)
Q (m³/s)
13,71
Tr 50 (años)
Q (m³/s)
15,47
Tr 100 (años)
101
Figura 7-1. Modelo en HEC-HMS de la cuenca de la Qda Casirba
Este modelo se presenta para dos escenarios teniendo en cuenta los hietogramas de
precipitación para la metodología de SCS y la de bloques alternos, con un CNII para el grupo B
de 57,4 y un periodo de retorno de 2,33 años. El resumen de los resultados se presentan en la
Tabla 7-2.
102
Figura 7-2. Modelo Metodología de SCS, con CNII grupo B de 57,4 y Tr de 2,33 años.
7.1.2 Metodología de bloques alternos (BA), con CNII grupo B de 57,4 y Tr de 2,33 años.
Los resultados de la modelación teniendo en cuenta estos parámetros se presentan a
continuación:
103
Figura 7-3. Modelo Metodología de BA, con CNII grupo B de 57,4 y Tr de 2,33 años.
Este modelo se presenta para dos escenarios teniendo en cuenta los hietogramas de
precipitación para la metodología de SCS y la de bloques alternos, con un CNII para el grupo C
de 71,2 y un periodo de retorno de 2,33 años. El resumen de los resultados se presentan en la
Tabla 7-2.
104
7.2.1 Metodología de SCS, con CNII grupo C de 71,2 y Tr de 2,33 años.
Los resultados de la modelación teniendo en cuenta estos parámetros se presentan a
continuación:
Figura 7-4. Modelo Metodología de SCS, con CNII grupo C de 71,2 y Tr de 2,33 años.
7.2.2 Metodología de bloques alternos (BA), con CNII grupo C de 71,2 y Tr de 2,33 años.
Los resultados de la modelación teniendo en cuenta estos parámetros se presentan a
continuación:
105
Figura 7-5. Modelo Metodología de BA, con CNII grupo C de 71,2 y Tr de 2,33 años.
Este modelo se presenta para dos escenarios teniendo en cuenta los hietogramas de
precipitación para la metodología de SCS y la de bloques alternos, con un CNII para el grupo B
de 57,4 y un periodo de retorno de 2,33 años. El resumen de los resultados se presentan en la
Tabla 7-2.
106
Figura 7-6. Modelo Metodología de SCS, con CNII grupo B de 57,4 y Tr de 50 años.
7.3.2 Metodología de bloques alternos (BA), con CNII grupo B de 57,4 y Tr de 50 años.
107
Figura 7-7. Modelo Metodología de BA, con CNII grupo B de 57,4 y Tr de 50 años.
Este modelo se presenta para dos escenarios teniendo en cuenta los hietogramas de
precipitación para la metodología de SCS y la de bloques alternos, con un CNII para el grupo B
de 57,4 y un periodo de retorno de 2,33 años. El resumen de los resultados se presentan en la
Tabla 7-2.
108
7.4.1 Metodología de SCS, con CNII grupo C de 71,2 y Tr de 50 años.
Los resultados de la modelación teniendo en cuenta estos parámetros se presentan a
continuación:
Figura 7-8. Modelo Metodología de SCS, con CNII grupo C de 71,2 y Tr de 50 años.
7.4.2 Metodología de bloques alternos (BA), con CNII grupo C de 71,2 y Tr de 50 años.
Los resultados de la modelación teniendo en cuenta estos parámetros se presentan a
continuación:
109
Figura 7-9. Modelo Metodología de BA, con CNII grupo C de 71,2 y Tr de 50 años.
Este modelo se presenta para dos escenarios teniendo en cuenta los hietogramas de
precipitación para la metodología de SCS y la de bloques alternos, con un CNII para el grupo B
de 57,4 y un periodo de retorno de 2,33 años. El resumen de los resultados se presentan en la
Tabla 7-2.
110
7.5.1 Metodología de SCS, con CNII grupo B de 57,4 y Tr de 100 años.
Los resultados de la modelación teniendo en cuenta estos parámetros se presentan a
continuación:
Figura 7-10. Modelo Metodología de SCS, con CNII grupo B de 57,4 y Tr de 50 años.
7.5.2 Metodología de bloques alternos (BA), con CNII grupo B de 57,4 y Tr de 100 años.
Los resultados de la modelación teniendo en cuenta estos parámetros se presentan a
continuación:
111
Figura 7-11. Modelo Metodología de BA, con CNII grupo B de 57,4 y Tr de 50 años.
Este modelo se presenta para dos escenarios teniendo en cuenta los hietogramas de
precipitación para la metodología de SCS y la de bloques alternos, con un CNII para el grupo B
de 57,4 y un periodo de retorno de 2,33 años. El resumen de los resultados se presentan en la
Tabla 7-2.
112
7.6.1 Metodología de SCS, con CNII grupo C de 71,2 y Tr de 100 años.
Los resultados de la modelación teniendo en cuenta estos parámetros se presentan a
continuación:
Figura 7-12. Modelo Metodología de SCS, con CNII grupo C de 71,2 y Tr de 100 años.
7.6.2 Metodología de bloques alternos (BA), con CNII grupo C de 71,2 y Tr de 100 años.
Los resultados de la modelación teniendo en cuenta estos parámetros se presentan a
continuación:
113
Figura 7-13. Modelo Metodología de BA, con CNII grupo C de 71,2 y Tr de 100 años.
114
Caudal
Modelo Tr (años) CNII Metodología máximo Volumen (m³)
(m³/s)
SCS 7,1
Modelo 1 2,33 B 246,4
Bloques Alternos 7,2
SCS 23,8
Modelo 2 2,33 C 414,1
Bloques Alternos 23,6
SCS 34,1
Modelo 3 50 B 794,8
Bloques Alternos 33,9
SCS 75,4
Modelo 4 50 C 1204,9
Bloques Alternos 73,9
SCS 41,8
Modelo 5 100 B 930,7
Bloques Alternos 41,5
SCS 88,2
Modelo 6 100 C 1390,2
Bloques Alternos 86,5
Los análisis hidráulicos del presente estudio se realizan para condiciones extremas, es
decir, teniendo un CNIII de 86 para los periodos de retorno de 2.33, 50 y 100 años, dado que
los caudales son más representativos con la zona en estudio. El resumen de resultados se
presenta a continuación en la Tabla 7-3.
115
7.7.1 Metodología de SCS, con CNIII y Tr de 2.33 años.
Los resultados de la modelación teniendo en cuenta estos parámetros se presentan a
continuación:
116
Figura 7-15. Modelo Metodología de SCS, con CNIII y Tr de 50 años.
117
Figura 7-16. Modelo Metodología de SCS, con CNIII y Tr de 100 años.
Dado que estos son los caudales más representativos para condiciones extremas, serán
los asumidos para desarrollar la modelación hidráulica que se presentará a continuación.
118
8 MODELACIÓN HIDRÁULICA
119
8.1 DETERMINACIÓN DEL COEFICIENTE DE RUGOSIDAD DE MANNING
120
Se resalta zonas con vegetación en la planicie de inundación, la cual se puede considerar
inefectiva para el flujo.
8.2 GÁLIBOS
Sobre corrientes de agua, relativamente limpias en toda época: mínimo dos metros
(2.00 m) por encima del Nivel de Aguas Máximas Extraordinarias (N.A.M.E.).
Sobre corrientes de agua que en algunos períodos transportan desechos, troncos y
otros objetos voluminosos: mínimo dos metros con cincuenta centímetros (2.50
m) por encima del N.A.M.E.
Sobre carreteras. Mínimo cinco metros (5.00 m), salvo que se estipule un valor
mayor para el caso particular.
Sobre vías férreas. Mínimo cinco metros con cincuenta centímetros (5.50 m),
salvo que se estipule un valor mayor para el caso particular.
Sobre cursos hídricos navegables. Dependerá del calado máximo de navegación,
por lo que el valor debe ser definido por el Ministerio de Transporte.”
121
Régimen del flujo: El programa HEC-RAS permite calcular de manera simultánea
flujo de tipo Mixto, es decir, considerándolo subcrítico o supercrítico.
Altura en la sección transversal inicial y final considerada es de: 0.078 m/m.
Caudal: correspondiente a los presentados en la Tabla 7-3 (condiciones
extremas).
Coeficientes de pérdidas de energía, “n” de Manning: correspondientes a 0.09 en
las márgenes y 0.05 en el cauce.
122
Figura 8-3. Sección hidráulica de la Quebrada Casirba (Tr=100años)
El puente sobre la Quebrada Casirba se encuentra a una altura de 5,40m desde el fondo
del cauce de acuerdo a las batimetrías obtenidas, es decir que el gálibo del puente teniendo en
cuenta un periodo de retorno de 100 años se encuentra a 2,70 m siendo mayor al mínimo
recomendado de 2.5 m por Manual de Diseño Geométrico de Carreteras (Instituto Nacional de
vías, 2009).
123
8.4 ANÁLISIS DE POSIBLE CONTROL HIDRÁULICO
Se debe tener en cuenta que la Quebrada Casirba es afluente al Río Casanare, por lo que
se debe considerar el río como un posible control hidráulico. Para esto, con un alto nivel de
incertidumbre y con respecto a información secundaria y estaciones limnigráficas como
(36027010) Puente Quemado y (36027020) La Cabuya operadas por el IDEAM (IDEAM,
Instituto de Hidrología, Meteorología y Estudios Ambientales, 2014) se estima que el Río
Casanare presenta un comportamiento en niveles aproximadamente de 6 m (Ver ANEXO 7
Modelación Hidráulica).
124
Figura 8-4. Perfil y sección hidráulica de la Quebrada Casirba. Teniendo en cuenta el Río
Casanare con nivel de 6 m como un posible control hidráulico
En vista que el nivel del Río Casanare es de 6m, la Quebrada Casirba busca empatarse
por energía, lo cual hace que se genere un área de inundación hasta el puente.
El resumen de resultados para los periodos de retorno de 2,33, 50 y 100 años teniendo
en cuenta el posible control hidráulico se presenta a continuación en la Tabla 8-2. Para
125
evidenciar los cambios de nivel, se muestran cinco registros aguas arriba y cinco registros
aguas abajo del puente en estudio.
Tabla 8-3. Resumen resultados modelos HEC-RAS. Con el posible control hidráulico
(Río Casanare)
Niveles (m)
Tr (años) Puente Quebrada Casirba
Aguas Arriba Aguas Abajo
2,33 1,17 1,32 1,33 2,15 2,08 2,01 2,00 2,01 2,01 2,04
Puente
50 1,57 1,63 1,89 2,67 2,59 1,38 1,97 2,04 2,03 2,05
100 1,64 1,71 2,00 2,78 2,70 1,45 1,93 2,05 2,03 2,05
El puente sobre la Quebrada Casirba se encuentra a una altura de 5,40m desde el fondo
del cauce de acuerdo a las batimetrías obtenidas, es decir que el gálibo del puente teniendo en
cuenta un periodo de retorno de 100 años y el posible control hidráulico con un nivel máximo de
2,78 m se encuentra a 2,62 m siendo mayor al mínimo recomendado de 2.5 m por Manual de
Diseño Geométrico de Carreteras (Instituto Nacional de vías, 2009).
126
9 CONCLUSIONES
Cuenta con una elevación media de la cuenca de 914,6msnm y una pendiente media de
cuenca del 25,6%, lo cual de acuerdo a la Resolución No. 2965 de septiembre 12 de 1995
emitida por el IGAC se clasifica como una pendiente Fuertemente Quebrado, es decir, son tierras
con diferentes formas: Inclinaciones y longitud de pendientes, las cuales pueden ser cortas o
largas con lomos redondeados, afilados, o ambos con pendientes
127
Aguada y (36010030) Sta Rita, de las cuales se concluye mediante polígonos de thiessen que las
estaciones que más influencia tienen sobre la cuenca en estudio son la estación La Cabuya (97%)
y estación Pte Quemado (3%).
128
Para conocer el comportamiento espacial de la precipitación, se elaboraron las isolineas
de precipitación total media anual para la zona del proyecto para el periodo 1994-2015
correspondientes al periodo común que tienen las cuatro estaciones anteriormente mencionadas y
que cubren todo el trazado en estudio, concluyendo que la precipitación media anual disminuye
en sentido oeste este variando desde los 3100 mm hasta los 1800 mm. Se observa que las
precipitaciones más altas se presentan en las estaciones Pte Quemado. La estación más cercana a
la zona de estudio es La Cabuya y presenta una precipitación total multianual de 2870 mm.
Debido a que sobre la Quebrada Casirba o cerca a ésta no hay estaciones de caudales que
permitan determinarlos de manera directa, se utiliza una metodología indirecta o método regional
de Transposición que utiliza relaciones precipitación - Área – Caudal y el modelo lluvia
escorrentía, los cuales permiten conocer los caudales base característicos de la zona de estudio.
La estación más cercana que registra parámetros de caudal es la estación 36017030 Pte
Carretera, la cual tiene un periodo de registro desde 1997 hasta 2015. Esta estación presenta un
comportamiento monomodal donde los meses con caudal por encima del promedio 9,73 m³/s
son de mayo a octubre. Posteriormente, se realizó un ajuste de distribución para estimar los
caudales esperados para los periodos de retorno 2,33 50 y 100 años para la cuenca Qda Casirba
los cuales se presenta a continuación.
Dentro del análisis de número de curva CN, se siguió el proceso descrito en el Manual
de drenaje de carreteras (Instituto Nacional de vías, 2009), necesario para la caracterización de
las cuencas del sistema hidrológico en estudio, se considera que la cuenca Qda. Casirba
pertenece al grupo de suelo B (CNII 57,4) y C (CNII 71,2) y para condiciones extremas CNIII 86
129
dado que la precipitación es mayor a 52 mm. Esto a través de promedio aritmético de los valores
teniendo en cuenta la clasificación del uso del suelo y su porcentaje de área dentro de la cuenca.
130
La modelación hidráulica desarrollada en el ítem 8 MODELACIÓN HIDRÁULICA, se
realizó en el software en HEC RAS. Para ello se cuenta con un alto porcentaje de incertidumbre
y con información limitada que consiste en 31 secciones batimétricas (ver Figura 8-1) del sitio
en estudio, 16 aguas arriba y 15 aguas abajo del puente con una separación aproximada de 20 m
entre cada sección, y algunas fotografías de la zona. El coeficiente de Manning para el tramo en
estudio se estimó de acuerdo con el procedimiento desarrollado por Cowan (Chow, 1994) 0.09
para las bancas izquierda y derecha y 0.05 para el fondo del canal. De acuerdo al registro
fotográfico se resalta zonas con vegetación en la planicie de inundación, la cual se puede
considerar inefectiva para el flujo.
Teniendo en cuenta que el puente sobre la Quebrada Casirba se encuentra a una altura de
5,40 m desde el fondo del cauce de acuerdo a las batimetrías obtenidas y para un periodo de
retorno de 100 años, el gálibo se encuentra a 2,70 m siendo mayor al mínimo recomendado de
2.5 m por Manual de Diseño Geométrico de Carreteras (Instituto Nacional de vías, 2009).
Para el análisis del posible control hidráulico, se observa que el nivel del Río Casanare es
de 6m, la Quebrada Casirba busca empatarse por energía, lo cual hace que se genere un área de
inundación hasta el puente. Este efecto se considera que no es permanente dado que las láminas
de flujo del Río no son homogéneas y son variables en el tiempo. En cuanto al gálibo del puente
y para un periodo de retorno de 100 años se encuentra a 2,62 m, mayor al mínimo recomendado
de 2.5 m por Manual de Diseño Geométrico de Carreteras (Instituto Nacional de vías, 2009).
131
precisamente para llegar a su estabilidad energética en la entrega al Río. También se puede
apreciar que si el nivel en el Río Casanare es muy bajo las velocidades en la Quebrada Casirba
no van a disminuir.
132
Zona de estudio
Zona de estudio
Figura 9-2. Zonas afectadas por inundación. Departamento del Casanare. (IDEAM,
Instituto de Hidrología, Meteorología y Estudios Ambientales, 2014)
133
en cuenta las modelaciones desarrolladas en las cuales se aprecia una distancia de entre 2,70 y
2,60 siendo mayor al mínimo recomendado de 2.5 m por Manual de Diseño Geométrico de
Carreteras (Instituto Nacional de vías, 2009).
Se podría asumir que el posible evento evidenciado por la comunidad donde la lámina de
agua llegaba al gálibo del puente fue para un periodo de retorno mayor de 100 años y es
aceptable en vista que como lo registra el informe hidrológico diario No. 134, 2016 (IDEAM,
Instituto de Hidrología, Meteorología y Estudios Ambientales, 2014) (ver Figura 9-3), el Río
Casanare ha llegado a presentar niveles de hasta 8 m.
Figura 9-3. Variación de niveles del Río Casanare. Informe hidrológico diario No. 134,
2016
134
10 RECOMENDACIONES
Es necesario instalar una mira en cualquier parte del estribo del puente que permita
registrar los niveles de la Quebrada Casirba, especialmente de los niveles máximos teniendo en
cuenta que no se cuenta con estaciones en la zona o en sus alrededores que registre y provea
información de dicho parámetro o si es posible instalar estaciones. El uso de metodologías de
indirectas de medición puede llegar a tener un porcentaje de incertidumbre asociada a la
precisión.
135
11 ANEXOS
136
12 BIBLIOGRAFÍA
(ONI), B. o. (April de 2019). El Niño and La Niña Years and Intensities. Obtenido de
https://ggweather.com/enso/oni.htm
137
Cudworth, A. G. (1992). Flood Hydrology Manual . Denver, Colorado: Superintendent of
Documents, U.S. Government Printing Office.
Hacienda, M. d. (29 de 12 de 2014). Contrato No. 237 de 2014, para la construcción del puente
vehicular "Lopeño" ubicado en la vía Sácama-La Cabuya, en el departamento de Arauca.
Bogotá, Colombia.
IDEAM. (2019). Estudio Nacional del Agua 2018. Bogotá: Panamericana Formas e Impresos
S.A.
138
Materón M. Hernán y Jiménez E. Henry. (1986). Hidrología Básica, Tomo I, II y III.
Universidad del Valle.
Vargas, R., & Díaz Granados, M. (1998). Curvas Sintéticas Regionalizadas Intensidad -
Duración - Frecuencia para Colombia. Bogotá, D.C.: Universidad de Los Andes.
139
World Meteorological Organization. (1994). Guide to Hidrological Practices, Data Acquisition
and processing, Analysis, Forecasting and Other Applications.
140