Puente San Juan-Vides V1 1

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 156

ESTUDIOS HIDROLÓGICOS E HIDRÁULICOS PARA

LA CONSTRUCCIÓN DEL PUENTE PEATONAL


SOBRE EL RÍO SAN JUAN, MUNICIPIOS DE PUERTO
CAICEDO Y ORITO, DEPARTAMENTO DEL
PUTUMAYO

I.A. Esp SIG Jhimmy CALVACHE FAJARDO


M.P. 00000-24847 AGR
CC 18130245 de Mocoa – Putumayo

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

El presente documento contiene el estudio de hidrología e hidráulica con el fin de proyectar la


construcción del puente peatonal a localizar en el paso de del río San Juan en las veredas San Juan
Vides y Damasco Vides en los municipios de Orito y Puerto Caicedo, respectivamente. El modelo de
informe se orienta de manera que se reconozca las condiciones físico-ambientales presentes en el
sector para apropiar los modos y métodos constructivos óptimos en la implantación de obra para el
diseño de la obra. La metodología del estudio presentado en este informe consiste en estimar la
intensidad, tiempo de duración y frecuencia de la precipitación dentro del área donde las estaciones
meteorológicas de IDEAM hacen referencia espacial. Del mismo modo, con los datos obtenidos de
caudal en el modelo hidrológico se realiza la simulación hidráulica de transporte de fluido en el vado
mediante el método de energía en presentación unidimensional. En el contexto de la localización, se
realizó el componente de dinámica fluvial por estudio multitemporal de la variación del componente
sedimentos y alineamiento de corriente.

El contenido empleado para determinar la condición geomorfológica fluvial se ha tenido en cuenta la


metodología de (Carvajal & Servicio Geológico Colombiano, 2011), en lo referente al componente
hidrológico e hidráulico se ha tomado el Protocolo de modelación hidrológica e hidráulica de (IDEAM,
2018) y el Código Colombiano de Puentes CCP-12 del (Ministerio de Transporte & AIS. INVIAS, 2014).

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.

1.2. OBJETIVOS ESPECÍFICOS

- Consultar, procesar y presentar la información cartográfica e hidrológica relacionada con el


alcance del estudio.
- Delimitar, estimar parámetros físicos-forma de la microcuenca de la escorrentía
- Calcular las lluvias para un periodo de retorno de diseño a partir de los registros de
precipitación máxima en 24 h con aplicación del modelo SCS.

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:

Figura 1. Localización general del área de proyecto

Fuente: Elaboración propia, 2022

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.

La obra se localiza al nor-occidente de la cabecera municipal de Puerto Caicedo, en una distancia


aproximada de 15 km de recorrido fluvial. Para llegar al sitio se parte del puerto fluvial de Puerto
Caicedo, donde se toma embarcación motorizada por aguas arriba del río San Juan, el cual, en un
transcurso de 15 km, se encuentra con la confluencia del río Vides, lugar donde así mismo se localiza
la población de San Juan Vides, los cuales serían los principales beneficiarios de la obra. El tránsito
por el río debe realizarse con embarcación de bajo calado (8 a 12 m) y con carga leve, pues existen
lugares que por procesos de imbricación y sobrelevación.

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

La geología es una de las características más predominantes al momento de plantear a meso-escala


la localización de obras con características permanentes y que cuya extensión sea de gran
envergadura. Es por ello que dentro del estudio se realizó una evaluación y documentación de las
características litológicas y tectónicas de la cuenca donde se desarrollará la obra.

La verificación compone la verificación en el contexto de la geología regional y las características


adyacentes de la geología tectónica que predomina en la transición de piedemonte a llanura
amazónica.

Página 6 de 156
Figura 2. Geología del área de proyecto

Fuente: Elaboración propia, 2022

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.

2.1.1.1. Grupo Orito (E3N1or)

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.

En algunos sectores se encuentran niveles de conglomerados finos, granosoportados, cuya


composición principal es chert y cuarzo, así como niveles arcillosos oscuros, con abundante materia
orgánica, delgadas láminas de yeso y capas centimétricas de carbón. Algunos horizontes fosilíferos
han sido reportados por diferentes autores (McGirk, 1949; Cáceres & Teatín, 1985; entre otros), en
los cuales se registra la presencia de foraminíferos arenáceos, ostrácodos, gastrópodos y fragmentos
de plantas.

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).

2.1.1.2. Terrazas medias (Qt2)

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.

2.1.1.3. Aluviones (Q2al)

Depósitos aluviales extensos se encuentran entre el piedemonte y la llanura amazónica,


especialmente en el extremo sureste de la Plancha 430 Mocoa, al sur y oriente de la Plancha 431
Piamonte y en el oriente de las planchas 449 Orito y 465 Churuyaco.

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.

El espesor de estas acumulaciones es muy variable y su edad es holocénica a reciente y están en


permanente proceso de acumulación y reorganización.

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

Fuente: (Carvajal & Servicio Geológico Colombiano, 2011)

Inicialmente, el proceso metodológico se basó en la recopilación de la información secundaria


necesaria para la interpretación del componente; ya fuera de tipo bibliográfico o cartográfico.

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.

2.2.1. Condiciones geomorfológicas del área de estudio

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

Fuente: (Instituto de Hidrología Meteorología y Estudios Ambientales, 2010)

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.

2.2.3. Jerarquización geomorfológica

Como se mencionó en la Figura 4, para la realización sistemática del levantamiento y jerarquización


geomorfológica, se tuvo en cuenta la metodología y nomenclatura planteada en la “Propuesta de
Estandarización de la Cartografía Geomorfológica en Colombia”, publicada por el Servicio Geológico
Colombiano en el año 2012.

A partir de la escala de trabajo, la definición de las categorías geomorfológicas, permitió el análisis


sistemático de la evolución geomorfológica del área de estudio. Esta concepción, es igualmente
necesaria con el fin de enfocar el pensamiento, hacia el análisis de las geoformas de manera sistémica
y coherente, con las operaciones funcionales de almacenamiento, consulta, análisis, generación y
generalización de información cartográfica de los Sistemas de Información Geográfica -SIG- (Carvajal,
2012).

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

Se refiere a grandes áreas geográficas o amplios espacios continentales o intra-continentales


caracterizados y definidos por estructuras geológicas y topográficas regionales que han tenido
deformación o basculamiento y posiblemente metamorfismo o intrusión ígnea.

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).

El área de estudio, se encuentra incluida dentro de la Provincia Geomorfológica de la cordillera


Occidental.

• Región geomorfológica (Ambiente morfogenético)

Es la agrupación de geoformas relacionadas genética y geográficamente. Están definidas por los


ambientes morfogenéticos y geológicos afectados por procesos geomórficos parecidos. El ambiente
morfogenético, hace alusión a las condiciones físicas, químicas, bióticas y climáticas bajo las cuales
se generaron las geoformas (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.

De manera general, los ambientes morfogenéticos se agrupan de la siguiente manera: Denudacional:


D; Fluvial y lagunar: F; Estructural: S; Glacial: G; Volcánico: V; Kárstico: K; Marino y costero: M; Eólico:
E; Antropogénico: A.

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.

• Unidad, subunidad y componente geomorfológico

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).

Tabla 3. Jerarquización geomorfológica para el área de estudio


GEOMORFO- PROVINCIA REGIÓN UNIDAD SUBUNIDAD COMPONENTE SÍMBOLO
ESTRUCTUR (AMBIENTE
A MORFOGENÉTICO)
Vertiente Escarpes de longitud
oriental, corta a larga, de forma
SISTEMA cordillera Escarpe de cóncavo convexa y
Ambiente de origen Escarpe de
OROGÉNICO Oriental erosión eventualmente recta, Deeme
Denudacional (D) Erosión
ANDINO (Peneplanacione menor con pendiente
s de la escarpada a muy
Amazonia) escarpada
Cauce activo En forma de U alargada FVAca
Vertiente
oriental,
SISTEMA cordillera Plano de De longitud corta a FTFpt1
Ambiente de origen Llanura de terraza T1 moderada y de forma
OROGÉNICO Oriental
Fluvial (F) inundación recta
ANDINO (Peneplanacione
s de la Plano de De longitud corta a FTFpt2
Amazonia) terraza T2 moderada y de forma
recta

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

2.2.3.1. Cauce activo (FVAca)

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.

2.2.3.2. Plano de terraza (FTFpt1, FTFpt2, FTFpt3)

Corresponde a las geoformas de pendiente ligeramente plana a ligeramente ondulada, de extensión


corta y configuración estrecha. Su origen está relacionado con procesos de erosión y acumulación
aluvial.

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.

Particularmente, los planos de terrazas son remanentes de anteriores niveles de sedimentación, en


los cuales la corriente ha realizado un proceso de incisión o socavación vertical como consecuencia
del rejuvenecimiento del paisaje.

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.

Figura 5. Geomorfología del área de proyecto

Fuente: Elaboración propia, 2022

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.

Tabla 4. Agentes y procesos geomorfológicos externos


PROCESOS GEOMORFOLÓGICOS
FUERZA O AGENTE
Denudativos Constructivos
AGUA, HIELO, TEMPERATURA, Meteorización de las rocas Principio de formación de suelos
ORGANISMOS
GRAVEDAD, MOVIMIENTOS Remoción en masa Sedimentación coluvial
SÍSMICOS
LLUVIA Erosión pluvial Sedimentación coluvial
GRAVEDAD + AGUA DEL SUELO Remoción torrencial, flujos de lodos y Sedimentación diluvial
escombros
ESCORRENTÍA Erosión fluvial Sedimentación aluvial
OLEAJE, CORRIENTES COSTERAS, Erosión marina o lacustre Sedimentación marina o lacustre
MAREAS
NIEVE Y HIELO EN MOVIMIENTO Erosión glaciar Sedimentación glaciar
(GLACIAR)
VIENTO Erosión eólica Sedimentación eólica
HOMBRE Degradación antrópica Conservación y recuperación de
tierras
Fuente: (Villota & Instituto Geográfico Agustin Codazzi (IGAC), 2005)

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.

2.2.4.1. Procesos erosivos

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.

A continuación, se menciona el procedimiento que se tuvo en cuenta para llevar a cabo la


identificación y caracterización de los procesos erosivos que predominan en el área de estudio:

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.

Información primaria: Atendiendo la referencia anterior, muchos de los procesos que no se


identificaron en los sensores remotos, fueron georreferenciados en campo; y aquellos que por tamaño
y área cumplían con los requisitos mínimos, fueron cartografiados.

En el área de estudio, se identificaron procesos erosivos activos, encaminados a la presencia de la


erosión laminar, en surcos y erosión en cárcavas.

• 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.

Figura 6. Morfología fluvial del área de proyecto

Fuente: Elaboración propia, 2022

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.

2.2.4.2. Procesos de movimientos en masa

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.

Información secundaria: Es importante mencionar, que la complementación a la interpretación de


los procesos morfodinámicos del área de estudio, fue aumentada con la información secundaria. El
uso de dicha base de datos, permitió georreferenciar los movimientos en masa relacionados con
Deslizamientos, Propagaciones laterales, Volcamientos, Deformaciones gravitacionales, Flujos,
Reptaciones y Caídas, que por escala y resolución no se identificaron en los sensores remotos.

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.

• Terraceo o pate vaca

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.

• Deslizamientos activos e inactivos

Corresponden a los procesos de remoción en masa violenta y dramática; en donde, se generan


grandes volúmenes de material de suelos residuales que se desprende cuesta abajo.

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

Fuente: Elaboración propia, 2022

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.

2.2.5. Diacronismo fluvial

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.

Figura 8. Parámetros estudiados en el análisis de la migración de los meandros

Fuente: (García et al., 2011)

2.2.5.1. Superficie movilizada

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.

Tabla 5. Superficie movilizada en el área de estudio


AÑO TIPO LONGITUD ANCHO (M) ÁREA (M2)
(M)
1976 Zona Activa 3473.8 2079.6 50.85
1985 Zona Activa 3536.2 1974.5 48.21
1999 Zona Activa 3537.3 1965.6 35.30
2008 Zona Activa 3490.4 1943.6 33.90
2022 Zona Activa 3490.4 1943.6 37.93
Fuente: Elaboración propia, 2022

Página 26 de 156
Figura 9. Dinámica traslacional de la zona activa entre 1976 y 2022

Fuente: Elaboración propia, 2022

2.2.5.2. Dinámica de los meandros

Los resultados obtenidos reflejan el aumento considerable de la magnitud de cambio en el meandro


del área del proyecto. La variación de la distancia es el parámetro que mejor refleja el aumento
progresivo de la dinámica, pero en el caso del proyecto este no es superior a 1 metro anual (promedio)
en la zona de influencia directa. Para entender la dinámica del meandro se ha tomado como referencia
la teoría de Lagasse et al. (2004) donde se reconoce los patrones de la dinámica del meandro por
medición directa de su eje.

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)

Tabla 6. Radios del meandro del proyecto


