Primera entrega proyecto de hidrología

Análisis cuenca Río Anaime (Tolima)

David Ricardo Villarreal Bermúdez - 25422357
Julián David Torres Bermúdez - 25422337
Lizan Teresa Peñaranda Chía - 25422545
Christian Camilo Pira Roa - 25422327

Presentado a Ingeniero Erasmo Rodríguez

Universidad Nacional de Colombia, Facultad de Ingeniería
Departamento de Ingeniería Civil y Agrícola
Bogotá D.C. 2016

I. Introducción
En el presente reporte se sustentará la primera entrega del proyecto final del curso de hidrología, el
cual tiene como objeto realizar la caracterización hidroclimatológia y la determinación de caudales
de diseño en la cuenca del rio Anaime.
En primer lugar se realizará una breve explicación sobre la obtención de la información
(cartográfica, de precipitación, caudales y temperatura). Posteriormente se mostrará la delimitación
manual de la cuenca hasta la estación hidrométrica Puente Luisa y se describirá el proceso.
En segundo lugar se determinará la ubicación de todas las estaciones hidrométricas usadas para
obtener los datos utilizados, a partir del glosario de la base de datos del IDEAM suministrado en el
curso. Luego a partir de una solicitud a la misma entidad, en la que se obtuvo la información de
interés (precipitación, caudales y temperatura) de las estaciones necesitadas, se seleccionará el
periodo de análisis, y finalmente se realizará un riguroso análisis estadístico de los datos de
precipitación, detección y remoción de puntos anómalos, llenado de datos faltantes y análisis de
consistencia de la información; para así poder realizar los mapas sobre la espacialización mensual
multianual de la precipitación para cada mes y para el año; también se construirán los diagramas de
caja y patilla de la precipitación promedio mensual multianual y caudales promedio en las
estaciones con la información consistente al final del análisis.

II. Metodología
El inicio del proyecto consistió en determinar y obtener las planchas y cartografía base necesarias
para cubrir toda la cuenca del Rio Anaime en escala 1:25000 a partir de la página virtual del IGAC,
esta cartografía debía principalmente cubrir la red de drenaje hacia las 2 estaciones hidrométricas:
Puente Luisa y El Chuzo. Mediante este proceso se determinaron 4 planchas es total [243-IV-A;
243-IV-C; 244-III-B; 244-III-D].
Utilizando el programa AutoCAD se realizó el proceso de delimitación manual de la cuenca. En
primer lugar se ubicó el punto de interés hasta donde se iba a estudiar el tramo del rio el cual
corresponde a la estación Pte Luisa. Posteriormente se iluminó los ríos y drenajes presentes en la
cuenca y sus alrededores para tener un mejor apoyo visual, con cierta dificultad debido a la gran
cantidad de afluentes o baja visibilidad en ciertas zonas. Finalmente se trazó la línea divisoria de
aguas que encierra a la cuenca hidrográfica de interés sabiendo que esta debe cortar las curvas de
nivel perpendicularmente y pasar por los puntos de mayor altura, en las planchas anteriormente
mencionadas. Asi mismo se utilizó AutoCAD para estimar el área de drenaje, obteniendo los
siguientes resultados. [Ver anexo 1].
2
A=75,0265 km
Luego a partir de la información del catálogo del IDEAM suministrada por el curso, se logró
determinar cuáles de las estaciones se encontraban en cercanías a la cuenca; se obtuvo la
información de las estaciones a través de la página virtual, se logró identificar en cuales de estas
estaciones se tenía información y de esta manera se determinaron las estaciones necesarias para
realizar el análisis. Se obtuvieron los registros de precipitación para las siguientes estaciones:

 Estación Hacienda Palogrande

 Estación Las Delicias
 Estación El Plan
 Estación La Cascada
 Estación Hacienda Cucuana
 Figura No. 1 Mapa de las estaciones hidrométricas utilizadas en el análisis

ESTACION LATITUD (GRADOS) LONGITUD (GRADOS)

