CAPITULO I

INTRODUCCION
El conocimiento de los procesos hidrodinámicos y morfológicos de una corriente aluvial,
dadas unas condiciones hidrológicas conocidas, permite desarrollar un adecuado plan de
manejo de los diseños hidráulicos que se lleven a cabo en dicho cauce. Entre los
procesos fluviales, se destacan tres procesos como elementales: la erosión, el transporte y la
depositación de sedimentos.
La erosión fluvial puede darse tanto en las bancas como en el fondo, dependiendo los grados
de libertad del cauce (Maza Álvarez & García Flores, 1996). La erosión natural del fondo en
cauces naturales producida por un incremento del caudal líquido, sea por crecidas naturales o
inducidas, es llamada socavación general, y la estimación de su profundidad es un tema de
gran importancia en ingeniería fluvial, especialmente útil para la toma de decisiones en el
diseño de obras de cruce a través del cauce, tales como puentes, tuberías, túneles
subfluviales, etc. (Farias, 2003).
La socavación general es un proceso físico complejo, cuyas bases teóricas aún no se
encuentran bien definidas y no es posible estimar con confianza los cambios en el lecho tras
el tránsito de un flujo. Los factores que influyen en la socavación general se pueden agrupar
en tres grupos: factores geomorfológicos, factores de transporte y el tipo de material que
conforma el lecho; todos factores únicos para cada río, lo cual dificulta la existencia de una
ley general (Melville & Coleman, 2000).
Uno de los factores fundamentales para estimar la profundidad de socavación en
cualquier tipo de cauce natural es el tipo de material del lecho que conforma el cauce y está
expuesto a las fuerzas erosivas de la corriente. El tipo de material que conforma el cauce
hace referencia al tamaño del mismo, e investigaciones han demostrado que en ríos
aluviales, es decir, aquellos que corren sobre materiales transportados por el propio río en el
pasado geológico, la granulometría del material del lecho está directamente relacionada con
las características geomorfológicas de un corriente (Autumn & Springer, 2006). Por lo tanto,
el estudio de socavación se puede acotar a ríos de montaña (característica geomorfológica),
en los cuales el tamaño medio de las partículas que conforman el lecho son superiores a 2
mm, y en general, corresponden a corrientes con buena capacidad de transporte, valles
encañonados y con pendientes longitudinales superiores a 0.1% (Leopold, Wolman & Miller,
1964). Este tipo de cauces son característicos en la geografía del Valle de Aburrá.
Los primeros pasos en la hidráulica fluvial para determinar la profundidad de socavación
partieron de la teoría del régimen, formulaciones completamente empíricas entre las que se
destacan la de Lacey (1930) y Blench (1939).
Posteriormente, a partir de la década de los 60, se le dio un nuevo enfoque a las
metodologías desarrolladas para calcular la profundidad de socavación, basadas en el
equilibrio entre la velocidad de flujo y la velocidad crítica para el material que conforma el
lecho, entre las que se destacan las formulaciones de Lischtvan & Lebediev (1959), Laursen
(1963), Maza Álvarez (1973) y Neill (1980), entre otras. Todas estas formulaciones semi-
empíricas, en su desarrollo conceptual, parten de hipótesis que simplifican el problema,
entre las que se destaca la asunción de flujo permanente, situación que durante una
creciente no ocurre y que no tiene en cuenta que la iniciación del movimiento de las
partículas de fondo ocurre luego de un denominado tiempo de “stress”. Asimismo, la
velocidad crítica es calculada a partir de coeficientes empíricos que representan las
condiciones de los cauces donde fueron calibrados y se debe verificar su aplicación en otras
cuencas, sin mencionar que el efecto de escondimiento de las partículas no es tenido en
cuenta por ninguna de estas metodologías (Schreider et al.,
2001. Yager, Kirchner & Drietrich, 2007).
Por otro lado, con el desarrollo de herramientas de modelación hidrodinámica que
incluyen módulos de transporte de sedimentos y cambios en el lecho (como el HEC-RAS o
el SRH-1D), es posible modelar procesos de agradación y degradación del lecho a partir
de las ecuaciones de transporte y de continuidad de sedimentos.
Bajo el panorama expuesto, este trabajo presenta una comparación entre las formulaciones
empíricas para estimar la profundidad de socavación general y los resultados
obtenidos mediante la modelación morfo-dinámica en un río de gravas, como el río Medellín,
y propone una metodología para el cálculo de la profundidad de socavación general ocurrida
durante una creciente a partir de modelos computacionales que tienen en cuenta el tránsito
de flujo no permanente, la gradación del lecho, el factor de escondimiento y condiciones
para la iniciación del movimiento (Melville & Coleman,
2000), como el HEC-RAS o el SRH-1D.
Este documento se distribuye de la siguiente forma, en el capítulo 1 se hace una
descripción de conceptos teóricos de la socavación, así como su formulación matemática
desde varios enfoques.
En el capítulo 2 se presentan las características geomorfométricas,
hidrométricas y sedimentológicas de la captacion analizada, y los escenarios de
modelación propuestos.
en el capítulo 3 se hace un a n á l i s i s d e sensibilidad de los parámetros de los modelos
evaluados, y posteriormente,
en el capítulo 4 se muestran los resultados de las modelaciones de socavación mediante
los diferentes métodos propuestos para su comparación, y así proponer la metodología más
apropiada para modelar el proceso de socavación, descrita en el capítulo 5. Finalmente se
presentan las limitaciones que se tuvieron para mejorar los resultados obtenidos, y las
conclusiones de este trabajo.
La investigación es aplicada – cuantitativa - de campo y laboratorio, donde la información de
campo se obtuvo de registros técnicos de proyectos generados en los últimos años, y la
determinada in situ en el desarrollo de la investigación; la aleatoriedad de las descargas
máximas y la falta de registros no limita la investigación en el sentido de que se utilizó como
caudal de prueba el máximo probable que pueda soportar la estructura con el mayor
estrechamiento.