AÑO RADIO MEANDRO ÁREA CIRCULO
(M) (HA)
1976 417 54.5
1985 283 25.1
1999 256 20.5
2008 159 8.0
2022 175 9.6
Fuente: Elaboración propia, 2022

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

Fuente: Elaboración propia, 2021

El ángulo de movimiento, por su parte es poco variable, incluso la geometría de localización se


conserva a pesar de las condiciones de frontera variable. Bajo esas instancias donde los ángulos
resultan tan poco variables resulta más viable proponer mecanismos de estabilización de corriente
para sostener la corriente. La distribución espacial de los círculos del meandro y las diferentes
posiciones de los centroides en el proceso dinámico se referencian en la Figura 11 y evidencian que
ese punto se sostendrá en periodo de evolución, y el ponteadero se considera viable alinearlo donde
anteriormente existía la obra de paso.

Página 29 de 156
2.2.6. Patrón morfológico actual

El reconocimiento del patrón hidrogeomorfológico de una corriente da cuentas del mecanismo de


implante, mantenimiento y perduración de las obras hidráulicas, para ello es fundamentar establecer
la edad, ciclo y momento de desarrollo de la corriente. Para (Universidad del Cauca, 2006), el ciclo de
un río se divide en tres etapas principales, comparables a la edad del hombre: juventud o zona de
montaña, madurez o valles intermedios y vejez o zonas planas. Cada una de estas etapas tiene sus
propias características, sin embargo, resulta difícil establecer límites precisos entre dichas etapas. El
ciclo de erosión brinda una descripción cualitativa de los valles de los ríos y de las áreas, pero cuando
se trata de asignarles valores cuantitativos surgen las dificultades, pero el concepto es un artificio de
gran valor, que permite describir con palabras el marco de desarrollo de un valle o de una región.

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.

Figura 13. Patrón de evolución del canal en la zona de influencia.

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.

Figura 14. Histórico de sismos en el área de influencia general de la microcuenca

Fuente: (Servicio Geológico Colombiano, 2022) 1 El círculo rojo indica el área de estudio.

• Factores del suelo

Aspectos como la permeabilidad, la textura, la estructura y la susceptibilidad a la erosión de los suelos


de las zonas de vertiente, son muy importantes en los procesos de cuantificación del aporte de
sedimentos y escorrentía superficial y subsuperficial que estas zonas hacen a los ríos. Por
consiguiente, deben ser tomados en cuenta en el análisis de la dinámica fluvial, puesto que influyen
no sólo en las condiciones hidrológicas sino también en la morfología del cauce, bien sea directa o
indirectamente.

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.

• Factores del uso y cobertura del suelo

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

Los cambios geomorfológicos en el río están asociados a procesos de degradación lateral,


degradación vertical y la agradación del cauce en el tiempo. También deben incluirse procesos que
tienen lugar a corto plazo como son las variaciones de las márgenes y el fondo del cauce, asociadas
usualmente con crecientes y las inundaciones por ella causadas.

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.

Tres clases de materiales se distinguen en un cauce natural considerando únicamente la resistencia


que ofrecen a ser transportados por una corriente: materiales no cohesivos o granulares, materiales
cohesivos y rocas.

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.

En particular, el río posee un complejo de socavación intensiva – progresiva, y finalmente estable, la


cual se encuentran a margen izquierda sobre el ponteadero. En este caso el suelo granular y
transportable se ha desgastado y transportado, realizando hoyas de profundización hasta llegar al
nivel de rodadura en la lodolita, lo que asemeja un coeficiente de Manning muy similar al concreto
pulido.

Página 35 de 156
Figura 15. Interacción geomorfológica del cauce de socavación.

Fuente: Elaboración propia, 2022

• 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.

El estudio hidráulico está orientado a la determinación de parámetros de interés a la hora de evaluar


los diferentes procesos fluviales, estando entre los más importantes las velocidades, los niveles de
flujo y el ancho de la superficie libre del flujo asociados a caudales para diferentes períodos de retorno.

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.

Procesos de degradación y agradación se acentúan durante crecientes, siendo evidente en el lugar


donde en la jerarquización geomorfológica indicó la presencia de un abanico antiguo. Si bien su
datación ya fue expuesta, en este proceso se debe considerar en el aspecto de diseño porque la
corriente actualmente se encuentra socavando sobre el talud ápice oriental, donde existen suelos muy
pobres de cohesión, cobertura compuesta por bosques de galería intervenidos. La zona anterior al
ponteadero presenta otro episodio de socavación lateral sobre talud compuesto por vegetación
secundaria en transición y pastos enmalezados.

2.2.7. Geometría hidráulica

El patronamiento de la corriente a nivel de planta – perfil permite reconocer variables de elementos


que delatan el probable comportamiento de la corriente a futuro, y cómo esta se servirá en relación
con la obra a implantar. Por esta razón Martín Vide (2004) menciona que si se le permiten todos sus
grados de libertad, un río es una corriente de agua que forma su cauce ajustando libremente la
anchura, el calado y la planta. El problema de deducir o predecir esta geometría del río ha ocupado a
muchos investigadores en geomorfología. Se cuenta con observaciones relativamente sencillas de las
características geométricas de los ríos que sugieren relaciones empíricas entre ellas. Dando un paso
más, la pregunta de porque adopta el río tales formas permanece todavía sin respuesta completa; este
es el caso del porqué de los meandros. Dada la situación, clasificar los ríos constituye un asunto
improrrogable en el diseño de obras de paso, puesto que de fondo el transporte de sedimentos
convierte esta fuente en un sistema de Flujo a dos fases, el cual es el movimiento simultáneo del agua
(fase líquida) y de los sólidos constituyentes del lecho (fase sólida) (Rocha Felices, 1998).

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í,

iQ0.44 < 0.0116 cauce meandriforme


Ec. 1
iQ0.44 > 0.0116 cauce trenzado

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.

Finalmente, λ es proporcional a Q95/D03

Figura 16. Parámetros de una planta sinuosa

Fuente: Martín Vide (2003)

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

Fuente: Elaboración propia, 2022

Otro parámetro geométrico de un río meandriforme es la sinuosidad ṣ, cociente de la longitud ᶴ a lo


largo del valle y la longitud de onda λ (como en la figura anterior). La geometría en planta de un
meandro regular y simple parece describirse bien mediante la ecuación: Ꝋ = Ꝋ0 sen (2π ṣ/ᶴ), donde
Ꝋ es el ángulo del eje con la dirección del valle y ṣ es la coordenada arco. Con estos meandros se
puede obtener meandros en cualquier grado de desarrollo, incluso al límite de su estrangulamiento
(Martín Vide, 2003), en ese sentido, la corriente se ha obtenido un grado de sinuosidad de 1.18
indicando que se trata de una corriente meándrica, en proceso de divagación conforme a las
definiciones de Martín (2017).

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.

La referencia asignataria corresponde a la siguiente descripción:

• Recto: cauce único con índice de sinuosidad inferior a 1.1.


• Sinuoso: cauce único con índice de sinuosidad entre 1.1 y 1.3.
• Meandriforme: cauce único con índice de sinuosidad superior a 1.3.
• Divagante: con sinuosidades o meandros, pero abundante material grueso en barras y alguna
subdivisión del cauce.
• Trenzado: de gravas, con división en subcauces móviles que se entrecruzan y que cambian
frecuentemente de posición.
• Anastomosado: de alta montaña en nuestra latitud, con sinuosidades e islas fijas.
• Rambla: con un cauce seco de caracteres próximos al trenzado pero característico de zonas
áridas.

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.

Figura 18. Descarga en el banco, en pies cúbicos por segundo

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.

Figura 19. Clasificación de los ríos aluviales

Fuente: Schumm (1960, 1963a, 1963b, 1985)

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.

2.2.8. Unidades de canal

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:

Figura 20. Formas fundamentales de piscina/rifle, reflejando la topografía de la cama

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.

• Adecuación piscícola (Ap): Son intervenciones que se realizan en el medio de áreas


hidrográficas donde para el abastecimiento del agua se aprovecha el medio superficial o
subterráneo con modificaciones de ambiente.

• 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.).

• Otras causas: si los cambios observados no responden a ninguna causa fácilmente


identificable, ya sea de origen natural o antrópico.

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.

2.2.10. Identificación de tramos modificados por acciones directas en el cauce

En el departamento del Putumayo es común que a menudo se desarrollen obras de lateralización,


descolmatación u modificaciones hidráulicas de los ríos cuando estos amenazan la infraestructura
existente sea obras públicas o viviendas particulares, por tal razón, una vez definido tipos morfológicos
en planta y caracterizados inicialmente, tanto para el trazado en situación de no alteración y en el
estado actual del curso fluvial, como se ha estudiado anteriormente, cabe la posibilidad de que existan
en la masa de agua tramos con el cauce claramente modificado (Ministerio para la transición
Ecológica, 2019), debiéndose identificar estos tramos afectados por cualesquiera de las acciones que
a continuación se detallan:

• 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.

• Acortado (Ac): cortas artificiales de meandros o actuaciones de reducción de curvatura. Se


incluyen todas las acciones realizadas sobre el cauce original que suponen una reducción en
la longitud del mismo, siempre que no supongan un desvío o cambio del 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.

2.2.11. Síntomas de dinámica vertical acelerada

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

Fuente: (Ministerio para la transición Ecológica, 2019)

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).

Teniendo en referencia la Figura 22, la estructura se encontraría en un área de incisión en lecho, y


perdida de espacio de rivera. En el primer caso no es crítico por lo ya citado (presencia de la formación
Orito), pero en el soporte derecho se deberá de proteger los suelos degradables con algún tipo de
obra continua de cobertura del talud.

Figura 22. Síntomas de dinámica vertical acelerada

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.

2.2.12.1. Tipo de sustrato

El sustrato del lecho se clasificará según los siguientes tipos:

• 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.

2.2.12.2. Tamaño dominante del sedimento

En la caracterización de este protocolo no se ha considerado necesario tomar muestras ni medir con


precisión los sedimentos, simplemente se propone la realización una rápida inspección visual tratando
de señalar cuál es el tamaño dominante en los sedimentos superficiales del lecho y de las posibles
barras o depósitos. Es importante señalar que el tamaño se estima sobre el mayor eje transversal del
elemento (el que cabría o no por un tamiz).

En el caso de que se detectase en el lecho la presencia de lodos procedentes de problemas presentes


o pasados de calidad de las aguas, también se procederá a señalar e identificar en esta
caracterización.

De acuerdo con lo anterior, las categorías de los sedimentos pueden ser las siguientes:

Tabla 8. Estructura y sustrato del lecho: tamaño dominante del sedimento


SIN SEDIMENTO SI EL CAUCE ES EN ROCA NO SE OBSERVARÁN
SEDIMENTOS SUPERFICIALES
ROCOSO Bloques más de 25,6 cm
GRUESO Cantos de 6,4 cm a 25,6 cm

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.

2.2.12.3. Clasificación de sedimentos

La clasificación de los sedimentos en el subtramo de muestreo proporciona mucha información acerca


del funcionamiento del cauce y su dinámica fluvial, es decir, de si el régimen de caudales existente en
el río es o no capaz de movilizar un determinado rango de sedimentos y/o los sedimentos existentes
o barras responden o no a unas condiciones morfológicas adecuadas.

Para caracterizar los sedimentos se comprobarán tres indicadores de fácil visualización:

• Imbricación: se da cuando cada partícula o clasto se encuentra bien colocada, inclinada y


sobre otra como las tejas de un tejado, en la dirección de la corriente. Indica una buena
clasificación de sedimentos y, por lo tanto, un buen trabajo fluvial. Si la imbricación es tenue
o inexistente puede ser síntoma de escasez de caudales generadores o de alteraciones
antrópicas.

• Acorazamiento: normalmente el material sedimentario en superficie (coraza) es de mayor


tamaño que el que hay inmediatamente debajo, debido a los propios procesos de movilización
condicionados por el caudal y las condiciones de flujo. Si esa coraza superior es muy marcada,
con partículas de tamaño mucho mayor (más de 3 o 4 veces mayor) que lo que está
depositado debajo, el acorazamiento se considera excesivo. También puede ser que el
acorazamiento sea débil o inexistente. La situación que marca una mayor salud en los
sedimentos fluviales es un acorazamiento moderado, de entre 1,5 a 3 veces mayores los
clastos de la coraza que el material inferior.

• 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.

Si el material sedimentario dominante es arena o finos, por su propia definición no existe el


acorazamiento y la imbricación se sustituye por la posible presencia de formas del lecho sedimentarias
que muestren una adecuada construcción de los depósitos, tales como los ripples o dunas. En estos

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:

- Clasificación efectiva: los sedimentos se encuentran imbricados, hay un acorazamiento moderado y


las barras se encuentran bien formadas.
- Clasificación limitada: falla alguno de los indicadores, no se observa imbricación clara o el
acorazamiento es muy alto o muy bajo o las barras son incipientes. Si el sedimento es arena o fino no
habrá una buena construcción de depósitos.
- Clasificación nula: fallan todos los aspectos. No se registra imbricación, no hay acorazamiento o este
es excesivo y no hay depósitos sedimentarios claros.

