EVALUACIÓN DE LA AFECTACIÓN DE LOS VERTIMIENTOS DE

AGUA RESIDUAL DE LOS CASCOS URBANOS DE LOS MUNICIPIOS

DE MIRAFLORES Y SAN EDUARDO SOBRE LA FUENTE
RECEPTORA DEL RÍO LENGUPA.

Willie Jonathan Borda Bastidas

Juana Valentina Camacho Sierra
EVALUACIÓN DE LA AFECTACIÓN DE LOS VERTIMIENTOS DE AGUA RESIDUAL DE
LOS CASCOS URBANOS DE LOS MUNICIPIOS DE MIRAFLORES Y SAN EDUARDO
SOBRE LA FUENTE RECEPTORA DEL RÍO LENGUPA.

WILLIE JONATHAN BORDA BASTIDAS

JUANA VALENTINA CAMACHO SIERRA

DOCENTE: Ingeniera CLAUDIA ROCÍO SUÁREZ CASTILLO

MODELAMIENTO Y SIMULACIÓN AMBIENTAL

UNIVERSIDAD SANTO TOMÁS
SECCIONAL TUNJA
2022
CONTENIDO
1. INTRODUCCIÓN ........................................................................................................................ 6
2. OBJETIVOS ................................................................................................................................ 6
2.1 Objetivo General: .............................................................................................................................................. 6
2.2 Objetivos Específicos: ..................................................................................................................................... 6
3. GENERALIDADES ..................................................................................................................... 7
3.1 Características generales de la cuenca Lengupa .......................................................................... 7
3.1.1 Localización ................................................................................................................ 7
3.1.2 Extensión .................................................................................................................... 8
3.1.3 Uso de suelo ............................................................................................................... 9
3.1.4 Morfometría ............................................................................................................. 11
3.2 Generalidades del municipio de Miraflores ................................................................................. 12
3.3 Generalidades del municipio de San Eduardo ............................................................................ 12
4. ANÁLISIS CLIMÁTICO ............................................................................................................ 12
4.1 Análisis y tratamiento de datos climáticos .................................................................................. 12
4.2 Precipitación............................................................................................................................................. 17
4.3 Identificación de periodos húmedos, moderados y secos ..................................................... 18
4.4 Brillo solar ................................................................................................................................................. 20
4.5 Humedad relativa ................................................................................................................................... 20
4.6 Temperatura............................................................................................................................................. 21
4.7 Tratamiento de datos............................................................................................................................ 22
4.8 Climogramas ............................................................................................................................................. 24
5. RESPUESTA HIDRLÓGICA ..................................................................................................... 25
5.1 Respuesta hidrológica estimada ...................................................................................................... 26
5.1.1 Balance hídrico ........................................................................................................ 26
6. CONCLUSIONES....................................................................................................................... 35
7. REFERENCIAS ......................................................................................................................... 36
8. ANEXOS .................................................................................................................................... 37
LISTA DE TABLAS
Tabla 1: Área de la cuenca mediante las herramientas ArcGIS, QGIS y Hec-HMS ....................... 8
Tabla 2: Litología de la cuenca Lengupá ......................................................................................... 9
Tabla 3. Índice Oceánico del Niño (ONI) ....................................................................................... 13
Tabla 4. Estaciones aferentes a la cuenca que contienen datos. ................................................. 16
Tabla 5. Correlación entre la precipitación y caudal anual ......................................................... 19
Tabla 6. Clasificación hidrológica de los suelos ........................................................................... 27
Tabla 7. Clasificación de clases antecedentes de humedad (AMC) por el método de
abstracciones de lluvia del SCS ....................................................................................................... 27
Tabla 8. Números de curva de escorrentía para tipos de uso del suelo determinados en la
cuenca Lengupa ............................................................................................................................... 28
Tabla 9. Factor de corrección para ETP según la latitud ............................................................. 31

LISTA DE FIGURAS
Figura 1. Delimitación de la cuenca................................................................................................. 8
Figura 2: Uso actual de suelo en la cuenca...................................................................................... 9
Figura 3: Formaciones y periodos con las planchas geológicas 191-Tunja y 210-Guateque ... 10
Figura 4. Drenajes de la cuenca Lengupá ..................................................................................... 11
Figura 5. Estaciones pluviométricas, climatológicas y Limnimétricas aferentes a la cuenca... 15
Figura 6. Polígonos de Thiessen con estaciones Pluviométricas y Climatológicas ................... 16
Figura 7. Mapa de brillo solar en Boyacá ...................................................................................... 20
Figura 8. Mapa de humedad relativa............................................................................................. 21

LISTA DE GRAFICAS
Grafica 1. Fenómeno de La Niña, El Niño y Neutral de acuerdo a precipitación y temperatura
........................................................................................................................................................... 14
Grafica 2. Coeficientes pluviométricos relacionados con el régimen de precipitaciones ......... 17
Grafica 3. Precipitación multianual mensual de 1999 a 2016...................................................... 18
Grafica 4. Hidrograma caudales anuales multianuales de la serie de años ............................... 19
Grafica 5. Temperatura mínima y máxima multianual mensual Miraflores 1999 a 2014 ....... 21
Grafica 6. Temperatura mínima y máxima multianual mensual Rondón 1999 a 2014 ............ 22
Grafica 7. Completado de datos faltantes por el método de regresión lineal para las dos
estaciones aferentes ........................................................................................................................ 23
Grafica 8. Curvas de doble masa, estación Miraflores y Rondón ................................................ 24
Grafica 9. Climogramas generados para la estación de Miraflores, periodo húmedo, seco y
moderado. ........................................................................................................................................ 25
Grafica 10. Curva de caudales e Hidrograma 1999 periodo húmedo ........................................ 33
Grafica 11. Curva de caudales e Hidrograma 2002 periodo seco ............................................... 34
Grafica 12. Curva de caudales e Hidrograma 2003 periodo moderado ..................................... 34
LISTA DE ECUACIONES
Ecuación 1. Ecuación para hallar el incremento conforme a las curvas de duración de
caudales y la precipitación estimada por medio de regresión lineal. .......................................... 22
Ecuación 2. Fórmula general para balance hídrico ..................................................................... 26
Ecuación 3. Fórmula explícita de balance hídrico ....................................................................... 26
Ecuación 4. Ecuación para hallar CN II ......................................................................................... 27
Ecuación 5. Ecuaciones para hallar CN I y CN III ......................................................................... 28
Ecuación 6. Ecuación para el cálculo de Pe y S............................................................................. 29
Ecuación 7. Ecuación para el cálculo de la Infiltración................................................................ 29
Ecuación 8. Ecuación para hallar el índice de calor mensual y un índice de calor anual ......... 30
Ecuación 9. Ecuación para hallar la constante de ajuste ............................................................. 30
Ecuación 10. Ecuación para hallar la evapotranspiración .......................................................... 30
Ecuación 11. Ecuación para hallar la reserva de humedad final ................................................ 31
Ecuación 12. Ecuación para hallar percolación ........................................................................... 31
Ecuación 13. Ecuación para hallar la escorrentía superficial ..................................................... 32
Ecuación 14. Ecuación para hallar la escorrentía total ............................................................... 32
Ecuación 15. Ecuación para hallar la probabilidad de excedencia ............................................. 32

LISTA DE ANEXOS
8.1 Anexo 1. Cierre de cuenca por QGIS .................................................................................. 37
8.2 Anexo 2. Cierre de cuenca por Hec-HMS........................................................................... 37
8.3 Anexo 3. Polígonos de Thiessen con estaciones pluviométricas y climatológicas ........ 39
8.4 Anexo 4. Mapa de pendientes cuenca Lengupá ................................................................ 40
 1. INTRODUCCIÓN

Los ríos constituyen un sistema de circulación de fluidos y sedimentos importantes a