1.2 OBJETIVOS
1.2.1 Objetivo General

1.2.2 Objetivos Específicos

 CAPÍTULO II
REVISIÓN DE LITERATURA
Marco Conceptual
ESTUDIO DE SOCAVACION EN LA CAPTACION-LOS DADOS
2.1 SOCAVACIÓN
La socavación es un tipo de erosión hídrica que hace referencia a la pérdida del material del
lecho y márgenes de un cauce, debido a la capacidad de transporte asociada a un evento
hidrológico. La reducción de este nivel respecto a un nivel de referencia es denominada
profundidad de socavación. La profundidad de socavación alcanzada depende del tipo y
tamaño de las partículas que conforman el lecho y la magnitud y duración del evento
hidrológico.
La socavación depende de muchos factores, que pueden agruparse en dos grupos
principales: los factores geomorfológicos y los factores de transporte.
Además, la socavación puede clasificarse en dos tipos según su naturaleza: socavación
general y socavación local.
2.1.1 FACTORES QUE INFLUYEN EN LA SOCAVACIÓN
FACTORES GEOMORFOLÓGICOS
Los factores geomorfológicos hacen referencia a las características de la cuenca y el río
analizado. Dentro de las características de la cuenca se incluyen los factores climáticos y los
usos y tipos de suelos, elementos de primera importancia para determinar las tasas de
erosión y transporte en el sitio de interés.
Las características de mayor interés son la pendiente, la geometría de la sección transversal,
su forma en planta, las características del material del lecho y las condiciones de borde del
canal.
La caracterización del material del lecho incluye la distribución del sedimento, la gradación, la
estratificación de capas en el lecho y la presencia de materiales cohesivos.
Los sedimentos no cohesivos son aquellos cuyo movimiento depende solamente de las
propiedades de las partículas que los componen y de la posición relativa entre las
mismas. Dentro de este grupo se encuentran las arenas y gravas. Los sedimentos
cohesivos, por el contrario, son aquellos cuyo movimiento inicial depende de la compleja
interacción físico-química entre las partículas coloidales y los efectos de la presión de poros.
La socavación en lechos de material cohesivo es un fenómeno mucho más complejo y no
puede ser evaluado según las características del tamaño de la partícula.
Otra característica importante son los controles geológicos, pues determinan los límites
probables de erosión.

FACTORES DE TRANSPORTE
 Los factores de transporte se relacionan con el transporte tanto de agua como de
sedimentos. Las características del flujo, como la velocidad, duración, caudal y
frecuencia, así como las tasas de transporte y tipo de sedimento transportado durante tales
eventos, son necesarias para estimar la profundidad de socavación

2.1.2 TIPOS DE SOCAVACIÓN

SOCAVACIÓN GENERAL
La socavación general es aquella disminución en el nivel base del lecho del cauce como
consecuencia de aumento en la velocidad y el esfuerzo cortante del flujo en el lecho, que
pone en movimiento las partículas de fondo y de las márgenes que se encuentran en
equilibrio, indistintamente de la presencia o no de cualquier estructura antropogénica.