Tabla 9. Clasificación de sedimentos

a) Imbricación en el lecho de gravas y b) Imbricación (ripples) en el lecho arenoso


cantos

Página 50 de 156
c) Acorazamiento excesivo d) Acorazamiento débil

e) Barras bien formadas

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.

2.2.12.4. Tipo de estructura vertical

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.

Se caracterizará el tipo de estructura longitudinal eligiendo entre los siguientes tipos:

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.

2.2.12.5. Movilidad de sedimentos

La movilidad de los sedimentos es un aspecto básico en la morfología fluvial. Para definirla se


comprobarán de forma rápida en el subtramo de muestreo cuatro indicadores:

• 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.

• El grado de colonización vegetal de las barras o depósitos sedimentarios. La vegetación dentro


del cauce suele indicar escasez de crecidas y caudales. La colonización vegetal reduce o
impide la movilización de sedimentos, por lo que puede generar déficits aguas abajo y llegar
a provocar procesos de incisión.

• 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 limitada (Ml): el sedimento muestra algunos síntomas de compactación o asiste a


una colonización vegetal moderada, o está cubierto por materiales finos o arcillosos que
impiden la movilidad.

• 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

Dentro de la hipsometría de la cuenca se puede identificar claramente las unidades de producción,


transporte y depósito, puesto que se trata de una cuenca larga de piedemonte donde los nacimientos
corresponden a fuentes permanente lóticas, se sugiere que la producción se encuentra sobre las
cotas 3551 a 2340 msnm, transporte hasta los 1280 msnm y el depósito hasta la desembocadura en
el río Putumayo sobre los 240 msnm.

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.

3.2. PARÁMETROS MORFOMÉTRICOS DE LA CUENCA

Las características morfométricas corresponden a la aplicación de procedimientos que a través del


estudio de la morfología y geomorfología caracterizan de manera cuantitativa los rasgos propios de
las cuencas hidrográficas en valores numéricos los cuales permiten comparar en forma exacta una
parte de la superficie terrestre con otra.

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.

Las características morfométricas corresponden a la aplicación de procedimientos que a través del


estudio de la morfología y geomorfología caracterizan de manera cuantitativa los rasgos propios de
las cuencas hidrográficas en valores numéricos los cuales permiten comparar en forma exacta una
parte de la superficie terrestre con otra.

Para la caracterización morfométrica se verificaron y se dibujaron las divisorias de aguas y se


calcularon las características morfométricas más relevantes, sobre la cartografía básica establecida
para el proyecto. Dichos parámetros fueron revisados por el SIG, con el fin de evitar inconsistencias
en las medidas. A continuación, se presenta el cálculo de los parámetros morfométricos del río San
Juan, a partir del modelo de elevación digital disponible.

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.

3.2.1. Parámetros de forma

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.

En términos básicos, la forma de la Cuenca Hidrográfica es importante pues se relaciona con el


Tiempo de Concentración (Tc), el cual es el tiempo necesario, desde el inicio de la precipitación, para
que toda la cuenca contribuya al cauce principal en estudio, es decir, el tiempo que toma el agua
precipitada en los límites más extremos de la cuenca para llegar al punto de salida de la misma, en
(Escuela superior de administración Pública - ESAP, 2012).

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

Fuente: Elaboración propia, 2022

3.2.1.1. Área de drenaje de la cuenca – Ac

Es la proyección horizontal del área de drenaje de un sistema de escorrentía dirigido, directa o


indirectamente, a un mismo cauce natural. El sitio que recoge toda la escorrentía que se produce en
una cuenca hidrográfica se denomina punto de concentración.

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).

La delimitación de una cuenca hidrográfica se realiza a partir de restituciones cartográficas y


fotogramétricas y para ello se deben seguir las reglas básicas que se mencionan a continuación:

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.

3.2.1.2. Perímetro de la cuenca

El perímetro de la cuenca o la longitud de la línea divisoria de la hoya es un parámetro importante,


pues en conexión con el área nos puede decir algo sobre la forma de la cuenca. Usualmente este
parámetro físico es simbolizado por la mayúscula P. El caso del río San Juan posee un perímetro de
144 km.

3.2.1.3. Ancho promedio – Ap

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
𝐴𝑝 =
𝐿𝑏

Donde A corresponde al área de la cuenca expresada en Km2, y Lb a la Longitud Axial calculada en


Km.

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.

3.2.1.4. Relación de elongación – Re

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.

3.2.1.5. Factor de forma – Kf

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.

3.2.1.6. Coeficiente de compacidad – Kc

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.

3.2.2. Parámetros hidrográficos

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.

3.2.2.1. Relación de bifurcación – Rb

Es la relación entre el número de corrientes de cualquier orden u (Nu) y el número de corrientes en


el siguiente orden superior u+1 (Nu+1)5:

𝑁𝑢 Ec. 7
𝑅𝑏 =
𝑁𝑢+1

La relación de bifurcación de la cuenca de análisis es de 4.2.

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
𝑃𝑚𝑐 =
𝐿

La pendiente media de la cuenca del proyecto es 0.054 m/m.

3.2.2.3. Patrón de drenaje

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).

Figura 27. Patrones de drenaje típicos

Fuente: Referencia al pie6

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.

# 𝐶𝑎𝑢𝑐𝑒𝑠 𝑜𝑟𝑑𝑒𝑛 1 Ec. 9


𝐶1 =
𝐴𝑟𝑒𝑎 𝑐𝑢𝑒𝑛𝑐𝑎 (𝐾𝑚2 )

La descripción de las características morfométricas, fisiográficas, hidrológicas e hidráulicas entre


otros, permite un mejor entendimiento del comportamiento del flujo de agua en la cuenca, (Bateman,
2007). La microcuenca de análisis posee una torrencialidad de 0.9, entendiéndose como una cuenca
con baja capacidad de generar torrentes.

3.2.2.5. Densidad de drenaje – Dd

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
𝐷𝑑 =
𝐴

La densidad de drenaje de la microcuenca de análisis es de 1.61, es decir, se trata de una cuenca


bien drenada. No obstante, de los cálculos individuales de la cuenca en general, se ha generado los
parámetros calculados para las 23 subcuencas del área del proyecto, los cuales se encuentran en el
anexo del proyecto (\\2022\10_Puente_Vides\02_HIDROLOGIA\01_MORFOMETRIA)

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)

3.2.4. Coberturas del suelo

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

Fuente: Elaboración propia, 2022

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.

4.1. MECANISMOS DE CÁLCULO

Con la información de la precipitación máxima de 24 horas se preparó una curva de Intensidad-


Frecuencia-Duración (IDF), de acuerdo a la siguiente metodología:

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

Tabla 12. Periodos de retorno para el diseño de obras de drenaje vial


OBRAS DE DRENAJE VIAL TR
DRENAJE SUPERFICIAL 2
CUNETAS 5
ZANJAS DE CORONACIÓN 10
ESTRUCTURAS DE CAÍDA 10
ALCANTARILLAS DE 0.90 M DE DIÁMETRO 10
ALCANTARILLAS MAYORES A 0.90 M DE DIÁMETRO 20
PUENTES MENORES (LUZ MENOR A 10 M) 25
PUENTES DE LUZ MAYOR O IGUAL A 10 M Y MENOR A 50
50 M
PUENTES DE LUZ MAYOR O IGUAL A 50 M 100
DRENAJE SUBSUPERFICIAL 2
Fuente: INVIAS & Ministerio de Transporte, 2009

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).

Valores del coeficiente de escorrentía en áreas urbanas

Tabla 13. Valores del coeficiente de escorrentía en áreas urbanas


TIPO DE ÁREA DE DRENAJE COEFICIENTE DE
ESCORRENTÍA C
PRADOS
SUELOS ARENOSOS, PLANOS, 2% 0.05 ˍ 0.10
SUELOS ARENOSOS, PROMEDIO, 2 -7 % 0.15 ˍ 0.20
SUELOS PESADOS (ARCILLOSOS), PLANOS, 2% 0.13 ˍ 0.17
SUELOS PESADOS (ARCILLOSOS), PROMEDIO, 2 0.18 ˍ 0.22
-7 %
SUELOS PESADOS (ARCILLOSOS), PENDIENTES, 0.25 ˍ 0.35
7%
DISTRITOS COMERCIALES
ÁREAS DE CENTRO DE CIUDAD 0.70 ˍ 0.95
ÁREAS VECINAS 0.50 ˍ 0.70
RESIDENCIAL
CASAS INDIVIDUALES SEPARADAS 0.30 ˍ 0.50
CASAS MULTIFAMILIARES SEPARADAS 0.40 ˍ 0.60
CASAS MULTIFAMILIARES UNIDAS 0.60 ˍ 0.75

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

Valores del coeficiente de escorrentía en áreas rurales:

Tabla 14. Valores del coeficiente de escorrentía en áreas rurales


VEGETACIÓN Y TEXTURA DEL SUELO
TOPOGRAFÍA Franco arenoso Franco limo arcilloso Arcilloso
BOSQUES
PLANO 0.10 0.30 0.40
ONDULADO 0.25 0.35 0.50
MONTAÑOSO 0.30 0.50 0.60
PASTOS
PLANO 0.10 0.30 0.40
ONDULADO 0.16 0.36 0.55
MONTAÑOSO 0.22 0.42 0.60
TIERRAS CULTIVADAS
PLANO 0.30 0.50 0.60
ONDULADO 0.40 0.60 0.70
MONTAÑOSO 0.52 0.72 0.82
NOTA: PLANO (PENDIENTE 0 - 5%); ONDULADO (PENDIENTE 5 - 10%); MONTAÑOSO (PENDIENTE
10 - 30%). PARA VALORES MAYORES AL 30 %, A FALTA DE DATOS, UTILIZAR LOS VALORES
PARA PENDIENTES ENTRE EL 10 Y EL 30 %.
Fuente: Ponce, 1989, citado por INVIAS & Ministerio de Transporte, 2009

4.2. ECUACIÓN DE INTENSIDAD

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).

Para determinar la intensidad de la lluvia se puede aplicar la Metodología Simplificada de cálculo de


las curvas intensidad – duración – frecuencia indicada en el Manual de Drenaje de Carreteras del

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

i: Intensidad de precipitación, en milímetros por hora (mm/h).


T: Periodo de retorno, en años.
M: Precipitación máxima promedio anual en 24 h a nivel multianual.
t: Duración de la lluvia, en minutos (min).

a, b, c, d: Parámetros de ajuste de la regresión. Estos parámetros dependen de la región del país


donde está ubicada la estación meteorológica. Las regiones y los valores de los parámetros se
muestran en la Tabla a continuación.

Tabla 15. Parámetros según regiones del país


REGIÓN A B C D
ANDINA 0.94 0.18 0.66 0.83
CARIBE 24.85 0.22 0.5 0.1
PACÍFICO 13.92 0.19 0.58 0.2
ORINOQUÍA 5.53 0.17 0.63 0.42
Fuente: Armenta Porras (2013)

4.3. CÁLCULOS HIDROLÓGICOS

De acuerdo con lo referenciado en el método de estimación de áreas aferentes de las estaciones


meteorológicas para la microcuenca, se dispuso a realizar los cálculos que implican: estimación de
lluvias máximas mensuales multianuales – con su serie histórica de picos, lluvias máximas multianuales
en 24 horas – con serie de picos, distribución mensual multianual de meses húmedos y temporada
seca.

Así las cosas, la estación aferente para la zona de estudio es la siguiente:

Tabla 16. Estaciones aferentes a la zona de estudio


MUNICIPIO CORRIENTE ESTACIÓN CÓDIGO TIPO
METEOROLÓGICA ESTACIÓN
PUERTO Putumayo Puerto Caicedo 47010110 Pluviométrica
CAICEDO
PASTO Nariño Santa Isabel 47010230 Pluviométrica
SANTIAGO Putumayo Balsayaco 47010020 Pluviográfica
SAN Putumayo Chungacaspi 47010050 Pluviométrica
FRANCISCO

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

Fuente: Elaboración propia, 2022

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:

Tabla 17. Porcentajes de estaciones aferentes a cada subcuenca


SUBCUENCA ESTACIÓN AFERENTE % DE AFERENCIA
C1 BALSAYACO [47010020] 45.5%
C1 CHUNGACASPI [47010050] 37.2%
C1 EL PEPINO - AUT [44015070] 1.0%
C1 PUERTO CAICEDO [47010110] 0.4%
C1 TORRE TV SAN FCO [47010180] 3.0%
C1 VILLAGARZON [44015010] 12.8%
C10 VILLAGARZON [44015010] 100.0%
C11 BALSAYACO [47010020] 74.1%
C11 SANTA ISABEL [47010230] 25.9%
C12 BALSAYACO [47010020] 1.8%
C12 PUERTO CAICEDO [47010110] 48.2%
C12 SANTA ISABEL [47010230] 50.0%
C13 PUERTO CAICEDO [47010110] 97.9%
C13 VILLAGARZON [44015010] 2.1%
C14 PUERTO CAICEDO [47010110] 26.2%
C14 VILLAGARZON [44015010] 73.8%
C15 PUERTO CAICEDO [47010110] 83.2%
C15 VILLAGARZON [44015010] 16.8%
C16 PUERTO CAICEDO [47010110] 96.8%
C16 VILLAGARZON [44015010] 3.2%
C17 PUERTO CAICEDO [47010110] 100.0%
C18 PUERTO CAICEDO [47010110] 100.0%
C19 PUERTO CAICEDO [47010110] 38.3%
C19 VILLAGARZON [44015010] 61.7%

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