través de las cuencas hidrográficas con el fin de ayudar a la conservación de los
ecosistemas y diferentes formas de vida, así pues, también ayudan con los procesos de
regulación sistemática del entorno como la captación y transporte de vertimientos,
gases y vapores para su tratamiento y purificación (Díaz & Torres, 2000).
Por otro lado, la presente investigación se basa en el análisis del comportamiento
hidrológico de la cuenca frente a temporadas húmedas, secas y moderadas mediante
los balances hídricos con las metodologías Soil Conservation Service (SCS) y ABCD
Thomas para poder compararlos con la realidad teniendo en cuenta la estación
limnimétrica (Chapasia). Esto, para que a través de un modelo, se logre analizar las
condiciones generadas en el agua, luego de realizar los vertimientos de los municipios
Miraflores y San Eduardo y permita representar los fenómenos físicos, químicos y
biológicos que se muestran en él (Valcárcel, y otros, 2010). Teniendo en cuenta lo
anterior, se emplea la modelación hidrológica ya que es una herramienta que ayuda a
realizar el análisis y la prevención de eventos extremos como las inundaciones, sin
embargo, en este caso se emplea con el fin de plantear pronósticos realistas para la
toma de decisiones referentes a la posible ubicación de la planta de tratamiento de
aguas residuales de los municipios y verificar la capacidad de autodepuración de la
fuente hídrica (Estrada & Pacheco, 2012).

Finalmente, a lo largo del documento se encuentran las generalidades de los municipios

de estudio y algunas características morfométricas y geológicas de la cuenca Lengupa,
seguido de esto, se encuentra el uso actual del suelo teniendo en cuenta la delimitación
de la cuenca. Así pues, se encuentra el análisis climático donde se menciona la relación
con el brillo solar, la temperatura, precipitación y la evapotranspiración que se presenta
para posteriormente aplicar el balance hídrico y comparar estos datos con la estación
limnimétrica.

2. OBJETIVOS

2.1 Objetivo General:

Evaluar la afectación de los vertimientos de agua residual de los cascos urbanos de los
municipios de Miraflores y San Eduardo sobre la fuente receptora del río Lengupa.
2.2 Objetivos Específicos:
 Estimar la respuesta hidrológica de la cuenca Lengupa generada durante tres
períodos: húmedos, secos y moderados influenciados por precipitaciones altas,
medias y bajas.
  Determinar el movimiento, comportamiento y dinámica del recurso hídrico
sobre la fuente receptora a partir de un modelo hidráulico para la ubicación de
las dos Plantas de Tratamiento de Aguas Residuales.
 Definir el estado de la fuente hídrica receptora de vertimientos de agua residual
mediante el uso de un modelo de calidad que garantice el análisis de las cargas
contaminantes.
 Proponer las acciones necesarias mitigando así los efectos ocasionados por los
vertimientos sobre la fuente receptora.

3. GENERALIDADES
3.1 Características generales de la cuenca Lengupa
3.1.1 Localización
La cuenca de Lengupá, pertenece al área hidrográfica del Orinoco, y así mismo, a la zona
hidrográfica del meta. Teniendo en cuenta lo anterior, el río Lengupá es uno de los más
importantes dentro del departamento de Boyacá, el cual nace a partir de la confluencia
de los ríos Fuche y Mueche en el municipio de Zetaquira para posteriormente
desembocar en el río Upía en el departamento de Casanare y finalmente unirse con el
río Meta (Corpoboyacá, 2020).
Así pues, esta cuenca se encentra en la región Andina y presenta cerca de 12 tipos de
coberturas siendo predominante el bosque denso alto de tierra firme con
aproximadamente el 24% del paisaje (Corpoboyacá, 2012). Cabe resaltar que la cuenca
Lengupá en su totalidad comprende hasta el municipio de Campohermoso, sin
embargo, con el fin de tener información real en cuanto al comportamiento de los
caudales del río Lengupá, ya que esta es la fuente receptora de los vertimientos de
Miraflores y San Eduardo, se identificó la estación limnimétrica más cercana, en este
caso la ubicada en el municipio de Miraflores, en la cual se realizó el cierre de la cuenca
mediante las herramientas, ArcGis, QGis y Hec-HMS con el fin de verificar los órdenes
de magnitud y la veracidad de los datos relacionados con el área de la cuenca. De esta
manera, el cierre por medio de la primera aplicación se muestra en la Figura 1, así
como el área de influencia de los municipios que la conforman.
 Figura 1. Delimitación de la cuenca.

Fuente: Elaboración propia autores

3.1.2 Extensión
Teniendo en cuenta lo anterior, en la Tabla 1 se muestran las áreas que se obtuvieron
en el cierre de la cuenca por medio de las tres herramientas, mostrando que no se
presentan variaciones considerables. Por otra parte, en el Anexo 1 y Anexo 2 se
muestra el cierre de la cuenca por QGis y Hec-HMS.
Tabla 1: Área de la cuenca mediante las herramientas ArcGIS, QGIS y Hec-HMS

Herramienta Área (Km2)

ArcGIS 915,78
QGIS 915,92
Hec-HMS 915,77
Fuente: Elaboración propia autores
Así mismo, presenta un perímetro aproximado de 153,9 Km, e incluye los municipios
de Rondón, Pesca, Tota, Zetaquirá, Aquitania, Ramiriquí, Miraflores, Berbeo y San
Eduardo del departamento de Boyacá.
3.1.3 Uso de suelo
Las rocas presentes en la cuenca, corresponde a unidades sedimentarias relativamente
consolidadas, así mismo, se desarrollan geoformas propiamente estructurales
denudacionales, y depósitos cuaternarios de tipo coluvial y aluvial, por esta razón, se
pueden encontrar escarpes elevados o topografías abruptas (Corpoboyacá, 2012). Así
pues, en la Figura 2 se muestran los usos de suelo actual de la cuenca.
Figura 2: Uso actual de suelo en la cuenca.

Fuente: Elaboración propia autores

Teniendo en cuenta lo anterior, los usos de suelo de la cuenca se basan en actividades
9.1
del campo, como lo son tipos de ganaderías, cultivos y algunos bosques y el paisaje que
predomina en la cuenca son montañas estructurales erosionales. De esta manera, su
litología varía a lo largo de la cuenca, como se muestra en la Tabla 2.
Tabla 2: Litología de la cuenca Lengupá

LITOLOGÍA DE LA CUENCA LENGUPÁ

Depósitos superficiales clásticos hidrogravigénicos
 Depósitos superficiales clásticos mixtos
Depósitos superficiales piroclásticos de ceniza volcánica
sobre rocas sedimentarias clásticas limo arcillosas
Rocas sedimentarias clásticas arenosas con intercalaciones
de limo arcillosas
Rocas sedimentarias clásticas limo arcillosas con
intercalaciones de arenosas
Rocas sedimentarias clásticas limo arcillosas y depósitos
superficiales clásticos gravigénicos mixtos
Rocas sedimentarias clásticas mixtas y depósitos
superficiales piroclásticos de ceniza volcánica
Adaptado de: Litología de suelos en Boyacá, IGAC.
Así pues, se observa que sobresalen los depósitos 10.1 superficiales piroclásticos sobre

rocas sedimentarias clásticas limo arcillosas, y en segundo ámbito se presentan rocas

sedimentarias clásticas arenosas. De igual manera, como se muestra en la Figura 3, se
emplearon las planchas 191-Tunja y 210-Guateque del servicio geológico colombiano,
con el fin de identificar las formaciones presentes y a qué periodo histórico pertenece,
de esta forma, se encontró que predominan las formaciones del periodo cretáceo
especialmente la formación UNE, Conejo, Chipaque grupo churuvita, Guaduas, entre
otras.
Figura 3: Formaciones y periodos con las planchas geológicas 191-Tunja y 210-
Guateque

Fuente: Elaboración propia autores

Por lo anterior, se define que la parte norte de la cuenca va a emplear un grupo
hidrológico B y la parte sur un grupo hidrológico C, debido a su composición.
 3.1.4 Morfometría
Por otro lado, algunas características morfométricas que hacen referencia la corriente
principal del río Lengupá son que presenta una pendiente promedio de 5,8 % al inicio
de la cuenca, sin embargo, a partir del paso por el municipio de Rondón, su pendiente
varía al rededor 1.23%. Así mismo, en el Anexo 4 se muestra el mapa de pendientes
que se presentan en la cuenca.
Además, sus cotas máximas y mínimas son:
 Cota superior: 3845
 Cota inferior: 1070