Los procesos fluviales y geomorfológicos que dan lugar a la socavación general ocurren
en un amplio rango de escalas espacio-temporales

Dentro de la socavación general se incluyen los procesos de degradación naturales como

incisión del lecho, migración lateral, cambios en las formas del lecho y socavación en
confluencias.

En términos de escala espacial, inestabilidades dentro de la cuenca pueden influir en los

procesos que contribuyen a la socavación general. Por ejemplo, cambios en los usos del
suelo pueden afectar la disponibilidad de sedimentos que transporta un río y por
consiguiente afectar directamente las tasas de socavación de algún un tramo en
particular.

Los cambios ocasionados por la socavación general a diferentes escalas espaciales son
proporcionales a la escala temporal en que ocurren, es decir, cambios en los patrones de
drenaje y perfiles longitudinales en una cuenca varían en escalas de tiempo geológicas.
A escala de tramo, el alineamiento de un cauce puede variar desde años a siglos,
cambios en la sección transversal pueden ocurrir de días a años, y cambios a escalas
locales, como las formas del lecho, pueden variar incluso a escala horaria.

Los cambios a largo plazo son consecuencia de una degradación progresiva y cuasi-
permanente debida a cambios hidrológicos y/o geomorfológicos a escala de cuenca, por
ejemplo, cambio climático, disminución del nivel base o construcción de una presa. A
corto plazo, la socavación general se desarrolla durante una creciente o crecientes
continuos durante un espacio temporal corto (días).
La socavación es un fenómeno complejo que para la comprensión del problema y su
solución requiere un análisis cualitativo inicial, enfocado en entender los procesos físicos
a escala de cuenca, de tramo y de sección.

La relación cualitativa propuesta por Lane (1955) ayuda a comprender las respuestas de
un río ante cambios en la cuenca que pueden producir agradación o degradación del
lecho. En la Figura 1-1 se representa el sistema fluvial como una balanza que permite
predecir la respuesta cualitativa de un cauce ante algún cambio que se presente en
alguno de los elementos característicos de la cuenca.

Figura 1 Analogía de la balanza de Lane

SOCAVACIÓN LOCAL

Contrario a la socavación general, la socavación local es el descenso abrupto en el lecho

fluvial, generado por la remoción del material de fondo debido a algún elemento físico
componente de alguna obra dentro del cauce o alguna singularidad natural que se
encuentre interpuesta dentro de la corriente. Dada la complejidad del movimiento del
fluido causante de la socavación localizada, su análisis se deriva en teorías diferentes y
en ocasiones contradictorias, dependiendo de las condiciones que se hayan tenido en
cuenta para su definición.

Este tipo de socavación ha sido más estudiado que la socavación general, pues los
análisis de socavación se realizan, en su mayoría, para conocer el diseño de
cimentaciones de obras que están en contacto directo con el cauce (puentes, espolones,
estribos, pilas, etc), ya que son obras de importancia en la subestimación de la
 socavación podría conllevar a la destrucción total de la estructura, o la sobreestimación
llevaría a adoptar profundidades de cimentación que hacen extremadamente costosa la
obra.

Sea cual sea el tipo de socavación estudiado, determinar la magnitud de la profundidad

de socavación es complicado por la naturaleza cíclica del proceso erosivo. La erosión
puede alcanzar su máximo valor durante una creciente, y volverse difícilmente visible tras
el paso de la misma al rellenarse con sedimento los “huecos” erosionados.

La socavación en lechos con fondo móvil se da cuando el río analizado transporta una
carga de sedimentos antes de entrar a la sección de control del tramo de interés, en
este caso el equilibrio en el balance de transporte determina la finalización del proceso
erosivo. Cuando el caudal sólido de ingreso al tramo de análisis es mínimo o nulo, se
denomina socavación en agua clara, y el proceso erosivo se detiene cuando el esfuerzo
cortante del cauce disminuye por debajo de los valores críticos de iniciación de
movimiento del material que conforma el lecho del tramo. También se considera erosión
en agua clara cuando el material del lecho es transportado en suspensión a través del
hueco de socavación hacia aguas abajo sin depósito en la sección de interés.

MODELAMIENTO HIDRÁULICO
MODELACIÓN EN HIDRÁULICA