Fuente: Elaboración propia, 2022

4.3.1. Análisis meteorológico

4.3.1.1. Precipitación

• Análisis de las series de tiempo

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

Fuente: Elaboración propia, 2022


Fuente: Elaboración propia, 2022
1/05/2010
1/05/2013
1/07/2011
1/09/2012 1/04/2014
1/11/2013 1/03/2015
1/01/2015 1/02/2016
1/03/2016 1/01/2017
1/05/2017 1/12/2017
1/07/2018 1/11/2018
1/09/2019 1/10/2019
1/11/2020 1/09/2020
1/01/2022 1/08/2021

0
0

20000
40000
60000
80000

-20000
100000
120000
140000
160000
180000
10000
20000
30000
40000
50000
60000

Precipitación acumulada multianual (mm) Precipitación acumulada multianual (mm)


Figura 30. Serie temporal de precipitación y curva de masa simple de la estación Santa Isabel

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

Fuente: Elaboración propia, 2022


Fuente: Elaboración propia, 2022
1/04/2008 1/08/2010
1/09/2009 1/10/2011
1/02/2011 1/12/2012
1/07/2012 1/02/2014
1/12/2013 1/04/2015
1/05/2015 1/06/2016
1/10/2016 1/08/2017
1/03/2018 1/10/2018
1/08/2019 1/12/2019
1/01/2021 1/02/2021
1/04/2022

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

Precipitación acumulada multianual (mm) Precipitación acumulada multianual (mm)


Figura 33. Serie temporal de precipitación y curva de masa simple de la estación Chungacaspi
Precipitación total mensual (mm) Precipitación total mensual (mm)

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

Fuente: Elaboración propia, 2022


Fuente: Elaboración propia, 2022
1/04/2008
1/09/2007
1/09/2009 1/03/2009
1/02/2011 1/09/2010
1/07/2012 1/03/2012
1/12/2013 1/09/2013
1/05/2015 1/03/2015
1/10/2016 1/09/2016
1/03/2018 1/03/2018
1/08/2019 1/09/2019
1/01/2021 1/03/2021

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

Precipitación acumulada multianual (mm) Precipitación acumulada multianual (mm)


Figura 34. Serie temporal de precipitación y curva de masa simple de la estación Villagarzón
Figura 36. Serie temporal de precipitación y curva de masa simple de la estación El Pepino
160 40000

140 35000

Precipitación acumulada multianual (mm)


120 30000
Precipitación total mensual (mm)

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

• Valores totales mensuales

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

Fuente: Elaboración propia, 2022

Figura 38. Distribución precipitación mensual media – Puerto Caicedo


400.0

350.0

300.0
Precipitación (mm)

250.0

200.0

150.0

100.0

50.0

0.0

Fuente: Elaboración propia, 2022

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

Fuente: Elaboración propia, 2022

Figura 40. Distribución precipitación mensual media – Chungacaspi


500.0
450.0
400.0
350.0
Precipitación (mm)

300.0
250.0
200.0
150.0
100.0
50.0
0.0

Fuente: Elaboración propia, 2022

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

Fuente: Elaboración propia, 2022

Figura 42. Distribución precipitación mensual media – Balsayaco


250.0

200.0
Precipitación (mm)

150.0

100.0

50.0

0.0

Fuente: Elaboración propia, 2022

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

Fuente: Elaboración propia, 2022

Tabla 18. Distribución temporal de la precipitación – Estación Santa Isabel


DATO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC MEDIA
PROMEDIO 60.6 64.8 81.6 92.2 107.9 104.0 118.3 99.5 70.1 64.8 68.5 68.0 1058.16
MÁXIMA 240.4 290.6 359.8 281.8 407.7 324.0 413.0 346.2 245.8 367.0 189.6 269.3 413.00
MÍNIMA 11.4 24.1 61 75.1 104.5 83 86.2 58.1 29.7 5.6 50 48 5.6
Fuente: Elaboración propia, 2022

Tabla 19. Distribución temporal de la precipitación – Estación Puerto Caicedo


DATO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC MEDIA
PROMEDIO 182.4 203.5 259.9 331.0 338.2 282.7 233.3 200.7 190.6 232.1 276.9 220.9 3060.20
MÁXIMA 520.1 842.0 677.0 779.5 719.4 605.0 485.0 714.8 450.8 497.5 775.1 620.2 842.00
MÍNIMA 5.3 23.5 17.1 48.3 39.7 50.5 18.3 36.3 24.6 28 46.6 28.6 5.3
Fuente: Elaboración propia, 2022

Tabla 20. Distribución temporal de la precipitación – Estación Torre TV


DATO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC MEDIA
PROMEDIO 168.8 185.5 239.0 303.4 409.6 528.2 600.2 410.1 257.7 201.8 173.6 174.0 3807.24
MÁXIMA 425.9 475.0 664.6 664.3 933.2 1046.1 1159.8 826.8 506.7 392.7 363.4 403.2 1159.80
MÍNIMA 75.2 44.7 136.4 160 276 298.2 375.3 178.2 152.6 155.3 87.9 65.1 44.7

Página 81 de 156
Fuente: Elaboración propia, 2022

Tabla 21. Distribución temporal de la precipitación – Estación Chungacaspi


DATO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC MEDIA
PROMEDIO 107.8 115.1 148.6 181.8 258.3 375.1 445.9 289.7 200.8 132.2 111.9 117.2 2591.00
MÁXIMA 242.8 302.2 438.6 435.0 570.0 814.2 773.3 600.2 389.5 377.2 209.9 273.7 814.20
MÍNIMA 16.3 11.5 25.2 15.3 30.3 27.5 34.5 20.2 10.6 9.9 15.9 17.6 9.9
Fuente: Elaboración propia, 2022

Tabla 22. Distribución temporal de la precipitación – Estación Villagarzón


DATO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC MEDIA
PROMEDIO 308.0 272.3 339.1 429.3 516.6 512.9 424.2 292.1 310.8 350.6 343.9 327.8 4512.52
MÁXIMA 801.0 709.3 574.0 795.0 854.9 878.3 720.9 665.4 681.5 588.2 881.4 600.3 881.40
MÍNIMA 12.4 62.3 185.8 125.2 178.4 187.9 212.8 95.9 0.4 4.6 78.6 52.3 0.4
Fuente: Elaboración propia, 2022

Tabla 23. Distribución temporal de la precipitación – Estación Balsayaco


DATO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC MEDIA
PROMEDIO 109.5 106.7 134.7 150.5 195.4 237.5 231.6 182.8 143.5 111.1 105.2 106.1 1804.90
MÁXIMA 213.5 294.3 288.8 315.6 415.7 365.6 393.4 286.2 231.7 195.9 173.3 245.4 415.70
MÍNIMA 11 8.1 11 14.2 22 19.8 25.8 18.8 11.4 12.8 15.5 11.1 8.1
Fuente: Elaboración propia, 2022

Tabla 24. Distribución temporal de la precipitación – Estación Pepino


DATO ENE FEB MAR ABR MAY JUN JUL AGO SEP OCT NOV DIC MEDIA
PROMEDIO 221.0 238.4 349.9 352.9 314.4 390.9 342.2 258.6 235.7 185.6 230.8 280.7 3400.92
MÁXIMA 401.5 430.0 567.0 715.1 560.1 884.7 569.8 585.0 470.7 370.3 732.0 672.0 884.70
MÍNIMA 40 75 107 307 223 341 55 153 106 93 4 130 4
Fuente: Elaboración propia, 2022

• Precipitación máxima en 24 horas

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

FRECUENCIAS DE PRECIPITACIÓN MÁXIMA ANUAL A NIVEL MULTIANUAL

300

Serie

Normal

Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)

Normal
Gumbel

200 Pearson

Log
Pearson
EV3

Abscisa

100

PROBABILIDAD P(X < Xi)

Fuente: Elaboración propia, 2022

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

FRECUENCIAS DE PRECIPITACIÓN MÁXIMA ANUAL A NIVEL MULTIANUAL

300

Serie

Normal

Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)

Normal
Gumbel

200 Pearson

Log
Pearson
EV3

Abscisa

100

PROBABILIDAD P(X < Xi)

Fuente: Elaboración propia, 2022

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

FRECUENCIAS DE PRECIPITACIÓN MÁXIMA ANUAL A NIVEL MULTIANUAL

300

Serie

Normal

Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)

Normal
Gumbel

200 Pearson

Log
Pearson
EV3

Abscisa

100

PROBABILIDAD P(X < Xi)

Fuente: Elaboración propia, 2022

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

Fuente: Elaboración propia, 2022

Figura 47. Frecuencia de precipitación máxima en 24 horas a nivel multianual estación Chungacaspi

FRECUENCIAS DE PRECIPITACIÓN MÁXIMA ANUAL A NIVEL MULTIANUAL

300

Serie

Normal

Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)

Normal
Gumbel

200 Pearson

Log
Pearson
EV3

Abscisa

100

PROBABILIDAD P(X < Xi)

Fuente: Elaboración propia, 2022

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

FRECUENCIAS DE PRECIPITACIÓN MÁXIMA ANUAL A NIVEL MULTIANUAL

300

Serie

Normal

Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)

Normal
Gumbel

200 Pearson

Log
Pearson
EV3

Abscisa

100

PROBABILIDAD P(X < Xi)

Fuente: Elaboración propia, 2022

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

FRECUENCIAS DE PRECIPITACIÓN MÁXIMA ANUAL A NIVEL MULTIANUAL

300

Serie

Normal

Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)

Normal
Gumbel

200 Pearson

Log
Pearson
EV3

Abscisa

100

PROBABILIDAD P(X < Xi)

Fuente: Elaboración propia, 2022

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

FRECUENCIAS DE PRECIPITACIÓN MÁXIMA ANUAL A NIVEL MULTIANUAL

300

Serie

Normal

Log
PRECIPITACIÓN MÁXIMA EN 24 HORAS (mm)

Normal
Gumbel

200 Pearson

Log
Pearson
EV3

Abscisa

100

PROBABILIDAD P(X < Xi)

Fuente: Elaboración propia, 2022

Página 89 de 156
4.4. CURVAS DE INTENSIDAD – DURACIÓN Y FRECUENCIA (IDF)

Para la determinación de estas curvas, inicialmente se consultó en el IDEAM la existencia de curvas


IDF de las estaciones analizadas, encontrando que no hay disponibilidad de esta información para la
integralidad, solo existe para la estación Villagarzón, pero los cálculos estimados son de una década
atrás. Para la estimación de curvas IDF, se empleó la formulación de “Curvas sintéticas de Intensidad
– Duración – Frecuencia para Colombia” de (Vargas & Diaz, 1998). Universidad de los Andes – 1998.
La ecuación definida en dicho estudio como Número 8 fue la empleada para la elaboración de las
curvas IDF. La zona en estudio se encuentra dentro de la zona Andina, por consiguiente, los datos de
las estaciones se analizan con base en los coeficientes definidos para la ecuación 11 en la región 1
(Andina) dentro de la referencia bibliográfica. No obstante, como se trata de un proyecto de
infraestructura donde la norma a seguir corresponde a las definidas por el INVIAS, se ha tomado las
curvas generadas con la implementación de la Ecuación 11 del presente estudio, consiguiendo las
siguientes curvas IDF:

Figura 51. Curvas IDF estación Santa Isabel


180
2 años 5 años 10 años 25 años
160
50 años 75 años 100 años
140

120
Intensidad (mm/hr)

100

80

60

40

20

0
0 20 40 60 80 100
Tiempo de duración (min)

Fuente Elaboración propia, 2022

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)

Fuente: Elaboración propia, 2022

Figura 53. Curvas IDF estación Torre TV


350
2 años 5 años 10 años 25 años
300 50 años 75 años 100 años

250
Intensidad (mm/hr)

200

150

100

50

0
0 20 40 60 80 100
Tiempo de duración (min)

Fuente: Elaboración propia, 2022

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

50 años 75 años 100 años


200
Intensidad (mm/hr)

150

100

50

0
0 20 40 60 80 100
Tiempo de duración (min)

Fuente: Elaboración propia, 2022

Figura 55. Curvas IDF estación Villagarzón


450
2 años 5 años 10 años 25 años
400
50 años 75 años 100 años
350

300
Intensidad (mm/hr)

250

200

150

100

50

0
0 20 40 60 80 100
Tiempo de duración (min)

Fuente: Elaboración propia, 2022

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)

Fuente: Elaboración propia, 2022