La longitud del cauce principal es de 70,85Km, de esta manera tambien se encontraron
diversos drenajes entre quebradas y ríos, los cuales se muestran en la Figura 4.
Figura 4. Drenajes de la cuenca Lengupá

Fuente: Elaboración propia autores

 3.2 Generalidades del municipio de Miraflores

El municipio de Miraflores se encuentra en el departamento de Boyacá a 1500 m.s.n.m,

ubicado en la región Andina con coordenadas 73° 8' 47'' W 5° 11' 49'' N. Limita por el
oriente con Berbeo por el sur con Campohermoso, por el suroriente con Páez, por el
norte con Zetaquira y por el occidente con Chinavita y Garagoa. Así mismo, se encuentra
en la provincia de Lengupá, de la cual es la capital (Alcaldía Municipal de Miraflores,
2022).
3.3 Generalidades del municipio de San Eduardo

El municipio de San Eduardo se encuentra en el departamento de Boyacá a una altura

aproximada de 1705 m.s.n.m, con coordenadas 73° 4' 41'' W 5° 13' 24'' N. Por otra parte,
limita por el norte con Aquitania, por el sur con Zetaquira, por el oriente con Aquitania
y por el occidente con Berbeo. Finalmente, al igual que Miraflores, hace parte de la
provincia de Lengupá (Alcaldía Municipal de San Eduardo, 2022).

4. ANÁLISIS CLIMÁTICO

El análisis climático del presente informe se centra en recolectar y procesar

información acerca del comportamiento de las variables climáticas, como lo son
precipitación, temperatura, brillo solar y humedad relativa, las cuales son obtenidas de
plataformas abiertas destacándose la del IDEAM. Por lo anterior, para hacer la descarga
de los datos meteorológicos se estableció un periodo de 22 años, comprendido desde el
año 1999 al 2021, por consiguiente, de acuerdo al tipo de análisis se requiere que los
valores sean descargados de forma diaria para así obtener una distribución
meteorológica mensual, anual y multianual, cuyo propósito sea organizar la
información de forma visual recurriendo a graficas que permitan entender la relación
que tienen dichas variables climáticas con el comportamiento y respuesta hídrica de la
cuenca Lengupá.
4.1 Análisis y tratamiento de datos climáticos
4.1.1 Fenómenos macro climáticos
La identificación de la variabilidad climática extrema, se refiere a las fluctuaciones
observadas en el clima durante periodos de tiempo relativamente cortos, es decir, es
temporal y transitorio (Montealegre, 2009), además, es importante definir los
fenómenos macro climáticos como la Oscilación del Sur (ENSO), refiriéndose al
fenómeno del Niño o de la Niña los cuales tienen incidencia sobre una extensión grande
de terreno afectando a todo un país y en ciertas ocasiones coincide con los mismos
periodos climáticos establecidos dentro de una cuenca. El ENSO depende de la
 temperatura, precipitación y de los vientos ocasionados en la franja tropical del océano
pacifico.
El fenómeno del Niño en Colombia hace referencia al calentamiento inusual de las aguas
del pacifico ecuatorial, generando una disminución de las precipitaciones y aumento de
la temperatura en el territorio colombiano, a diferencia del fenómeno de la Niña que
produce un aumento sobre lo normal de la lluvia y para la Región Andina se presenta
en los meses de octubre-noviembre y es menos frecuente que el fenómeno del Niño
(Valiente, 1999).
Por otro parte los índices, como el Índice Oceánico del Niño (ONI), se encarga de
cuantificar a través de la media de tres meses consecutivos las anomalías en la
temperatura superficial del mar, esta es valorada en un umbral de +/- 0,5°C, donde si el
valor es consecutivo por más de 5 veces y es mayor o igual a 0,5 significa que es un
periodo de sequía asociada al fenómeno del Niño, por el contrario, si es menor o igual a
-0,5°C indica fenómeno de la Niña, y finalmente, si se encuentra entre -0,5°C < X < 0,5°C,
quiere decir que épocas moderadas. De este modo, la relaciona los valores de los años
para la serie de 1999 a 2021, según los valores determinados por el ONI el fenómeno
de La Niña se presentó en 1999, 2000, 2010, 2011 y 2021, para la fase neutral se
destacan los años 2001, 2003 y 2013 y para el evento de tipo El Niño están los años
2002, 2004 y 2015.
Tabla 3 relaciona los valores de los años para la serie de 1999 a 2021, según los valores
determinados por el ONI el fenómeno de La Niña se presentó en 1999, 2000, 2010, 2011
y 2021, para la fase neutral se destacan los años 2001, 2003 y 2013 y para el evento de
tipo El Niño están los años 2002, 2004 y 2015.
Tabla 3. Índice Oceánico del Niño (ONI)
Year DJF JFM FMA MAM AMJ MJJ JJA JAS ASO SON OND NDJ
1999 -1.5 -1.3 -1.1 -1.0 -1.0 -1.0 -1.1 -1.1 -1.2 -1.3 -1.5 -1.7
2000 -1.7 -1.4 -1.1 -0.8 -0.7 -0.6 -0.6 -0.5 -0.5 -0.6 -0.7 -0.7
2001 -0.7 -0.5 -0.4 -0.3 -0.3 -0.1 -0.1 -0.1 -0.2 -0.3 -0.3 -0.3
2002 -0.1 0.0 0.1 0.2 0.4 0.7 0.8 0.9 1.0 1.2 1.3 1.1
2003 0.9 0.6 0.4 0.0 -0.3 -0.2 0.1 0.2 0.3 0.3 0.4 0.4
2004 0.4 0.3 0.2 0.2 0.2 0.3 0.5 0.6 0.7 0.7 0.7 0.7
2005 0.6 0.6 0.4 0.4 0.3 0.1 -0.1 -0.1 -0.1 -0.3 -0.6 -0.8
2006 -0.9 -0.8 -0.6 -0.4 -0.1 0.0 0.1 0.3 0.5 0.8 0.9 0.9
2007 0.7 0.2 -0.1 -0.3 -0.4 -0.5 -0.6 -0.8 -1.1 -1.3 -1.5 -1.6
2008 -1.6 -1.5 -1.3 -1.0 -0.8 -0.6 -0.4 -0.2 -0.2 -0.4 -0.6 -0.7
2009 -0.8 -0.8 -0.6 -0.3 0.0 0.3 0.5 0.6 0.7 1.0 1.4 1.6
2010 1.5 1.2 0.8 0.4 0.2 -0.7 -1.0 -1.3 -1.6 -1.6 -1.6 -1.6
2011 -1.4 -1.2 -0.9 -0.7 -0.6 -0.4 -0.5 -0.6 -0.8 -1.0 -1.1 -1.0
2012 -0.9 -0.7 -0.6 -0.5 -0.3 0.0 0.2 0.4 0.4 0.3 0.1 -0.2
2013 -0.4 -0.4 -0.3 -0.3 -0.4 -0.4 -0.4 -0.3 -0.3 -0.2 -0.2 -0.3
2014 -0.4 -0.5 -0.3 0.0 0.2 0.2 0.0 0.1 0.2 0.5 0.6 0.7
2015 0.5 0.5 0.5 0.7 0.9 1.2 1.5 1.9 2.2 2.4 2.6 2.6
2016 2.5 2.1 1.6 0.9 0.4 -0.1 -0.4 -0.5 -0.6 -0.7 -0.7 -0.6
2017 -0.3 -0.2 0.1 0.2 0.3 0.3 0.1 -0.1 -0.4 -0.7 -0.8 -1.0
2018 -0.9 -0.9 -0.7 -0.5 -0.2 0.0 0.1 0.2 0.5 0.8 0.9 0.8
2019 0.7 0.7 0.7 0.7 0.5 0.5 0.3 0.1 0.2 0.3 0.5 0.5
2020 0.5 0.5 0.4 0.2 -0.1 -0.3 -0.4 -0.6 -0.9 -1.2 -1.3 -1.2
2021 -1.0 -0.9 -0.8 -0.7 -0.5 -0.4 -0.4 -0.5 -0.7 -0.8 -1.0 -1.0
Adaptado de: Índice Oceánico del Niño (ONI)
De la misma forma, utilizando la plataforma Climate Data Library, se puede establecer
con mayor precisión las variaciones que ocurren en el municipio de Miraflores
(municipio de referencia) al correlacionar la precipitación y temperatura con los
fenómenos de La Niña, El Niño o Neutral (moderado), evidenciado en la Grafica 1 para
los años de 1999 a 2016.
Grafica 1. Fenómeno de La Niña, El Niño y Neutral de acuerdo a precipitación y
temperatura