HDA, PALOGRANDE 4°21'0'' N 75°25'0'' W
LAS DELICIAS 4°22'48'' N 75°30'40,5'' W
EL PLAN 4°21'4,5'' N 75°30'44,9'' W
LA CASCADA 4°17'5'' N 75°32'31'' W
HDA, CUCUANA 4°20'28,9'' N 75°31'6,8'' W
PTE, LUISA 4°20'28,4'' N 75°31'12,4'' W
EL CHUZO 4°24'27,7'' N 75°25'55,4'' W
Tabla No1. Coordenadas estaciones

 Determinación del periodo de análisis común

Al analizar los registros de información sobre precipitación mensual para las 5 estaciones, se logró
determinar que el periodo de análisis común observando la uniformidad en la información será de
enero de 1987 a diciembre de 2015
Estación Periodos de registros
La Cascada 1983 – 2001 ; 2010 - 2015
Hacienda Cucuana 1987 - 2015
El Plan 1983 - 2015
Las Delicias 1984 - 2015
Hacienda Palogrande 1986 - 2015
Tabla No. 1 Periodos de análisis para las estaciones seleccionadas

III. Análisis de datos de precipitación

 Detección y remoción de puntos anómalos en las series

FALTA DESCRIPCION DEL METODO – FALTA ANEXO 2
Se debe recordar que a lo largo de la historia, Colombia se ha visto afectada por una serie de ciclos
del fenómeno “ENSO”, ocasionando aumentos o decrementos significativos de las lluvias durante
su época de influencia. Por lo tanto, fue necesario recurrir a algunas fuentes bibliográficas para
obtener un histórico de periodos en los cuales el fenómeno ENSO afectó la hidrología del país,
siendo esta referencia usada en el análisis de los datos anómalos, en donde, dependiendo del
periodo en el que se presentó el fenómeno, se decidió aceptar o rechazar el dato. Sin embargo, si
algún dato se encontraba en el periodo afectado pero era bastante exagerado, se consideró como
anómalo. Se obtuvieron los resultados mostrados en el Anexo 2.
 Llenado de datos faltantes
Luego de determinar y remover los valores atípicos, se determinó el llenado de los datos a partir de
la utilización del Software “DataFit”, por medio del análisis estadístico. DataFit es un programa de
análisis estadístico, el cual realiza regresiones no lineares a partir de datos que pueden tener de 1 a
20 variables independientes, el programa a partir de los datos realiza iteraciones buscando cual es el
modelo de regresión que mejor se ajusta, obteniendo resultados de regresiones tipo polinómica,
logarítmica, polinomio, polinomio inverso, etc. Posteriormente, el programa arroja el coeficiente de
correlación para la regresión elegida, de acuerdo a lo establecido para el desarrollo del trabajo, un
buen coeficiente R² es que sea mayor a 0,6, de este modo, se elegía un modelo de regresión de
acuerdo al valor del coeficiente de correlación R².
Finalmente, una vez elegido el modelo estadístico que mejor se ajustara a los datos analizados,
DataFit obtiene la ecuación que define el modelo, con las diferentes variables; La ecuación es usada
para obtener los diferentes datos faltantes. En el anexo 3 se pueden apreciar los resultados obtenidos
mediante DataFit, se registraron los datos de los coeficientes de relación para los diferentes casos
analizados, y los resultados finales.
Aquellos valores que evaluándolos por el método anterior, no cumplieron con un coeficiente de
correlación mayor a 0.6; fueron hallados mediante el método IDW (Inverso de la distancia al
cuadrado), que consiste en una relación entre el dato de precipitación faltante en una estación
especifica y las precipitaciones captadas por las otras estaciones teniendo en cuenta las distancias de
separación entre ellas. Las tablas con el llenado de los datos faltantes se pueden hallar en el anexo 3.
Las distancias fueron halladas ubicando las estaciones en “Google Earth” con las coordenadas
geográficas y midiendo las diferentes distancias. Las tablas con las distancias se pueden hallar en el
Anexo 4.

 Análisis de consistencia