Figura 57. Curvas IDF estación Pepino


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)

Fuente: Elaboración propia, 2022

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).

La distribución espacial de la lluvia tiene gran incidencia en el comportamiento de caudales de una


cuenca hidrográfica, teniendo en cuenta la variación de características físicas de la cuenca como el
relieve, tipo de suelo, tipo y densidad de cobertura vegetal, uso del suelo, además de la intervención
humana (actividades), destacando la urbanización, que conlleva a la impermeabilización del suelo en
gran porcentaje del área, por lo que el comportamiento hidrológico de las corrientes que componen
el sistema hídrico, está claramente determinado, tanto espacial como temporalmente, por la
ocurrencia de la precipitación a lo largo de su territorio.

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.

Si bien a continuación, se describe la metodología general, al interior del texto, se describe o


desarrollan los elementos de la misma, incluyendo componentes nuevos que fueron necesario
desarrollar para el logro de uno de los objetivos del estudio.

4.5.1. Determinación del área hidrológica del estudio

La delimitación de la cuenca hidrográfica, se realizó a partir de la base cartográfica en formato digital


del Instituto Geográfico Agustín Codazzi escala 1:25.000 o una mayor escala o con información de un
modelo de elevación digital MDT obtenido de imágenes ALOS-PALSAR, y se delimita usando la
aplicación ArcHydro de ArcGIS.

4.5.2. Estimación del índice de capacidad de infiltración y escurrimiento del suelo

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.

4.5.3. Determinación del programa de cómputo hidrológico

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.

4.5.4. Análisis estadístico de la información de la estación seleccionada

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.

4.5.5. Cálculo de caudales

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.

Este método se desarrolló empíricamente por el departamento estadounidense de conservación de


suelos, en el año 1972, con el fin de estimar las abstracciones (las cuales se refieren principalmente

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.

Figura 58. Variables en el método de abstracciones del SCS.

Fuente: Chow, Maidment, & Mays (1994)

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

Fuente: Chow et al. (1994)

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.

Primero se clasifica hidrológicamente el suelo y se determina su potencial de escurrimiento. Según la


tasa de infiltración, la distribución de finos y gruesos y condiciones de drenaje, se puede establecer si
el suelo tiene un potencial bajo, moderadamente bajo, moderadamente alto y alto.

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

Fuente: Chow et al. (1994)

4.5.5.1. Tipo de suelo hidrológico

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.

Página 100 de 156


4.5.5.2. Uso y tratamiento del suelo

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

Página 101 de 156


del suelo. Dentro de los cultivos en hileras rectas se incluyen los cultivos que siguen la mayor
pendiente y los orientados según las curvas de nivel en hileras rectas (a media ladera). En las áreas
con 1 o 2% de pendiente, los cultivos transversales en hileras rectas son casi equivalentes a los
cultivos siguiendo las curvas de nivel. En lo referente a cultivos por líneas de nivel, se considera que
los surcos que se usan para plantar los granos pequeños y legumbres generalmente son muy
pequeños y tienden a desaparecer debido al efecto del clima. Los surcos y bancales en curvas de
nivel como se usan en los cultivos por hileras son generalmente grandes; se pueden hacer al sembrar
y después se reducen en tamaño por los cultivos; las reducciones que sufre el escurrimiento debido
a los cultivos a lo largo de las curvas de nivel son mayores al disminuir la pendiente del terreno; las
terrazas pueden tener cierta pendiente, pueden estar a nivel con los extremos abiertos, o a nivel con
los extremos cerrados; los efectos de las terrazas de pendiente con los extremos abiertos se
consideran en las evaluaciones y se incluyen los efectos de los surcos a nivel y de las salidas cubiertas
de pasto. Las terrazas a nivel con extremos cerrados deben considerarse como surcos a nivel.

4.5.5.3. Combinaciones hidrológicas de suelo y vegetación

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.

4.5.5.4. Condiciones de humedad antecedente a la lluvia de diseño

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.

Página 102 de 156


Tabla 32. Condiciones de humedad antecedente a la lluvia de diseño
PERÍODO VEGETATIVO LLUVIA ANTECEDENTE
PERÍODO INACTIVO DE LAS PLANTAS Lluvia acumulada de los 5 días previos al
aguacero en consideración
I 00,0 – 12,7 mm
II 12,7 – 28,0 mm
III más de 28,0 mm
PERÍODO DE CRECIMIENTO DE LAS
PLANTAS
I 00,0 – 35,5 mm
II 35,5 – 53,3 mm
III más de 53,3 mm
Fuente: Chow et al. (1994)

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.

4.5.5.5. Determinación de la lluvia efectiva

Este parámetro se calcula aplicando la relación:

(𝑃𝑒 )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).

Página 103 de 156


Tabla 33. Selección del valor CN” considerando diferentes combinaciones hidrológicas de suelo y vegetación
COBERTURA Y USO DEL TRATAMIENTO CONDICIÓN GRUPOS HIDROLÓGICOS DE
SUELO O PRÁCTICA HIDROLÓGICA SUELOS
A B C D
NÚMERO DE CURVA (CN)
RASTROJO Hileras Rectas 77 86 91 94
CULTIVOS EN HILERAS Hileras Rectas Mala 72 81 88 91
Hileras Rectas Buena 67 78 85 89
Curvas de Nivel Mala 70 79 84 88
Curvas de Nivel Buena 65 75 82 86
Cur/Niv y Mala 66 74 80 82
Terrazas
Cur/Niv y Buena 62 71 78 81
Terrazas
CULTIVOS EN HILERAS Hileras Rectas Mala 65 76 84 86
ESTRECHAS Hileras Rectas Buena 63 75 83 87
Curvas de Nivel Mala 63 74 82 85
Curvas de Nivel Buena 61 73 81 84
Cur/Niv y Mala 61 72 79 82
Terrazas
Cur/Niv y Buena 59 70 78 81
Terrazas
LEGUMBRES EN HILERAS Hileras Rectas Mala 66 77 85 89
ESTRECHAS O PASTOS EN Hileras Rectas Buena 58 72 81 85
ROTACIÓN 1/
Curvas de Nivel Mala 64 75 83 85
Curvas de Nivel Buena 55 69 78 83
Cur/Niv y Mala 63 73 80 83
Terrazas
Cur/Niv y Buena 51 67 76 80
Terrazas
PASTOS DE PASTOREO Mala 68 79 86 89
Regular 49 69 79 84
Buena 39 61 74 80
Curvas de Nivel Mala 47 87 81 88
Curvas de Nivel Regular 25 59 75 83
Curvas de Nivel Buena 6 35 70 79
PASTOS DE CORTE Buena 30 58 71 78
BOSQUE Mala 45 66 77 83
Regular 36 60 73 79

Página 104 de 156


COBERTURA Y USO DEL TRATAMIENTO CONDICIÓN GRUPOS HIDROLÓGICOS DE
SUELO O PRÁCTICA HIDROLÓGICA SUELOS
A B C D
NÚMERO DE CURVA (CN)
Buena 25 55 70 77
SOLARES ---------- 59 74 82 86
VÍAS DESTAPADAS 2/ ---------- 72 82 87 89
PAVIMENTOS ---------- 74 84 90 92
1/ Siembra tupida o al volé
2/ Incluyendo derecho de vía
Fuente: (Chow et al., 1994)

La condición hidrológica que hace referencia al grado de infiltración corresponde a la siguiente


interpretación:

Buena Para alta capacidad de infiltración, bajo potencial de escorrentía.


Regular Para regular capacidad de infiltración, regular potencial de escorrentía.
Mala Para baja capacidad de infiltración, alto potencial de escorrentía.

Para la condición hidrológica se hace referencia al grado de cobertura vegetal y puede establecerse
de la siguiente manera:

Buena Para cobertura de más del 75% del área.


Regular Con cobertura entre 50% y 75%.
Mala Cobertura menor del 50%.

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:

Tabla 34. Valores de CN para diferentes condiciones de humedad antecedente


CN PARA CONDICIONES CN PARA CONDICIONES
CONTINUACIÓN
II I III II I III
100 100 100 74 55 88
99 97 100 73 54 87
98 94 99 72 53 86
97 91 99 71 52 86
96 89 99 70 51 85
95 87 98 69 50 84

Página 105 de 156


94 85 98 68 48 84
93 83 98 67 47 83
92 81 97 66 46 82
91 80 97 65 45 82
90 78 96 64 44 81
89 76 96 63 43 80
88 75 95 62 42 79
87 73 95 61 41 78
86 72 94 60 40 78
85 70 94 59 39 77
84 68 93 58 38 76
82 66 92 56 36 75
81 64 92 55 35 74
80 63 91 54 34 73
79 72 91 53 33 72
78 60 90 52 32 71
77 59 89 51 31 70
76 58 89 50 31 70
75 57 88 49 30 69
Fuente: (Chow et al., 1994)

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:

[𝐶𝑁((𝑃/𝑅) + 2) − 200]2 200 Ec. 14


𝑃𝑒 = ;𝑃 ≥ 𝑅 ∗( ) − 2; 𝑅 = 2.54
𝐶𝑁 [𝐶𝑁((𝑃⁄𝑅 ) − 8) + 800] 𝐶𝑁

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. 15


𝐶𝑁𝑐𝑜𝑚𝑝 =
𝐴𝑡𝑜𝑡𝑎𝑙

Página 106 de 156


Las pérdidas iniciales de la lluvia se estiman a partir de la ecuación de hipótesis que define el método
del SCS:

𝐹𝑎 𝑃𝑒 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.

El tiempo de rezago corresponde al tiempo comprendido entre el centroide del histograma de


precipitación efectiva y el pico del hidrograma de escorrentía directa correspondiente. La experiencia
del SCS ha mostrado la siguiente relación típica para estimar el tiempo de rezago en una cuenca:

𝑡𝑙 = 0.6 · 𝑡𝑐 Ec. 18

Donde tl = tiempo de desfase de la hoya hidrográfica en horas y tc es el tiempo de concentración.


El hidrograma unitario de una cuenca es el hidrograma de escurrimiento directo, producido por 1
unidad (mm, cm, etc.) de lluvia en exceso, que cae con una intensidad uniforme sobre toda la cuenca
durante un tiempo determinado; es un hidrograma típico de la cuenca, que se relaciona con las
características de la cuenca (área, forma, pendiente, etc.) que son relativamente constantes.

El hidrograma unitario está basado en el supuesto de que precipitaciones de características similares


(en duración, intensidad, distribución, etc.) generan hidrogramas unitarios parecidos. Por lo tanto, las
ordenadas del hidrograma son proporcionales al volumen de escurrimiento directo.

Si no se cuenta con registros simultáneos de lluvia y escurrimientos, el hidrograma unitario se realiza


a partir de alguna de las características fisiográficas de la cuenca. Para esta situación se utilizan los
denominados hidrogramas unitarios sintéticos, como el hidrograma Unitario Adimensional del Servicio
de Conservación de Suelos (SCS).

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

Página 107 de 156


Unidos para diferentes condiciones climatológicas y áreas de drenaje. Para definir el hidrograma
unitario adimensional se calcula el caudal pico (qp) y el tiempo al pico (tp) con las expresiones:

A
q p = 0.208 ∗ t = √t c + 0.6 ∗ t c Ec. 19
tp p

Dónde: q pes el caudal pico, en m3/s/mm.


A es el área de la cuenca, en km2.
t c es el tiempo de concentración, en horas.
t p es el tiempo pico, en horas.

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.

Figura 61. Hidrograma Unitario Adimensional del SCS.

Fuente: Chow et al. (1994)

4.5.6. Modelación hidrológica

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

Página 108 de 156


cada subcuenca, para lo cual se emplearon diferentes escenarios, asociados a distintos períodos de
retorno, a la simultaneidad de eventos, a la probabilidad de ocurrencia y a las condiciones de frontera
y así determinar los caudales máximos de las crecientes correspondientes a los diferentes períodos
de retorno.

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).

El método L-Moments es una gran herramienta para la regionalización de datos hidrológicos,


permitiendo de esta manera obtener caudales de diseños en una zona de proyecto. Es un método
fundamentado en la estadística y es considerado en la actualidad como el método más preciso, porque
involucra no solo la supervisión de los datos, sino también la identificación de zonas homogéneas,
comprende además la selección de la distribución de frecuencia regional y la estimación de los
caudales de diseños (Ayros & Salazar, 2011).

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.

Página 109 de 156


Figura 62. Envolvente de creciente máxima Región 16 – Tr 100 años

Fuente: (Núñez Galeano, 2015)

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.

4.6. RESULTADOS HIDROLÓGICOS

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.

Página 110 de 156


4.6.1. Análisis de frecuencia de lluvias máximas en 24 horas y determinación de aguacero
de diseño

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]

