Estimación de Caudales Maximos para La Microcuenca Del Municipio de San Francisco de Sales Cundinamarca
Estimación de Caudales Maximos para La Microcuenca Del Municipio de San Francisco de Sales Cundinamarca
Estimación de Caudales Maximos para La Microcuenca Del Municipio de San Francisco de Sales Cundinamarca
1
DETERMINACIÓN DE LOS CAUDAL MÁXIMO Y CARACTERIZACIÓN
MORFOMÉTRICA DE LA MICROCUENCA DEL CAUCE EN EL SECTOR
HIDRODEPENDIENTE DEL MUNICIPIO DE SAN FRANCISCO - CUNDINAMARCA
2
TABLA DE CONTENIDO
INTRODUCCIÓN ............................................................................................................ 9
1. GENERALIDADES ................................................................................................ 10
3
3.1 METODOLOGÍAS ........................................................................................... 19
3.1.1 Limitaciones .................................................................................................. 31
3.2 INSTRUMENTOS O HERRAMIENTAS UTILIZADAS ............................... 37
3.3 ALCANCES Y LIMITACIONES .................................................................... 37
5. CONCLUSIONES Y RECOMENDACIÓNES.................................................... 67
BIBLIOGRAFÍA ............................................................................................................. 68
APÉNDICE ...................................................................................................................... 69
ANEXOS .......................................................................................................................... 69
4
LISTA DE IMÁGENES
5
LISTA DE TABLAS
6
LISTA DE GRAFICAS
7
RESUMEN
El Agua, es prioridad para todos los que habitan este planeta, por esta razón se han constituido en
todo el mundo organizaciones que velan por la gestión adecuada del recurso hídrico, sin embargo,
la gestión por parte de las organizaciones gubernamentales encargadas es insuficiente, generando
escasez y desinformación hidrológica en gran parte de las cuencas del país, este es el caso de San
Francisco Cundinamarca, debido a ello la alcaldía no cuenta con un estudio hidrológico de su cauce
principal, por tanto, el objetivo principal de este documento es presentar la modelación hidrológica
de la microcuenca que abastece al municipio, mediante el uso del software HEC-HMS así como
las herramientas de los Sistema de Información Geográfica (SIG) e informes hidrológicos de una
estación del IDEAM -Instituto de Hidrología, Meteorología y Estudios Ambientales-, además del
modelo de elevación digital de Vertex, portal de datos de la Instalación de Satélites de Alaska para
imágenes de la Tierra con sensores remoto, con el fin de estimar los caudales máximos para los
periodos de retorno de 2, 10, 25, 50 y 100 años bajo la modelación del HEC-HMS con los métodos
de los Hidrograma unitarios de scs, snyder y Clark, el método de estimación por módulo de
rendimiento que acoge la función de distribución log Pearson III y el método racional.
ABSTRAC
The water is priority for all people from this planet, for this reason have been constituted in around
the world organizations that keep vigil to the proper management of water resources, however, the
management for part of government organizations is insufficient, that generates scarcity and
desinformation hydrological in large part of the basins from the country, that is the case of San
Francisco – Cundinamarca, because of that the mayor´s office haven’t an hydrological study of its
main channel, for that, the principal objective of this document is to present the hydrological model
of the micro basin that caters to the people by using the software HEC-HMS, as well as the tools
of geographic information system (GIS) and hydrological reports from the IDEAM station -
hydrology institute, meteorology and studies environmental institute- also the elevations modeling
digital of Vertex -data portal of the Alaska Satellite Installation for Earth images with remote
sensors – in order to estimate the max flows for the return periods of 2, 10, 25, 50 and 100 years
using the modeling of HEC-HMS with the methods of unity hydrographs SCS, SNYDER and
CLARK, the method of estimate for performance module that receive the function of distribution
Log Pearson III and the rational method.
8
INTRODUCCIÓN
Es un hecho, que las cuencas hidrográficas constituyen la base que sostiene la estabilidad
social actual, aportando un recurso considerado no renovable y del cual depende la
continuidad de la vida en el planeta, tal como se conoce. De allí la importancia del cuidado
del recurso hídrico por lo que surge la necesidad de conocer a mayor profundidad parámetros
como las características, el comportamiento y cuidados que a la comunidad tanto le ha
costado y que no ha sabido apreciar particularmente a la población colombiana ocasionando
pérdidas considerables a nivel económico, productivo y social.
Una de las poblaciones afectadas por esta problemática se encuentra al norte del
departamento de Cundinamarca en el municipio de San Francisco donde se evidencia escasez
del recurso agua debido posiblemente a su exceso de habitantes en temporada turística y a
una insuficiencia en la oferta del cuerpo de agua que abastece al municipio, lo que lleva al
caso en cuestión, en el proyecto de grado se efectuó una metodología que tuvo contemplado
la obtención de un modelo digital de elevación (DEM) de resolución de 12.5 metros, el cual
se usó para la delimitación de la cuenca hidrológica usando el software de Sistema de
información geográfica ARCGIS, con esta delimitación se seleccionó la estación más cercana
al área de influencia. A si mismo se desarrolló el procesamiento de los datos obtenidos de la
estación Climatológica Sabaneta operada por el Instituto de Hidrología, Meteorología y
Estudios Ambientales IDEAM, para rastrear datos atípicos y obtener datos de precipitación
máxima de 24 horas.
Además de ello se realizó la identificación del tipo de suelo y cobertura vegetal que compone
el área en estudio, con el fin de obtener la curva número, crear las curvas IDF sintéticas para
diferentes tiempos de retorno, elaborar el hietograma por bloque alterno y el hietograma de
precipitación neta e infiltración, permitiendo la elaboración de los hidrogramas unitarios con
ayuda del Software de modelado Hidrológica HEC-HMS y por último la estimación del
cálculo del caudal máximo mediante diferentes métodos.
9
1. GENERALIDADES
¿Cuáles son los caudales máximos del río cañas hasta el municipio San Francisco
Cundinamarca?
1.3 JUSTIFICACIÓN
La planeación adecuada del recurso hídrico debe ser considerada prioridad para el país así
como para las instituciones educativas, desafortunadamente a lo largo de los años las
actividades antrópicas han afectado este recurso debido a esto las consecuencias de estas
prácticas y a la planeación y manejo inadecuado han sido más evidentes durante los últimos
tiempos, en donde zonas en las que las comunidades contaban con disposición del recurso en
este momento presentan problemas de desabastecimiento o carencia de infraestructura
adecuada para tal fin.
Así como el desarrollo del presente trabajo de grado nos ofrece herramientas y experiencia
en dicha modelación del mismo modo permite desarrollar competencias para futuros
proyectos en la vida laboral.
10
1.4 OBJETIVOS
Estimar los caudales máximos de la subcuenca del río Cañas aplicando la herramienta HEC-
HMS
1.5 CRONOGRAMA
FASE ACTIVIDAD
1 Recolección de información
2 Procesamiento de datos
3 Modelación de la cuenca
4 Elaboración del documento final
5 Entrega del documento final
MAY
ABR
ENE
FEB
DIC
FASES/
Recolección de información
Procesamiento de datos
Modelación de la cuenca
Elaboración del documento final
Entrega del documento final
11
2. MARCOS DE REFERENCIA
Las características físicas de una cuenca tienen una relación estrecha con el comportamiento
de los caudales que transitan por ella. Los parámetros Morfométricos de una cuenca integran
un conjunto de estimaciones realizadas, en la mayoría de los casos, al iniciar un estudio
hidrológico, con fines de aprovechamiento o control. Según (Gaspari a, y otros, 2012). El
análisis morfométrico es el estudio de un conjunto de variables lineales, de superficie, de
relieve y drenaje; que permite conocer las características físicas de una cuenca, lo cual
permite realizar comparaciones entre varias cuencas, así como ayuda a la interpretación de
la funcionalidad hidrológica y en la definición de las estrategias para la formulación de su
manejo.
12
La obtención de las curvas IDF debe realizarse con información pluviográfica de estaciones
ubicadas en la localidad, derivando las curvas de frecuencia correspondientes mediante
análisis puntuales de frecuencia de eventos extremos máximos. La distribución de
probabilidad de Gumbel se recomienda para estos análisis, aunque otras también pueden ser
ajustadas. Eventualmente, es posible hacer análisis regionales de frecuencia en caso de
disponer de más de una estación pluviográfica. Si no existe información en la población, debe
recurrirse a estaciones localizadas en la zona lo más cercanas a la población. Si esto no
permite derivar curvas IDF aceptables para el proyecto, deben ajustarse curvas IDF por
métodos sintéticos, preferencialmente derivados con información pluviográfica colombiana.
De acuerdo con el nivel de complejidad del sistema, la manera mínima permitida de
obtención de las curvas IDF se define en la siguiente Tabla 3.
13
Imagen 1. Componentes básicos del flujo en el hidrograma de escorrentía
“Está diseñado para simular los procesos hidrológicos completos de los sistemas de cuencas
dendríticas. El software incluye muchos procedimientos tradicionales de análisis
hidrológico, como la infiltración de eventos, los hidrogramas unitarios y el enrutamiento
hidrológico. HEC-HMS también incluye los procedimientos necesarios para la simulación
continua, incluida la evapotranspiración, el deshielo y la contabilidad de la humedad del
suelo. También se proporcionan capacidades avanzadas para la simulación de escorrentía
en cuadrícula utilizando la transformación de escorrentía cuasi-distribuida lineal
(ModClark). Se proporcionan herramientas de análisis complementarias para la
optimización del modelo, el flujo de flujo de predicción, la reducción del área de
profundidad, la evaluación de la incertidumbre del modelo, la erosión y el transporte de
sedimentos, y la calidad del agua.
El software presenta un entorno de trabajo completamente integrado que incluye una base
de datos, utilidades de entrada de datos, motor de computación y herramientas de informe
de resultados. Una interfaz gráfica de usuario permite al usuario un movimiento sin
interrupciones entre las diferentes partes del software. Los resultados de simulación se
almacenan en HEC-DSS (sistema de almacenamiento de datos) y se pueden usar junto con
otros programas para estudios de disponibilidad de agua, drenaje urbano, pronóstico de
flujo, impacto de urbanización futura, diseño de aliviaderos de embalses, reducción de daños
por inundación, regulación de llanuras de inundación y operación de sistemas.” (US ARMY
CORPS OF ENGINEERS, s.f.)
14
2.2 MARCO GEOGRÁFICO
Este municipio cuenta con 10 veredas (Arrayán, El Peñón, Juan de Vera, La Laja, Muña,
Pueblo Viejo, Sabaneta, San Antonio, San Miguel y Tóriba.) Ver Imagen 3
15
Imagen 3. Mapa político del municipio de San francisco, Cundinamarca. (Alcaldía
de San Francisco de Sales - Cundinamarca, 2013)
La cabecera del municipio de San Francisco cuenta con 2851 habitantes, de acuerdo a los
resultados consolidados en el Censo del 2005, realizado por el DANE (DANE, 2005, pág.
43), la población total en este municipio es 8187 habitantes de los cuales el 65% de la
población vive en la zona rural y el 35% corresponde a la población de la cabecera municipal.
Así como el 50.4% corresponde a hombres y el 49.6% a mujeres.
16
2.4 ESTADO DEL ARTE
Con el fin de conocer las investigaciones que han sido desarrolladas acerca de la estimación
de caudales máximos y caracterización mormetrica de cuencas en un ámbito nacional e
internacional, se realizó una búsqueda en diferentes bases de datos y en la web con el fin de
recopilar la suficiente información teórica y metodológica planteada para abordar este tema.
A continuación, se sintetizan los aportes más importantes de algunas investigaciones que se
encaminan al tema principal de este trabajo.
Boris Eduardo Carías Juárez, Egly Tatiana Chacón Novoa y Miguel Ángel Márquez, en su
trabajo de grado: Validación de metodologías para el cálculo de caudales máximos en El
Salvador (cita), calculan los tiempos de concentración validando la metodología de estos y
calculan caudales máximos por las metodologías hidrometeorológicas de Formula Racional,
Soil Conversation Service, Hidrogramas Unitarios Complejos e Hidrogramas sintéticos,
realizando a su vez un estudio estadístico de los caudales máximos determinando la función
de distribución que mayor se ajuste, comparan los datos con metodologías puntuales y
regionales y validan las metodologías determinando los factores de ajuste correspondiente.
Luego de la obtención de resultados utilizadas las metodologías propuestas para determinar
tiempos de concentración (como la Kirpich, Agencia Federal de Aviación, Bureau of
Reclamation y California), concluyeron que es posible utilizar cualquiera de las metodologías
utilizadas, pues los resultados obtenidos se encuentran dentro de rangos preestablecidos para
las regiones hidrográficas estudiadas (1-3 m/s). Para la determinación de caudales máximos,
el Método Racional subestima los valores de caudales máximos para las áreas analizadas en
la investigación, la metodología del hidrograma Sintético de Sneyder, determina muy bien
los caudales máximos siempre se conozcan en primera instancia detalladamente las
características físicas de la cuenca, las metodologías HS triangular y SCS arrojan valores
similares entre sí, al igual que los caudales máximos, pero con volúmenes de escurrimiento
diferentes. La ventaja de un método estadístico puntual es que los caudales máximos se
calculan en base a eventos reales y la desventaja es que, a registros largos de caudales, la
extrapolación de caudales para periodos de retorno mayores a las longitudes de registro debe
tomarse con criterio y reversa para ser utilizados y por último la función de probabilidad que
mejor se ajustó a los datos de las estaciones tomadas para el estudio fue la distribución Log-
Normal de dos parámetros.
José David Gamba Mejía, José Gabriel González Hernández en su trabajo de grado: Curvas
De Rendimiento De Caudales Máximos Instantáneos Asociados A Diferentes Periodos
De Retorno En Colombia, Caso De Estudio Río Magdalena. En este estudio se realizó la
construcción de las curvas de rendimiento hídrico para caudales máximos instantáneos y la
17
determinación de la distribución que mejor representa los mismos para facilitar la aplicación
de estos estudios en actividades de ingeniería hidrológica e hidráulica en el rio Magdalena,
con este fin se proporcionó parámetros de diseño y prevención a la hora de realizar
construcciones dentro de la cuenca, y a su vez, utilizarse como criterio hidrológico para
estimar caudales en aquellos sitios donde no se posea información de los mismos. Se
concluyo en base al análisis realizado, la distribución de probabilidad que mejor se ajusta a
la cuenca del Magdalena es la distribución hidrológica GEV (método de momentos
ponderados) debido a que fue la que mejor asimilo el comportamiento de los datos con un
valor promedio de 2.43 y una desviación estándar de 1.941, cabe resaltar que se pueden
utilizar las distribuciones de probabilidad PEARSON TIPO III (Método de momentos) y
LOG-PEARSONTIPO (Método Sam) cuyas desviaciones estándar fueron 2.187 y 2.152
respectivamente, bastante cercanas a la distribución probabilística que mejor se ajustó.
Ubaldo Andres Diaz Arias, Daniel Felipe Guacaneme Prieto, Francy Lorena Alonso Nore en
su trabajo de grado: Estimación De Caudales Máximos Para El Sector Hidrológico De
La Quebrada El Diamante, Mediante Modelación Con Hec-Hms (cita). se efectuó la
modelación hidrológica de la quebrada el Diamante para efectos de estimar el caudal máximo
mediante HEC-HMS, donde La metodología empleada para realizar el proyecto de grado
consistió la obtención de un DEM de 10 metros, con el cual se obtuvo la delimitación y los
parámetros de la cuenca utilizando el software ArcGIS, posteriormente se realizó la selección
de las estaciones de influencia según los polígonos de Thiessen. Posteriormente, a los datos
de cada estación se le realizo la búsqueda de datos atípicos y se encontró la precipitación
máxima en 24 horas para los ultimo 10 años disponibles. Luego, se identificó la cobertura
vegetal y el tipo de suelo presente en la cuenca, para determinar el número de curva, además
se crearon las curvas IDF sintéticas por el método regionalizado para Colombia, con las
intensidades para un tiempo de retorno de 100 años, se elaboró el hidrograma de bloque
alterno, el hietograma de precipitación neta e infiltración, los hidrogramas de caudales con
HEC-HMS y de forma manual, además de hallar el caudal máximo por el método racional.
Como resultado de este estudio se obtuvo que: El método del hidrograma unitario triangular
(utilizado para realizar el proceso de forma manual) y el hidrograma unitario curvilíneo
(utilizado por el software HEC-HMS) arrojan resultados similares. Con el análisis de método
racional, se determina que el caudal máximo encontrado, es mayor que los otros caudales
máximos encontrados con los otros medios evaluados (manual y HEC HMS). Lo anterior
porque no se cumple con la restricción del área establecida para el método racional de que el
área de la cuenca debe ser menor a 80Ha.
18
3. FUNDAMENTO METODOLÓGICO
3.1 METODOLOGÍAS
Existen metodologías empíricas y algunas con antigüedad para estimar el cálculo de caudales
de crecientes máximos instantáneos anuales de crecientes dentro de estas se encuentran las
siguientes:
Si se cuenta con una estación hidrométrica sobre una corriente de agua de la hoya
hidrográfica hasta el sitio de proyecto de cruce con la vía que se estudia, con datos históricos
de caudal máximo instantáneo, se toman sus valores máximos instantáneos anuales, se
ajustan a distribuciones probabilísticas conocidas, y se hacen inferencias estadísticas, como
se indica en el numeral 2.3.4. Se debe tener en cuenta que los caudales máximo instantáneos
son los valores de caudal máximo que cruzan por una sección hidrométrica en una corriente,
y son registrados en estaciones limnigráficas a través de datos de lecturas de miras máximas
instantáneas, integradas con la respectiva curva de calibración de caudales líquidos. En las
estaciones limnimétricas solamente se pueden calcular valores puntuales de caudal en las
horas del día en las cuales se leen las lecturas de mira de las mismas.
Es posible que la estación hidrométrica sobre la corriente de agua que cruza la vía que se
estudia no se encuentre exactamente en el sitio de proyecto de la obra de drenaje vial, sino
en la misma hoya hidrográfica, en otra ubicación.
Donde:
19
AEH= Área hoya hidrográfica hasta la estación hidrométrica, en kilómetros cuadrados
(Km2).
El exponente x es un valor que fluctúa usualmente entre 0.5 y 0.75. A falta de datos de
investigación, se acostumbra tomar un valor igual a 0.5.
En todo caso, para la aplicación de esta metodología es imperativo que las áreas de drenaje
no sean muy diferentes. En general, esta diferencia no podrá ser mayor o menor al 50 % del
valor original del área de drenaje.
La misma metodología se puede aplicar para hoyas hidrográficas que sean hidrológica y
climatológicamente homogéneas.
Esta metodología se desarrolla cuando se cuenta con varias estaciones hidrométricas con
datos de caudales máximos instantáneos anuales históricos en una región de cierta magnitud
e hidrológicamente homogénea.
1. Se calculan los valores de caudal diferentes periodos de retorno en cada una de las
estaciones hidrométricas.
2. Se calcula la relación generalizada para la región, entre el caudal con periodo de
retorno de 2.33 años Q2.33 (caudal promedio en la distribución Gumbel) y el área de
drenaje A,con base en los datos de cada una de las estaciones hidrométricas.
Q2.33=f(A)
3. Para cada estación, se calcula la relación entre el caudal para cada periodo de retorno
T definido (dando valores de periodos de retorno iguales a 2, 5, 10, 20, 50 y 100 años)
y el caudal con periodo de retorno de 2.33 años QT/Q2.33).
4. Se debe realizar un ensayo (test) de homogeneidad estadístico, con el fin de asegurar
que los datos obtenidos correspondan a una región hidrológicamente homogénea. Si
el resultado de una de las estaciones queda por fuera de la franja de confianza del
referido ensayo, se debe eliminar tal valor. La pregunta de si los datos en un grupo de
estaciones son homogéneos, se puede contestar en un sentido estadístico
determinando si ellos difieren de uno a otro por cantidades que no pueden ser
explicadas al azar. Por otro lado, donde esas diferencias no son más que debidas a la
casualidad, se puede concluir que los datos representan aspectos diferentes de la
misma entidad y, por consiguiente, pueden ser tenidos como homogéneos.
5. Se calcula el valor promedio o el valor mediano de la relación QT/Q2.33 para la región
para los diferentes periodos de retorno considerados, con base en los valores
individuales para cada estación.
6. Para una estación no aforada en la región con un área de drenaje A, el valor del caudal
con un periodo de retorno determinado QT se calcula estableciendo, en primer lugar,
20
el valor del Q2.33 con la relación generalizada Q2.33= f (A) para la región y, luego,
utilizando la relación generalizada QT/Q2.33 para la referida región.
Así mismo encontramos otras metodologías para determinar el cálculo de los caudales y
crecientes máximos instantáneos anuales de crecientes.
El uso del método racional solo se puede usar para áreas de drenaje con áreas no
mayores a 200 hectáreas. (Fuente….)
El método del hidrograma en cuencas con áreas muy grandes, es posible desarrolla
una división en subcuencas homogéneas.
En el método donde se usa ecuaciones Log Pearson III el análisis del caudal máximo
en una cuenca se debe contar con datos de por lo menos 10 años de registro continuo
por 10 años, esto para estimaciones de cálculo con tiempos de retorno de 25 y 100
años.
DESCRIPCIÓN MÉTODO
DE Descarg
Análisis de Ecuaciones Hidrogra Raciona
INSTALACION a pico
frecuencia de regresión ma SCS l
ES de SCS
Instalacion
Corriente de flujo, es
Método
canales, puentes y Método Método complejas
alternati
alcantarillas habitual alternativo o
vo
mayores 48" hidrograma
necesitada
Instalacion
Método
Alcantarillas 48" y Método es Método
alternati
más pequeño alternativo complejas habitual
vo
o
21
hidrograma
necesitada
Instalacion
Desagües de
es
tormentas,
complejas Método
carretera zanjas y
o habitual
alcantarillas
hidrograma
laterales
necesitada
Evaluaci
ón
prelimin
Cuencas de Método
ar y
detención habitual
revisión
de
permiso
Fuente:xxxxx
Es necesario mencionar cada una de als metodologías existentes para el cálculo de caudal
máximo de una cuenca, estos aportan información académica al desarrollo del proyecto de
grado, sin embargo, no se aplican debido a las restricciones que tiene cada una de estas.
- Método Racional:
Ahora bien, este método se implementa si la cuenca cumple con un parámetro morfométrico
el cual corresponde a que la cuenca a modelar no debe tener un área superior a 1 km2.
Sin embargo, mencionamos que el método racional tiene variables tales como la escorrentía
superficial, la intensidad de precipitación y el área de la cuenca. La siguiente ecuación es la
usada en este método:33.
𝑄 = 0.278𝑥 𝐶𝐼𝐴
Donde:
22
(33OSPINA GARCÍA, Michael Jeffrey. MANCIPE GIRONZA, Juan Sebastián.
Estudio Morfométrico y Estimación de caudal de creciente de la quebrada Canauchita
del municipio de Otanche-Boyacá utilizando el método racional. Trabajo de grado
para optar al título de Tecnólogo en Construcciones Civiles. Universidad Distrital
Francisco José de Caldas. Facultad tecnológica. 2015.29 p)
A continuación, se denotan los fundamentos teóricos para determinar las variables del
método racional. [Schmidth, 1986: p.356]
𝐾𝑥 𝑇 𝑚
𝑖=
(𝑑 + 𝑐)𝑛
Donde k, m, n y c son constantes que se calculan mediante un análisis de correlación lineal
múltiple; T: Período de retorno en años, d: Duración en minutos e i: Intensidad en mm/h.
Corresponde al tiempo en que una gota de agua tarda en recorrer la longitud desde el punto
más distante del rio principal hasta el punto de cierre de la modelación. Estos tiempos de
concentración se calculan a partir de las características físicas de la cuenca, los cuales
corresponde a las pendientes, longitudes, elevaciones medias y el área de la cuenca. Es
necesario mencionar que todas las formulas tienen factores de corrección que dependen de
la cobertura de la cuenca [German Monsalve, 1999: p.180]
Desarrollada a partir de la información del SCS en siete cuencas rurales en Tennessee con
canales bien definidos y pendientes empinadas (de 3% a 10%)
En donde
23
L: Longitud máxima del río desde aguas arriba hasta la salida o punto de cierre, (en
metros)
S: Pendiente del cauce o H/L (m/m) donde H es la diferencia de elevación entre el
punto más elevado y el punto de interés.
En donde
Dicha ecuación se desarrolla para cuencas montañosas en California, formula conocida como
la Formula de California Culverts Practice. [Ven Te Chow, 1994: p. 96]
0.285
L3
Tc = 0.94788 ( ) (𝐻𝑜𝑟𝑎𝑠
H
En donde
24
Sin embargo, la mayoría de las cuencas en Colombia cuentan con estaciones hidrométricas,
así como de registros pluviográficos mínimos. Lo anterior genero la necesidad de desarrollar
métodos con los cuales sea posible obtener hidrogramas unitarios los cuales usen datos de
las características generales de la cuenca modelada. Estos hidrogramas se denominan
hidrogramas sintéticos. [Francisco Aparicio, 2001: p.228]
Esta metodología se utiliza para áreas de drenaje mayores a 2.5 Km2, es el hidrograma de
escorrentía superficial total resultante de un volumen unitario de lluvia neta, uniformemente
distribuido en espacio y tiempo. La altura de la lluvia neta o efectiva correponde con la altura
de escorrentía superficial total del hidrograma unitario.
Donde:
d= Lluvia neta o efectiva total, en milímetros (mm).
A= Área de drenaje de la hoya hidrográfica, en Kilómetros cuadrados (Km2).
Q0= Escorrentía superficial total, en metros cúbicos por segundo (m3/s).
T= Tiempo base del hidrograma unitario de la hoya, en horas (h).
Con estos principios, Sherman introdujo el denominado hidrograma unitario, el cual es una
herramienta útil en la transformación de datos de lluvia a caudal.
Los hidrogramas de escorrentía superficial en una hoya se pueden calcular por dos
metodologías, a saber: con datos históricos de creciente o mediante metodologías
relacionadas con hidrogramas unitarios sintéticos.
25
La metodología es la siguiente:
El hidrograma unitario sintético triangular fue desarrollado en 1958 por Mockus, el cual
busca un caudal pico Qp, en función del área de la cuenca (km2) y en el tiempo base (horas):
Mockus concluye que el tiempo base y el tiempo pico se relacionan mediante la siguiente
expresión:
En donde:
26
(Horas)
Donde:
Además, la duración:
𝐶𝐴
𝑄𝑝 =
𝑇𝑝
C= 2.08 (283.4 en el sistema inglés de unidades).
A= Área de drenaje en kilómetros cuadrados.
Tp= tiempo de ocurrencia del pico.
𝑡𝑟
𝑇𝑝 = + 𝑡𝑝
2
El hidrograma unitario sintético del Soil Conservation Service se representa de la siguiente
manera:
27
Hidrograma Unitario Sintético de Snyder. [German Monsalve, 1999: p.219].
La estimación del hidrograma Unitario de Snyder es procedente si no se cuenta con los datos
históricos de caudal y precipitación para el cálculo del hidrograma unitario de una cuenca.
𝑡𝑟 = 𝑡𝑝 /5.5
Donde:
tr: Duración de la lluvia neta (horas).
0.275 𝐶𝑝 𝐴
𝑞𝑝 =
𝑡𝑝
Donde:
28
El hidrograma unitario sintético obtenido representa a un milímetro de precipitación neta
sobre toda la cuenca.
𝑡𝑝
𝑇 = 3+
8
Para otra lluvia neta de duración tR diferente a tP el tiempo de retardo correspondiente a tpR
el caudal pico qpR y el caudal base TR según Linsley son:
Adicionalmente,
Calculando este hidrograma con valores de ll neta de duración tr o con una duración tR en
donde:
W75: Ancho del HU, horas, un valor de las ordenadas igual al 75% del caudal pico
qP o qPR.
W50: Ancho del HU, horas, un valor de las ordenadas igual al 50% del caudal pico
qP o qPR.
29
El valor del tiempo de retardo de la cuenca tp presenta valores iguales al tiempo entre el
centro geométrico de la duración de la lluvia neta tr y el pico del hidrograma unitario (horas)
basados en las características tP, tr,qP, W75, W50 y T (ó tPR, trR, qPR, W75R, W50R, y TR).
Si no se cumple se debe ajustar el hidrograma (ajustando el tiempo base T o TR) hasta cumplir
con la anterior ecuación.
En donde:
El valor del tiempo promedio de ascenso del hidrograma unitario, Ta, se calcula de acuerdo
con la siguiente ecuación para hoyas hidrográficas con uso del suelo predominantemente
forestado:
30
Y de acuerdo con la siguiente ecuación, para las hoyas hidrográficas con uso del suelo
predominantemente agrícola:
Donde
Donde:
3.1.1 Limitaciones
Método racional
La suposición del método racional que limite su uso a las áreas de drenaje pequeñas de hasta
80.9 hectáreas incluyen:
31
La fracción de lluvia que se convierte en escorrentía [C] depende de la intensidad o
volumen de la lluvia. Esta suposición es razonable para áreas impermeables, como
calles, tejados y estacionamientos. Para áreas permeables, la fracción de escorrentía
varía con la intensidad de la lluvia y el volumen acumulado de lluvia. El coeficiente
de escorrentía seleccionado debe ser apropiado para las condiciones de tormenta,
suelo y uso del suelo.
Método SCS
32
Cuando el CN ponderado es de menos de 40 o más de 98, utilice otro procedimiento
para calcular el gasto pico.
Si el área de drenaje es mayor que 2589.99 hectáreas, la precipitación de punto debe
ajustarse usando el procedimiento del capítulo 4, NEH-4 (SCS, 1985) o TP-40
(Hershfield, 1961).
Las características analizadas de la cuenca deben ser dentro de los límites de los que se
utilizan para desarrollar las ecuaciones. Cada informe del USGS tiene una tabla que muestra
el límite superior e inferior para la característica de cuencas utilizada en el análisis de
regresión; se debe tener cuidado usando ecuación fuera de los límites especificados. Las
ecuaciones de regresión se pueden utilizar en secuencia que no son significativamente
afectados por la regulación artificial, urbanización o desviación. Las ecuaciones no se deben
utilizar inmediatamente aguas abajo de un lago o cuerpo de agua, sino usados para determinar
el flujo pico se enruten a través de un embalse.
33
METODOLOGÍAS PARA CAUDAL MÁXIMO
MÉTODO MÉTODO DIRECTO MÉTODO CREAGER MÉTODO DEL H.U. MÉTODO
RACIONAL (PERÚ) TRIANGULAR DEL
SISTEMA
DIPEO
Q=CIA/360 Q=(1/n)*(AR)^(2/3)*S^(1/2) Q=(C1+C2)*LOG(T)*A^(mA^- Qp=(0,208)*(A*Pe/tp) Q=P(34H)*
n) A* Cr* Cp/r *
Cf/Tc
C: n:coeficiente de rugosidad A: área de la cuenca aportante tb=2.67tp P:
coeficiente (km2) precipitación
de máxima en 24
escorrentía horas (mm)
I: intensidad A:area hidráulica promedio T: periodo de retorno (años) tp=(de/2)+tr A: área de la
máxima de (m2) cuenca (km2)
lluvia
PARÁMETROS
(mm/hr)
A: área de la R: radio hidraulico promedio C1,C2: coeficientes por Cr: coeficiente
cuenca (ha) (m) regiones hidráulicas de reducción
S:pendiente m/m m,n: exponentes por regiones tr=0.6tc = 0.005* Cr/p:
hidráulicas (L/(raiz(s))^0.64 coeficiente de
P=P24h*(D/1440)^0.25 reducción
según
permeabilidad
y pendiente
Longitudinal y
transversal
Pe=(P-0.20S)^2/(P+0.80S) Cf: coeficiente
(in) que considera
la forma de la
curva del
caudal
34
METODOLOGÍAS PARA CAUDAL MÁXIMO
MÉTODO MÉTODO DIRECTO MÉTODO CREAGER MÉTODO DEL H.U. MÉTODO
RACIONAL (PERÚ) TRIANGULAR DEL
SISTEMA
DIPEO
s=(1000/CN)-10 tc: tiempo de
concentración
(sg)
A: área de la cuenca (km2)
Pe: altura de precipitación en
exceso (mm)
tp: tiempo pico (hr)
tb: tiempo base (hr)
de: duración de exceso (hr)
tc: tiempo de concentración
(hr)
tr: tiempo de retraso (hr)
L: longitud de cauce principal
(m)
S: pendiente del cauce (%)
P: precipitación de diseño (in)
s: atracción inicial
CN: curva numero
cuencas >1 se debe contar con los solo para Perú 200 KM2 >
LIMITACIONES
35
METODOLOGÍAS PARA CAUDAL MÁXIMO
MÉTODO MÉTODO DIRECTO MÉTODO CREAGER MÉTODO DEL H.U. MÉTODO
RACIONAL (PERÚ) TRIANGULAR DEL
SISTEMA
DIPEO
36
Por otro lado se desarrolla un análisis adicional la cual se plasma en el Anexo 1.
Metodologías para cálculo de caudal máximo en cuencas.
Las herramientas utilizadas para este trabajo fueron los software en los que se modelaron
la cuenca y los programas básicos del office mencionados a continuación: ARCGIS,
HEC-HMS, WMS, y los componentes del office como Word y Excel.
37
4. ANÁLISIS DE RESULTADOS
Siguiendo con este análisis es importante mencionar que esta cuenca tiene # corrientes,
donde el rio cañas (corriente principal) tiene una longitud de 12.62 km,
38
De igual manera la cuenca objeto del presente estudio se caracteriza por tener una
densidad de corriente (Ds) 0.3 ríos/km y una densidad de drenajes (Dd) 0.49 siendo la
cuenca pobremente drenada.
Sumado a lo anterior la forma de la cuenca tiende a ser alargada esto debido a que el
factor forma de esta corresponde a 0.38, así mismo de acuerdo al índice de alargamiento
corresponde a 2.63 lo que indica que es una cuenca alargada de este modo se deduce que
los tiempos de viaje son largos. Lo que nos permite deducir que la cuenca tiene menor
susceptibilidad a la formación de crecientes
Por otro lado de acuerdo al Anexo 2 el tipo de red de drenaje tipo detrítico esto de acuerdo
a la guía básica para la caracterización morfométricos de cuencas hidrográficas en la cual
menciona “los tributarios se unen a la corriente principal formando ángulos. Su
presencia indica suelos homogéneos” (Reyes Trujillo, Barroso, & Carvajal Escobar,
2010, pág. 62)
PORCENTAJE PORCENTAJE
DEL TOTAL SOBRE EL LIMITE
(%) INFERIOR (%)
0.965 100
2.229 99.035
4.415 96.806
4.849 92.391
4.680 87.542
4.540 82.862
7.264 78.322
6.943 71.058
8.655 64.116
10.994 55.460
12.952 44.466
14.010 31.515
10.564 17.505
4.742 6.941
2.199 2.199
100
39
Grafica 1. Curva Hipsométrica. Fuente: Autores
40
De acuerdo a la Grafica 1 se analiza que la cuenca se halla en una etapa joven tendiendo
a presentar más procesos erosivos, confirmando esto los diferentes procesos erosivos
presentados a lo largo del municipio de San Francisco.
subcuenca TC (hr) TR ( TC TR
retardo) (min) (min)
SUB1 2.066499 0.723275 123.99 43.39649
SUB2 1.597608 0.559163 95.85648 33.54977
SUB3 1.405044 0.491766 84.30267 29.50593
SUB4 1.148427 0.40195 68.90564 24.11697
CUENCA 5.289245 1.851236 317.3547 111.0741
Tabla 7. Cálculos obtenidos para cada una de las subdivisiones de la subcuenca.
Fuente: Autores
4.3 PRECIPITACIÓN
42
Francisco. Fuente: http://www.ideam.gov.co/solicitud-de-informacion
PRECIPITACIÓN
AÑO
(mm)
1987 44.0
1988 100.1
1989 68.0
1990 92.0
1991 52.0
1992 75.0
1993 68.2
1994 105.0
1995 83.0
1996 36.8
1997 51.7
1998 54.8
1999 67.0
2000 69.0
2001 57.0
2002 81.1
2003 55.0
2004 49.3
43
2005 69.2
2006 70.0
2007 71.5
2008 69.5
2009 57.6
2010 85.1
2011 79.5
2012 51.3
2013 72.5
2015 49.3
Número de
28
Registros
Media 67.30
Mediana 68.60
Máximo 105.00
Mínimo 36.80
Desviación
16.66
Estándar
Tabla 8. Datos de precipitación máxima en 24 horas. Fuente: Autores
DATOS SABANETA
ESTADÍSTICOS
Min (mm) 36.8
Q1 (mm) 52.7000
Q2 (mm) 68.6000
Q3 (mm) 78.375
Max (mm) 105.0
Rango RI (mm) 25.6750
a -24.325
b 14.1875
c 116.8875
d 155.4
Tabla 9. Datos y cálculos estadísticos de la estación Sabaneta, Municipio San
Francisco, Cundinamarca. Fuente: Autores
44
Grafica 2. Caja de bigotes. Fuente: Autores
WATER RESOURCES
COUNCIL
Media (mm) 1.8
DesvEst (mm) 0.109
Coef Asi (mm) -0.174
Kn (mm) 2.534
Yh (Umbral) (mm) 2.091
Xh (mm) 123.334
45
Xh sin desface 122.334
(mm)
Tabla 10. Metodo water resources Council. Fuente: Autores
Para la modelación de la cuenca es necesario obtener las curvas IDF con tiempos de
retorno de 2, 5, 10, 25, 50 y 100 años. A continuación se presenta la gráfica resultante
del cálculo realizado.
46
47
Grafica 4. Curva IDF para la subcuenca hidrográfica hasta el municipio San Francisco, Cundinamarca. Fuente: Autores
48
Es necesario mencionar que la subcuenca hasta el municipio San Francisco,
Cundinamarca pertenece a la región andina por tanto se adoptan valores para el cálculo
de esta curva. De acuerdo a los Valores de los coeficientes a, b, c y d para el cálculo de
las curvas intensidad-duración-frecuencia, IDF, para Colombia.
Coeficiente de la
Ecuación para la
región Andina (R1)
a 0.94
b 0.18
c 0.66
d 0.83
Tabla 11. Valores de los coeficientes a, b, c y d para el cálculo de las curvas
intensidad-duración-frecuencia, IDF, para Colombia. (INVIAS, INSTITUTO NACIONAL
DE VIAS, 2009, págs. 2 - 71)
Una vez obtenidos los datos de intensidad para diferentes periodos de retornos, se
procedió a realizar el cálculo de la lluvia de diseño por bloque alterno.
hietograma
hietograma
49
dividió en cuatro subdivisiones, a continuación se presenta la figura en donde se puede
evidenciar esta subdivisión.
SUBDIVISIÓN A(Km2)
sub_1 38.266
sub_2 30.496
sub_3 8.830
sub_4 5.515
Tabla 13. Área de las subdivisiones de la subcuenca modelada. Fuente: Autores
50
C 0.5064 0.5342 0.5248 0.519
CN 88.7692 89.4716 89.38723 91.75909
D 6.980132 6.231229 3.35309 2.649847
S(m/km) 105.7 125.6 134 170.5
H 0.916313 0.820168 0.761388 0.818912
Tabla 14. Datos de cada una de las subdivisiones de la subcuenca para la
modelación del hidrograma unitario SCS. Fuente: Autores
51
Imagen 6. Resultado del programa HEC-HMS para el hidrograma unitario SCS,
tiempo de retorno 2 años. Fuente: Autores.
52
Imagen 7. Resultado del programa HEC-HMS para el hidrograma unitario SCS,
tiempo de retorno 10 años. Fuente: Autores.
53
Imagen 8. Resultado del programa HEC-HMS para el hidrograma unitario SCS,
tiempo de retorno 25 años. Fuente: Autores.
54
Imagen 9. Resultado del programa HEC-HMS para el hidrograma unitario SCS,
tiempo de retorno 50 años. Fuente: Autores.
Grafica 9. Hidrograma unitario SCS para tiempo de retorno 100 años. Fuente:
Autores
55
Imagen 10. Resultado del programa HEC-HMS para el hidrograma unitario SCS,
tiempo de retorno 100 años. Fuente: Autores.
56
PEAK DISCHARGE: 22.5 M3/S
Grafica 10. Hidrograma unitario Snyder para tiempo de retorno 2 años. Fuente:
Autores
57
Grafica 11. Hidrograma unitario Snyder para tiempo de retorno 10 años. Fuente:
Autores
58
Grafica 12. Hidrograma unitario Snyder para tiempo de retorno 25 años. Fuente:
Autores
59
Grafica 13. Hidrograma unitario Snyder para tiempo de retorno 50 años. Fuente:
Autores
60
Grafica 14. Hidrograma unitario Snyder para tiempo de retorno 100 años.
Fuente: Autores
Grafica 15. Hidrograma unitario Clark para tiempo de retorno 2 años. Fuente:
Autores
62
Tiempo de retorno 10 años
Grafica 16. Hidrograma unitario Clark para tiempo de retorno 10 años. Fuente:
Autores
63
Grafica 17. Hidrograma unitario Clark para tiempo de retorno 25 años. Fuente:
Autores
64
Grafica 18. Hidrograma unitario Clark para tiempo de retorno 50 años. Fuente:
Autores
65
Grafica 19. Hidrograma unitario Clark para tiempo de retorno 100 años. Fuente:
Autores
66
5. CONCLUSIONES Y RECOMENDACIÓNES
Los cálculos máximos obtenidos para cada uno de los tiempos de retorno así como
para cada uno de los métodos se denota que no presentan desviaciones
significativas en los resultados de un método a otro, lo que se concluye que la
relación de estos métodos.
Como es evidente los resultados del cálculo de los caudales máximos para cada
uno de los tiempos de retorno los datos del método del SCS presentan valores
mayores comparando de los métodos de Snyder y de Clark.
67
BIBLIOGRAFÍA
68
Press.
Ministerio de Vivienda, Ciudad y Territorio. (200). Reglamento Técnico del Sector de Agua
Potable y Saneamiento Básico - RAS. BOGOTÁ D.C.
Reyes Trujillo, A., Barroso, F. U., & Carvajal Escobar, Y. (2010). Guía basica para la
caracterizacion morfometrica de cuencas hidrograficas. Cali, Colombia.
UNIVERSIDAD DE CUNDINAMARCA. (2017). FORMULACIÓN DEL PLAN
MUNICIPAL DE GESTIÓN DEL RIESGO DE DESASTRES MUNICIPIO DE SAN
FRANCISCO CUNDINAMARCA. BOGOTA.
US ARMY CORPS OF ENGINEERS. (s.f.). http://www.hec.usace.army.mil. Obtenido de
http://www.hec.usace.army.mil/software/hec-hms/:
http://www.hec.usace.army.mil/software/hec-hms/
VÁSQUEZ RANGEL, C., HERRERA LÓPEZ, D., & GUTIÉRREZ REY, Y. (2014).
Caracterización morfométrica de la cuenca de la quebrada Tatamaco, del municipio
Villavieja del departamento del Huila, mediante el uso de la herramienta HEC-
GEOHMS. Bogotá D.C.
APÉNDICE
ANEXOS
69