El análisis de consistencia se realizó sobre los datos completos, después de haber removido los
datos atípicos y realizado el llenado de los datos faltantes mediante los procedimientos descritos
anteriormente. Se realizaron las curvas de doble masa entre la precipitación acumulada de cada año
durante el periodo de análisis para cada estación de estudio versus la precipitación acumulada de
forma anual para una estación índice en el análisis, la cual para este caso, es el promedio de
precipitación anual de todas las estaciones. Utilizando las herramientas de Excel se realizaron las
curvas de doble masa que se muestran a continuación.
 Curva de doble masa - Estación Hacienda Cucuana
1800.00
1600.00 f(x) = 0.845269453008827 x + 341.74932243702
1400.00 R² = 0.693482315036158
1200.00
Pacum (mm)

1000.00
800.00
600.00
400.00
200.00
0.00
600.00 800.00 1000.00 1200.00 1400.00 1600.00 1800.00
Pacum E.I. (mm)

Grafica No. 1 Curva de doble masa para la estación Hacienda Cucuana

Curva de doble masa - Estación La cascada

1800.00
1600.00 f(x) = 0.56116129717673 x + 494.542444172714
R² = 0.719696261955434
1400.00
1200.00
Pacum (mm)

1000.00
800.00
600.00
400.00
200.00
0.00
600.00 800.00 1000.00 1200.00 1400.00 1600.00 1800.00 2000.00 2200.00

Pacum E.I. (mm)

Grafica No. 2 Curva de doble masa para la estación La Cascada

 Curva de doble masa - Estación El Plan
1800.00
1600.00
f(x) = 0.994852773141444 x + 154.570653043409
1400.00 R² = 0.807702230545724
Pacum (mm)
1200.00
1000.00
800.00
600.00
400.00
200.00
0.00
600.00 700.00 800.00 900.00 1000.00 1100.00 1200.00 1300.00 1400.00 1500.00 1600.00
Pacum E.I. (mm)

Grafica No. 3 Curva de doble masa para la estación El Plan

Curva de doble masa - Estación Las Delicias

1800.00
1600.00
P acum. (mm)

1400.00 f(x) = 0.598402050864432 x + 523.215841711892

R² = 0.60549802268042
1200.00
1000.00
800.00
600.00
400.00
200.00
0.00
200.00 400.00 600.00 800.00 1000.00 1200.00 1400.00 1600.00 1800.00 2000.00
P acum. E.I. (mm)

Grafica No. 4 Curva de doble masa para la estación Las Delicias

Curva de doble masa - Estación Hda. Palogrande

1800.00
1600.00
1400.00 f(x) = 0.621624750434278 x + 252.071782335162
R² = 0.667202790634596
P acum. (mm)

1200.00
1000.00
800.00
600.00
400.00
200.00
0.00
600.00 800.00 1000.00 1200.00 1400.00 1600.00 1800.00 2000.00 2200.00 2400.00
P acum. E.I. (mm)

Grafica No. 5 Curva de doble masa para la estación Hacienda Palogrande

Se observa como los diagramas de doble masa (Precipitación acumulada anual de cada estación vs
precipitación acumulada en la estación índice) no poseen una tendencia lineal, por lo tanto, se deben
aplicar factores de corrección a la información. Este comportamiento no homogéneo se puede deber
a cambios en los aparatos, cambios en los operarios y cambio de cobertura alrededor. Para realizar
la corrección por no homogeneidad se halla la pendiente So en la ecuación
Rantiguo∗So
Rnuevo=
S1
Esta So se multiplica por la precipitación acumulada en la estación índice hallando la precipitación
acumulada corregida (Las tablas con los datos completos se encuentran en el Anexo 5);
posteriormente se halla la diferencia entre la precipitación acumulada corregida y la precipitación
acumulada inicial, dando como resultado los deltas de precipitación, que se reparten
ponderadamente entre los meses de los años correspondientes a cada periodo de análisis, las tablas
de precipitación corregida se encuentran en el Anexo 6. De esta manera se obtuvo el diagrama de
doble masa corregido y consistente, con una pendiente aproximadamente constante y uniforme,
como se observa a continuación:

Curva de doble masa (corregida) - Estación Hacienda Cucuana

2000.00
1800.00
f(x) = 1.097 x
1600.00 R² = 1
1400.00
Pacum (mm)