Página 111 de 156


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
PUERTO
C15 8.20 CAICEDO 6.82 0.83 52.90 44.00
[47010110] 53.94
VILLAGARZON
C15 8.20 1.38 0.17 59.10 9.94
[44015010]
PUERTO
C16 30.85 CAICEDO 29.87 0.97 30.01 29.05
[47010110] 30.12
VILLAGARZON
C16 30.85 0.98 0.03 33.53 1.07
[44015010]
PUERTO
C17 27.95 CAICEDO 27.95 1.00 40.57 40.57 40.57
[47010110]
PUERTO
C18 7.70 CAICEDO 7.70 1.00 66.98 66.98 66.98
[47010110]
PUERTO
C19 13.57 CAICEDO 5.19 0.38 41.47 15.87
[47010110] 44.47
VILLAGARZON
C19 13.57 8.37 0.62 46.33 28.60
[44015010]
EL PEPINO - AUT
C2 157.35 34.08 0.22 25.83 5.59
[44015070]
PUERTO
C2 157.35 CAICEDO 2.20 0.01 25.89 0.36
[47010110] 27.45
TORRE TV SAN
C2 157.35 16.33 0.10 21.56 2.24
FCO [47010180]
VILLAGARZON
C2 157.35 104.74 0.67 28.93 19.26
[44015010]
PUERTO
C20 15.36 CAICEDO 15.36 1.00 48.40 48.40 48.40
[47010110]
PUERTO
C21 15.29 CAICEDO 14.09 0.92 47.00 43.31
[47010110] 47.43
VILLAGARZON
C21 15.29 1.20 0.08 52.51 4.13
[44015010]
VILLAGARZON
C22 11.69 11.69 1.00 49.50 49.50 49.50
[44015010]
PUERTO
C23 6.31 CAICEDO 6.31 1.00 80.07 80.07 80.07
[47010110]
SANTA ISABEL
C3 71.03 71.03 1.00 24.51 24.51 24.51
[47010230]
CHUNGACASPI
C4 82.88 4.67 0.06 22.72 1.28
[47010050]
EL PEPINO - AUT
C4 82.88 10.93 0.13 37.89 5.00
[44015070]
35.48
TORRE TV SAN
C4 82.88 40.22 0.49 31.62 15.34
FCO [47010180]
VILLAGARZON
C4 82.88 27.06 0.33 42.44 13.86
[44015010]
BALSAYACO
C5 29.06 17.50 0.60 39.74 23.93
[47010020]
38.38
SANTA ISABEL
C5 29.06 11.56 0.40 36.32 14.45
[47010230]
SANTA ISABEL
C6 38.63 38.63 1.00 32.18 32.18 32.18
[47010230]

Página 112 de 156


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
PUERTO
C7 29.33 CAICEDO 2.08 0.07 63.02 4.48
[47010110] 69.88
VILLAGARZON
C7 29.33 27.24 0.93 70.41 65.40
[44015010]
BALSAYACO
C8 41.59 12.84 0.31 29.94 9.24
[47010020]
PUERTO
C8 41.59 CAICEDO 22.72 0.55 58.84 32.15
[47010110] 46.87
SANTA ISABEL
C8 41.59 4.39 0.11 27.36 2.89
[47010230]
VILLAGARZON
C8 41.59 1.64 0.04 65.74 2.59
[44015010]
PUERTO
C9 39.63 CAICEDO 7.37 0.19 39.68 7.38
[47010110] 43.47
VILLAGARZON
C9 39.63 32.26 0.81 44.34 36.09
[44015010]
Fuente: Elaboración propia, 2022

4.6.2. Determinación de las crecientes para el periodo de retorno de diseño

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.

En detalle, el procedimiento utilizado para construir el hidrograma de la creciente de diseño es el


siguiente:

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.

Página 113 de 156


5. Debido a que la distribución temporal de los aguaceros en cualquier estación es aleatoria, se
adopta la siguiente curva que son los términos medios de acuerdo a las curvas de probabilidad
del 50% entre los Cuartiles I y IV los más comúnmente empleados (Back, 2011), las cuales son
usualmente empleadas para este tipo de proyectos para obtener una distribución temporal de la
precipitación. En la Tabla 47 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 y en la Figura 63 aparece la curva de distribución
adoptada en los presentes análisis.
Tabla 36. Distribución temporal de la lluvia adoptada.
FRACCIÓN FRACCIÓN
TIEMPO PRECIPITACIÓN
0 0
0.1 0.04
0.2 0.13
0.3 0.28
0.4 0.49
0.5 0.7
0.6 0.84
0.7 0.92
0.8 0.97
0.9 0.995
1 1,000
Fuente: (Chow et al., 1994)

Página 114 de 156


Figura 63. Curva de distribución del aguacero seleccionado en el presente estudio
DISTRIBUCCIÓN DE LA LLUVIA

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 [%]

Fuente: Chow et al. (1994)

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.

Página 115 de 156


Figura 64. Hietograma 100 años asociado a la Subcuenca 1
16

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)

Fuente: Elaboración propia, 2022

Figura 65. Hietograma 100 años asociado a la Subcuenca 2


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)

Fuente: Elaboración propia, 2022

Página 116 de 156


Figura 66. Hietograma 100 años asociado a la Subcuenca 3
14

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)

Fuente: Elaboración propia, 2022

Figura 67. Hietograma 100 años asociado a la Subcuenca 4


20

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)

Fuente: Elaboración propia, 2022

Página 117 de 156


Figura 68. Hietograma 100 años asociado a la Subcuenca 5
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)

Fuente: Elaboración propia, 2022

Figura 69. Hietograma 100 años asociado a la Subcuenca 6


18

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)

Fuente: Elaboración propia, 2022

Página 118 de 156


Figura 70. Hietograma 100 años asociado a la Subcuenca 7
35

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)

Fuente: Elaboración propia, 2022

Figura 71. Hietograma 100 años asociado a la Subcuenca 8


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)

Fuente: Elaboración propia, 2022

Página 119 de 156


Figura 72. Hietograma 100 años asociado a la Subcuenca 9
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)

Fuente: Elaboración propia, 2022

Figura 73. Hietograma 100 años asociado a la Subcuenca 10


30

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)

Fuente: Elaboración propia, 2022

Página 120 de 156


Figura 74. Hietograma 100 años asociado a la Subcuenca 11
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)

Fuente: Elaboración propia, 2022

Figura 75. Hietograma 100 años asociado a la Subcuenca 12


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)

Fuente: Elaboración propia, 2022

Página 121 de 156


Figura 76. Hietograma 100 años asociado a la Subcuenca 13
35

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)

Fuente: Elaboración propia, 2022

Figura 77. Hietograma 100 años asociado a la Subcuenca 14


40

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)

Fuente: Elaboración propia, 2022

Página 122 de 156


Figura 78. Hietograma 100 años asociado a la Subcuenca 15
30

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)

Fuente: Elaboración propia, 2022

Figura 79. Hietograma 100 años asociado a la Subcuenca 16


Hietograma precipitación 100 años - Eslabón
18

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)

Fuente: Elaboración propia, 2022

Página 123 de 156


Figura 80. Hietograma 100 años asociado a la Subcuenca 17
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)

Fuente: Elaboración propia, 2022

Figura 81. Hietograma 100 años asociado a la Subcuenca 18


35

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)

Fuente: Elaboración propia, 2022

Página 124 de 156


Figura 82. Hietograma 100 años asociado a la Subcuenca 19
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)

Fuente: Elaboración propia, 2022

Figura 83. Hietograma 100 años asociado a la Subcuenca 20


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)

Fuente: Elaboración propia, 2022

Página 125 de 156


Figura 84. Hietograma 100 años asociado a la Subcuenca 21
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)

Fuente: Elaboración propia, 2022

Figura 85. Hietograma 100 años asociado a la Subcuenca 22


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)

Fuente: Elaboración propia, 2022

Página 126 de 156


Figura 86. Hietograma 100 años asociado a la Subcuenca 23
40

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)

Fuente: Elaboración propia, 2022

8. Se caracteriza la cobertura vegetal de la cuenca hasta el sitio de estudio y el tipo de suelos


superficiales de la cuenca hasta esta sección para evaluar el parámetro CN asociado a la
cobertura.
9. El método de selección de las estaciones para las microcuencas, se partió de la información
disponible en la zona de estudio, que para este caso fue las estaciones Churuyaco y San Miguel.
10. Con los datos anteriores se selecciona el número de curva CN que representa el complejo
hidrológico suelo-cubierta vegetal de la cuenca, de acuerdo con el procedimiento mostrado. El
valor de curva número CN para la condición de humedad antecedente II se determinó en función
del tipo de suelo, cobertura vegetal y relieve. Estas características fueron determinadas mediante
el análisis de la información mostrada anteriormente, apoyada con visitas de campo efectuadas
con este propósito.
11. Se estima la lluvia efectiva, es decir, el escurrimiento directo con las siguientes fórmulas, usando
el valor “CN” seleccionado anteriormente, que representa el complejo hidrológico suelo-cubierta
vegetal para la condición más apropiada. El detalle de los cálculos realizados con el programa
HEC-HMS se muestra en el Anexo 02 del componente Hidrología.
12. Se calcula la longitud del cauce principal (cauce más largo) desde su nacimiento hasta el extremo
inferior de cada cuenca, la pendiente longitudinal media y el desnivel de la corriente entre su
nacimiento y la sección mencionada.

Página 127 de 156


13. Se determina el tiempo de concentración para la cuenca utilizando la pendiente longitudinal y los
desniveles correspondientes.
14. Se calculan los hidrogramas triangulares para cada incremento de escurrimiento en horas,
Para el intervalo de tiempo “D” seleccionado, se calcula el tiempo al pico, tp, el tiempo base, tb y
el caudal máximo, qp para 1 mm de escorrentía, es decir, de lluvia efectiva.
15. Se determinan los caudales pico o máximos de los hidrogramas triangulares para cada incremento
de escurrimiento. Estos máximos se obtienen multiplicando el caudal máximo de 1 mm de
escurrimiento (calculado en el paso anterior) por el incremento de lluvia efectiva (escurrimiento).
Se parte de los incrementos de tiempo “D” (en horas) definidos anteriormente y de los
incrementos de la lluvia efectiva (escurrimiento) en mm. Se consideran los valores de los caudales
máximos (caudales pico) de cada uno de los hidrogramas unitarios en m3/s/mm, es decir, caudales
producidos por 1 mm de lluvia efectiva, determinados de manera idéntica para cada intervalo de
tiempo con la ecuación:

0,2083 ⋅ 𝐴 Ec. 20
𝑞𝑃 =
𝑡𝑃

qP: caudal pico del hidrograma triangular, en m3/s/mm


A: área de la cuenca aferente, en km2
tP: tiempo al pico del hidrograma unitario triangular, en horas

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.

Página 128 de 156


Figura 87. Esquema hidrológico generado para el modelo de la cuenca San Juan

Fuente: Elaboración propia, 2022

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.

Página 129 de 156


Figura 88. Hidrograma de creciente del nodo Juntion-11 para periodo de 100 años

Fuente: Elaboración propia, 2022

5. SIMULACIÓN HIDRÁULICA

La simulación hidráulica es la materialización espacial y numérica del comportamiento de las variables


físicas de la corriente en el cauce para un caudal de diseño. Este modelo comprende la espacialización
de las zonas inundables, la socavación general del Talweg, la capacidad de cizallamiento, velocidad,
energía, y tirante hidráulico.

Página 130 de 156


Aprovechando las potencialidades GIS del desarrollo metodológico propuesto, la simulación de los
caudales obtenidos en el aparte anterior se realizó mediante las herramientas RAS Mapper de HEC-
RAS 6.1 (Ackerman, Wre, Jensen, & Brunner, 2010). Según (US Army Corps of Engineers, 2016),
HEC-RAS resuelve el flujo gradualmente variado a partir de la ecuación de balance de energía
(trinomio de Bernoulli) entre dos secciones dadas, excepto en los casos en los que se simules
estructuras como puentes, vertederos o tramos cortos entubados (culverts). En tales casos HEC-RAS
resuelve la ecuación de conservación de la cantidad de movimiento, así como ciertas ecuaciones de
carácter empírico establecidas ad hoc para estas estructuras.

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

Página 131 de 156


ser ortogonal, sin embargo, si fuera así la discretización numérica es más simplificada y
eficiente.
• Tabla detallada de propiedades hidráulicas para celdas y contornos de celdas computaciones
2D. Cada celda y contorno de celda es pre-procesada con la finalidad de obtener tablas de
propiedades hidráulicas basadas en el terreno subyacente empleado en la modelización.
Básicamente el pre-proceso, calcula una relación detallada de Elevación-Volumen para cada
celda; y para cada contorno de celda calcula la relación Elevación-Perímetro Mojado,
Elevación-Área, Elevación-Rugosidad, y demás propiedades hidráulicas. Estas relaciones que
crea el programa permiten al usuario crear celdas computacionales grandes conservando los
detalles del terreno; lo cual resulta ventajoso porque hace más rápido los tiempos de cálculo,
pues genera mayores detalles hidráulicos a nivel de cada celda. La elección del tamaño de
celda se basa en el nivel de detalle que se quiere obtener.
• Mapas detallados de inundación y animaciones de los mismos. HEC-RAS, a través de su
herramienta RAS Mapper, ofrece la posibilidad de visualizar mapas de zonas inundadas, así
como la animación del flujo de agua cuando acontece la inundación. Dicho proceso se basa
en función del terreno subyacente y no en el tamaño de celda computacional de la malla
generada.
• Solución del esquema numérico basado en multiprocesadores. La solución obtenida de la
modelización, ha sido programada para que aproveche los procesadores que tiene el
ordenador, lo que permite que se ejecute más rápido si solo empleáramos uno.
• Permite trabajar con procesadores de 64-Bit y 32-Bit HEC-RAS puede trabajar ahora en
ordenadores que tengan procesadores de 32 y 64 bit, siendo conocido que un procesador de
64-Bit se ejecutará más rápido que el de 32-Bit y puede manejar gran cantidad de datos.

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:

Conservación de la masa: Suponiendo que el flujo es incompresible, la forma diferencial inestable de


la ecuación de conservación de masa (continuidad) es:

𝜕𝐻 𝜕(ℎ𝑢) 𝜕(ℎ𝑣) Ec. 21


+ + +𝑞 =0
𝜕𝑡 𝜕𝑥 𝜕𝑧

Página 132 de 156


donde t es el tiempo, u y v son los componentes de velocidad en la dirección x e y respectivamente y
q es un término de flujo fuente / sumidero.
En forma vectorial, la ecuación de continuidad toma la forma: ∂H.

𝜕𝐻 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

Página 133 de 156


de análisis y su aplicación al problema directo. Caso contrario son los modelos estocásticos,
que a cierto conjunto de datos de entrada le corresponden diferentes conjuntos de salida con
una probabilidad de ocurrencia, por lo que la variable de estado es de carácter aleatorio
• El segundo criterio se centra en la capacidad del modelo en determinar las variaciones
temporales (evolución) del proceso de análisis, de ser así estaríamos hablando de un modelo
dinámico o no estacionario, en los que el tiempo es una variable independiente. Si el modelo
carece de esa propiedad hablamos de un modelo estático o estacionario.
• El tercer criterio analiza el número de dimensiones en las que se soluciona la variable de
estado y los parámetros del modelo. Si la estructura del modelo matemático solo tiene
presente los cambios en una dimensión, como es el caso del tiempo, estaremos hablando de
modelos aglutinados, mientras que, si las soluciones tienen en cuenta dos o más dimensiones,
estableciendo los cambios de las variables en el tiempo en una coordenada o dos, definiendo
esto como modelos distribuidos.
Centrando el interés a modelaciones de la propagación de avenidas en ríos la estructura del modelo
se centra en resolver las ecuaciones del flujo variable en lámina libre o ecuaciones de Saint Venant,
las cuales se deducen las ecuaciones de las leyes físicas de conservación de un fluido (Bladé, 2009).
Dichos modelos parten de la simplificación del problema y eso repercute en las ecuaciones de
gobierno que tienen una precisión adecuada y minimizan el costo computacional, ya que dichas
ecuaciones parten de las establecidas por Navier – Stokes, las cuales describen el movimiento
instantáneo y las variables medias del teorema de Reynolds y su discretización es tridimensional a lo
largo del dominio de estudio. Las simplificaciones se centran primordialmente en integrar la
profundidad de las ecuaciones Reynolds, eliminado la dimensión vertical, obteniendo las ecuaciones
en dos dimensiones (2D); el siguiente paso es llevarlas a una dimensión (1D), siendo ecuaciones más
simples que pueden representar el movimiento de la onda en no permanente.

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

Página 134 de 156


grande sea el paso de tiempo es más sencilla la descripción de los procesos hidrológicos. Al contrario,
la modelación de los procesos a un paso de tiempo corto, necesita una descripción de los procesos
más complejos (Gnouma 2006). La selección del modelo está condicionada principalmente a los
objetivos del estudio y a la disponibilidad de la información. Otros criterios que inciden en la selección
del modelo son: los recursos económicos, técnicos, computacionales y el tiempo disponible para llevar
a cabo el estudio.

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

Consiste en definir y precisar la representación del sistema fluvial en el modelo. En él se seleccionan


tanto el espaciamiento como las secciones transversales del cauce y la planicie que mejor representan
la batimetría y la topografía, respectivamente, se define el tamaño de las celdas de la malla
computacional para representar convenientemente la topografía de toda el área de estudio. Para ello
primero se definen las longitudes de los cauces a modelar y se delimita de manera precisa la extensión
del área a modelar. Es conveniente señalar aquí la necesidad de ubicar los extremos o fronteras.

También se definen en la etapa de esquematización las condiciones hidrodinámicas en las fronteras


del modelo; igualmente se definen los valores de la rugosidad de los cauces y de la llanura de
inundación, la viscosidad turbulenta o alguna ecuación para determinarla, según los recursos
disponibles en el mismo modelo.

5.1.2. Implementación del modelo

La implementación del modelo consiste en suministrar al modelo la siguiente información:


1 Información física esquematizada anteriormente según los requerimientos y herramientas
disponibles en el mismo.
2 Definir la topografía en el modelo hidráulico mediante el modelo digital de terreno (MDT).
3 Los valores de los parámetros físicos (rugosidad, viscosidad turbulenta en todos los tramos y
celdas del modelo.
4 Las condiciones hidrodinámicas en las fronteras del modelo, generalmente las hidrógrafas de
caudales en las fronteras superiores de los diferentes cauces en el área de estudio, y la curva de

Página 135 de 156


calibración nivel-caudal o las hidrógrafas de niveles en las fronteras inferiores o de salida del
modelo.
5 Finalmente se seleccionan los parámetros numéricos del modelo. Principalmente el intervalo de
tiempo computacional y el espaciamiento máximo en los cauces. Estos parámetros se establecen
teniendo en cuenta diferentes aspectos, tales como, el grado de detalle o resolución requerido
en los resultados del modelo, el tiempo computacional (no debe ser excesivamente grande) y
especialmente, garantizar la estabilidad numérica del modelo y la precisión de los resultados. En
este sentido, cada modelo matemático define, dependiendo del tipo de esquema numérico
adoptado para resolver las ecuaciones hidrodinámicas, las condiciones que deben cumplirse
para garantizar su estabilidad numérica; por ello, siempre debe chequearse el estricto
cumplimiento de las mismas. Una vez concluido este proceso se realizan la primeras corridas o
simulaciones en el modelo, siguiendo, preferiblemente y en su orden, las etapas de análisis de
sensibilidad, calibración, verificación y explotación o aplicaciones del modelo para generar
finalmente los niveles de agua y las condiciones hidrodinámicas para los eventos hidrológicos de
interés del estudio, es decir, las crecientes correspondientes a los diferentes periodos de retorno
seleccionados.

5.1.3. Análisis de sensibilidad

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.

5.1.4. Verificación del modelo

El termino validación no necesariamente denota el establecimiento de la verdad. Preferiblemente


denota el establecimiento de legitimidad típicamente dada en términos de argumentos, métodos y
convenciones (Merriam 1963). La validación permite corroborar el desempeño del modelo en un
contexto espacio – temporal diferente al de la calibración; ésta puede ser espacial (en un punto
diferente de donde se calibró), temporal (empleando series con un periodo de tiempo diferente al
calibrado) y espacio – temporal (en un punto y en un periodo diferente a los empleados en la
calibración); en todos los casos anteriores, se corre el modelo ya calibrado (IDEAM 2017). En síntesis,
el modelo se calibrará y verificará para todos los periodos y condiciones hidrológicas que se
consideren importantes en el estudio.

Página 136 de 156


5.1.5. Aplicaciones del modelo

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.

5.1.6. Condiciones de frontera

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.

5.2. IMPLEMENTACIÓN DEL MODELO

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.

5.2.1. Configuración del dominio 1D

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.

Página 137 de 156


Posteriormente, se generó el Modelo De Elevación digital (DEM) creado por el equipo de SIG, el
modelo posee una resolución de 0.10x0.10m que fue generado con información primaria de la
agregación de la topobatimetría.

5.2.2. Condiciones de borde e iniciales del modelo

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.

Figura 89. Esquema de condiciones de borde de la simulación hidráulica

El esquema en gris indica la estructura propuesta


Fuente: Elaboración propia, 2022

Página 138 de 156


5.2.3. Coeficientes de rugosidad

Cuando se requiere la definición de los condicionantes hidráulicos de un modelo hidráulico, ya sea


unidimensional o bidimensional, una de las variables que más influyen en el comportamiento del
modelo es precisamente el coeficiente de rugosidad.

Afortunadamente, desde hace muchos años, se ha estado estudiando y analizando diferentes


situaciones de flujo obteniendo de este modo una cierta relación entre unas características físicas de
cauces y márgenes y unos valores de coeficiente de Manning representativos.

Para seleccionar el coeficiente de rugosidad de Manning en cauces naturales intervenidos y no


intervenidos, se dispone de mucha literatura donde se puede encontrar toda información previa de
valores ya establecidos, verificados y tabulados, es realmente de mucha ayuda, pero es claro, esos
valores responden a unas características muy concretas que seguramente presentará similitudes con
las que tiene el cauce que debes modelizar, pero al final de cada caso es único. Por ello es importante
entender los factores que afectan el valor n para adquirir un conocimiento básico del problema y
disminuir la incertidumbre.

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,

Página 139 de 156


tablas tomadas de Chow (2004) e imágenes de Barnes (1967) y proceder o a estimar los valores de
la n de Manning para cada sector.

El coeficiente de resistencia se determinará con la siguiente expresión, similar a la de Cowan (Ochoa


Tomás. 2002):

𝑛 = (𝑛0 + 𝑛1 + 𝑛2 + 𝑛3 + 𝑛4 + 𝑛5 + 𝑛6 )𝑚1 ⋅ 𝑚2 Ec. 23

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.

5.3. RESULTADOS DE LA MODELACIÓN

Los resultados obtenidos en la modelación hidráulica unidimensional del proyecto, se exponen a


continuación por cada uno de los ítems de análisis:

5.3.1. Esquema de planta de inundación

Esta categoría simplifica el alcance y zonas de desborde la inundación sobre el área de análisis.

Página 140 de 156


Figura 90. Esquema en planta de la inundación

El esquema en gris indica la estructura propuesta


Fuente: Elaboración propia, 2022

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.

La altura máxima de la inundación es de 5.98 m, y esta se alcanza precisamente sobre la solera de


socavación lateral derecha axial al ponteadero. No obstante, de acuerdo con la proyección del nuevo
ponteadero, se tendría 1.86 m de altura efectiva después del N.A.M.E., localizándose en las cotas de
implante, es decir, sobre 297 msnm.

Página 141 de 156


Figura 91. Perfil hidráulico de la inundación por laterales y canal

C: Canal, L: Izquierda, R: Derecha


Fuente: Elaboración propia, 2022

5.3.2. Esquema de velocidades

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.

Página 142 de 156


Figura 92. Esquema en planta de velocidades

El esquema en gris indica la estructura propuesta


Fuente: Elaboración propia, 2022

5.3.3. Tirante hidráulico

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.

Página 143 de 156


Figura 93. Esquema del NAME vs Puente propuesto - Vista aguas arriba
SanJuan
RS = 73 BR

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)

Fuente: Elaboración propia, 2022

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)

Fuente: Elaboración propia, 2022

Página 144 de 156


5.4. ANÁLISIS DE SOCAVACIÓN

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:

• Suelos cohesivos (m/s)

𝑉𝑒 = 0.60 𝑔𝑑 1.18 𝑏 𝐻𝑠 𝑥 Ec. 24

• Suelos no cohesivos (m/s)

𝑉𝑒 = 0.68 𝑑𝑚 0.28 𝑏 𝐻𝑠 𝑥 Ec. 25

Donde,

Ve = Velocidad media suficiente para degradar el cauce en m/seg.


gd = Peso volumétrico del material seco que se encuentra a una profundidad Hs, medida desde la
superficie del agua (Ton/m3)
b = Coeficiente que depende de la frecuencia con que se repite la avenida que se estudia.
x = Es un exponente variable que está en función del peso volumétrico gs del material seco (Ton/m3)

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í.

Página 145 de 156


( 1 ) - Perfil antes de la erosión.
( 2 ) - Perfil después de la erosión

Cálculo de la profundidad de la socavación en suelos homogéneos:

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

Página 146 de 156


granulométrica de suelos no cohesivos, por tanto, el cálculo del diámetro de material seleccionado
para la presente socavación responde a los resultados se indica a continuación:

Tabla 37. Granulometría para análisis de socavación


LUGAR D60 D50 ESTE NORTE
SONDEO DERECHA 0.37 0.33 1042166.3 573886.6
SONDEO IZQUIERDA NA NA NA NA
CAUCE (RECOLECCIÓN DE MUESTRAS) --- 0.87 1042224.1 573863.3

Fuente: Elaboración propia, 2022

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:

Tabla 38. Valores de socavación obtenidos en el modelo hidráulico


