Balance Hídrico
Balance Hídrico
Balance Hídrico
1
ESTUDIO HIDROLÓGICO Y BALANCE HÍDRICO PARA DETERMINAR
LA OFERTA Y LA DEMANDA DE AGUA DE LA CUENCA DE LA
QUEBRADA NISCOTA PARA UN ACUEDUCTO INTERVEREDAL EN
NUNCHÍA, CASANARE
2
3
Nota de aceptación
______________________________________
______________________________________
______________________________________
______________________________________
Presidente del Jurado
______________________________________
Jurado
______________________________________
Jurado
4
Agradecimientos
A nuestras familias, por hacer nuestro sueño posible y ser el motor para nosotros.
A los profesores de cada uno de los módulos, quienes con sus conocimientos nos guiaron
5
TABLA DE CONTENIDO
RESUMEN 15
ABSTRACT 16
INTRODUCCIÓN 17
2. MARCOS DE REFERENCIA 22
6
2.1.14 Tiempo de concentración 25
2.1.15 Transpiración 25
2.2 MARCO TEÓRICO 25
2.2.1 Balance hídrico 25
2.2.2 Área de la microcuenca 27
2.2.3 Perímetro 27
2.2.4 Densidad de corrientes 27
2.2.5 Densidad de drenaje (Dd) 27
2.2.6 Longitud del cauce principal (Lcp) 28
2.2.7 Pendiente de la corriente principal 28
2.2.8 Coeficiente de compacidad 28
2.2.9 Factor de forma (Kf) 29
2.2.10 Índice de alargamiento 29
2.2.11 Índice de asimetría 30
2.2.12 Elevación media de la cuenca 30
2.2.13 Pendiente media de la cuenca 31
2.2.14 Escorrentía 31
2.2.15 Evaporación 31
2.2.16 Evapotranspiración 31
2.2.17 Evapotranspiración potencial (ETP) 32
2.2.18 Evapotranspiración real (ETR) 33
2.3 MARCO JURIDICO 33
2.3.1 Constitución política de Colombia 33
2.3.2 Leyes aplicables 33
2.3.3 Decretos 34
2.3.4 Resoluciones 34
2.3.5 Acuerdos 35
2.4 MARCO GEOGRÁFICO 35
2.4.1 Localización del proyecto 36
2.4.2 Hidrografía de la zona de estudio 36
2.4.3 Sitio de cierre de la cuenca 38
2.4.4 Geología de la zona de estudio 38
2.4.5 Estratigrafía 39
2.4.6 Geología estructural 40
7
2.4.7 Climatología del municipio de Nunchía 42
2.4.7.1 Estaciones climáticas y pluviométricas cercanas a la zona de estudio 42
3. METODOLOGÍA 48
4. CARACTERÍSTICAS DE LA CUENCA 51
5. BALANCE HIDRICO 61
8
6.4 DEMANDA HIDRICA TOTAL 67
6.4.1 Método aritmético 68
6.4.2 Método geométrico 69
6.4.3 Método exponencial 69
6.4.4 Calculo de la dotación bruta 72
6.5 INDICE DE ESCASEZ 73
CONCLUSIONES 99
RECOMENDACIONES 102
CITAS 103
BIBLIOGRAFIA 104
9
LISTA DE FIGURAS
10
FIGURA. 21 COBERTURAS DEL SUELO EN LA CUENCA DE LA QUEBRADA NISCOTA. 78
FIGURA. 22 ÁREA SUPERVISADA POR LA ESTACIÓN LA CHAPARRERA Y ÁREA DE LA QUEBRADA
NISCOTA. 84
FIGURA. 23 HIDROGRAMA UNITARIO DE ACUERDO A SCS, PARA LA QUEBRADA NISCOTA. 88
FIGURA. 24 HIDROGRAMA UNITARIO DE ACUERDO A SCS, PARA LA QUEBRADA NISCOTA. 89
FIGURA. 25 ÁREA SUPERVISADA POR LA ESTACIÓN EL PLAYÓN Y ÁREA DE LA QUEBRADA NISCOTA.
90
FIGURA. 26 CURVA DE DURACIÓN DE CAUDALES Q. NISCOTA – SITIO DE CAPTACIÓN. 93
FIGURA. 27 RESULTADOS DEL AFORO REALIZADO PARA LA QUEBRADA NISCOTA 93
FIGURA. 28 GEOMETRÍA QUEBRADA NISCOTA SECTOR DE LA BOCATOMA. 95
FIGURA. 29 CONDICIONES DE FLUJO QUEBRADA NISCOTA SECTOR DE LA BOCATOMA. 96
FIGURA. 30 INUNDACIÓN CON CAUDAL MÁXIMO DE LA QUEBRADA NISCOTA SECTOR DE LA
BOCATOMA. 96
FIGURA. 31 SECCIÓN TRANSVERSAL DE LA INUNDACIÓN CON CAUDAL MÁXIMO DE LA QUEBRADA
NISCOTA SECTOR DE LA BOCATOMA. 97
11
LISTA DE TABLAS
12
TABLA. 21 PERIODOS DE RETORNO SEGÚN DISTRIBUCIÓN DE GUMBEL PARA PRECIPITACIÓN
MÁXIMA ANUAL EN 24HR DE LA ESTACIÓN LA CHAPARRERA, SERIE DE DATOS 1996 2016. 82
TABLA. 22 PRECIPITACIÓN EXCEDENTE PARA DIFERENTES PERIODOS DE RETORNO. 83
TABLA. 23 CÁLCULO DE LOS TIEMPOS DE CONCENTRACIÓN Y DE RETARDO PARA LA CUENCA DE LA
QUEBRADA NISCOTA. 83
TABLA. 24 HIDROGRAMA UNITARIO DATOS. 88
TABLA. 25 HIDROGRAMA UNITARIO DATOS. 89
TABLA. 26 CAUDALES MÁXIMOS MENSUALES ESTACIÓN EL PLAYÓN. 91
TABLA. 27 CAUDALES MÁXIMOS MENSUALES SEGÚN DISTRIBUCIÓN DE GUMBEL PARA DIFERENTES
PERIODOS DE RETORNO DE LA ESTACIÓN LA EL PLAYÓN, SERIE DE DATOS 1996 2016. 91
TABLA. 28 VARIACIÓN DEL ÍNDICE DE ESCASEZ PARA CADA MES DEL AÑO 99
13
LISTA DE ANEXOS
14
RESUMEN
Con este proyecto se busca realizar el balance hídrico de la cuenca de la Quebrada Niscota,
con el fin de conocer su uso potencial como fuente abastecedora de agua para un acueducto que
beneficie seis veredas de la parte baja del municipio de Nunchía en el departamento de Casanare,
esto con el fin de evaluar la oferta hídrica de acuerdo a la resolución 865 de 2004 emanada del
Ministerio de Ambiente, Vivienda y Desarrollo Territorial. Para ello se buscará la información
necesaria concerniente a censos poblacionales, datos meteorológicos oficiales, imágenes satelitales
y mapas geológicos de la cuenca.
Una vez recolectada esta información se procede a su respectivo análisis, para esto se
utilizarán herramientas de sistemas de información geográfica SIG para encontrar los diferentes
parámetros morfológicos de la cuenca y el Software Hec Ras, para determinar los caudales
mínimos y máximos que se pueden esperar para el sitio de interés que corresponde a la captación.
Posterior al análisis de la información se realizará el balance hídrico involucrando las condiciones
morfológicas y climatológicas de la cuenca, el cual nos permitirá conocer la oferta hídrica y
compararla con la demanda hídrica requerida para abastecer a la población beneficiada actual y
futura, durante todos los meses del año. Así mismo es de gran importancia comprobar que se
cumpla con el caudal ecológico luego de captar agua de esta fuente.
15
ABSTRACT
This project seeks to carry out the water balance of the Niscota basin, in order to know
its potential use as a source of water for an exploitation that benefits six trails in the lower part of
the municipality of Nunchía in the Department of Casanare, this in order to evaluate the water
supply according to resolution 865 of 2004 emanating from the Ministry of Environment, Housing
and Territorial Development. For this purpose, it can search for the necessary information for
population censuses, official meteorological data, satellite images and geological maps of the
basin.
Once this information has been collected, its respective analysis can be used, for which
GIS information system tools are used for the different morphological parameters of the basin and
the Hec Ras software, to determine the minimum and maximum flows that can be expected for the
site of interest that corresponds to the catchment. After the analysis of the information, the water
balance is achieved by involving the morphological and climatological conditions of the basin,
which allows us to know the water supply and the comparison with the water demand required to
supply the beneficiary population at present and in the future, during all months of the year. It is
also of great importance to verify that it can be seen with the ecological flow after capturing water
from this source.
16
INTRODUCCIÓN
17
Finalmente se obtiene el índice de escases y caudal ecológico. Como información adicional
se realiza una evaluación de los caudales máximos y mínimos.
18
1. GENERALIDADES DEL TRABAJO DE GRADO
De acuerdo con el contexto en que se desarrolla este trabajo de grado, se puede clasificar
dentro de la línea de investigación de “Saneamiento de Comunidades”.
En el año 2004 se inició la construcción de un acueducto para surtir de agua a las veredas:
Cañas, Corea, Guacharacas, Macuco, Redención, Romero, Sirivana y La Yopalosa del municipio
de Nunchía en el departamento de Casanare; para ello se planeó captar el agua de la Quebrada
Niscota y se construyeron las siguientes obras: bocatoma, aducción, desarenador, tanque de
almacenamiento y parte de la línea de conducción; sin embargo, este proyecto nunca se puso en
funcionamiento por que la bocatoma fue socavada y destruida por la corriente del río. Dada esta
problemática las entidades gubernamentales perdieron el interés en el proyecto; no obstante,
debido a la necesidad de la población de abastecerse de agua potable, especialmente durante la
época de verano, vale la pena evaluar la posibilidad de completar este sistema de acueducto.
Por otra parte; para otorgar el permiso de captación de aguas superficiales, Corporinoquia
exige el estudio completo del sistema de acueducto y el estudio hidrológico de la cuenca, el cual
debe contener el cálculo del índice de escases de acuerdo a la metodología presentada en la
Resolución 865 de 2004.
19
1.2.2 Pregunta de investigación
1.3 JUSTIFICACIÓN
Las autoridades ambientales cumplen una misión muy importante al controlar el uso del
recurso hídrico; por tanto, es deber de los diseñadores de un sistema de acueducto revisar los
caudales necesarios para abastecer a la población beneficiada sin que ello signifique un perjuicio
para el medio ambiente.
20
1.5 OBJETIVOS
1.5.1 Objetivo general
21
2 MARCOS DE REFERENCIA
2.1.2 Caudal
Es el volumen de fluido que atraviesa una superficie dada en la unidad de tiempo. Se usa
para expresar el caudal de un rio, de un fluido por una conducción, etc.
Zona delimitada topográficamente que drena mediante una red hidrográfica, es decir, la
superficie total de tierras que drenan en un cierto punto de un rio o curso de agua.
Curva que muestra el porcentaje de tiempo durante el cual el caudal en un curso de agua
superior a determinados valores, con independencia de su continuidad en el tiempo.
2.1.5 Hidrograma
Gráfica que muestra las variaciones de caudal de un rio o del nivel del agua subterránea en
una perforación respecto a un nivel de referencia frente al tiempo. El Hidrograma unitario es el
22
nombre con el que se conoce un método de cálculo que permite transformar la precipitación en
escorrentía superficial, sirve para facilitar la predicción de cómo responderán cuencas fluviales
concretas ante diferentes supuestos de precipitación. En el Hidrógrafa de caudales se representan
las variaciones de caudal respecto al tiempo.
2.1.6 Hidrología
Ciencia que estudia el ciclo del agua y su evolución sobre la superficie de la tierra y bajo
el suelo, constituyendo una de sus ramas principales el análisis del volumen hídrico de una región
determinada y su distribución en el espacio y el tiempo.
2.1.8 IDEAM
2.1.9 Infiltración
23
2.1.10 Isoterma
Es una curva que une los vértices, en un plano cartográfico, que presentan las
mismas temperaturas en la unidad de tiempo considerada. Así, para una misma área, se pueden
diseñar un gran número de planos con isotermas, por ejemplo: Isotermas de la temperatura media
de largo periodo del mes de enero, de febrero, etc., o las isotermas de las temperaturas medias
anuales.
2.1.11 Isoyeta
Línea curva que une los puntos, en un mapa, que presentan las mismas precipitaciones en
la unidad de tiempo considerada. Así, para una misma área, se pueden diseñar un gran número de
mapas con isoyetas, por ejemplo: isoyetas de la precipitación media del mes de enero, febrero, etc.,
o las isoyetas de las precipitaciones medias anuales, etc.
Los sistemas de captación de agua tienen como propósito recolectar el agua proveniente de
una fuente superficial para su utilización en un acueducto.
Son sistemas que permiten relacionar cualquier tipo de dato con una localización
geográfica. Esto quiere decir que en un solo mapa el sistema muestra la distribución de recursos,
edificios, poblaciones, entre otros datos de regiones o todo un país. Este es un conjunto que mezcla
hardware, software y datos geográficos, y los muestra en una representación gráfica. Los SIG
representan una herramienta muy útil en el procesamiento de datos hidrológicos y climáticos de
una cuenca o región.
24
2.1.14 Tiempo de concentración
Se define como el tiempo que pasa desde el final de la lluvia neta, hasta el final de la
escorrentía directa. Representa el tiempo que tarda en llegar al sitio de aforo la última gota de
lluvia que cae en el extremo más alejado de la cuenca y que circula por escorrentía directa.
2.1.15 Transpiración
Donde,
P = Precipitación.
ES = Escorrentía Superficial.
ETR = Evapotranspiración.
I = Infiltración.
25
El balance hídrico de una cuenca es fundamental para el conocimiento del afluente y su
capacidad de abastecimiento a las poblaciones aledañas que se verán beneficiadas con el recurso,
desde el punto de vista teórico el análisis de este permitirá realizar juicios para el uso del afluente
como el principal y conseguir un uso más racional de los recursos de agua en el espacio y en el
tiempo con la construcción de obras de captación para posterior distribución. El balance hídrico
brindará información acerca de la capacidad del rio y los meses del año en el cual el afluente
aumenta su caudal o por el contrario lo disminuye como también pondrá permitir realizar
predicciones hidrológicas, basado en datos que son registrados por las estaciones meteorológicas
cercanas al área de estudio, información que se solicita al IDEAM.
Existen varios métodos para hallar el balance hídrico de una cuenca y sus componentes, ya
que algunos no tienen en cuenta las posibles variaciones ecológicas ni las características naturales
del medio, como método general se tiene la ecuación del balance hídrico pero “la técnica del
balance hídrico implica mediciones de ambos aspectos, almacenamientos y flujos de agua: sin
embargo, algunas mediciones se eliminan en función del volumen y periodo de tiempo utilizados
para el cálculo del balance” (UNESCO, 1971) es decir, que se reducirá a una sola ecuación en su
forma más general y dependiendo de su complejidad esta podrá reducirse o ampliarse.
Por otro lado, el método para el cálculo de los componentes requiere datos de precipitación
y caudales que son básicos los cuales aportan entidades como el IDEAM ya mencionada con
anterioridad.
26
2.2.2 Área de la microcuenca
2.2.3 Perímetro
Definida como la cantidad de ríos por kilómetro cuadrado, que presenta una cuenca.
Se conoce como densidad de drenaje a la longitud de los ríos y quebradas que llegan o
tributan al río principal dentro del área de la cuenca. Este es un parámetro revelador del régimen y
de la morfología de la cuenca, porque relaciona la longitud de los cursos de agua con el área total.
De esta manera, los valores altos reflejan un fuerte escurrimiento. La longitud total de los cauces
dentro de una cuenca hidrográfica (L), dividida por la superficie total de la Cuenca (A), define la
densidad de drenaje o longitud de cauces por unidad de área. Este parámetro se expresa en
Km/Km².
27
= Ecuación (2).
Como su nombre lo indica, representa el cambio de elevación del perfil del cauce principal,
este parámetro se expresa normalmente en grados y/o en porcentaje. Esta característica contribuye
a definir la velocidad de la escorrentía superficial de la corriente de agua, en donde a mayor
pendiente, mayor velocidad del agua; este parámetro se asocia al régimen hidráulico de la corriente
y a su torrencialidad.
28
unidad, para cuencas cuya forma se asemeja a una circunferencia, el valor de este coeficiente será
muy cercano a uno (1.0). Se obtiene aplicando la ecuación (3).
= . Ecuación (3).
√
Donde:
A = área de la cuenca.
P = perímetro de la cuenca.
= Ecuación (4).
Donde:
A = área de la cuenca.
Lc = Longitud axial de la cuenca.
29
2.2.11 Índice de asimetría
= Ecuación (5).
∑ ( ∗ )
= Ecuación (6).
Donde:
Em = Elevación media de la cuenca.
Ei = Elevación promedio de la i-ésima franja entre curvas de nivel.
A = Área total de la cuenca.
Ai = Área de la i-ésima franja.
30
2.2.13 Pendiente media de la cuenca
Para hallar la pendiente media de la cuenca se utiliza el principio de Alvord, con el cual se
analiza la pendiente media que existente entre curvas de nivel, teniendo en cuenta el área entre las
curvas, su separación y la suma de la longitud de dichas líneas, delimitadas dentro del área de la
cuenca y definida mediante la siguiente expresión:
∗
= Ecuación (7).
Donde:
D = Equidistancia entre curvas de nivel.
Lc = longitud total de las curvas de nivel dentro del área de la cuenca.
A = Área total de la cuenca.
2.2.14 Escorrentía
Es la lámina de agua que circula en una cuenca de drenaje, es decir la altura en milímetros
de agua de lluvia escurrida y extendida dependiendo la pendiente del terreno.
2.2.15 Evaporación
2.2.16 Evapotranspiración
31
2.2.17 Evapotranspiración potencial (ETP)
= Ecuación (9).
= − + + . Ecuación (10).
ETPsc = ETP mensual en mm/mes para mese de 30 días y 12 horas de sol. (Teóricas)
= ∗ ∗ Ecuación (11).
Donde,
ETPsc = Evapotranspiración Potencial sin corregir.
ETP = Evapotranspiración Potencial corregida.
i = Índice de temperatura mensual.
32
a = Factor adimensional.
N = Factor de corrección por latitud y longitud.
33
ambiente. Contiene los principios fundamentales sobre prevención y control de la
contaminación del aire, del agua y del suelo.
Ley 99 de 1993: Por el cual se crea el Ministerio del Medio Ambiente, se reordena
el Sector Público encargado de la Gestión y Conservación del medio ambiente y los
Recursos Naturales Renovables.
2.3.3 Decretos
El decreto 1076 del 26 de mayo de 2015, es el Decreto Único Reglamentario del Sector
Ambiente y Desarrollo Sostenible. Compila en un solo cuerpo normativo todos los decretos
reglamentarios vigentes expedidos hasta la fecha, que desarrollan las leyes en materia ambiental.
Teniendo en cuenta esta finalidad este decreto no contiene ninguna disposición nueva, ni modifica
las existentes.
2.3.4 Resoluciones
34
Resolución 0330 del 8 de junio de 2017: Por el cual se adopta el Reglamento
Técnico para el Sector de Agua Potable y Saneamiento Básico – RAS y se deroga
la Resolución 1096 de 2000, 424 de 2001, 0668 de 2003, 1459 de 2005 y 2320
de 2009.
2.3.5 Acuerdos
35
2.4.1 Localización del proyecto
36
Grande, Salitre, la Rumbita, Negra, Aguablanca, Oreja Marrana, Cueta, San Martín, Castañal, las
Cañas, Niscota, las Coloradas y del Caño los Mangos.
Al occidente del municipio se encuentra la microcuenca del río Pauto el cual nace en el
cerro Romeral en la zona de páramo del municipio de Socotá, en el Departamento de Boyacá, a
3.700 msnm, en Casanare baña los municipios de Támara, Nunchía, San Luis, Pore y Trinidad.
El río Tocaría entrega sus aguas al río Cravo Sur y este junto con el río Pauto entregan sus
aguas al río Meta el cual es uno de los principales afluentes del río Orinoco, uno de los sistemas
hídricos más importantes del mundo.
El municipio de Nunchía cuenta con varias quebradas dentro de las cuales las más
importantes son:
37
Quebrada Chicharra: Nace en la cuchilla la cofradía a 425 m.s.n.m. desemboca
en el río Pauto; presenta una longitud de 7 km.
Quebrada Guevareña: Nace a los 375 msnm, desemboca en la quebrada la
Chicharra, su longitud es de 4.65 Km.
Quebrada Jaboncillo: Nace a los 550 m.s.n.m. desembocando en el margen
izquierdo del río Tocaría. Es intermitente, lo que indica que en épocas de verano
puede secarse en su totalidad.
Quebrada la Picadura: Nace en la cota 675 m.s.n.m. y recorre 2.4 kilómetros.
Quebrada Quiquimare: Nace en las estribaciones de la cordillera oriental en los
límites entre los municipios de Nunchía y paya, desemboca en la quebrada Niscota,
aporta el agua para el acueducto de la vereda La Capilla del municipio de Nunchía.
Quebrada Niscota: Nace en el municipio de Paya en el departamento de Boyacá.
Esta quebrada entrega sus aguas al río Nunchía y este a su vez al río Tocaría.
Quebrada Tamuría: La micro cuenca de la quebrada Tamuría se localiza sobre
territorio de las veredas Alto Tamuría, Barbacoas y cofradía del municipio de
Nunchía, desemboca en el río Pauto. La quebrada Tamuría tiene una longitud de
11.35 Km.
Quebrada el Zaque: Abastece el acueducto del casco urbano, tiene una pequeña
cuenca de 150 hectáreas. Su intervención antrópica ha sido significativa causando
gran deterioro de la cuenca, en épocas de verano su caudal se ha tornado insuficiente
para abastecer a la población.
38
Depósitos fluviales actuales
2.4.5 Estratigrafía
39
2.4.6 Geología estructural
Falla de Guaicáramo
40
Cuellar et al., 1988 (en CARRILLO y VERGARA 1981), mediante estudios de sismicidad
instrumental le atribuyen a esta falla actividad sísmica actual.
Falla de Yopal
Sinclinal de Nunchía
Dueñas y Van Der Hammen (2007), atribuyen el origen de esta estructura a los pulsos
compresivos de las placas tectónicas durante el mioceno, con los cuales se incrementa el
levantamiento de la cordillera oriental y se define el Piedemonte Llanero, producto de ello, se
genera un cinturón deformado estrecho y alargado que daría lugar al sinclinal de Nunchía.
Esta estructura posee una extensión de 60km, presenta una morfología asimétrica y se
encentra afectada en su flanco oeste por el sistema de fallas de Guaicáramo y en su flanco este por
la Falla de Yopal.
41
2.4.7 Climatología del municipio de Nunchía
42
Las estaciones: El Morro, La Chaparrera, Desecho El Hda. Son estaciones pluviométricas
es decir solo registran datos pluviométricos, en las figuras 5, 6 y 7 se pueden apreciar los
pluviogramas realizados con base en los datos del IDEAM para estas estaciones.
300
250
200
150
100
50
0
1 2 3 4 5 6 7 8 9 10 11 12
PREC 2 42.5 63.9 183 343 295 276 234 217 198 137 36.6
43
Figura. 7 Pluviograma estación La Chaparrera Casanare.
Fuente: Autores – Información IDEAM.
44
Figura. 9 Climograma estación Támara, Casanare
Fuente: Autores – Información IDEAM.
45
Figura. 10 Temperatura y humedad relativa de la estación Támara.
Fuente: Autores – Información de IDEAM.
46
Figura. 12 Temperatura y humedad relativa de la estación Corinto.
Fuente: Autores – Información de IDEAM.
Dado que la cuenca de la Quebrada Niscota posee elevaciones del terreno que oscilan entre
350 y 1300 m.s.n.m. presenta diferentes intensidades de lluvia y temperatura para un período de
tiempo dado. Teniendo en cuenta que los datos recolectados en las estaciones meteorológicas son
puntuales es necesario distribuirlos en toda el área de la cuenca, para ello se emplean métodos
gráficos y matemáticos, para la realización de este estudio se utiliza las herramientas de sistemas
de información geográfica SIG para elaborar mapas de isolíneas y encontrar los parámetros
característicos de la cuenca, teniendo como base la información de las estaciones meteorológicas
de IDEAM.
47
3 METODOLOGÍA
mediante isolíneas.
48
3.1.2 FASE 2: Análisis de la zona en estudio.
evapotranspiración de la cuenca.
De los mapas geológicos obtener la geología regional y local del sitio de estudio.
49
3.1.4 FASE 4: Resultados:
Análisis del área de inundación ocasionada por caudales máximos para un periodo
de retorno de 100 años.
Verificar los resultados del aforo realizado al final del verano frente a los datos de
caudal encontrados mediante métodos indirectos para la Quebrada Niscota.
enfocándose en el balance hídrico, Para el desarrollo de este proyecto se realizarán visitas al sitio
se llevará a cabo el proyecto con el fin de describir la zona de estudio y registrar la información
Excel, Hec RAS y ArcGIS, los cuales serán las herramientas que se utilizarán para el
50
4 CARACTERÍSTICAS DE LA CUENCA
51
un área total de 6987 hectáreas equivalentes a 69.87 Km2 que aportan aguas superficiales en los
periodos lluviosos y de infiltración acumulada para los periodos secos.
La quebrada Niscota en su cauce principal presenta una longitud total de 16.5 Kilómetros
desde su nacimiento hasta el sitio de la captación, siguiendo la sinuosidad del cauce.
Para el caso de la quebrada Niscota se obtiene una pendiente media de 5.9% desde el sitio
de captación hasta su nacimiento, como se observa en la siguiente Figura.
52
4.1.3 Coeficiente de Gravelius o de compacidad (kc).
Se obtiene un coeficiente de Gravelius igual a 1.44 para la cuenca de interés, lo que indica
una forma de oval redonda a oval oblonga con propensión a crecientes moderadas.
Por medio de sistemas de información geográfica se determinó que la cuenca tiene una
longitud axial máxima de 13036m.
Para la cuenca de interés el coeficiente de forma igual a 0.41, lo cual indica una forma
alargada con propensión a crecientes moderadas más aun teniendo en cuenta que el área de interés
es la parte inicial de la cuenca, donde el volumen captado no alcanza a producir crecientes
importantes.
53
4.1.5 Índice de alargamiento
54
Tabla. 5 Cálculo de la elevación media por el método de las franjas de área.
Fuente: Autores.
Intervalo entre Ei Ai EixAi %A Curva Hipsométrica
2 2
curvas de nivel (m) (Km ) Km % Aacum Elevación (m)
1800 1883,5 1841,8 0,23 431,0 0,34 0,34 1841,8
1600 1800 1700 1,15 1955,0 1,65 1,98 1700
1400 1600 1500 4,63 6946,5 6,63 8,61 1500
1200 1400 1300 8,56 11130,6 12,25 20,86 1300
1000 1200 1100 12,34 13571,8 17,66 38,52 1100
800 1000 900 16,73 15060,6 23,95 62,47 900
600 800 700 15,13 10592,4 21,66 84,12 700
400 600 500 8,75 4374,0 12,52 96,64 500
400 335 367,5 2,35 861,8 3,36 100,00 367,5
69,87 64923,7 100
Elevación media (m) = 929,2
55
4.1.8 Pendiente media de la cuenca – Método de Alvord
∗
=
Dada su pendiente media se puede afirmar que la cuenca posee un relieve fuertemente
accidentado.
56
4.2 PARÁMETROS RELACIONADOS CON EL DRENAJE DE LA MICROCUENCA
4.2.1 Orden de drenaje
Esta microcuenca presenta una densidad de corrientes de 168 ríos / 69.87Km² = 2.40
ríos/Km².
57
4.2.3 Densidad de drenaje (Dd)
Dd = 2.14 Km/Km2
Teniendo como base la información de las precipitaciones medias anuales de las estaciones
cercanas al proyecto y con las ayudas de los SIG, se elabora el mapa de isoyetas para caracterizar
la cuenca con una precipitación media anual de 3120 mm/año.
58
Figura. 17 Mapa realizado con herramientas de SIG para las Isoyetas de la
microcuenca, precipitación media anual.
Fuente: Autores.
59
Figura. 18 Mapa realizado con SIG para las Isotermas de la microcuenca,
temperatura media anual.
Fuente: Autores.
60
5 BALANCE HÍDRICO
El balance hídrico es una forma de estimar la cantidad de agua que puede presentar una
región en un periodo de tiempo determinado. Usando un método numérico se comparan las
entradas de agua lluvia con las salidas (evaporación y escorrentía).
Por otra parte, se realizó la extracción de la información mensual dentro del área de la
cuenca hidrográfica a partir del geo-procesamiento de las Capas Ráster, esto con el fin de
cuantificar la distribución media de los parámetros climatológicos. Para este caso, se utilizó el
módulo de estadísticas zonales incluidas dentro de las herramientas de análisis espacial de
ArcMap.
61
Tabla. 9 Valores de precipitaciones y temperaturas medias mensuales obtenidos
mediante el procesamiento de mapas de SIG.
Fuente: Autores.
62
Teniendo los valores mensuales de evapotranspiración corregida se procede a elaborar el
balance hídrico de la cuenca.
63
Figura. 19 Representación gráfica del balance hídrico mensual
Fuente: Autores.
64
6 OFERTA, DEMANDA E INDICE DE ESCASEZ DE LA CUENCA
Para este caso la reducción de la oferta hídrica por caudal ecológico es:
65
6.1.1.2 Reducción por calidad del agua
La oferta del recurso hídrico se ve afectada por la calidad del agua, esta puede ser alterada
por las actividades humanas o por el escurrimiento de bosques o por el escurrimiento de aguas en
zonas de producción agrícola y ganadera, también el aumento de carga sólida en el agua por
descarga de sedimentos se manifiesta con mayor frecuencia en las corrientes de régimen torrencial
y depende en esencia de la intensidad de la lluvia en la parte alta de las cuencas, que interactúa con
el sistema cobertura vegetal y suelo en dichas cuencas. Esta alteración está también directamente
relacionada con las actividades humanas.
Una vez se conozca el estado de la calidad del agua de las fuentes de abastecimiento como
de los cuerpos de agua, la oferta hídrica de estos sistemas se debe afectar por el 25%,
correspondiendo a la condición de calidad del agua.
En este caso la reducción de la oferta hídrica por calidad del agua es: 0.25*110038263m3/año =
27509566 m3/año.
Para obtener la oferta hídrica neta se debe tener en cuenta la reducción por calidad del agua
(25%) y la reducción por caudal ecológico (25%); por lo tanto, la oferta hídrica neta de la cuenca
se puede expresar como:
Demanda Potencial de Agua: volumen potencial de agua, que se requiere para atender las
actividades socioeconómicas en un espacio y tiempo determinado. (Incluye la cantidad de agua
abastecida y contabilizada, agua usada no contabilizada y el agua requerida para desarrollar
actividades socioeconómicas deprimidas y/o no abastecidas de manera efectiva). (IDEAM, 2008).
66
6.2 INDICE DE ESCASEZ
El Índice de Escasez de agua es la relación porcentual entre la demanda potencial de agua
del conjunto de actividades sociales y económicas con la oferta hídrica disponible, luego de aplicar
factores de reducción por régimen de estiaje y fuentes frágiles. (IDEAM, 2008)
= % Ecuación (12).
Donde:
Ie = Índice de escasez (%)
Dh = Demanda hídrica potencial de agua (Mm3)
On = Oferta hídrica superficial neta (Mm3)
Fr = factor de reducción.
Donde:
67
DUS = Demanda de Agua para el Sector Servicios.
DUA = Demanda de Agua para Uso Agrícola.
DUP = Demanda de agua para uso pecuario.
Como demanda doméstica se tiene la demanda necesaria para suplir de agua a la población
beneficiada por el acueducto Inter veredal de la Quebrada Niscota, con el cual se espera abastecer
a 2320 habitantes censados para el escenario actual. A continuación, se realiza una proyección de
la población actual con un horizonte de 25 años, para ello se utilizan tres métodos recomendados
por el Reglamento Técnico para el Sector de Agua Potable y Saneamiento Básico – RAS.
Este método sugiere que la población aumenta a una tasa constante de crecimiento
aritmético, representada de la siguiente forma:
− −
Tasa de crecimiento: = = Ecuación (14).
−
= + Ecuación (15).
Donde,
m = Período intercensal entre los censos P1 y Po
t1 = Año del censo P1
t2 = Año del censo Po
k = Tasa de crecimiento aritmético
P = Población al final del período de diseño
n = Período comprendido
Para proyectar la población por medio de este método se tienen los datos de los censos
realizados por el DANE para el año 1993 con 6894 habitantes y para el año 2005 con 8420
habitantes para todo el municipio de Nunchía. Tomando como dato base la población actual a
beneficiarse con el acueducto Interveredal igual a 2320 habitantes, es posible proyectarla por el
método aritmético con la misma tasa de crecimiento.
68
De esta forma se obtiene una tasa de crecimiento de 1.51% para el municipio de Nunchía,
aplicando este método a la población contabilizada, beneficiada con el proyecto en el año 2018 de
2320 habitantes se obtiene para el año 2043 la población crecerá a 2776 habitantes.
Pf = Puc * (1 + r) (Tf-Tuc)
Donde,
Pf = Población proyectada
Puc = Población del censo inicial
Tf = Año de la proyección
Tuc = Año del último censo
r = Tasa de crecimiento anual
Tuc - Tci
Donde,
Ka = Pendiente de la recta
Pci = Población censo inicial
Puc = Población del último censo
Tci = Año del censo inicial
69
Tuc = Año del último censo
r = Tasa de crecimiento anual
70
Tabla. 14 Proyección de la población beneficiada con el acueducto Interveredal, teniendo
como fuente abastecedora la Quebrada Niscota.
Fuente: SELMAC S.A.S. Consultores.
PROYECCIÓN DE CRECIMIENTO EXPONENCIAL Y GEOMÉTRICO DE
LA POBLACIÓN
Tasas de Crecimiento Fuente DANE para el Departamento de Casanare (%)
(Periodo 2015-2020)
POBLACION INICIAL CENSO (2320
hab) USUARIOS 464
AÑO METODO Población de saturación 0
GEOMETRICO EXPONENCIAL T. Geométrico T. Exponencial
2018 2320 2320 2320 2320
2019 2359 2360 2359 2360
2020 2399 2400 2399 2400
2021 2440 2441 2440 2441
2022 2481 2483 2481 2483
2023 2523 2526 2523 2526
2024 2566 2569 2566 2569
2025 2609 2613 2609 2613
2026 2653 2658 2653 2658
2027 2698 2704 2698 2704
2028 2744 2750 2744 2750
2029 2790 2797 2790 2797
2030 2837 2845 2837 2845
2031 2885 2894 2885 2894
2032 2934 2944 2934 2944
2033 2984 2994 2984 2994
2034 3034 3045 3034 3045
2035 3085 3097 3085 3097
2036 3137 3150 3137 3150
2037 3190 3204 3190 3204
2038 3244 3259 3244 3259
2039 3299 3315 3299 3315
2040 3355 3372 3355 3372
2041 3412 3430 3412 3430
2042 3470 3489 3470 3489
2043 3529 3549 3529 3549
71
DANE Casanare indicadores
r=1,69% anual r=1,70% anual
demográficos (cuadro A81)
SELECCIÓN 3549 HABITANTES PARA EL 2043
Tabla. 15 Dotación neta máxima por habitante según la altura sobre el nivel del mar de la
zona atendida.
Fuente: Decreto 330 de 2017 Ministerio de Ambiente, Vivienda y Desarrollo Territorial.
Por lo tanto, la demanda por uso doméstico corresponde a la demanda del acueducto Inter
veredal, la cual es la siguiente:
72
dbruta Dneta
1 % P Ecuación (18).
Los términos DUI, DUS, DUA y DUP de la ecuación (13), tienen un valor de cero debido
a que no existe demanda hídrica para los sectores industrial, servicios, el uso de agua para los
sectores agrícola y pecuario de esta fuente son muy bajos.
= % Ecuación (19).
= . %
73
7 ANÁLISIS DE CAUDALES MÁXIMOS MEDIOS Y MINIMOS
Para el estudio de caudales que transitan por el punto de estudio de la cuenca de la Quebrada
Niscota, se tiene en cuenta la información de las estaciones más cercanas, las cuales son: Estación
El Playón código 35217060 de tipo limnigráfica y la estación La Chaparrera código 35210050 de
tipo pluviográfica, las cuales se encuentran dentro de la cuenca del río Tocaría de la cual es afluente
la quebrada Niscota. Debido a que la estación La Chaparrera registra precipitaciones y la estación
El Playón registra caudales, se realiza una transformación de la lluvia registrada en la estación La
Chaparrera en escorrentía superficial y luego por el método del hidrograma unitario del SCS, se
obtienen los caudales máximos y mínimos de la cuenca los cuales se pueden comparar con los
observados en la estación limnigráfica El Playón, para finalmente realizar una modelación de
caudal máximo utilizando el programa HEC RAS, 4.2.1 para analizar las inundaciones en la zona
aledaña al sitio de captación.
Con los caudales mínimos mensuales encontrados se evaluará la oferta y la demanda más
crítica durante el año.
Para emplear este método de transformación se debe estimar; la infiltración del agua en el
suelo además de los tiempos de concentración tc y de retardo tL de la cuenca
74
7.2 INFILTRACIÓN
La infiltración se define como el proceso por el cual el agua penetra el suelo desde la
superficie hasta las capas inferiores del suelo. La capacidad de infiltración es la capacidad máxima
con que el suelo en una condición dada puede absorber agua.
Las condiciones del suelo afectan la capacidad de infiltración, es así como el uso del suelo
y la humedad del mismo, afectan directamente su estado.
75
7.3.1 Clasificación Hidrológica de los Suelos
Suelo tipo A: Potencial de escurrimiento bajo. Suelos con altas capacidades de infiltración
cuando están completamente húmedos, principalmente arenas y gravas muy bien
ordenadas. Suelos con alta transmisión de agua.
Suelo tipo C: Suelos con capacidades de infiltración baja cuando están completamente
húmedos, principalmente suelos que contiene una capa que impide el movimiento hacia
abajo o suelos con textura fina o moderadamente fina. Estos suelos tienen baja transmisión
de agua.
Suelos tipo D: Suelos con capacidades de infiltración muy bajas cuando están
completamente húmedos. Suelos que se expanden significativamente cuando se mojan,
arcillas altamente plásticas y ciertos suelos salinos. Suelos con transmisión del agua muy
baja. Tomado de la Resolución 865 de 2004, emanada por el Ministerio de Ambiente,
Vivienda y Desarrollo Territorial.
76
De acuerdo a la clasificación del SCS, los suelos de la cuenca se pueden clasificar como: Suelo
tipo C, con esta clasificación se eligen los valores del número de curva (CN) de la siguiente tabla.
77
Con la ayuda de ArcGIS se elabora un mapa de coberturas del área de la cuenca en estudio,
a partir de una imagen satelital de 12.5X12.5 de resolución, se identifica cada cobertura y se
clasifica de acuerdo a la metodología Corine Land Cover a escala 1:25.000, luego de la tabla de
atributos se extrae el área de cada cobertura, la cual se presenta en términos de porcentajes en la
Tabla 17, con la cual se calcula el valor del CNII ponderado.
78
Tabla. 17 Valor del número de curva CNII ponderado para el área de la cuenca.
Fuente: Autores.
COBERTURA AREA PORCENTAJE CN CN*Ai
Herbazal denso de tierra firme no arbolado 17,19 24,61 86 1478,3
Herbazal denso de tierra firme con arbustos 4,75 6,80 86 408,5
Pastos limpios 0,44 0,63 70 30,8
Pastos arbolados 0,14 0,20 70 9,8
Pastos enmalezados 0,2 0,29 70 14,0
Mosaico de cultivos, pastos y espacios naturales 15,76 22,57 70 1103,2
Bosque denso alto de tierra firme 2,1 3,01 60 126,0
Arbustal denso 0,12 0,17 70 8,4
Tierras desnudas y degradadas 0,2 0,29 80 16,0
Zonas pantanosas 6,49 9,29 70 454,3
Red vial y territorios asociados 0,2 0,29 82 16,4
Rios 0,44 0,63 60 26,4
Mosaico de pastos y espacios naturales 21,81 31,23 75 1635,8
Σ 69,84 100,00 5327,9
NUMERO DE CURVA CN = 76,3
I 0 - 36
II 36.1 – 52.5
79
Para la estación La Chaparrera código 35210050, el día más lluvioso para la serie de datos
estudiada (1996 a 2016) se presentó el día 24 de mayo de 2016, como se puede apreciar en la Tabla
19.
= . .
Ecuación (20).
76.3
CN = = 88
0.43 + 0.0057 ∗ 76.3
80
P = Pe + Ia + Fa
Dónde:
P= precipitación total
Pe = exceso de precipitación;
Ia = (abstracción inicial antes de encharcamiento);
Ia = 0.2S
Fa = abstracción continuada,
Desarrolló la ecuación de la precipitación efectiva, la cual se expresa de la siguiente forma:
( . )
= ( . )
Ecuación (21).
[ ( ) ]
= [ (( ∗ ) ) ]
≥ ∗ (200 ) − 2; = 2.54 Ecuación (22).
Dónde:
= − Ecuación (23).
81
Tabla. 20 Restitución de datos de precipitación máxima en 24hr estación la Chaparrera.
Fuente: Autores – Datos de IDEAM.
AÑO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
1996 2 0 19,1 59,5 55,2 134,4 44,6 97,7 78 43,6 34 52 134,4
1997 0 34,8 39,5 47 45,3 36,5 91 94,3 69,9 40,9 41 0 94,3
1998 0 34,5 35,7 79,7 132,5 57 52,6 130 127 55 44,5 6,5 132,5
1999 5,5 79,6 20,5 67 68,9 36,6 107 117 60 45 54 27,5 117
2000 7,2 7,2 84,6 63,5 87,5 73 84,8 64 57,6 51 67,5 130 130
2001 0 0 9 74,5 81,4 73,1 73 89,1 74,6 52,9 47,9 41 89,1
2002 0 0 34,3 65,5 105 94 61 51 73 40 39 39,8 105
2003 3,2 22,6 33,9 66,5 87,5 62 46 97,5 47 34,5 28,5 31 97,5
2004 0 56,8 33,6 97 90 54 91,5 81 82 58,4 51,3 39,7 97
2005 17 33,5 11,2 105 72 95,8 87 84,5 66 75,7 28,5 6 105
2006 10,5 0 48,8 33,5 116,5 60 80 56 43,2 103,2 137 0 137
2007 0 15 39,5 32 41 31 67 28,5 85,5 68 36,5 39,5 85,5
2008 0 0 21,3 49 80 141 132 38,5 118,5 55 20 5,5 141
2009 5 4 58 46,5 86 70 43 112 55,2 48,5 30 7,1 112
2010 0 10 40 39,6 122 108,9 49 106 36 44,1 119 17 122
2011 9 42 11 123 108 108,5 69,5 97,5 48,5 143,8 60 37 143,8
2012 0 17,9 56,5 75,4 60 72 121 140 84 62,7 42,5 60 140
2013 0 5,5 33,9 51 62 56 88 57 76,5 47 34,5 11 88
2014 0 33,5 25,5 130 89,3 68 79 47 27,9 79 65,5 29 130
2015 5,6 77 20 23,5 70 82,4 85 53 67 102 68,5 123 123
2016 0 4,2 14 130 145 57,5 130 48,5 83 53 76 39 145
MEDIOS 3,10 22,77 32,85 69,46 85,96 74,84 80,10 80,48 69,54 62,06 53,60 35,31 117,58
MAXIMOS 17,0 79,6 84,6 130,0 145,0 141,0 132,0 140,0 127,0 143,8 137,0 130,0 145,0
MINIMOS 0,0 0,0 9,0 23,5 41,0 31,0 43,0 28,5 27,9 34,5 20,0 0,0 0,0
82
Aplicando las ecuaciones (21) y (22) se obtienen los siguientes valores para la precipitación
excedente.
83
Figura. 22 Área supervisada por la estación La Chaparrera y área de la Quebrada
Niscota.
Fuente: Autores.
84
7.5.2 Hidrogramas
Un hidrograma es una gráfica que representa la variación del caudal con relación al tiempo;
el área bajo la curva representa el volumen de agua que pasa por la sección del río en un intervalo
de tiempo del hidrograma. Monsalve. Op.cit., p. 194-195
Los hidrogramas sintéticos se utilizan para desarrollar hidrogramas unitarios para otros
puntos en la corriente dentro de la misma cuenca o para cuencas cercanas de características
similares. Existen tres tipos de hidrogramas unitarios sintéticos: aquellos que relacionan las
características del hidrograma (tasa de flujo pico, flujo base, etc.) con las características de la
cuenca (Snyder, 1938; Gray, 1961), aquellos basados en hidrogramas unitarios adimensionales
(SCS, 1972), y aquellos basados en modelos de almacenamiento en la cuenca (Clark, 1943).
85
.
= ∗( ∗ ) Ecuación (25).
L= longitud del cauce principal desde la divisoria de aguas hasta la salida (m).
Lc = longitud del cauce principal desde el punto más cercano al centroide de la cuenca hasta la
salida (m).
Ct = coeficiente asociado al almacenamiento de la cuenca.
La fórmula de Snyder para caudal pico es:
∗
= Ecuación (26).
Suponiendo por simplicidad una lluvia efectiva uniforme, el tiempo al pico del hidrograma
unitario es igual a la mitad de la duración de la lluvia efectiva más el tiempo de desfase. Por
consiguiente, el tiempo al pico en términos del tiempo de desfase es:
= Ecuación (27).
Este hidrograma unitario sintético fue desarrollado por Víctor Mockus* en 1950. Este
Hidrograma fue estudiado con base en el análisis de un gran número de hidrogramas unitarios
naturales de un amplio rango de tamaños de cuencas hidrográficas y sitios geográficos.
En el método del SCS se calcula el tiempo de desfase por medio de la siguiente formula:
. ( . ) .
= ∗ . . Ecuación (28).
Donde,
86
S = Pendiente del cauce principal en (m/m).
= Ecuación (29).
Donde,
Para el tiempo al pico se tiene en forma experimental la relación tp/tr =5, el SCS teniendo
en cuenta estas ecuaciones propone que el caudal pico se puede encontrar con la siguiente fórmula:
. ∗
= Ecuación (30).
Donde,
87
Figura. 23 Hidrograma Unitario de acuerdo a SCS, para la quebrada Niscota.
Fuente: Autores.
88
Figura. 24 Hidrograma Unitario de acuerdo a SCS, para la quebrada Niscota.
Fuente: Autores.
89
7.6 CAUDAL MÁXIMO OBSERVADO EN LA ESTACIÓN EL PLAYÓN
90
Tabla. 26 Caudales máximos mensuales estación El Playón.
Fuente: Autores – Datos de IDEAM.
AÑO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC
1996 35 61 35 227 555 659 532 340 699 1249 364 518 1249
1997 44 42 37 227 474 1176 863 460 358 426 82 19 1176
1998 12 30 156 615 600 695 1125 790 1125 311 305 344 1125
1999 15,4 81,6 24,2 212,6 302 681 834,1 927,2 749,4 753 232,2 105 927,2
2000 19,4 21 517,6 158 436,2 278 717 681 827 314 232,2 97,8 827
2001 31,9 12,6 18,6 52,5 570 827 481,6 446,4 364,4 284 336,8 357,5 827
2002 31,9 12,7 63 378,2 287 753 385,1 844,8 532 323 681 58,5 844,8
2003 13,6 10,3 34,8 235 949,1 485,2 439,6 532 314,2 375,8 382,6 162,7 949,1
2004 23 26 71 1103 570 547 670 532 597 366 224 101 1103
2005 22 98 12 801 1008 938 266 1004 521 1154 416 34 1154
2006 207 170,3 226,6 266,2 1100 644,5 693,7 687,5 555 595,3 420,4 108,6 1100
2007 55,5 52,5 235 691,8 474,4 485,2 470,8 1154 670 1118 254,6 76,5 1154
2008 42,5 32,6 54 204,6 532 702,6 1008 409 923,5 532 790 85,2 1008
2009 37 35,9 40,3 162,9 830,5 738,6 702,6 1060 521,2 1103 600,4 54 1103
2010 31,5 42,5 282,2 938,2 1444 1474 1356 1186 558,6 830,5 827 147 1474
2011 34,8 35,9 99,6 808,5 862,5 855,4 295 532 551 666,4 1118 207 1118
2012 41,4 34,8 108,6 558,6 680,9 729 715,2 246,2 478 485,2 112,2 43,6 729
2013 10,72 14,29 49,85 162,8 284,3 230 451,3 434,8 256,2 409,7 248,3 72,91 451,3
2014 15,57 36,11 39,96 312,3 403,4 793,4 873,1 535,4 203,8 292,7 342,1 86,49 873,1
2015 33,24 51,92 39,17 242,4 352,5 745,1 959,6 907,1 459,6 352,5 313,6 81,35 959,6
MEDIOS 37,87 45,10 107,22 417,88 635,79 721,85 691,94 685,47 563,20 597,06 414,12 138,01 1007,61
MAXIMOS 207,0 170,3 517,6 1103,0 1444,0 1474,0 1356,0 1186,0 1125,0 1249,0 1118,0 518,0 1474,0
MINIMOS 10,7 10,3 12,0 52,5 284,3 230,0 266,0 246,2 203,8 284,0 82,0 19,0 451,3
Según los reportes de IDEAM del año 2010 el caudal máximo presentado en el río Tocaría
medido en la en la estación 35217060 El Playón en el municipio de Nunchía, registra un caudal
máximo absoluto de 1474 m3/s para el mes de junio. Calculando los caudales para diferentes
periodos de retorno por medio de la distribución de Gumbel, se encuentran los valores de la Tabla
27, donde para un periodo de retorno de cien años el caudal máximo probable que puede transitar
por el río Tocaría en el sitio de ubicación de la estación El Playón es de 2085m 3/s.
Teniendo en cuenta que la cuenca de la Quebrada Niscota de 69.87 Km2 hace parte de la
cuenca del río Tocaría, cuyo caudal es monitoreado por la estación El playón, sitio en el cual
91
confluyen las aguas captadas en un área de 1393 km2, es posible relacionar el caudal máximo
medido en la cuenca del río Tocaría con el caudal máximo posible para la cuenca de la Quebrada
Niscota para el mismo intervalo de tiempo, relacionando el caudal con la raíz cuadrada del área,
de la siguiente manera:
Q1 √ 1
Q2 √ 2,
De la relación anterior se obtiene que el caudal máximo probable registrado para la cuenca
de la Quebrada Niscota para un periodo de retorno de cien años es 466.9 m 3/s.
Esta curva expresa el porcentaje de tiempo en el cual el caudal que fluye por un curso de
agua en un sitio determinado es superior a un determinado valor.
Teniendo en cuenta los datos de caudales medios observados en la estación El Playón para
el periodo 2010 a 2015, se construye la curva de duración de caudales para el río Tocaría (cuyo
cálculo se encuentra consignado en el Anexo C). En forma análoga, se obtiene la curva de duración
de caudales para la Quebrada Niscota, relacionando el caudal y el área de las dos cuencas de la
siguiente forma:
Q1 A1
Q2 A2
De esta manera se obtiene la curva de duración de caudales para la quebrada Niscota, como se
observa en la Figura 26.
92
Figura. 26 Curva de Duración de Caudales Q. Niscota – Sitio de Captación.
Fuente: Autores.
93
duración de caudales de la Figura 26 se observa que durante el 98% del tiempo se presentan
caudales que superan este valor.
En resumen los caudales máximos encontrados para la Quebrada Niscota son: En primer
lugar 385.3 m3/s, obtenidos mediante la transformación de la precipitación en escorrentía y la
construcción de un hidrograma unitario según el SCS, a partir de los datos de precipitación máxima
en 24Hr, medidos en la estación la Chaparrera, en segundo lugar 466.9 m3/s obtenidos de los datos
registrados directamente en la estación El Playón y una relación entre el área de la cuenca de la
quebrada y el área vigilada por la estación. De los dos caudales encontrados se utiliza el mayor
para realizar una modelación de la inundación del área aledaña a la bocatoma, sitio de interés del
presente estudio, con la ayuda del software HEC RAS 4.2.1
94
Figura. 28 Geometría Quebrada Niscota sector de la bocatoma.
Fuente: Autores, tomado de HEC-RAS
Posteriormente se definen las condiciones de flujo, para este caso realizamos la modelación
con flujo estacionario dado que contamos con el caudal máximo posible para un periodo de retorno
de 100 años igual a 466.9 m3/s, los cuales empiezan a fluir aguas arriba desde la estación 3600,
como se observa en la Figura 30.
95
Figura. 29 Condiciones de Flujo Quebrada Niscota sector de la bocatoma.
Fuente: Autores, tomado de HEC-RAS
96
Figura. 31 Sección transversal de la inundación con caudal máximo de la Quebrada Niscota
sector de la bocatoma.
Fuente: Autores, tomado de HEC-RAS
Abriendo la tabla standard se pueden apreciar los valores numéricos de los parámetros
utilizados y generados por la ejecución de HEC-RAS como son: Caudal utilizado, cota mínima del
terreno, elevación de la lámina de agua, elevación del agua para flujo crítico, pendiente del tramo,
velocidad, área de flujo, ancho superior y número de Froude.
97
Tabla. 27 Parámetros obtenidos con HecRas, para caudal máximo.
Fuente: Autores – Datos de IDEAM.
98
CONCLUSIONES
Tabla. 28 variación del índice de escasez para cada mes del año
Fuente: Autores
Parámetro Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic Anual
Precipitación (mm) 12,8 52,9 98,1 294,1 419,3 391,0 374,7 338,4 314,4 289,0 167,0 54,9 2806,5
Area de la cuenca Km2 69,87 69,87 69,87 69,87 69,87 69,87 69,87 69,87 69,87 69,87 69,87 69,87
Oferta hídrica mensual m3 891541,2 3696821,7 6854247,0 20547369,6 29292997,5 27319868,7 26180987,7 23645405,4 21964333,2 20189635,2 11671084,8 3836561,7 196090854
Demanda hídrica mensual m3 20529,5 18542,8 20529,5 19867,3 20529,5 19867,3 20529,5 20529,5 19867,3 20529,5 19867,3 20529,5 241719
Reducción por caudal ecológico 222885,3 924205,4 1713561,8 5136842,4 7323249,4 6829967,2 6545246,9 5911351,4 5491083,3 5047408,8 2917771,2 959140,4 49022713
Reducción por calidad 222885,3 924205,4 1713561,8 5136842,4 7323249,4 6829967,2 6545246,9 5911351,4 5491083,3 5047408,8 2917771,2 959140,4 49022713
Índice de escasez % 4,61 1,00 0,60 0,19 0,14 0,15 0,16 0,17 0,18 0,20 0,34 1,07
De acuerdo con el balance hídrico realizado para la cuenca de la Quebrada Niscota, desde
febrero hasta abril, todos los valores de Δ corresponden a déficits, dado que la evapotranspiración
supera el valor de las precipitaciones medias mensuales; siendo enero el mes que presenta un
mayor déficit de agua en la cuenca, lo cual confirma los valores observados para el índice de
escasez mensual presentados en la Tabla 28.
99
En abril, la precipitación es mayor que la evapotranspiración (P>ETP) y se comienza a reponer
humedad en el suelo, sin embargo, el suelo aún posee capacidad de almacenar agua y durante los
meses de mayo a noviembre el suelo ya está saturado por tanto el Δ de Almacenamiento es cero;
para el mes de diciembre tenemos, 100 mm de variación de almacenamiento menos 68.5mm de
diferencia para un total de 31.5mm; por lo tanto, diciembre no presenta Déficit, debido a que se
utilizan 68mm acumulados durante el mes de noviembre.
Al comparar el caudal máximo encontrado para la Quebrada Niscota por medio del
hidrograma unitario de la SCS involucrando el método del número de curva CN de 385.3 m3/s y
el caudal encontrado mediante los caudales medidos en la estación más cercana y relación de áreas
de cuencas de 466.9 m3/s, se observa una diferencia del 17%, esto refleja el grado de incertidumbre
con el que se cuenta dado que se están utilizando métodos indirectos para el encontrar estos
caudales dado que no se cuenta con una estación hidrológica instrumentada dentro de la cuenca de
estudio.
100
encontrados indirectamente para la quebrada, además se observa que aún en la época más seca del
año el agua que escurre por el cauce principal de la cuenca es suficiente para abastecer el proyecto
del acueducto interveredal sin experimentar presiones sobre la cuenca incluso cuando se presentan
los caudales mínimos.
101
RECOMENDACIONES
Para una mejor comprensión y para mayor precisión del comportamiento de caudales se
deben utilizar información de estaciones hidrometeorológicas con series de datos de muchos más
años.
No hay inconveniente para que la corporación autorice la concesión de captación de agua
superficial, dado que no se está afectando en forma significativa el hábitat de la fauna y la flora
con la captación del agua que demanda el acueducto Inter veredal de Nunchía.
Realizar el estudio de socavación del sitio de ubicación de la captación, para tener una
mejor comprensión de la profundidad de la cimentación de una estructura fundada en el lecho del
río, para ello se pueden usar las metodologías propuestas por varios autores como, por ejemplo:
Lischtvan-Lebediev (1959), Laursen (1963), Maza Álvarez & García Flores (1978) y Neill (1980),
entre otros.
Otra opción es realizar una modelación con el Software HEC-RAS, utilizando el módulo
de sedimentos el cual permite calcular el transporte de sedimentos mediante diferentes métodos, y
para condiciones de concentración, distribución granulométrica y velocidad de caída de las
partículas.
102
CITAS
(heras, n.d.)
(UNESCO, 1971)
(Ramirez, 2003)
103
BIBLIOGRAFIA
Colombia. MAVDT. PNGIRH. (2010). Política Nacional para la Gestión Integral del Recurso
Hídrico. Bogotá, D.C. - Colombia: Nuevas Ediciones Ltda.
Chow Ven Te, Hidrología Aplicada, Editorial Mc Graw Hill 1994.
Delgado, W., Charry, A., & Ortiz, J. (2015, 06 12). Universidad Católica de Colombia. Retrieved
from https://www.ucatolica.edu.co/portal/wp-content/uploads/2015/12/Documento-balance-
hidrico.-U.-Cat%C3%B3lica.-Julio-2015.pdf
Manual de drenaje para carreteras, Ministerio de Transporte, Instituto Nacional de Vías 2009.
IDEAM. (2010). Estudio Nacional del Agua 2010. En M. y. Instituto de Hidrología, Estudio
Nacional del Agua 2010. Bogotá D.C.
104
Reyes Trujillo, A., Ulises Barroso, F., & Carvajal Escobar, Y. (2010). Guía Básica para la
Caracterización Morfométrica de Cuencas Hidrográficas. Cali: Universidad del Valle.
105
Anexo A - INFORMACION SOLICITADA A IDEAM
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
Anexo B - ISOYETAS E ISOTERMAS MENSUALES ELEBORADAS CON LOS DATOS
DE LAS ESTACIONES CERCANAS A LA CUENCA DE LA QUEBRADA NISCOTA
122
Anexo V Mayo Anexo VI Junio
123
Anexo XI Noviembre Anexo XII Diciembre
124
Anexo XV Marzo Anexo XVI Abril
125
Anexo XXI Septiembre Anexo XXII Octubre
126
Anexo C - CURVA DE DURACIÓN DE CAUDALES RIO TOCARÍA, Y QUEBRADA
NISCOTA, ELABORADAS CON CAUDALES MEDIOS MENSUALES PERIODO 2010 -
2015
33 117,1 50
34 94,3
35 81,25 0
0% 20% 40% 60% 80% 100%
36 28,88 Frecuencia
37 12,02
38 8,56
39 11,67 **Nota: El año 2012 se descarta por tener datos insuficientes.
127
Anexo D - CATEGORIAS PARA EL INDICE DE ESCASEZ
Tabla D-1 Categorías para el índice del uso del agua o índice de escasez.
128
Anexo E - GRANULOMETRÍA DEL MATERIAL DEL LECHO DE LA QUEBRADA
NISCOTA D50=42MM.
129