Fuente: Climate Data Library ([Link]

A partir de la Grafica 1 para la serie de años 1999-2016, se evidencia que los periodos
húmedos son 1999, 2000, 2008 y 2011 identificados por medio de las barras de color
azul, para los secos están los años 2003, 2005, 2010, 2015 y 2016 (barras de color rojo),
los demás son periodos moderados reconocidos por las barras de color gris.
4.1.2 Identificación de estaciones susceptibles a análisis
Para cuantificar la incidencia que tienen las variables climáticas de la cuenca y así
determinar acertadamente la respuesta hídrica de la cuenca Lengupá, se parte por
identificar las estaciones limnimétricas, pluviométricas y climatológicas que son
aferentes a la cuenca en mención. A continuación, en la Figura 5 se visualizan todas las
estaciones aferentes a la cuenca sin discriminar las que tienen ausencia de datos
completos para la serie de tiempo establecida.
Figura 5. Estaciones pluviométricas, climatológicas y Limnimétricas aferentes a la
cuenca

Fuente: Elaboración propia autores

Así mismo, en la Tabla 4 se pueden observar las estaciones que contienen datos
completos, luego de identificar estas estaciones se procede al descargue y análisis de la
información para lograr determinar los 9 años objeto de estudio para las temporadas
húmedas, secas y moderadas.
 Tabla 4. Estaciones aferentes a la cuenca que contienen datos.

Municipio Código Tipo de Coordenadas decimales Altura

estación estación Latitud Longitud (msnm)
Rodón 35085020 PM y CL 5,35841667 -73,2036111 2120
Zetaquira 35080010 PM 5,28297222 -73,1692778 1720
Berbeo 35080030 PM 5,18372222 -73,0867222 1370
Miraflores 35085040 PM y CL 5,19255556 -73,1447778 1640
Miraflores 35087060 LM 5,16411111 -73,0915278 1210
Fuente: Elaboración propia autores
A partir de la Tabla 4, se realiza la división de la cuenca mediante la metodología de
Polígonos de Thiessen con el fin de determinar el área de influencia de cada estación,
de esta manera, se obtuvo la Figura 6 en la cual se evidencian dichas áreas con relación
a las estaciones tanto climatológicas como pluviométricas de la cuenca.
Figura 6. Polígonos de Thiessen con estaciones Pluviométricas y Climatológicas

Fuente: Elaboración propia autores

Sin embargo, por diversos factores de la asignatura, en este caso se analizarán los datos
obtenidos a partir de dos estaciones tanto pluviométricas (PM) como climatológicas
(CL) y una limnimétrica (LM), con el fin de generar los polígonos de Thiessen, los cuales
se pueden observar en el Anexo 3. Dichas estaciones se encuentran en el municipio de
Rondón y de Miraflores, así mimo, los datos que contienen estas estaciones son de
calidad por presentar un mínimo de valores faltantes y de esta manera poder realizar
el tratamiento de datos respectivo, cuya finalidad sea definir los años húmedos,
moderados y secos presentes en la cuenca Lengupá.
4.2 Precipitación
La precipitación es el producto de la condensación del vapor de agua atmosférico que
se deposita como líquido o sólido en la superficie del suelo o de los océanos (IDEAM,
2005). De igual manera, según la estacionalidad de la precipitación y dependiendo de
cuantas veces se presenten periodos secos y húmedos, el régimen puede ser: 1)
monomodal, caracterizado por tener un periodo seco y uno húmedo muy marcado
durante el año en ciertas épocas de este o 2) bimodal donde se registran dos
temporadas secas alternadas con dos lluviosas (Guzmán, Ruíz, & Cadena, 2014).
Una forma de determinar la estacionalidad presente en un lugar es por medio de la
realización de coeficientes pluviométricos generados a partir de la mensual anual
multianual, los cuales si el valor presentado en un mes es menor a 1,00 significa que
hace parte de un periodo seco, sin embargo, si es superior a 1,00 se caracteriza por ser
un periodo húmedo, debido a lo mencionado anteriormente, la Grafica 2 representa los
coeficientes pluviométricos obtenidos para las dos estaciones pluviométricas aferentes
a la cuenca de Lengupá.
Grafica 2. Coeficientes pluviométricos relacionados con el régimen de precipitaciones

Fuente: Elaboración propia autores

Según el estudio de Guzmán, Ruíz y Cadenas en 2014, establecieron que el tipo de
régimen que presenta el occidente de Boyacá es Bimodal 3 (Bm3), teniendo dos fases
secas y dos de mayores lluvias, identificando que la primera temporada seca es menos
prolongada y las épocas de mayores lluvias son similares (Guzmán, Ruíz, & Cadena,
2014). Sin embargo, al visualizar la Grafica 3, en la estación del municipio de Miraflores
se aprecia un régimen monomodal conforme a que el primer trimestre del año y los dos
últimos meses son caracterizados por bajas precipitaciones y para el municipio de
Rondón es bimodal porque se diferencian dos escenarios donde el coeficiente
pluviométrico es superior a 1,00 durante los meses de abril a agosto y de octubre a
noviembre.
De igual modo, para observar el comportamiento de la precipitación media mensual-
multianual para los años 1999 al 2016, se recurrió nuevamente a la plataforma Climate
Data Library, lo cual permite identificar un comportamiento monomodal en el
municipio de Miraflores, debido a que solo se presenta un pico en las precipitaciones,
contrario a ello, en el municipio de Rondón se aprecia dos picos significando un régimen
bimodal.
Grafica 3. Precipitación multianual mensual de 1999 a 2016