1200.00
1000.00
800.00
600.00
400.00
200.00
0.00
600.00 800.00 1000.00 1200.00 1400.00 1600.00 1800.00
Pacum E.I. (mm)

Grafica No. 6 Curva de doble masa corregida para la estación Hacienda Cucuana

Curva de doble masa (corregida) - Estación La Cascada

1400.00
1200.00
f(x) = 0.701790000000001 x
Pacum (mm)

1000.00
R² = 1
800.00
600.00
400.00
200.00
0.00
600.00 800.00 1000.00 1200.00 1400.00 1600.00 1800.00
Pacum E.I. (mm)

Grafica No. 7 Curva de doble masa corregida para la estación La Cascada

 Curva de doble masa (corregida) - Estación El Plan
2000.00
1800.00
f(x) = 1.10860000000001 x − 6.59383658785373E-12
1600.00 R² = 1
1400.00
P acum. (mm)

1200.00
1000.00
800.00
600.00
400.00
200.00
0.00
600.00 800.00 1000.00 1200.00 1400.00 1600.00 1800.00
Pacum E.I. (mm)

Grafica No. 8 Curva de doble masa corregida para la estación El Plan

Curva de doble masa (corregida) - Estación Las Delicias

1600.00
1400.00 f(x) = 0.889450000000002 x
1200.00
P acum. (mm)

R² = 1
1000.00
800.00
600.00
400.00
200.00
0.00
600.00 800.00 1000.00 1200.00 1400.00 1600.00 1800.00
P acum. E.I. (mm)

Grafica No. 9 Curva de doble masa corregida para la estación Las Delicias

Curva de doble masa (corregida) - Estación Hda. Pa-

logrande
1500.00
f(x) = 0.838490000000002 x
R² = 1
P acum. (mm)

1000.00

500.00

0.00
600.00 800.00 1000.00 1200.00 1400.00 1600.00 1800.00
P acum. (mm)

Grafica No. 10 Curva de doble masa corregida para la estación Hacienda Palogrande
  Espacializacion precipitación mensual multianual FALTA 12 mapas
 Precipitación areal promedio anual FALTA 1 MAPA
 Diagrama de cajas y patillas precipitación promedio mensual multianual FALTA

IV. Análisis de datos de caudal

Los análisis estadísticos se aplican a los datos obtenidos de las estaciones Puente Luisa y El Chuzo

Figura 2. Mapa de las estaciones limnimétricas utilizadas en el análisis

 Detección y corrección de puntos anómalos FALTA

V. Conclusiones
En cuanto a consistencia de los datos a través de los diagramas de doble masa, se observa que en
general los datos de precipitación no son consistentes, tal como se comentó en el análisis de
consistencia, y por lo tanto fue necesario realizar la corrección por no homogeneidad en todas las
gráficas.

El proceso de delimitación manual de una cuenca hidrográfica puede llegar a ser dispendioso,
motivo por el cual existen varios programas y ayudas electrónicas que facilitan y dan precisión al
momento de trabajar con una cuenca hidrográfica. Sin embargo, el resultado obtenido al delimitar la
cuenca de forma manual puede llegar a ser más claro visualmente que al delimitar la cuenca de
forma automática, ya que la delimitación manual trabaja con mapas reales elaborados por el IGAC,
mientras que la delimitación automática se basa en los datos proporcionados en los “DEM” que son
de modelos digitales del terreno lo cuales no poseen gran precisión de cada detalle y su utilidad
varía mucho con la resolución
VI. BIBLIOGRAFIA
[1] Periodos del ENSO en Colombia
www.ggweather.com/enso/oni
“El niño and La niña Years and Intensites”
[2] Figuras 1 y 2
Tomadas y modificadas de Google Earth

VII. TABLA DE ANEXOS

Anexo 1. Delimitación manual de la cuenca

Anexo 2. Datos anómalos en cada estación

Anexo 3. Llenado de datos faltantes

Anexo 4. Distancias entre estaciones

Anexo 5. Precipitación acumulada corregida

Anexo 6. Precipitación mensual corregida