Para los ingenieros de diseño es esencial predecir en forma exacta el comportamiento de las
estructuras hidráulicas bajo las condiciones de diseño (por ejemplo, la capacidad máxima de
descarga), las condiciones de operación y las situaciones de emergencia (por ejemplo, la
creciente máxima para un vertedero). Durante las etapas de diseño, los ingenieros necesitan
una “herramienta” de predicción confiable para comparar los comportamientos de diferentes
opciones de diseño.

En mecánica de fluidos e hidráulica, muchos problemas no pueden resolverse mediante la

teoría o la referencia a datos empíricos. En la práctica, un diseño apropiado de estructuras
hidráulicas se deduce utilizando un modelo bien definido como un sistema que opera en
forma similar a otros sistemas y el cual suministra una predicción acertada de las
características de estos.

En Hidráulica, se utilizan muy a menudo dos tipos de modelos: modelos computacionales y

modelos físicos. Los modelos numéricos son programas de computador que resuelven las
ecuaciones básicas de la mecánica de fluidos: por ejemplo, la integración numérica de la
ecuación de remanso para predecir el perfil longitudinal de la superficie libre de un flujo
gradualmente variado en un canal abierto. Sin embargo, su aplicación está restringida a
situaciones de flujo y condiciones de frontera simples para las cuales las ecuaciones básicas
puedan integrarse de forma numérica y tengan significado. La calibración y la validación de
modelos computacionales es en extremo difícil, y la mayoría de los modelos de computador
es aplicable solo en rango muy específico de condiciones de flujo y condiciones de frontera.
Con mayor frecuencia, se debe utilizar los modelos físicos, incluyendo casos para la
validación de modelos computacionales. (Chanson 2004).

2.2 MODELACIÓN DE LA SOCAVACIÓN

La estimación de la profundidad de socavación generalmente se basa en relaciones
empíricas o modelación computacional, que aproxima numéricamente la solución de las
ecuaciones de continuidad y transporte de sedimentos.

Dentro de las relaciones empíricas se destacan dos enfoques, el primero se basa en la teoría
del régimen y el segundo en el equilibrio entre la velocidad del flujo y la velocidad crítica para
el movimiento de las partículas que conforman el lecho.

2.2.1 CONDICIONES DE SIMILITUD

Para satisfacer las condiciones de similitud de todo Prototipo-Modelo, se requiere cumplir

con las siguientes condiciones (Gherardelli, 2007):
1) Similitud Geométrica, que se cumple cuando la relación de todas las dimensiones
de longitud homólogas (geométricamente equivalentes) son iguales; de donde
implica que solo la similar dad de forma se toma en cuenta.

la relación de longitud

Lp …………………………………………………… (01)
Lr=
Lm

Dónde: Lr: Relación de longitudes, que debe ser constante para todas
las partes del modelo y prototipo.

2) Similitud Cinemática, la relación de todos los puntos geométricamente homólogos,

velocidades y aceleraciones entre modelo y prototipo en cualquier espacio de
tiempo son constantes.

Lr ………………………………………………… (2)
vr=
Tr
vr Lr
ar = = ………………………………………………… (3)
Tr Tr 2
En el que la relación de tiempo, Tr, es:
1
Tp ………………………………………………… (4)
Tr= =Lr 2
Tm
 3) Similitud Dinámica, se debe satisfacer la similitud cinemática y considerar una
relación constante de la acción de las fuerzas sobre las partículas del fluido,
tales como fuerza de fricción, de gravedad, de Coriolis y de inercia. Lo anterior
implica que el polígono de fuerzas sobre los puntos homólogos debe ser
geométricamente similar.

Fr=Mr . ar ………………………………………………… (5)

La relación se puede escribir de la siguiente forma:
………………………………………….. (6)
Fr=ρr . Lr 2 .vr 2

4. Semejanza Sedimentológica, fundamentada básicamente las condiciones de flujo de las

aguas en el cauce y su interrelación con los sedimentos del lecho móvil, específicamente en la
investigación para cauces no cohesivos, no se considera efecto de vegetación u otros obstáculos
que genere resistencia al flujo. La similitud, se fundamenta en el equilibrio de una partícula en
el lecho, la cual debe ser semejante a la que será utilizada en laboratorio.

Yalin y Karahan (1979) citado en Martínez (2001), son quienes ensayan la zona baja del
diagrama de Shields, en el que se observa que en la curva de separación del movimiento baja
más el valor asintótico pasa a valores de 0.45 (Frs) alcanzándose dicho valor para R*>70. En
donde el número de Reynolds ya no tiene incidencia en el movimiento.