Fuente: Climate Data Library ([Link]

4.3 Identificación de periodos húmedos, moderados y secos
Para establecer de forma asertiva los tres periodos en cuestión, se hizo necesario
recurrir a tres variables que son los fenómenos macro climáticos (ENSO), los valores de
caudales medios anuales y la precipitación media anual, con ello se hace el debido
análisis y contraste de la información permitiendo definir de forma precisa los nueve
años. A continuación, se relacionan la Grafica 4 y Tabla 5.
 Grafica 4. Hidrograma caudales anuales multianuales de la serie de años

Fuente: Elaboración propia autores

La Grafica 4 presenta el comportamiento de los caudales a medida del tiempo,
permitiendo observar valores altos, moderados, bajos y sin valor de acuerdo a los
reportados por la estación limnimétrica Chapasia, también, es de resaltar que para la
mayoría de años de la serie se obtuvo una cantidad significativa de datos faltantes a
excepción de nueve años, esto son 1999 al 2004, 2007, 2017 y 2018, destacando que
para 1999 se presenta el mayor valor de caudal y el menor se generó en 2001 conforme
a lo registrado en la estación y por el proceso de cálculos realizados para obtener los
caudales medios anuales.
Por otra parte, correlacionando las precipitaciones anuales totales para los nueve años
mencionados anteriormente se detalla la Tabla 5, donde de acuerdo al caudal alto o
bajo se organiza cada año en húmedo, moderado y seco.
Tabla 5. Correlación entre la precipitación y caudal anual
Estación Periodos
pluviométrica Húmedo Moderado Seco
1999 2004 2007 2000 2003 2018 2001 2002 2017
Q Chapasia 60 57 56 48 44 49 35 39 43
(m3/s)
Miraflores 1583,8 1686,4 1082,2 1731,4 1878,1 1639,6 1497 1564,6 1679,7
P (mm)

Rondón 1929,2 2101,7 1773,9 2113,4 2076 2095,3 1926,9 2529,4 2261,6
Suma 3513 3788,1 2856,1 3844,8 3954,1 3734,9 3423,9 4094 3941,3

Fuente: Elaboración propia autores

De esta forma, al comparar los valores sumados entre las dos estaciones pluviométricas
de Miraflores y Rondón con los caudales medio anuales, no siempre coincide la relación
de que si uno es alto (caudal) el otro también (precipitación), aunque lo adecuado sería
 que fuese así, conforme a ello, la variable principal para establecer las épocas climáticas
es mediante el valor de los caudales medios anuales.
4.4 Brillo solar
El brillo solar es el número de horas de sol efectivo que se relaciona con la cantidad de
tiempo durante el cual la superficie de la cuenca es irradiada por el sol directo. De esta
manera, con ayuda del atlas interactivo del IDEAM se realizó la descarga del mapa
nacional multianual de brillo solar en el territorio, dando como resultado la Figura 7.
Figura 7. Mapa de brillo solar en Boyacá

5-6
horas
4-5
horas
3-4
horas

Adaptado de: Atlas de Radiación Solar IDEAM. 2014

Teniendo en cuenta lo anterior, y como se muestra en la Figura 7 la cuenca está
influenciada por un brillo solar que comprende de 4 a 5 horas. Por otra parte, se obtuvo
una gráfica de brillo solar en el municipio de Tunja, con el fin de contrastarla con la
incidencia de son en la cuenca Lengupa debido a su cercanía. De esta manera, se
encontró que los meses en los que aumenta el brillo solar son enero, febrero y
diciembre.
4.5 Humedad relativa
Se entiende como el porcentaje de humedad que contiene el aire, respecto al total de lo
que es capaz de contener en función de la temperatura y presión. Esta va
incrementando a medida que el aire se enfría, o cuando se agrega vapor de agua al aire
(IDEAM, 2005).
 Figura 8. Mapa de humedad relativa

500-1000mm
1000-1500mm
1500-2000mm
2000-2500mm
2500-3000mm

Fuente: Atlas de Humedad Relativa IDEAM. 2014

Así pues, al igual que se hizo con el brillo solar, se tomó la gráfica del municipio de Tunja,
en este caso si es posible observar algunas diferencias, esto puede deberse a la
temperatura ambiente o la concentración de lluvias en la cuenca.
4.6 Temperatura
La temperatura ambiente de la cuenca del río Lengupá fluctúa entre los 11 y 22°C, como
se puede observar en las Grafica 5 y Grafica 6, mostrando la temperatura multianual
mensual de los dos lugares (Miraflores y Rondón) donde se encuentran las estaciones
climatológicas de aferencia con la cuenca, destacando que la temperatura media
multianual mensual es alrededor de los 16,5°C.
Grafica 5. Temperatura mínima y máxima multianual mensual Miraflores 1999 a
2014

Tomado: Climate Data Library ([Link]

 Grafica 6. Temperatura mínima y máxima multianual mensual Rondón 1999 a 2014

Tomado: Climate Data Library ([Link]

4.7 Tratamiento de datos
El tratamiento de datos se realiza con el objetivo de: 1) comparar la consistencia de los
datos por diferentes métodos como por ejemplo la curva de doble masa, la cual es una
representación de la altura de precipitación acumulada a través del tiempo, desde el
inicio de una tormenta o serie de años hasta su terminación y 2) completar los datos
faltantes mediante algún método, para este caso, se emplea la regresión lineal que
consiste en el cálculo de los valores inexistente a través de la ecuación de la recta de
regresión lineal, derivada para este caso mediante la distribución o dispersión formada
entre dos estaciones.
Ecuación 1. Ecuación para hallar el incremento conforme a las curvas de duración de
caudales y la precipitación estimada por medio de regresión lineal.
△𝑃
𝐼𝑛𝑐 = ; 𝑦 = 𝑎 + 𝑏𝑥
△𝑡
Donde:
Inc = Incremento
△ 𝑃 = Incremento de precipitación en el intervalo (mm)
△ 𝑡 = Incremento en el tiempo
y = Precipitación estimada
x = Precipitación de la estación que contiene datos completos
a y b = son constantes
4.7.1 Completado de datos faltantes por regresión lineal
Los datos faltantes o inexistentes entre las dos estaciones son de aproximadamente un
total de 205 datos para los nueve años escogidos, estos datos inexistentes están
distribuidos en 5 años (1999, 2003, 2004, 2007 y 2018); La Grafica 7 se realiza para el
año 2007, esto debido a que es el año con mayor cantidad de valores inexistentes (92),
repartidos en 31 para la estación de Rondón y 61 para Miraflores, razón por la cual, se
hizo necesario utilizar las dos estaciones más cercanas a estas que son la estación de
Zetaquira (Rondón) y Berbeo (Miraflores), así mismo, la Grafica 7 muestra la
dispersión generada con la estación completa e incompleta y sus respectivas
ecuaciones para la línea de tendencia de tipo lineal.
Grafica 7. Completado de datos faltantes por el método de regresión lineal para las
dos estaciones aferentes

Fuente: Elaboración propia autores

Con base en la Grafica 7, se obtiene la ecuación necesaria en cada uno de los casos para
completar los datos faltantes, donde x se reemplaza por el valor de la precipitación
diaria de la estación con valores completos. El procedimiento se realizó para cada uno
de los años generando las gráficas correspondientes.
4.7.2 Curvas de doble masa
Las curvas de doble masa son obtenidas a partir del cálculo de precipitaciones
acumuladas en las dos estaciones pluviométricas de Miraflores y Rondón, así mimo, se
escogen los nueve años establecidos con anterioridad, para establecer el
comportamiento y las posibles fluctuaciones generadas y así lograr al realizar una
comparación entre las dos gráficas obtenidas. Razón de lo anterior, se muestran las
siguientes graficas (Ver Grafica 8).
 Grafica 8. Curvas de doble masa, estación Miraflores y Rondón

Fuente: Elaboración propia autores

4.8 Climogramas
Los climogramas o diagramas ombrotermicos son representaciones anuales de
temperatura y precipitación media, son importantes para determinar la influencia que
tienen estas dos variables en relación a periodos secos y húmedos. Por consiguiente, la
Grafica 9, tiene en cuenta los años 1999 (húmedo),2002 (seco) y 2003 (moderado), los
cuales son los que generan menos sesgos para ser analizados, de igual forma se
presentan los climogramas de la estación de Miraflores debido a que estas tienen
aferencia en los dos municipios de estudio (Miraflores y San Eduardo). Al visualizar la
Grafica 9 se establece que el pico más alto de precipitación se presenta en el periodo
seco con un valor superior a 300 mm para el mes de agosto, por el contrario, el periodo
húmedo arroga una precipitación inferior a los 250 mm.
 Grafica 9. Climogramas generados para la estación de Miraflores, periodo húmedo,
seco y moderado.

Fuente: Elaboración propia autores

A nivel general el comportamiento de los tres periodos se puede inferir que la

precipitación y la temperatura son variables inversamente proporcionales, donde si
aumenta una, disminuya la otra, es decir, si la precipitación media mensual es elevada,
la temperatura media mensual va a disminuir en una proporción aproximada de 3 a 4
°C.

5. RESPUESTA HIDRLÓGICA

De acuerdo a los valores obtenidos en el análisis climático se pudo determinar cuáles

son los nueve años a tener en cuenta para realizar la cuantificación de las diferentes
variables que intervienen para obtener las curvas de duración de caudales por medio
de la estación limnimétrica (Chapasia), pluviométricas y climatológicas identificadas
con anterioridad, conducen al modelo del sistema hidrológico general (balance hídrico)
y al método ABCD Thomas, cuya finalidad sea comparar cuál de las dos metodologías
empleadas se aproxima a los valores observados (reales) la estación limnimétrica.
Conocer la respuesta hídrica de la cuenca es fundamental para predecir futuros
comportamientos que se puede generar en la cuenca y así poder tomar decisiones
acertadas relacionadas con las intervenciones a realizar, como por ejemplo la ubicación
de Plantas de Tratamiento de Aguas Residuales.
 5.1 Respuesta hidrológica estimada
5.1.1 Balance hídrico
Consiste en una serie de pasos articulados que permite cuantificar y analizar
parámetros morfológicos y climáticos de una zona de estudio, cuya finalidad radica en
suministrar información valiosa acerca del estado y comportamiento del recurso
hídrico a través de los años. El método más empleado para realizar el balance es el de
lluvia-escorrentía, el cual se basa en utilizar un valor de Número de Curva - CN
ponderado o CNII (valor establecido según las condiciones del suelo), este método es
considerado de tipo indirecto.
[Link] Ecuación de balance hídrico
Se basa en establecer las variables que conectan las entradas (ganancias) y salidas
(perdidas) dentro de un sistema, en este caso una cuenca hidrológica, con la ayuda de
ecuaciones se establece el principio de conservación de masas haciendo referencia a la
ecuación de continuidad, por lo anterior la ecuación general se basa en la sumatoria de
entradas menos salidas para así poder establecer ese cambio de almacenamiento que
ocurre durante un periodo de tiempo y sobre un área en particular.

Ecuación 2. Fórmula general para balance hídrico

∑ 𝐸𝑛𝑡𝑟𝑎𝑑𝑎𝑠 − ∑ 𝑆𝑎𝑙𝑖𝑑𝑎𝑠 =△ 𝐴𝑙𝑚𝑎𝑐𝑒𝑛𝑎𝑚𝑖𝑒𝑛𝑡𝑜

Por lo anterior la formula general de balance es:

Ecuación 3. Fórmula explícita de balance hídrico

𝑃 − 𝑅𝑜 − 𝐼 − 𝐸𝑉𝑃 =△ 𝐴𝑙𝑚𝑎𝑐𝑒𝑛𝑎𝑚𝑖𝑒𝑛𝑡𝑜
𝑃 = 𝑅𝑜 + 𝐼 + 𝐸𝑉𝑃
Donde:
P = Precipitación
Ro= Escorrentía superficial
I = Infiltración
EVP = Evapotranspiración
[Link] Calculo del Número de Curva – CN
Este es un número adimensional que va desde 0 a 100, entre mayor sea el CN significa
que hay superficies impermeables permitiendo un porcentaje de escorrentía alto. Este
tipo de CN se establece conforme a las características del suelo como el tipo, uso y
formaciones del suelo, asignando el grupo hidrológico (ver Tabla 6) y el valor
correspondiente, además de la pendiente del terreno contribuye a ese porcentaje de
escurrimiento, como se menciona a partir de la Figura 2.
Tabla 6. Clasificación hidrológica de los suelos
Grupo hidrológico Descripción

Grupo A Bajo potencial de Altas tasa de infiltración: arena profunda, suelos

escorrentía profundos depositados por el viento, limos agregados.

Grupo B Moderadamente bajo Tasa de infiltración moderada: suelos poco profundos

potencial de depositados por el viento, marga arenosa.
escorrentía

Grupo C Moderadamente alto Infiltración lenta: marga arcillosa, margas arenosas

potencial de poco profundas, suelo con bajo contenido orgánico y
escorrentía con altos contenidos de arcilla.

Grupo D Alto potencial de Baja tasa de infiltración: suelos que se expanden

escorrentía significativamente cuando se mojan, arcillas altamente
plásticas y ciertos suelos salinos.

Fuente: (Chow, Maidment, & Mays, 1994)

Así mismo, existen tres tipos de CN que son clasificado de acuerdo a las condiciones
antecedentes de humedad (AMC), como se muestra en la Tabla 7.
Tabla 7. Clasificación de clases antecedentes de humedad (AMC) por el método de
abstracciones de lluvia del SCS
Condición de humedad Precipitación acumulada de los 5 días previos al
antecedente (AMC) evento en consideración (cm)

Condición seca – AMC I 0 – 3,5

Condición normal – AMC II 3,5 – 5,25

Condición húmeda – AMC III Más de 5,25

Fuente: Soil Conservation Service, 1972, tabla 4.2 (Chow, Maidment, & Mays, 1994)
Para realizar el cálculo, en primera instancia se utiliza la ecuación para hallar el CN II o
ponderado, y se detalla a continuación:
Ecuación 4. Ecuación para hallar CN II
(𝐶𝑁1 ∗ 𝐴1 ) + (𝐶𝑁2 ∗ 𝐴2 ) + (𝐶𝑁𝑛 + 𝐴𝑛 )
𝐶𝑁 𝐼𝐼 =
𝐴𝑟𝑒𝑎 𝑡𝑜𝑡𝑎𝑙
Los números de curva equivalentes pueden calcularse por:
 Ecuación 5. Ecuaciones para hallar CN I y CN III
𝐶𝑁 𝐼𝐼 𝐶𝑁 𝐼𝐼
𝐶𝑁 𝐼 = ; 𝐶𝑁 𝐼𝐼𝐼 =
2,3 − 0,012 𝐶𝑁 𝐼𝐼 0,43 + 0,0057 𝐶𝑁 𝐼𝐼
Donde:
𝐶𝑁1, 𝐶𝑁2 hasta 𝐶𝑁𝑛 = Son los valores asignados de acuerdo a los componentes del suelo
𝐴1 , 𝐴2 hasta 𝐴𝑛 = So las áreas generadas de acuerdo al 𝐶𝑁1, 𝐶𝑁2 o 𝐶𝑁𝑛
𝐶𝑁 𝐼 = Número de curva para una condición seca
𝐶𝑁 𝐼𝐼 = Número de curva para una condición normal
𝐶𝑁 𝐼𝐼𝐼 = Número de curva para una condición húmeda
De acuerdo a las características del suelo explicadas con anterioridad se determinan
dos tipos de CN II uno con valor de 64,66 (Estación Rondón) y 76,23 (Estación
Miraflores), debido a que hay dos estaciones pluviométricas que tienen aferencia con la
cuenca, los cuales son relacionados de acuerdo a la Tabla 8.
Tabla 8. Números de curva de escorrentía para tipos de uso del suelo determinados
en la cuenca Lengupa
Thiessen Estación Tipo de suelo CN II Área Total
(m2)
Área (m2) Grupo Condición CN
hidrológico
Thiessen 1 Rondón Ganadería extensiva 64,66 521219387
493670 B Regular 69

217155694 B Regular 69

60826439 B Buena 67

Cultivo de papa, maíz, cebada

54714169 B Buena 78

Bosque natural

188029415 B Buena 55

Thiessen 2 Miraflores Cultivo de arveja, hortalizas 76,23 254128368

28625496 C Buena 85

Ganadería extensiva

988598 C Regular 73

65708091 C Regular 79
 82864895 C Buena 74

Bosque natural

75941288 C Regular 73

El área de estudio corresponde a dos tipos uno grupo hidrológico B, moderadamente

bajo potencial de escorrentía y el otro C, moderadamente alto potencial de escorrentía,
esto se definió por la presencia de arenas y arcillas en la zona.
[Link] Precipitación efectiva método SCS para abstracciones
El Soil Conservation Service en 1972 desarrollo un método para calcular las pérdidas
de la precipitación de una tormenta, la precipitación efectiva es la diferencia entre la
precipitación total y la infiltración, correspondiendo a una cantidad de lluvia que
produce escorrentía superficial, está dada por la siguiente ecuación:
Ecuación 6. Ecuación para el cálculo de Pe y S
(𝑃 − 0,2𝑆)2 1000
𝑃𝑒 = ;𝑆 = − 10
𝑃 + 0,8𝑆 𝐶𝑁
Donde:
Pe = Precipitación efectiva en milímetros
P = Precipitación en milímetros
S = Retención potencial máxima en pulgadas o milímetros
CN = Numero de Curva
El valor de la precipitación efectiva siempre será menor o igual a la profundidad de
precipitación y tiene que ser calculado de acuerdo al valor de humedad antecedente
(AMC I, II o III) de los 5 días de precipitación acumulada.
[Link] Infiltración
Proceso por el cual el agua se transfiere desde la superficie del terreno hacia las
profundidades, muchos factores influyen en la tasa de infiltración como la morfología,
composición del terreno y su cubierta vegetal (Chow, Maidment, & Mays, 1994). Por
tanto, la ecuación es la diferencia entre la precipitación y la precipitación efectiva dada
en milímetros.
Ecuación 7. Ecuación para el cálculo de la Infiltración

𝐼 = 𝑃 − 𝑃𝑒
 [Link] Calculo de evapotranspiración, estimación con método de
Thornthwaite
La evapotranspiración es la combinación de evaporación (proceso por el cual el agua
líquida pasa al estado de vapor en condiciones naturales) desde la superficie del suelo
y la transpiración de la vegetación y su valor está dominado por el suministro de
energía (luz solar) y el transporte de vapor (Ven Te Chow & Mays, 1994).
El método de Thornthwaite, es el más utilizado para el cálculo de evapotranspiración
(ETP) de un lugar para cada mes del año a partir de parámetros básicos como la
temperatura media mensual y el número teórico de horas del sol establecido por el
rango de latitud (norte-sur). Mediante la utilización de las siguientes ecuaciones se
obtiene el valor de la ETP potencial media para cada mes del año.
Ecuación de índice de calor mensual y un índice de calor anual (sumatoria de los índices
mensuales)
Ecuación 8. Ecuación para hallar el índice de calor mensual y un índice de calor anual
12
𝑡 1,514
𝑖=( ) ; 𝐼 = ∑𝑖
5
𝑖=1

Se procede a calcular la ETP mensual sin corregir, mediante la siguiente ecuación,

estableciendo el coeficiente experimental de ajuste y la ETP sin corregir, es decir, sin
ninguna calibración que varía en función de la latitud.
Ecuación 9. Ecuación para hallar la constante de ajuste

𝑎 = 0,675𝑥10−6 ∗ 𝐼 3 − 0,771𝑥10−4 ∗ 𝐼 2 + 1,792𝑥10−2 ∗ 𝐼 + 0,49239

Ecuación 10. Ecuación para hallar la evapotranspiración
10 ∗ 𝑡 𝑎
𝐸𝑇𝑃 = ( ) ∗ 1,6
𝐼
Donde:
ETP = Evapotranspiración potencial media mensual (mm/mes)
i = índice de calor mensual
I = índice de calor anual
t = Temperatura media mensual (°C)
a = Coeficiente experimental de ajuste
Por último, se procede hacer la calibración de la ETP que consiste en multiplicar el valor
obtenido de ETP por un factor de correlación en este caso y para donde están ubicadas
las estaciones corresponde a una latitud 5°N (Colombia) para cada año expuestos en la
Tabla 9.
 Tabla 9. Factor de corrección para ETP según la latitud

Lat Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
5°N 1,02 0,93 1,03 1,02 1,06 1,03 1,06 1,05 1,01 1,03 0,99 1,02
Fuente: (Método de Thornthwaite)
De esta forma se halla el valor de la ETP para cada mes del año.
[Link] Reserva de humedad del suelo
La humedad del suelo corresponde al volumen de agua que se infiltra y antes de
alcanzar el nivel freático es removido por las raíces de plantas o evaporación. Por tanto,
existen dos tipos de humedad del suelo, una inicial y otra final, las cuales son
determinadas por las características de uso de suelo y vegetación presente, influyendo
en que se presente un menor o mayor porcentaje de inflación y ETP.
La reserva de humedad inicial (RHi), se determina con base a las características del
suelo o a estudios detallados que indican un valor que varía entre 0 a 100, por ende, al
no encontrar el valor real del suelo se define que sea 100 para iniciar y para continuar
con el día siguiente siempre constituye el valor obtenido del día o mes anterior
establecido por la reserva de humedad final.
La reserva de humedad final (RHf), es calculada a partir de la influencia de tres
parámetros los cuales son RHi, infiltración y ETP, donde las dos primeras aportan en la
humedad del suelo evidenciado por medio a la siguiente ecuación.
Ecuación 11. Ecuación para hallar la reserva de humedad final
𝑅𝐻𝑓 = 𝑅𝐻𝑖 + 𝐼 − 𝐸𝑇𝑃

Dependiendo del resultado se siguen tres condiciones para obtener un valor acertado,
donde si el resultado de RHf supera los 100 se establece como valor verdadero 100
(RHi+I-ETP>=100 =100) o por el contrario si el valor es negativo inferior a cero se
establece como valor real cero (RHi+I-ETP<0=0), si no se cumplen las dos condiciones
anteriores y el resultado esta entre 0 y 100 se deja el valor correspondiente.
[Link] Percolación
Hace referencia al proceso mediante el cual el agua atraviesa el suelo, está determinado
por la infiltración, las abstracciones o perdidas por diferencias de humedades (RHi y
RHf) y la ETP, representado conforme a la siguiente ecuación.
Ecuación 12. Ecuación para hallar percolación

𝑃𝑒𝑟𝑐 = 𝐼 − (𝑅𝐻𝑖 − 𝑅𝐻𝑓 ) − 𝐸𝑇𝑃

De igual forma si el resultado obtenido es inferior a cero, se establece como valor real
cero, y de no ser así se deja el valor positivo superior a cero.
 [Link] Escorrentía superficial y total
La escorrentía superficial es el proceso por el cual el agua fluye por la superficie del
terreno hacia los cauces, esto debido a una abundancia de precipitación, lo que ocasiona
un exceso de infiltración o de saturación del suelo. Esta es determinada de acuerdo a
los antecedentes de escorrentía del mes o día anterior y percolación, por lo cual, se
puede recurrir a dos ecuaciones, según el caso de que se tenga o no el valor antecedente
de escorrentía:
Ecuación 13. Ecuación para hallar la escorrentía superficial

𝐸𝑠𝑐 𝑆𝑢𝑝 = 𝑃𝑒𝑟𝑐 ∗ (1 − 0,99)

𝐸𝑠𝑐 𝑆𝑢𝑝 = 𝐸𝑠𝑐 𝑑𝑖𝑎 𝑜 𝑚𝑒𝑠 𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟 ∗ 0,99 + 𝑃𝑒𝑟𝑐 ∗ (1 − 0,99)
La escorrentía total hace referencia al resultado generado por la suma entre la
precipitación efectiva y la escorrentía superficial significando el aporte a un caudal.
Ecuación 14. Ecuación para hallar la escorrentía total

𝐸𝑠𝑐 𝑡𝑜𝑡𝑎𝑙 = 𝑃𝑒 + 𝐸𝑠𝑐 𝑆𝑢𝑝

[Link] Caudal
Finalmente, para poder obtener las curvas de duración de caudales es necesario
convertir la escorrentía total (mm/día) o unidad de volumen segundo (m3/s), por lo
cual se hace la conversión de unidades según corresponda y de este modo se multiplica
por el valor de área establecido que es aferente a cada estación (Miraflores y Rondón).
[Link] Curvas de duración de caudales
La curva de duración es un procedimiento grafico para el análisis de los datos de
caudales y representan la frecuencia acumulada de ocurrencias de un caudal (diario,
semanal, mensual o anual) determinado. Este método es comúnmente utilizado para
validar modelos hidrológicos o para comparar y modelar la respuesta hidrológica de
una cuenca. La ecuación que representa esa probabilidad de ocurrencia es la siguiente:
Ecuación 15. Ecuación para hallar la probabilidad de excedencia
𝑚
𝑃 = ( ) ∗ 100
𝑁
Donde:
P = Probabilidad de excedencia u ocurrencia
m = Número de veces que se presente en ese tiempo el caudal
N = Número de datos
De manera análoga, el hidrograma es la representación gráfica de las variables del
caudal de una corriente con respecto al tiempo. Este es usado debido a que permite
contribuir en el análisis obtenido entre la estación limnimétrica (Chapasia) y el
balance hidrológico realizado, para así estimar la relevancia del método de lluvia-
escorrentía.
[Link] Análisis de curvas de duración de caudales

Teniendo en cuenta lo mencionado anteriormente, y tomando los resultados obtenidos

a partir del balance hídrico para los años 1999, 2004 y 2007, se realiza la comparación
con la estación limnimétrica para el año 1999, debido a que el comportamiento del
caudal medio-multianual es el mayor entre los periodos húmedos.
Grafica 10. Curva de caudales e Hidrograma 1999 periodo húmedo

En la Grafica 10 se logra observar una variación notable entre lo obtenido con el

balance y la estación limnimétrica, debido a que no se presentaron valores superiores
a los 100 m3/s en comparación con la estación limnimétrica que registra valores
superiores a los 800 m3/s. Sin embargo, al ser un periodo húmedo, la mayoría de datos
con caudales altos deberían encontrarse en un porcentaje de probabilidad menor al
10% para la curva de balance, no obstante, estos tienden a estar ubicados por encima
del 50% de probabilidad. Por otra parte, en el hidrograma se aprecia de mejor manera
la incertidumbre y poca precisión que tiene el balance para este año evidenciando
tendencias distorsionadas en la distribución entre las dos curvas.

De acuerdo a los tres años (2001, 2002 y 2017) seleccionados para los periodos secos,
el 2002 contiene datos completos para precipitación y además cuenta con la
particularidad de ser un año consecutivo, logrando representar una adecuada precisión
en la secuencia de datos trabajados.
 Grafica 11. Curva de caudales e Hidrograma 2002 periodo seco

Respecto a la Grafica 11 se puede observar que se presenta una intersección entre las
dos gráficas cuando la probabilidad es igual al 20%, así mismo, la diferencia de caudales
máximos mensuales es de aproximadamente 100 m3/s de acuerdo a las curvas
generadas, logrando suponer que el balance para este año coincide parcialmente con el
generado a partir de la estación limnimétrica, Sin embargo, es necesario realizar una
calibración más rigurosa en función al CN asignado. En el hidrograma se visualiza una
sobrevaloración en los caudales mensuales al inicio del año, a pesar de ello, se observa
una mínima correlación acorde al aumento y disminución de los caudales para las dos
líneas de tendencia. Es decir, si está aumentando los valores de caudales en la estación
limnimétrica, los valores de caudales del balance también aumentan.
Finalmente, para los periodos moderados se analizaron los años 2000, 2003 y 2018, de
los cuales se identificó para el 2003 un mejor comportamiento gráfico. Además de ser
el año en donde el balance hídrico se acercó a los datos reales.
Grafica 12. Curva de caudales e Hidrograma 2003 periodo moderado

Con base a la Grafica 12 se puede inferir que la mayoría de valores se encuentran

situados previos al 10% de la probabilidad, de igual manera, la diferencia entre lo
observado y calculado es de aproximadamente 55 m3/s con lo cual se evidencia una
disminución en el rango de variación con respecto a los demás años analizados. Por
consiguiente, se establece que es el año con mayor correlación y precisión con respecto
al comportamiento obtenido a partir de la estación limnimétrica. Igualmente, el
hidrograma representa un acercamiento en el valor de los caudales para los meses de
abril a octubre. Por otra parte, aunque imitó el comportamiento de la estación
limnimétrica, no fue posible alcanzar en ningún mes el valor máximo registrado por la
estación.

6. CONCLUSIONES
A manera de conclusión, los tres cierres generados a la altura de la estación
limnimétrica (chapasia) conforma la representación gráfica y extensión de la cuenca
aproximadamente de 915,78 Km2, obtenido conforme a la utilización de tres distintas
herramientas de procesamiento, mediante los cuales se estableció un margen de error
de +/- 0.14 Km2, corroborando la superficie total de la cuenca.
Por otra parte, se encuentra que existe una moderada diferencia en los fenómenos
meteorológicos que se presentan en la cuenca, con respecto a los fenómenos macro
climáticos como el ENSO y el ONI, debiéndose a diversos factores principalmente a las
características meteorológicas de la zona de estudio, además de la presencia de
diversos microclimas a nivel local.
Así mismo, mediante el análisis visual del balance hídrico se encontró que al momento
de comparar con los valores de la estación limnimétrica (Chapasia) hay un mayor
acercamiento a la realidad en las épocas de sequías y moderadas, esto puede deberse al
valor asignado del CN II, puesto que refleja las condiciones edafológicas de la cuenca, y
por ende va a afectar ciertas variables implícitas en el cálculo del balance hídrico, como
lo son la tasa de infiltración, el porcentaje de humedad y la escorrentía. Además, debido
a que los años seleccionados no son consecutivos, no determinan adecuadamente los
valores que relacionan factores acumulativos principalmente asociado con las
características intrínsecas del suelo, presentando una subvaloración en el aporte de
caudal por escorrentía total a la cuenca.
Posteriormente, las curvas de duración de caudales son las que realmente evalúan el
acercamiento del método indirecto de lluvia-escorrentía con la realidad medida por la
estación limnimétrica, puesto que son el resultado obtenido tras la ejecución de varios
cálculos cuya calibración, validación, precisión y comprobación son importante en la
gestión integral del recurso hídrico, también con respeto a las curvas de duración se
hace indispensable establecer el aporte de caudal en diferentes épocas, por ejemplo el
aporte de caudal más alto para el año de 1999, 2002 y 2003 por estación limnimétrica
fue de 816, 246 y 191 m3/s durante un solo día de esos años y el aporte por balance
hídrico fue 87, 189 y 143 m3/s respectivamente, de modo que el modelo de sistema
hidrológico general es demasiado objetivo porque solo constituye las tasas de flujo de
entradas y de salidas, sin tener en cuenta variables de estado o una adecuada
calibración de parámetros para la evaluación hidrológica.
Finalmente, las variaciones extremas de temperatura influyen en la
evapotranspiración, debido a que existen diversos pisos térmicos dentro de la cuenca
desde la cota superior a 3845 hasta 1070 m.s.n.m, así pues, esta temperatura puede
llegar a influir en la duración e intensidad de las precipitaciones a lo largo de años
susceptible a estudio.

7. REFERENCIAS

Alcaldía Municipal de Miraflores. (15 de 02 de 2022). Municipal de Miraflores. Obtenido de

Municipal de Miraflores: [Link]

Alcaldía Municipal de San Eduardo. (14 de 02 de 2022). Municipio de San Eduardo. Obtenido de
Municipio de San Eduardo: [Link]

Chow, V. T., Maidment, D., & Mays, L. (1994). Agua superficial. En V. T. Chow, D. Maidment, & L.
Mays, Hidrología aplicada (pág. 156). McGRAW-HILL.

Corpoboyacá. (2012). Plan de ordenamiento y manejo de la cuenca hidrgráfica del río Upía. Fase
de prospectiva.

Corpoboyacá. (2020). Información previa para el establecimiento de la meta de carga global

contaminante para el segundo quinquenio de la cuenca del río Lengupá. Boyacá:
Subdirección de Ecosistemas y Gestión Ambiental.

Díaz, H. C., & Torres, J. G. (2000). Importancia de los ríos en el entorno ambiental. Revista Del
Instituto De investigación De La Facultad De Minas, Metalurgia Y Ciencias geográficas, 57-
63.

Estrada, V., & Pacheco, R. (2012). Modelación hidrológica con HEC-HMS en cuencas montañosas
de. Ingeniería hidráulica y ambiental, 94-105.

Guzmán, Ruíz, & Cadena. (2014). Regionalización de Colombia según la estacionalidad de la

precipitación media mensual, a través del análisis de componentes principales (ACP).
Bogotá: IDEAM.

IDEAM. (2005). Atlas Climatológico de Colombia. Bogotá: IDEAM.

Montealegre, J. (2009). Informe final de las actividades desarrolladas en el marco del contrato
prestación de servicios No IDEAM 022 - 2009. Bogotá: IDEAM.

Valcárcel, L., Borroto, J., Alberro, N., Griffith, J., Derivet, M., Flores, P., . . . Domínguez, J. (2010).
Modelación de la calidad del agua en el segmento. Scielo.

Valiente, Ó. (1999). Evolución en el estudio del fenómeno ENSO (El Niño-Oscilante del Sur): de
anomalía «local» a la predicción de variaciones climáticas globales. Investigaciones
Geográficas, 5-20.
Ven Te Chow, M. D., & Mays, L. (1994). Agua atmosfeérica. En M. D. Ven Te Chow, & L. Mays,
Hidrología aplicada (pág. 93). McGRAW-HILL.

8. ANEXOS
8.1 Anexo 1. Cierre de cuenca por QGIS

8.2 Anexo 2. Cierre de cuenca por Hec-HMS.

8.3 Anexo 3. Polígonos de Thiessen con estaciones pluviométricas y
climatológicas
 8.4 Anexo 4. Mapa de pendientes cuenca Lengupá

Powered by TCPDF ([Link])

 Índice de comentarios

9.1 poner porcentajes, esta ambiguo el análisis

10.1 no se, no entiendo

Powered by TCPDF ([Link])