Socavación por contracción
Variable Izquierda Canal Derecha
Ys (m): 0 0 0
Vc (m/s): 10.6 0.86 0.53
Ecuación Clear Clear Clear
Socavación de muelle
Todos los muelles Ys (m): 1.77
Froude (#): 0.08
Ecuación: CSU
Socavación de pilar
Variable Izquierda Derecha
Pilar Ys (m): 0.45 1.84
Froude (#): 0.06 0.18
Ecuación: CSU
Socavaciones combinadas
Muelle + Contracción (m):
Muelle + Contracción Izquierda (m): 3.61
Muelle + Contracción Derecha (m): 0.45
Fuente: Elaboración propia, 2022

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.

Página 147 de 156


La ampliación del cauce se daría en 6.08 m, estos datos son los requeridos para el diseño de las obras
de protección y limitación de impacto – divagación de la estructura en la dinámica fluvial, en ese
sentido, la luz efectiva será de por lo menos 103 metros como lo expone el diseño hidráulico.

Figura 95. Esquema frontal de socavación lateral y general de la estructura del puente

Fuente: Elaboración propia, 2022

5.5. DISEÑO FLUVIAL

La estabilidad de la obra en las condiciones de diseño presentadas por el componente geotécnico y


estructural no permanece sustentable sin que se realicen intervenciones directas en el cauce que
garanticen la limitación del efecto de socavación y migración del cauce.

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.

Página 148 de 156


• Protección de estribos

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.

Figura 96. Esquema conceptual de protección de estribos

Fuente: Geomatrix, 2021

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:

Figura 97. Arreglos tipo de Hidromalla Geoestera 300

Fuente: Geomatrix, 2021

La propuesta se protección se orienta básicamente a la torre derecha y el muerto de la misma


orientación, pues son aquellos que mantienen condiciones se socavación, por como se describió en
la sección anterior. El tipo de arreglo dado es el siguiente:

Página 149 de 156


Figura 98. Arreglo de hidromallas para protección de estribo y torre derecha

Fuente: Adaptado de Geomatrix, 2022

• Protección por migración

Para la migración se ha propuesto la implementación de protección anti-socavación utilizando el


producto de colchacreto articulado.

El diseño propuesto para la protección se encuentra adicionado en el Anexo


\\Puente_Vides\06_INFORME\ Proteccion PVides Derecha.XLSM

Figura 99. Esquema conceptual de protección de migración por socavación

Fuente: Geomatrix, 2021

La guía de instalación se anexa al folder de cálculo de socavación.

Página 150 de 156


CONCLUSIONES Y RECOMENDACIONES

Los criterios de diseño que debe satisfacer la infraestructura son los siguientes:

• Periodo de retorno: 100 años


• Caudal: 725.5 m3/s
• Régimen de flujo: Mixto
• Tirante crítico sobre tablero: 1.89 m
• Distancia crítica del tirante a la base del puente: 1.86m
• Velocidad de flujo: 1.4 m/s
• Socavación derecha: 0.45 m
• Socavación izquierda: 3.61 m
• Socavación de lecho: 1.77 m

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:

Tabla 39. Centrado de obras para replanteo


ÍTEM ESTE NORTE ELEVACIÓN
PILA DERECHA 1042164.7339 573883.1847 267
PILA IZQUIERDA 1042266.4118 573856.9814 267
ESTRIBO DERECHO 1042155.6760 573885.5482 267-Criterio Estructural
ESTRIBO IZQUIERDO 1042278.6907 573853.8171 267-Criterio Estructural
Fuente: Elaboración propia, 2022

El esquema de implante se presenta en la siguiente Figura 100:

Página 151 de 156


Figura 100. Implante del Puente en el Fotomosaico de referencia

Fuente: Elaboración propia, 2022

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.

Figura 101. Implante de protección anti-socavación local

Fuente: Elaboración propia, 2022

Página 152 de 156


c) Se recomienda realizar la descolmatación de la barra de acceso al asentamiento, de manera
que el material presente se pueda agradar en la conformación del talud que soportará la
medida de protección anti-socavación lateral. Propuesta que se hace mediante el uso de
colchacreto articulado.

Figura 102. Implante de protección anti-socavación lateral

Fuente: Elaboración propia, 2022

d) Se recomienda al geotecnista desarrollar análisis condicentes a la verificación de la estabilidad


al desgaste o meteorización de la lodolita del implante izquierdo de la infraestructura.

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.

Página 153 de 156


h) El gestor del proyecto deberá realizar el trámite respectivo ante la autoridad ambiental para
obtener el permiso de ocupación del cauce por la obra a construir.

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.

Página 154 de 156


Gonzales Fletcher, A., Proyecto Radargrametrico del Amazonas, B. (Colombia) spa, Rivera Hermida, H.,
Diazgranados, D. A., & Jaramillo Berrio, G. (n.d.). La Amazonia colombiana y sus recursos; estudio
elaborado por el Grupo Tecnico del Proyecto Redargrametrico del Amazonas durante el periodo junio
1974-diciembre 1979. Bogota (Colombia) PRORADAM.
Guzman, D., Ruíz, J. ., & Cadena, M. (2014). Regionalización de Colombia según la estacionalidad de la
precipitación media mensual, a través análisis de componentes principales. In Grupo de
Modelamiento de Tiempo, Clima y Escenarios de Cambio Climático (Vol. 1).
https://doi.org/10.1016/s1013-7025(09)70018-1
Horton, R. E. (1945). Geological Society of America Bulletin. Geological Society Of America Bulletin, 56(1),
151–180. https://doi.org/10.1130/0016-7606(1945)56
IDEAM. (2010a). Leyenda nacional de coberturas de la tierra. Metodología CORINE Land Cover Adaptada
para Colombia. Escala 1:100.000. (N. Martinez & U. Murcia, Eds.). Retrieved from
http://siatac.co/c/document_library/get_file?uuid=a64629ad-2dbe-4e1e-a561-
fc16b8037522&groupId=762
IDEAM. (2010b). Sistemas Morfogénicos del Territorio Colombiano. 7 El dominio amazónico. Bogotá D.C.
IDEAM. (2018). Protocolo De Modelación Hidrológica E Hidráulica. Retrieved from
http://documentacion.ideam.gov.co/openbiblio/bvirtual/023833/Protocolo_Modelacion_HH.pdf
IGAC. (2014). Estudio general de suelos y zonificación de tierras Departamento del Putumayo Escala 1:
100.000. In Imprenta Nacional de Colombia (Ed.), Instituto Geográfico Agustín Codazzi (Vol. 1).
https://doi.org/10.1016/B978-84-458-1898-5.50004-3
INGEOMINAS. (2002). Memoria explicativa de la Geología de la Plancha 303, Colombia, Departamentos de
Huila, Tolima y Meta. Ingeominas Bogotá, 88.
Instituto de Hidrología Meteorología y Estudios Ambientales. (2010). Sistemas morfogénicos del territorio
colombiano. Bogotá D.C.
INVIAS, I. N. de V., & Ministerio de Transporte. (2009). Manual de Drenaje para Carreteras (J. Flores & A.
Bolaños, Eds.). Bogotá.
Karalis, S., Karymbalis, E., Valkanou, K., Chalkias, C., Katsafados, P., Kalogeropoulos, K., … Bofilios, A.
(2014). Assessment of the Relationships among Catchments ’ Morphometric Parameters and
Hydrologic Indices. International Journal of Geosciences, 2014(December), 1571–1583.
Keller, E. A., & Rockwell, T. K. (1984). Tectonic Geomorphology , Quaternary Chronology , and
Paleoseismicity. (November 2015). https://doi.org/10.1007/978-3-642-69759-3
Kondolf, G. M., Montgomery, D. R., Piégay, H., & Schmitt, L. (2005). Geomorphic Classification of Rivers
and Streams. In Tools in Fluvial Geomorphology (Vol. 9). https://doi.org/10.1002/0470868333.ch7
Lagasse, P. F., Spitz, W. J., Zevenbergen, L. W., & Zachmann, D. W. (2004). Handbook for Predicting Stream
Meander Migration and Supporting Software. https://doi.org/10.17226/23346
Lane, E. W. (1957). A study of the shape of channels formed by natural streams flowing in erodible material.
U.S. Army Corps of Engineers Report, (9), 141.
Leopold, L. B., Wolman, M. G., & Miller, J. P. (1964). Fluvial processes in geomorphology (Primera; W.F.
Freemanand Company, Ed.). Retrieved from www.doverpublications.com
Martín, M. (2017). ÍNDICE SINUOSIDAD DEL CAUCE | Cursos de Ingeniería, Medio Ambiente y Calidad.
Retrieved April 24, 2020, from http://eimaformacion.com/indice-de-sinuosidad-de-un-cauce/
Martín Vide, J. P. (2003). Ingeniería de Rios (p. 333). p. 333.
Martín Vide, J. P. (2004). Introducción a la ingeniería fluvial. 3o IOP. ETSECCP.
https://doi.org/10.1157/13068212

Página 155 de 156


Ministerio de Transporte, & AIS. INVIAS. (2014). Sección 2: Caracteristicas generales de diseño y
ubicación. Norma Colombiana de Diseño de Puentes, CCP 14, 25.
Ministerio de Transporte, & INVIAS, I. N. de V. (2015). Guía de diseño de pavimentos con placa-huella. 101.
Retrieved from
ttps://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&ved=0ahUKEwjZ-
ua32vvVAhWEPCYKHUFABZYQFggmMAA&url=https://www.mintransporte.gov.co/descargar.php?id
File=15033&usg=AFQjCNHVSOkDmudk0TOtlLyk8ZDYRfSMgw
Ministerio para la transición Ecológica. (2019). Protocolo de caracterización hidromorfológica de masas de
agua de la categoría ríos (M-R-HMF-20). Madrid, España.
Núñez Galeano, L. (2015). Adaptación del metodo L-Moments para la regionalizacion de eventos máximos
para las cuencas de Colombia. Pontificia Universidad Javeriana.
Robertson, K. G., Jaramillo, O., Castiblanco, M. A., & Instituto de Hidrología Meteorología y Estudios
Ambientales. (2013). Guía Metodológica Para La Elaboración De Mapas Geomorfológicos a Escala
1:100.000. Instituto de Hidrología, Meteorología y Estudios Ambintales - IDEAM, 88. Retrieved from
http://fs03eja1.cormagdalena.com.co/nuevaweb/Niveles/Definiciones.pdf
Rocha Felices, A. (1998). Introducción a la Hidráulica Fluvial. https://doi.org/10.1007/s13398-014-0173-7.2
Rojo, J. (1997). Morfometria de Cuencas.
Salvatierra-Colindres, F. (2006). ESTUDIO DE LA CUENCA DEL RÍO CAJOLÁ CON FINES DE MANEJO DE
SUELOS, MUNICIPIO DE CAJOLÁ, DEPARTAMENTO DE QUETZALTENANGO, GUATEMALA (pp. 1–
19). pp. 1–19. Guatemala: Universidad de San Carlos de Guatemala.
Schumm, S. A. (1960). The shape of alluvial channels in relation to sediment type. Office, 24.
Schumm, S. A. (1963a). A tentative classification of alluvial river channels. Geological Survey Circular 477,
16. Retrieved from
http://scholar.google.com/scholar?hl=en&btnG=Search&q=intitle:A+Tentative+Classification+of+Allu
vial+River+Channels#1
Schumm, S. A. (1963b). Sinuosity of Alluvial Rivers on the Great Plains. (September 1963), 1–11.
Schumm, S. A. (1985). Patterns of alluvial rivers. 5–28.
Trujillo, A. R., Barroso, F. U., & Escobar, Y. C. (2010). Guía básica para la caracterización morfométrica de
cuencas hidrográficas (1st ed.). Retrieved from http://www.jstor.org/stable/j.ctt1rfsqf5
UNAD. (2012). Patrones de Drenaje. Retrieved June 6, 2016, from
http://datateca.unad.edu.co/contenidos/30172/MODULO
HIDROLOGIA/leccin_12_patrones_de_drenaje.html
Universidad del Cauca. (2006). Capítulo 1. Geomorfología Fluvial.
US Army Corps of Engineers. (2016). HEC-RAS 5.0 Reference Manual (p. 960). p. 960.
Vargas, R., & Diaz, M. (1998). Curvas sintéticas regionalizadas de Intensidad-Duración- Frecuencia para
Colombia. 200. https://doi.org/10.1017/CBO9781107415324.004
Verstappen, H. T., Zuidam, R. A. van., Meijerink, A. M. J., & Nossin, J. J. (1991). The ITC system of
geomorphologic survey : a basis for the evaluation of natural resources and hazards. Enschede, the
Netherlands: International Institute for Aerospace Survey and Earth Sciences (ITC).
Villota, H., & Instituto Geográfico Agustin Codazzi (IGAC). (2005). Geomorfología Aplicada a Levantamientos
Edafológicos y zonificación Física de Tierras. Instituto Geográfico Agustín Codazzi, 184.
https://doi.org/10.1080/1350176032000046967

Página 156 de 156

También podría gustarte