Puente San Juan-Vides V1 1
Puente San Juan-Vides V1 1
Puente San Juan-Vides V1 1
Abril de 2022
Página 1 de 156
CONTENIDO
INTRODUCCIÓN .................................................................................................................................................. 4
1. OBJETIVOS ................................................................................................................................................... 4
1.1. OBJETIVO GENERAL ......................................................................................................................... 4
1.2. OBJETIVOS ESPECÍFICOS ............................................................................................................... 4
2. GENERALIDADES ........................................................................................................................................ 5
2.1. GEOLOGÍA ........................................................................................................................................... 6
2.1.1. Litología ........................................................................................................................................ 7
2.1.2. Tectónica ....................................................................................................................................10
2.2. GEOMORFOLOGÍA ...........................................................................................................................10
2.2.1. Condiciones geomorfológicas del área de estudio .............................................................12
2.2.2. Evolución geomorfológica .......................................................................................................13
2.2.3. Jerarquización geomorfológica ..............................................................................................14
2.2.4. Morfodinámica del área de estudio........................................................................................19
2.2.5. Diacronismo fluvial ....................................................................................................................24
2.2.6. Patrón morfológico actual ........................................................................................................30
2.2.7. Geometría hidráulica ................................................................................................................37
2.2.8. Unidades de canal ....................................................................................................................42
2.2.9. Caracterización inicial de los cambios de tipo morfológico en las últimas décadas .....44
2.2.10. Identificación de tramos modificados por acciones directas en el cauce .......................45
2.2.11. Síntomas de dinámica vertical acelerada .............................................................................46
2.2.12. Estructura y sustrato del lecho ...............................................................................................48
3. HIDROGRAFÍA ............................................................................................................................................56
3.1. HIPSOMETRÍA DE LA CUENCA .....................................................................................................57
3.2. PARÁMETROS MORFOMÉTRICOS DE LA CUENCA.................................................................57
3.2.1. Parámetros de forma ................................................................................................................58
3.2.2. Parámetros hidrográficos ........................................................................................................62
3.2.3. Tiempo de concentración ........................................................................................................65
3.2.4. Coberturas del suelo ................................................................................................................65
4. HIDROLOGÍA ..............................................................................................................................................66
4.1. MECANISMOS DE CÁLCULO.........................................................................................................67
4.2. ECUACIÓN DE INTENSIDAD ..........................................................................................................69
4.3. CÁLCULOS HIDROLÓGICOS .........................................................................................................70
4.3.1. Análisis meteorológico .............................................................................................................73
Página 2 de 156
4.4. CURVAS DE INTENSIDAD – dURACIÓN Y FRECUENCIA (iDF) ...............................................90
4.5. METODOLOGÍA MODELACIÓN HIDROLÓGICA ........................................................................94
4.5.1. Determinación del área hidrológica del estudio ..................................................................94
4.5.2. Estimación del índice de capacidad de infiltración y escurrimiento del suelo ................94
4.5.3. Determinación del programa de cómputo hidrológico .......................................................95
4.5.4. Análisis estadístico de la información de la estación seleccionada ..................................95
4.5.5. Cálculo de caudales .................................................................................................................96
4.5.6. Modelación hidrológica ......................................................................................................... 108
4.5.7. Calibración del modelo lluvia – escorrentía para validar el valor de las crecientes
calculadas con el método S.C.S. .......................................................................................................... 109
4.6. RESULTADOS HIDROLÓGICOS ................................................................................................. 110
4.6.1. Análisis de frecuencia de lluvias máximas en 24 horas y determinación de aguacero de
diseño 111
4.6.2. Determinación de las crecientes para el periodo de retorno de diseño ....................... 113
5. SIMULACIÓN HIDRÁULICA .................................................................................................................. 130
5.1. PARAMETRIZACIÓN...................................................................................................................... 133
5.1.1. Esquematización .................................................................................................................... 135
5.1.2. Implementación del modelo ................................................................................................. 135
5.1.3. Análisis de sensibilidad ......................................................................................................... 136
5.1.4. Verificación del modelo......................................................................................................... 136
5.1.5. Aplicaciones del modelo ....................................................................................................... 137
5.1.6. Condiciones de frontera ....................................................................................................... 137
5.2. IMPLEMENTACIÓN DEL MODELO ............................................................................................. 137
5.2.1. Configuración del dominio 1D ............................................................................................. 137
5.2.2. Condiciones de borde e iniciales del modelo ................................................................... 138
5.2.3. Coeficientes de rugosidad .................................................................................................... 139
5.3. RESULTADOS DE LA MODELACIÓN ........................................................................................ 140
5.3.1. Esquema de planta de inundación ...................................................................................... 140
5.3.2. Esquema de velocidades ...................................................................................................... 142
5.3.3. Tirante hidráulico ................................................................................................................... 143
5.4. ANÁLISIS DE SOCAVACIÓN ....................................................................................................... 145
5.5. DISEÑO FLUVIAL ........................................................................................................................... 148
CONCLUSIONES Y RECOMENDACIONES ............................................................................................... 151
BIBLIOGRAFÍA ................................................................................................................................................. 154
Página 3 de 156
INTRODUCCIÓN
De esta manera, los capítulos a continuación describen los aspectos que relacionan la revisión de
información secundaria, el análisis de hidrología, los detalles de la hidrografía general a nivel de
microcuenca, las actividades de reconocimiento - operación en campo y las conclusiones finales de
los resultados obtenidos.
1. OBJETIVOS
1.1. OBJETIVO GENERAL
El estudio tiene como principal objetivo el estudio hidrológico e hidráulico para la construcción del
puente localizado en el paso del río San Juan en las veredas San Juan Vides y Damasco Vides en los
municipios de Orito y Puerto Caicedo.
Página 4 de 156
- Determinar el sitio óptimo de localización de la estructura a través de un análisis de dinámica
fluvial.
2. GENERALIDADES
La zona de estudio se localiza en los Municipios de Orito y Puerto Caicedo, sobre el perímetro de
rural, en las veredas San Juan Vides y Damasco Vides. Cartográficamente se localiza en la plancha
IGAC (1:25.000) 449-IV-A, en una extensión de análisis de 40 hectáreas y 250 metros de recorrido.
Las coordenadas que delimitan el área de estudio específico de la obra se indican en la Tabla 1:
Página 5 de 156
Tabla 1. Coordenadas del área de estudio
ORIENTACIÓN MAGNA OESTE WGS1984
Este Norte Longitud Latitud
NORTE 1042106 574176 76°41'57.217"W 0°44'42.939"N
ESTE 1042548 573914 76°41'42.918"W 0°44'34.418"N
SUR 1042142 573901 76°41'56.052"W 0°44'22.08"N
OESTE 1041664 573860 76°42'11.518"W 0°44'32.658"N
Fuente: Elaboración propia, 2022
Actualmente el sector que encuentra conformado por predios cuya cubierta del suelo predominante
son pastos limpios, vegetación secundaria en transición, y áreas con espacios naturales y cultivos -en
mosaico-, adicionalmente, existen áreas importantes conformadas por bosques productivos donde se
obtiene maderas de calidad comercial. En el marco general, la cuenca está constituida en mayor
proporción, por vegetación secundaria en transición, bosques densos en tierra firme y bosques
abiertos.
De acuerdo con las observaciones propias y entrevistas no estructuradas desarrolladas con los
morados del lugar, la actividad predominante en el área es la agrícola, seguida de la actividad de
aprovechamiento forestal de madera valiosa y tercero la ganadería doble propósito.
2.1. GEOLOGÍA
Página 6 de 156
Figura 2. Geología del área de proyecto
2.1.1. Litología
La microcuenca de la referencia dentro del área propuesta para la construcción del puente ésta se
encuentra localizada dentro de una unidad litológica correspondiente al periodo Paleógeno-Neógeno,
entendiéndose como Grupo Orito, muy aferente a las terrazas secundarias del periodo Cuaternario.
Dentro de la memoria descriptiva de la plancha geológica 449 el (INGEOMINAS, 2002) refiere que el
Grupo Orito fue definido por Miley & McGirk (1948) en el río Orito; las primeras referencias, en la
Cuenca del Putumayo, se remontan a reportes de Texas Petroleum Company, en los cuales el Grupo
Página 7 de 156
Orito está conformado por las formaciones Orteguaza y Belén (McGirk, 1949). Otros reportes internos
de la Texas indican que el grupo está constituido por las formaciones Orteguaza, Belén y Orito, aunque
algunas veces nombran la unidad como Formación Orito solamente. Reportes de Shell subdividen el
Grupo Orito en las formaciones Orteguaza y Belén - Orito.
Localización: Los afloramientos del Grupo Orito están restringidos a la llanura amazónica, en donde
forma un relieve ondulado, debido fundamentalmente a su composición litológica que no permite
generar relieves sobresalientes. Aunque el Grupo Orito ocupa gran parte de la zona suroriental del
área cartografiada (planchas 430 Mocoa, 431 Piamonte, 449 Orito y 465 Churuyaco), las exposiciones
por lo regular son pobres y están enmascaradas por depósitos cuaternarios; adicionalmente, salvo los
afloramientos cercanos al piedemonte, la estratificación de la unidad tiene un buzamiento muy bajo a
subhorizontal y horizontal, lo que imposibilita, para la escala de presentación del mapa geológico,
determinar con detalle la secuencia estratigráfica y los límites de las formaciones que conforman el
grupo.
Descripción: El Grupo Orito está conformado por estratos delgados a gruesos de lodolitas, limolitas
y lodolitas arenosas de color gris, rojo, morado, amarillo y moteadas, con intercalación de estratos
gruesos lenticulares a ondulosos no paralelos y no continuos de litoarenitas de textura variable,
inmaduras a submaduras, con ondulitas y laminación inclinada en artesa, plana paralela y localmente
macizas. Las litoarenitas presentan micas, laminillas y restos de materia orgánica, óxidos de hierro,
yeso y localmente bioperturbación.
Espesor: El espesor del Grupo Orito, medido en la sección río Fragua - río Congor de la Plancha 431
Piamonte, es de 1.830 m. Una secuencia afectada por fallamiento, en el río San Juan (Plancha 449
Orito), dio un espesor de 1.546 m, mientras que en la quebrada La Chorrera (Plancha 465 Churuyaco),
se midieron 946 m, sin que se hubiera podido observar el tope de la secuencia. Con anterioridad se
tenía el reporte de McGirk (1949), que indicaba un rango entre 1.000 y 1.300 m de espesor.
Contactos: El contacto inferior con el segmento superior de la Formación Pepino fue identificado
como neto, paraconforme y gradual; en la base del Grupo Orito las capas de arcillolitas arenosas
grises marcan un cambio transicional con el Miembro Superior de la Formación Pepino, situación que
fue observada en la quebrada La Chorrera (Plancha 465 Churuyaco).
El límite superior del Grupo Orito con los depósitos cuaternarios que lo suprayacen es una
discordancia, comúnmente de tipo angular. También se reportan límites paraconforme y disconforme.
Página 8 de 156
Edad: Las asociaciones fosilíferas encontradas en las sedimentitas del Grupo Orito corresponden al
intervalo Oligoceno - Mioceno medio (McGirk, 1949). En el área cercana al Municipio de Orito (Plancha
449 Orito), las unidades correspondientes a la secuencia que se considera como Grupo Orito
(formaciones La Paloma, Ospina y San Miguel) tienen edad que va del Oligoceno medio al Plioceno
según Miley & McGirk (1948).
Las terrazas medias son el nivel intermedio de terrazas presentes en la zona; regularmente ocupan
las laderas de los valles de las corrientes de segundo y tercer orden; los constituyentes son guijos,
guijarros y bloques de tamaño medio a grueso, composición heterogénea y matriz de arena media a
gruesa.
En la Plancha 449 Orito los depósitos de terraza se encuentran en la llanura amazónica drenada por
los ríos Putumayo, Vides, Conejo, Orito, Caldero y San Juan, mientras que en la Plancha 465
Churuyaco están asociadas al curso bajo de los ríos Churuyaco, Rumiyaco y San Miguel.
Estos depósitos aluviales se encuentran siempre relacionados a los cauces actuales y activos de los
diferentes ríos que drenan la zona y descargan su caudal a la llanura amazónica. Los depósitos de
este tipo más extensos se encuentran al oriente de la ciudad de Mocoa, en el sector de Villa Garzón -
Puerto Limón en la Plancha 430 Mocoa; el río Putumayo y sus afluentes Vides, Orito, San Juan, Conejo
y otros, menores en la Plancha 449 Orito y el río San Miguel con sus tributarios Churuyaco, Rumiyaco
y Gavilanes en la Plancha 465 Churuyaco. En varios sitios estos depósitos no fueron diferenciados de
las terrazas más bajas (Qt3).
Las acumulaciones aluviales están conformadas por partículas tamaño bloque, canto, grava, arena,
limo y arcilla, de muy variada composición litológica, con o sin matriz, y refleja la constitución geológica
de las áreas drenadas.
Página 9 de 156
2.1.2. Tectónica
El proyecto no se localiza en ningún plano de falla descrito ni aparente. El desarrollo del sector solo
evidencia la prominencia de la formación Orito con las particulares lodolitas características.
2.2. GEOMORFOLOGÍA
De acuerdo con Carvajal & Servicio Geológico Colombiano (2011) la geomorfología como la ciencia ,
trata de la descripción y el estudio de la génesis, clasificación, procesos y evolución de las formas
antiguas y actuales del terreno (tanto terrestres como submarinas) y su relación con las estructuras
infrayacentes. Esta ciencia está enfocada al análisis de la superficie de la tierra en donde interactúan
la litósfera, la hidrósfera, la atmósfera y la biosfera. Igualmente las geoformas por definición son la
expresión superficial del terreno, de la interacción dependiente de los materiales constituyentes de
las mismas y su disposición estructural, de los procesos geomórficos que interactúan según el
ambiente morfogenético específico donde se desarrollan, y el tiempo de duración de la acción de los
mencionados procesos (Keller & Rockwell, 1984).
De acuerdo con los objetivos del presente estudio, la cartografía temática del componente
geomorfológico, se representa en un mapa a escala 1:2.000 (áreas urbanas) y 1:5.000 (áreas rurales),
bajo la estructuración de una base estructurada por un Sistema de Información Geográfica -SIG-,
tanto, se siguieron los lineamientos y modificaciones establecidas por Carvajal & Servicio Geológico
Colombiano (2011) donde, cada elemento cartográfico, está representado, por símbolos,
convenciones y colores, definidos en la metodología mencionada y basada en la metodología del
International Institute for Aerospace Survey and Earth Sciences -I.T.C.- propuesta por Verstappen,
Herman Th. and A van Zuidam, Robert (Verstappen, Zuidam, Meijerink, & Nossin, 1991).
La Figura 3 plantea el enfoque del análisis geomorfológico propuesto (Carvajal & Servicio Geológico
Colombiano, 2011); donde, el proceso analítico del terreno comprende desde lo regional hasta lo local,
relacionando la escala aplicada con las jerarquías geomorfológicas, donde la base regional está
fundamentada en la génesis geológica de las geoformas y los ambientes morfogenéticos, y el detalle
basado en los ambientes morfogenéticos, expresión morfológica, litología y procesos morfodinámicos.
Página 10 de 156
Figura 3. Esquema de jerarquización geomorfológica propuesto para INGEOMINAS
Para el caso de la revisión bibliográfica existente, se consultaron las bases documentales de los
centros especializados, tales como el Servicio Geológico Colombiano -SGC-; Instituto Geográfico
Agustín Codazzi -IGAC-, y la Facultad de Geología de la Universidad Nacional de Colombia -Sede
Bogotá D.C.
Junto a lo mencionado, se consultaron, vía Web, los centros de documentación de aquellas entidades
u organismos que indirectamente desarrollan el componente geomorfológico, tales como: Instituto de
Hidrología, Meteorología y Estudios Ambientales -IDEAM-, Instituto de Investigación de Recursos
Biológicos Alexander von Humboldt, Ministerio de Medio Ambiente, Autoridad Nacional de Licencias
Ambientales -ANLA-, Corporación Autónoma Regional de la Amazonía -CORPOAMAZONÍA-.
Cartográficamente, se contó con el límite de las cuencas hidrográficas, cartografía base (curvas de
nivel), cartografía temática de geológica, cartografía temática de estructuras geológicas, cartografía
temática del patrón de drenaje y cartografía temática de pendientes. Además, se contó con los
sensores remotos (Ortofotos); Modelo Digital de Sombras (DSM); y el Modelo Digital de Elevación
(DEM) levantado en este proyecto mediante Vehículo Aéreo Tripulado Remotamente – RPAS.
La Tabla 2 a continuación presenta las características de los sensores remotos empleados para la
generación del componente geomorfológico.
Página 11 de 156
Tabla 2. Características de los sensores remotos utilizados en la interpretación
SENSOR ESCALA TAMAÑO DEL FECHA DE
PIXEL PRODUCCIÓN
LANDSAT 2 1:100.000 60 m 14-Jul-1969
LANDSAT 4 1:100.000 60 m 01-Feb-1976
LANDSAT 5 1:100.000 60 m 28-Abr-1985
LANDSAT 5 1:100.000 60 m 19-Ago-1985
IKONOS 1:5.000 1M 14-Ene-2008
GEOEYE 1:5.000 1m 11-Sep-2014
WORLDWIDE-2 1:1000 50 cm 11-Feb-2022
PHANTOM 4 RTK 1:200 7 cm 13-Abr-2022
Fuente: Elaboración propia, 2022
A partir del análisis del DEM, se generaron variables geométricas (expuestas en la metodología
aplicada), que en diversos casos apoyaron la sectorización de unidades homogéneas; variables que
en un sistema de información geográfica son evaluadas para una mejor definición de las unidades
geomorfológicas. El DEM empleado se obtuvo de la clasificación de datos LAZ de la fotogrametría
La estructuración del mapa preliminar, consistió en definir las unidades geomorfológicas que
cumplieran con las características de forma, escala de mapeo (área mínima cartografiable) y
descripción de la unidad. Esta fase, está fundamentada en la interpretación de los sensores remotos
disponibles para el área de las cuencas (sensores remotos), y de las variables nombradas
anteriormente; diferenciando las geoformas por su origen (ambiente morfogenético) y las relaciones
entre las mismas.
A partir de la interpretación temática y una vez definidas las unidades de geomorfología, se procedió
a ajustar la geometría de las mismas a través del patrón de drenaje, base cartográfica y su
correspondencia con el modelo de sombras.
Como ya se describió, el área de estudio, se localiza al nor-occidente del municipio de Puerto Caicedo;
donde, el río San Juan transita en sentido de noroeste - sureste -, hasta desembocar en el río
Putumayo, a la altura del municipio de Puerto Caicedo.
El área de estudio del río San Juan, se desarrolla sobre los Ambientes de origen Fluvial (F) y
Denudacional (D), de la vertiente oriental de la cordillera occidental del Sistema orogénico andino.
En el lugar de la obra el relieve es llanura; con superficies planas o suavemente inclinadas con
contraste de terrazas ligeramente pronunciadas; de forma cóncavas o convexas, dependiendo del
ambiente morfogenético, del mismo modo, el paisaje está acompañado de la sobreelevación del
Grupo orito que sobresale en la vertiente izquierda del río. El gradiente altitudinal, oscila entre los
intervalos de 260 a 310 m.s.n.m.
Página 12 de 156
2.2.2. Evolución geomorfológica
Instituto de Hidrología Meteorología y Estudios Ambientales (2010) indica que cuando nos referimos
a los sistemas aluviales como aquellos modelados y formaciones superficiales producto de la
sedimentación por las corrientes mayores del drenaje. En este trabajo, la caracterización se basa en
el conjunto de las formas, formaciones superficiales y los procesos morfogénicos que integran el
modelado.
La característica general de los sistemas aluviales es la inundabilidad periódica por el desborde de los
ríos autóctonos y alóctonos y que definen las llanuras de inundación. En relación con la lluvia, también
el encharcamiento de estos espacios es otra de las características y especialmente en el clima
superhúmedo de la Amazonia.
Por la pendiente suave de la mayoría del espacio amazónico, y a pesar de los frecuentes controles
estructurales que diferentes tipos de lineamientos tectónicos ejercen sobre las corrientes, una de las
principales características de la red de drenaje es la divagación con formación de meandros, cauces
abandonados e incluso, desplazamientos laterales de los ríos.
La diferenciación que se presenta está ligada en primer lugar a la altura relativa de las formaciones y
secundariamente a la relación con el origen andino o amazónico de los ríos.
Esta parte de los sistemas aluviales se tomó básicamente de las unidades fisiográficas definidas para
el mapa de suelos hecho por (Gonzales Fletcher, Proyecto Radargrametrico del Amazonas, Rivera
Hermida, Diazgranados, & Jaramillo Berrio, n.d.), además de las caracterizaciones publicadas por
Botero (1980), Khobzi et ál. (1980) y Van Zuidam (1974), ubicando al río San Juan y sus afluentes de
la zona de depósito dentro de la zona de disección modera en transición a zonas inundadas, aunque
no son permanentes, podría citarse como zonas de inundación transicional y temporal por épocas
continuas.
Página 13 de 156
Figura 4. Perfil ríos Caquetá-Putumayo-Amazonas. Estructura en bloques, pendiente y grados de disección
Así las cosas, la organización de los sistemas aluviales y de interfluvio se aprecia en la Figura 4. Los
depósitos aluviales proceden, en su mayoría, de la Cordillera Oriental; sin embargo, se reconocen los
aportes por disección y escurrimiento hídrico que funcionan en las plataformas y en los afloramientos
del escudo, aunque estos aportes parecen ser menores que los de la cordillera.
Para definir las diferentes unidades geomorfológicas, fue necesario clasificar y definir cada unidad de
acuerdo con las características homogéneas, tales como relieve, litología y ambientes morfogenéticos.
El resultado, fue un mapa geomorfológico en el cual se pueden apreciar los ambientes
geomorfológicos que dan lugar a las diferentes unidades geomorfológicas.
A continuación, se describen las características geomorfológicas que hacen parte del área de estudio,
enmarcadas dentro de los conceptos de morfoestructura, morfogénesis, morfometría y
morfodinámica.
Página 14 de 156
• Geomorfoestructura
El área de estudio, se encuentra incluida dentro del Sistema de Orogénico Andino (Peneplanaciones
de la Amazonia).
• Provincia geomorfológica
Corresponde a un conjunto de regiones con geoformas parecidas y definidas por un macro relieve y
una génesis geológica similar. Localmente se relacionan con las regiones naturales y con los terrenos
geológicos de Colombia los cuales están delimitados por el trazo de megafracturas y suturas definidas
o inferidas (Carvajal, 2012).
Según este mismo autor, se determina con base en la interpretación de los procesos geomorfológicos
registrados (origen tanto endógeno como exógeno), que dieron lugar a la formación, evolución y
modificación de las mismas.
Para el área de estudio, se determinó solo dos ambientes morfogenéticos o regiones geomorfológicas:
Ambiente de origen denudacional (D): Determinado por la actividad de los procesos erosivos
hídricos y pluviales, y principalmente producto de procesos de meteorización, erosión y remoción en
masa, sobre geoformas pre existentes.
Página 15 de 156
Ambiente de origen fluvial (F): Corresponde a las geoformas generadas por procesos (erosión –
sedimentación), generadas por corrientes de agua tales como caños, arroyos, ríos, lagos o lagunas.
Esta categoría, se propone como la unidad básica de la cartografía geomorfológica. El término unidad
geomorfológica, se define como una geoforma individual genéticamente homogénea; mientras que la
subunidad, es una subdivisión de la unidad geomorfológica. Esta última, se encuentra determinada
fundamentalmente por los contrastes morfológicos y morfométricos, que relacionan el tipo de material
o la disposición estructural de éstos con la correspondiente topografía del terreno.
Por último, el componente o elemento geomorfológico está establecido por los rasgos del relieve
determinados por la morfometría detallada del terreno (Formas de ladera, escarpes estructurales,
pendientes anómalas, perfiles de suelos residuales o remanentes, escarpes naturales o antrópicos,
relieves internos de laderas o flancos, crestas, formas del valle) (Carvajal, 2012).
Este mismo concepto, lo sustenta Robertson y Otros (2013), mencionando que el componente o
elemento geomorfológico corresponde al máximo nivel de detalle asociado a escarpes o cambios de
pendiente (Taludes, orillares, escarpes de falla).
Con base tanto en los objetivos del estudio como en lo descrito anteriormente, la tabla siguiente
muestra la jerarquización realizada a partir tanto de la interpretación de los sensores remotos como
del reconocimiento de campo, para el área de estudio del río San Juan se presenta la siguiente
jerarquización geomorfológica (Tabla 3).
Página 16 de 156
GEOMORFO- PROVINCIA REGIÓN UNIDAD SUBUNIDAD COMPONENTE SÍMBOLO
ESTRUCTUR (AMBIENTE
A MORFOGENÉTICO)
Plano de De longitud moderada y FTFpt3
terraza T3 de forma recta a
cóncava
Fuente: Elaboración propia, 2022
Esta geoforma de origen fluvial, corresponde al cauce permanente sobre la cual se desarrolla
actualmente el tránsito del río San Juan, y su tributario inmediato el río Vides, así como los afluentes
identificados en la cartografía. La base de conformación corresponde a arenas bien lavadas, con
alguna presencia de arcillas no consolidadas, y otro estrato prominente son las arcillolitas del grupo
Orito. Este cauce ha sido permanente por lo menos los últimos 60 años, por como se puede datar en
las imágenes producto de sensores remotos.
Están constituidas por materiales clásticos gruesos como: limos y arcillas y en menor proporción
gravas o arenas, la disminución del tamaño de grano varía a lo largo del curso de la corriente y
lateralmente a medida que se aleja del cauce activo actual.
Los niveles más altos de terrazas, son los más antiguos y normalmente contienen los suelos más
evolucionados, a menos que hubiesen sido disectados. Cada nivel de terraza está separado de otros
por escarpes verticales a subverticales en los que afloran las capas de sedimentación.
El plano de terraza intermedio de comporta como un remanente erosivo de una terraza superior, que
fue lavada por el despojo de coberturas y socavación general dada por la migración de los cauces de
anidamiento.
Las terrazas más antiguas se localizan en las zonas donde actualmente la presencia de corrientes
activas y cauces abandonados están más alejados, por tanto, su moldeado solo se verá modificado
de manera somera por las aguas de escorrentía e infiltración.
Página 17 de 156
2.2.3.3. Escarpe de erosión (Deeme)
Corresponde a la forma de terreno convexa que se erige desde las terrazas más recientes hasta las
zonas donde se conforman morfoestructuras de origen estructural o hasta el recambio de hoyas
hidrográficas. Aunque la mayor parte poseen procesos de erosión activos, que particularmente
agregan clastos a terrazas inferiores, algunos son de garantía estable, y su agregación se ofrece por
bloques desprendibles de lodolitas.
Página 18 de 156
2.2.4. Morfodinámica del área de estudio
De acuerdo con Robertson, Jaramillo, Castiblanco, & Instituto de Hidrología Meteorología y Estudios
Ambientales (2013), la morfodinámica, es la parte de la geomorfología que trata de los procesos
geodinámicos externos (principalmente denudativos), tanto antiguos como recientes que han
modelado y continúan modelando el relieve y son los responsables del estado actual de las geoformas
o unidades de terreno. Consiste en determinar el tipo de proceso y su intensidad en la conformación
de las geoformas (procesos erosivos y movimientos de remoción en masa.
De manera general, la tabla a continuación, resume los principales agentes geomorfológicos, agentes
de la meteorización y fuerzas de desplazamiento, así como los procesos exógenos a que dan lugar.
El área específica del proyecto posee en la margen derecha procesos activos de erosión fluvial,
acompañados de erosión pluvial dada por la degradación antrópica presente en el intersticio de la
obra. Si bien los procesos naturales citados son de carácter general y genérico en todas las corrientes,
su aumento en la magnitud está completamente direccionada y relacionada con la alteración del
componente ecosistémico forestal del lugar.
Corresponde a las diversas formas de cómo los agentes móviles (agua, vientos, glaciares, etc.),
desprenden y transportan los productos de la meteorización y de la sedimentación (Villota & Instituto
Geográfico Agustin Codazzi (IGAC), 2005).
Página 19 de 156
El transporte de dichas partículas, se presenta en forma de erosión geológica; erosión acelerada
(antrópica); erosión por escurrimiento (laminar, en surcos y en cárcavas); erosión fluvial; erosión eólica
o erosión diferencial.
Interpretación sensores remotos: Con base en los sensores remotos, la identificación de los
procesos erosivos, se realizó a partir de realces, refinamientos, contrates de color y todas aquellas
combinaciones que se pudieron obtener de las herramientas del Sistema de Información Geográfica
-SIG-. Es importante mencionar, que los diferentes tratamientos realizados por medio del SIG,
dependieron tanto de la escala de los sensores remotos como de la resolución de los mismos.
• Erosión laminar
Corresponde a la remoción uniforme de una lámina delgada de suelo de una superficie inclinada, sin
que se formen claramente canales de desagüe. Por lo general, se detecta cuando el horizonte A del
suelo es removido totalmente, aflorando el horizonte B o el C.
• Erosión en surcos
Se desarrolla a partir de la erosión laminar; donde, la remoción del suelo ocurre en mayor cantidad a
lo largo de pequeños canales, que se han generado tanto por la actividad antrópica (uso del suelo y
ganado) como por variables naturales (pendiente y clima).
• Erosión en cárcavas
Corresponde a la erosión por escurrimiento muy avanzado; donde, se generan surcos profundos y
amplios, los cuales están regulados tanto en forma como en tamaño, por las características del suelo,
relieve, clima, cobertura vegetal y uso actual del suelo.
La morfología fluvial en el proyecto se comporta de manera activa en todas las modalidades de pérdida
de material originario de los horizontes primarios; en principio, podrá notarse que las terrazas medias
y altas presentan erosión laminar horizontal con conclusiones por cárcavas para entrega final en el
Página 20 de 156
cauce permanente. La erosión laminar está definida en las vegas de divagación, y son el producto del
lavado permanente de estas zonas en épocas en las cuales la fuente supera los umbrales del canal.
La presencia de estas huellas erosivas sugiere que la obra de paso deba contar con plataformas de
desarrollo-acceso que garanticen el aislamiento o contacto de la superficie afectada por laminación
con arranque de tablero. Por esa consecuencia, se deberá además desarrollar obras de drenaje
complementarias que limiten la acción de laminación, preferiblemente en modo espina de pescado,
sea superficial o subsuperficial. Por parte de caminos de cruce en el puente, se deberá procurar que
estos sean de carácter peatonal únicamente, pues el terraceo por bestias o ganado, genera alta
presión de poros que traduce a formar líneas de socavación.
Página 21 de 156
En materia de los cultivos, se debe procurar que los drenajes sean de carácter lateral en los colectores,
y sembrar a contra-pendiente para apoyar la sujeción de terreno por raíces.
Hace referencia a todo movimiento ladera abajo del material de una pendiente bajo la influencia de la
gravedad, usualmente lubricado por agua (lluvia, agua subterránea, deshielo). El movimiento puede
ser muy lento, lento o rápido.
El procedimiento que se tuvo en cuenta para llevar a cabo la identificación y caracterización de los
procesos morfodinámicos que predominan en la zona de estudio, corresponden a:
Interpretación sensores remotos: Con base en los sensores remotos, la identificación de los
procesos morfodinámicos, se realizó a partir de realces, refinamientos, contrates de color y todas
aquellas combinaciones obtenidas de las herramientas del Sistema de Información Geográfica -SIG.
Al interior del área de estudio, se identificaron movimientos en masa relacionados con el “terraceo o
pate vaca”; las cuales son clasificadas por (Villota & Instituto Geográfico Agustin Codazzi (IGAC),
2005), como un flujo lento de suelo en estado plástico; deslizamientos activos e inactivos con
reptación.
Este movimiento del terreno, se define como un desplazamiento lento; donde actúa la gravedad, el
agua del suelo y las sobrecargas o presiones externas. El suelo se comporta como un flujo en estado
plástico (Villota, 2005). A partir de estas sobrecargas (ganadería), el suelo se deforma presentado las
“terracetas” o “pate vaca”. Este proceso, conlleva a un desgarramiento de la parte alta de los suelos
de la ladera, lo que genera un deslizamiento y, posteriormente, la generación de procesos erosivos,
tales como surcos y cárcavas.
Página 22 de 156
Es importante mencionar, que el criterio utilizado para diferenciar un deslizamiento activo de uno
inactivo, corresponde a la ausencia o presencia de vegetación, respectivamente.
• Desplomes
Corresponde a los procesos asociados a la dinámica fluvial, la cual genera fuerzas laterales sobre los
taludes de las márgenes o riberas. Generalmente, se presenta si la pendiente de las márgenes es
vertical a subvertical (escarpes de terraza); en donde, la caída del material de suelo es libre, o sea sin
contacto con la superficie que hace parte del talud.
Es importante mencionar, que, debido a la escala de trabajo, no todos los procesos morfodinámicos
han sido cartografiados. La Figura 7 representa los escasos y leves procesos que posee el cauce,
principalmente como deslizamiento vertical activo sobre estratos arenosos con bases arcillosas y
limosas de limitada cohesión. De igual forma, se identifica una línea de desplome vertical en la zona
inferior de donde quedaría el apoyo izquierdo (sugerido) de la infraestructura. Ese desplome trata de
los procesos morfodinámicos dados por las presiones de la roca sobre el compuesto inferior de la
lodolita, lo que provoca que se desarrollen desplomes en bloques y se acumulen en la rivera para dar
de cuenta la evidencia del fenómeno erosivo.
El terraceo, más que una consecuencia natural, se trata de una consecuencia antrópica del desarrollo
del sector, pues la población es concentrada en menos de 2 hectáreas, con dinámicas amplias y
densas, lo que provoca que el paso y movimiento desarrolle cambios significantes en la tipología de
terreno.
Los deslizamientos son subverticales por la pérdida progresiva del suelo de soporte por la corriente
del río, fenómeno completamente natural, e incrementado por la modificación de las coberturas de
galería -como se indicó arriba-.
Página 23 de 156
Figura 7. Procesos de movimientos en masa en el área de proyecto
Los movimientos en masa por desplome sugieren una evaluación puntual de la estabilidad en lo que
confiere al desplome de bloques de lodolita, pues si bien en el lugar se puede observar la continuidad
de la masa de roca, sus bases por la presión dada llegan a fracturarse y desintegrarse.
Para García et al. (2011) los ríos aluviales de llanura son sistemas de gran complejidad geomorfológica
y ambiental que destacan por sus activos procesos de dinámica fluvial. Su ajuste permanente y las
fluctuaciones de caudal líquido y sólido se traducen en movilidad lateral y vertical (Werritty, 1997). Sin
embargo, diversas actividades humanas sobre la cuenca, la llanura de inundación y el cauce modifican
Página 24 de 156
esta movilidad, que es la que garantiza la riqueza y diversidad de esos sistemas naturales (Malavoi et
al., 1998).
Para el desarrollo del componente de dinámica fluvial se ha sustentado con metodologías de análisis
vectorial de la cartografía fluvial a través de los SIG. El uso de los Sistemas de Información Geográfica
(SIG) supone un avance en los estudios de dinámica fluvial, facilitando el proceso de análisis y
aportando datos cuantitativos sobre la variación temporal y espacial de los cursos fluviales. En este
contexto, a nivel internacional se han realizado varios acercamientos a esta temática (Lagasse, Spitz,
Zevenbergen, & Zachmann, 2004).
El río San Juan indica que los patrones de la evolución temporal de la dinámica geomorfológica en el
tramo estudiado son muy similares a lo largo del tiempo en los ingresos a las zonas de control. Esta
evolución está estrechamente ligada a las perturbaciones que sufren desde la última mitad del siglo
XX por la intervención antrópica, las cuales han desencadenado una intensa simplificación de los
cauces. En este sentido, estos cursos medios y altos se caracterizan por reducir sus divagaciones y
estrechar o dos corrientes permanentes por varios años.
Para conocer la evolución de los cursos fluviales se ha realizado un análisis diacrónico mediante
ortoimágenes e imágenes satelitales de diferentes fuentes. Se han seleccionado y utilizado fotogramas
aéreos verticales de las siguientes fechas: 1976 1976/85, 1985/99, 1999/08 y 2018/22. La ortoimagen
del año 2022 se obtuvo mediante vuelo no tripulado con RPAS (Drone) fue referenciada mediante el
sistema de Referencia MAGNA Sirgas origen Oeste, mientras que las otras imágenes fueron obtenidas
en UTM 18N, se proyectaron al sistema de referencia de la ortoimagen.
El análisis de la dinámica fluvial se ha llevado a cabo mediante dos metodologías diferentes, de cara
a reforzar los resultados. Ambas metodologías se realizan mediante el uso de técnicas y herramientas
de SIG vectorial.
La primera metodología tiene como objetivo conocer la movilización general de los cauces en cada
período de estudio. Se basa en el cálculo de la superficie movilizada, que es la superficie modificada
entre la orilla del cauce en la fotografía de inicio del período y la misma orilla del cauce en la fotografía
final del período (Downward & Brookes, 1994). También se calculó la superficie activa, que incluye el
trazado del cauce de cada par de años analizado y la superficie que queda entre ambos trazados. Esta
superficie resulta muy útil para conocer cuantitativamente el nivel de reducción de la superficie de
cambio a lo largo del tiempo.
De forma paralela se realizó un análisis de la migración del meandro donde se localizará la obra, así
como de la tendencia de su futura evolución a partir de la metodología de Lagasse et al. (2004). Esta
metodología se basa en la edición del círculo circunscrito sobre la orilla cóncava de cada meandro en
los diferentes años estudiados. Posteriormente se han de editar los centroides de cada círculo, que
son la referencia para calcular el radio de curvatura en cada año y la variación de la distancia y del
ángulo de movimiento entre diferentes años, como ilustra la figura a continuación. Con los datos
Página 25 de 156
obtenidos y según la metodología propuesta por Lagasse et al. (2004), se puede calcular una
predicción o tendencia de la migración de los meandros, teniendo como referencia la variación de los
parámetros (ángulo, distancia y radio) en el último período.
El análisis integral muestra de manera general que no existe mayores cambios en la geometría del
cauce en los diferentes tiempos a nivel de superficie y parámetros de esbeltez, lo que pone de
manifiesto la estabilidad de planta del cauce, es decir; el espejo de agua se mantiene en medida que
la socavación lateral aumenta hasta encontrarse con la lodolita y control natural de la formación orito,
a partir de la cual ya es completamente improbable que evolucione. En ese sentido, la superficie
hidráulica de movilización no estaría netamente influenciada por sus áreas, sino por la disposición de
las mismas, sus ángulos y radio del meandro en el ataque.
Página 26 de 156
Figura 9. Dinámica traslacional de la zona activa entre 1976 y 2022
Página 27 de 156
Figura 10. Modos de desarrollo del bucle del meandro.
a) extensión, b) traslación, c) rotación, d) conversión a un bucle compuesto, e) corte del cuello por cierre, f) corte
diagonal por vertedero y g) corte del cuello por vertedero.
Fuente: (Lagasse et al., 2004)
Por su parte, la rata de cambio del radio de curvatura tampoco sufre incrementos exponenciales de
reducción, con comportamiento generalizado, donde se traduce todo a un fenómeno de consolidación
de corte de cuello por vertedero, que se sostiene por la desembocadura del río Vides y por el control
geológico de la formación Orito.
Tabla 7. Variación de la distancia de cambio (m) y el ángulo (º) en el área del proyecto
Tiempos Orientación Azimut Longitud (m) Rata (m/año) Desarrollo
1976 - 1985 67 160 17.7 Extensión
Corte del cuello
1985 - 1999 63 51 3.6
por vertedero
Corte del cuello
1999 - 2008 83 90 10
por vertedero
Página 28 de 156
Tiempos Orientación Azimut Longitud (m) Rata (m/año) Desarrollo
Corte del cuello
2008 - 2022 249 26 1.8
por vertedero
Fuente: Elaboración propia, 2022
Figura 11. Variación de la distancia de cambio (m) y el ángulo (º) en el área del proyecto
Página 29 de 156
2.2.6. Patrón morfológico actual
El rasgo más característico de una corriente joven es la actividad de esa corriente y la rápida erosión
de su cauce. En esta etapa, las laderas del valle llegan directamente (o casi) a las orillas del río y
adoptan un perfil transversal en forma de V. No existe planicie de inundación, o sólo una muy angosta.
El gradiente de la corriente es pronunciado, marcado por cascadas y rápidos y los tributarios tienden
a ser pocos y pequeños.
La juventud de un río corresponde a la parte más alta de la cuenca hidrográfica en donde se originan
el caudal y los sedimentos. Está caracterizada por tener fuertes pendientes, velocidades altas y
caudales bajos. El cauce transcurre por relieves escarpados y estratos rocosos principalmente. La
energía del río se consume en profundizar el cauce.
En un río maduro, la velocidad de corte hacia abajo es más lenta; el gradiente más suave y la mayoría
de las cascadas y rápidos han sido eliminados. A medida que avanza el ensanchamiento del valle con
mayor rapidez que su profundización, comienza a formarse una planicie de inundación; y la corriente
comienza a divagar o a hacer meandros a través de la planicie. En esta etapa ha alcanzado el valle su
mayor profundidad. Es la de transferencia o transporte de agua y sedimentos de la zona de montaña
a la zona baja. La energía del río se consume en profundizar y ampliar el cauce. El río forma meandros
y trenzamientos.
En la etapa de vejez, el ensanchamiento del valle, aunque lento, todavía predomina sobre el corte
hacia abajo y la planicie de inundación es más ancha que la faja de meandros. Los lagos en “media
luna”, los meandros abandonados y los bordos naturales son ahora más comunes que en la madurez.
Corresponde a la parte baja en donde el sedimento se deposita. Se caracteriza por tener pendientes
bajas, velocidades bajas y altos caudales. El cauce transcurre en estratos aluviales de gran espesor.
La tendencia del cauce es a ampliarse. El patrón que representa a la corriente se comporta de manera
adolescente, donde la variabilidad se influye por la cesión de zonas libres y fácilmente erosionables.
Página 30 de 156
2.2.6.1. Patrones que afectan el paisaje fluvial
Existen factores a nivel de cuenca, de valle y de cauce, como son los geomorfológicos, geológicos,
sismológicos, suelo, uso y cobertura de suelo, sedimentos, hidrológicos, hidráulicos, eco-sistémicos y
el nivel base, que afectan la estabilidad del cauce. En forma general, se hace a continuación una breve
descripción de los factores que de una u otra manera afectan la dinámica de un río.
• Factores ecológicos
Los depósitos superficiales existentes y las clases de rocas con sus características de meteorización,
erosión, discontinuidades, rasgos topográficos y drenaje, influyen sobre el paisaje fluvial. Es
importante conocer y definir los contactos litológicos y las fallas geológicas existentes en la cuenca
del río pues pueden inducir movimientos en masa de las laderas que al caer a un cauce crean
eventualmente represamientos del flujo con consecuencias de magnitud impredecible hacia aguas
abajo.
Figura 12. Evolución de un canal.
Incisión inicial (a, b), pasando por el ensanchamiento (c, d), hasta la agradación (d, e) y eventual estabilidad. Fuente:
Lagasse, O. F., Schall J. D., Richardson, E. V. (2001)., citado por (Universidad del Cauca, 2006)
En el área de estudio se logró encontrar prácticamente todos los patrones de evolución local, siendo
como particular la referencia que indica que el puente actual se sitúa a la izquierda en el
ensanchamiento homogéneo activo, mientras que a la izquierda se encuentra en un proceso de
Página 31 de 156
ensanchamiento heterogéneo estable. La obra proyectada dentro de estas dinámicas, lo cual ofrece
estabilidad por la misma razón que el ensanchamiento lateral izquierdo ya no es posible que
evolucione.
Agradación heterogénea activa – Fahta, Agradación heterogénea estable – Fahte, Agradación homogénea activa –
Faha, Ensanchamiento heterogéneo activo – Fehta, Ensanchamiento homogéneo estable – Fehe
Fuente: Elaboración propia, 2022
• Factores sísmicos
La sismicidad de una zona se considera un factor detonante que induce por ejemplo movimientos en
masa que pueden llegar a afectar la dinámica fluvial de un cauce. Es importante tener en cuenta
aspectos como: datos históricos e instrumentales de localización y magnitud y frecuencia o período
de retorno de los sismos.
Página 32 de 156
En el plan municipal de gestión de riesgo de desastres del municipio de Orito no se reporta
deslizamientos o efectos desencadenantes de deslizamientos, de hecho, se reconoce que en el río no
se desarrolla inundaciones de tipo torrencial, su comportamiento es laminar de largo alcance en zona
de llanura.
Fuente: (Servicio Geológico Colombiano, 2022) 1 El círculo rojo indica el área de estudio.
http://sish.sgc.gov.co/visor/sesionServlet?metodo=irAIntensidadesTodas&idDepartamento=&idMunicipio=&cuadranteXMin=&cuadra
nteXMax=&cuadranteYMin=&cuadranteYMax=
Página 33 de 156
En ese sentido, sobre las márgenes en el área del proyecto se encuentran suelos de características
asociadas a la unidad cartografiable VUAai (lado derecho), la cual en el Estudios General de Suelos
del Departamento del Putumayo (IGAC, 2014), referencia que corresponde a la subclase 5hs-U, que
presenta condiciones de suelos que se han desarrollado a partir de depósitos aluviales heterogéneos;
son superficiales a profundos, bien y pobremente drenados, texturas moderadamente finas y
moderadamente gruesas, reacción moderada y fuertemente ácida, alta a baja saturación de aluminio,
fertilidad baja a alta. Estas tierras para uso y manejo presentan limitaciones por inundaciones
frecuentes, exceso de humedad debido a las altas precipitaciones, drenaje deficiente sectorizado, alta
saturación de aluminio en algunos sectores y fertilidad baja, mientras que al lado izquierdo se
encuentra la unidad LUBd2, con la subclase 7pec-U, que representa suelos con material parental de
limolitas, lodolitas y arcillolitas. Los suelos son profundos a, bien drenados, de texturas muy finas y
finas, fuertemente ácidos, con alta saturación de aluminio, de fertilidad baja. Los suelos están
afectados por erosión y otros por movimientos en masa tipo pata de vaca.
El uso y cobertura del suelo es el factor que más fácilmente da cuenta del grado de intervención
antrópica sobre el medio y, por consiguiente, permite explicar algunos de los cambios que se
presentan en los ríos.
La variación del uso y cobertura puede en algunas áreas favorecer los procesos de inundación,
degradación, agradación y movilidad de los cauces y en otras, propender por la recuperación de áreas
degradadas. Sin embargo, la revegetación de antiguos cauces o madre viejas, el secamiento de
planicies de inundación o la consolidación de centros poblados en áreas de antigua divagación de un
río, pueden llegar a generar un problema puesto que para una creciente.
En el área de estudio se pudo observar que el uso del suelo principal es la agricultura por unidades
sembradas en hileras, existen pastos limpios y pasturas enmalezadas. No obstante, se puede observar
que existe un relicto de bosques de galería en la zona posterior a la obra, estos se componen por
arbustos y árboles que no superan los 5 metros de altura, con madurez muy baja, lo que habla de un
estado de preservación no sustentable, esto por dos razones, una se debe a los caudales
extraordinarios de los procesos fluviales y la frecuencia con la que los eventos extraordinarios tienen
lugar (Ballesteros et al., 2017), y otro trata del proceso acelerado de intervención de vegas agrícolas
sin proveer una barrera de sustento ante degradación del uso del suelo.
• Factores geomorfológicos
Página 34 de 156
El estudio de las formas del sistema fluvial es indispensable para conocer e identificar los controles
naturales que ofrece la cuenca al paso de una creciente. Dichos controles pueden tener un carácter
permanente como por ejemplo las terrazas aluviales (llanuras de inundación abandonadas) o pueden
tener un carácter transitorio como son los complejos de orillares, los cuales pueden ser modificados
con el paso de una creciente extrema.
En la corriente del sector se logra reconocer pocas formas específicas que se puedan considerar
como formadoras o actuadoras de la corriente, caso de ello es el rápido que se identifica en todos los
sets de imágenes citados en el estudio. El factor formador más determinante es el complejo orillar de
desembocadura que se localiza a 260 m del ponteadero, que en trasfondo permite garantizar una
superficie estable para mantener la formación del paso del puente.
• Sedimentos
Las características del material que transporta una corriente están asociadas con la configuración
geológica de la cuenca de drenaje. Los cauces son usualmente susceptibles de ser degradados, al
presentarse un desequilibrio entre la capacidad de transporte de sedimentos y el suministro de los
mismos, alterando la geometría hidráulica de la sección transversal. Sin embargo, para simplificar la
modelación hidráulica, en muchos casos se recurre a considerar la geometría constante en el tiempo.
La interacción entre el flujo y el material granular aluvial ha sido más ampliamente estudiada debido a
que es el caso más frecuente asociado con problemas en la hidráulica de ríos.
Página 35 de 156
Figura 15. Interacción geomorfológica del cauce de socavación.
• Factores hidráulicos
Los factores hidráulicos involucrados en la estabilidad del cauce son numerosos e incluyen magnitud
y frecuencia de las crecientes, velocidades del flujo, niveles del agua, gradiente hidráulico, formas de
lecho, transporte de sedimentos, resistencia al flujo, tipo de flujo, contracciones al flujo, basuras y
cuerpos flotantes.
Página 36 de 156
En los factores hidráulicos que se destacan, está la formación de dunas de contracción de corriente y
desvío de masas, es decir, barras sumergidas que actúan aumentando la velocidad de la corriente
sobre la margen de menor resistencia topográfica, generando nuevas barras en la sección más lenta
y socavando en la otra sección de la superficie, en la microcuenca no hay cambio de barras, y
actualmente se comienza un proceso de formación por agradación cóncava y socavación convexa
precisamente antes y posterior al puente. Dicho esto, se trata de un elemento deposicional que deberá
intervenirse constantemente o de lo contrario deberá diseñarse con lecho estable.
La primera relación empírica de interés es en qué condiciones un río forma un cauce trenzado o un
cauce único meandriforme. En la relación i·Q0.44 donde i es la pendiente en (m/m), Q el caudal en
(m³/seg) se puede establecer una frontera entre una y otra morfología (Martín Vide, 2003). Así,
De acuerdo con lo anterior, teniendo en cuenta el caudal medio de la fuente determinado por
(Calvache-Fajardo & Cuesvas, 2016) es de 723.2 m3/s, la relación u expresión adquiere un valor de
0.18, hablaríamos de un cauce trenzado.
Página 37 de 156
No obstante, de acuerdo con los autores de la cita matemática, la expresión no es determinante para
establecer la forma del cauce de manera general, porque es muy frecuente que un río sea trenzado
en su tramo alto y luego pase a ser meandriforme, o en la mayoría de los casos de llanura, que estos
ríos sean trenzados y sinuosos, la cual sería la forma más apropiada de reconocer el cauce de la
quebrada.
En los ríos con meandros se aceptan como significativas, pese a la lógica dispersión de los datos de
campo, algunas relaciones empíricas entre la longitud de onda de la forma λ (Figura 16), su amplitud
α, el caudal Q y la anchura del cauce en la superficie libre β. La longitud λ es de 7 a 11 veces la
anchura del río β (la proporción más común es 10; entre λ y β se supone una relación potencial pero
el exponente de β resulta prácticamente la unidad) y la amplitud α es unas tres veces la anchura β,
aunque esta última relación tiene menos fundamento, debido a que cambia más durante la evolución
de un meandro, mientras que por el contrario λ puede permaneces invariable.
De acuerdo con las relaciones anteriores se tiene que el río en la zona de estudio posee una longitud
de onda de 610 m, el ancho del cauce es 103.2 m, lo cual establece una relación de 4 a 1, mientras
que la razón de amplitud 1309 m, es 10 veces el ancho del cauce. Por lo anterior se concluye que se
trata de un cauce meándrico en proceso de transformación desde trenzas residuales, del mismo, que
la fuente posee condiciones de desarrollo bastante homogéneas que indican un río con una madurez
segura.
Página 38 de 156
Figura 17. Parámetros de sinuosidad del cauce
Por otro lado, el (Ministerio para la transición Ecológica, 2019) propone el cálculo del índice de
sinuosidad teniendo en cuenta solo la longitud real de la corriente medida desde el eje del espejo,
Página 39 de 156
sobre la longitud axial (línea recta que une el inicio con el final del tramo) de la corriente misma. El
término matemático es el siguiente:
𝐿𝑟𝑒𝑎𝑙 Ec. 2
𝐼𝑠𝑖𝑛𝑢𝑜𝑠𝑖𝑑𝑎𝑑 =
𝐿𝑖−𝑓
Dada la expresión anterior, el índice de sinuosidad sería de 1.32, teniendo en cuenta que la longitud
real es de 1730 m, sobre el valor de λ que ya se trató con la otra metodología.
Los umbrales propuestos pueden modificarse levemente, si se observa la existencia de otros atributos
morfológicos que indiquen con claridad la vinculación de la masa o tramo con un determinado tipo
morfológico.
En materia de definición de este patrón de drenaje muchos autores han desarrollado mecanismo de
clasificaciones de ríos, que puede dividirse a grandes rasgos en dos enfoques: 1) relaciones
cuantitativas (que pueden ser empíricas o teóricas) y 2) marcos conceptuales (Kondolf, Montgomery,
Piégay, & Schmitt, 2005).
Relaciones cuantitativas: En los estudios de (Lane, 1957; Leopold, Wolman, & Miller, 1964) se observó
que para una descarga determinada, los canales trenzados se producen en pendientes más
pronunciadas que los ríos serpenteantes. Los estudios reconocieron un continuo del patrón del canal,
pero Leopold y Wolman (1957) propusieron un umbral entre los ríos serpenteantes y los ríos trenzados
(Figura a continuación), proporcionando un medio para predecir los cambios en el patrón del canal
en función de la alteración de la descarga o de la pendiente del canal. Ambos estudios también
reconocieron que hay factores adicionales que afectan el patrón del canal, como el tamaño de los
granos, la carga de sedimentos, la vegetación ribereña, la rugosidad del canal, el ancho y la
profundidad. Los investigadores posteriores modificaron el marco del esquema para incluir el tamaño
de los granos (que altera la ubicación de los límites entre los diferentes patrones de los canales) y
Página 40 de 156
para distinguir los canales anastomosados y los canales errantes; los ríos errantes son de transición
entre morfologías serpenteantes y trenzadas, mientras que los ríos anastomosados son canales
multihilos separados por islas cortadas de la planicie de inundación y se distinguen de los canales
trenzados que se forman por deposición de barras y división del flujo en el canal sub-secuencial.
El patrón del canal (serpenteante, recto, trenzado) como una función de la pendiente del canal y la descarga del
banco. Fuente: Leopold et al. (1964)
Desde este análisis ubican al río en una transición de ríos trenzados a serpenteantes, lo cual agrega
más al criterio de clasificación ya definido.
Marcos conceptuales: De acuerdo con Schumm (1960, 1963a) los trabajos sobre ríos de lecho de
arena y grava en las grandes llanuras del oeste de los Estados Unidos destacaron que el patrón y la
estabilidad del cauce están fuertemente influenciados por la carga impuesta por el río (tamaño del
sedimento y modo de transporte) y el contenido de arcilla-limo de la llanura de inundación (que
proporciona la cohesión necesaria para el desarrollo de los meandros del río). Basándose en estas
observaciones, Schumm (1963b, 1985) planteó un marco conceptual para clasificar los ríos aluviales
que relacionaba el patrón del cauce y la estabilidad con (1) el contenido de arcilla-limo de los bancos,
(2) el modo de transporte de los sedimentos (carga suspendida, carga mixta, carga del lecho), (3) la
relación entre la carga del lecho y la carga total (una función de la potencia de la corriente, el tamaño
de los sedimentos y el suministro), y (4) la pendiente y la relación anchura-profundidad del cauce
Página 41 de 156
(figura seguida). Estudios posteriores señalaron el papel de la vegetación ribereña y la fuerza de las
raíces para afectar la cohesión de los bancos, el ancho del cauce y el patrón del cauce.
Bajo los conceptos anteriores, la corriente se clasificaría localmente en la categoría 3b para la zona
de estudio, es decir, Meandriforme por gradiente de la relación ancho-profundo, estabilidad relativa,
del mismo tenor es para el coeficiente de carga (15% apróx), tamaño de sedimento, carga de
sedimentos, velocidad de flujo y energía de contacto.
Las unidades de canal son unidades morfológicas de escala inferior al alcance (por ejemplo, diferentes
tipos de piscinas, barras, escalones, canales de rifle) que forman los componentes básicos de las
morfologías de mayor alcance, como las piscinas escalonadas o los canales de rifle de piscina. Bisson
y otros (1982) elaboraron una clasificación descriptiva y detallada de las unidades del canal en los
arroyos del noroeste del Pacífico para cuantificar los diferentes tipos de hábitat físico de los
salmónidos. Trabajos posteriores realizados por otros examinaron la hidráulica de las unidades de
canal, su respuesta a la tala de madera y a la eliminación de grandes desechos leñosos, y las
Página 42 de 156
características físicas y biológicas de las unidades de canal en los arroyos empinados. Además, se
han elaborado clasificaciones jerárquicas de las unidades de canal y su hidráulica, pero estos
enfoques utilizan principalmente descripciones cualitativas de la corriente (agua rápida frente a agua
lenta) y de la rugosidad de la superficie del agua ("turbulenta" frente a "no turbulenta") (Kondolf et al.,
2005).
Bisson y otros (2006) recomienda que estos últimos términos sean reemplazados por flujo "áspero"
versus flujo "suave", ya que la mayoría de los flujos de los ríos son muy turbulentos, sensu stricto. La
mayoría de las clasificaciones de unidades de canal se centran en el canal húmedo, generalmente
excluyendo las barras, pero una clasificación detallada de los tipos de barras y los procesos físicos
asociados ha sido desarrollada por Church y Jones (1982). Las figuras asociadas a las unidades de
canal son las siguientes:
Plunge Pool-PP: Piscina de inmersión, Lateral Scour Pool-LSP: Piscina de socavación lateral, Backwater Pool-BP:
Piscina de remanso, Dammed Pool-DP: Piscina represada, Glide-G: Deslizamiento, Riffle-R: Regata, Rapid-Ra: Rápido,
Falls-F: Caida, Cascade-C: Cascada, Side Channel-SC: Canal lateral. Fuente: (Frissell, Liss, Warren, & Hurley, 2014)
Página 43 de 156
Con la clasificación anterior se desarrolló el reconocimiento de unidades de canal en el tramo de
estudio, siendo reconocida la piscina de socavación lateral, un rápido antecesor y un remanso de
salida que por topografía se identifica con profundidades de 2 a 6 metros.
2.2.9. Caracterización inicial de los cambios de tipo morfológico en las últimas décadas
Teniendo en cuenta los tiempos comprendidos entre el año 2014 y 2022, sobre los tramos
hidromorfológicos presentes en la masa de agua, tanto en la situación de no alteración o natural como
en la actual o modificada se procederá a caracterizar el posible origen de los cambios detectados en
los últimos años. Estos cambios morfológicos responden a las distintas presiones existentes sobre las
masas de agua y son la consecuencia de las alteraciones hidrológicas y morfológicas de las variables
que determinan la dinámica fluvial (Ministerio para la transición Ecológica, 2019).
En caso de haber identificado cambios en el tipo morfológico se deberá definir el origen o causa de
los mismos como:
• Existencia de regulación aguas arriba (Or): la presencia de embalses modifica los volúmenes
de aportación y lamina las avenidas. Esto conlleva procesos de simplificación y estrechamiento
de los cauces.
• Obras civiles de movilidad (Ocm): Construcción de vías, caminos o adecuaciones por cable u
otros mecanismos – modos de transporte.
• Cambios de los usos de suelo en la cuenca vertiente (Cus): los cambios de los usos del suelo
modifican las condiciones de infiltración/escorrentía a partir de las precipitaciones. Esto altera
la hidrología superficial y en consecuencia los volúmenes y caudales circulantes por el cauce
receptor, lo que origina procesos de simplificación geomorfológica.
• Acciones directas sobre el cauce (Adc): cambios antrópicos sobre el trazado, sección
transversal o pendiente longitudinal de un cauce (desvíos, canalizaciones, acortamientos,
etc.).
Si bien en la zona de estudio no se han desarrollado acciones permanentes que influyan de manera
compleja en la modificación hidromorfológica de la corriente, si se logra vislumbrar desde la década
Página 44 de 156
de los 80’, algunas alteraciones paisajísticas que por manera directa no pueden ser masivas, pero por
el uso y manejo de las mismas pueden sucumbir en los cambios del curso del río.
En la corriente se puede notar como existen modificaciones en el cambio del uso de suelo, donde los
bosques de galería con uso de preservación fueron sustituidos por pastos y cultivos transitorios
aislados.
• Desviado (De): creación de un cauce nuevo. Se incluyen todas las acciones realizadas sobre
el cauce original que suponen un cambio total o parcial de su trazado en planta.
• Estrechado (Es): reducción artificial de la anchura del cauce activo, por eliminación de brazos,
barras/islas activas, o por el desarrollo de intervenciones con capacidad para reducir la
anchura.
• Canalizado (Ca): obra de encauzamiento que fija totalmente las márgenes. Se incluyen todas
las acciones llevadas a cabo sobre el cauce que supongan su canalización mediante la
regularización de su trazado, sección transversal o pendiente, con el objetivo de encauzar el
río. La canalización puede ser en tierras, sobre el material original del cauce, o bien mediante
la aportación de materiales externos (piedra, escollera, hormigón, etc.).
• Cubierto (Cu): cauce artificialmente enterrado en una longitud superior a los 200 metros, cuyas
secciones pueden presentar diferente geometría (circular, abovedada, etc.), que se ubica (en
su totalidad o sensiblemente) sobre el trazado del cauce original y al que reemplaza desde el
punto de vista de la conducción del flujo líquido y sólido. Está construido a partir de materiales
como el hormigón, el acero, el aluminio o el plástico de alta densidad.
Página 45 de 156
Para realizar este análisis se ha buscado sobre incidentes que permitan reconocer si en el municipio
en los últimos años se han desarrollado actividades de ese tipo, dando a lugar que no se han
contratado acciones de descolmatación, del mismo modo, los lugareños refieren que en el sitio se no
se adecúa ni se realiza ningún tipo de intervención para mantener la línea de la corriente.
De acuerdo con (Ministerio para la transición Ecológica, 2019) como complemento al estudio en
gabinete de la alteración del régimen hidrológico, tanto en caudales líquidos como sólidos, se
analizarán en los tramos hidromorfológicos de la masa de agua los movimientos verticales de la fuente
(dinámica vertical) que pueden ser nulos o escasos, o bien puede definirse por dos procesos
opuestos, el de incisión o encajamiento del lecho y el de acreción o elevación del lecho por
acumulación de sedimentos.
Es normal que un cauce confinado con pendiente alta asista a un proceso natural progresivo de
incisión. También es normal que en un cauce de llanura con muy poca pendiente los sedimentos se
acumulen y eleven ligeramente algunas zonas de depósito. La mayor parte de los cursos fluviales
asisten a una situación de equilibrio, con zonas de incisión y otras de acumulación y con procesos
muy lentos, apenas perceptibles a escala temporal humana. El problema surge cuando esta dinámica
vertical resulta acelerada, es decir, se manifiesta muy rápidamente acompañada de síntomas
evidentes: acumulaciones excesivas o extensas por acreción y descalzamiento de las márgenes o
desnudez de las raíces vegetales por incisión.
Figura 21. Efecto de incisión en el lecho como consecuencia del estrechamiento del cauce por ocupación de las
márgenes del río
Página 46 de 156
Para la caracterización de este componente se identificará en el área de estudio por cada sección,
que tipo de efecto se genera en el río, empezando por; Disminución del nivel freático (Dnf), Pérdida
de espacio de ribera (Per) e Incisión del lecho (Ile).
Disminución del nivel freático (Dnf), Pérdida de espacio de ribera (Per) e Incisión del lecho (Ile)
Fuente: Elaboración propia, 2022
Página 47 de 156
2.2.12. Estructura y sustrato del lecho
La estructura y sustrato del lecho son elementos esenciales en el ecosistema fluvial, puesto que de
ellos depende el desarrollo de una gran parte de los ciclos biológicos asociados al cauce.
Esta se caracterizará de acuerdo con (Ministerio para la transición Ecológica, 2019), a través de los
trabajos de campo realizados en cada uno de los subtramos correspondientes a los tramos
hidromorfológicos en que ha quedado dividida la masa de agua en estudio.
• En roca: el cauce está labrado en lecho rocoso, aunque puede haber algunos sedimentos.
• Coluvial: el lecho y las orillas están dominados por material procedente de las laderas cercanas
por acción de la gravedad.
• Aluvial: todo el lecho y las orillas están conformados por sedimentos del río.
• Mixto: Combinación de alguno de los anteriores.
A lo largo de la fuente se observa que su depósito está constituido por sustratos mixtos, principalmente
aluviales de lecho, otros pocos se tratan de arenas mal gradadas y con compuestosde arcillas con
limos.
De acuerdo con lo anterior, las categorías de los sedimentos pueden ser las siguientes:
Página 48 de 156
Gravas de 2 mm a 64 mm
LODO Arenas de 0,063 mm a 2 mm
Limos y Arcillas menos de 0,063 mm
FINO Lodos antropogénicos Lodos procedentes de problemas de calidad de las
aguas
Fuente: Ministerio para la transición Ecológica (2019)
Por la clasificación anterior, se puede concluir que el lecho de la corriente está compuesto por arenas.
Por otro lado, los laterales se componen por gravas y arenas.
• Formación de barras: observando si están bien formadas (completas, con formas nítidas,
generalmente apuntadas hacia aguas abajo, con huellas de que las crecidas las han elevado
o reconstruido) o bien son incipientes o mal estructuradas. Para este indicador es importante
poder reconocer el río en aguas bajas y comprobar también la morfología en fotografía aérea.
Página 49 de 156
casos se prestará especial atención a la formación de depósitos, es decir, a que el material esté bien
distribuido, con aspecto de naturalidad en el lecho y en los depósitos laterales.
Como conclusión final, se procederá a establecer las siguientes categorías para la clasificación de los
sedimentos del subtramo:
Página 50 de 156
c) Acorazamiento excesivo d) Acorazamiento débil
Página 51 de 156
f) Barras incipientes, mal estructuradas
Fuente: (Ministerio para la transición Ecológica, 2019)
Con la correspondiente definición anterior, se puede indicar claramente que el lecho posee una
clasificación efectiva, donde el lecho se encuentra imbricado en arena, barras laterales bien
conformadas con un acorazamiento debil.
Es la secuencia de elementos del fondo del cauce, en el sentido de la pendiente dominante a visualizar
en el subtramo de muestreo. En grandes ríos es factible ver en fotografía aérea esta estructura, pero
en la mayor parte de los cursos fluviales será necesario registrarla directamente en campo, siendo
muy fácil de identificar. Se deberán seleccionar las tipologías principales de las que aparecen en el
listado adjunto, especificándose sus características. En caso de que la estructura longitudinal del
subtramo de muestreo sea fruto de las alteraciones humanas deberá señalarse este aspecto.
Página 52 de 156
Figura 23. Tipo de estructuras longitudinales
Poza/marmita de gigante (Pmg), Salto/poza (Spz), Rápido/poza (Rpz), Rápido/remanso (Rre), Rápido contínuo (Rc),
Grada (Gr), Rampa (Ra), Tabla (Ta), Indeterminada (In)
Fuente: (Ministerio para la transición Ecológica, 2019)
Para el análisis del componente se recurrió a la utilización del eje perfil levantado en la batimetría del
proyecto, donde el puente actual se cruza a 87 grados con la corriente actual en la abscisa propuesta.
Página 53 de 156
Figura 24. Tipo de estructuras longitudinales presentes en el proyecto
Poza/marmita de gigante (Pmg), Salto/poza (Spz), Rápido/poza (Rpz), Rápido/remanso (Rre), Rápido contínuo (Rc),
Grada (Gr), Rampa (Ra), Tabla (Ta), Indeterminada (In)
Fuente: Elaboración propia, 2022
Con la imagen anterior se puede deducir que la proyección de la nueva estructura se localiza en un
Rápido-Poza sucedida por un Remanso. Esta situación es muy importante en la definición del diseño,
porque garantiza la movilidad de los fluidos en razón de masas, se aquieta, y finalmente se distribuye
en el espejo final. Sin embargo, el proceso de imbricación de salida sugiere intervenciones para evitar
que el cauce de constriña y disminuya el área hidráulica.
• Si el sedimento del propio lecho mojado o de las barras de las orillas está suelto y es fácilmente
removible manualmente, o bien está compactado o encostrado.
• Si el material sedimentario aparece cubierto de una pátina o capa de finos, que puede ser
síntoma de problemas de inactividad en los procesos hidromorfológicos.
• La presencia de madera muerta transportada (troncos, ramas y otros restos) y/o arribazones,
son buenos indicadores de dinámica geomorfológica y efectividad de transporte de las
crecidas, además de intervenir en los procesos de sedimentación generando efectos trampa,
sombra o pantalla.
Página 54 de 156
La movilidad de los sedimentos está íntimamente relacionada con el funcionamiento de la cuenca
hidrográfica y en especial, con el régimen de caudales. Cuencas con importante regulación, donde la
generación de sedimentos estará alterada, así como su régimen de caudales, especialmente en
avenidas, tendrán pocos sedimentos y estarán muy fijados por la vegetación riparia, tanto macrófitos
como vegetación leñosa.
Como conclusión de este apartado, la movilidad de sedimentos se clasificará en estas tres categorías:
• Movilidad efectiva (Me): el sedimento está suelto y es fácilmente movible, no está cubierto por
una capa de finos, no hay colonización vegetal o ésta es muy débil, hay madera muerta
transportada y/o arribazones integrados con los sedimentos.
• Movilidad mediana (Mm): el sedimento está suelto, pero por tamaño, proceso de imbricación
o acorazamiento no es fácilmente transportable por la masa de agua corriente.
• Movilidad nula (Mn): el sedimento está encostrado superficialmente o bien cubierto de una
capa continua de material fino o está totalmente colonizado por vegetación.
Como se indicó en el componente de estructuras verticales, en la zona del puente existe un proceso
de socavación de fondo con conformación de pozas, pero que viene dinamizado por un rápido
antecesor, y una descarga de otro río que intenta limitarlo, peor por las condiciones de calado, este
se supera y no reviste mayor dificultad en la entrega. La movilidad mediana se da por el proceso de
agregación de suelos de salida en la poza por la rampa ascendente.
Página 55 de 156
Figura 25. Tipo de movilidad de los sedimentos
Movilidad efectiva (Me), Movilidad mediana (Mm), Movilidad limitada (Ml), Movilidad nula (Mn)
Fuente: Elaboración propia, 2021
3. HIDROGRAFÍA
Debido a la localización geográfica y las condiciones el paisaje en el cual se localiza los municipios de
Orito y Puerto Caicedo, el biotipo de lomas y llanuras predominantes. Este paisaje está acompañado
por el tránsito inicial de grandes ríos, que por condición topográfica los condiciona a buscar
pendientes de menor esfuerzo cuando estos llegan a las geoformas fluviales en paisaje de llanura.
La fuente en particular, es una fuente que recibe muchas corrientes aportantes desde zona de
montaña y piedemonte. Su nacimiento se localiza sobre los 3551 msnm y desemboca en el río
Putumayo sobre los 240 msnm.
Página 56 de 156
3.1. HIPSOMETRÍA DE LA CUENCA
Las alturas con mayor frecuencia de área son las comprendidas entre los 1800 y 1500 msnm, es decir,
la zona de mayor ocupación en la cuenca, es el área de transporte, donde se localizan los afleuntes
tributarios permanentes de mayor proporción, y por consecuencia de la distribución, el puente se
localiza en la zona de depósito.
Estos parámetros están relacionados con el régimen hidrológico de una cuenca o microcuenca ya
que es una función de numerosos factores, entre los que predomina el clima y la forma del terreno en
el cual se desarrolla el fenómeno. Las formas de la superficie terrestre, y en particular su situación en
altitud, tienen influencia decisiva sobre los más importantes factores condicionantes del régimen
hidrológico, como precipitación, escorrentía, infiltración y formación de depósitos y sedimentos. Los
valores morfométricos son fundamentales para documentar la analogía territorial y establecer
relaciones hidrológicas de generalización y expresan en valores simples las características de paisajes
complejos.
Página 57 de 156
Estos elementos físicos proporcionan la más conveniente posibilidad de conocer la variación en el
espacio de los elementos del régimen hidrológico; las características físicas o morfométricas
calculadas para la corriente.
Para los efectos de plantear ajustes de ingeniería en diseño de obras de paso, la Hidrología orientada
a la estimación de los caudales de Escorrentía Superficial que se genera sobre una Cuenca
Hidrográfica ante la ocurrencia de determinada precipitación de diseño. Si bien la magnitud del caudal
depende en gran medida de la magnitud de la Precipitación, es muy cierto que la forma de la cuenca
juega un papel predominante en la respuesta de ésta a dicha precipitación2.
2
Publicado el 19/06/2015 en http://ingenieriacivil.tutorialesaldia.com/algunos-parametros-de-forma-y-drenaje-de-la-cuenca-
hidrografica-y-su-relacion-con-el-tiempo-de-concentracion/ Recuperado el 17/04/2016
Página 58 de 156
Figura 26. Microcuencas de la cuenca del río San Juan
De acuerdo con su extensión, existen algunas sub-áreas que se clasifican como se presenta en la
tabla enseguida:
Tabla 10. Definición de cuencas por tamaño
AC (KM2) NOMBRE
<5 Unidad
5-20 Sector
Página 59 de 156
AC (KM2) NOMBRE
20-100 Microcuenca
100-300 Subcuenca
>300 Cuenca
Fuente: (Trujillo, Barroso, & Escobar, 2010)
Para el caso de la cuenca San Juan, el área de drenaje es de 1205 km2, es decir, se denomina como
Cuenca en nivel de la clasificación de Trujillo et al., (2010).
a) La divisoria de aguas pasa por los puntos más altos de las cordilleras, cruzando los valles que
estas delimitan.
b) Su delimitación comienza en el punto de concentración y se continúa a cada lado de este
punto con líneas siempre perpendiculares a las curvas de nivel.
c) La divisoria de aguas nunca debe interceptar los cauces naturales.
No obstante, y teniendo en cuenta que el área de análisis no requiere determinar los parámetros de
las microcuencas aportantes hasta el sitio del proyecto, porque la conformación hidrográfica ofrece
líneas mínimas de tributación.
El ancho promedio es la relación existente entre el área superficial y la longitud axial de la corriente
principal. Esta característica es muy importante para determinar el factor forma de la microcuenca y
la susceptibilidad a la torrencialidad (Escuela superior de administración Pública - ESAP, 2012).
𝐴 Ec. 3
𝐴𝑝 =
𝐿𝑏
Página 60 de 156
Tabla 11. Relación de ancho promedio con factor de forma y características de torrencialidad
ANCHO PROMEDIO FORMA CARACTERÍSTICAS
≤ 0.26 Oval, oblonga, rectangular Menor susceptibilidad a la
torrencialidad
0.26 - 0.52 Oval redonda a oval oblonga Susceptibilidad media a la
torrencialidad
˃ 0.52 Oval redonda a redonda Mayor susceptibilidad a la
torrencialidad, por consiguiente, se
debe priorizar como una microcuenca
de manejo especial. Uso de cobertura
Fuente: Horton (1945)
Para el caso de la microcuenca el ancho promedio es de 20.4, es decir, una forma Oval redonda.
Se define como la relación entre el diámetro de un círculo que posea la misma área de la cuenca y
cuyo diámetro sea igual la longitud de la cuenca.
El valor de la relación de elongación se acerca a la unidad cuando la cuenca es muy plana y circular,
cuando la cuenca es plana con porciones accidentales, la relación de elongación está entre 0.5 y 0.8.
√𝐴 Ec. 4
𝑅𝑒 = 1.128
𝐿
Donde A corresponde al área de la cuenca expresada en Km2, y L a la Longitud del cauce principal
calculada en Km.
La relación de elongación de la cuenca de análisis es 0.64, es decir, plana con algunas proporciones
accidentadas.
Es la relación entre el ancho medio de la cuenca y la longitud de su cauce principal. El ancho medio
se obtiene cuando se divide el área de la cuenca por la longitud del cauce principal, por lo tanto, el
Coeficiente de Forma queda definido así3:
𝐴 Ec. 5
𝐾𝑓 =
𝐿2
3
Ibid. 5
Página 61 de 156
Según Salvatierra-Colindres (2006) citado por Karalis et al., (2014), la relación de forma da alguna
indicación de la tendencia de las avenidas en el cauce, porque una cuenca con un factor de forma
bajo, tiene menos tendencias a concentrar las intensidad de lluvias que una cuenca de igual área.
Para la cuenca se encontró un factor de forma de 0.32, es decir, una microcuenca con bajas
capacidades de concentrar lluvias.
El índice propuesto por Graveluis (1914) citado por UNAD (2012), compara la forma de la cuenca con
la de una circunferencia, cuyo círculo inscrito tiene la misma área de la cuenca en estudio. Se define
como la razón entre el perímetro de la cuenca que es la misma longitud del parteaguas o divisoria que
la encierra y el perímetro de la circunferencia. El grado de aproximación de este índice a la unidad
indicará la tendencia a concentrar fuertes volúmenes de aguas de escurrimiento, siendo más
acentuado cuanto más cercano a uno sea, es decir mayor concentración de agua 4.
𝑃 Ec. 6
𝐾𝑐 = 0.28 ( )
√𝐴
Para la cuenca, el coeficiente de compacidad es 1.15, es decir, no concentraría fuertes volúmenes de agua.
Los parámetros hidrográficos corresponden a cuyos factores que por condiciones de red hídrica co-
localizada en un espacio geográfico determina comportamientos de flujo en la cuenca hidrográfica.
𝑁𝑢 Ec. 7
𝑅𝑏 =
𝑁𝑢+1
4
Publicado el 12/10/2012 en http://es.slideshare.net/nuevvo/4-geomorfologia Recuperado el 17/04/2016
5
Ibid., 6
Página 62 de 156
3.2.2.2. Pendiente media del cauce principal – Pmc
Es la diferencia total de elevación del cauce principal (cota máxima – cota mínima), dividida por su
longitud total.
𝐻 Ec. 8
𝑃𝑚𝑐 =
𝐿
Dependiendo de la forma que presente el o los Talweg que posee la cuenca, se puede hablar de
patrones de drenaje, entendido esto como la forma que tiene el o los cauces encargados de drenar o
extraer el agua de los terrenos circundantes (vertientes), en (Rojo, 1997).
Comparando la geometría del drenaje de la cuenca con lo descrito en la Figura 27, la cuenca de
análisis tiene un patrón de drenaje Dendrítico.
6
Tomado de http://4.bp.blogspot.com/_w1kycNNBkOE/TByjswg6VdI/AAAAAAAADeg/7o76RtwkyKA/s1600/t2.jpg Recuperado el
05/06/2016
Página 63 de 156
3.2.2.4. Coeficiente de torrencialidad
Índice que mide el grado de torrencialidad de la cuenca, por medio de la relación del número de
cauces de orden uno con respecto al área total de la misma. A mayor magnitud, mayor grado de
torrencialidad presenta una cuenca.
Con el fin de catalogar una cuenca como bien o mal drenada, analizando su densidad de drenaje, se
puede considerar que valores de densidad de drenaje próximo a 0.5 km/km 2 o mayores indican la
eficiencia de la red de drenaje7.
Valores bajos de Dd generalmente están asociados con regiones de alta resistencia a la erosión, muy
permeables y de bajo relieve. Valores altos fundamentalmente son encontrados en regiones de suelos
impermeables, con poca vegetación y de relieve montañoso8.
𝐿𝑡 Ec. 10
𝐷𝑑 =
𝐴
7
Ibid 6
8
Ibíd. 8
Página 64 de 156
3.2.3. Tiempo de concentración
Se define como el tiempo mínimo necesario para que todos los puntos de una cuenca estén aportando
agua de escorrentía de forma simultánea al punto de salida, punto de desagüe o punto de cierre. Está
determinado por el tiempo que tarda en llegar a la salida de la cuenca el agua que procede del punto
hidrológicamente más alejado, y representa el momento a partir del cual el caudal de escorrentía es
constante, al tiempo que máximo; el punto hidrológicamente más alejado es aquél desde el que el
agua de escorrentía emplea más tiempo en llegar a la salida (Ibáñez, Moreno, & Gisbert, 2010).
El cálculo del tiempo de concentración de las microcuencas de análisis fue estimado a través de las
diferentes metodologías aplicables y validadas en diferentes escenarios similares, así mismo se
determinó el valor promedio, desviación standard, la mediana ponderada y mayor de datos aportantes
para escogencia del valor más adecuado. Los datos asociados a todas las subcuencas están en el
anexo (\\10_Puente_Vides\02_HIDROLOGIA\01_MORFOMETRIA)
Como el elemento que permite determinar el coeficiente de escorrentía es la cobertura del suelo,
haciendo uso de la metodología de CORINE Land Cover adaptada para Colombia (IDEAM, 2010a) en
nivel de interpretación III, se determinó la composición de cubierta de la microcuenca. La cobertura
se obtuvo del Sistema de Información Ambiental de Colombia9 desarrolladas por el SINCHI para el
año 2020.
Así las cosas, puede notarse que se trata de una cuenca cuya predominancia está dada por la
vegetación secundaria en transición, bosques densos en tierra firme, pastos enmalezados, cultivos
con espacios naturales y bosques abiertos. El esquema espacial de la cuenca en términos de
cobertura se espacializa en la siguiente Figura.
9
http://www.siac.gov.co/catalogo-de-mapas
Página 65 de 156
Figura 28. Coberturas del suelo de la microcuenca
4. HIDROLOGÍA
El estudio hidrológico tiene como objeto determinar la intensidad máxima de precipitación o lluvia de
diseño, con la cual se estima las condiciones de localización, espaciamiento y fijación de las unidades
de descole. Para estimar dicha condición de diseño se hizo uso de los datos IDEAM para los registros
máximo históricos multianuales de cada una de las estaciones pluviométricas aferentes al trazado de
diseño. Dicha aferencia fue determinada a través de la construcción de polígonos de Thiessen
(diagramas de Voronoi - teselación de Dirichlet) (Bateman, 2007).
Página 66 de 156
Debido a diferentes razones, estas cuencas han sido poco estudiadas, la información hidrológica
existente se remonta a datos limnimétricos con poca abundancia de registros que no simplifican mayor
representatividad, y en cuanto a la regulación a nivel de Planes de Ordenamiento de Cuencas, Plan
de Ordenamiento de Ronda Hídrica y Planes de Ordenación del Recurso Hídrico, son inexistentes,
con referencia a estudios de hidráulica, éstos son muy mínimos, solo algunas aproximaciones para el
desarrollo de infraestructura, principalmente de mitigación de amenazas.
En ese sentido, cuando la ausencia de datos hidroclimatológicos son ausentes, se pretende satisfacer
esa información con una serie de operaciones y cálculos para estimar los valores de escorrentía de
una cuenca, bajo ese contexto, el cálculo de parámetros físicos y morfométricos guardan con la
magnitud del caudal una relación directa, en particular con la respuesta de ésta a máximas
precipitaciones10. Es por eso que se deben adelantar estudios hidrológicos eficientes, iniciando con
los cálculos morfométricos de las cuencas, dado que los aspectos físicos que determinan la forma de
una cuenca están relacionados con la fisiografía, geomorfología del terreno y particularmente el
comportamiento hidráulico.
Es así como se obtiene los valores hidrológicos del área de estudio, basado en la estimación de la
intensidad, tiempo y duración de la precipitación, retorno, el coeficiente de escorrentía, impluvium y
parámetros hidrológicos regionales para aplicación de la Metodología Simplificada de cálculo de las
curvas intensidad – duración – frecuencia indicada en el Manual de Drenaje de Carreteras del INV
2009.
La información que permitió obtener el presente cálculo proviene de los registros de valores máximos
de precipitación en 24 horas para las diferentes estaciones.
Identificar, para cada año de la serie histórica, el valor máximo de precipitación registrado en
veinticuatro (24) horas.
La intensidad de la lluvia calcularla para diferentes duraciones de aguacero y para cada año de la
serie histórica. Generalmente se utilizan duraciones de aguacero de 5, 10, 15, 20, 25 y 30 minutos.
10
Publicado el 19/06/2015 en http://ingenieriacivil.tutorialesaldia.com/algunos-parametros-de-forma-y-drenaje-de-la-cuenca-
hidrografica-y-su-relacion-con-el-tiempo-de-concentracion/ Recuperado el 17/04/2016
Página 67 de 156
Aplicar los valores de los periodos de retorno (T) recomendados en el Manual de Drenaje Vial del
INV 2009, que se indican en la Tabla 12 a continuación. Estimación de la precipitación máxima
probable
De acuerdo con lo anterior, se usará el Tiempo de Retorno de 100 años, teniendo como argumentos
la longitud de la onda del meando en el cauce (103.2 m).
Página 68 de 156
TIPO DE ÁREA DE DRENAJE COEFICIENTE DE
ESCORRENTÍA C
SUBURBANA 0.25 ˍ 0.40
ÁREAS DE APARTAMENTOS DE VIVIENDA 0.50 ˍ 0.70
CALLES
ASFALTADAS 0.70 ˍ 0.95
CONCRETO 0.80 ˍ 0.95
LADRILLO 0.70 ˍ 0.85
CALZADAS Y ALAMEDAS 0.75 ˍ 0.85
TECHOS 0.75 ˍ 0.95
Adaptado de: MONSALVE, 1999, citado por INVIAS & Ministerio de Transporte, 2009
La ecuación de intensidad es un modelo matemático iterativo donde se obtienen los las relaciones o
cocientes a la lluvia de 24 horas que se emplean para duraciones de varias horas, el presente cálculo
se estima a partir del concepto de hidrosistema mixto, donde se procesa los factores de confianza y
se contrasta con los resultados pluviométricos de la estación meteorológica. El modelo aplica a las
zonas geográficas disímiles como la que trata la guía del (Armenta Porras, 2013).
Página 69 de 156
INVIAS 2009, en el cual se recomienda la utilización de la precipitación máxima promedio anual en 24
horas y aplicar la expresión que a continuación se indica:
𝑎 × 𝑇 𝑏 × 𝑀𝑑 Ec. 11
𝑖=
𝑡𝑐
60
Página 70 de 156
MUNICIPIO CORRIENTE ESTACIÓN CÓDIGO TIPO
METEOROLÓGICA ESTACIÓN
SAN Putumayo Torre TV 47010180 Pluviométrica
FRANCISCO
MOCOA Putumayo El Pepino 44015070 Climática Principal
VILLAGARZÓN Putumayo Villagarzón 44015010 Agroclimatológica
De acuerdo con lo anterior y la imagen asociada a la referencia, las estaciones aferentes a la cuenca
son siete, se debe tener en cuenta que todas las estaciones meteorológicas están influyendo en la
cuenca, pero para los cálculos individuales se subcuencas se realizará una intersección entre
polígonos de Thiessen con las hoyas hidrográficas:
Página 71 de 156
SUBCUENCA ESTACIÓN AFERENTE % DE AFERENCIA
C2 EL PEPINO - AUT [44015070] 21.7%
C2 PUERTO CAICEDO [47010110] 1.4%
C2 TORRE TV SAN FCO [47010180] 10.4%
C2 VILLAGARZON [44015010] 66.6%
C20 PUERTO CAICEDO [47010110] 100.0%
C21 PUERTO CAICEDO [47010110] 92.1%
C21 VILLAGARZON [44015010] 7.9%
C22 VILLAGARZON [44015010] 100.0%
C23 PUERTO CAICEDO [47010110] 100.0%
C3 SANTA ISABEL [47010230] 100.0%
C4 CHUNGACASPI [47010050] 5.6%
C4 EL PEPINO - AUT [44015070] 13.2%
C4 TORRE TV SAN FCO [47010180] 48.5%
C4 VILLAGARZON [44015010] 32.7%
C5 BALSAYACO [47010020] 60.2%
C5 SANTA ISABEL [47010230] 39.8%
C6 SANTA ISABEL [47010230] 100.0%
C7 PUERTO CAICEDO [47010110] 7.1%
C7 VILLAGARZON [44015010] 92.9%
C8 BALSAYACO [47010020] 30.9%
C8 PUERTO CAICEDO [47010110] 54.6%
C8 SANTA ISABEL [47010230] 10.6%
C8 VILLAGARZON [44015010] 3.9%
C9 PUERTO CAICEDO [47010110] 18.6%
C9 VILLAGARZON [44015010] 81.4%
Fuente: Elaboración propia, 2022
Página 72 de 156
Figura 29. Polígonos de Thiessen asociados al proyecto
4.3.1.1. Precipitación
Considerando las series de tiempo originales, se procede a realizar las curvas de masas simple para
cada estación seleccionada.
Se observa que para la mayoría de las series de precipitación la curva de masas presenta una
pendiente constante (sin saltos ni cambios bruscos). Sin embargo, las estaciones de piedemonte y
zona de llanura albergan mayor concentración de intensidades y mayor precipitación acumulada.
Página 73 de 156
Precipitación total mensual (mm) Precipitación total mensual (mm)
100
120
0
20
40
60
80
100
150
200
250
0
50
1/11/1978 1/07/1989
1/01/1980 1/06/1990
1/03/1981 1/05/1991
1/05/1982 1/04/1992
1/07/1983 1/03/1993
1/09/1984 1/02/1994
1/11/1985 1/01/1995
1/01/1987
1/12/1995
1/03/1988
1/11/1996
1/05/1989
1/10/1997
1/07/1990
1/09/1998
1/09/1991
1/08/1999
1/11/1992
1/01/1994 1/07/2000
1/03/1995 1/06/2001
1/05/1996 1/05/2002
1/07/1997 1/04/2003
1/09/1998 1/03/2004
1/11/1999 1/02/2005
1/01/2001 1/01/2006
1/03/2002 1/12/2006
Página 74 de 156
1/05/2003 1/11/2007
1/07/2004 1/10/2008
1/09/2005 1/09/2009
1/11/2006
1/08/2010
1/01/2008
1/07/2011
1/03/2009
1/06/2012
0
0
20000
40000
60000
80000
-20000
100000
120000
140000
160000
180000
10000
20000
30000
40000
50000
60000
Figura 31. Serie temporal de precipitación y curva de masa simple de la estación Puerto Caicedo
Precipitación total mensual (mm) Precipitación total mensual (mm)
100
120
140
160
180
0
20
40
60
80
100
120
140
160
0
20
40
60
80
1/06/1971 1/02/1979
1/11/1972 1/04/1980
1/04/1974 1/06/1981
1/09/1975 1/08/1982
1/02/1977 1/10/1983
1/07/1978 1/12/1984
1/12/1979 1/02/1986
1/05/1981 1/04/1987
1/10/1982 1/06/1988
1/03/1984 1/08/1989
1/08/1985 1/10/1990
1/01/1987 1/12/1991
1/02/1993
1/06/1988
1/04/1994
1/11/1989
1/06/1995
1/04/1991
1/08/1996
1/09/1992
1/10/1997
1/02/1994
1/12/1998
1/07/1995
1/02/2000
1/12/1996
1/04/2001
1/05/1998
1/06/2002
Página 75 de 156
1/10/1999
1/08/2003
1/03/2001
1/10/2004
1/08/2002
1/12/2005
1/01/2004 1/02/2007
1/06/2005 1/04/2008
1/11/2006 1/06/2009
0
0
20000
40000
60000
80000
50000
100000
120000
140000
160000
100000
150000
200000
250000
Figura 32. Serie temporal de precipitación y curva de masa simple de la estación Torre TV
100
100
150
200
250
300
0
0
10
30
40
50
60
70
80
90
50
20
1/01/1970 1/09/1965
1/06/1971 1/03/1967
1/11/1972 1/09/1968
1/04/1974 1/03/1970
1/09/1975 1/09/1971
1/02/1977 1/03/1973
1/07/1978 1/09/1974
1/12/1979 1/03/1976
1/05/1981 1/09/1977
1/10/1982 1/03/1979
1/03/1984 1/09/1980
1/08/1985 1/03/1982
1/01/1987 1/09/1983
1/06/1988 1/03/1985
1/11/1989 1/09/1986
1/04/1991 1/03/1988
1/09/1992 1/09/1989
1/03/1991
1/02/1994
1/09/1992
1/07/1995
1/03/1994
1/12/1996
1/09/1995
1/05/1998
Página 76 de 156
1/03/1997
1/10/1999
1/09/1998
1/03/2001
1/03/2000
1/08/2002
1/09/2001
1/01/2004
1/03/2003
1/06/2005
1/09/2004
1/11/2006
1/03/2006
0
0
20000
40000
60000
80000
50000
100000
120000
100000
150000
200000
250000
300000
Figura 35. Serie temporal de precipitación y curva de masa simple de la estación Balsayaco
140 35000
100 25000
80 20000
60 15000
40 10000
20 5000
0 0
1/01/2002
1/05/2002
1/09/2002
1/01/2003
1/05/2003
1/09/2003
1/01/2004
1/05/2004
1/09/2004
1/01/2005
1/05/2005
1/09/2005
1/01/2006
1/05/2006
1/09/2006
1/01/2007
1/05/2007
1/09/2007
1/01/2008
1/05/2008
1/09/2008
1/01/2009
1/05/2009
1/09/2009
1/01/2010
1/05/2010
1/09/2010
1/01/2011
1/05/2011
1/09/2011
1/01/2012
1/05/2012
1/09/2012
1/01/2013
1/05/2013
Fuente: Elaboración propia, 2022
Las lluvias de la zona alta de la orogenia son en mayor parte de distribución unimodal, mientras que
la parte baja se comportan de manera tendiente a generar dos momentos pluviométricos. Las
precipitaciones de mayor magnitud se presentan dentro del segundo cuatrimestre del año en todas
las secciones, es decir, entre abril y agosto, mientras que los periodos de menor intensidad mensual
se al final e inicio de año. Las lluvias máximas promedio anual se manifiestan superiores a los 3000
mm/año, en el piso altitudinal comprendido entre los 1500 a 200 msnm.
Página 77 de 156
Figura 37. Distribución precipitación mensual media - Santa Isabel
140.0
120.0
100.0
Precipitación (mm)
80.0
60.0
40.0
20.0
0.0
350.0
300.0
Precipitación (mm)
250.0
200.0
150.0
100.0
50.0
0.0
Página 78 de 156
Figura 39. Distribución precipitación mensual media – Torre TV
700.0
600.0
500.0
Precipitación (mm)
400.0
300.0
200.0
100.0
0.0
300.0
250.0
200.0
150.0
100.0
50.0
0.0
Página 79 de 156
Figura 41. Distribución precipitación mensual media – Villagarzón
600.0
500.0
400.0
Precipitación (mm)
300.0
200.0
100.0
0.0
200.0
Precipitación (mm)
150.0
100.0
50.0
0.0
Página 80 de 156
Figura 43. Distribución precipitación mensual media – El Pepino
450.0
400.0
350.0
300.0
Precipitación (mm)
250.0
200.0
150.0
100.0
50.0
0.0
Página 81 de 156
Fuente: Elaboración propia, 2022
Para la zona en estudio, se realizó el análisis de lluvias máximas para diferentes períodos de retorno
a las estaciones seleccionadas. Para el análisis de las lluvias máximas en 24 horas, se recopilaron los
máximos valores anuales de precipitación en 24 horas registrada, mediante las distribuciones de
probabilidad Normal, Gumbel, Lognormal III, Pearson, Log Pearson III y EV3 se hicieron ajustes y
empleando la prueba de Chi-cuadrado se determina cual distribución de probabilidad presenta mejor
ajuste. Las Tablas a continuación presentan los valores obtenidos por los diferentes métodos para
períodos de retorno de 5,10,15,25, 50 y 100 años.
Página 82 de 156
Tabla 25. Análisis de Frecuencia de Lluvia Máxima en 24 horas - Estación Santa Isabel
TR NORMAL GUMBEL PEARSON LOG- PEAR LOG-NOR EV3
AÑOS (mm) (mm) (mm) (mm) (mm) (mm)
5 30.8 30.9 30.8 30.9 30.5 NA
10 33.4 34.8 33.4 33.6 33.6 NA
15 34.7 37.0 34.7 35.0 35.3 NA
25 36.1 39.7 36.1 36.5 37.3 NA
50 37.9 43.3 37.9 38.3 39.9 NA
100 39.5 47.0 39.5 39.9 42.4 NA
CHI 2 1.201 2.479 1.201 1.173 1.994 NA
Fuente: Elaboración propia, 2022
Figura 44. Frecuencia de precipitación máxima en 24 horas a nivel multianual estación Santa Isabel
300
Serie
Normal
Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)
Normal
Gumbel
200 Pearson
Log
Pearson
EV3
Abscisa
100
Página 83 de 156
Tabla 26. Análisis de Frecuencia de Lluvia Máxima en 24 horas - Estación Puerto Caicedo
TR NORMAL GUMBEL PEARSON LOG- PEAR LOG-NOR EV3
AÑOS (mm) (mm) (mm) (mm) (mm) (mm)
5 79.2 79.0 79.3 80.1 78.1 78.9
10 86.5 89.8 83.5 82.1 87.2 83.6
15 90.2 95.9 85.4 82.5 92.1 85.7
25 94.3 103.5 87.2 82.7 98.0 87.9
50 99.3 113.6 89.0 82.8 105.8 90.5
100 103.9 123.7 90.4 82.8 113.2 92.6
CHI 2 21.971 46.108 6.952 5.666 38.588 7.698
Fuente: Elaboración propia, 2022
Figura 45. Frecuencia de precipitación máxima en 24 horas a nivel multianual estación Puerto Caicedo
300
Serie
Normal
Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)
Normal
Gumbel
200 Pearson
Log
Pearson
EV3
Abscisa
100
Página 84 de 156
Tabla 27. Análisis de Frecuencia de Lluvia Máxima en 24 horas - Estación Torre TV
TR NORMAL GUMBEL PEARSON LOG- PEAR LOG-NOR EV3
AÑOS (mm) (mm) (mm) (mm) (mm) (mm)
5 58.1 58.0 58.0 58.1 57.9 NA
10 61.1 62.5 61.3 61.2 61.3 NA
15 62.6 65.0 63.1 62.7 63.1 NA
25 64.4 68.2 65.1 64.5 65.2 NA
50 66.4 72.4 67.6 66.6 67.9 NA
100 68.3 76.5 69.9 68.5 70.3 NA
CHI 2 1.538 1.832 1.467 1.555 1.461 NA
Fuente: Elaboración propia, 2022
Figura 46. Frecuencia de precipitación máxima en 24 horas a nivel multianual estación Torre TV
300
Serie
Normal
Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)
Normal
Gumbel
200 Pearson
Log
Pearson
EV3
Abscisa
100
Página 85 de 156
Tabla 28. Análisis de Frecuencia de Lluvia Máxima en 24 horas - Estación Chungacaspi
TR NORMAL GUMBEL PEARSON LOG- PEAR LOG-NOR EV3
AÑOS (mm) (mm) (mm) (mm) (mm) (mm)
5 44.1 43.8 43.8 44.6 43.1 43.9
10 48.8 50.7 45.9 44.7 49.2 46.8
15 51.1 54.6 46.7 45.1 52.5 48.2
25 53.8 59.3 47.4 46.4 56.5 49.6
50 57.0 65.8 48.0 49.5 61.8 51.2
100 59.9 72.2 48.4 54.9 67.0 52.6
CHI 2 26.989 43.591 13.387 29.371 39.713 15.00
Figura 47. Frecuencia de precipitación máxima en 24 horas a nivel multianual estación Chungacaspi
300
Serie
Normal
Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)
Normal
Gumbel
200 Pearson
Log
Pearson
EV3
Abscisa
100
Página 86 de 156
Tabla 29. Análisis de Frecuencia de Lluvia Máxima en 24 horas - Estación Villagarzón
TR NORMAL GUMBEL PEARSON LOG- PEAR LOG-NOR EV3
AÑOS (mm) (mm) (mm) (mm) (mm) (mm)
5 86.8 86.3 86.9 87.4 86.1 86.5
10 93.2 95.6 90.8 90.2 93.7 90.6
15 96.4 100.9 92.5 91.1 97.8 92.4
25 100.0 107.4 94.2 91.8 102.6 94.4
50 104.4 116.1 96.0 92.3 108.8 96.6
100 108.4 124.7 97.4 92.5 114.7 98.4
CHI 2 14.185 32.873 3.981 4.083 23.791 3.753
Fuente: Elaboración propia, 2022
Figura 48. Frecuencia de precipitación máxima en 24 horas a nivel multianual estación Villagarzón
300
Serie
Normal
Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)
Normal
Gumbel
200 Pearson
Log
Pearson
EV3
Abscisa
100
Página 87 de 156
Tabla 30. Análisis de Frecuencia de Lluvia Máxima en 24 horas - Estación Balsayaco
TR NORMAL GUMBEL PEARSON LOG- PEAR LOG-NOR EV3
AÑOS (mm) (mm) (mm) (mm) (mm) (mm)
5 33.8 33.6 33.8 33.4 33.5 33.7
10 36.4 37.4 35.2 33.6 36.6 35.3
15 37.6 39.5 35.8 33.6 38.2 36.0
25 39.1 42.1 36.4 33.6 40.2 36.8
50 40.8 45.6 37.0 33.7 42.6 37.7
100 42.4 49.1 37.4 34.0 45.0 38.4
CHI 2 6.347 11.439 3.686 6.472 8.925 3.518
Fuente: Elaboración propia, 2022
Figura 49. Frecuencia de precipitación máxima en 24 horas a nivel multianual estación Balsayaco
300
Serie
Normal
Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)
Normal
Gumbel
200 Pearson
Log
Pearson
EV3
Abscisa
100
Página 88 de 156
Tabla 31. Análisis de Frecuencia de Lluvia Máxima en 24 horas - Estación Pepino
TR NORMAL GUMBEL PEARSON LOG- PEAR LOG-NOR EV3
AÑOS (mm) (mm) (mm) (mm) (mm) (mm)
5 77.9 80.5 78.1 78.8 77.0 77.6
10 84.6 92.2 83.7 84.4 85.2 81.8
15 88.0 98.7 86.4 87.0 89.6 83.8
25 91.8 106.9 89.2 89.6 95.0 85.8
50 96.4 117.8 92.6 92.4 101.9 88.2
100 100.5 128.6 95.4 94.5 108.5 90.1
CHI 2 2.937 4.437 2.499 2.347 4.539 2.938
Fuente: Elaboración propia, 2022
Figura 50. Frecuencia de precipitación máxima en 24 horas a nivel multianual estación Pepino
300
Serie
Normal
Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)
Normal
Gumbel
200 Pearson
Log
Pearson
EV3
Abscisa
100
Página 89 de 156
4.4. CURVAS DE INTENSIDAD – DURACIÓN Y FRECUENCIA (IDF)
120
Intensidad (mm/hr)
100
80
60
40
20
0
0 20 40 60 80 100
Tiempo de duración (min)
Página 90 de 156
Figura 52. Curvas IDF estación Puerto Caicedo
400
2 años 5 años 10 años 25 años
350
50 años 75 años 100 años
300
Intensidad (mm/hr)
250
200
150
100
50
0
0 20 40 60 80 100
Tiempo de duración (min)
250
Intensidad (mm/hr)
200
150
100
50
0
0 20 40 60 80 100
Tiempo de duración (min)
Página 91 de 156
Figura 54. Curvas IDF estación Chungacaspi
250
2 años 5 años 10 años 25 años
150
100
50
0
0 20 40 60 80 100
Tiempo de duración (min)
300
Intensidad (mm/hr)
250
200
150
100
50
0
0 20 40 60 80 100
Tiempo de duración (min)
Página 92 de 156
Figura 56. Curvas IDF estación Balsayaco
200
2 años 5 años 10 años 25 años
180
50 años 75 años 100 años
160
140
Intensidad (mm/hr)
120
100
80
60
40
20
0
0 20 40 60 80 100
Tiempo de duración (min)
300
Intensidad (mm/hr)
250
200
150
100
50
0
0 20 40 60 80 100
Tiempo de duración (min)
Página 93 de 156
4.5. METODOLOGÍA MODELACIÓN HIDROLÓGICA
El clima regional de Colombia está definido por las propiedades relevantes de los valores medios de
los campos, que están forzados por las condiciones de frontera de la atmósfera, de la distribución de
continente–océano y las variaciones de la radiación solar sobre la superficie, en el espacio y tiempo.
Estos fenómenos atmosféricos son fuertemente influenciados por la distribución de la topografía y la
vegetación en la superficie del territorio colombiano (Guzman, Ruíz, & Cadena, 2014).
De acuerdo a lo anterior, se aplicó un modelo lluvia-escorrentía para definir las crecientes, el cual se
inicia con un estudio detallado de los aguaceros que producen los caudales máximos y se
complementa con una caracterización completa de la cuenca aferente en los seis ríos tributarios y de
los principales canales de aguas lluvias para culminar con la construcción de los hidrogramas
correspondientes.
La estimación del índice de capacidad de infiltración y escurrimiento del suelo se adelantó por medio
de un reconocimiento de campo, con ayuda de información existente, que contengan la
caracterización de la cobertura vegetal de las cuencas de estudio y el tipo de suelos superficiales de
la cuenca hasta esta sección para evaluar el parámetro CN.
La selección del valor CN en cada caso, se realizó considerando los siguientes cuatro parámetros:
Página 94 de 156
• Tipo de suelo hidrológico
• Uso y tratamiento del suelo
• Combinaciones hidrológicas de suelo y vegetación, y
• Condiciones de humedad antecedente a la lluvia de diseño,
La clasificación de los suelos en grupos hidrológicos principales se realizó tomando como base el tipo
hidrológico a partir de planos de suelos agrícolas. Donde no se dispone de esos planos, se hizo uso
de un buen criterio. Usualmente en una misma cuenca se reportan diferentes tipos de suelos. En ese
caso se consideraron por separado, calculando el área ocupada por cada uno de ellos y determinando
el valor CN individualmente para luego ponderarlos y obtener un CN representativo.
Para efectuar una buena clasificación y definir la respuesta hidrológica de los diferentes suelos, en el
documento original del S.C.S. se presenta una lista extensa de suelos agrológicos pertenecientes a
los diferentes grupos hidrológicos. En nuestro medio se deben utilizar criterios más generales, dado
que aún no se han elaborados listas de este tipo. Por lo anterior, lo usual es aplicar criterios
conservadores para compensar la falta de información y evitar errores por defecto.
En la actualidad existe una gran cantidad de software que ayuda a realizar modelación hidrológica
continua y/o de eventos; así mismo es importante tener en cuenta que cada uno de ellos tiene una
aplicación (para una gama de casos determinada), lo cual debe ser tenido en cuenta antes de su
selección.
La aplicación del modelo lluvia-caudal para definir las crecientes se inicia con un estudio detallado de
los aguaceros que producen los caudales máximos y se complementa con una caracterización
completa de las cuencas aferentes, para luego culminar con la construcción de los hidrogramas
correspondientes.
Para determinar el caudal generado por la lluvia se aplicó el modelo lluvia-caudal del Soil Conservation
Service de los E.U.A. Para el efecto, se definen las áreas aferentes, se selecciona el período de retorno
del aguacero de diseño, se calcula el valor de la tormenta con el período de retorno requerido, se
evalúan los tiempos de concentración, se selecciona el Número de Curva CN, el cual es un concepto
similar al coeficiente de escorrentía y, finalmente, se calculan los caudales de diseño que serán el
insumo principal para la modelación hidráulica.
El análisis estadístico de la información de las estaciones con datos de caudales máximos anuales en
el cauce, se realizó para obtener los caudales máximos para diferentes períodos de retorno
Página 95 de 156
seleccionados mediante un análisis estadístico, considerando las diferentes funciones de distribución
de probabilidades (Gumbel, LogNormal, Pearson III, Log-Pearson III,, etc.) aplicadas a las series de
caudales máximos, con el propósito de seleccionar la que mejor represente el comportamiento de los
datos disponibles partiendo de la prueba de bondad con el Chi Cuadrado.
En ausencia de datos de caudal del río y/o canales, se realizó un estudio de identificación y
caracterización de la distribución espacial y temporal de la precipitación máxima en 24 horas en la
zona de estudio, de acuerdo a las estaciones meteorológicas encontradas dentro y/o cercanas a cada
una de las cuencas, dado que el modelo Lluvia - Caudal del Soil Conservation Service requiere, como
insumo, el valor de la precipitación máxima diaria con un período de retorno (2, 5, 10, 15, 25, 50, 75 y
100 años) predeterminado. Para la cuenca se produjo una serie de precipitaciones máximas en 24
horas. La obtención de las lluvias con diferentes períodos de retorno se logró mediante un análisis de
frecuencia utilizando las distribuciones Normal, Log Normal, Pearson Tipo III, Log Pearson tipo III,
Valores Extremos y Gumbel, con el propósito de seleccionar la que mejor represente el
comportamiento de los datos disponibles partiendo de pruebas de bondad con el Chi Cuadrado.
El modelo lluvia-caudal del Soil Conservation Service requiere, como insumo, el valor de la
precipitación máxima diaria con un período de retorno predeterminado. Para el efecto, se
seleccionaron las estaciones de precipitación, como representativa de la cuenca desde su nacimiento
hasta el sector de estudio. Estas son las estaciones que se encuentran dentro y/o cerca de la zona
estudiada.
Para la cuenca se produjo una serie de precipitaciones máximas en 24 horas. La obtención de las
lluvias con diferentes períodos de retorno se logró mediante un análisis de frecuencia utilizando las
distribuciones Normal, Log Normal, Pearson tipo III, Log Pearson tipo III y Gumbel, con el propósito
de seleccionar la que mejor represente el comportamiento de los datos disponibles.
El cálculo del caudal de aguas lluvias para los diferentes períodos de retorno se llevó a cabo utilizando
un método de generación de hidrogramas haciendo uso del modelo lluvia- escorrentía, para lo cual
se obtuvieron los hietogramas de precipitación efectiva, a través del uso de modelos de infiltración,
con el fin de determinar el hidrograma de escorrentía directa.
Debido a las limitaciones que se tienen en las mediciones hidrológicas directas, es necesario
implementar modelos lluvia – escorrentía, los cuales proporcionan una predicción o extrapolación
cuantitativa que puede ser útil en la toma de decisiones, especialmente en cuencas que no han sido
calibradas. Para el cálculo de los hidrogramas se utilizó el modelo del Soil Conservation Service
(S.C.S.) de los Estados Unidos de América.
Página 96 de 156
a la intercepción, detención superficial e infiltración) necesarias a la lluvia total que se precipita sobre
la cuenca para obtener la precipitación de excesos, que luego conformó la escorrentía superficial
directa sobre la cuenca.
Con base en lo anterior, se puede decir que el exceso de precipitación o la escorrentía superficial
directa es siempre menor o igual a la precipitación total; también, que existe una cierta cantidad de
precipitación para la cual no ocurrirá escorrentía y, por tanto, si se le sustrae esta cantidad a la
precipitación total, es posible obtener la escorrentía potencial.
La capacidad de infiltración que tiene el suelo va disminuyendo a medida que pasa el tiempo, debido
a lo cual, al momento de separar la parte de la precipitación de excesos en el hietograma de
precipitación total, se sugiere seguir una curva descendente que refleja la disminución natural de la
capacidad de infiltración.
Página 97 de 156
Figura 59. Disminución de la capacidad de infiltración
Dicho método parte de un número de curva de escorrentía (CN), a partir del cual se calculan las
pérdidas en la cuenca de una tormenta, donde se incluyen la intercepción, la detención superficial y
la infiltración. El número de curva varía en un rango de 1 a 100 y depende de la caracterización de
cuatro propiedades productoras de escorrentía de una cuenca hidrográfica.
Igualmente, el método incluye una clasificación del suelo según su uso y tratamiento. El uso del suelo
pertenece a la cobertura de la cuenca asociado a los cuerpos de agua (tipo de vegetación, tierras en
descanso, usos no agrícolas, cuerpos de agua y superficies impermeables). Por otro lado, el
tratamiento del suelo aplica principalmente a los usos agrícolas y las prácticas mecánicas y de manejo
que se utilicen.
La selección del valor CN en cada caso, se realiza considerando los siguientes cuatro parámetros:
• Tipo de suelo hidrológico.
• Uso y tratamiento del suelo.
• Combinaciones hidrológicas de suelo y vegetación.
• Condiciones de humedad antecedente a la lluvia de diseño.
La aplicación gráfica de este parámetro se presenta en la figura enseguida:
Página 98 de 156
Figura 60. Solución gráfica para determinar el valor de la escorrentía directa a partir de la lluvia, utilizando el número
de curva CN
Desde el punto de vista de su capacidad de infiltración y para definir el valor CN junto con los otros
tres parámetros mencionados, el Soil Conservation Service ha clasificado los suelos en cuatro grupos:
A, B, C y D, ordenados según su potencial de escurrimiento (tasa de infiltración para suelos
“desnudos” luego de un humedecimiento prolongado, profundidad del nivel freático de invierno,
infiltración y permeabilidad del suelo luego del humedecimiento prolongado y profundidad hasta un
estrato de permeabilidad muy lenta):
El tipo de material en cada caso particular se obtiene a partir de la información de un estudio de suelos,
considerando una clasificación que va desde los tipos “A” de bajo potencial de escorrentía y con altas
infiltraciones (suelos arenosos) hasta los de alto potencial de escorrentía o suelos con infiltración muy
lenta cuando están muy húmedos (arcillosos) tipos “C”. A continuación, se presenta la clasificación
establecida por el S.C.S.
“A”. (Potencial de escorrentía mínimo). Suelos que tienen alta tasa de infiltración aun cuando estén
muy húmedos. Consistentes en arenas o gravas profundas bien drenadas o excesivamente drenadas
con poco limo y arcilla, lo mismo que loess muy permeables. Esos suelos tienen una alta rata de
transmisión de agua.
“B”. (Moderadamente bajo potencial de escorrentía). Suelos con tasas de infiltración moderadas
cuando están muy húmedos. Suelos moderadamente profundos a profundos, moderadamente bien
Página 99 de 156
drenados a bien drenados, suelos con texturas moderadamente finas a moderadamente gruesas y
permeabilidad moderadamente lenta a moderadamente rápida. Son suelos con tasas de transmisión
de agua moderadas. La mayor parte de los suelos arenosos son menos profundos o menos compactos
que los del grupo A. El grupo, en conjunto, posee una infiltración media superior después de haberse
mojado completamente.
“C”. (Moderadamente alto potencial de escorrentía). Suelos con infiltración lenta cuando están muy
húmedos. Consisten en suelos con estratos que impiden el movimiento del agua hacia abajo; suelos
de textura moderadamente fina a fina; suelos con infiltración lenta debido a sales o suelos con lentes
de agua moderada. Esos suelos pueden ser pobremente drenados, o bien, moderadamente drenados
con estratos de permeabilidad lenta a muy lenta. Comprende los suelos poco profundos y los que
contienen mucha arcilla y coloides, aunque menos que los del grupo D. El grupo tiene una infiltración
inferior a la promedia después de saturación.
“D”. (Potencial de escorrentía máximo). Suelos con infiltración muy lenta cuando están muy húmedos.
Son suelos arcillosos con alto potencial de expansión; suelos con nivel freático alto permanente;
suelos con “claypan” o estrato arcilloso superficial; suelos con infiltración muy lenta debido a sales o
alkali y suelos poco profundos ubicados sobre materiales casi impermeables. Estos suelos tienen una
tasa de transmisión de agua muy lenta. Incluyen a la mayor parte de arcillas expansivas y algunos
suelos poco profundos con subhorizontes casi impermeables cerca de la superficie.
La clasificación de los suelos en grupos hidrológicos principales se puede realizar tomando como
base el tipo hidrológico a partir de planos de suelos agrícolas. Si no se dispone de esos planos, se
debe hacer uso de un buen criterio. Usualmente en una misma cuenca se reportan diferentes tipos
de suelos. En ese caso se pueden considerar por separado, indicando el área ocupada por cada uno
de ellos y determinar el valor CN individualmente para luego ponderarlos y obtener un CN
representativo.
Para efectuar una buena clasificación y definir la respuesta hidrológica de los diferentes suelos, en el
documento original del S.C.S. se presenta una lista extensa de suelos agrológicos pertenecientes a
los diferentes grupos hidrológicos. En nuestro medio se deben utilizar criterios más generales, dado
que aún no se han elaborados listas de este tipo. Por lo anterior, lo usual es aplicar criterios
conservadores para compensar la falta de información y evitar errores por defecto.
La clasificación anterior utiliza la premisa de que los suelos con perfiles semejantes en lo que se refiere
a su espesor, textura, estructura, contenido de materia orgánica y expansividad por saturación,
producen una respuesta semejante ante un aguacero de duración considerable. En la comparación
se supone que los suelos están desnudos y que ya han presentado el máximo cambio de volumen por
humedad y que la lluvia es superior a la infiltración potencial.
El segundo parámetro a considerar para definir el valor CN, posee dos componentes: el uso y el
tratamiento. Estos dos parámetros se clasifican en relación con las crecientes que puede producir el
escurrimiento. Mientras mayor sea la retención total ocasionada por un determinado uso y tratamiento,
menor será el caudal pico del hidrograma. El uso del suelo se refiere al tipo de vegetación existente,
presencia de cultivos y usos no agrícolas tales como cuerpos de agua representados por lagunas,
ciénagas y otros. El tratamiento del suelo se aplica fundamentalmente a los usos agrícolas incluyendo
todo tipo de laboreo con rastrillos, nivelaciones, terraceos, control de pastoreo y rotación de cultivos.
Para definir el valor CN, junto con los demás parámetros, se diferencia entre suelos cultivados, pastos
y bosques. En los primeros se identifican tierras en descanso, pastos, cultivos en hilera, granos,
cultivos de alta densidad por la escasa distancia entre ellos, rotaciones, cultivos a lo largo de las curvas
de nivel y cultivos en terrazas.
Rotación de cultivos: La sucesión de los cultivos en una cuenca posee ciertos efectos hidrológicos.
Las rotaciones pueden ser malas o buenas, en función de la cantidad de vegetación densa en la
rotación. Las malas rotaciones son las que se hacen plantando un cultivo en hileras o de granos
pequeños todos los años, consecutivamente. Por otro lado, una mala rotación puede combinar de
varias maneras dos cultivos, granos pequeños o dejar en barbecho durante varias temporadas. Las
buenas rotaciones contienen alfalfa, legumbres que se siembran muy juntas o pastos para mejorar la
textura de la tierra y aumentar la infiltración. Por ejemplo, una rotación de 2 años de trigo (o maíz) y
barbecho puede ser buena para las cosechas cuando existen lluvias reducidas y actúan como un
factor limitante, aunque hidrológicamente es una mala rotación.
Para las praderas naturales y pastizales: Se consideran tres condiciones, que se basan en
consideraciones hidrológicas más que en la producción de forraje., Las praderas malas o pastizales
tienen exceso de pastoreo, o tienen una cubierta vegetal en menos del 50% del área. Las praderas o
pastizales regulares tienen una cubierta vegetal entre el 50 y el 75% del área y no tienen exceso de
pastoreo. Las praderas buenas o pastizales tienen 75% de cubierta con pastoreo ligero.
Los lotes de bosque se clasifican con respecto a los factores hidrológicos y no por la producción de
madera. Los lotes de bosque malos son los que presentan pastoreo excesivo, los que se queman
regularmente, lo que destruye la protección superficial, y los que poseen árboles pequeños. Los lotes
de bosque regulares, poseen algo de pastoreo, pero no se queman; estos lotes pueden tener algo de
cubrimiento superficial, pero generalmente no están protegidos. Los lotes de bosque bueno están
protegidos contra el pastoreo, de manera que el suelo está cubierto superficialmente y posee
arbustos. Para el bosque comercial la clasificación hidrológica se determina en función del espesor y
calidad del cubrimiento superficial, del humus y de la compacidad del humus.
En lo referente a tipos varios de terrenos, se sabe que, generalmente, las cubiertas de las casas, las
carreteras y las áreas urbanas ocupan sólo partes muy pequeñas de la cuenca (en las zonas rurales).
En estos casos, las áreas se consideran dentro de otros tipos de cubiertas debidas a diferentes usos
En la Tabla 32 se combinan los grupos de suelos, el uso del suelo y las clases de tratamiento formando
complejos hidrológicos del suelo y la vegetación. Esta tabla se utiliza para seleccionar el número de
curva CN en función de los parámetros anteriores y considerando una determinada humedad
antecedente, como se mostrará adelante, Los números muestran el valor relativo de los complejos
como productores directos de escurrimiento. Cuanto más elevado sea el número, mayor es el volumen
de escurrimiento directo que puede esperarse de un aguacero. Esta tabla se preparó usando datos
de cuencas aforadas con suelo y vegetación conocidos. El valor de los aguaceros se utilizó para
construir una curva en función del escurrimiento directo para las crecientes anuales y otras avenidas
importantes.
Para seleccionar el valor CN más apropiado en cada caso, se dividió las cuencas de las quebradas en
cuencas homogéneas por concepto de suelos, de vegetación, de usos y de tratamientos y este valor
se determinó individualmente para cada una de ellas. El valor final de CN es el resultado de ponderar
los anteriores, considerando el porcentaje de área ocupada por cada zona homogénea.
Para seleccionar el valor CN se tiene en cuenta el estado de humedad del suelo previo a la lluvia que
se estudia, considerando la cantidad de lluvia caída en el período de los cinco (5) días anteriores,
estableciéndose tres (3) clases de relaciones con dicha cantidad de lluvia. En la Tabla 40 se presentan
las condiciones y los valores correspondientes. Las condiciones medias quedan incluidas en la
condición II, que corresponde al cuadro de curvas de escorrentía para los complejos suelo cobertura
(CN) que se presenta más adelante.
La probabilidad de que la lluvia conjunta de los cinco días anteriores a un aguacero fuerte supere los
53,3 mm es baja y, además, se considera que aún en ese caso, la saturación no es completa y no
alcanza a impermeabilizar completamente la superficie de la cuenca. Por ese motivo se aplica la
condición II, la cual es media en lo referente a su capacidad de escorrentía directa.
(𝑃𝑒 )2 Ec. 12
𝑆′ = − 𝑃𝑒
𝑄
Dado que la escorrentía acumulada “Q” es la misma lluvia efectiva, inicialmente se determina la
infiltración potencial “S” por medio de la siguiente expresión, para lo cual se debe conocer el número
de curva CN.
1.000 Ec. 13
𝑆= − 10,0
𝐶𝑁
La siguiente tabla se utiliza para determinar los valores de CN en función de las prácticas agrícolas y
condición hidrológica promedia (condición II).
Para la condición hidrológica se hace referencia al grado de cobertura vegetal y puede establecerse
de la siguiente manera:
Al determinar el valor CN, se debe considerar la pendiente media de la cuenca, diferenciando entre
alta, media o baja. Después de haber determinado el valor “CN” para la condición II, se puede definir
el correspondiente “CN” para la condición I o para la condición III, de acuerdo con las correlaciones
de la Tabla 34:
Una vez establecidas estas propiedades, es posible determinar el número CN siguiendo los cuadros
guía de esta metodología. El número CN permite, entonces establecer la escorrentía directa a partir
de las siguientes expresiones:
Dónde:
P: Profundidad de precipitación en cm.
Pe: Profundidad de exceso de precipitación en cm.
Debido a las distintas condiciones del suelo y coberturas, se puede ponderar el valor de CN respecto
al área:
𝐹𝑎 𝑃𝑒 Ec. 16
=
𝑆 𝑃 − 𝐼𝑎
Donde Ia se refiere a las pérdidas iniciales antes del encharcamiento, que a su vez se definen como:
25400 Ec. 17
𝐼𝑎 = 0.2 · 𝑆 = 0.2 · ( − 254)
𝑁𝐶
El tiempo de concentración se define como el tiempo mínimo en el cual la totalidad de la cuenca está
aportando flujo de escorrentía al punto de descarga, determinado por el tiempo transcurrido desde
que el agua que se encuentra en el punto más alejado llega hasta el punto de control.
𝑡𝑙 = 0.6 · 𝑡𝑐 Ec. 18
El SCS de los Estados Unidos propone usar el hidrograma unitario adimensional de la figura mostrada
a continuación, basado en la combinación de hidrogramas unitarios en diferentes cuencas en Estados
A
q p = 0.208 ∗ t = √t c + 0.6 ∗ t c Ec. 19
tp p
La forma del hidrograma unitario queda definida al multiplicar los valores de las ordenadas y las
abscisas, que aparecen en la tabla del modelo, por q p y t p. Para la elaboración del hidrograma se
escoge un valor de t/tp y con la tabla se obtiene q/qp. Una vez obtenido este valor y conociendo qp
se halla el valor de q, de la relación t/tp, conocido tp se despeja el valor de t. Este procedimiento se
repite las veces que sea necesario para hallar los valores de q y t, con lo cual se puede obtener la
forma del hidrograma unitario.
Con los datos obtenidos se realizó la modelación hidrológica por medio del software HEC-HMS para
obtener los hidrogramas para los diferentes períodos de retorno, analizando el comportamiento de
4.5.7. Calibración del modelo lluvia – escorrentía para validar el valor de las crecientes
calculadas con el método S.C.S.
Debido a que el modelo lluvia-caudal mediante el cual se calcularon las crecientes en todas las
microcuencas, se adelanta su calibración por medio de un análisis regional utilizando la Adaptación
del Método L-Moments para la Regionalización de Eventos Máximos para las Cuencas de Colombia
(Núñez Galeano, 2015).
La Adaptación del Método L-Moments para la Regionalización de Eventos Máximos para las Cuencas
de Colombia, es el resultado del análisis de los caudales máximos instantáneos en cada una de las
regiones homogéneas encontradas en la investigación. Con el estudio se obtuvieron curvas que
presentan las crecientes con diferentes períodos de retorno en función del área de la cuenca.
En la Figura 62 siguiente aparecen las envolventes de crecientes máximas en función del área de
drenaje para la región 16 (Parte alta Caquetá, Caguán, Putumayo) y aplicadas a crecientes con un
período de retorno de 100 años.
Como puede notarse en el nomograma anterior el caudal determinado por el método hidrológico
(1035 m3/s) S.C.S. es muy similar al generado con el cruce de los datos de las áreas (1100 m3/s), por
tanto, puede considerarse que el modelo hidrológico está acorde con el principio de calibración
regional.
La aplicación del modelo lluvia-caudal para definir las crecientes se inicia con un estudio detallado de
los aguaceros que producen los caudales máximos y se complementa con una caracterización
completa de las cuencas aferentes para luego culminar con la construcción de los hidrogramas
correspondientes.
El río San Juan, se puede considerar como un sistema lótico con modificaciones considerables en su
morfología, teniendo en cuenta que ha tenido alteraciones morfológicas considerables en su forma
natural del cauce permanente por la ausencia de vegetación en diferentes etapas, y por la intervención
antrópica en todas sus expresiones.
El cálculo de las lluvias máximas con diferentes períodos de retorno, se efectuó mediante una hoja
electrónica en Excel (Anexo - Hidrología), la cual, contiene todo el procedimiento descrito
anteriormente en el numeral 4.5.
Para cada distribución de las lluvias en cada una de las subcuencas, se aplicó el principio de áreas
aferentes por índice de precipitación individual de cada polígono de Thiessen, es decir, se realizó una
operación SIG de intersección, donde se cruzó cada una de las subcuencas con los polígonos de
Thiessen y de distribuyó la intensidad en función del área, encontrando así el valor de precipitación
asociado al periodo de retorno de 100 años para agregar al modelo hidrológico. En la Tabla 35 se
presentan los resultados de frecuencia de las precipitaciones máximas.
Tabla 35. Precipitación Máxima en 24 Horas para cada una de las subcuencas
AREA ÁREA % DE
SUBCU NOMBRE FACTOR DE PRECIPITACIÓN PRECIPITACIÓN
TOTAL INTERSECADA PRECIPITA
ENCA ESTACIÓN DISTRIBUCIÓN POR TC (MM) SUBCUENCA (MM)
(KM2) (KM2) CIÓN
BALSAYACO
C1 140.33 63.88 0.46 20.15 9.17
[47010020]
CHUNGACASPI
C1 140.33 52.15 0.37 23.69 8.80
[47010050]
EL PEPINO - AUT
C1 140.33 1.47 0.01 39.51 0.41
[44015070]
PUERTO 25.23
C1 140.33 CAICEDO 0.60 0.00 39.61 0.17
[47010110]
TORRE TV SAN
C1 140.33 4.22 0.03 32.97 0.99
FCO [47010180]
VILLAGARZON
C1 140.33 18.02 0.13 44.25 5.68
[44015010]
VILLAGARZON
C10 18.07 18.07 1.00 52.57 52.57 52.57
[44015010]
BALSAYACO
C11 54.28 40.21 0.74 44.94 33.29
[47010020]
43.94
SANTA ISABEL
C11 54.28 14.08 0.26 41.07 10.65
[47010230]
BALSAYACO
C12 51.48 0.92 0.02 19.79 0.35
[47010020]
PUERTO
C12 51.48 CAICEDO 24.83 0.48 38.90 18.76 28.16
[47010110]
SANTA ISABEL
C12 51.48 25.73 0.50 18.09 9.04
[47010230]
PUERTO
C13 8.78 CAICEDO 8.60 0.98 66.98 65.57
[47010110] 67.14
VILLAGARZON
C13 8.78 0.18 0.02 74.83 1.57
[44015010]
PUERTO
C14 6.22 CAICEDO 1.63 0.26 69.74 18.29
[47010110] 75.77
VILLAGARZON
C14 6.22 4.59 0.74 77.92 57.48
[44015010]
En general, el cálculo del hidrograma de la creciente posee un carácter aproximado. Por tratarse de
un modelo lluvia-caudal, inicialmente se determina el valor de la lluvia de diseño, su duración y su
distribución en el tiempo y en el espacio. A partir de ella se calcula la lluvia efectiva y, por último, se
construye el hidrograma de la creciente producida, considerando las características de la cuenca y
aplicando las propiedades del hidrograma unitario.
1. Se calcula el área de la cuenca considerada. Su determinación se efectuó por medio del programa
ArcGIS. Para valorar el área se delimitó en dos subcuencas tomando en consideración las líneas
divisorias.
2. Se utiliza como insumo el valor de las lluvias máximas en 24 horas con diferentes períodos de
retorno tal como se muestra en la Tabla 35.
3. La duración de la lluvia de diseño, como se trata de una cuenca relativamente grande, se adopta
una lluvia con una duración de 6,0 horas, la cual, fue calibrado el método de L-Moments para las
cuencas hidrográficas de Colombia.
4. A partir del área de la cuenca, se aplica un factor de reducción de 0,85 para representar la
distribución espacial de la lluvia en la zona estudiada. Para seleccionar este valor se parte del
concepto que una lluvia puntual abarca hasta unos 25 km2, la cual, fue calibrado el método de L-
Moments para las cuencas hidrográficas de Colombia.
1.00
0.90
0.80
0.70
0.60
Lluvia[%]
0.50
0.40
0.30
0.20
0.10
0.00
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
Tiempo [%]
6. Se selecciona el intervalo “D” para todo el período del aguacero menor a la quinta parte del tiempo
de concentración, Así, la lluvia anteriormente seleccionada se divide en 10 intervalos de tiempo
de duración “D” cada uno.
7. Se definen los valores secuenciales del aguacero de diseño. Inicialmente se seleccionan los
incrementos de tiempo “D”, en horas; se registran los porcentajes de cada incremento respecto
a la duración de la lluvia total, se escriben los incrementos de la lluvia frente a cada incremento
de tiempo, en mm y se calculan los valores acumulados de la lluvia, en mm. En las siguientes 23
Figuras, se presenta los hietogramas estimados para la estación tenida en cuenta en el modelo
hidrológico de la cuenca.
14
12
Precipitación (mm)
10
0
5
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
15
25
35
45
55
65
75
85
95
Tiempo (min)
14
12
Precipitación (mm)
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
12
10
Precipitación (mm)
0
5
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
15
25
35
45
55
65
75
85
95
Tiempo (min)
18
16
14
Precipitación (mm)
12
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
20
Precipitación (mm)
15
10
0
5
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
15
25
35
45
55
65
75
85
95
Tiempo (min)
16
14
12
Precipitación (mm)
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
30
25
Precipitación (mm)
20
15
10
0
5
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
15
25
35
45
55
65
75
85
95
Tiempo (min)
20
Precipitación (mm)
15
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
20
Precipitación (mm)
15
10
0
5
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
15
25
35
45
55
65
75
85
95
Tiempo (min)
25
20
Precipitación (mm)
15
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
20
Precipitación (mm)
15
10
0
5
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
15
25
35
45
55
65
75
85
95
Tiempo (min)
14
12
Precipitación (mm)
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
30
25
Precipitación (mm)
20
15
10
0
5
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
15
25
35
45
55
65
75
85
95
Tiempo (min)
35
30
Precipitación (mm)
25
20
15
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
25
20
Precipitación (mm)
15
10
0
5
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
15
25
35
45
55
65
75
85
95
Tiempo (min)
16
14
12
Precipitación (mm)
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
20
Precipitación (mm)
15
10
0
5
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
15
25
35
45
55
65
75
85
95
Tiempo (min)
30
25
Precipitación (mm)
20
15
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
20
Precipitación (mm)
15
10
0
5
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
15
25
35
45
55
65
75
85
95
Tiempo (min)
20
Precipitación (mm)
15
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
20
Precipitación (mm)
15
10
0
5
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
15
25
35
45
55
65
75
85
95
Tiempo (min)
20
Precipitación (mm)
15
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
35
30
Precipitación (mm)
25
20
15
10
0
205
105
115
125
135
145
155
165
175
185
195
215
225
235
245
255
265
275
285
295
305
315
325
335
345
355
5
15
25
35
45
55
65
75
85
95
Tiempo (min)
0,2083 ⋅ 𝐴 Ec. 20
𝑞𝑃 =
𝑡𝑃
Se tabulan los caudales pico producidos por cada incremento de lluvia efectiva en m 3/s multiplicando
los incrementos de la lluvia efectiva por sus caudales. Se indica la hora en la cual se inicia cada uno
de los hidrogramas incrementales, en horas. Se inicia con la hora cero y, posteriormente, se va
sumando el incremento de tiempo “D”. También se anota la hora en la cual se presenta el caudal pico
de cada uno de los hidrogramas incrementales, lo cual se hace sumando a la hora inicial el tiempo al
pico en horas. Se registra la hora en la cual termina cada uno de los hidrogramas incrementales, para
lo cual se suma a los tiempos iniciales el tiempo base en horas. En el Anexo 2 se muestran los cálculos
de los hidrogramas utilizando el programa HEC-HMS.
16. Se construye el hidrograma final de la creciente resultante de la cuenca, con el siguiente
procedimiento:
• Se dibujan los hidrogramas triangulares de los incrementos. El programa HEC-HMS efectúa
esta operación.
• Se suman las ordenadas de los hidrogramas triangulares en varias verticales para obtener
diferentes puntos que, unidos, originan el hidrograma final. En muchos casos, sólo es
necesario sumar las ordenadas en las horas que representan el principio, el máximo y el final
de cada hidrograma de los incrementos.
Así las cosas, el cálculo de caudales para el proyecto se realiza a la descarga del río Vides y el
principal cauce del río San Juan. El caudal pico a tiempo de retorno de 100 años es de 725.5 m3/s,
mientras que la cuenca integralmente posee un caudal de 1035 m3/s.
5. SIMULACIÓN HIDRÁULICA
HEC-RAS en su versión 6.1, resulta novedosa por su análisis en 2D. Este tipo de modelización ha
permitido que desarrolle nuevas capacidades para la simulación de flujo en 2D; así por ejemplo puede
ejecutar modelos en 1D, 2D o una combinación de ambos. Estas capacidades serán descritas a
continuación (la información fue tomada del manual del usuario de HEC-RAS 2D).
• Puede realizar modelización en 1D, 2D o una combinación de ambos. La capacidad de ejecutar
modelos combinados, permite al usuario trabajar en amplios sistemas de ríos, así, por ejemplo,
se puede usar el modelo 1D sobre el cauce de un río y el modelo 2D sobre las zonas
adyacentes al mismo, en donde se requiere mayor detalle los resultados hidráulicos.
• Emplea las ecuaciones de Saint-Venant y Onda Difusa en 2D. El programa permite elegir entre
las ecuaciones de Saint-Venant u Onda Difusa en 2D para llevar a cabo la modelización. De
forma general las ecuaciones de Onda Difusa en 2D permiten que el software procese la
información rápidamente y tiene mayores propiedades de estabilidad, sin embargo, las
ecuaciones de Saint-Venant 2D son aplicables a una gama más amplia en la resolución de
problemas.
• Esquema numérico empleado: Volumen Finito Implícito. El método de volúmenes finitos,
incrementa la mejora en la estabilidad y robustez sobre las técnicas de diferencias finitas y
elementos finitos; pudiendo manejar adecuadamente situaciones de flujo en regímenes
Subcrítico, Supercrítico y Mixto.
• Algoritmo de solución para el acoplamiento de modelos 1D y 2D. Este algoritmo permite la
retroalimentación directa en cada paso de tiempo entre los elementos de flujo 1D y 2D. Por
ejemplo, el caso del río modelado en 1D que se conecta a cualquiera área (modelada en 2D)
mediante algún dique (estructura lateral). Si consideramos que el flujo se desplaza por encima
del dique, o por efecto de la ruptura del dique, del modelo 1D hacia el 2D, entonces el
programa emplea la ecuación de vertedero para resolver el cómputo del flujo. Por cada paso
de tiempo, la ecuación de vertedero emplea los resultados de los modelos 1D y 2D permitiendo
la contabilización exacta del vertedero sumergido.
• Mallas computacionales estructuradas y no estructuradas. HEC-RAS 5.0 fue diseñado para
trabajar con mallas no estructuradas, pero también puede trabajar con mallas estructuradas.
Las celdas computacionales de una determinada malla pueden ser triángulos, cuadrados,
rectángulos o polígonos de hasta 8 lados como máximo. La malla computacional no necesita
El software requiere como datos de entrada el régimen de flujo, las condiciones de borde, los niveles
de la sección a banca llena y el valor de la rugosidad de Manning, y soluciona la ecuación para flujo
gradualmente variado en canales naturales y artificiales en régimen permanente. Las profundidades
de flujo se calculan a partir del método estándar por pasos y como resultado del modelo se obtiene la
superficie libre del agua asociada a una condición de caudal impuesta que puede variar a lo largo del
cauce, según US Army Corps of Engineers, (2016, citado por IDEAM & Universidad Nacional de
Colombia, 2013).
Los mecanismos de cálculo para el modelo se indican en las siguientes expresiones matemáticas:
𝜕𝐻 Ec. 22
+ ∇ ∗ ℎ𝑉 + 𝑞 = 0
𝜕𝑡
donde V = (u, v) es el vector de velocidad y el operador diferencial del (∇) es el vector de los
operadores de derivada parcial dado por
∇ = ∂ ∂x ∂∂ /, (/ y).
5.1. PARAMETRIZACIÓN
Los modelos hidrodinámicos se clasifican como modelaciones matemáticas, los cuales constituye un
método cognoscitivo en el que el objeto de estudio o cognitivo (río, humedales, quebradas, entre
otros) es reemplazado con otro, llamado modelo (imagen matemática) que cumple con relación al
primero unas condiciones de analogía y semejanza, es decir que hay una equivalencia entre los
resultados del modelo y el tramo del objeto de estudio (Domínguez, 2011), lo anterior se implementa
mediante un algoritmo lógico – numérico, ejecutado por un ordenador, permitiendo así estudiar la
cualidades del proceso original (Domínguez, 2011; Pena, 2011).
Partiendo de lo anterior, los modelos hidrodinámicos guardan esa analogía y semejanza con los
cuerpos de agua superficial tanto loticos como lenticos y que permiten estimar la delimitación de la
amenaza de inundación en áreas adyacentes a los mismos y susceptibles a este fenómeno,
relacionando los caudales circundantes (producto del modelo hidrológico establecido) al eje ecológico
principal y planicies de inundación, que asocian a su vez un campo de alturas de lámina de agua y un
campo de velocidades, (Ashley et al., 2005). Dichos modelos cuando se usan para la gestión del riesgo
de inundación y en general el entendimiento hidrodinámico del tramo de simulación, se establece una
serie de aproximaciones sucesivas mediante la modificación del dominio computacional (diferentes
áreas de inundación), variaciones en el refinamiento de malla y geometría, relacionando lo anterior a
diferentes periodos de retorno y planteamiento de escenarios futuros, asociados con cambios en las
variables ambientales y geométricas del medio. (Ven T Chow, 1994; Ashley et al., 2005; Middelkoop
et al., 2004; Gonzalo, 2011; Pena, 2011).
Los modelos matemáticos pueden clasificarse en primera instancia según tres criterios, (Domínguez,
2011), como se menciona a continuación:
• El primer criterio se base en la capacidad del modelo en tener presente la incertidumbre del
proceso, diferenciándolos en modelos determinista, el cual sus valores de entrada y
parámetros definirán un conjunto finito de valores de salida, asumiendo un control del proceso
Los modelos desarrollados por lo general tienen una relación en la permanencia de los caudales
modelados en función del tiempo (permanente y no permanente), dado la discretización del esquema
numérico que utilicen, dando una versatilidad de los mismos en el segundo criterio de selección
descrito, ya que los análisis pueden simplificar el problema (permanente) o aproximando al mismo, al
momento de una modelación de crecientes, ya que las corrientes con una pendiente y velocidad alta,
como es el caso de las corrientes lóticas o de montaña, pueden conservar la consistencia en pequeños
intervalos de tiempo, (Pena, 2011; Bladé, 2009).
Paralelamente, el tercer criterio de selección es el que mayor incidencia tiene a la hora de diferenciar
este tipo de modelos, en donde el objetivo final es la resolución conjunta de las ecuaciones de
gobierno en 1D o 2D, encontrando el balance entre la lógica del problema a describir, su precisión y
costo.
Teniendo en cuenta que existe una variedad de procesos según el tipo de modelo, la selección
depende de la escala espacial y temporal de la modelación. En cuanto a la escala temporal, entre más
Para desarrollar el proceso de modelación hidráulica, el cual se incluirá: esquematización del sistema
fluvial (cauces y planicie), implementación en el modelo, análisis de sensibilidad, calibración,
verificación y manejo de inundaciones para diferentes escenarios hidrológicos. Una vez calibrado y
verificado el modelo hidráulico, se procede a realizar las diferentes simulaciones que permitirán
generar las condiciones hidrodinámicas en los cauces y en la planicie, especialmente los niveles y
velocidades del agua y la extensión de las áreas inundadas, para los diferentes eventos hidrológicos
asociados a los periodos de retorno preestablecidos. Se va determinar la velocidad, dirección,
profundidad del flujo y las llanuras de inundación. El desarrollo del proceso va a incluir:
5.1.1. Esquematización
Previo a la calibración del modelo, se hace un análisis de sensibilidad que consiste en identificar los
parámetros más sensibles ante cambios pequeños en sus valores; con esto se logra ser más eficiente
en la calibración. En esta etapa se busca encontrar los parámetros que tienen mayor peso o relevancia
en la calibración del modelo (IDEAM, 2018). Los cuales se realizarán las pruebas con cambios en los
datos de entrada, y establecer qué tanto cambia los niveles de agua, las profundidades, las
velocidades del flujo y los caudales que te arrojando el modelo hidráulico.
Una vez calibrado y verificado el modelo hidráulico, se procede a realizar las diferentes simulaciones
que permitirán generar las condiciones hidrodinámicas en los cauces y en la planicie, especialmente
los niveles y velocidades del agua y la extensión de las áreas inundadas, para los diferentes eventos
hidrológicos asociados a los periodos de retorno preestablecidos. Para cada periodo de retorno se
introducen al modelo las hidrógrafas de caudales correspondientes.
También se debe realizar la modelación hidráulica entre los empalmes de dos caudales sean de ríos
o canales, considerando el efecto de remanso en dicho empalme, esto con el objeto de además de
considerarlo en la modelación, sirva para aportar la información necesaria en las recomendaciones
finales.
Se realizará la definición de las condiciones de frontera para las cuencas en estudio, en especial la
condición de frontera en los diferentes sitios de entrega de los afluentes al río o quebrada,
considerando la pendiente de fondo del cauce. Por otro lado, en las fronteras de inicio se deberá
definir los caudales que llegan a estos puntos de acuerdo a los resultados del modelo hidrológico.
Para desarrollar la modelación hidráulica del río San Juan fue necesario desarrollar ciertos pasos
metodológicos con el propósito de obtener las variables de diseño hidráulico del puente a construir.
A continuación, se realiza una descripción detallada de los pasos metodológicos desarrollados para
realizar la modelación hidráulica.
Para poder calcular los patrones de flujo en el dominio 1D se hace necesario la discretización de esta
área en elementos más pequeños.
Se ha digitalizado las bancas del canal a distancia de 2 metros de cada lado, estos valores de distancia
fueron seleccionados teniendo en cuenta que un espaciado más pequeño sin detrimento a los altos
costos computacionales, éste favorecerá al proceso de calibración y verificación del modelo
hidrodinámico, el cual consume tiempo y está íntimamente ligado a los tiempos de cálculo del
escenario planteado.
Para alimentar el modelo, se fijaron una serie de puntos a lo largo del dominio 1D por donde entraría
el flujo de agua proveniente de la parte alta. En la siguiente figura se presentan todas las condiciones
de borde e iniciales del modelo hidráulico para cada uno de los escenarios propuestos.
Para el río objeto de estudio, los factores que más influencia tienen en la determinación del valor de
la n de Manning son: por la textura del lecho, por efectos del relieve del lecho creado por las
microformas, por el tamaño y a la forma de la sección transversal, por las irregularidades de la sección
transversal, por las variaciones de tamaño y forma de la sección, asociado a las obstrucciones al flujo,
por la vegetación, relacionado con el transporte de sedimentos y por la presencia de curvas o
meandros en planta y las características morfológicas del mismo. Teniendo en cuenta estos
condicionantes, disponemos de un método de obtención del valor del coeficiente de rugosidad de
Manning, se utilizó el método de Cowan, la cual, es más representativo para todo lo largo del cauce.
Así, con este método, partiendo de un canal teórico, uniforme y liso de un material dado al cual le
corresponde un valor de coeficiente de rugosidad de Manning determinado, a éste se le van
añadiendo valores que representan la presencia de elementos y características que condicionan el
paso del flujo.
Evidentemente, existen unas tablas orientativas de rangos de valores de cada uno de esos factores
de n de Manning, y que nos ayudaron a poder determinar el valor de coeficiente de rugosidad que
tendrá nuestro cauce para el tramo de estudio.
Con el fin de representar lo mejor posible la resistencia hidráulica del río se adelantará un trabajo de
campo consistente en analizar la granulometría del material del lecho, los posibles obstáculos en el
propio cauce, las condiciones de las laderas, la presencia de vegetación y todos los demás parámetros
que influyen en este parámetro, el cual se expresa mediante el coeficiente de rugosidad de Manning
y de la observación de una imagen de satélite geo-referenciada, de los apuntes y fotos de campo,
n0 - Valor de n para un río recto y prismático que considera sólo la textura del lecho
n1 - Valor de n por efectos del relieve del lecho creado por las microformas
n2 - Valor de n debido al tamaño y a la forma de la sección transversal
n3 - Valor de n debido a las irregularidades de la sección transversal
n4 - Valor de n debido a las variaciones de tamaño y forma de la sección
n5 - Valor de n asociado a las obstrucciones al flujo
n6 - Valor de n por vegetación
m1 - Factor correctivo relacionado con el transporte de sedimentos
m2 - Factor correctivo por la presencia de curvas o meandros en planta
Para el propósito de estudio se seleccionó para el cauce un Manning de 0.028, para la solera derecha
de 0.031 y para la solera izquierda de 0.02.
Esta categoría simplifica el alcance y zonas de desborde la inundación sobre el área de análisis.
De esta categoría puede indicarse que completamente em todos los flancos del sector de análisis la
corriente supera los niveles de contención naturales en el paso actual. El estribo y muerto derecho es
completamente cubierto por la masa de agua, mientras que anclaje izquierdo no llega a ser
intervenido.
El esquema de velocidades permite reconocer como el río se encuentra interviniendo en los procesos
morfohidráulicos de la corriente, principalmente en los puntos de agregación y socavación, y en una
medida menor, donde existe pérdida de la vegetación nativa que actúa como salvamento del material
deposicional.
El esquema de velocidad indica que el punto máximo corresponde a 4.42 m/s, marco en el punto de
activación de torrente por desembocadura del río Vides, es decir 180 metros antes de la obra
proyectada. Bajo la estructura propuesta la velocidad se encuentra catalogada en 1.4 m/s, la cual es
de condiciones aceptables para el diseño hidráulico.
Las velocidades más bajas se encuentran entre los 0.4 y 0.7 m/s, y se presentan en la rampa-remanso
que existe posterior al ponteadero.
El Tirante hidráulico o distancia de flujo es la distancia vertical desde el punto más bajo de la base del
canal del río hasta la superficie libre del líquido transportado. Es el principal factor de análisis a evaluar,
puesto que éste condiciona la altura final que deben tener el puente.
Dicha evaluación arroja que, entre el puente propuesto y el Nivel de Aguas Máximas Extraordinarias
– NAME, existe una distancia de 1.86 m.
300
Legend
WS PF 1
Crit PF 1
290
0.0 m/s
0.2 m/s
0.4 m/s
280
Elevation (m)
0.6 m/s
0.8 m/s
1.0 m/s
270
Ground
Bank Sta
260
250
0 20 40 60 80 100 120 140 160
Station (m)
Figura 94. Esquema del NAME vs Puente propuesto - Vista aguas abajo
SanJuan
RS = 73 BR
310
Legend
WS PF 1
300 Crit PF 1
0.0 m/s
0.2 m/s
290
0.4 m/s
Elevation (m)
0.6 m/s
280 0.8 m/s
1.0 m/s
Ground
270
Bank Sta
260
250
0 20 40 60 80 100 120 140 160
Station (m)
Para estimar el descenso que se puede producir en el lecho del río San Juan en la zona de la autopista
hidráulica -antes y contacto de pontón- se ha tomado como referencia la metodología de socavación
general de Lischtvan – Levediev (Guitelman & Pérez, 1999), la cual tiene la siguiente definición
matemática:
Velocidad erosiva que es la velocidad media que se requiere para degradar el fondo está dado por
las siguientes expresiones:
Donde,
Hs = Tirante considerado, a cuya profundidad se desea conocer que valor de Ve se requiere para
arrastrar y levantar al material (m)
dm = Es el diámetro medio (mm) de los granos del fondo obtenido según la expresión.
dm = 0.01 S di pi
En el cual,
di = Diámetro medio, en mm, de una fracción en la curva granulométrica de la muestra total que se
analiza.
pi = Peso de esa misma porción, comparada respecto al peso total de la muestra. Las fracciones
escogidas no deben ser iguales entre sí.
Suelos cohesivos
1
𝛼 ∙ 𝐻05⁄3 1+𝑥 Ec. 26
𝐻𝑠 = ( )
0,60 ∙ 𝑔𝑑 1.18
Suelos no cohesivos
1
𝛼 ∙ 𝐻05⁄3 1+𝑥 Ec. 27
𝐻𝑠 = ( )
0,68 ∙ 𝑑𝑚 0.18
Donde,
A = Qd / (Hm5/3 Be m)
Qd = caudal de diseño (m3/seg)
Be = ancho efectivo de la superficie del líquido en la sección transversal
m = coeficiente de contracción.
Hm = profundidad media de la sección = Área / Be
x = exponente variable que depende del diámetro del material
dm = diámetro medio (mm)
Teniendo en consideración las variables que implican el cálculo de la socavación, siendo esta
particular donde la construcción de los cimientos del puente se realiza sobre anclaje sobre suelo
aluvial reciente y terraza reciente, se ha tenido en cuenta el material caracterizado en la fracción
Aplicando la metodología descrita, en ejecución del módulo de diseño hidráulico del puente, se
obtiene los siguientes resultados de la socavación:
De acuerdo con lo anterior, se puede indicar que en 100 años la margen del pilar derecho tendrá una
socavación de 3.61 m, mientras que la izquierda será de 0.45 m. Por su parte, el canal del río por
efecto de contracción socavará 1.77 m.
Figura 95. Esquema frontal de socavación lateral y general de la estructura del puente
Como la condición del cauce no es torrencial, y los efectos de perdida lateral de suelo se manifiestan
mediante desplomes verticales, se propone intervención en los puntos críticos señalados como
complejos en la geomorfología fluvial.
De manera integral se debe proteger las zonas de estribos mediante cubierta con Hidromalla
Geoestera (Nombre comercial casa Geomatrix), esta protección deberá de realizarse con un manto
que cubra 10 metros a cada torre se deberá perfilar el terreno a máximo 12° y acuñar el proceso
constructivo sugerido por el diseñador geotécnico.
La protección posee diversas tipologías de acuerdo al elemento que se desea resguardar, los
esquemas típicos de arreglo en las mallas son cuatro, los cuales se exponen enseguida:
Los criterios de diseño que debe satisfacer la infraestructura son los siguientes:
a) Para los efectos que representa la correcta localización del puente, se deberá realizar
amojonamiento estático de puntos geodésicos. Las coordenadas centroides para la
localización de la estructura son las siguientes:
b) El constructor deberá aplicar obra de protección anti-socavación local, para lo cual aquí se
recomienda el uso de geomallas, el esquema pictórico de implante se encuentra en el plano
de diseño hidráulico.
e) En ambas márgenes deberá fortalecerse con la siembra de material vegetal dominante y nativo
de las orillas, principalmente los citados en la evaluación ecosistémica. Para ello se
recomienda un espaciado de siembra de 6x6 m.
f) Las actividades de descolmatación deberán realizarse por lo menos cada 1.5 a 2 años, o cada
ocasión que la fuente realice acumulaciones significativas en el área diseñada de intervención
ésta descolmatación deberá ser manual, puesto que por las condiciones de acceso no es
posible desarrollarse mediante maquinaria amarilla.
g) La localización del puente es sugerida por criterios geomorfológicos e hidráulicos, por tal
razón, se deberá chequear las condiciones de suelo y viento por parte de los especialistas
específicos para validar el esquema sugerido.
i) El gestor deberá elaborar el plan de inversión del 1% de la obra como mecanismo ambiental
para obtener el permiso de ocupación, y en ese contexto, se propone que se desarrolle
procesos de revegetarización del bosque de galería.
BIBLIOGRAFÍA
Ackerman, C. T., Wre, D., Jensen, M. R., & Brunner, G. W. (2010). Geospatial capabilities of HEC-RAS for
model development and mapping. Second Joint Federal Interagency Conference.
Armenta Porras, G. E. (2013). Análisis detallado del efecto foehn generado por la cordillera Oriental en el
alto Magdalena (Huila y Tolima). 100. Retrieved from http://www.bdigital.unal.edu.co/11389/
Ballesteros, J., Bodoque, J. ., Díez-Herrero, A., Génova, M., Gutiérrez, E., Moya, J., … Saz, M. . (2017).
Dendrogeomorfología Los árboles, fuente del conocimiento de los procesos y desastres naturales
(IPL). La Pobla Llarga.
Bateman, A. (2007). Hidrologia basica y aplicada. Grupo de Investigación En Transporte de Sedimetos, 70.
Calvache-Fajardo, J., CISFA S.A.S, & Emerald Energy Colombia. Estudio de régimen hídrico de la
subcuenca hidrográfica del Río Pescado. , (2016).
Calvache-Fajardo, J., & Cuesvas, E. (2016). MODELACIÓN DE MÁXIMA AVENIDA DEL RÍO PUTUMAYO
PARA LA IDENTIFICACIÓN DE LA AMENAZA POR INUNDACIÓN DE LA POBLACIÓN DE PUERTO
CAICEDO, DEPARTAMENTO DEL PUTUMAYO. Universidad de Manizales.
Carvajal, J. H., & Servicio Geológico Colombiano. (2011). Propuesta de Estandarización de la Cartografía
Geomorfológica en Colombia. Ingeominas, 71. Retrieved from
https://www2.sgc.gov.co/Publicaciones/Cientificas/NoSeriadas/Documents/Propuesta de
estandarización.PDF
Chow, V., Maidment, D., & Mays, L. (1994). Hidrología aplicada. Hidrologia Aplicada, p. 575 pp. Retrieved
from http://bases.bireme.br/cgi-
bin/wxislind.exe/iah/online/?IsisScript=iah/iah.xis&src=google&base=REPIDISCA&lang=p&nextActio
n=lnk&exprSearch=158911&indexSearch=ID%5Cnhttp://www.sidalc.net/cgi-
bin/wxis.exe/?IsisScript=BINAI.xis&method=post&formato=2&cantidad=
Downward, S. R., & Brookes, A. (1994). A methodology for quantifying river channel planform change using
GIS. (224), 449–456.
Escuela superior de administración Pública - ESAP. (2012). Características Morfométricas de Cuencas
Hidrográficas. Centro de Documentación e Información Municipal, 4–6.
Frissell, C. A., Liss, W. J., Warren, C. E., & Hurley, M. D. (2014). A Hierarchical Framework for Stream Habitat
Classification: Viewing Streams in a Watershed Context. (March 1986).
https://doi.org/10.1007/BF01867358
García, D. G., Naverac, V. A., Bea, E. D., Ibisate, A., Matauco, G. De, & Ojeda, O. (2011). Análisis diacrónico
de la migración de cauces fluviales mediante técnicas de SIG vectorial. (2006), 18–20.