Para Fuentes, R (2002). El número de Reynolds pierde su influencia si sus valores son muy
grandes. En este caso, esto ocurre para un valor R*≈70 (Yalin, 1971), Por lo que se obviaría la
similitud de Reynolds; a este tipo de modelos se los conoce con el nombre de modelos
Eulerianos a gran número de Reynolds.

Segùn French (1988), los investigadores {Bogardi, 1959, o (Henderson, 1966) y Raudkivi,
1976} han aseverado que el requisito de que Fs y R* tengan exactamente el mismo valor en el
modelo y prototipo puede y debe relajarse un poco, cuando R*>=100 ya que el flujo alrededor
de los granos es totalmente turbulento por tanto no es necesario que se cumpla

Rm=R p ……………..……………………….. (6)

Maza y Alvarez (1968) indican que para que exista similitud entre modelo y prototipo, en
cuanto al número de Reynolds no se puede conseguir la similitud numérica, por lo que para la
similitud debería cumplirse que ambos puntos caigan en el mismo tipo de fondo para el
diagrama de Shields.

2.2.2 Teoría del régimen

La teoría del régimen fue iniciada en 1895 por Kennedy, quién para construir una red de
canales no revestidos para riego, observo y midió las dimensiones de 22 canales del sistema
Alto Bari Doab (India) que ya estaban en operación, y su sección se había ajustado a unas
dimensiones estables para los caudales, líquidos y sólidos, que transportaban. Kennedy
obtuvo que la velocidad media de flujo está en función de la profundidad de equilibrio de los
canales, y con esa relación dimensionó futuros canales (Maza Álvarez & García Flores,
1996).

……………………………………..(7)

Lindley (1919) empleó el término régimen para definir que un canal se encuentra en
régimen cuando su sección y pendiente están en equilibrio para el caudal transportado,
de tal manera que aumentos o disminuciones de él, hacen que el ancho y la profundidad
se modifiquen (Maza Flores & García Flores, 1996). Numerosos autores propusieron
relaciones empíricas de este tipo, basados en datos de canales de riego de la India,
Pakistan, Egipto y Estados Unidos.

En este trabajo se describen los métodos de Lacey (1930) y Blench (1939), el primero
por estar basado en la recopilación de muchos métodos anteriores a él, y el segundo por
ser uno de los más difundidos.

LACEY

El autor continuó y amplio los estudios de Lindley, quien ya había observado 4345 km de
canales en la India. En 1930 presentó las ecuaciones para obtener los principales
parámetros hidráulicos y geométricos de canales estables.

La ecuación de Lacey (y todas las formulaciones de la teoría del régimen) no es

dimensionalmente correcta, por tanto debe ingresarse en las unidades adecuadas. La
ecuación aquí mostrada se aplica utilizando el sistema métrico y el d50 en mm.
 ……………………………………..(8)

BLENCH

Basado en las observaciones de varios autores, Blench en 1939 y 1941, así como en
publicaciones posteriores, presentó sus fórmulas básicas y de diseño.

Para gravas (d50>2mm) la ecuación de Blench se puede escribir como

………………….……………..(9 )

Donde el d50 está expresado en mm y las demás variables en el sistema métrico

modelos más comerciales de éste tipo se destacan el HEC-6 (actual HEC-RAS), GSTAR
(actual SRH) y el MIKE11, desarrollados por el U.S. Army Corps of Engineers, el U.S.
Bureau of Reclamation y el Danish Hydraulics Institute, respectivamente. Los resultados
de éstos y de cualquier otro modelo unidimensional para predecir cambios en el lecho
depende en gran medida de las ecuaciones de transporte adoptadas en la modelación
(Melville & Coleman, 2000).

Asimismo, Fan (1994) concluye en sus estudios que los modelos unidimensionales son
altamente dependientes de los datos de entrada, incluso más que los modelos
bidimensionales, siendo esenciales para su desarrollo, calibración, validación e
implementación. El autor encontró que para los mismos datos, diferentes modelos
pueden predecir resultados significativamente diferentes (Melville & Coleman, 2000).

Para el presente trabajo, se seleccionaron los dos modelos unidimensionales más

comerciales en el medio y que han tenido un mayor y continuo desarrollo en modelación
de cambios en el lecho, el HEC-RAS (antiguo HEC-6) y el SRH-1D (antiguo GSTAR). A
continuación se presenta una descripción de sus principales características, capacidades
y limitaciones.