Science & Mathematics > Earth Sciences">
Tesis - Diana Violeta Díaz Llatas
Tesis - Diana Violeta Díaz Llatas
Tesis - Diana Violeta Díaz Llatas
FACULTAD DE INGENIERÍA
ESCUELA ACADÉMICO PROFESIONAL DE INGENIERÍA HIDRÁULICA
TESIS
Para optar el título profesional de:
INGENIERO HIDRÁULICO
Asesorada por:
Ph.D. OSWALDO ORTIZ VERA
CAJAMARCA – PERU
2018
AGRADECIMIENTOS
A las instituciones que apoyaron con diversas actividades como lo son la Universidad
Nacional de Cajamarca, E.P.S. SEDACAJ con su gerente general Ing. Oscar Santisteban
Kaneko, equipo de logística y operadores y el Instituto Geofísico del Perú (IGP); quienes
complementaron las labores necesarias para el proyecto.
A todo el equipo de trabajo que conforma el proyecto anteriormente mencionado, por su
apoyo en diversas etapas del proyecto de investigación.
A mi familia quienes has sido un gran apoyo desde el inicio hasta la actualidad, a mis padres
Gladis y Leoncio por su apoyo incondicional, por incentivarme siempre a lograr mis metas
conservando la humildad.
ii
DEDICATORIA
La presente tesis va dedicada primeramente a Dios, porque es quien día a día guía mi
camino y me permite enfrentar nuevos retos. A mi madre Gladis y mi padre Leoncio, por
su cariño y apoyo incondicional en el logro de mis metas. A mis hermanos Javier, Henry,
Werner, Marco, Lenin y Roxana quienes con su ejemplo me siguen enseñando que se
pueden lograr grandes cosas.
iii
ÍNDICE
I. INTRODUCCIÓN ............................................................................................................. 1
v
ÍNDICE DE TABLAS
vi
ÍNDICE DE FIGURAS
vii
Figura 27. Porcentaje de material retenido (%) sub – superficial para cada punto de
muestreo, respecto a su diámetro (mm). ............................................................................ 66
Figura 28. Cuatro puntos de muestreo seleccionados en el lecho superficial del río
Ronquillo............................................................................................................................ 69
Figura 29. Curvas granulométricas de lecho superficial de río Ronquillo. ....................... 70
Figura 30. Material retenido (%) en cada punto de muestreo del material lecho superficial
............................................................................................................................................ 71
Figura 31. Concentraciones de sedimentos en suspensión (SSC) y caudales líquidos (Q)
medidos en la estación Ronquillo. ..................................................................................... 76
Figura 32. Curva hipsométrica y de frecuencia de altitudes de la cuenca Ronquillo. ...... 78
Figura 33. Eventos de caudales extremos que llenaron el barraje. ................................... 79
Figura 34. Serie temporal de concentraciones de sedimentos Vs. Caudales líquidos en la
estación Ronquillo. ............................................................................................................ 81
Figura 35. Correlación entre concentraciones de sedimentos en suspensión (SSC) y
caudales líquidos (Q) medidos en la estación Ronquillo. .................................................. 82
Figura 36. Volumen de sedimento acumulado (m3) para cada evento de crecida utilizando
las fórmulas empíricas. ...................................................................................................... 83
Figura 37. Volúmenes acumulados de sedimentos de fondo Estimados y Observados. ... 85
Figura 38. Coeficiente de proporcionalidad entre volumen estimado (Ve) y volumen
medido en campo (Vm). .................................................................................................... 86
Figura 39. (a) Levantamiento topográfico del cauce de río Ronquillo y (b) batimetría del
vaso formado por el barraje de la captación Ronquillo. .................................................... 93
Figura 40. (a) Material del lecho superficial y (b) material del lecho sub-superficial del río
Ronquillo............................................................................................................................ 93
Figura 41. (a) Tamizado de material sub-superficial, (b) Material para densidad de gruesos,
(c) y (d) Ensayo para determinar densidad de material fino. ............................................. 94
Figura 42. (a) Filtrado d muestras de sedimentos en suspensión, (b) Pesado de muestras
filtros secos para obtensión de concentración de sedimentos en suspensión..................... 95
viii
RESUMEN
Palabras claves:
Sedimento de fondo.
Sedimento en suspensión.
Río Ronquillo.
Volumen de sedimento de fondo.
Coeficiente de proporcionalidad.
ix
ABSTRACT
Regarding the bedload sediment, a sediment monitoring strategy was established, using the
dam formed by the weir of the catchment structure as an accumulator or totalizer of
sediments of known storage capacity, the accumulator filling periods, flows and time
involved were recorded. With this strategy, 4 volumes of bedload sediments measured
during 2016 were recorded.
Four empirical formulas were selected: Parker, Klingteman and McLean 1982 (P-K-Mc),
Smart and Jaeggi 1983 (S-J), Bathurst (1987) and Meyer Peter and Muller 1948 (M-P-M),
the results of each of them were compared with the volume of sediment observed in the
field, all formulas used over-estimated the volume of actual bottom sediment, the one that
in smaller amount over-estimates the formula of Parker, Klingteman and McLean 1982 (P-
K-Mc).
Keywords:
Bedload Sediment.
Suspended Sediment.
Río Ronquillo.
Bedload Sediment Volume.
Proportionality coefficient.
x
I. INTRODUCCIÓN
Los ríos de montaña están típicamente caracterizados como ríos que tienen gradientes
longitudinales que oscilan entre 0.1% y 20% (PAPANICOLAOU, A. et al., 2004), sus
lechos están compuestos de sedimentos gruesos y tienden a desarrollar procesos de
acorazamiento o blindaje, los cuales ocasionan una alta variabilidad espacial y temporal en
el arrastre de fondo (GRAY, J. et al., 2010).
La medición de sedimentos es una tarea compleja, ya que realizar la elección del sistema
de medición más preciso es a menudo difícil porque depende de muchos parámetros
(HABERSACK, H. et al., 2016), estos parámetros están relacionados con la gran
variabilidad espacial y temporal del proceso de arrastre de fondo, limitaciones de los
muestreadores, técnicas de muestreo y recursos para recolectar datos precisos y
cuantificables de carga de lecho (GRAY, J. et al., 2010).
1
1.1. PLANTEAMIENTO DEL PROBLEMA
Las investigaciones realizadas hasta el momento muestran que los modelos hidráulicos, en
este caso de transporte de sedimentos desarrollados sobre la base de la evidencia empírica
en ríos de arena de suave pendiente no son aplicables a ríos de montaña de sedimento
grueso, en consecuencia su aplicación en ríos de grava con lecho acorazado plantea diversas
incertidumbres (en cuanto al diámetro característico del sedimento más adecuado para
sustituir en las fórmulas, así como en lo concerniente al umbral de rotura o alteración de la
coraza, y a las diversas fases de transporte).
2
1.3. OBJETIVOS
Objetivo general
Objetivos específicos
3
II. MARCO TEÓRICO
4
la fórmula asume que el incremento de las tasas de transporte es debido al
movimiento promedio de las partículas, relacionado directamente a la velocidad del
flujo líquido. Lo cual contrasta con lo observado durante el periodo de monitoreo.
A demás, las constantes empíricas de la fórmula se desarrollaron en canales de
aforo, conteniendo diversos tamaños de partículas de sedimento (0.78 a 28.6 mm)
similares al tipo de granulometría obtenida en ambas cuencas. Sin embargo, este
resultado todavía posee gran incertidumbre durante ocurrencia de grandes avenidas
y la influencia de eventos extremos como El Niño.
5
rendimiento relativamente bueno del enfoque del esfuerzo de Yager (2006) para las
corrientes con mayores concentraciones de roca indica que esta corrección basada
en la física para la resistencia al flujo adicional es un paso adelante para una mejor
caracterización de tales condiciones de corriente desde un punto de vista teórico.
6
2.2. BASES TEÓRICAS
Un sistema fluvial se puede dividir en tres zonas, cada una de las cuales es un sistema de
proceso de respuesta con sus propias entradas y salidas de sedimentos (Figura 1). Según
(CHARLTON, R. 2007) estas zonas son:
7
- Área (A):
Es la relación entre el área “A” de la cuenca y el cuadrado del máximo recorrido “L”. Este
parámetro mide la tendencia de la cuenca hacia las crecidas, rápidas y muy intensas a lentas
y sostenidas, según que su factor de forma tienda hacia valores extremos grandes o
pequeños respectivamente (VILLÓN, M. 2002).
𝐴 (1)
𝐹=
𝐿2
Donde:
A = Área de la cuenca.
L = longitud del máximo recorrido.
8
Las cuencas redondeadas tienen tiempo de concentración cortos con gastos pico muy
fuertes y recesiones rápidas, mientras que las alargadas tienen gastos pico más atenuados y
recesiones más prolongadas.
𝐾𝑐 = 0.28𝑃𝐴−0.5 (2)
Donde:
P = Perímetro de la cuenca.
A = Área de la cuenca.
Las formas de la cuenca, en concordancia con los valores que adopte los índices de
compacidad, de muestra en la Tabla 1.
Tabla 1
Formas de la Cuenca de acuerdo al Índice de Compacidad Kc.
9
𝑛
1 (3)
𝐻 = ∑ 𝐻𝑖 ∗ 𝐴𝑖
𝐴
𝑖=1
Donde:
𝑛
𝐶 (4)
𝑆𝑐 = ∑ 𝑙𝑖
𝐴
𝑖=1
Donde:
Sc = Pendiente de la cuenca.
C = Equidistancia entre curvas de nivel.
A = Áreas de la cuenca.
Li = Longitud de cada curva de nivel.
La clasificación de las cuencas según la pendiente de laderas, se aprecia en la Tabla 2.
10
Tabla 2
Clasificación de las cuencas de acuerdo a la pendiente promedio de laderas.
0–3 Plano P1
3–7 Suave P2
7 – 12 Mediano P3
20 – 35 Accidentado P4
35 – 50 Fuerte P5
50 – 75 Muy Fuerte P6
50 – 75 Escarpado P7
75 Muy escarpado P8
Este es un parámetro importante, ya que permite ver el comportamiento del recurso hídrico,
como por ejemplo, para la determinación de las características óptimas de su
aprovechamiento hidroeléctrico, o en la solución de problemas de inundaciones (VILLÓN,
M. 2002).
∑𝑛𝑖=1 𝑙𝑖
𝑆𝑜 = (5)
𝑙
∑𝑛 √ 𝑖
[ 𝑖=1 𝑆𝑖 ]
Donde:
𝑆𝑜 = Pendiente del cauce principal.
n = Número de tramos de similar pendiente.
𝑙𝑖 = Longitud del tramo i.
𝑆𝑖 = Pendiente del tramo i.
11
- Longitud del máximo recorrido (L)
Es la medida de la mayor trayectoria de las partículas del flujo comprendida entre el punto
más bajo del colector común conocido como punto emisor, y el punto más alto o inicio del
recorrido sobre la línea de divortio aquarum. Este parámetro tiene relación directa con el
tiempo de concentración de la cuenca, el mismo que depende de la geometría de la cuenca,
de la pendiente del recorrido y de la cobertura vegetal (ORTIZ, O. 2004).
Este parámetro, llamado también tiempo de equilibrio, es el tiempo que toma la partícula,
hidráulicamente más lejana, en viajar hasta en punto emisor. Para ello se supone que el
tiempo de duración de la lluvia es de por lo menos igual al tiempo de concentración y que
se distribuye uniformemente en toda la cuenca. Este parámetro tiene estrecha relación con
el gasto pico y con el tiempo de recesión de la cuenca, tiempos de concentración muy cortos
tienen gastos pico intensos y recesiones muy rápidas, en cambio los tiempos de
concentración más largos determinan gastos pico más atenuados y recesiones mucho más
sostenidas (VILLÓN, M. 2002).
0.38
𝐿 ∗ 𝐿𝐶
𝑇𝑐 = 𝐶 ( ) (6)
√𝑆
Donde:
𝑇𝑐 = Tiempo de concentración, en horas.
L= Longitud de máximo recorrido, en Km.
𝐿𝐶 = Longitud de centroide, en Km.
S = Pendiente media del máximo recorrido.
C = Coeficiente de escorrentía.
12
- Coeficiente Orográfico (Co):
𝐻2
𝐶𝑜 =
𝐴 (7)
Donde:
𝐶𝑜 = Coeficiente Orográfico, adimensional.
H = Altitud media del Relieve.
A = Área de la cuenca.
b. CUENCA RONQUILLO
- Red hidrográfica
La cuenca Ronquillo está conformada por los ríos Balconcillo, Manzana, Cushunga, Tres
Ríos, a continuación se describen con mayor detalle.
La unión de los ríos Balconcillo y Cushunga y Manzana forman el río Tres Ríos, el cual
metros más abajo, toma el nombre de Río Ronquillo, debido a que adopta el nombre de la
zona (Ronquillo).
El río Balconcillo se forma después de recibir los aportes de la quebrada Conga. El río
Cushunga se origina de la unión de las quebradas Portachuelo y Poroporo, y corriente abajo
recibe los aportes de Garbanzo Pampa y Sacchayoc (BENAVIDES, I. et al., 2007).
13
El río Manzana se forma de la confluencia de las quebradas Choro Mayo y Carhuacasha
Corral Blanco (BENAVIDES, I. et al., 2007).
El río Ronquillo es la fuente que abastece del recurso hídrico a la captación Ronquillo, el
agua captada se conduce hasta la Planta de Tratamiento de Agua Potable Santa Apolonia,
la cual abastece al 30% de la población Cajamarquina.
14
- Demografía en cuenca ronquillo.
Dentro de la cuenca Ronquillo se ubican 10 centros poblados los cuales son: Carhuaquero,
Coñorcucho, Chamis, Corisorgona, Yuraccirca, Cushunga, Merinillo, Capillania,
Secsemayo (lote I y II) y Carhuaconga. Figura 3.
15
En la cuenca Ronquillo se encuentran las siguientes formaciones geológicas que en
(VÁSQUEZ, G. & RODRÍGUEZ, M. 2011) las describen de la siguiente manera.
Esta formación está formada por areniscas blancas de grano medio a grueso, tiene un grosor
promedio de 500 m. en algunos lugares se observa estratificación cruzada y marcas de
oleaje.
Esta formación consiste en una secuencia fosilífera de calizas arenosas, lutitas calcáreas y
margas, las que por intemperismo adquieren un color crema amarillento. Su aspecto terroso
amarillento es una característica para distinguirla en el campo. Generalmente los bancos de
margas se presentan muy nodulosos y las calizas frescas muestran colores gris parduzcos
algo azulados. Su grosor varía de 200 a 250 m.
Esta formación yace concordantemente sobre la formación Chúlec e infrayace con suave
discordancia a la formación Yumagual, relación observable en la cuenca de Pulluicana, en
la carretera Cajamarca- La Encañada, al este de los Baños del Inca. Edad y correlación.- la
formación Pariatambo contiene restos de moluscos, estas especies son típicamente
16
pelágicas del Albiano medio. Se correlaciona con la parte superior de la formación
Crisnejas y con la formación Yacu Ushco. Se extiende ampliamente en los andes centrales
del Perú.
Consiste en una secuencia de margas y calizas gris parduzcas en bancos más o menos
uniformes, destacando un miembro medio lutáceo margoso, amarillento, dentro de un
conjunto homogéneo presenta escarpas debido a su dureza uniforme. En algunos horizontes
se observan nodulaciones calcáreas. Tiene un grosor aproximado de 700 m.
Esta formación consiste de calizas gris oscuras o azuladas, con delgados lechos de lutitas y
margas. Las calizas se presentan en bancos gruesos con escasos fósiles. Esta formación
yace concordantemente sobre la formación Quilquiñán y con la misma relación infrayace
a la formación Celendín. Su grosor varía entre los 600 y 700 m.
Esta unidad consiste en gruesos estratos de rocas volcánicas, intercaladas en la base con
areniscas rojizas y en la parte superior de una espesa secuencia de aglomerados y
piroclásticos bien estratificados. Alcanza un espesor de 900 m.
El volcánico san pablo yace con suave discordancia erosional al volcánico Chilete e
infrayace al volcánico Huambos en igual relación.
En casi todo el sector andino del departamento y en ambos lados de la divisoria continental
se encuentra aflorando una unidad de tobas ácidas, en los alrededores del pueblo del mismo
17
nombre. La topografía del volcánico Huambos es bastante característica, generalmente
forma llanuras delimitadas por farallones o escarpas donde se aprecia la estratificación,
sobre estas superficies se ha desarrollado un sistema dendrítico de drenaje
Respecto a la litología y grosor, el volcánico Huambos está compuesto por tobas y brechas
de composición ácida. La litología común muestra fragmentos de cuarzo hasta de 3 mm.
de diámetro y cristales de biotita en una matriz feldespática que probablemente corresponde
a una toba dacítica. Se encuentran brechas dacíticas compuestas por bloques grandes de
toba envueltos por una matriz tobácea, esta litología se observa en el cuadrángulo de Chota.
En la mayoría de los casos los piroclásticos del volcánico Huambos están bien estratificados
en capas medianas a gruesas, parcialmente compactadas con escasos niveles de toba
soldada o ignimbrita. El volcánico Huambos tiene un grosor promedio de 50 a 100 m. y un
máximo de aproximadamente 400 en la localidad de Huambos.
Fuente: INGEMET.
18
- Actividades socioeconómicas.
Ganadería
Los pobladores de la cuenca en estudio, en mayoría crían ganado ovino del cual provechan
su lana para sus tejidos y su carne para su alimentación, en menor cantidad podemos
encontrar ganado vacuno y avícola.
En 60 % del total de las familias poseen sus yuntas para la preparación de los terrenos, las
cuales venden posteriormente para solventar los gastos de alimentación y educación de sus
hijos u otras necesidades. La ganadería no es una actividad de mucha importancia por falta
de pasturas (CHUQUIRUNA, M. et al., 2008).
Agricultura
- Geomorfología
El relieve característico de la zona del sector Urubamba ‐ Ronquillo, está vinculado a los
procesos tectónicos y a las modificaciones bioclimática que se han generado desde el inicio
de su aparición. Las condicionantes morfológicas como la inestabilidad, vulnerabilidad y
riesgo, siempre han condicionado el uso y ocupación del territorio. En términos generales,
la geomorfología es muy favorable a la concentración de flujos asociados a sólidos, por lo
que es frecuente los procesos de geodinámica externa, esto es ayudado por la presencia de
pendientes empinadas, presencia del canal de escurrimiento principal, la frecuencia de
litologías poco consistentes y una cierta proporción de terrenos sometidos a condiciones
semiáridas, climática y ausencia de cobertura vegetal, estas condiciones aumentan las
acciones erosivas de las cuencas y la gran presencia de fenómenos de remoción de masas
de suelos y de rocas.
19
Como parte del relieve andino de Cajamarca es su elevada tasa de erosión actual, con
predominio de acciones erosivas ligadas al escurrimiento difuso y concentrado en surcos y
cárcavas, motivadas tanto por causas naturales como por labores agrícolas en terrenos
altamente susceptibles a la erosión. Según lo apreciado, los Peligros Naturales y Antrópico
que inciden sobre el núcleo urbano de zona del sector Urubamba – Ronquillo y su entorno
inmediato son de origen Geológico (Sismicidad), origen climático (Inundaciones generadas
por la acción pluvial) y de origen geológicos‐climáticos (Deslizamientos) (MPC, 2012).
Fuente: (CRUZADO)
La morfología dominante son las terrazas en la parte central. Más hacia el frente la
pendiente es más fuerte y se ubican escarpes viejos (parcialmente cubiertos por
vegetación) de hundimientos.
20
c. CONCEPTOS GENERALES
La mayoría de sedimentos suministrados a los cauces del río varía mucho en tamaño, desde
las partículas microscópicas de arcilla a grandes rocas. Como es de conocimiento, existen
diferentes mecanismos implicados en el transporte de sedimentos gruesos y finos dentro de
los cauces fluviales. Estos se ilustran en la Figura 6 y tienen importantes implicaciones para
la forma en que los sedimentos de diferentes tamaños se transfieren a través del sistema.
(CHARLTON, R. 2007).
El material grueso, típicamente arenas gruesas, gravas finas y las partículas más grandes,
se mueve a lo largo del lecho del canal como carga de fondo. Las partículas están en
contacto continuo o regular al lecho del canal y se mueven por rodamiento, deslizamiento
o en una serie de movimientos de salto llamados saltación. (CHARLTON, R. 2007).
21
La acción del agua sobre el fondo puede caracterizarse por una tensión cortante en el fondo
“”. La resistencia de la partícula a ser movida puede relacionarse con su peso sumergido,
el cual es función de peso específico sumergido (s-), y del diámetro “D” que caracteriza
el volumen de la partícula. Con estas tres variables puede formarse el parámetro
adimensional de Shield o tensión córtate 𝜏̅ (VIDE, J. 2007):
𝜏𝑜
𝜏̅ =
(𝛾𝑠 − )𝐷 (8)
Esta relación compara la fuerza que tiende a producir el movimiento (acción de arrastre
proporcional a 𝜏𝑜 𝐷2 ) con la fuerza que procura estabilizarlo o mantenerlo en reposo (acción
del peso proporcional a (𝛾𝑠 − )𝐷3 .
Para un cauce prismático, el peso del líquido “W” es igual al producto del peso específico
por volumen de control W= γ.A.dx donde “A” es el área transversal del cauce. La
componente en la dirección del flujo Wx = g.A.dx.sen θ. Si el ángulo “θ” es pequeño senθ
≅ tgθ ≅S siendo “S” la pendiente del tramo de cauce analizado (OLAZÁBAL, J. 2015).
22
Por otro lado, existe una fuerza de fricción aplicada en todo el contorno del cauce que se
opone al escurrimiento τo.P.dx, donde “P” es el perímetro mojado del cauce. Analizando
el equilibrio de ambas fuerzas τo.P.dx = γ.Α.S.dx resulta la Ec. 16 (OLAZÁBAL, J. 2015).
(9)
o = RhS
Donde:
: Tensión en el fondo (Kgf.m-3)
R: Radio hidráulico (m)
S: Pendiente media de la cuenca. -
𝛾𝑠 : Peso específico del suelo (Kgf.m-3)
𝛾: Peso específico del agua (Kgf.m-3)
Tamaño de partículas
Las partículas que conforman los suelos varían en un amplio rango de tipos como son:
grava, arena, limo o arcilla, dependiendo del tamaño predominante de partículas.
Varias organizaciones han desarrollado límites de las partículas según su tamaño los cuales
podemos observar en la Tabla 3.
23
Tabla 3. Clasificación de Partículas Según su Tamaño en milímetros (mm).
Grava
Arena Limo Arcilla
Nombre de (Haddadchi,
(Haddadchi et (Haddadchi et (Haddadchi
organización Omid, &
al.) al.) et al.)
Dehghani)
Instituto tecnológico de
Massachusetts (Smith & 2 mm 2 a 0.06 mm 0.06 a 0.002 mm 0.002 mm
Owens)
Departamento de
agricultura de Estados 2 mm 2 a 0.05 mm 0.05 a 0.002 mm 0.002 mm
Unidos (USDA)
Asociación Americana
de Funcionarios del
0.075 a 0.002
76.2 a 2 mm 2 a 0.075 mm 0.002 mm
Transporte y Carreteras mm
Estatales (AASHTO)
Sistema unificado de
clasificación de suelos
(U.S. Army Corps of
Engineers; U.S. Burneau 4.75 a 0.075 Finos (es decir, limos y arcillas
76.2 a 4.75 mm
mm 0.075 mm )
of Reclamation;
American Society for
Testing and Materials)
24
- Distribución granulométrica de las partículas.
Primero el suelo se seca en el horno, y luego todos los grumos se disgregan en partículas
pequeñas antes de ser pasados por la malla. Una vez que pasan las mallas, se determina la
masa del suelo retenida en cada malla. Cuando se analizan suelos cohesivos, resulta difícil
disgregar los grumos en partículas individuales. En tal caso, el suelo se mezcla con agua
para formar una lechada que luego se lava a través de las mallas. Las porciones retenidas
en cada malla se recolectan por separado y se secan en horno antes de que la masa retenida
en cada malla sea determinada.
Una vez terminado el análisis por cribado, sus resultados se presentan generalmente en
graficas semi-logarítmicas como curvas de distribución granulométrica (o de tamaño de
grano). Los diámetros de las partículas se grafican en escala logarítmica y el porcentaje
correspondiente de finos en escala aritmética. Por ejemplo, las curvas de distribución de
suelos se muestran en la Figura 8.
25
- Diámetros característicos de las partículas
D35: Diámetro de la malla por donde pasan el 35 por ciento de los sedimentos de la muestra.
Es el diámetro característico propuesto por Einstein para representar el diámetro de la
muestra.
D50: Diámetro que en muchos casos representa el diámetro medio. Shields utilizo este
diámetro para su estudio del inicio del movimiento.
∆1 𝐷1 + ∆2 𝐷2 + ∆3 𝐷3 +. . +∆𝑛 𝐷𝑛 𝐴𝑖 𝐷𝑖
𝐷𝑚 = =∑ (10)
100 100
Dónde:
Δi: Representa una porción del porcentaje del gráfico de distribución granulométrica
Di: Valor medio del diámetro que corresponde a Δi.
La densidad s es la masa que se tiene por unidad de tiempo. Existen tablas con valores
aproximados de los sólidos según su composición mineralógica, pero los valores reales se
determinan directamente en laboratorio. (DAS, B. 2001).
26
𝑀
𝜌𝑠 =
𝑉 (11)
Donde:
𝑃
𝛾𝑠 = (12)
𝑉
𝛾𝑠 = 𝜌𝑠 ∗ 𝑔
(13)
Donde:
𝛾𝑠 : Peso específico de la partícula. (Kg.s-m2) o (N.m-2)
𝑃: Peso de la partícula. (N)
𝑔: Aceleración de la gravedad. (m2.s-1)
- Entender los factores que afectan al valor de “n” con el fin de adquirir conocimientos
básicos del problema y disminuir el rango de incertidumbres.
- Consultar una tabla de valores comunes de n para diferentes tipos de canales o cauces.
- Examinar y familiarizarse con la apariencia de canales o cauces comunes cuyos
coeficientes de manning se conocen.
27
Método de Cowan:
Donde:
𝑛𝑜 : Es un valor básico de “n” para un canal recto, uniforme y liso en los materiales naturales
involucrados.
28
Tabla 4. Valores utilizados por Cowan para la estimación del “n” de manning.
29
Método de Strickler:
Cuando los lechos de los cauces naturales están constituidos por material pedregoso, donde
el sedimento es representado por un diámetro medio, se recomienda el uso de la ecuación
de Strickler para la estimación de no.
𝑛𝑜 = 0.038𝐷1/6 (15)
Donde:
D: Diámetro representativo de la rugosidad superficial (m)
El diámetro “D” es equivalente al diámetro D65, D90 o D95 dependiendo del acorazamiento
del lecho. Particularmente cuando los sedimentos ofrecen una granulometría gruesa y
extendida, el diámetro medio de la coraza es cercano al D90 o D95 obtenido de la curva
granulométrica original del lecho. (CHOW, V.T. 1983)
Existen muchos conceptos para definir el caudal dominante que a continuación veremos.
30
Figura 9. Vista transversal de un río que muestra el nivel de desborde del mismo.
Muchos autores han propuesto diversas fórmulas para poder estimar este caudal entre los
cuales, siendo La que mayor se ajusta a las características del río Ronquillo es la de
Williams 1978, esta fórmula se encuentra dentro del libro de (ALONSO, R. 2004).
Williams propuso una ecuación para estimar la capacidad máxima del cauce utilizando 233
datos de caudal y realizando una ecuación de regresión múltiple obtuvo la Ec. 16.
Los parámetros de aplicación que menciona el autor corresponden a 0.5 ≤ Q (m3/s) ≤ 28320,
0.000041 ≤ S ≤ 0.0810 y 0.00019 ≤ d50 (m) ≤ 0.19.
31
d. SEDIMENTOS DE FONDO Y MÉTODOS DE MEDICIÓN.
El gasto sólido que lleva un río en un momento determinado puede ser menor que su
capacidad de transporte. Así, un río que escurre sobre un fondo constituido por piedras de
gran tamaño puede tener como consecuencia de su velocidad una gran capacidad de
transporte, pero en realidad pudiera no haber transporte sólido (debido al gran tamaño de
las piedras que constituyen el lecho). En un caso como éste se dice que el río se encuentra
en estado de erosión latente. Este caso es frecuente en ciertos torrentes. Cuando un río
transporta su máxima capacidad de sólidos se dice que se halla en estado aluvional o de
saturación (ROCHA, A. 1998).
- Medición directa:
32
Figura 10. Trampa de Arrastre de Fondo
- Muestreadores:
El cálculo del arrastre de fondo se puede efectuar a partir de muestras recogidas por un
dispositivo que está situado por debajo del lecho de la corriente durante un tiempo
determinado y que luego son extraídas para pesarlas. Se han utilizado numerosos
dispositivos y su variedad demuestra la dificultad que existe para tomar una muestra exacta
y representativa. Los problemas que plantean los muestreadores del arrastre de fondo son:
La forma más sencilla es una cesta de alambre con una aleta estabilizadora como se ilustra
en la Figura 11. El material que toma este modelo es reducido porque interfiere con la
corriente y cierta cantidad de materia cada vez mayor a medida que la cesta se llena, se
desvía en torno al muestreador.
33
Figura 11. Muestreador Tipo Canasta
34
e. SEDIMENTO EN SUSPENSIÓN Y MÉTODOS DE MEDICIÓN.
El material fino en los ríos de transporta en dos formas, como material sólido y en solución
, el límite para diferenciarlos es el diámetro de sus partículas, siendo la partículas con
diámetros menores a 0,45 µm (micras) para el material en solución (OWENS, P. et al.,
2005).
- Muestreo manual
Para poder realizar este tipo de muestreo se tienen que tener las siguientes consideraciones:
Dentro de este tipo de muestreo se presentan diversas formas de realizarlo como pueden
ser:
35
Muestreador de profundidad integrado.
Este muestreador consta de un recipiente, una boquilla, un orificio para permitir la salida
de aire o una válvula de control, por el que una única muestra de material compuesto de
descarga, se recoge moviendo el dispositivo de muestreo verticalmente a través de la
corriente (VANONI, V. 2006).
Punto único
Si bien es la mejor práctica para tomar muestras a través de la sección transversal del flujo,
se reconoce que esto no siempre es posible, especialmente cuando el río está crecido.
36
- Muestreo automático
Pasivo
Figura 14. Muestreador pasivo que toma muestras de agua mediante el sistema de sifonaje.
37
Bomba automatizada
38
En este muestreador la velocidad de flujo dentro de la cámara principal de la toma de
muestras se reduce por un factor de alrededor de 600 en comparación con el flujo de
ambiente (PHILLIPS, J. et al., 2000).
Para poder estimar el transporte de fondo se presentan las siguientes fórmulas empíricas:
Las fórmulas que mayormente se utilizadas para la estimación del gasto sólido de fondo,
son las siguientes:
Meyer-Peter (1883-1969) fue director del antiguo Laboratorio de Hidráulica del Politécnico
de Zurich. En este laboratorio se han realizado durante muchos años importantes
investigaciones y estudios en modelo hidráulico sobre diversos problemas de transporte
sólido en ríos de fuerte pendiente, entre ellos, varios para el Perú.
Es una ecuación empírica obtenida en suiza a partir de ensayos de laboratorio y muy usada
en Europa. Cubre hasta pendientes de 2% y hasta tamaños de material de 30 mm (VIDE, J.
2007) y tirante de 0.01 a 1.20 m. Los experimentos se condujeron con materiales naturales
y sintéticos cuyo diámetro osciló entre 0.4 y 30 mm.
𝛾𝑠 3 0.5
𝑛𝑠 1.5
(17)
𝑔𝐵 = 8𝛾𝑠 ((𝑔( 𝐷𝑚 )) [( )1.5 𝜏∗ − 0.047]
𝛾𝑠 − 𝛾 𝑛
1/6
𝐷
𝑛𝑠 = 90 (18)
26
Donde:
39
𝐷𝑚 : Diámetro medio (m)
𝜏∗ : Número adimensional de Shields.
𝑛𝑠 : Coeficiente de manning debido a la partícula.
𝑛: Coeficiente de manning.
g: Aceleración debido a la gravedad (m2/s)
D90: Diametro de la partícula (D=D90) de sedimento superficial.
La fórmula de (PARKER, G. et al., 1982) difiere del resto de fórmulas en que se basa única
y exclusivamente en datos de campo. Estas fórmulas son de las más completas y fiables
(GARCÍA,C. & SALA, M. 1998).
1.- Todas los tamaños disponibles de las partículas, gruesas y finas, son aproximadamente
igual de móviles en el movimiento inicial y durante el transporte.
2. Las distribuciones de tamaño de partículas de la carga de fondo son los mismos que se
encuentran en el material del lecho y no cambian con la etapa de flujo.
3. Una superficie de material gravoso está presente durante flujo alto, así como durante
flujo bajo.
4. La cantidad de transporte de grava se predicen por el tamaño medio de las partículas del
sedimento sub-superficial.
5. Las partículas de varios tamaños son arrastradas a aproximadamente la misma tensión
de cizallamiento.
6. La distribución del tamaño no es importante.
7. Este modelo puede ser adecuado para pequeñas y medianos canales que tienen tamaños
de partícula entre 0,6 a 102 mm, con pendientes medias, y que no están dominados por la
carga de arena
40
Esta ecuación presenta las siguientes 3 componentes:
= 0.0025∅14.2
50 0.95 < ∅50
Donde ∅50 es Parámetro adimensional en base a tensión de corte necesario para movilizar
∗
partículas de tamaño D50 y 𝜏50 es la tensión de referencia de Shield que toma un valor de
0.0876.
∗
𝜏50
∅50 = ∗
𝜏𝑟50 (20)
∗
𝑢∗2
𝜏50 = (21)
(𝑠 − 1) ∙ 𝑔 ∙ 𝐷50𝑠𝑢𝑏
Todos los parámetros se relacionan para calcular la tasa de transporte de sedimento por
medio de la ecuación:
𝑊 ∗ ∙ 𝐵 ∙ 𝑢∗3 ∙ 𝜌𝑠
𝑄𝑏 = (22)
(𝑠 − 1)𝑔
También presentada por Smart (1984), basada en experimentos en canales con 77 datos
propios del laboratorio y 137 de Meyer Peter y Müller (1948) (LÓPEZ, R. 2003).
41
esfuerzo cortante adimensional sobre el valor crítico, las relaciones propuestas por Smart y
Jaeggi corresponden a una versión más completa de la ecuación de transporte de Meyer
Peter y Müller (AGUIRRE, J. et al., 2000), válidos para diámetros mayores a 0.4 mm y
pendientes entre 0.4% y 20% (SMART, G. 1984), en la forma:
∗
𝐷90 0.2 0.6 ∗ ∗0.5 ∗ ∗
Φ = 4 [ ] 𝑆 𝐶 𝜏 [𝜏 − 𝜏𝑐𝑗 ] (23)
𝐷30
Donde:
D90 y D30 representan los diámetros del 90% y 30% del material respectivamente y *cj
es el esfuerzo cortante crítico dado por:
∗
𝑆
𝜏𝑐𝑗 = 𝜏𝑐∗ cos(arctan 𝑆) [1 − ]
𝑡𝑎𝑛 𝜙 (24)
- Bathurst (1987).
∗
2.5 𝑆 3/2
Φ = (𝑞 − 𝑞𝑐 ) (25)
(Δ + 1)𝐷50 √(𝑔Δ𝐷50 )
De la cual D50 es el diámetro del 50% del material sólido, q es el caudal de agua por
unidad de ancho y qc es el caudal crítico y se obtiene de la siguiente expresión:
42
𝑞𝑐
= 0.15 𝑆 −1.12
3
√𝑔𝐷50 (26)
Para D16, este deberá ser tomado de la capa superficial del lecho, y es de la forma:
𝑞𝑐
= 0.21 𝑆 −1.12
3
√𝑔𝐷16 (27)
𝑞𝑠
Φ∗ = (28)
𝐷√𝑔∆𝐷
𝑡
(29)
𝑉𝑠 = ∫ 𝑄𝑠(𝑡) ∗ 𝑑(𝑡)
0
𝑞𝑠 (𝑡) ∗ 𝐿 (30)
Qs(t) =
𝜌𝑠
Donde:
𝑉𝑠 : Volumen de sedimento en un intervalo de tiempo (L3).
𝑄𝑠: Volumen de sedimento por unidad de tiempo (L3T-1).
𝑞𝑠 (𝑡): Gasto de sedimento por unidad de longitud (M L-1T-1)
43
L: Ancho de cauce (L).
𝜌𝑠 : Densidad de material depositado (M L3)
𝑉𝑠 (31)
𝑉𝑠𝑐 =
(1 − 𝑛)
Donde:
44
- Curva de gasto
En este método se deriva una relación entre los valores de concentración y el caudal de los
ríos en el momento del muestreo por una transformada logarítmica y regresión de mínimos
cuadrados (WALLING, D. & WEBB, B. 1981, 1988; WEBB, B. et al., 1997).
𝐶 = 𝑎𝑄 𝑏 (32)
45
Palabras claves:
- Río Ronquillo.- Es el tramo final de río Tres Ríos, conocido con el nombre de río
Ronquillo ya que adopta el nombre de la zona por el que pasa.
46
III. MATRIALES Y MÉTODOS
3.2. METODOLOGÍA
Para lograr los objetivos mencionados, se realizó la siguiente metodología
47
METODOLOGÍA
48
a. CARTOGRAFÍA DE LA CUENCA
Para obtener los parámetros geomorfológicos de la cuenca se utilizaron imágenes del radar
Alos Palsar de topografía con una resolución de 12.5 m correspondiente a la zona de interés,
las imágenes se procesaron utilizando el software QGis.
49
En la Tabla 5 se presentan algunos parámetros geomorfológicos de la cuenca después de realzar
la delimitación.
Luego de haber realizado el análisis de consistencia del registro de niveles de agua, se obtiene
los caudales para la estación, el cálculo de los caudales se realizó mediante la aplicación de la
fórmula de Poleni (PATT, H. & JÜPNER, R. 2013)(Ecuación 33).
2 (33)
𝑄= × 𝜇 × √2𝑔 × 𝑤 × ℎ3⁄2
3
50
De la ecuación anterior, “Q” representa el caudal (m³.s-1), “µ” representa un valor empírico
adimensional para describir la forma de la presa (µ = 0,75), g representa la gravedad terrestre
(m/s²), “w” representa el ancho de la cresta de la presa (w = 6,95 m) y “h” representa el nivel
del agua (m) sobre la cresta de la presa.
Figura 18. Piezómetro ubicado en estación Ronquillo, el cual registra los niveles de agua
cada un periodo de tiempo determinado.
51
VOLÚMENES DE SEDIMENTOS DE FONDO MEDIDO EN CAMPO.
Se registraron los eventos que llenaron la presa, el caudal que se tomó como base para
considerar que se inició el transporte de fondo fue el caudal dominante estimado, a partir de
ello se registró el tiempo de duración del evento y los caudales líquidos involucrados.
La compuerta ubicada en el corte C-C corresponde a la del canal de captación, es utilizada para
la regulación del caudal captado.
Las compuertas que se observan en los cortes A-A y B-B se utilizan para hacer la limpieza del
vaso de la presa, estas se abren cuando el vaso de la presa se llena de sedimentos para poder
realizar la limpieza, cuando la presa se ha quedado sin sedimento, las compuertas se cierran
hasta un próximo evento.
52
Figura 19. Vista en planta de captación Ronquillo y cortes de compuertas.
En las Figuras 20 y 21 se muestran los resultados de las barimetrías del totalizador o vaso de
la presa vacío y lleno respectivamente, los colores en cada una de la figuras representan las
elevaciones en metros sobre el nivel del mar (msnm). En las tablas de sus leyendas los colores
van enumerados desde el uno hasta seis, siendo uno el punto con sedimento más bajo y seis el
más alto.
54
Figura 20. Vista en planta de batimetría en época de estiaje de vaso de presa.
55
Figura 21. Vista en planta de batimetría de vaso de presa lleno.
56
- Volumen de sedimentos de fondo observado.
La capacidad de almacenamiento cuando se llena corresponde a 68.2 m3, hay que tener en
cuenta que el volumen estimado con las batimetrías incluyen sedimentos y poros.
Fueron 4 los eventos que llenaron el vaso de la presa, es decir, para cada uno de los cuatro
eventos se observó 68.2 m3 de sedimento acumulado. A continuación se presenta un gráfico de
frecuencias acumuladas, el cual nos permite observar el llenado del vaso de la presa y la
dinámica de los caudales involucrados en cada evento.
Figura 22. Cuatro eventos monitoreados que llenaron el vaso de la presa durante el año 2016.
50
45
40
35
30 Evento
25 09/02/2016
20 11/02/2016
09/03/2016
15
14/03/2016
10
5
0
2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5 10.5 11.5 12.5 13.5 14.5 15.5 16.5 17.5
Ql (m3/s)
57
En la Figura 23, el eje “y” representa al llenado de la presa 100% representa que la presa ya se
llenó (68.2 m3) y en el eje “x” se presenta los caudales involucrados en cado evento.
Para obtener la longitud del tramo con el que se trabajará se realizó una visita de campo en la
zona, observando las características del río Ronquillo se decidió trabajar con los 200 m del río
ubicados aguas arriba de la captación, ya que se consideró un tramo representativo tanto para
estimar el ancho promedio del río, pendiente y granulometría.
La línea azul central representa el alineamiento trazado para obtener la pendiente del cauce,
las progresivas del alineamiento tienen su punto de inicio en la progresiva 0+00 ubicada en la
captación Ronquillo y van aumentando aguas arriba de la captación Ronquillo finalizando en
la progresiva 0+200 Km.
58
Figura 23. Vista en planta y topografía de 200 m de río Ronquillo
59
Figura 24. Perfil longitudinal del cauce del río Ronquillo
60
La figura 24 corresponde al perfil longitudinal de río Ronquillo realizado a partir del
alineamiento mostrado en la Figura 26, el cual tiene un ancho promedio de 6.60 m, se
presentan las progresivas cada 10 m en el eje “x” y las pendientes de tramos entre
progresivas, además las cotas de elevación cada 1 metro en el eje “y”.
Donde:
𝐻𝑚á𝑥: Altura máxima.
𝐻𝑚𝑖𝑛: Altura mínima.
𝐿𝑜𝑛𝑔𝑖𝑡𝑢𝑑: Longitud del tramo.
2 (35)
𝑛
𝑆𝑓 =
1 1 1
+ + ⋯+
[√𝑆𝑓1 √𝑆𝑓2 √𝑆𝑓𝑛]
61
Donde:
n= Número de intervalos iguales.
Sfn= Pendiente en intervalo “n”.
Tabla 6. Resumen de resultados para pendiente del cauce del río Ronquillo.
Método Pendiente
Pendiente uniforme 3.161 %
Método de Taylor y
3.163 %
Schwarz
Promedio ponderado 3.162 %
Pendiente promedio 3.162 %
Para la obtención de los parámetros hidráulicos del río como son el área, perímetro mojado y
tirante del agua, se utilizó el software HEC-RAS.
62
Una vez realizada la calibración de la información, se ingresaron los caudales de los 4 eventos
monitoreados y se obtuvieron sus respectivos parámetros.
Se utilizó el método propuesto por Cowan (1956) con el cual se estimó un valor inicial del
coeficiente de rugosidad, considerando los siguientes valores:
63
Para calibrar el valor del coeficiente de manning se utilizó el programa HEC RAS, teniendo
caudales medidos y simulados con la topografía del río se realizaron variaciones de coeficiente de
manning, obteniéndose un coeficiente de rugosidad mostrado en la Tabla 7.
Tabla 7.
Coeficiente de rugosidad para cauce de río Ronquillo.
Coeficiente de rugosidad
n= 0.052
El muestreo para este tipo de material se realizó utilizando el método volumétrico, este método
es el más completo de todos y consiste en tomar del lecho un cierto volumen de material sub-
superficial, con ese propósito se retira primero la capa superficial en un espesor comparable al
tamaño de la mayor partícula observada en la superficie. El volumen que se toma a
continuación debe ser representativo del material granular del cauce, para lo cual debe seguirse
el criterio de la mayor partícula observada la cual no represente más del 1% en peso total de la
muestra utilizando la siguiente fórmula (Vide, 2007).
𝑃 = 𝐷𝑚á𝑥 ∗ 10 (36)
Donde:
Teniendo en cuenta las consideraciones anteriores, se observó que el diámetro mayor fue de
47 cm aproximadamente, lo cual nos indica que se deben extraer 470 Kg del lecho sub-
superficial.
64
Para la presente tesis se ubicaron 3 puntos de muestreo a lo largo de los 200 m del río Ronquillo
para la cual se realizó su respectiva curva granulométrica.
El método utilizado para este análisis fue el de tamizado, se utilizaron los tamices
recomendados por la Norma ASTM D-422.
65
Para el análisis del material extraído del cauce del río Ronquillo, se utilizaron los tamices desde
el de 3” hasta el tamiz Nº 200 como se observa en el eje secundario “x”.
Figura 27. Porcentaje de material retenido (%) sub – superficial para cada punto de
muestreo, respecto a su diámetro (mm).
Tabla 8.
Porcentaje de material del cauce.
Material Porcentaje
Finos 0.40%
Arenas 37.90%
Suelo Gravoso
66
Tabla 9.
Clasificación del material del lecho del río Ronquillo.
El análisis nos da como resultado que el lecho del río Ronquillo es un lecho gravoso con arenas
y material bien gradado. Predominan en mayor cantidad las gravas con diámetros desde 2 mm
hasta los 25 mm y en un menor porcentaje las gravas de diámetro mayor a 37.5 mm. Respecto
a las arenas, existe la mayor predominancia de arenas de diámetro de 0.850 mm y en un menor
porcentaje de arenas menores a 0.355 mm.
Tabla 10.
Diámetros característicos de material del lecho sub - superficial del río Ronquillo.
Diámetros
Valor
característicos
D16 16% 2.1 mm
D50 50% 8.8 mm
D84 84% 31.6 mm
D90 90% 39.8 mm
Dm 10.95 mm
67
Densidad y peso específico:
Para determinar la densidad del material fino y grueso (material que pasa la malla número 4 y
el material que se retiene en la malla número 4 respectivamente) se realizó tomando en cuenta
los procedimientos recomendados en la norma ASTM C 127, obteniéndose los siguientes
resultados:
Tabla 11.
Densidad y peso específico del material del lecho del río Ronquillo.
Para el muestreo del material superficial se utilizó el método fotográfico (RAMÍREZ, C. et al.,
2011), el cual es uno de los métodos utilizados cuando el material presente corresponde a
diámetros muy grandes y se hace complicado su traslado hacie el laboratorio.
68
Figura 28. Cuatro puntos de muestreo seleccionados en el lecho superficial del río Ronquillo
(a) (b)
(c) (d)
Las figuras (a), (b), (c) y (d) corresponden a los puntos de muestreo seleccionados de los 6
realizados a los cuales llamaremos punto de muestreo 1, 2, 3 y 4 respectivamente.
69
Análisis granulométrico del material del lecho superficial.
Para la clasificación de este tipo de material se consideraron los siguientes diámetros: menores
4.75 mm, 6.30 mm, 7.50 mm, 12.50 mm, 19.00 mm, 25 mm, 35.70 mm, 50 mm, 75 mm, 100
mm, 150 mm, 200 mm, 250 mm, 400 mm y mayores a 500 mm. Después del cálculo de las
áreas del material analizado se obtuvo los siguientes resultados.
Se observa que el punto de muestreo 3 presenta mayor tamaño de diámetros mientras que los
puntos 1,2 y 5 tienen diámetros parecidos.
70
Tabla 12.
Diámetros característicos de material superficial del lecho del río Ronquillo.
Diámetros
Valor
característicos
D10 10% 20.34 mm
D 16 16% 29.59 mm
D50 50% 94.73 mm
D84 84% 343.84 mm
D90 90% 446.22 mm
Dm 129.48 mm
Se observa que el diámetro promedio del material del lecho superficial de lecho es de 130 mm
aproximadamente.
Figura 30. Material retenido (%) en cada punto de muestreo del material lecho superficial
71
En la Figura 31 se observa que el material predominante en el lecho superficial del río varía
entre los 35 mm hasta los 150 mm, es por ello que para este tipo de sedimento se seleccionó el
método fotográfico el cual nos permitió realizar su clasificación.
Una vez obtenida la información mediante los trabajos de campo y de gabinete se procedió a
la aplicación de las fórmulas empíricas para la estimación del sedimento de fondo.
Como primer paso se estimaron los caudales sólidos para cada uno de los cuatro eventos
utilizando las siguientes fórmulas:
72
Tabla 13.
Fórmulas utilizadas para la estimación del volumen de sedimento acumulado.
FÓRMULA SIGLAS DIÁMETROS RESTRICCIONES
Sub-sustrato basado P-K-Mc D50sub Adecuado para pequeñas y
en la ecuación de medianos canales que
Parker, Klingteman y tienen tamaños de partícula
McLean 1982 entre 0,6 a 102 mm, con
pendientes medias
Smart y Jaeggi 1983 S-J D90 sub, D30 sub Válida para diámetros
mayores a 0.4 mm y
pendientes entre 0.4% y
20%
Bathurst 1987 Bathurst D50sub Ó D16sup Válida para pendientes
desde 0.4% hasta 4% y
cumple la condición 0.3 ≤
d/D84 ≤ 50
Meyer –Peter y M-P-M Dm, D90sup Válida para diámetro
Muller medio entre 0.4 y 30 mm
Donde:
73
Para estimar la porosidad del material se realizó en campo 3 puntos de muestreo a partir de los
cuales se obtuvo un valor promedio correspondiente al 30%, el cual se utilizó para realizar la
corrección por porosidad como se mencionó en el párrafo anterior.
d. SEDIMENTO EN SUSPENSIÓN.
MONITOREO
El monitoreo se realizó en todos los puntos en donde están ubicadas la estación hidrométrica
Ronquillo, el trabajo de muestreo fue realizado por cada observador designado y debidamente
capacitado.
PROCESAMIENTO DE MUESTRAS
El método que se aplicó es el de “Curva de gasto”, en este método se deriva una relación entre
los valores de concentración y el caudal del río en el momento del muestreo por una
transformada logarítmica y regresión de mínimos cuadrados (WALLING, D. & WEBB, B.
1981, 1988; WEBB, B. et al., 1997).
𝐶 = 𝑎𝑄 𝑏 (37)
A partir de este método se han desarrollado diversas variantes que consisten básicamente en
el agrupamiento de los datos teniendo en cuenta la estación, año, meses, etc.
75
Figura 31. Concentraciones de sedimentos en suspensión (SSC) y caudales líquidos (Q)
medidos en la estación Ronquillo.
10000
1000
y = 126.62x1.3129
R² = 0.42
SSC (mg/L)
100
10
0
0.01 0.10 1.00 10.00
Q (m3/s)
76
IV. ANÁLISIS Y DISCUCIÓN DE RESULTADOS
4.1. CARTOGRAFÍA
La cuenca en estudio tiene un área correspondiente a 41.7 km2, con una variación de alturas
desde los 3508 msnm hasta los 3986 msnm, tiene una pendiente de laderas promedio de 29%
y la pendiente promedio del cauce del río Ronquillo corresponde es de 3.2%.
Tabla 14.
Principales parámetros geomorfológicos de la cuenca Ronquillo.
Parámetros Valor
Área de la cuenca, Km2 41.7
Perímetro de la cuenca, Km 34.5
Pendiente promedio de la cuenca, % 29.0
Elevación máxima, msnm 3986.0
Elevación media, msnm 3508.0
Elevación mínima, msnm 2835.0
Tiempo de concentración, mín 235.0
Longitud del máximo recorrido, Km 13.2
Pendiente media del río ubicado aguas arriba
3.2
de la captación, %
Precipitación anual, mm 1280
Descarga máxima medida, m3 s-1 33.7
Descarga media medida, m3 s-1 0.2
Descarga mínima medida, m3 s-1 0.1
- Según la pendiente del cauce de río ubicado aguas arriba de la captación Ronquillo y
tomando en cuenta que los ríos de montaña que tienen gradientes longitudinales que
77
oscilan entre 0.1% y 20% (PAPANICOLAOU, A. et al., 2004), se puede considerar al
río Ronquillo como un río de montaña.
3800
Elevación (m.s.n.m)
3600
3400
3200
3000
2800
0 20 40 60 80 100
Porcentaje de Área (%)
- La curva hipsométrica presentada es una curva sinusoidal, según este tipo de curva se
puede interpretar y caracterizar a la cuenca como una cuenca que se encuentra en
equilibrio y se encuentra en fase de madures.
a. CAUDALES LÍQUIDOS.
Los 4 eventos monitoreados son considerados como eventos extremos en el año 2016, ya que
lograron llenar el vaso de la presa (totalizador) con el sedimento transportado durante su
ocurrencia.
78
Figura 33. Eventos de caudales extremos que llenaron el barraje.
10.5
10.0
9.5
9.0
8.5
8.0
7.5
7.0
6.5
6.0
5.5
5.0
4.5
4.0
3.5
3.0
2.5
2.0
00 60 120 180 240 300 360 420
Duración (minutos)
79
De los cuatro eventos presentados se observa que los que ocurrieron en las fechas de 11/02/16
y 14/03/16 presentan los picos más altos, mientras los ocurridos el 09/02/16 y 09/03/16
presentan picos menores y comportamientos similares.
Conocer los aportes de sedimentos en suspensión de la cuenca obteniendo una relación entre
las concentraciones de sedimentos en suspensión y los caudales líquidos es un gran aporte, es
un gran reto, debido a la aleatoriedad de las variables involucradas. Sin embargo en la cuenca
Ronquillo se han obtenido los primeros resultados del monitoreo puntual en la estación
Ronquillo.
80
Figura 34. Serie temporal de concentraciones de sedimentos Vs. Caudales líquidos en la estación Ronquillo.
Concentración de sedimento en suspensión (SSC) Vs Caudale líquido (Q)
Estación Ronquillo
1800.00 3.60
1700.00 3.40
1600.00 3.20
1500.00 3.00
1400.00 2.80
1300.00 2.60
1200.00 2.40
1100.00 2.20
SSC (mg/l)
Q m3/s
1000.00 2.00
900.00 1.80
800.00 1.60
700.00 1.40
600.00 1.20
500.00 1.00
400.00 0.80
300.00 0.60
200.00 0.40
100.00 0.20
0.00 0.00
15/02/2016 15/03/2016 15/04/2016 15/05/2016 15/06/2016 15/07/2016 15/08/2016
Fechas
81
Figura 35. Correlación entre concentraciones de sedimentos en suspensión (SSC) y caudales
líquidos (Q) medidos en la estación Ronquillo.
1000
y = 41.315x1.7211
R² = 0.45
y = 0.2409x-1.011
100 R² = 0.02
SSC (mg/L)
10
0
0.01 0.10 1.00 10.00
Q (m3/s)
Para los caudales iguales o mayores a 0.10 m3.s-1se presenta un coeficiente de determinación
de R2= 0.45 con coeficientes a= 41,315 y b= 1,72; el coeficiente “a” es la constante de la
ecuación, mientras que “b” representa la intensidad de las fuerzas erosivas, para el presente
cado el exponente “b” es mayor que 1, ello indica que las concentraciones crecen cuando las
descargas se incrementan (efecto de purga).
82
concentraciones decrecen rápido con el incremento en la descarga, por lo tanto el transporte
declina cuando la descarga aumenta.
Los volúmenes obtenidos utilizando las fórmulas empíricas para cada evento de crecida se
pueden observar en la Figura 36.
Figura 36. Volumen de sedimento acumulado (m3) para cada evento de crecida utilizando las
fórmulas empíricas.
1500
1000
500
Ve (m3)
0
09/03/2016 14/03/2016
2000
1500
1000
500
0
P-K-Mc S-J Bathurst M-P-M
Método
Los resultados obtenidos muestran que la formula con la que se obtienen los valores más bajos
de volúmenes de sedimentos en los cuatro eventos monitoreados es la de P-K-Mc “Sub-
sustrato basado en la ecuación de Parker, Klingteman y McLean”.
83
En cuanto a las fórmulas con la que se estimó los valores más altos, varían de acuerdo al evento.
Se presentan 2 grupos bien marcados, para los eventos ocurridos el 9 de febrero y 9 de marzo
los volúmenes estimados más altos corresponden a las fórmulas de S-J (Smart y Jaeggi) y M-
P-M (Meyer Peter y Muller), mientras que en los eventos ocurridos el 11 de febrero y 14 de
marzo, los valores más altos se obtuvieron con la fórmula de Bathurst.
La variación de resultados obtenidos con cada fórmula se debe a que los eventos monitoreados
son distintos (caudales, tiempo de duración), lo que nos permite observar la sensibilidad de
cada fórmula aplicada.
Los valores de los volúmenes acumulados obtenidos tanto estimados como observados,
muestras en la Tabla 15, a partir de ella se aprecia la gran variabilidad de los resultados
obtenidos mediante fórmulas empíricas respecto al valor medido en campo.
Tabla 15.
Tabla comparativa de volúmenes acumulados de sedimentos de fondo estimados Vs
volúmenes observados.
Evento Volumen de Sedimento de fondo m3
09/02/2016 11/02/2016 09/03/2016 14/03/2016
Método
P-K-Mc 641.96 463.78 597.43 700.21
S-J 1797.73 1174.21 1780.46 1883.29
Bathurst 1623.87 1673.92 1617.24 2061.96
M-P-M 1795.90 1290.16 1669.16 1954.18
Observado 68.20 68.20 68.20 68.20
Gráficamente (Figura 37) se aprecia que para los eventos monitoreados las fórmulas están
sobre-estimando los volúmenes reales.
84
Figura 37. Volúmenes acumulados de sedimentos de fondo Estimados y Observados.
c J rst do
-M S- -M
hu -P va
P-
K
Ba
t M ser
b
O
Método
Para poder estimar cuanto es la cantidad que sobre estima cada fórmula empírica se analizó los
coeficientes de proporcionalidad para cada una de las 4 fórmulas (Figura 38).
La fórmula empírica que sobre estima en menor cantidad al volumen de sedimento de fondo
observado corresponde a la fórmula de P-K-Mc “Sub-sustrato basado en la ecuación de Parker,
Klingteman y McLean”, con una mediana de 9.1, es decir que esta fórmula está estimando 9.1
veces el volumen de sedimento de fondo observado.
En las fórmulas de S-J (Smart y Jaeggi) y Bathurst se están presentando valores extremos, ello
representa la sensibilidad de estas fórmulas en cuanto a la variación de caudales (caudales
extremos y duración del evento).
La fórmula de M-P-M (Meyer Peter y Muller), presenta una mediana de 26.3 veces la sobre-
estimación del sedimento de fondo y una amplia distribución de los valores obtenidos
85
Figura 38. Coeficiente de proporcionalidad entre volumen estimado (Ve) y volumen medido
en campo (Vm).
86
V. CONCLUCIONES Y RECOMENDACIONES.
5.1. CONCLUCIONES.
- La cuenca ronquillo tiene un área de 41.7 km2, con una variación de alturas desde los
2835 msnm hasta los 3986 msnm, la pendiente de laderas promedio es de 27.92% y
pendiente promedio del tramos del cauce del río Ronquillo seleccionado es de 3.16%.
Por la magnitud de la pendiente promedio del río, este es considerado como un río de
montaña y según la pendiente promedio de ladera se considera que es una cuenca con
relieve accidentado.
87
5.2. RECOMENDACIONES.
- Se debe tener en consideración que con eventos que presenten caudales mayores al
caudal dominante (2.8 m3.s-1) pueden presentarse problemas de colmatación por
arrastre de sedimentos de fondo en la captación Ronquillo.
- Se debe tener en cuenta que los valores de sobre estimación de las fórmulas empíricas
utilizadas en el presente estudio solo son válidas para eventos similares a los mostrados,
ya que para eventos extremos los resultados pueden variar.
88
VI. BILBIOGRAFÍA
6.1. LIBROS
Alonso, R. L. (2004). Problemas resueltos de ingeniería hidráulica forestal (Vol. 47): Universitat
de Lleida.
Charlton, R. (2007). Fundamentals of fluvial geomorphology: Routledge.
Chow, V. T. (1983). Hidráulica de Canales Abiertos: ISBN 968-13-1327.
Das, B. M. (2001). Fundamentos de ingeniería geotécnica: Thomson Learning.
Juan Pedro Martín, V., & Upc Edicions, U. (2007). Ingeniería de ríos (Vol. 9): Univ. Politéc. de
Catalunya.
Mijares, A., & Mijares, F. J. F. J. A. (1989). Fundamentos de hidrología de superficie: Limusa.
Miller, J. R., Mackin, G., & Miller, S. M. O. (2015). Application of geochemical tracers to
fluvial sediment: Springer.
Reyes, A., Barroso, F. U., & Carvajal, Y. (2010). Guía Básica para la caracterización
morfométrica de cuencas hidrográficas: Programa Editorial UNIVALLE.
Rocha, A. (1998). Introducción a la hidráulica fluvial: Universidad Nacional de Ingeniería.
Facultad de Ingeniería Civil.
Rubio, T. O. (2011). Hidráulica de ríos y procesos morfológicos: Ecoe ediciones.
Vide, J. P. M. (2007). Ingeniería de ríos (Vol. 9): Univ. Politéc. de Catalunya.
Villón, M. (2002). Hidrología. Instituto Tecnológico de Costa Rica.
Hudson, N. (1997). Medición sobre el Terreno de la Erosión del Suelo y de la Escorrentía (Vol.
68): Food & Agriculture Org.
6.2. ARTÍCULOS
Aguirre-Pe, J., Olivero, M., & Moncada, A. (2000). Transporte de sedimentos en cauces de alta
pendiente. Ingeniería del agua, 7(4), 353-365.
Barnes, H. H. (1967). Roughness characteristics of natural channels: US Govt. Print. Off.
Barry, J. J., Buffington, J. M., Goodwin, P., King, J. G., & Emmett, W. W. (2008). Performance
of bed-load transport equations relative to geomorphic significance: Predicting effective
89
discharge and its transport rate. Journal of Hydraulic Engineering, 134(5), 601-615. doi:
https://doi.org/10.1061/(ASCE)0733-9429(2008)134:5(601)
Cowan, W. L. (1956). Estimating hydraulic roughness coefficients. Agricultural Engineering,
37(7), 473-475.
Elsitdié, L. G. C., Vide, J. P. M., & Martín, M. D. M. (Writers). (2009). Coeficientes de
resistencia, transporte de sedimentos y caudal dominante en regiones semiáridas, Todo
sobre Coeficientes de resistencia, transporte de sedimentos y caudal dominante en
regiones semiáridas.
Farias, H. D. (2005). Análisis del concepto de caudal dominante en ríos aluviales y evaluación de
metodologías de cálculo. Paper presented at the Segundo Simposio Regional sobre
Hidráulica de Ríos.
García, C., & Sala, M. (1998). Aplicación de fórmulas de transporte de fondo a un río de gravas:
comparación con las tasas reales de transporte obtenidas en el río Tordera. Ingeniería del
agua, 5(1), 59-72.
Glysson, G. D. (1987). Sediment-transport curves: US Geological Survey.
Gray, J. R., Laronne, J. B., Osterkamp, W., & Vericat, D. (2010). Bed Load Research
International Cooperative–BRIC. US Geological Survey Scientific Investigations Report,
5091, 208-217.
Gray, J. R., & Simões, F. J. (2008). Estimating sediment discharge Sedimentation engineering:
Processes, measurements, modeling, and practice (pp. 1067-1088).
Habersack, H., Kreisler, A., Rindler, R., Aigner, J., Seitz, H., Liedermann, M., & Laronne, J. B.
(2016). Integrated automatic and continuous bedload monitoring in gravel bed rivers.
Geomorphology, 291, 80-93. doi: http://dx.doi.org/10.1016/j.geomorph.2016.10.020
Haddadchi, A., Omid, M. H., & Dehghani, A. A. (2013). Bedload equation analysis using bed
load-material grain size. Journal of Hydrology and Hydromechanics, 61(3), 241-249.
López, R. (2003). Fórmulas para el cálculo aproximado de la capacidad de transporte de
sedimentos en ríos de montaña. Cimbra, 2003, núm. 351, p. 32-36.
López, R., Vericat, D., & Batalla, R. J. (2015). Evaluación de fórmulas de transporte de fondo en
un río de gravas acorazado. Tecnología y ciencias del agua, 6(2), 5-20.
Nitsche, M., Rickenmann, D., Turowski, J. M., Badoux, A., & Kirchner, J. W. (2011).
Evaluation of bedload transport predictions using flow resistance equations to account for
90
macro‐roughness in steep mountain streams. Water Resources Research, 47(8). doi:
10.1029/2011WR010645
Ortiz, M., Huerta-Tamayo, L., & Hinojosa, A. (2003). Transporte de sedimento por tracción de
marea en el Estero de Punta Banda, Baja California, México. GEOS, 23(3), 283-294.
Owens, P., Batalla, R., Collins, A., Gomez, B., Hicks, D., Horowitz, A., Peacock, D. (2005).
Fine-grained sediment in river systems: environmental significance and management
issues. River research and applications, 21(7), 693-717. doi: 10.1002/rra.878
Papanicolaou, A., Bdour, A., & Wicklein, E. (2004). One-dimensional hydrodynamic/sediment
transport model applicable to steep mountain streams. Journal of Hydraulic Research,
42(4), 1-19. doi: 10.1080/00221686.2004.9728402
Parker, G., Klingeman, P. C., & McLean, D. G. (1982). Bedload and size distribution in paved
gravel-bed streams. Journal of the Hydraulics Division, 108(4), 544-571.
Patt, H., & Jüpner, R. (2013). Hochwasser-Handbuch. Auswirkungen und Schutz, Berlin,
Heidelberg. doi: https://doi.org/10.1007/978-3-642-28191-4
Perks, M. T. (2014). Suspended Sediment Sampling. British Society for Geomorphology.
Phillips, J., Webb, B., Walling, D., & Leeks, G. (1999). Estimating the suspended sediment loads
of rivers in the LOIS study area using infrequent samples. Hydrological processes, 13(7),
1035-1050.
Phillips, J. M., Russell, M. A., & Walling, D. E. (2000). Time-integrated sampling of fluvial
suspended sediment: a simple methodology for small catchments. Hydrological
Processes, 14(14), 2589-2602.
Ramírez, C., García, J. L., Ramirez, C. A., & Bocanegra, R. A. (2011). Caracterización de los
materiales del lecho del río Cauca. Ingeniería de Recursos Naturales y del Ambiente,
1(1), 10-17.
Rascher, E., Baewert, H., Schmidt, K.-H., & Morche, D. (2012). Comparing the results of bed
load transport equations to field measurements in an Alpine river. Paper presented at the
EGU General Assembly Conference Abstracts.
Smart, G. M. (1984). Sediment transport formula for steep channels. Journal of Hydraulic
Engineering, 110(3), 267-276. doi: 10.1061/(ASCE)0733-9429(1984)110:3(267)
91
Smith, T. B., & Owens, P. N. (2014). Flume and field based evaluation of a time-integrated
suspended sediment sampler for the analysis of sediment properties. Earth Surface
Processes and Landforms, 39(9), 1197-1207.
Vanoni, V. A. (2006). Sedimentation Engineering: Theory. Measurements, Modeling, and
Practice, Manuals and Reports on Engineering Practice (54).
Wren, D. G., Barkdoll, B. D., Kuhnle, R. A., & Derrow, R. W. (2000). Field techniques for
suspended-sediment measurement. Journal of Hydraulic Engineering, 126(2), 97-104.
6.3. OTROS
Benavides, I., Ángeles, I., Salazar, E., & Abásolo, J. (2007). Inventario de Fuentes de Agua
Superficial de la Cuenca del Mashcon.
Cruzado, G. Proyecto Gestión de Riesgo de Deslizamiento del Ronquillo, Urubamba y
Corisorgona.
Chuquiruna , M., Cotrina , S., & Sánchez , I. (2008). Diagnóstico Socio Ecológico de Nuestra
Jalca en Microcuenca en Ronquillo. Proyecto Páramo Andino.
Municipalidad Provincial de Cajamarca. (2012). Estudio de Evaluación de Riesgo de desastres
del Sector Crítico Urbano: Urubamba-Ronquillo. Municipalidad Provincial de
Cajamarca.
Ortiz, O. (2004). Evaluación hidrológica. Revista Hidrored, 2-7.
Olazábal, J. B. Q. (2015). Estudio esperimental del transporte de sedimentos en suspensióm y
fondo y comparación con fórmulas empíricas en los ríos Puyango-Tumbes y Zarumilla.
[TESIS].
Vásquez, G. C., & Rodríguez, M. E. C. (2011). Estudio de Geología. Gobierno Regional de
Cajamarca.
92
VII. ANEXOS.
(a) (b)
Figura 39. (a) Levantamiento topográfico del cauce de río Ronquillo y (b) batimetría del
vaso formado por el barraje de la captación Ronquillo.
(a) (b)
Figura 40. (a) Material del lecho superficial y (b) material del lecho sub-superficial del río
Ronquillo.
93
ANEXO 2: TRABAJOS DE LABORATORIO.
(a) (b)
(c) (d)
Figura 41. (a) Tamizado de material sub-superficial, (b) Material para densidad de gruesos,
(c) y (d) Ensayo para determinar densidad de material fino.
94
(a) (b)
Figura 42. (a) Filtrado d muestras de sedimentos en suspensión, (b) Pesado de muestras
filtros secos para obtención de concentración de sedimentos en suspensión.
95
ANEXO 3: PRESUPUESTO.
Tabla 16
Presupuesto para elaboración de tesis.
96